• <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
    觀察力
    精品一区二区三卡| 亚洲欧美一区二区三区黑人| 曰老女人黄片| 久9热在线精品视频| 久久毛片免费看一区二区三区| 真人做人爱边吃奶动态| 欧美精品一区二区大全| 黑人巨大精品欧美一区二区蜜桃| 精品国产乱子伦一区二区三区 | 午夜福利在线免费观看网站| 亚洲精品成人av观看孕妇| 一个人免费在线观看的高清视频 | 欧美日韩黄片免| 美女扒开内裤让男人捅视频| 狠狠狠狠99中文字幕| 亚洲国产毛片av蜜桃av| 国产精品国产三级国产专区5o| 精品一区在线观看国产| 国产色视频综合| 大陆偷拍与自拍| 欧美变态另类bdsm刘玥| 亚洲精品第二区| 咕卡用的链子| 91大片在线观看| 黄色片一级片一级黄色片| 在线观看免费日韩欧美大片| 美女高潮喷水抽搐中文字幕| av天堂在线播放| 一二三四在线观看免费中文在| 人人澡人人妻人| 免费不卡黄色视频| 日韩大码丰满熟妇| 女性生殖器流出的白浆| 久久狼人影院| 老司机靠b影院| 色老头精品视频在线观看| 久久国产亚洲av麻豆专区| 法律面前人人平等表现在哪些方面 | 男女无遮挡免费网站观看| 每晚都被弄得嗷嗷叫到高潮| 欧美日韩av久久| 真人做人爱边吃奶动态| av国产精品久久久久影院| 女警被强在线播放| 他把我摸到了高潮在线观看 | 91大片在线观看| 99热网站在线观看| 俄罗斯特黄特色一大片| 美女高潮喷水抽搐中文字幕| 一个人免费看片子| 秋霞在线观看毛片| 我要看黄色一级片免费的| 国产日韩一区二区三区精品不卡| 波多野结衣一区麻豆| 成在线人永久免费视频| 国产精品二区激情视频| 一区福利在线观看| av在线app专区| 黄网站色视频无遮挡免费观看| 91精品国产国语对白视频| 国产日韩一区二区三区精品不卡| 中文字幕人妻熟女乱码| 国产欧美日韩精品亚洲av| 午夜福利视频在线观看免费| 午夜激情av网站| 久久青草综合色| 老司机福利观看| 在线看a的网站| 亚洲第一青青草原| 国产欧美日韩一区二区三区在线| 在线 av 中文字幕| 正在播放国产对白刺激| 男人操女人黄网站| 国产有黄有色有爽视频| 久久av网站| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品1区2区在线观看. | 成人三级做爰电影| 大型av网站在线播放| 精品乱码久久久久久99久播| 嫩草影视91久久| 久久久久国产精品人妻一区二区| 亚洲欧美一区二区三区久久| 人人澡人人妻人| 大码成人一级视频| 国产av又大| 国内毛片毛片毛片毛片毛片| 成年动漫av网址| 女人高潮潮喷娇喘18禁视频| 亚洲精品一卡2卡三卡4卡5卡 | 久久精品亚洲熟妇少妇任你| 9色porny在线观看| 欧美日韩亚洲国产一区二区在线观看 | 亚洲第一青青草原| 精品人妻一区二区三区麻豆| 狠狠精品人妻久久久久久综合| 亚洲五月色婷婷综合| 久久女婷五月综合色啪小说| 久久精品aⅴ一区二区三区四区| 欧美日韩黄片免| 亚洲欧美精品综合一区二区三区| 日韩免费高清中文字幕av| 真人做人爱边吃奶动态| 青青草视频在线视频观看| 可以免费在线观看a视频的电影网站| 极品人妻少妇av视频| 国产亚洲精品久久久久5区| 真人做人爱边吃奶动态| 女人久久www免费人成看片| 人妻人人澡人人爽人人| 久久久精品94久久精品| 国产精品久久久久久人妻精品电影 | 国产亚洲欧美精品永久| 精品一区二区三区四区五区乱码| 久9热在线精品视频| 成年人黄色毛片网站| 成年美女黄网站色视频大全免费| 亚洲人成77777在线视频| 欧美日本中文国产一区发布| www.999成人在线观看| 纵有疾风起免费观看全集完整版| 亚洲九九香蕉| 国产一区二区激情短视频 | 国产一区二区三区av在线| 中文字幕另类日韩欧美亚洲嫩草| www.999成人在线观看| 亚洲成人国产一区在线观看| 亚洲少妇的诱惑av| 99久久人妻综合| 国产无遮挡羞羞视频在线观看| 欧美日韩精品网址| 国产成人精品久久二区二区免费| 精品高清国产在线一区| 亚洲情色 制服丝袜| 国产亚洲欧美精品永久| 国产欧美日韩综合在线一区二区| 亚洲人成电影观看| 久久精品亚洲av国产电影网| 美女大奶头黄色视频| 久久人妻福利社区极品人妻图片| 黄色视频,在线免费观看| 亚洲精品美女久久av网站| 一二三四社区在线视频社区8| 人妻一区二区av| 1024香蕉在线观看| 亚洲精品一二三| 日韩电影二区| 精品国产一区二区三区四区第35| 午夜免费成人在线视频| 99国产极品粉嫩在线观看| 国产深夜福利视频在线观看| 丝袜在线中文字幕| 少妇粗大呻吟视频| 91字幕亚洲| 国产精品久久久久久精品古装| 亚洲五月色婷婷综合| 天堂俺去俺来也www色官网| 国产日韩欧美视频二区| 国产精品熟女久久久久浪| 欧美黄色淫秽网站| 国产精品自产拍在线观看55亚洲 | 啦啦啦啦在线视频资源| 成年女人毛片免费观看观看9 | 午夜福利一区二区在线看| 亚洲成国产人片在线观看| 91麻豆精品激情在线观看国产 | 法律面前人人平等表现在哪些方面 | 最近最新免费中文字幕在线| 国产一区二区三区在线臀色熟女 | 99九九在线精品视频| 成年人免费黄色播放视频| 一个人免费看片子| 精品人妻1区二区| 热re99久久精品国产66热6| 久久人人爽人人片av| 91九色精品人成在线观看| 大片电影免费在线观看免费| 一本—道久久a久久精品蜜桃钙片| 波多野结衣一区麻豆| 最近最新中文字幕大全免费视频| 99国产精品一区二区蜜桃av | 欧美日韩福利视频一区二区| 麻豆av在线久日| 亚洲一区二区三区欧美精品| 亚洲av国产av综合av卡| 天天操日日干夜夜撸| 伊人亚洲综合成人网| 天天操日日干夜夜撸| 黄频高清免费视频| 国产免费福利视频在线观看| 久久中文字幕一级| 丝袜美腿诱惑在线| 纵有疾风起免费观看全集完整版| 狠狠狠狠99中文字幕| 香蕉丝袜av| 国产伦人伦偷精品视频| 国产精品一区二区在线观看99| 五月开心婷婷网| av网站免费在线观看视频| 国产福利在线免费观看视频| 老司机影院毛片| av有码第一页| 国产欧美亚洲国产| 亚洲精品一二三| 亚洲精品美女久久久久99蜜臀| 在线精品无人区一区二区三| 在线永久观看黄色视频| a级片在线免费高清观看视频| 一二三四社区在线视频社区8| 老汉色av国产亚洲站长工具| 成年人免费黄色播放视频| 成在线人永久免费视频| 黄频高清免费视频| 亚洲激情五月婷婷啪啪| 成年人黄色毛片网站| bbb黄色大片| 国产男人的电影天堂91| 久久热在线av| 亚洲 欧美一区二区三区| 91九色精品人成在线观看| 国产精品久久久久成人av| 欧美日韩国产mv在线观看视频| 热re99久久国产66热| 美国免费a级毛片| 色老头精品视频在线观看| 韩国高清视频一区二区三区| 久久久久久久国产电影| 在线观看人妻少妇| av视频免费观看在线观看| 99re6热这里在线精品视频| 一个人免费在线观看的高清视频 | 精品免费久久久久久久清纯 | 精品国内亚洲2022精品成人 | 激情视频va一区二区三区| 天堂俺去俺来也www色官网| 久久久久久久大尺度免费视频| 丝瓜视频免费看黄片| 精品人妻一区二区三区麻豆| 久久亚洲精品不卡| 日本精品一区二区三区蜜桃| 亚洲av日韩精品久久久久久密| 狠狠狠狠99中文字幕| 后天国语完整版免费观看| 色婷婷av一区二区三区视频| 国产精品一区二区在线不卡| 啦啦啦视频在线资源免费观看| 人妻人人澡人人爽人人| 在线观看www视频免费| tube8黄色片| 国产主播在线观看一区二区| 国产精品二区激情视频| 欧美激情高清一区二区三区| 视频在线观看一区二区三区| 日本撒尿小便嘘嘘汇集6| 久久影院123| 欧美乱码精品一区二区三区| 久久人人爽人人片av| 女人高潮潮喷娇喘18禁视频| 国产精品国产av在线观看| 大香蕉久久网| 欧美在线一区亚洲| 别揉我奶头~嗯~啊~动态视频 | 亚洲精品中文字幕在线视频| 涩涩av久久男人的天堂| 人妻一区二区av| 国产福利在线免费观看视频| 亚洲av国产av综合av卡| 俄罗斯特黄特色一大片| 大型av网站在线播放| 亚洲人成电影观看| 爱豆传媒免费全集在线观看| 日韩大片免费观看网站| 亚洲精品国产色婷婷电影| 狂野欧美激情性xxxx| 久久九九热精品免费| 午夜福利免费观看在线| 久久精品国产亚洲av香蕉五月 | 一本大道久久a久久精品| 搡老乐熟女国产| 一本—道久久a久久精品蜜桃钙片| 欧美成狂野欧美在线观看| 亚洲精品中文字幕一二三四区 | 欧美乱码精品一区二区三区| 亚洲精品国产一区二区精华液| 69精品国产乱码久久久| 亚洲黑人精品在线| 99精国产麻豆久久婷婷| 国产有黄有色有爽视频| 欧美精品人与动牲交sv欧美| 久久这里只有精品19| 91精品国产国语对白视频| 亚洲一区中文字幕在线| 亚洲精品在线美女| 啦啦啦免费观看视频1| 最近最新中文字幕大全免费视频| 首页视频小说图片口味搜索| av天堂久久9| 男男h啪啪无遮挡| 黄色视频在线播放观看不卡| 12—13女人毛片做爰片一| 欧美在线黄色| 日本精品一区二区三区蜜桃| 黄频高清免费视频| 国产成人精品无人区| 亚洲精品自拍成人| 汤姆久久久久久久影院中文字幕| 在线观看免费午夜福利视频| 欧美黄色淫秽网站| 亚洲情色 制服丝袜| 欧美国产精品va在线观看不卡| 亚洲国产成人一精品久久久| 成年动漫av网址| 女性生殖器流出的白浆| 伊人亚洲综合成人网| 免费在线观看黄色视频的| 青春草亚洲视频在线观看| 欧美国产精品va在线观看不卡| 亚洲精品国产av蜜桃| 精品一品国产午夜福利视频| 女人久久www免费人成看片| 一级,二级,三级黄色视频| 波多野结衣一区麻豆| 一级片'在线观看视频| 日本av手机在线免费观看| av不卡在线播放| 国产欧美日韩一区二区三区在线| 丝袜喷水一区| 久久国产精品影院| 女警被强在线播放| 69精品国产乱码久久久| 亚洲熟女精品中文字幕| 国产成人免费无遮挡视频| 一区在线观看完整版| 日本wwww免费看| 久久久精品94久久精品| 在线十欧美十亚洲十日本专区| 欧美日韩国产mv在线观看视频| 日韩中文字幕欧美一区二区| 久久人人爽人人片av| 一个人免费看片子| 国产一区二区激情短视频 | 啦啦啦视频在线资源免费观看| 一边摸一边抽搐一进一出视频| www.自偷自拍.com| 91精品三级在线观看| 国产精品成人在线| 丰满迷人的少妇在线观看| 老熟妇乱子伦视频在线观看 | 久久狼人影院| 老司机影院毛片| 亚洲伊人久久精品综合| 国产成人av激情在线播放| 秋霞在线观看毛片| 国产精品国产av在线观看| 交换朋友夫妻互换小说| 操美女的视频在线观看| 黄色毛片三级朝国网站| 韩国高清视频一区二区三区| 黄片小视频在线播放| 十八禁人妻一区二区| 中国国产av一级| 国产亚洲一区二区精品| 精品乱码久久久久久99久播| 91国产中文字幕| 伦理电影免费视频| 天堂俺去俺来也www色官网| 少妇被粗大的猛进出69影院| 久久人妻熟女aⅴ| 欧美日韩福利视频一区二区| 侵犯人妻中文字幕一二三四区| 精品亚洲成国产av| 欧美激情高清一区二区三区| 精品国产一区二区久久| 99久久综合免费| 国产精品1区2区在线观看. | 亚洲精品粉嫩美女一区| 99精品久久久久人妻精品| av国产精品久久久久影院| 男人舔女人的私密视频| 色94色欧美一区二区| 精品一区二区三区四区五区乱码| 精品国产国语对白av| 国产野战对白在线观看| 久久中文字幕一级| e午夜精品久久久久久久| 91精品伊人久久大香线蕉| 国产高清视频在线播放一区 | 国产精品 欧美亚洲| 97精品久久久久久久久久精品| 日韩欧美一区视频在线观看| 51午夜福利影视在线观看| 夜夜夜夜夜久久久久| 一本久久精品| 在线观看一区二区三区激情| 国产成人精品久久二区二区91| 2018国产大陆天天弄谢| 国产亚洲av高清不卡| 午夜福利一区二区在线看| 90打野战视频偷拍视频| 国产一区二区在线观看av| 人妻人人澡人人爽人人| videos熟女内射| 久久中文字幕一级| 亚洲成国产人片在线观看| 热99国产精品久久久久久7| 青春草视频在线免费观看| 精品国产一区二区三区久久久樱花| 一区二区av电影网| 免费在线观看视频国产中文字幕亚洲 | 日本欧美视频一区| 久久女婷五月综合色啪小说| 99热网站在线观看| 女性被躁到高潮视频| av视频免费观看在线观看| 一区二区三区精品91| 国产精品一区二区免费欧美 | 亚洲性夜色夜夜综合| 久久香蕉激情| 搡老熟女国产l中国老女人| 久久精品熟女亚洲av麻豆精品| 纵有疾风起免费观看全集完整版| 欧美日韩亚洲国产一区二区在线观看 | 丰满迷人的少妇在线观看| 在线观看www视频免费| 日韩熟女老妇一区二区性免费视频| 国产成人av教育| 中文字幕精品免费在线观看视频| 久久亚洲精品不卡| 老汉色av国产亚洲站长工具| 午夜福利免费观看在线| 各种免费的搞黄视频| 成年人午夜在线观看视频| 精品乱码久久久久久99久播| 精品一区二区三区四区五区乱码| 丝袜脚勾引网站| 久久久久精品人妻al黑| 99热网站在线观看| 久久国产精品影院| 女人久久www免费人成看片| 亚洲精品乱久久久久久| 久久人人爽人人片av| 90打野战视频偷拍视频| 亚洲avbb在线观看| 免费在线观看影片大全网站| 国内毛片毛片毛片毛片毛片| 搡老乐熟女国产| 十八禁高潮呻吟视频| 亚洲精品自拍成人| 国产免费视频播放在线视频| av欧美777| 99热国产这里只有精品6| 中文字幕精品免费在线观看视频| 国产91精品成人一区二区三区 | 一区在线观看完整版| 两性午夜刺激爽爽歪歪视频在线观看 | 19禁男女啪啪无遮挡网站| 老司机午夜十八禁免费视频| a 毛片基地| 国产欧美日韩精品亚洲av| 国产成人欧美在线观看 | 国产福利在线免费观看视频| 国产成人av激情在线播放| 久久久国产欧美日韩av| 日韩大码丰满熟妇| 午夜91福利影院| 人妻 亚洲 视频| 别揉我奶头~嗯~啊~动态视频 | 国产三级黄色录像| 欧美另类亚洲清纯唯美| 黄色 视频免费看| 51午夜福利影视在线观看| 久久99热这里只频精品6学生| 国产日韩欧美在线精品| tocl精华| 国产老妇伦熟女老妇高清| 19禁男女啪啪无遮挡网站| 成人国产av品久久久| av有码第一页| 免费不卡黄色视频| 亚洲精品第二区| 国内毛片毛片毛片毛片毛片| 久久热在线av| 丁香六月天网| 大香蕉久久成人网| 女人爽到高潮嗷嗷叫在线视频| 久久亚洲国产成人精品v| 最新在线观看一区二区三区| 99热全是精品| 狠狠婷婷综合久久久久久88av| 精品第一国产精品| 国产日韩欧美在线精品| av福利片在线| 97精品久久久久久久久久精品| 成在线人永久免费视频| av在线播放精品| 人妻 亚洲 视频| 99热网站在线观看| 999久久久精品免费观看国产| 欧美xxⅹ黑人| 91精品伊人久久大香线蕉| 国产91精品成人一区二区三区 | 亚洲国产av新网站| 韩国高清视频一区二区三区| 成年美女黄网站色视频大全免费| 日日夜夜操网爽| 国产又色又爽无遮挡免| 狂野欧美激情性xxxx| 三上悠亚av全集在线观看| 亚洲av国产av综合av卡| 久久精品亚洲熟妇少妇任你| 女性生殖器流出的白浆| 亚洲欧洲精品一区二区精品久久久| 精品人妻一区二区三区麻豆| 欧美精品一区二区免费开放| 老司机亚洲免费影院| 亚洲av日韩精品久久久久久密| 久久久国产一区二区| 日韩欧美国产一区二区入口| 十八禁人妻一区二区| 成年人午夜在线观看视频| 大香蕉久久网| svipshipincom国产片| 天天躁夜夜躁狠狠躁躁| 亚洲激情五月婷婷啪啪| 老司机影院成人| 成在线人永久免费视频| 午夜免费鲁丝| 国产伦理片在线播放av一区| 亚洲国产精品一区三区| 少妇被粗大的猛进出69影院| 女人精品久久久久毛片| 18禁观看日本| 一区在线观看完整版| 黄网站色视频无遮挡免费观看| 国产成人精品无人区| netflix在线观看网站| 精品人妻熟女毛片av久久网站| 国产激情久久老熟女| tube8黄色片| 性少妇av在线| 999久久久国产精品视频| 青草久久国产| 精品福利观看| 免费观看人在逋| 亚洲成人手机| 在线天堂中文资源库| 精品国内亚洲2022精品成人 | 69av精品久久久久久 | 一区二区三区乱码不卡18| 亚洲精品粉嫩美女一区| 人人澡人人妻人| kizo精华| 夫妻午夜视频| 午夜福利乱码中文字幕| 热99re8久久精品国产| 国产精品熟女久久久久浪| www.av在线官网国产| 亚洲欧美日韩高清在线视频 | 欧美人与性动交α欧美软件| 日韩熟女老妇一区二区性免费视频| 人人澡人人妻人| 宅男免费午夜| av一本久久久久| 午夜福利乱码中文字幕| 久久久国产一区二区| 免费少妇av软件| 日韩中文字幕视频在线看片| 午夜福利视频在线观看免费| 日韩一卡2卡3卡4卡2021年| 成人国产一区最新在线观看| 欧美人与性动交α欧美精品济南到| 999久久久精品免费观看国产| 免费高清在线观看视频在线观看| 国产精品免费视频内射| av超薄肉色丝袜交足视频| 狂野欧美激情性bbbbbb| 国产成人精品久久二区二区91| 久久精品亚洲熟妇少妇任你| 美女脱内裤让男人舔精品视频| 真人做人爱边吃奶动态| 国产免费现黄频在线看| 大片免费播放器 马上看| 欧美激情 高清一区二区三区| av又黄又爽大尺度在线免费看| 男人舔女人的私密视频| 亚洲人成电影观看| 亚洲欧美精品综合一区二区三区| 国产福利在线免费观看视频| 亚洲欧美清纯卡通| 女性生殖器流出的白浆| 精品国产乱子伦一区二区三区 | 美女主播在线视频| 国产亚洲av高清不卡| 伊人亚洲综合成人网| 国产在视频线精品| 免费高清在线观看日韩| 国产成人免费无遮挡视频| 亚洲男人天堂网一区| 夫妻午夜视频| 日韩 欧美 亚洲 中文字幕| 首页视频小说图片口味搜索| 久久精品国产a三级三级三级| 在线观看免费午夜福利视频| 蜜桃国产av成人99| 亚洲美女黄色视频免费看| 麻豆乱淫一区二区| 麻豆av在线久日| 成年人午夜在线观看视频| 亚洲精品中文字幕一二三四区 | 视频区图区小说| 建设人人有责人人尽责人人享有的| 99国产精品一区二区蜜桃av | 日韩免费高清中文字幕av| 亚洲av电影在线进入| 亚洲精品av麻豆狂野|