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

    直墻前矩形月池水波振蕩研究分析

    2016-10-12 06:55:20鄭碧芳張洪生楊小巖
    海洋工程 2016年5期
    關(guān)鍵詞:直墻月池勢函數(shù)

    鄭碧芳,張洪生,楊小巖,2

    (1.上海海事大學(xué) 海洋科學(xué)與工程學(xué)院,上海 201306;2.上海交通大學(xué) 船舶海洋與建筑工程學(xué)院,上海 200240)

    直墻前矩形月池水波振蕩研究分析

    鄭碧芳1,張洪生1,楊小巖1,2

    (1.上海海事大學(xué) 海洋科學(xué)與工程學(xué)院,上海 201306;2.上海交通大學(xué) 船舶海洋與建筑工程學(xué)院,上海 200240)

    月池內(nèi)流體存在活塞和晃蕩兩類振蕩現(xiàn)象?;诰€性勢流理論,推導(dǎo)了波浪斜向入射下,直墻前矩形月池輻射和繞射問題的解析解。通過分離變量法和特征函數(shù)展開法求解了速度勢函數(shù),根據(jù)邊界條件來確定速度勢函數(shù)中的未知系數(shù),由速度勢函數(shù)計(jì)算斜向波與矩形月池相互作用的水動力系數(shù)和波浪激勵力,對它們的變化規(guī)律進(jìn)行了分析討論,研究了底部開口大小、波浪入射角度對矩形月池水動力特性的影響以及直墻遠(yuǎn)近對波浪力的影響。結(jié)果表明,月池底部開口大小對流體水平作用的影響較小,而對流體垂直作用的影響較大;波浪入射角度的變化對矩形月池橫蕩和橫搖運(yùn)動時(shí)的水動力特性有一定的影響;在一定條件下,直墻的存在會使得月池在水平方向所受到的波浪力比開敞水域中的要大。

    矩形月池;直墻;輻射;繞射;解析解;水動力特性

    Abstract:Two types of oscillation piston and sloshing,always exist in the fluid of moon pool.In the ocean,the wave is often obliquely incident.Based on the linear potential theory,an analytical solution is presented for the radiation and diffraction of a rectangular structure with a moon pool in front of a vertical wall in oblique seas.The expressions for the potentials are obtained by the method of separation of variables and the eigenfunction expansion,and the unknown coefficients are determined by the boundary conditions.The hydrodynamic coefficients and wave excitation forces are investigated,and the effects of the width of opening,the angle of incidence and the distance between the vertical wall and the structure on the hydrodynamic characteristics of the rectangular structure with a moon pool are studied mainly.The analytical solution indicates that the opening has little effect on the horizontal force but a large effect on the vertical force.And the angle of incidence has some impacts on the sway and roll motion of the floating body.The horizontal force on the structure in the presence of a wall is larger than those in the absence of the vertical wall due to the reflection of the vertical wall.

    Keywords:rectangular structure with a moon pool; vertical wall; radiation; diffraction; analytical solution; hydrodynamic characteristic

    船舶與海洋工程中經(jīng)常應(yīng)用矩形月池結(jié)構(gòu),所以研究矩形月池與波浪的相互作用是一個(gè)很重要的水動力問題。關(guān)于波浪對開敞水域中結(jié)構(gòu)物的作用,學(xué)者們已進(jìn)行了一定的研究。黃磊等[1]研究了具有矩形月池的鉆井船在微幅橫蕩簡諧運(yùn)動激勵下的輻射問題,討論了月池參數(shù)對船舶水動力系數(shù)的影響;ZHENG等[2]采用分離變量法分析了自由表面上的無限長矩形浮體的波浪輻射和繞射問題;ZHOU等[3]研究了矩形月池的波浪輻射和繞射問題;MOLIN[4]研究了月池中的活塞和晃蕩現(xiàn)象;TEUNG和SEAH[5]采用特征函數(shù)展開的方法調(diào)查了月池中的共振現(xiàn)象,重點(diǎn)研究了共振頻率附近水動力參數(shù)的變化情況;康莊和姚熊亮[6]研究了月池內(nèi)流體的活塞振蕩特性;戴偉[7]研究了均勻流條件下月池的水動力特性;AALBERS[8]研究了波浪作用下月池內(nèi)的共振問題,并分析了月池對船舶垂蕩運(yùn)動的影響;FALTINSON等[9]研究了流體在由兩個(gè)矩形船體形成的月池中的二維活塞運(yùn)動。鑒于實(shí)際海洋中波浪多為斜向入射,一些學(xué)者研究了斜向波入射條件下,波浪與浮式結(jié)構(gòu)物的相互作用問題。ZHENG等[10-11]研究了波浪斜向入射下,矩形浮體以及淹沒在有限水深的矩形結(jié)構(gòu)的波浪輻射和繞射問題;YANG和ZHANG[12]研究了波浪斜向入射下,矩形月池的波浪輻射和繞射問題;EVANS和MORRIS[13]研究了波浪斜向入射條件下,位于深水處的垂直固定板的波浪反射和透射問題,文獻(xiàn)[14-18]研究了斜向波與柱體結(jié)構(gòu)物相互作用的問題。

    直墻的存在會影響波動特性,而且結(jié)構(gòu)物上的波浪力也將發(fā)生變化。部分學(xué)者采用解析法研究了波浪與直墻前結(jié)構(gòu)物的作用問題。HSU 和 WU[19]采用分離變量法研究了直墻前二維矩形浮標(biāo)的垂蕩問題,發(fā)現(xiàn)水動力系數(shù)的一系列峰值會隨浮標(biāo)與直墻間距離的變化而有規(guī)律地出現(xiàn);ZHENG 等[20]進(jìn)一步研究了直墻前二維矩形浮標(biāo)的三個(gè)自由度的輻射問題,重點(diǎn)研究了吃水、浮標(biāo)與直墻間距離等對水動力系數(shù)的影響;滕斌和寧德志[21-22]應(yīng)用映像原理,將直墻前單個(gè)圓柱體對波浪的繞射問題變換為雙柱對雙向波浪的繞射問題,利用特征函數(shù)展開法,建立了直墻前圓柱體波浪繞射的解析求解過程;還采用同樣的方法研究了正交直墻前圓柱的振蕩問題;ZHANG和ZHOU[23]研究了直墻前矩形月池的波浪輻射和繞射問題,并分析了開口位置和直墻的存在對月池水動力特性的影響。

    本文推導(dǎo)了有限水深中、波浪斜向入射條件下直墻前無限長矩形月池的輻射和繞射問題的解析解。首先,給出入射波速度勢函數(shù)的表達(dá)式,把整個(gè)水域劃分為若干個(gè)合適的子區(qū)域;然后,根據(jù)邊界條件,應(yīng)用分離變量法和特征函數(shù)展開法求解各子區(qū)域的速度勢函數(shù);最后,分析討論底部開口大小、波浪入射角度以及直墻遠(yuǎn)近對矩形月池水動力系數(shù)和波浪激勵力的影響。

    圖1 直墻前矩形月池示意Fig.1 The sketch of a rectangular structure with a moon pool in front of a vertical wall

    1 數(shù)學(xué)模型

    一個(gè)寬度為2b,吃水為h1的無限長的矩形結(jié)構(gòu)物漂浮在深度為h的水中,右側(cè)D處有一在y方向上無限長的直墻。在矩形結(jié)構(gòu)物的底部中心處有一個(gè)寬度為2a的開口。在未受擾動的靜水面上建立如圖1所示的笛卡爾坐標(biāo)系O-xyz。z軸坐標(biāo)向上為正,x軸方向向右。假定矩形月池在y軸方向上無限長且入射波與x軸形成一個(gè)夾角θ(0°<θ<90°)。

    假定流體無粘、不可壓縮且運(yùn)動無旋,入射波和結(jié)構(gòu)物的振幅很小,總的速度勢Φ可以表示為[12]:

    (1)

    (2)

    和邊界條件:

    (3)

    (4)

    (5)

    (6)

    式中:ν=ω2/g,g是重力加速度,Sb是月池的濕表面,δjl為Kronecker 函數(shù),其定義為:

    (7)

    頻散關(guān)系為:

    ω2=gkthkh

    (8)

    φ(0)可表示為:

    (9)

    物面邊界條件展開為[3]:

    (10)

    2 求解過程

    2.1速度勢函數(shù)的表達(dá)式

    用特征函數(shù)展開法求解上述邊值問題的解析解,把整個(gè)水域劃分為圖1所示的七個(gè)子區(qū)域,分別是I (b≤x≤b+D,-h≤z≤0)、II (a≤x≤b,-h≤z≤-h1)、III (-a≤x≤a,-h≤z≤0)、IV (-b≤x≤-a,-h≤z≤-h1)、V (-≤x≤-b,-h≤z≤0)、VI (a≤x≤b,-h1≤z≤0)和VII (-b≤x≤-a,-h1≤z≤0)。七個(gè)子區(qū)域內(nèi)的輻射速度勢和繞射速度勢分別記為:和

    根據(jù)分離變量法和各區(qū)域之間的速度和壓力連續(xù)邊界條件,每個(gè)子區(qū)域的速度勢函數(shù)可以表示為:

    (13)

    (14)

    (15)

    (16)

    (17)

    (18)

    (19)

    其中:

    (20)

    (21)

    (22)

    βn=(n-1)π/(h-h1)n=1,2,3,...

    (23)

    (24)

    α1=-ik1,k1th(k1h1)=ω2/gn=1
    αnth(αnh1)=-ω2/gn=2,3,…

    (25)

    (26)

    k0=ksinθ

    (27)

    (28)

    (29)

    其中:

    (30)

    (31)

    2.2未知系數(shù)的求解

    在x=±b和x=±a斷面,運(yùn)用速度連續(xù)和壓力連續(xù)條件,可以得到:

    (32)

    (33)

    (34)

    (35)

    (36)

    (37)

    (38)

    (39)

    (40)

    (41)

    (42)

    (43)

    根據(jù)文獻(xiàn)[3],將含有未知數(shù)的速度勢函數(shù)表達(dá)式代入式(32)~(43)中,利用三角函數(shù)的正交性,可將以上邊值問題轉(zhuǎn)換為一個(gè)無窮維線性方程組進(jìn)行求解。為得到未知系數(shù),需把所有的無窮序列進(jìn)行截?cái)嗵幚?,取無窮級數(shù)的前N項(xiàng),即產(chǎn)生一個(gè)12N階的線性方程組。通過高斯消元法求解方程組,即可得到速度勢函數(shù)的表達(dá)式。

    3 水動力計(jì)算

    3.1輻射作用力

    在得到速度勢函數(shù)的表達(dá)式后,由Bernoulli方程可得到月池濕表面的壓強(qiáng)分布,沿月池表面對壓強(qiáng)進(jìn)行積分,可得到月池上的輻射波浪力為:

    (44)

    式中:ρ是密度,τij是在模態(tài)j作單位振幅的運(yùn)動時(shí)所受到的i方向的輻射作用力。

    由波浪輻射力可以得到附加質(zhì)量μij和阻尼系數(shù)λij的表達(dá)式,可將它們無量綱化為:

    (45)

    (46)

    3.2波浪激勵力

    矩形月池所受到的波浪激勵力可以由入射和繞射速度勢求得,垂直分力為:

    (47)

    水平分力為:

    (48)

    以及關(guān)于(0,0)的傾覆力矩為:

    (49)

    將波浪激勵力(矩)無量綱化為:

    fi=Fi/(2ρgη(0)bbδ3i)

    (50)

    4 分析與討論

    通過計(jì)算斜向波與直墻前矩形月池相互作用的水動力系數(shù)和波浪激勵力,對得到的速度勢函數(shù)的正確性進(jìn)行驗(yàn)證,然后定性分析波浪入射角度、底部開口大小對月池水動力特性的影響以及直墻遠(yuǎn)近對波浪力的影響。

    4.1結(jié)果驗(yàn)證

    4.1.1 收斂性

    在速度勢函數(shù)的計(jì)算中,已有學(xué)者通過數(shù)值試驗(yàn)的方法考察了取不同N值時(shí)水動力系數(shù)和波浪激勵力計(jì)算結(jié)果的收斂情況。文獻(xiàn)[3]分別討論了N=50和N=100時(shí)附加質(zhì)量、阻尼系數(shù)曲線和波浪激勵力的變化曲線,發(fā)現(xiàn)在這兩種情況下得到的曲線吻合良好,可知N取50時(shí)即有足夠的運(yùn)算精度。下文中所有的水動力系數(shù)和波浪力都是在N=50時(shí)計(jì)算得到的。

    4.1.2 水動力系數(shù)的對稱性驗(yàn)證

    由格林公式知,附加質(zhì)量和阻尼系數(shù)均具有對稱性[24]。利用這一特性可對所得速度勢函數(shù)的正確性進(jìn)行驗(yàn)證。圖2給出了橫搖和橫蕩模態(tài)兩種組合下求得的附加質(zhì)量和阻尼系數(shù)曲線,由圖可知這兩種組合下的水動力系數(shù)曲線吻合良好。

    圖2 θ=30°,h/h1=3,b/h1=0.5,a/b=0.5和D/h1=3時(shí),橫蕩和橫搖運(yùn)動組合下附加質(zhì)量和阻尼系數(shù)的對稱性Fig.2 Symmetry of the added mass and damping coefficients between coupled sway and roll motions withθ=30°,h/h1=3,b/h1=0.5,a/b=0.5and D/h1=3

    4.2入射角度對水動力特性的影響

    考慮入射角θ分別取15°、45°和75°時(shí),入射角度的大小對直墻前矩形月池的水動力系數(shù)的影響。在垂蕩模態(tài)下,不同入射角度引起的附加質(zhì)量和阻尼系數(shù)曲線如圖3(a)和3(b)所示。在圖3(a)中θ分別取15°、45°和75°時(shí),μ11分別在2kb≈ 0.575、0.617和0.662處取得峰值,圖3(a)和3(b)表明在垂蕩模態(tài)下,不同入射角度引起的附加質(zhì)量和阻尼系數(shù)的變化幾乎是相同的。圖3(c)~3(d),3(e)~3(f)表明在橫蕩和橫搖模態(tài)下,不同入射角度引起的附加質(zhì)量和阻尼系數(shù)的變化與垂蕩模態(tài)下的變化有很大不同:入射角度不同附加質(zhì)量也不相同,且主要表現(xiàn)在2kb> 2.25時(shí);而阻尼系數(shù)的區(qū)別相對很小,在不同入射角度下曲線基本保持一致,只在峰值處有略微的區(qū)別。更具體地說,從圖3(c)和3(e)可以看到入射角θ分別取15°、45°和75°時(shí)附加質(zhì)量曲線波動變化很大,不同入射角下其峰值也有很大的差別。所有這些都表明入射角度的不同會對橫蕩和橫搖運(yùn)動的共振頻率產(chǎn)生重要的影響。

    圖3 h/h1=3,b/h1=0.5,a/b=0.5和D/h1=3時(shí)入射波角度對水動力系數(shù)的影響Fig.3 Effect of the angle of incidence on hydrodynamic coefficients with h/h1=3,b/h1=0.5,a/b=0.5and D/h1=3

    不同入射角度下垂直分力的實(shí)部和虛部曲線分別如圖4(a)和4(b)所示。在圖4(a)中當(dāng)入射角θ分別取15°、45°和75°時(shí),峰值的絕對值是0.906、0.660和0.234,依次下降了27.2%和64.5%,由圖4(b)可知峰值的絕對值是0.852、0.573和0.155,依次下降了32.7%和72.9%。圖4(c)~4(d),4(e)~4(f)表明入射角度不同時(shí)波浪激勵力實(shí)部和虛部的曲線變化也有很大的不同。圖4(c)和4(e)表明,當(dāng)入射角θ分別取15°、45°和75°時(shí),水平分力和傾覆力矩實(shí)部峰值的絕對值在2kb= 0.535時(shí)分別是1.891、1.607、1.186和1.527、1.298、0.958,分別依次下降了15.0%、26.2%和15.0%、26.2%。圖4(d)和4(f)表明,水平分力和傾覆力矩虛部的第一個(gè)峰值所對應(yīng)的頻率分別與它們實(shí)部的第一個(gè)峰值所對應(yīng)的頻率是一樣的;而虛部曲線隨著角度的增加則出現(xiàn)很大的差異。

    圖4 h/h1=3,b/h1=0.5,a/b=0.5和D/h1=3時(shí)入射波的角度對波浪激勵力的影響Fig.4 Effect of the angle of incidence on wave excitation forces with h/h1=3,b/h1=0.5,a/b=0.5and D/h1=3

    圖5給出了直墻存在的情況下,矩形月池所受波浪力的垂向分力和水平分力隨入射角度的變化規(guī)律。從圖中可以看出,盡管2kb取不同的值,但垂向波浪力和水平方向的波浪力都隨著入射角度的增加而逐漸減小;入射角度在小于45°時(shí)對波浪力的影響相對不太明顯,但在45°~90°范圍內(nèi),波浪力隨著角度的增加而迅速減小。

    圖5 h/h1=3,b/h1=0.5,a/b=0.5和D/h1=3時(shí)波浪力隨入射角度的變化Fig.5 Variation of the wave excitation forces with the angle of incidence in the case of h/h1=3,b/h1=0.5,a/b=0.5and D/h1=3

    4.3底部開口大小對水動力特性的影響

    對a/b取不同值如0.2、0.4、0.6和0.8時(shí)分別進(jìn)行計(jì)算。a/b取不同值時(shí)垂蕩運(yùn)動引起的附加質(zhì)量力和阻尼系數(shù)的結(jié)果如圖6(a)和6(b)所示。由圖6(a)可知,當(dāng)a/b取0.2、0.4、0.6和0.8時(shí),μ11分別在2kb約為0.498、0.535、0.617和0.617處取得峰值。圖6(b)表明,雖然開口寬度不同,但垂蕩運(yùn)動引起的阻尼系數(shù)的變化幾乎是相同的。圖6(c)和6(e)表明,當(dāng)取不同的a/b值時(shí),橫蕩和橫搖模態(tài)下的附加質(zhì)量變化有很大不同,在2kb=3~4之間,曲線出現(xiàn)不同的峰值。圖6(d)和6(f)表明,不同開口大小對阻尼系數(shù)的影響不大。

    圖6 θ=30°,h/h1=3,b/h1=0.5和D/h1=3時(shí)底部開口寬度對水動力系數(shù)的影響Fig.6 Effect of the width of opening on hydrodynamic coefficients with θ=30°,h/h1=3,b/h1=0.5and D/h1=3

    a/b取不同值時(shí),垂直分力的實(shí)部和虛部曲線分別如圖7(a)和7(b)所示。因?yàn)檫@個(gè)力是來自水平底部的壓力,所以底部開口的寬度對其有很大的影響。由圖7(a)可知,當(dāng)a/b分別為0.2、0.4、0.6和0.8時(shí),Re(f1)對應(yīng)峰值的絕對值分別是2.864、1.111、0.518和0.138,依次下降了 61.2%、53.3%和73.4%,且從圖7(b)可知Im(f1)峰值的絕對值分別是1.684、1.020、0.473和0.127,依次下降了39.4%、53.6%和73.2%。圖7(c)~7(d),7(e)~7(f)表明,水平分力、傾覆力矩的實(shí)部曲線在取不同的a/b時(shí)基本都是一致的;但虛部由于底部開口大小不同,曲線峰值的大小也明顯不同。

    圖7 θ=30°,h/h1=3,b/h1=0.5和D/h1=3時(shí)底部開口寬度對波浪激勵力的影響Fig.7 Effect of the width of opening on wave excitation forces with θ=30°,h/h1=3,b/h1=0.5andD/h1=3

    4.4直墻遠(yuǎn)近對波浪力的影響

    圖8是在波浪正向入射條件下,矩形月池與直墻間取不同距離(即l/h1取不同值)時(shí),作用于月池水平方向波浪力的分布,以及和開敞水域中矩形月池所受波浪力的對比情況。由圖可以看到,存在直墻時(shí)矩形月池所受的最大波浪力相較于開敞水域中矩形月池所受的最大波浪力增大了兩倍左右。由于直墻的存在,結(jié)構(gòu)物上的波浪力隨著波數(shù)的變化而發(fā)生振蕩,且矩形月池與直墻間距離較大時(shí),反射作用減小,其所受到的水平方向波浪力的峰值減小。圖8(a)和8(b)表明波浪的水平分力曲線在有無直墻時(shí)的整體走勢大致相同,但在直墻結(jié)構(gòu)存在的情況下,當(dāng)2kb<1 時(shí),結(jié)構(gòu)所受的波浪水平力曲線出現(xiàn)多個(gè)波動。

    圖8 h/h1=3,b/h1=0.5,a/b=0.5時(shí)直墻距離對波浪力的影響Fig.8 Effect of the distance between the vertical wall and the structure on wave excitation forces with h/h1=3,b/h1=0.5,a/b=0.5

    5 結(jié) 語

    本文通過分離變量法和特征函數(shù)展開法,推導(dǎo)了有限水深中、波浪斜向入射條件下直墻前矩形月池輻射和繞射問題的解析解。利用所推導(dǎo)的解析解,計(jì)算了作用在矩形月池上的波浪激勵力,且主要研究了波浪入射角度、底部開口寬度對矩形月池水動力特性的影響以及直墻遠(yuǎn)近對波浪力的影響。研究表明入射角度對阻尼系數(shù)以及垂蕩運(yùn)動的附加質(zhì)量影響不大,但對橫蕩和橫搖運(yùn)動的附加質(zhì)量有很大的影響;入射角度對波浪激勵力的垂直分力影響不明顯,隨著入射角度的增加水平分力和傾覆力矩會逐漸變小,因此入射角度對水平分力和傾覆力矩的影響很大。底部開口寬度對阻尼系數(shù)的影響不明顯,但對附加質(zhì)量的影響很大,不同的開口寬度引起的附加質(zhì)量峰值明顯不同;由于月池底部開口,月池內(nèi)外流體的運(yùn)動可視為整體豎向的運(yùn)動,因此底部開口寬度對月池的水平力和晃動頻率基本沒有影響,但對于月池的垂直力和活塞頻率有明顯的影響。月池所受到的水平波浪力的峰值由于直墻結(jié)構(gòu)的存在而增大,當(dāng)2kb<1時(shí)水平向波浪力出現(xiàn)多個(gè)波動。由于本文是基于線性勢流理論采用解析方法而進(jìn)行計(jì)算求解的,因此還需對計(jì)算結(jié)果、尤其是對于入射角大于45°以及a/b小于0.2時(shí)的計(jì)算結(jié)果,通過物理模型實(shí)驗(yàn)或者數(shù)值解進(jìn)行進(jìn)一步的驗(yàn)證。本文得到的矩形月池的水動力系數(shù)和波浪力的變化規(guī)律可為月池的設(shè)計(jì)提供重要參考。

    [1] 黃磊,劉利琴,唐友剛,等.水平激勵下鉆井船矩形月池的水動力特性研究[J].天津大學(xué)學(xué)報(bào):自然科學(xué)與工程技術(shù)版,2015,48(11):1 001-1 008.(HUANG Lei,LIU Liqin,TANG Yougang,et al..Hydrodynamic characteristics of rectangular moonpool of drilling unit under horizontal excitation [J].Journal of Tianjin University,2015,48(11):1 001-1 008.(in Chinese))

    [2] ZHENG Y H,YOU Y G,SHEN Y M.On the radiation and diffraction of water waves by a rectangular buoy [J].Ocean Engineering,2004,31(8):1 063-1 082.

    [3] ZHOU H W,WU G X,ZHANG H S.Wave radiation and diffraction by a two-dimensional floating rectangular body with an opening in its bottom [J].Journal of Engineering Mathematics,2013,83(1):1-22.

    [4] MOLIN B.On the piston and sloshing modes in moonpools [J].Journal of Fluid Mechanics,2001,430(1):27-50.

    [5] YEUNG R W,SEAH R K M.On Helmholtz and higher-order resonance of twin floating bodies[J].Journal of Engineering Mathematics,2007,58 (1-4):251-265.

    [6] 康莊,姚熊亮.浪流條件下月池內(nèi)流體活塞振蕩特性實(shí)驗(yàn)研究[J].哈爾濱工程大學(xué)學(xué)報(bào),2008,29(3):209-216.(KANG Zhuang,YAO Xiongliang.Experimental research on piston oscillating behavior of fluid in a moon pool with wave-flow [J].Journal of Harbin Engineering University,2008,29(3):209-216.(in Chinese))

    [7] 戴偉.均勻流條件下月池水動力特性研究[D].哈爾濱:哈爾濱工程大學(xué),2005.(DAI Wei.Hydrodynamic performance research on moon pool under uniform flow condition [D].Harbin:Harbin Engineering University,2005.(in Chinese))

    [8] ALBERS A B.The water motions in a moonpool[J].Ocean Engineering,1984,11(6):557-579.

    [9] FALTINSEN O M,ROGNEBAKKE O F,TIMOKHA A N.Two-dimensional resonant piston-like sloshing in a moon pool [J].Journal of Fluid Mechanics,2007,575(9):359-397.

    [10] ZHENG Y H,SHEN Y M,YOU Y G,et al.Wave radiation by a floating rectangular structure in oblique seas [J].Ocean Engineering,2006,33(1):59-81.

    [11] ZHENG Y H,LIU P F,SHEN Y M,et al.On the radiation and diffraction of linear water waves by an infinitely long rectangular structure submerged in oblique seas [J].Ocean Engineering,2007,34(3-4):436-450.

    [12] YANG X Y,ZHANG H S.Wave radiation and diffraction by a floating rectangular structure with an opening at its bottom in oblique seas [J].Journal of Hydrodynamics,ser.B.2015,Accepted.

    [13] EVANS D V,MORRIS C A N.The effects of a fixed vertical barrier on obliquely incident surface waves in deep water [J].Journal of Institute of Mathematics and its Applications,1972,9(2):198-204.

    [14] POLITIS C G,PAPALEXANDRIS M V,ATHANASSOULIS G A.A boundary integral equation method for oblique water-wave scattering by cylinders governed by the modied Helmholtz equation [J].Applied Ocean Research,2002,24(4):215-233.

    [15] BOLTON W E,URSELL F.The wave force on an infinitely long circular cylinder in an oblique sea [J].Journal of Fluid Mechanics,1973,57(2):241-256.

    [16] BAI K J.Diffraction of oblique waves by an infinite cylinder [J].Journal of Fluid Mechanics,1975,68(3):513-535.

    [17] LEONARD J W,HUANG M C,HUDSPETH R T.Hydrodynamic interference between floating cylinders in oblique seas [J].Applied Ocean Research,1983,5(3):158-166.

    [18] GARRISON C J.Interaction of oblique waves with an infinite cylinder [J].Applied Ocean Research,1984,6(1):4-15.

    [19] HSU H H,WU Y C.The hydrodynamic coefficients for an oscillating rectangular structure on a free surface with sidewall.Ocean Engineering,1997,24 (2):177-199.

    [20] ZHENG Y H,SHEN Y M,YOU Y G,et al.On the radiation and diffraction of water waves by a rectangular structure with a sidewall[J].Ocean Engineering,2004,31 (17):2 087-2 104.

    [21] 滕斌,寧德志.波浪對直墻前垂直圓柱的繞射[J].海洋工程,2003,21(4):48-52.(TENG Bin,NING Dezhi.Wave diffraction from a uniform cylinder in front of a vertical wall [J].The Ocean Engineering,2003,21(4):48-52.(in Chinese))

    [22] 寧德志,滕斌.正交直墻前直立圓柱振蕩問題的研究[J].中國工程科學(xué),2003,5(5):84-91.(NING Dezhi,TENG Bin.Study on the oscillation of a uniform cylinder in front of two verticd walls intersecting mormally[J].Engineering Sciences,2003,5 (5):84-91.(in Chinese))

    [23] ZHANG H S,ZHOU H W.Wave radiation and diffraction by a two-dimensional floating body with an opening near a side wall [J].China Ocean Engineering,2013,27(4):437-450.

    [24] MARUO H.The drift of a body floating on waves [J].Journal of Ship Research,1960,4 (03):1-5.

    The oscillation of water waves in a rectangular structure with a moon pool in front of a vertical wall

    ZHENG Bifang1,ZHANG Hongsheng1,YANG Xiaoyan1,2

    (1.College of Ocean Science and Engineering,Shanghai Maritime University,Shanghai 201306,China; 2.School of Naval Architecture,Ocean and Civil Engineering,Shanghai Jiao Tong University,Shanghai 200240,China )

    P731.2;O353.2

    A

    10.16483/j.issn.1005-9865.2016.05.005

    張洪生。 E-mail:hszhang@shmtu.edu.cn

    1005-9865(2016)05-0040-11

    2015-12-08

    國家自然科學(xué)基金資助項(xiàng)目(51079082);上海市自然科學(xué)基金資助項(xiàng)目(14ZR1419600);上海市教委一流學(xué)科建設(shè)項(xiàng)目;2013年上海市研究生教育創(chuàng)新計(jì)劃實(shí)施項(xiàng)目(第二批)

    鄭碧芳(1992-),女,山西應(yīng)縣人,碩士研究生,主要從事水波動力學(xué)的研究。

    猜你喜歡
    直墻月池勢函數(shù)
    航天器姿態(tài)受限的協(xié)同勢函數(shù)族設(shè)計(jì)方法
    “直墻平底”之于碑刻·篆刻·書法的藝術(shù)價(jià)值
    中國篆刻(2023年12期)2024-01-22 10:53:18
    大洋勘探船月池阻力特性研究
    數(shù)學(xué)理論與應(yīng)用(2022年1期)2022-04-15 09:03:32
    直墻半圓拱巷道圍巖應(yīng)力分布解析
    深埋直墻拱形隧道穩(wěn)定性研究
    金屬鎢級聯(lián)碰撞中勢函數(shù)的影響
    不同月池形狀對鉆井船附加阻力的影響
    鉆井船月池區(qū)波面升高運(yùn)動分析
    直墻半圓拱可縮性U型鋼支架卡纜臨界約束力理論研究
    99热国产这里只有精品6| 十八禁网站网址无遮挡| 一区二区三区激情视频| 国产成人精品久久二区二区免费| 夜夜夜夜夜久久久久| 汤姆久久久久久久影院中文字幕| 啦啦啦 在线观看视频| 国产xxxxx性猛交| 一区二区三区激情视频| 老司机深夜福利视频在线观看| 欧美久久黑人一区二区| 无人区码免费观看不卡 | 一本色道久久久久久精品综合| 久久精品国产亚洲av高清一级| 精品乱码久久久久久99久播| 亚洲美女黄片视频| 亚洲精品在线美女| 777久久人妻少妇嫩草av网站| 国产在视频线精品| 男女高潮啪啪啪动态图| 精品国产一区二区久久| 国产黄色免费在线视频| 久久中文字幕一级| kizo精华| 涩涩av久久男人的天堂| 九色亚洲精品在线播放| 夜夜骑夜夜射夜夜干| 91大片在线观看| 人人妻人人澡人人看| 桃红色精品国产亚洲av| 性色av乱码一区二区三区2| 亚洲av日韩精品久久久久久密| 国产av一区二区精品久久| 91字幕亚洲| 亚洲五月色婷婷综合| 亚洲欧美色中文字幕在线| 精品少妇久久久久久888优播| 国产极品粉嫩免费观看在线| 一级,二级,三级黄色视频| 国产精品一区二区在线观看99| kizo精华| 男女无遮挡免费网站观看| 伊人久久大香线蕉亚洲五| 欧美日韩精品网址| 在线观看舔阴道视频| 国产成人啪精品午夜网站| 操美女的视频在线观看| 成年人黄色毛片网站| 久久99一区二区三区| 99热网站在线观看| 1024视频免费在线观看| 久久精品国产亚洲av高清一级| 黄色视频不卡| av视频免费观看在线观看| 国产日韩欧美亚洲二区| 一区二区三区激情视频| 欧美性长视频在线观看| 啦啦啦在线免费观看视频4| 亚洲第一欧美日韩一区二区三区 | 99精品久久久久人妻精品| 又大又爽又粗| 精品福利永久在线观看| 免费观看人在逋| 成人亚洲精品一区在线观看| 亚洲第一青青草原| 黑人巨大精品欧美一区二区蜜桃| 午夜视频精品福利| 国产精品亚洲一级av第二区| 久久精品亚洲熟妇少妇任你| 成人精品一区二区免费| 亚洲欧美日韩高清在线视频 | 男人舔女人的私密视频| 亚洲欧美日韩另类电影网站| 国产精品国产高清国产av | 成人三级做爰电影| 黑人欧美特级aaaaaa片| 淫妇啪啪啪对白视频| 少妇的丰满在线观看| 国产成人av教育| 日韩大码丰满熟妇| 久久久久久久精品吃奶| 90打野战视频偷拍视频| 国产精品九九99| 欧美黄色片欧美黄色片| 一区福利在线观看| 久久久国产欧美日韩av| 999久久久国产精品视频| 亚洲欧美日韩另类电影网站| 国产精品熟女久久久久浪| 午夜日韩欧美国产| 91大片在线观看| 精品午夜福利视频在线观看一区 | 欧美成人免费av一区二区三区 | 亚洲欧美一区二区三区黑人| 19禁男女啪啪无遮挡网站| 午夜福利,免费看| 亚洲精品中文字幕一二三四区 | av不卡在线播放| 老司机在亚洲福利影院| 欧美精品亚洲一区二区| 国产精品99久久99久久久不卡| 欧美性长视频在线观看| 欧美日韩国产mv在线观看视频| 91精品国产国语对白视频| 午夜福利视频精品| 亚洲专区字幕在线| 岛国在线观看网站| 国产亚洲av高清不卡| 亚洲专区国产一区二区| 亚洲精品成人av观看孕妇| 动漫黄色视频在线观看| 12—13女人毛片做爰片一| 欧美日韩亚洲高清精品| 成人av一区二区三区在线看| 亚洲精品在线美女| 美女视频免费永久观看网站| av线在线观看网站| 亚洲少妇的诱惑av| 男女无遮挡免费网站观看| h视频一区二区三区| 午夜福利乱码中文字幕| 国产成人影院久久av| 大香蕉久久成人网| 两性夫妻黄色片| 日韩一卡2卡3卡4卡2021年| 十分钟在线观看高清视频www| 黄片大片在线免费观看| 国产伦人伦偷精品视频| 免费观看av网站的网址| 香蕉丝袜av| 少妇裸体淫交视频免费看高清 | 在线十欧美十亚洲十日本专区| 亚洲色图综合在线观看| 真人做人爱边吃奶动态| 国产成人免费无遮挡视频| 又黄又粗又硬又大视频| 麻豆乱淫一区二区| 伊人久久大香线蕉亚洲五| 久久久久久久国产电影| 日韩制服丝袜自拍偷拍| 天天躁狠狠躁夜夜躁狠狠躁| 91成人精品电影| 精品视频人人做人人爽| 国产精品香港三级国产av潘金莲| 国产欧美日韩一区二区三| 国产精品久久久人人做人人爽| 国产一区有黄有色的免费视频| 欧美+亚洲+日韩+国产| 女人被躁到高潮嗷嗷叫费观| 黄色 视频免费看| 日本五十路高清| av欧美777| 亚洲视频免费观看视频| 成年人免费黄色播放视频| 欧美日本中文国产一区发布| 女性生殖器流出的白浆| 亚洲七黄色美女视频| 亚洲黑人精品在线| 亚洲国产精品成人综合色| 色综合站精品国产| 最近视频中文字幕2019在线8| 午夜亚洲福利在线播放| 国产精品久久久久久人妻精品电影| 狠狠狠狠99中文字幕| 99riav亚洲国产免费| 精品熟女少妇八av免费久了| 很黄的视频免费| 精品日产1卡2卡| 中文字幕最新亚洲高清| 真实男女啪啪啪动态图| 99久久精品一区二区三区| 757午夜福利合集在线观看| 亚洲熟女毛片儿| 国产一区在线观看成人免费| 亚洲精品美女久久av网站| 国产精品av久久久久免费| 亚洲黑人精品在线| 亚洲人成电影免费在线| 亚洲成a人片在线一区二区| 欧洲精品卡2卡3卡4卡5卡区| 高潮久久久久久久久久久不卡| 国产亚洲精品一区二区www| 久久香蕉精品热| 午夜影院日韩av| 国产蜜桃级精品一区二区三区| 88av欧美| 99国产综合亚洲精品| 午夜福利在线观看吧| 无人区码免费观看不卡| 熟女少妇亚洲综合色aaa.| 99热这里只有精品一区 | 国产精品永久免费网站| 久久久久久久久免费视频了| 亚洲在线自拍视频| 欧美绝顶高潮抽搐喷水| 全区人妻精品视频| 国产精品影院久久| av欧美777| 亚洲精品粉嫩美女一区| 久久久精品欧美日韩精品| 欧美日韩乱码在线| 最近在线观看免费完整版| 国产成人av教育| 日本在线视频免费播放| 免费一级毛片在线播放高清视频| 一区二区三区国产精品乱码| 久久精品亚洲精品国产色婷小说| 日韩欧美 国产精品| 欧美在线黄色| 国产黄a三级三级三级人| 叶爱在线成人免费视频播放| 精品久久久久久久末码| 1024手机看黄色片| 欧美中文日本在线观看视频| 免费在线观看日本一区| 久久久久免费精品人妻一区二区| 亚洲国产高清在线一区二区三| 亚洲精品在线观看二区| 国产主播在线观看一区二区| 亚洲av电影在线进入| 久久久久国内视频| 国产成人一区二区三区免费视频网站| 成人三级做爰电影| 久久久精品大字幕| 国产免费男女视频| 变态另类成人亚洲欧美熟女| 九九在线视频观看精品| 国产精品影院久久| 免费av毛片视频| 舔av片在线| 午夜福利在线观看吧| 国产精品美女特级片免费视频播放器 | 99精品欧美一区二区三区四区| 国产精品av久久久久免费| 一边摸一边抽搐一进一小说| 欧美一级a爱片免费观看看| 欧美激情久久久久久爽电影| 精品电影一区二区在线| 在线观看一区二区三区| 19禁男女啪啪无遮挡网站| 亚洲中文av在线| 国产欧美日韩精品亚洲av| 一边摸一边抽搐一进一小说| 欧美日韩一级在线毛片| 亚洲第一欧美日韩一区二区三区| 亚洲人成网站高清观看| 黄片大片在线免费观看| tocl精华| 亚洲男人的天堂狠狠| 久久久久久久久中文| 精品不卡国产一区二区三区| 国产男靠女视频免费网站| 午夜亚洲福利在线播放| 日韩欧美国产在线观看| 日本黄色视频三级网站网址| 舔av片在线| 狂野欧美激情性xxxx| 99视频精品全部免费 在线 | 黑人欧美特级aaaaaa片| 午夜免费成人在线视频| 天天躁日日操中文字幕| 99久久综合精品五月天人人| 日韩国内少妇激情av| 婷婷六月久久综合丁香| 亚洲av电影在线进入| 国产成人aa在线观看| 国产伦一二天堂av在线观看| 97人妻精品一区二区三区麻豆| 国产成人精品无人区| 国产97色在线日韩免费| 少妇的逼水好多| 两个人的视频大全免费| 一进一出好大好爽视频| 欧美成人性av电影在线观看| 99久久精品国产亚洲精品| 国产熟女xx| 亚洲美女黄片视频| 亚洲精品国产精品久久久不卡| 可以在线观看的亚洲视频| 国产亚洲精品av在线| 欧美日本亚洲视频在线播放| 久久久久久国产a免费观看| a级毛片在线看网站| 亚洲第一欧美日韩一区二区三区| 搡老熟女国产l中国老女人| 国产精品久久久久久久电影 | 在线a可以看的网站| 无遮挡黄片免费观看| 亚洲精品粉嫩美女一区| 国模一区二区三区四区视频 | 好男人电影高清在线观看| 97超视频在线观看视频| 亚洲黑人精品在线| 日本黄色片子视频| 怎么达到女性高潮| 亚洲精品在线美女| 国产精品久久久人人做人人爽| 欧美成人性av电影在线观看| 国产日本99.免费观看| 天堂网av新在线| 一级黄色大片毛片| 色噜噜av男人的天堂激情| 欧美性猛交╳xxx乱大交人| 美女cb高潮喷水在线观看 | www.自偷自拍.com| 美女高潮喷水抽搐中文字幕| 高清毛片免费观看视频网站| 国产精品永久免费网站| 国产亚洲精品久久久久久毛片| 热99在线观看视频| 91在线精品国自产拍蜜月 | 禁无遮挡网站| 午夜久久久久精精品| 日本熟妇午夜| 午夜视频精品福利| 岛国视频午夜一区免费看| 嫩草影院入口| 欧美一级毛片孕妇| 国产亚洲精品一区二区www| 成人高潮视频无遮挡免费网站| 精品福利观看| 亚洲国产精品999在线| 亚洲狠狠婷婷综合久久图片| 欧美日韩一级在线毛片| 黑人欧美特级aaaaaa片| 又紧又爽又黄一区二区| 午夜福利欧美成人| 我要搜黄色片| xxxwww97欧美| 小蜜桃在线观看免费完整版高清| 夜夜看夜夜爽夜夜摸| 国产精品,欧美在线| 国产精品 国内视频| 欧美极品一区二区三区四区| 在线国产一区二区在线| 午夜福利高清视频| 哪里可以看免费的av片| 给我免费播放毛片高清在线观看| 中文在线观看免费www的网站| av天堂中文字幕网| 国产高清激情床上av| 精品免费久久久久久久清纯| 夜夜夜夜夜久久久久| 一级毛片高清免费大全| ponron亚洲| 床上黄色一级片| 男人和女人高潮做爰伦理| 最新美女视频免费是黄的| 日本成人三级电影网站| 九九久久精品国产亚洲av麻豆 | 久久午夜亚洲精品久久| 久久热在线av| 不卡一级毛片| 90打野战视频偷拍视频| 亚洲欧美激情综合另类| 黑人操中国人逼视频| 91字幕亚洲| 亚洲av电影在线进入| 很黄的视频免费| 成人高潮视频无遮挡免费网站| 黄片小视频在线播放| 国产成人一区二区三区免费视频网站| 久久久久国产精品人妻aⅴ院| 香蕉av资源在线| 国产精华一区二区三区| 久久久精品大字幕| 国内精品美女久久久久久| 黄片小视频在线播放| 十八禁人妻一区二区| 亚洲在线观看片| 我的老师免费观看完整版| 1024香蕉在线观看| 国模一区二区三区四区视频 | 亚洲中文字幕日韩| 色吧在线观看| 久久九九热精品免费| 99久久精品一区二区三区| 午夜免费成人在线视频| 99国产精品一区二区蜜桃av| www.精华液| 三级男女做爰猛烈吃奶摸视频| 一a级毛片在线观看| 色播亚洲综合网| 熟女少妇亚洲综合色aaa.| 久久香蕉精品热| 成人特级av手机在线观看| 成年免费大片在线观看| 超碰成人久久| 午夜免费观看网址| 欧美色欧美亚洲另类二区| 欧美高清成人免费视频www| 亚洲成人精品中文字幕电影| 欧美zozozo另类| 亚洲在线观看片| 婷婷精品国产亚洲av在线| 一夜夜www| 老司机福利观看| 免费无遮挡裸体视频| 精品乱码久久久久久99久播| 亚洲人成电影免费在线| 黑人操中国人逼视频| 日韩欧美精品v在线| 国产高清有码在线观看视频| 中出人妻视频一区二区| 国产视频一区二区在线看| 19禁男女啪啪无遮挡网站| 色视频www国产| 亚洲欧美精品综合一区二区三区| 色吧在线观看| 国产精品美女特级片免费视频播放器 | 一级毛片精品| 男女午夜视频在线观看| cao死你这个sao货| 免费观看精品视频网站| 久久久久亚洲av毛片大全| 亚洲七黄色美女视频| 特大巨黑吊av在线直播| 国产精品1区2区在线观看.| 亚洲精品色激情综合| 国产成人aa在线观看| 黑人操中国人逼视频| 久久亚洲精品不卡| 给我免费播放毛片高清在线观看| 天堂av国产一区二区熟女人妻| 亚洲美女黄片视频| 欧美zozozo另类| 日本一本二区三区精品| 国产爱豆传媒在线观看| 免费观看人在逋| 久久午夜综合久久蜜桃| 亚洲 国产 在线| 日本a在线网址| 久久国产精品人妻蜜桃| 成人国产一区最新在线观看| 国产成人精品无人区| 国产在线精品亚洲第一网站| 一级a爱片免费观看的视频| 久久久久九九精品影院| 一级毛片女人18水好多| 国产精品久久久av美女十八| 国产精华一区二区三区| 好看av亚洲va欧美ⅴa在| 国产91精品成人一区二区三区| 久久久久性生活片| 午夜两性在线视频| 亚洲性夜色夜夜综合| 级片在线观看| 亚洲真实伦在线观看| 午夜免费激情av| 国产亚洲精品综合一区在线观看| 午夜福利在线观看免费完整高清在 | 欧美一区二区精品小视频在线| 18禁美女被吸乳视频| 免费在线观看日本一区| 国产精品久久久av美女十八| 精品一区二区三区视频在线 | 美女被艹到高潮喷水动态| 少妇人妻一区二区三区视频| 十八禁人妻一区二区| 看黄色毛片网站| 亚洲国产看品久久| 亚洲国产色片| 在线十欧美十亚洲十日本专区| 99精品欧美一区二区三区四区| 成年免费大片在线观看| 亚洲中文字幕日韩| 久久久精品欧美日韩精品| 伦理电影免费视频| 亚洲精品在线观看二区| 精品久久久久久久末码| 欧美高清成人免费视频www| 久久精品国产亚洲av香蕉五月| 青草久久国产| 亚洲专区字幕在线| 日本精品一区二区三区蜜桃| 欧美中文日本在线观看视频| 亚洲精品在线美女| 国产伦人伦偷精品视频| 亚洲在线自拍视频| 少妇人妻一区二区三区视频| 午夜精品在线福利| 日韩国内少妇激情av| 哪里可以看免费的av片| 欧美日韩亚洲国产一区二区在线观看| 免费在线观看影片大全网站| 性色av乱码一区二区三区2| 99热精品在线国产| 色噜噜av男人的天堂激情| 亚洲av第一区精品v没综合| 欧美日韩亚洲国产一区二区在线观看| av中文乱码字幕在线| 哪里可以看免费的av片| 人妻久久中文字幕网| 1000部很黄的大片| 一个人观看的视频www高清免费观看 | 久99久视频精品免费| 国产精品98久久久久久宅男小说| 亚洲av日韩精品久久久久久密| 日韩欧美国产在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 麻豆成人午夜福利视频| 婷婷精品国产亚洲av| 亚洲av电影在线进入| 国产在线精品亚洲第一网站| 亚洲国产精品sss在线观看| 欧美乱妇无乱码| 岛国在线观看网站| 久久久久久久午夜电影| 精品久久久久久久久久久久久| 午夜福利18| 成人av在线播放网站| 国产精品一区二区精品视频观看| 淫秽高清视频在线观看| 九色国产91popny在线| 搡老岳熟女国产| 午夜激情欧美在线| 精品一区二区三区视频在线观看免费| 两个人视频免费观看高清| 午夜日韩欧美国产| 天堂动漫精品| 国产精品 欧美亚洲| 亚洲精品一区av在线观看| av女优亚洲男人天堂 | 亚洲人成伊人成综合网2020| 日本黄色视频三级网站网址| 亚洲aⅴ乱码一区二区在线播放| netflix在线观看网站| 免费高清视频大片| 亚洲第一电影网av| 日韩成人在线观看一区二区三区| 亚洲人成电影免费在线| 国产99白浆流出| 精品久久久久久久末码| 757午夜福利合集在线观看| 国产一区二区三区在线臀色熟女| 亚洲在线自拍视频| 男女做爰动态图高潮gif福利片| 最新在线观看一区二区三区| 久久久久久九九精品二区国产| 久久久久九九精品影院| 成年版毛片免费区| 一夜夜www| 黑人操中国人逼视频| 国产精品,欧美在线| 蜜桃久久精品国产亚洲av| 男女午夜视频在线观看| 日韩有码中文字幕| 精品国产超薄肉色丝袜足j| 母亲3免费完整高清在线观看| 天天一区二区日本电影三级| 在线观看美女被高潮喷水网站 | 人妻夜夜爽99麻豆av| 黄片大片在线免费观看| av福利片在线观看| 精品免费久久久久久久清纯| 91九色精品人成在线观看| 琪琪午夜伦伦电影理论片6080| 亚洲欧美日韩高清在线视频| 黄片小视频在线播放| 级片在线观看| 九九久久精品国产亚洲av麻豆 | 成人精品一区二区免费| 天堂动漫精品| 中文在线观看免费www的网站| 国产精品av视频在线免费观看| 一卡2卡三卡四卡精品乱码亚洲| 精品久久久久久,| 桃色一区二区三区在线观看| 成人国产一区最新在线观看| 麻豆国产av国片精品| 精品一区二区三区视频在线观看免费| 国产精品久久视频播放| 网址你懂的国产日韩在线| 99国产精品一区二区蜜桃av| 一个人看视频在线观看www免费 | 九色成人免费人妻av| 啦啦啦观看免费观看视频高清| 日本精品一区二区三区蜜桃| 久久人妻av系列| 国产91精品成人一区二区三区| 午夜久久久久精精品| 丰满的人妻完整版| 日韩欧美在线二视频| 精品国产三级普通话版| 亚洲男人的天堂狠狠| 男人舔女人下体高潮全视频| 色综合欧美亚洲国产小说| 亚洲成人中文字幕在线播放| 亚洲第一欧美日韩一区二区三区| 亚洲自拍偷在线| 久久久国产成人精品二区| 色在线成人网| 成年女人毛片免费观看观看9| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品亚洲av一区麻豆| 亚洲自拍偷在线| www国产在线视频色| 搡老岳熟女国产| 欧美日韩亚洲国产一区二区在线观看| 欧美在线黄色| 蜜桃久久精品国产亚洲av| 久久精品aⅴ一区二区三区四区| 在线观看免费视频日本深夜| 男女视频在线观看网站免费| 日本成人三级电影网站| 丝袜人妻中文字幕| 中文亚洲av片在线观看爽| a在线观看视频网站| 成人特级黄色片久久久久久久| 禁无遮挡网站| 免费在线观看影片大全网站| 国产成人aa在线观看| 国产av在哪里看| 国产精品亚洲av一区麻豆| 亚洲成人中文字幕在线播放|