• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    波流作用下圓柱體入水特性的三維數(shù)值模擬研究

    2015-04-26 02:56:33孫龍泉姚熊亮
    船舶力學(xué) 2015年10期
    關(guān)鍵詞:波流圓柱體質(zhì)心

    楊 衡,孫龍泉,劉 瑩,姚熊亮

    (哈爾濱工程大學(xué) 船舶工程學(xué)院,哈爾濱 150001)

    波流作用下圓柱體入水特性的三維數(shù)值模擬研究

    楊 衡,孫龍泉,劉 瑩,姚熊亮

    (哈爾濱工程大學(xué) 船舶工程學(xué)院,哈爾濱 150001)

    結(jié)構(gòu)在實際海況中入水受到多種載荷的共同作用,同時還伴隨波浪作用的影響,因此該過程是一個強非線性的過程。針對結(jié)構(gòu)在波、流中入水過程的特點,將入射波(波、流)引入非線性雙漸進法,研究三維剛體圓柱體在波、流及波流聯(lián)合作用下入水過程運動響應(yīng)及姿態(tài)的變化,計算結(jié)果與試驗結(jié)果符合得較好,非線性雙漸進法適用于分析三維剛體波、流中入水問題。結(jié)果表明在近波面附近,結(jié)構(gòu)受波浪作用明顯,入水相位、浪級、流速及波流速度矢量差異對結(jié)構(gòu)入水運動速度及軌跡影響明顯。

    入水;波流;非線性;雙漸進法(DAA)

    0 引 言

    結(jié)構(gòu)入水過程中會產(chǎn)生砰擊現(xiàn)象,在真實的海況條件下,還要受到波浪、洋流的影響,結(jié)構(gòu)入水過程中運動響應(yīng)及姿態(tài)與波、流的傳播方向、浪級和入水相位[1]有關(guān),對結(jié)構(gòu)入水過程的安全性及入水軌跡的穩(wěn)定性非常重要。結(jié)構(gòu)在波、流中入水,除受到浮力、結(jié)構(gòu)自身重力、砰擊載荷外,還要考慮到波、流載荷對結(jié)構(gòu)的影響,同時還伴隨波、流的運動及波面的非線性變化,因此該過程是一個強非線性瞬態(tài)流固耦合過程,包括濕表面、自由面邊界條件和伯努利方程的非線性。

    目前對入水砰擊問題研究較多,研究對象以二維楔形體為主,自由液面多為自由平面。在國內(nèi)外有關(guān)物體入水問題的研究文獻中,真正考慮波浪因素的并不多,其中,F(xiàn)altinsen[2]采用附加質(zhì)量法計算入水物體所受水動力,基于勢流邊界元理論得出附加質(zhì)量隨浸深變化的關(guān)系曲線,并考慮入射波的作用,但未考慮入水物體對波浪的反作用。王文華等[3]采用二維CFD方法數(shù)值模擬了圓柱在規(guī)則波中的入水過程。儲慧林[4]建立了二維數(shù)值水槽,重點研究了波浪相位及波高對魚雷砰擊壓力、彈道和速度的影響。袁緒龍等[5]研究了具有細長前錐段外形的超空泡導(dǎo)彈在高速入水時彈道和流體動力的情況,通過流場-彈道耦合方法,分析了不同預(yù)置舵角下入水彈道與流體動力的變化規(guī)律。對于波浪因素的研究,權(quán)曉波等[6]基于波浪理論得到了二階Stokes波的速度表達式,并通過給定入射邊界速度分布形式的方法實現(xiàn)波浪的數(shù)值模擬,利用動網(wǎng)格技術(shù)、UDF技術(shù)等,探討了水下航行體模型出水時波浪對其流場、力學(xué)特性等的影響。Wang和Wu[7]基于FEM方法求解了波浪與垂直固定圓柱體全非線性相互作用過程。劉云龍、張阿漫等[8]將波浪作用以入射勢方式引入水下爆炸氣泡計算中,得到氣泡脈動與波浪面相互作用及波浪面的變形。對于潛射導(dǎo)彈出水、魚類入水、水下爆炸氣泡脈動等瞬態(tài)流固耦合領(lǐng)域,考慮流體的可壓縮性是非常必要的,王詩平、孫士麗等[9]對可壓縮流場中氣泡脈動特性進行數(shù)值模擬,得出可壓縮性對氣泡射流速度的影響。王文華[10]利用液面捕捉法和直角切割網(wǎng)格系統(tǒng)解決入水過程中瞬時移動的自由液面和動邊界問題。

    本文的基于二階雙漸進法[11-13](DAA2),推導(dǎo)了考慮入射波作用下二階DAA方程,通過非線性伯努利方程[14],將二階DAA方程與結(jié)構(gòu)運動方程耦合求解,討論了剛體圓柱入水過程中運動特性,與模型試驗進行對比,驗證本文方法可靠性。著重探討了波、流作用下三維結(jié)構(gòu)入水特性,為進一步研究海洋結(jié)構(gòu)物在波浪中入水問題提供了方法。

    1 數(shù)值計算方法

    水中的結(jié)構(gòu)物與流場的瞬態(tài)耦合問題有這樣的特點:早期瞬態(tài)(高頻)響應(yīng)主要是聲輻射問題,即所謂平面波或曲面波近似解問題;后期瞬態(tài)(低頻)響應(yīng)主要是“虛質(zhì)量”問題。Geers等人(1975-1980),根據(jù)這些特點提出了雙漸近法(DAA)。雙漸近法(DAA)之所以能廣泛應(yīng)用于水中結(jié)構(gòu)物瞬態(tài)流固耦合領(lǐng)域,其中一個重要原因是它考慮了流體可壓縮性,假設(shè)流場為各向同性、無粘、無旋但可壓縮的理想流體,則考慮可壓縮性的勢函數(shù)微分方程式:

    建立結(jié)構(gòu)在波、流中入水?dāng)?shù)值計算模型如圖1所示,結(jié)構(gòu)以速度V在重力作用下垂直入水,波、流向右傳播。

    結(jié)構(gòu)運動的邊界條件為:

    式中:V為結(jié)構(gòu)運動速度;n是法向矢量坐標(biāo);Φi為波浪產(chǎn)生的速度勢。

    無窮遠的邊界條件為:

    當(dāng)考慮自由液面效應(yīng)時,流場在滿足控制方程以及物面不可穿透條件的同時,還需滿足自由液面邊界條件:

    運動學(xué)邊界條件:

    圖1 數(shù)值模型Fig.1 Calculate model

    動力學(xué)邊界條件:

    式中:g為重力加速度,ξ為自由液面的波幅。

    在延遲勢法基礎(chǔ)上,滿足上述邊界條件,求解控制方程(1),得到速度勢方程如下:

    in標(biāo)r0的距離。

    引入了一個線性假設(shè),P=(ρΦ˙r,)t,得到考慮入射波作用下雙漸近方程如下:

    其中:P=Ps+Pb,Ps為流體中的散射壓力,Pb為入射波壓力。

    將方程(6)對時間進行一次積分

    在得到流場速度勢的基礎(chǔ)上,可以通過非線性伯努利方程[11]得到流場壓力計為Pd。

    將計算所得的動壓力加載到結(jié)構(gòu)上,得結(jié)構(gòu)動響應(yīng)。

    結(jié)構(gòu)運動速度ω˙n以及結(jié)構(gòu)加速度ω¨n作為下一次計算的初值,從而完成下一次循環(huán)計算。

    2 波浪作用下結(jié)構(gòu)入水特性

    2.1 波浪基本理論

    波浪是海面受風(fēng)壓擾動后,在重力作用下產(chǎn)生的周期性運動,粘性的作用很小。當(dāng)模型尺度相對于波浪波長較小,不規(guī)則波浪可簡化為規(guī)則的三維平面進行波。波浪運動采用空間固定坐標(biāo)系OXYZ,波浪向X軸正向傳播,如圖2所示。

    波數(shù)k=2π/λ,頻率ω=2π/T,T為波浪周期。

    線性波速度勢為:

    圖2 三維波浪示意圖Fig.2 The 3D wave

    2.2 模型試驗結(jié)果對比分析

    應(yīng)用本文方法計算,選取與文獻[11]具有相同的質(zhì)量,質(zhì)心位置,轉(zhuǎn)動慣量計算模型,計算該圓柱縮比模型入水過程中運動參數(shù)變化,并與試驗值進行對比分析。圖3給出結(jié)構(gòu)入水過程中運動速度及姿態(tài)計算值與試驗值對比結(jié)果。

    圖3是圓柱入水過程中三個方向的速度變化曲線,V1是橫向速度、V2法向速度、V3是垂向速度,從圖中可以看出,結(jié)構(gòu)入水過程中受到的垂向載荷變化劇烈,結(jié)構(gòu)的垂向加速度變化劇烈,相比于垂向速度和法向速度,橫向速度的變化曲線最為平緩。

    圖4、5是結(jié)構(gòu)入水過程中俯仰角、偏航角的變化曲線,從中可以看出,在圓柱入水過程中,在砰擊載荷和波浪載荷的作用下,圓柱俯仰角減??;偏航角增大,即圓柱在入水過程中首部受力和尾部受力的不平衡產(chǎn)生了偏轉(zhuǎn)力矩,引起了圓柱運動姿態(tài)的改變。

    圖3 速度時歷曲線Fig.3 The history curve of the velocities

    圖5 偏航角時歷曲線Fig.5 The history curve of the yaw angle

    圖6 節(jié)點分布Fig.6 The distribution of note

    2.3 波浪作用下圓柱入水運動特性分析

    選取波浪浪級為5級,取上述計算模型,垂直入水速度為25m/s,水平速度為0,入水位置為波谷,波浪向右傳播。取圓柱體三個特征點:頭部、尾部中心點、質(zhì)心進行,如圖6所示。

    計算過程中忽略波面變形及氣墊效應(yīng)對圓柱入水特性影響,計算結(jié)果如圖7和圖8所示。

    由圖7可以看出,圓柱在波浪中入水過程三個節(jié)點的運動速度呈現(xiàn)明顯的差異,頭部中心點(節(jié)點5)首先入水,受到波浪作用,頭部質(zhì)點速度增加,隨著侵水深度的增加,波浪作用力減小,頭部速度增加變緩。在圓柱質(zhì)心入水前,圓柱頭部處于波浪作用,產(chǎn)生相對于質(zhì)心逆時針彎矩,使圓柱尾部產(chǎn)生反向速度;圓柱質(zhì)心入水后,一方面圓柱尾部逐漸侵水,受到波浪作用,另一方面圓柱頭部遠離波浪面,波浪作用力減小,兩者共同作用結(jié)果產(chǎn)生繞圓柱順時針彎矩,圓柱尾部速度逐步恢復(fù)至零值。圓柱質(zhì)心速度與圓柱頭部質(zhì)心速度變化趨勢相近。從圖8可以看出,波浪作用對圓柱垂向速度影響較小,在圓柱頭入水和尾入水瞬間,垂向速度曲線產(chǎn)生明顯的波動。由以上分析,可將圓柱體在波浪中入水過程分為三個明顯的過程:

    (1)頭入水—質(zhì)心入水:圓柱頭部中心、質(zhì)心速度正向增加,尾部速度負向增加;

    (2)質(zhì)心入水—尾入水:圓柱頭部中心、質(zhì)心速度正向增加,尾部速度負向減??;

    (3)尾入水—:尾入水時刻,圓柱三點的速度出現(xiàn)波動,隨著圓柱繼續(xù)下沉,波浪的作用減弱,各點速度變化平穩(wěn)。由于波浪引起的水質(zhì)點運動主要集中在近水面的區(qū)域內(nèi),在距水面較遠處波浪對圓柱運動的影響基本可以忽略。

    圖7 圓柱在波浪中入水各點橫向速度對比Fig.7 The comparison of the horizontal velocities of water entry in wave

    圖8 圓柱在波浪中入水各點垂向速度對比Fig.8 The comparison of the vertical velocities of water entry in wave

    2.4 不同波浪相位下圓柱入水運動特性分析

    對于無限水深平面進行波,波浪在同一豎直線上的水質(zhì)點做同相位的圓軌線運動,處于波峰A點、波節(jié)B點和波谷C點的水質(zhì)點的速度如圖9所示,調(diào)整波面方程中相位角ε可使圓柱處于波峰、波節(jié)和波谷三個入水位置,計算不同相位下,圓柱體入水過程運動特性。

    圖9 不同入水位置示意圖 Fig.9 The sketch of different entry location

    圖10 俯仰角及偏航角定義Fig.10 The definition of the pitch angle and the yaw angle

    本文對結(jié)構(gòu)入水姿態(tài)的描述采用與出水問題相近的偏航角與俯仰角,其定義如圖10所示,偏航角定義為結(jié)構(gòu)縱軸與大地坐標(biāo)系X軸正向之間的夾角θ,偏航角定義為結(jié)構(gòu)質(zhì)量中心與初始點之間連線與Z軸負向之間夾角α,結(jié)構(gòu)向X軸正向偏轉(zhuǎn)時俯仰角和偏航角為正向。

    由圖11可以看出,由于波浪波峰與波谷處水質(zhì)點運動速度相反,圓柱體在波峰與波谷入水時產(chǎn)生相反的運動軌跡,相同的入水速度及計算時間范圍內(nèi),波谷入水時,圓柱體橫向運動距離大于波峰入水情況,波節(jié)處入水情況下圓柱體橫向運動距離最小,其運動規(guī)律與波谷入水情況相似。圖12給出不同相位處圓柱體入水過程橫向速度變化,其決定質(zhì)心橫向運動的距離,三個入水相位中,波谷入水時橫向速度峰值最大。圖12中第一個速度峰值時刻對應(yīng)圓柱體頭部完全入水時刻,第二個速度峰值時刻對應(yīng)圓柱質(zhì)心垂向運動位移等于半個波幅值時,當(dāng)圓柱體質(zhì)心垂向運動位移等于波幅后,波浪影響逐步較小,圓柱體質(zhì)心橫向速度變化平緩。

    圖11 質(zhì)心運動軌跡Fig.11 The motion track of the center of mass

    圖12不同相位處入水質(zhì)心橫向速度比較Fig.12 The comparison of the horizontal velocities of water entry in different phase

    圖13、14給出圓柱在不同相位入水時,圓柱入水及水中運動過程俯仰角與偏航角變化情況,圓柱在三個相位處入水時,均產(chǎn)生明顯的偏航角,在近波面,圓柱體受到波浪的影響較大,俯仰角變化速率較大,隨入水深度增加,波浪的影響逐漸較小,圓柱體在表面壓差的作用下,俯仰角繼續(xù)減小甚至發(fā)生傾覆。偏航角表征圓柱體入水過程結(jié)構(gòu)質(zhì)心偏離初始位置情況,表明入水過程橫向載荷的變化。由圖14可以看出,圓柱體在不同相位處入水偏航角變化不同,波峰入水與波節(jié)、波谷入水產(chǎn)生相反的偏航角,波峰入水時的偏航角最大值大于波峰入水,波節(jié)入水過程中偏航角變化最小。隨入水深度增加,波浪的影響逐漸較小,偏航角逐漸趨于平衡。

    圖13 不同相位處入水俯仰角比較Fig.13 The comparison of the pitch angle of water entry in different phase

    圖14 不同相位處入水偏航角比較Fig.14 The comparison of the yaw angle of water entry in different phase

    2.5 浪級對圓柱入水運動特性影響

    為了研究不同浪級對結(jié)構(gòu)入水影響,選取波高、周期如圖15所示規(guī)則波,圓柱接觸水面時的波浪正處于波谷位置。

    圖16給出圓柱在不同浪級作用下入水過程橫向?qū)Ρ惹€。隨著浪級的增大,圓柱入水過程中質(zhì)心橫向速度明顯增加。

    圖17給出圓柱在不同浪級下俯仰角對比,隨著浪級的增大,圓柱入水過程俯仰角增加,波浪誘導(dǎo)結(jié)構(gòu)產(chǎn)生的俯仰力矩大于橫向載荷,浪級變化作用對結(jié)構(gòu)入水姿態(tài)影響明顯。

    圖15 不同浪級示意圖Fig.15 The sketch of different wave scale

    圖16 不同浪級入水橫向速度對比Fig.16 The comparison of the horizontal velocities of water entry in different wave scale

    圖17 不同浪級入水俯仰角對比Fig.17 The comparison of the pitch angle of water entry in different wave scale

    3 海流作用下結(jié)構(gòu)入水特性

    海流又稱洋流,是海水在地轉(zhuǎn)偏向力、引潮力等作用而產(chǎn)生的相對穩(wěn)定的大規(guī)模穩(wěn)定的流動,因此,可將海流簡化為均勻流。計算時取垂直入水速度為25 m/s,均勻流向右傳播,流速度為2.0 m/s。下圖給出均勻洋流作用下圓柱體入水過程三個節(jié)點橫向速度和垂向速度比較。

    圖18 橫向速度對比Fig.18 The comparison of the horizontal velocities

    圖19 垂向速度對比Fig.19 The comparison of the vertical velocities

    由圖18可以看出,在均勻流作用下圓柱體入水過程也存在三個明顯的拐點,即頭部(節(jié)點5)入水、質(zhì)心(節(jié)點743)和尾部(節(jié)點1)入水三個時刻。在尾部入水時刻,三個特征節(jié)點的橫向速度出現(xiàn)明顯的波動,因為在圓柱體全部入水后,附加質(zhì)量變化率在此時刻忽然變?yōu)榱?。在整個入水過程中,頭部質(zhì)心(節(jié)點5)較尾部質(zhì)心(節(jié)點1)變化劇烈,圓柱質(zhì)心的橫向速度在整個入水過程中逐漸趨于洋流速度2.0 m/s。圓柱體三個特征節(jié)點的垂向速度在尾部入水前并未出現(xiàn)明顯變化,在尾部入水時刻出現(xiàn)波動,之后繼續(xù)下降。在圓柱體全部入水后,在自身重力和浮力作用下,且此時圓柱體偏轉(zhuǎn)角度較大,三個節(jié)點的垂向速度發(fā)生較大的差異。

    圖20 俯仰力矩變化曲線 Fig.20 The history curve of the moment

    圖21 水平合力變化曲線Fig.21 The history curve of the horizontal force

    圖20給出圓柱體入水過程中繞圓柱體質(zhì)心俯仰力矩變化,由圖可以看出,俯仰力矩在尾部入水前為負值,隨著入水深度的增加,圓柱體質(zhì)心上部和下部壓力分布逐步均勻,俯仰力矩并逐步減小。在尾入水時刻,俯仰力矩波動較大。完全入水后,由于圓柱體表面僅受到靜水壓力的影響,上下載荷分布較均勻,俯仰力矩變化較平緩。圖21給出圓柱入水過程中水平方向合力的變化,由圖可以看出,圓柱體入水深度增加,水平合力逐漸增大,在尾部入水時刻,水平合力出現(xiàn)較大波動。在圓柱體完全入水后,由于迎、背流面壓差僅為靜水壓力所造成的壓差,其相差較小,因此,合外力變化較小且趨于零值。

    針對不同流速下圓柱體入水影響,分別選取流速為0.5 m/s、1.0 m/s、1.5 m/s及2.0 m/s,計算不同流速下圓柱體的運動狀態(tài)與載荷特性。圖22給出不同流速下圓柱體入水過程中水平速度對比。

    由圖22可以看出,隨洋流速度增大,圓柱體入水質(zhì)心水平方向速度增大,且逐漸趨于流速,圓柱體完全入水后,質(zhì)心水平速度變化趨于平緩。下圖給出不同流速下圓柱體入水過程中俯仰角和偏航角對比。

    圖22 水平速度對比Fig.22 The comparison of the horizontal velocities

    圖23 俯仰角對比Fig.23 The comparison of the pitch angle

    圖24 偏航角對比Fig.24 The comparison of the yaw angle

    由圖23、24可以看出,隨洋流速度增大,圓柱入水過程的俯仰角和偏航角越大,即俯仰力矩和水平合力越大。在整個入水過程中,俯仰角和偏航角隨流速增加變化越劇烈。

    4 波與流共同作用下結(jié)構(gòu)入水分析

    在實際海況下,波與流往往是同時存在的,在計算實際海況下結(jié)構(gòu)入水問題,要考慮波與流共同對結(jié)構(gòu)物的作用。相對于船舶與海洋結(jié)構(gòu)物的尺度,海洋中海流的尺度較廣,且沿水深變化較小,因而可將波、流共同作用下結(jié)構(gòu)入水簡化為波浪與均勻流共同作用下結(jié)構(gòu)入水問題。

    4.1 波與流速度矢量共線

    分別選取洋流速度與波浪傳播方向同向和逆向兩種形式,波浪浪級取為3級、洋流速度取為1.5 m/s,計算波流共線時圓柱體入水過程運動特性。圖25給出不同流速下圓柱體入水過程中水平速度對比。

    圖25 水平速度對比Fig.25 The comparison of the horizontal velocities

    由圖25可以看出,在波流共同作用下,圓柱體入水過程運動特性與姿態(tài)發(fā)生較大變化。三級海浪的作用效果大于1.5 m/s流速,波流同向時增大了圓柱體質(zhì)心的橫向速度,波流反向時,抑制了圓柱體橫向速度變化,由于波浪的作用大于流速,圓柱體質(zhì)心橫向速度仍為正向。下圖給出不同流速下圓柱體入水過程中俯仰角和偏航角對比。

    圖26 俯仰角對比Fig.26 The comparison of the pitch angle

    圖27 偏航角對比Fig.27 The comparison of the yaw angle

    由圖26、27可以看出,波浪、洋流、波流共同作用四種情況下,圓柱均按順時針方向繞質(zhì)心做俯仰運動,波流同向時俯仰角和偏航角最大,反向時最小。在給定波浪等級和洋流速度情況下,僅有波浪作用時圓柱體的俯仰角、偏航角大于僅有洋流作用時俯仰角、偏航角。

    下圖給出不同流速下圓柱體入水過程中俯仰角速度和偏航角速度對比。

    28俯仰角速度對比Fig.28 The comparison of the pitch angle velocity

    圖29 偏航角速度對比Fig.29 The comparison of the yaw angle velocity

    由圖28、29可以看出,波流的共同作用下,俯仰角速度增加到最大值的時刻提前,在相同的入水時間條件下,波流共同作用下圓柱體俯仰角速度逐漸減小為零,而偏航角速度變化規(guī)律相同。

    4.2 波與流速度矢量相差90°

    實際海洋環(huán)境下,波與流的傳播方向有時并不同向,本節(jié)選取波浪傳播方向沿X軸正向,流傳播方向沿Y正向,計算波與流傳播方向不同時結(jié)構(gòu)入水運動特性,并與僅存在波浪和僅存在洋流作用下圓柱入水特性進行對比。

    圖30給出僅存在波浪、僅存在洋流及波流共同作用下圓柱質(zhì)心位置變化。

    由圖30可以看出,僅存在波浪時,質(zhì)心偏移方向沿X軸正向,僅存在洋流時,質(zhì)心偏移方向沿Y正向,波流共同作用時,質(zhì)心運動偏轉(zhuǎn)方向根據(jù)浪級的大小和流速大小確定。

    圖30 質(zhì)心位置對比Fig.30 The motion track of the center of mass

    5 結(jié) 論

    本文采用非線性雙漸進法對圓柱非定常入水過程進行模擬,計算值與試驗值吻合較好,本方法適用于三維入水問題求解。通過對波、流及波流載荷共同作用下三維剛體結(jié)構(gòu)入水運動、姿態(tài)變化特性得到以下結(jié)論:

    (1)圓柱入水過程中,存在兩個關(guān)鍵時刻,尾部入水瞬間和頭部入水瞬間,對結(jié)構(gòu)運動速度產(chǎn)生明顯的影響。在這兩個時刻,入水砰擊力存在峰值,其本質(zhì)原因是在入水過程中附加質(zhì)量的劇烈變化,也就是濕表面的劇烈變化。

    (2)不同波浪相位入水時,圓柱體的運動軌跡與姿態(tài)差異較大,波峰入水與波節(jié)、波谷入水俯仰與偏航運動相反;波浪浪級增大加劇圓柱體的運動的偏移。

    (3)在近波面,波浪對結(jié)構(gòu)作用效果明顯,結(jié)構(gòu)完全入水后,波浪作用效果逐漸減弱。同時波浪載荷對結(jié)構(gòu)軸向載荷影響較橫向載荷小。

    (4)同流速對結(jié)構(gòu)入水產(chǎn)生明顯影響,流速越大,圓柱體入水速度和姿態(tài)變化越劇烈。

    (5)波流共同作用下,結(jié)構(gòu)入水特性改變,波流同向加劇圓柱入水姿態(tài)變化,波流反向抑制圓柱體入水姿態(tài)變化,波流流向不共線時,圓柱體產(chǎn)生兩方向俯仰與偏航運動。

    [1]殷崇一.潛射導(dǎo)彈發(fā)射與出水載荷研究[D].西安:西北工業(yè)大學(xué),2004.Yin Chongyi.Research on launch and effluent load of submarine-launched missile[D].Xian:Northwestern Polytechnical University,2004.

    [2]Faltinsen O M.Sea loads on ships and offshore structures[M].UK:Cambridge University Press,1990.

    [3]王文華,王言英.圓柱在波浪中入水的數(shù)值模擬[J].上海交通大學(xué)學(xué)報,2010,44(10):1393-1399. Wang Wenhua,Wang Yanying.Numerical study on cylinder entering water in wave[J].Journal of Shanghai Jiao Tong University,2010,44(10):1393-1399.

    [4]儲慧林.魚雷在波浪中入水的數(shù)值模擬[D].鎮(zhèn)江:江蘇科技大學(xué),2014.Chu Huilin.Numerical simulation of torpedo entering water in wave[D].Zhenjiang:Jiangsu University of Science and Technology,2014.

    [5]袁緒龍,朱 珠.預(yù)置舵角對高速入水彈道和流體動力的影響[J].應(yīng)用力學(xué)學(xué)報,2015,33(1):11-16. Yuan Xulong,Zhu Zhu.Influence of preset rudder angle on trajectory and hydrodynamic at high-speed water-entry[J]. Chinese Journal of Applied Mechanics,2015,33(1):11-16.

    [6]權(quán)曉波,孔德才,李 巖.波浪模擬及其對水下航行體出水過程影響[J].哈爾濱工業(yè)大學(xué)學(xué)報,2011,43(3):140-144. Qüan Xiaobo,Kong Decai,Li Yan.Wave simulation and its effects on the exceeding water process of the underwater vehicle[J].Journal of Harbin Institute of Technology,2011,43(3):140-144.

    [7]Wang C Z,Wu G X.Interactions between fully nonlinear water waves and cylinder arrays in a wave tank[J].Ocean Engineering,2010,37:400-417.

    [8]劉云龍,張阿漫,王詩平,田昭麗.基于邊界元法的氣泡同波浪相互作用研究[J].物理學(xué)報,2012,61(22):1-8. Liu Yunlong,Zhang Aman,Wang Shiping,Tian Zhaoli.Research on interaction between bubble and surface waves based on BEM[J].Acta Phys.Sin.,2012,61(22):1-8.

    [9]王詩平,孫士麗,張阿漫,陳 玉.可壓縮流場中氣泡脈動數(shù)值模擬[J].力學(xué)學(xué)報,2012,44(3):513-520. Wang Shiping,Sun Shili,Zhang Aman,Chen Yu.Numercal simulation of bubble dynamics in compressible fluid[J].Acta Mechanica Sinica,2012,44(3):513-520.

    [10]王文華,黃 一,王言英,翟鋼軍,黃亞南.彈性楔形體各狀態(tài)參數(shù)對入水運動性能的影響[J].船舶力學(xué),2014,18(11): 1320-1330.Wang Wenhua,Huang Yi,Wang Yanying,Zhai Gangjun,Huang Yanan.Effect of status parameters for elastic wedge on dynamic performance of water-entry[J].Journal of Ship Mechanics,2014,18(11):1320-1330.

    [11]Geers T L.Residual potential and approximation methods for three dimensional fluid-structure interaction problems[J].J Acoust.Soc.Amer.,1971,49:1505-1510.

    [12]Geers T L.Doubly asympotic approximation for transient motions of submerged structures[J].J Acoust.Soc.Ame.,1978, 64:1500-1508.

    [13]Geers T L,Felippa C A.Doubly asymptotic approximations for vibration analysis of submerged structures[J].J Acoust.Soc. Amer.,1980,73:1152-1159.

    [14]孫士麗.具有自由液面效應(yīng)的高頻流固耦合問題[D].哈爾濱:哈爾濱工程大學(xué),2010. Sun Shili.The high frequency fluid structure interaction problems with free surface effect[D].Harbin:Harbin Engineering University,2010.

    3D numerical simulation on water entry of cylindrical under wave and stream action

    YANG Heng,SUN Long-quan,LIU Ying,YAO Xiong-liang
    (Harbin Engineering University,Haerbin 150001,China)

    During water entry of the structure in actual sea states,it will receive the interaction of a variety of loads,as well as the influence of wave and stream action,so this is a strong non-linear process.According to the characteristics of water entry of the structure,the incident wave and stream will be introduced to the nonlinear doubly asymptotic approximation method,so as to study the motion response and attitude change of 3D rigid cylindrical under wave and stream action during water entry.The calculation results fit with the test results well,this method is available for the analysis of water entry of 3D rigid body.The results show that in the vicinity of wave surface,the wave affects the structure distinctly;phases of water entry,wave scale,stream velocity and the differences of wave and stream velocities vector influence the speed, motion trajectory of the structure in water notably.

    water entry;wave and stream;non-linear;DAA(doubly asymptotic approximation)

    U671.5

    A

    10.3969/j.issn.1007-7294.2015.10.003

    1007-7294(2015)10-1186-11

    2015-05-03

    船舶預(yù)研支撐技術(shù)基金資助項目(13J1.5.1);國家自然科學(xué)基金資助項目(51309060)

    楊 衡(1986-),男,博士研究生,E-mail:yanghengmvp@163.com;孫龍泉(1983-),男,講師。

    猜你喜歡
    波流圓柱體質(zhì)心
    波流耦合下樁周珊瑚砂沖刷機理研究
    海洋通報(2022年5期)2022-11-30 12:04:02
    重型半掛汽車質(zhì)量與質(zhì)心位置估計
    基于GNSS測量的天宮二號質(zhì)心確定
    波流聯(lián)合作用下海上輸油漂浮軟管動力響應(yīng)分析
    偶感
    槽道內(nèi)渦波流場展向渦的分布特征
    找出圓柱體
    娃娃畫報(2016年6期)2016-07-29 19:31:34
    圓柱體上的最短路徑
    一種海洋測高衛(wèi)星質(zhì)心在軌估計算法
    航天器工程(2014年5期)2014-03-11 16:35:53
    觀察力
    少妇被粗大猛烈的视频| 王馨瑶露胸无遮挡在线观看| 精品午夜福利在线看| 亚洲伊人久久精品综合| 超色免费av| 成人毛片60女人毛片免费| 在线精品无人区一区二区三| 男女啪啪激烈高潮av片| 男女国产视频网站| 亚洲精品日本国产第一区| 国产精品蜜桃在线观看| 99re6热这里在线精品视频| 亚洲国产日韩一区二区| 亚洲精品久久久久久婷婷小说| 两性夫妻黄色片| 99久久中文字幕三级久久日本| 久久99热这里只频精品6学生| 久久人人爽人人片av| 精品亚洲成国产av| 大话2 男鬼变身卡| 亚洲内射少妇av| 国语对白做爰xxxⅹ性视频网站| 纯流量卡能插随身wifi吗| 国产毛片在线视频| 如日韩欧美国产精品一区二区三区| 国产片内射在线| 亚洲国产色片| 国产又色又爽无遮挡免| 国产极品天堂在线| 国产高清国产精品国产三级| 免费黄色在线免费观看| 18+在线观看网站| 边亲边吃奶的免费视频| 亚洲欧洲精品一区二区精品久久久 | 日韩一卡2卡3卡4卡2021年| 黄色毛片三级朝国网站| 伊人久久国产一区二区| 男女免费视频国产| 亚洲av福利一区| 国产精品久久久久成人av| 新久久久久国产一级毛片| 久久女婷五月综合色啪小说| 看非洲黑人一级黄片| 精品久久久久久电影网| 成人国语在线视频| 色哟哟·www| 卡戴珊不雅视频在线播放| 国产成人免费无遮挡视频| 国产成人免费无遮挡视频| 中文字幕av电影在线播放| 久久久久久久久免费视频了| 欧美日韩亚洲国产一区二区在线观看 | 久久综合国产亚洲精品| 最近的中文字幕免费完整| 欧美人与性动交α欧美软件| 亚洲一区二区三区欧美精品| freevideosex欧美| 两性夫妻黄色片| 亚洲av成人精品一二三区| 亚洲综合色网址| www.精华液| 人人澡人人妻人| 亚洲欧美一区二区三区久久| 欧美精品亚洲一区二区| 欧美日韩av久久| 在线亚洲精品国产二区图片欧美| 午夜影院在线不卡| 国产午夜精品一二区理论片| 巨乳人妻的诱惑在线观看| 国产精品香港三级国产av潘金莲 | 中国三级夫妇交换| 亚洲欧洲国产日韩| 精品久久久久久电影网| 日本wwww免费看| 中文精品一卡2卡3卡4更新| 欧美少妇被猛烈插入视频| 成人漫画全彩无遮挡| 成人国语在线视频| 多毛熟女@视频| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲国产av影院在线观看| 大香蕉久久成人网| 最新的欧美精品一区二区| 欧美日韩精品成人综合77777| 午夜福利一区二区在线看| 亚洲视频免费观看视频| 18+在线观看网站| 久久精品国产鲁丝片午夜精品| 国产一区有黄有色的免费视频| 宅男免费午夜| 国产伦理片在线播放av一区| 国产亚洲精品第一综合不卡| 晚上一个人看的免费电影| 亚洲精品中文字幕在线视频| 午夜精品国产一区二区电影| 久久久久久人妻| 精品亚洲成国产av| 欧美黄色片欧美黄色片| 亚洲精品久久成人aⅴ小说| 涩涩av久久男人的天堂| 欧美亚洲日本最大视频资源| 免费黄网站久久成人精品| www.自偷自拍.com| 三上悠亚av全集在线观看| 亚洲av综合色区一区| 午夜老司机福利剧场| 国产片特级美女逼逼视频| 免费日韩欧美在线观看| 午夜福利影视在线免费观看| 亚洲欧美一区二区三区国产| 中文字幕精品免费在线观看视频| 色播在线永久视频| 欧美bdsm另类| 久久99蜜桃精品久久| 人人妻人人澡人人爽人人夜夜| 久久久久精品人妻al黑| 国产在线视频一区二区| 亚洲,欧美,日韩| 亚洲久久久国产精品| 亚洲欧美成人精品一区二区| 日本黄色日本黄色录像| 久久精品人人爽人人爽视色| 久久99精品国语久久久| 日韩三级伦理在线观看| 青青草视频在线视频观看| 免费黄色在线免费观看| 在线观看免费视频网站a站| 成人亚洲精品一区在线观看| av网站免费在线观看视频| 国产精品一国产av| 人人妻人人澡人人看| 日本wwww免费看| 9色porny在线观看| 又黄又粗又硬又大视频| 丰满少妇做爰视频| 亚洲欧美精品综合一区二区三区 | 在线 av 中文字幕| 国产1区2区3区精品| 国产精品免费视频内射| 韩国高清视频一区二区三区| 新久久久久国产一级毛片| 十八禁高潮呻吟视频| 18禁国产床啪视频网站| 91精品三级在线观看| 欧美日韩综合久久久久久| 国产成人精品久久久久久| 美女脱内裤让男人舔精品视频| 国产男女超爽视频在线观看| 成人免费观看视频高清| 亚洲精品美女久久久久99蜜臀 | 99热网站在线观看| 亚洲一级一片aⅴ在线观看| 搡女人真爽免费视频火全软件| 少妇被粗大猛烈的视频| av网站在线播放免费| 欧美精品国产亚洲| 欧美另类一区| 日韩免费高清中文字幕av| www日本在线高清视频| 777久久人妻少妇嫩草av网站| 少妇人妻精品综合一区二区| 国产精品国产三级国产专区5o| 国产日韩欧美视频二区| 老汉色∧v一级毛片| 亚洲精品第二区| 天天躁夜夜躁狠狠躁躁| 欧美成人午夜免费资源| 一级爰片在线观看| 国产一区二区三区综合在线观看| 久久毛片免费看一区二区三区| 亚洲精品一二三| 看免费成人av毛片| 不卡视频在线观看欧美| 久久久久久久久久久免费av| 欧美xxⅹ黑人| 麻豆av在线久日| av免费在线看不卡| 少妇人妻久久综合中文| 一本—道久久a久久精品蜜桃钙片| 亚洲欧洲精品一区二区精品久久久 | 国产精品无大码| 亚洲,欧美,日韩| 黄片小视频在线播放| 夫妻性生交免费视频一级片| 午夜日本视频在线| 亚洲熟女精品中文字幕| av片东京热男人的天堂| 2018国产大陆天天弄谢| 看免费av毛片| 欧美日韩国产mv在线观看视频| 国产精品嫩草影院av在线观看| 侵犯人妻中文字幕一二三四区| 999精品在线视频| 亚洲图色成人| 国产xxxxx性猛交| 久久亚洲国产成人精品v| 久久av网站| 人人妻人人澡人人爽人人夜夜| 一级,二级,三级黄色视频| 成人国产av品久久久| 丰满饥渴人妻一区二区三| 男女午夜视频在线观看| 1024香蕉在线观看| 国产av码专区亚洲av| 欧美日韩视频精品一区| 少妇人妻 视频| 欧美另类一区| 亚洲精品国产色婷婷电影| 在线观看美女被高潮喷水网站| xxxhd国产人妻xxx| 最近2019中文字幕mv第一页| 久久人人97超碰香蕉20202| 九九爱精品视频在线观看| 亚洲欧美色中文字幕在线| 在线精品无人区一区二区三| 夫妻性生交免费视频一级片| 黄频高清免费视频| 哪个播放器可以免费观看大片| 菩萨蛮人人尽说江南好唐韦庄| 日韩大片免费观看网站| 国产精品女同一区二区软件| 春色校园在线视频观看| 午夜精品国产一区二区电影| 久久久欧美国产精品| 亚洲激情五月婷婷啪啪| 一区二区av电影网| 91在线精品国自产拍蜜月| 国产1区2区3区精品| 又黄又粗又硬又大视频| 黄网站色视频无遮挡免费观看| 狠狠精品人妻久久久久久综合| 欧美97在线视频| 成年人免费黄色播放视频| 免费高清在线观看视频在线观看| 成人国产av品久久久| 午夜影院在线不卡| 久久亚洲国产成人精品v| 国产精品国产三级专区第一集| 一级毛片 在线播放| 91久久精品国产一区二区三区| 一级毛片我不卡| 狠狠精品人妻久久久久久综合| 免费在线观看完整版高清| 久久久国产精品麻豆| 一级,二级,三级黄色视频| 2018国产大陆天天弄谢| 久久99热这里只频精品6学生| 亚洲精品第二区| 丝袜喷水一区| 伦精品一区二区三区| 大香蕉久久网| 日韩人妻精品一区2区三区| 中文字幕色久视频| videossex国产| 久久精品国产亚洲av天美| 99re6热这里在线精品视频| 一级片免费观看大全| 色哟哟·www| 99久国产av精品国产电影| 老汉色∧v一级毛片| 久热这里只有精品99| 香蕉丝袜av| 性少妇av在线| 日日爽夜夜爽网站| 侵犯人妻中文字幕一二三四区| 老女人水多毛片| 大香蕉久久网| 五月伊人婷婷丁香| 久久ye,这里只有精品| 免费日韩欧美在线观看| 丰满饥渴人妻一区二区三| 久久久久久久亚洲中文字幕| 亚洲男人天堂网一区| 亚洲av中文av极速乱| 青春草亚洲视频在线观看| 精品久久久精品久久久| av线在线观看网站| 在线精品无人区一区二区三| 老鸭窝网址在线观看| 久久久久国产一级毛片高清牌| 国产精品久久久久久av不卡| 婷婷成人精品国产| 国产日韩欧美在线精品| 国产白丝娇喘喷水9色精品| 两性夫妻黄色片| 国产人伦9x9x在线观看 | 91精品三级在线观看| 最新的欧美精品一区二区| 国产成人精品婷婷| 成年人午夜在线观看视频| 免费少妇av软件| 国产女主播在线喷水免费视频网站| xxx大片免费视频| 老司机亚洲免费影院| 精品亚洲成a人片在线观看| 国产男女超爽视频在线观看| 日韩不卡一区二区三区视频在线| 久久午夜福利片| 国产无遮挡羞羞视频在线观看| 不卡视频在线观看欧美| 26uuu在线亚洲综合色| 9191精品国产免费久久| 热99久久久久精品小说推荐| 91国产中文字幕| 亚洲婷婷狠狠爱综合网| 999久久久国产精品视频| 日本-黄色视频高清免费观看| 欧美日韩精品成人综合77777| 青青草视频在线视频观看| 亚洲国产成人一精品久久久| 亚洲国产日韩一区二区| 久久精品国产自在天天线| 一区二区三区乱码不卡18| 老熟女久久久| 国产在线免费精品| 国产成人91sexporn| 亚洲四区av| 久久毛片免费看一区二区三区| av网站免费在线观看视频| 国产综合精华液| 久久久久久久大尺度免费视频| 99国产精品免费福利视频| 亚洲av国产av综合av卡| av卡一久久| 女人高潮潮喷娇喘18禁视频| 伊人亚洲综合成人网| 久久久久久久精品精品| 一级黄片播放器| 欧美成人精品欧美一级黄| 九草在线视频观看| 亚洲激情五月婷婷啪啪| 亚洲av.av天堂| 成年人午夜在线观看视频| 岛国毛片在线播放| 亚洲一区二区三区欧美精品| 自拍欧美九色日韩亚洲蝌蚪91| 国产成人a∨麻豆精品| 亚洲精品国产av蜜桃| 久久99一区二区三区| xxx大片免费视频| 欧美成人午夜精品| 国产成人精品福利久久| 七月丁香在线播放| 久久 成人 亚洲| 久久精品熟女亚洲av麻豆精品| 亚洲人成电影观看| 啦啦啦中文免费视频观看日本| 自线自在国产av| a级毛片在线看网站| 激情五月婷婷亚洲| 国产成人欧美| 精品国产乱码久久久久久男人| 国产成人免费观看mmmm| 老司机影院成人| 欧美激情 高清一区二区三区| 欧美bdsm另类| 99九九在线精品视频| 免费黄频网站在线观看国产| 纵有疾风起免费观看全集完整版| 欧美国产精品一级二级三级| 最近的中文字幕免费完整| 亚洲视频免费观看视频| 免费人妻精品一区二区三区视频| 久久人人爽av亚洲精品天堂| 国产深夜福利视频在线观看| 欧美xxⅹ黑人| 国产男人的电影天堂91| 亚洲三区欧美一区| 亚洲,欧美,日韩| av福利片在线| 男男h啪啪无遮挡| 女人高潮潮喷娇喘18禁视频| 国产精品 国内视频| 精品第一国产精品| 精品99又大又爽又粗少妇毛片| 欧美激情高清一区二区三区 | 最新的欧美精品一区二区| 在线观看美女被高潮喷水网站| 久热久热在线精品观看| 亚洲国产日韩一区二区| 日本午夜av视频| 中文字幕精品免费在线观看视频| 亚洲第一青青草原| 精品午夜福利在线看| 久久精品国产a三级三级三级| 黄片无遮挡物在线观看| 久久久久网色| 成人亚洲欧美一区二区av| 精品少妇久久久久久888优播| 一级爰片在线观看| 欧美bdsm另类| 最近2019中文字幕mv第一页| 久久女婷五月综合色啪小说| 制服人妻中文乱码| tube8黄色片| 久久99一区二区三区| 国产日韩一区二区三区精品不卡| 97精品久久久久久久久久精品| 美女脱内裤让男人舔精品视频| 国产日韩一区二区三区精品不卡| 国产一区有黄有色的免费视频| 久久久久精品久久久久真实原创| 一本久久精品| 久久久久久久久久久免费av| 最近手机中文字幕大全| 一级片免费观看大全| 亚洲精品成人av观看孕妇| 在线观看免费日韩欧美大片| 狠狠精品人妻久久久久久综合| 欧美人与性动交α欧美精品济南到 | 亚洲精品国产色婷婷电影| 亚洲综合色网址| 丝袜在线中文字幕| 纵有疾风起免费观看全集完整版| av在线播放精品| 亚洲欧美成人精品一区二区| 黑人猛操日本美女一级片| 精品一区二区三卡| 日本vs欧美在线观看视频| 飞空精品影院首页| freevideosex欧美| 9色porny在线观看| 丝袜脚勾引网站| 国产午夜精品一二区理论片| 美女国产视频在线观看| 国产精品国产三级国产专区5o| 欧美激情高清一区二区三区 | 多毛熟女@视频| 青春草国产在线视频| 在线天堂最新版资源| 在线观看美女被高潮喷水网站| 美女午夜性视频免费| 中文欧美无线码| 亚洲欧美一区二区三区国产| 亚洲成av片中文字幕在线观看 | 九色亚洲精品在线播放| 人成视频在线观看免费观看| 亚洲av免费高清在线观看| 人人妻人人澡人人爽人人夜夜| av在线播放精品| 不卡视频在线观看欧美| 精品第一国产精品| 欧美在线黄色| xxx大片免费视频| 伊人久久大香线蕉亚洲五| 欧美日韩综合久久久久久| 一级毛片我不卡| 日日爽夜夜爽网站| 国产av一区二区精品久久| 激情五月婷婷亚洲| 久久久欧美国产精品| 久久午夜综合久久蜜桃| 亚洲国产精品999| 色视频在线一区二区三区| 欧美黄色片欧美黄色片| 久久综合国产亚洲精品| 国产又色又爽无遮挡免| 欧美精品国产亚洲| 亚洲欧美色中文字幕在线| 大码成人一级视频| 免费少妇av软件| 亚洲内射少妇av| 少妇的逼水好多| 国产人伦9x9x在线观看 | 久久久久视频综合| 国产日韩欧美亚洲二区| 国产成人精品久久二区二区91 | 国产精品偷伦视频观看了| 最近2019中文字幕mv第一页| 亚洲精品自拍成人| 亚洲欧美精品综合一区二区三区 | 亚洲精品国产av蜜桃| 亚洲国产看品久久| 国产欧美日韩综合在线一区二区| 一二三四在线观看免费中文在| 永久网站在线| 99九九在线精品视频| 日韩熟女老妇一区二区性免费视频| 校园人妻丝袜中文字幕| 在线观看免费视频网站a站| 可以免费在线观看a视频的电影网站 | 欧美激情极品国产一区二区三区| 中文欧美无线码| 国产精品免费视频内射| 最近最新中文字幕大全免费视频 | 大话2 男鬼变身卡| 激情视频va一区二区三区| 国产一区二区在线观看av| 国产xxxxx性猛交| 久久久久久久久免费视频了| 国产av码专区亚洲av| 狂野欧美激情性bbbbbb| 国产欧美日韩综合在线一区二区| 亚洲国产毛片av蜜桃av| 在线观看国产h片| 久久精品人人爽人人爽视色| 欧美日韩精品网址| 欧美日韩视频精品一区| 欧美精品高潮呻吟av久久| 激情五月婷婷亚洲| 色94色欧美一区二区| 9热在线视频观看99| 欧美人与性动交α欧美精品济南到 | 大香蕉久久网| 曰老女人黄片| av免费观看日本| 久久久久久久亚洲中文字幕| 人人澡人人妻人| 熟女电影av网| 久久久精品区二区三区| 中国三级夫妇交换| 欧美黄色片欧美黄色片| 日韩免费高清中文字幕av| 男女边摸边吃奶| 国产精品国产三级国产专区5o| 亚洲美女黄色视频免费看| 美女视频免费永久观看网站| 亚洲精品一二三| 免费人妻精品一区二区三区视频| 日韩 亚洲 欧美在线| 一级片免费观看大全| 亚洲欧美中文字幕日韩二区| 亚洲av欧美aⅴ国产| 日韩免费高清中文字幕av| 欧美日韩一区二区视频在线观看视频在线| 高清欧美精品videossex| 日本-黄色视频高清免费观看| 久久精品国产综合久久久| 不卡视频在线观看欧美| 欧美日韩视频高清一区二区三区二| 欧美激情 高清一区二区三区| 成人18禁高潮啪啪吃奶动态图| 中文字幕人妻丝袜一区二区 | 热re99久久精品国产66热6| 成人国语在线视频| 国产成人一区二区在线| 天天躁日日躁夜夜躁夜夜| 日韩av不卡免费在线播放| 亚洲图色成人| 制服人妻中文乱码| 中文字幕另类日韩欧美亚洲嫩草| 国产无遮挡羞羞视频在线观看| av免费在线看不卡| 亚洲久久久国产精品| 在现免费观看毛片| 欧美日韩一级在线毛片| 久久精品aⅴ一区二区三区四区 | 亚洲欧美成人精品一区二区| 午夜福利影视在线免费观看| 国产精品偷伦视频观看了| 少妇精品久久久久久久| 看免费av毛片| 男女免费视频国产| 欧美成人精品欧美一级黄| 亚洲精品国产一区二区精华液| 五月开心婷婷网| 国产黄色视频一区二区在线观看| 日韩一卡2卡3卡4卡2021年| 黄片小视频在线播放| 免费在线观看视频国产中文字幕亚洲 | 超碰成人久久| av福利片在线| 少妇被粗大的猛进出69影院| 久久精品国产综合久久久| 国产黄频视频在线观看| 国产xxxxx性猛交| 高清不卡的av网站| 亚洲精品视频女| av在线app专区| 我的亚洲天堂| 男女免费视频国产| 自线自在国产av| 国产精品一区二区在线不卡| 免费观看在线日韩| 国产精品蜜桃在线观看| 国产亚洲午夜精品一区二区久久| 国产伦理片在线播放av一区| xxxhd国产人妻xxx| 久久久久久久大尺度免费视频| 国产精品av久久久久免费| 日韩免费高清中文字幕av| 满18在线观看网站| av女优亚洲男人天堂| 亚洲精品久久久久久婷婷小说| 国产精品久久久久久久久免| 精品久久久久久电影网| 亚洲三级黄色毛片| 成人18禁高潮啪啪吃奶动态图| 亚洲国产日韩一区二区| 午夜福利影视在线免费观看| 日产精品乱码卡一卡2卡三| 欧美亚洲日本最大视频资源| 视频区图区小说| 国产午夜精品一二区理论片| 又黄又粗又硬又大视频| 精品国产乱码久久久久久小说| 人妻 亚洲 视频| 亚洲欧美清纯卡通| 免费在线观看黄色视频的| 国产成人午夜福利电影在线观看| 久久国内精品自在自线图片| 老司机影院成人| 免费播放大片免费观看视频在线观看| 日韩在线高清观看一区二区三区| 热99国产精品久久久久久7| 欧美bdsm另类| 久久亚洲国产成人精品v| 咕卡用的链子| 亚洲精品视频女| 男人舔女人的私密视频| av网站免费在线观看视频| 男女边吃奶边做爰视频| 亚洲视频免费观看视频| 日本爱情动作片www.在线观看| 色网站视频免费|