• <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型鋼支架卡纜臨界約束力理論研究
    亚洲av国产av综合av卡| 日韩大码丰满熟妇| 视频区欧美日本亚洲| 成人免费观看视频高清| 国产精品自产拍在线观看55亚洲 | 亚洲情色 制服丝袜| 可以免费在线观看a视频的电影网站| 久久香蕉激情| 久9热在线精品视频| 精品福利永久在线观看| 亚洲一区二区三区欧美精品| 69av精品久久久久久 | 男女之事视频高清在线观看| 成人亚洲精品一区在线观看| 视频区欧美日本亚洲| 叶爱在线成人免费视频播放| 日韩熟女老妇一区二区性免费视频| 91成年电影在线观看| 成人亚洲精品一区在线观看| 香蕉丝袜av| 51午夜福利影视在线观看| 国产精品久久久久久人妻精品电影 | 久久久久久久大尺度免费视频| 成人av一区二区三区在线看 | 国产黄频视频在线观看| 欧美精品人与动牲交sv欧美| 精品福利永久在线观看| 亚洲国产欧美一区二区综合| 午夜福利在线观看吧| 国产黄频视频在线观看| 国产av又大| 99精品久久久久人妻精品| 国产麻豆69| 亚洲人成电影免费在线| 最黄视频免费看| 99热国产这里只有精品6| 精品国产一区二区久久| 欧美一级毛片孕妇| 欧美少妇被猛烈插入视频| 青春草亚洲视频在线观看| 高清在线国产一区| 狂野欧美激情性xxxx| 十分钟在线观看高清视频www| 国产一区二区在线观看av| 成人国语在线视频| 精品久久蜜臀av无| 最新在线观看一区二区三区| 各种免费的搞黄视频| 久久久久视频综合| 日韩中文字幕视频在线看片| 伦理电影免费视频| 日日夜夜操网爽| 91九色精品人成在线观看| 欧美 亚洲 国产 日韩一| 秋霞在线观看毛片| 欧美日韩精品网址| 亚洲成人免费电影在线观看| 国产深夜福利视频在线观看| 精品国产一区二区三区四区第35| 国产av一区二区精品久久| 欧美变态另类bdsm刘玥| 成人手机av| 十八禁人妻一区二区| 纵有疾风起免费观看全集完整版| 老鸭窝网址在线观看| 一区二区三区精品91| 国产欧美日韩综合在线一区二区| 国产成人免费观看mmmm| 国产成人欧美在线观看 | 日本91视频免费播放| 巨乳人妻的诱惑在线观看| 久久久欧美国产精品| 亚洲中文av在线| 免费观看av网站的网址| 中国美女看黄片| 久久精品久久久久久噜噜老黄| 少妇的丰满在线观看| 青春草亚洲视频在线观看| 亚洲中文日韩欧美视频| 亚洲精品av麻豆狂野| 亚洲av片天天在线观看| 久久国产亚洲av麻豆专区| 99久久综合免费| 不卡av一区二区三区| 欧美成人午夜精品| 丝袜脚勾引网站| 老熟妇仑乱视频hdxx| av天堂在线播放| 欧美精品高潮呻吟av久久| 久久国产精品男人的天堂亚洲| 精品亚洲成国产av| 天天躁狠狠躁夜夜躁狠狠躁| 精品亚洲成a人片在线观看| 日本欧美视频一区| 考比视频在线观看| 99精品欧美一区二区三区四区| 人妻一区二区av| 深夜精品福利| 美女高潮到喷水免费观看| videosex国产| 午夜福利免费观看在线| 亚洲国产精品一区二区三区在线| 男女边摸边吃奶| 国产黄色免费在线视频| 这个男人来自地球电影免费观看| 国产熟女午夜一区二区三区| 99热国产这里只有精品6| 两人在一起打扑克的视频| 十分钟在线观看高清视频www| av一本久久久久| 日日爽夜夜爽网站| 亚洲专区字幕在线| 亚洲一区中文字幕在线| 9热在线视频观看99| 一区二区三区激情视频| 亚洲精品国产av蜜桃| 黑人巨大精品欧美一区二区蜜桃| 国产一区有黄有色的免费视频| h视频一区二区三区| 亚洲久久久国产精品| 欧美激情高清一区二区三区| 交换朋友夫妻互换小说| 亚洲av国产av综合av卡| 亚洲天堂av无毛| 日韩一卡2卡3卡4卡2021年| av又黄又爽大尺度在线免费看| 国产亚洲午夜精品一区二区久久| 久久国产精品人妻蜜桃| 91精品国产国语对白视频| 欧美激情久久久久久爽电影 | 丁香六月欧美| 少妇粗大呻吟视频| 男女国产视频网站| 久久精品亚洲av国产电影网| 天天添夜夜摸| 国产欧美亚洲国产| 91大片在线观看| 三级毛片av免费| 亚洲人成电影免费在线| 免费看十八禁软件| 狠狠狠狠99中文字幕| 国产精品久久久人人做人人爽| 一级黄色大片毛片| 正在播放国产对白刺激| 首页视频小说图片口味搜索| av超薄肉色丝袜交足视频| 欧美激情 高清一区二区三区| 精品一区二区三区av网在线观看 | 日韩大片免费观看网站| 久久精品国产亚洲av香蕉五月 | 水蜜桃什么品种好| 老司机福利观看| 欧美成狂野欧美在线观看| 在线观看免费视频网站a站| 啦啦啦视频在线资源免费观看| 超色免费av| 日韩三级视频一区二区三区| 亚洲少妇的诱惑av| 日日爽夜夜爽网站| 国产成人免费观看mmmm| 视频在线观看一区二区三区| 岛国毛片在线播放| www.自偷自拍.com| 美女国产高潮福利片在线看| 亚洲欧美一区二区三区久久| 亚洲成人手机| 我要看黄色一级片免费的| 老司机影院成人| 久久国产精品影院| 动漫黄色视频在线观看| 自线自在国产av| 久久久久久久大尺度免费视频| 国产高清视频在线播放一区 | 久久天堂一区二区三区四区| 极品少妇高潮喷水抽搐| 久久人人97超碰香蕉20202| 国产成人欧美| 纯流量卡能插随身wifi吗| 亚洲欧美激情在线| 亚洲三区欧美一区| 天天躁夜夜躁狠狠躁躁| 亚洲第一av免费看| 波多野结衣av一区二区av| 精品少妇一区二区三区视频日本电影| 久久久久久久精品精品| 精品一区在线观看国产| 两个人看的免费小视频| 丰满迷人的少妇在线观看| 别揉我奶头~嗯~啊~动态视频 | 中亚洲国语对白在线视频| 伊人亚洲综合成人网| 岛国毛片在线播放| www.av在线官网国产| 美女中出高潮动态图| 老熟女久久久| 成人av一区二区三区在线看 | 欧美大码av| 伦理电影免费视频| 亚洲自偷自拍图片 自拍| 老司机靠b影院| 久久精品aⅴ一区二区三区四区| 视频区欧美日本亚洲| 精品久久久久久电影网| 亚洲av成人一区二区三| 水蜜桃什么品种好| 亚洲全国av大片| 人妻 亚洲 视频| 两个人免费观看高清视频| 午夜免费观看性视频| 视频区欧美日本亚洲| cao死你这个sao货| 国产黄色免费在线视频| 一区二区三区精品91| 亚洲欧美成人综合另类久久久| 爱豆传媒免费全集在线观看| 伦理电影免费视频| 91精品伊人久久大香线蕉| 免费日韩欧美在线观看| 久久人人爽人人片av| 韩国精品一区二区三区| h视频一区二区三区| 国产精品九九99| 51午夜福利影视在线观看| 久久精品成人免费网站| 久久九九热精品免费| 精品国产乱码久久久久久男人| 国产成人精品无人区| 后天国语完整版免费观看| 狂野欧美激情性xxxx| 久久综合国产亚洲精品| av不卡在线播放| 亚洲熟女毛片儿| 精品人妻1区二区| 三上悠亚av全集在线观看| 日韩欧美国产一区二区入口| 丁香六月天网| 日韩 欧美 亚洲 中文字幕| 老汉色av国产亚洲站长工具| 国产精品自产拍在线观看55亚洲 | 精品久久蜜臀av无| 大片电影免费在线观看免费| 精品国内亚洲2022精品成人 | 热re99久久精品国产66热6| 悠悠久久av| 青春草视频在线免费观看| 人人妻人人爽人人添夜夜欢视频| 国产精品一区二区在线观看99| 久久免费观看电影| 日本wwww免费看| 成人手机av| 色94色欧美一区二区| 亚洲,欧美精品.| videosex国产| 1024视频免费在线观看| 国产av又大| 在线av久久热| 久久久久久久久免费视频了| 丝袜脚勾引网站| 中文字幕制服av| 欧美日本中文国产一区发布| 老鸭窝网址在线观看| 欧美另类一区| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲成av片中文字幕在线观看| 男女免费视频国产| 在线av久久热| 王馨瑶露胸无遮挡在线观看| 久久久久国产精品人妻一区二区| 一个人免费看片子| 80岁老熟妇乱子伦牲交| 人人妻人人澡人人爽人人夜夜| 捣出白浆h1v1| 黄色视频,在线免费观看| 麻豆av在线久日| 最黄视频免费看| 中国国产av一级| 亚洲 欧美一区二区三区| av一本久久久久| 国产精品一区二区在线不卡| 青青草视频在线视频观看| 精品国产乱子伦一区二区三区 | 亚洲欧美精品综合一区二区三区| 久久久久久免费高清国产稀缺| 国产一级毛片在线| 欧美在线黄色| 在线观看免费视频网站a站| 亚洲精品美女久久av网站| 欧美日本中文国产一区发布| 99精国产麻豆久久婷婷| 日韩人妻精品一区2区三区| 国产极品粉嫩免费观看在线| 国产片内射在线| 丝袜美足系列| 亚洲精品一区蜜桃| 国产欧美日韩一区二区三 | 男女床上黄色一级片免费看| 日韩,欧美,国产一区二区三区| 91av网站免费观看| 亚洲av电影在线进入| 1024视频免费在线观看| 视频区图区小说| 婷婷丁香在线五月| 国产成人一区二区三区免费视频网站| 麻豆av在线久日| 欧美日韩一级在线毛片| 欧美国产精品一级二级三级| 亚洲精品国产av成人精品| 曰老女人黄片| 狠狠精品人妻久久久久久综合| 日韩视频在线欧美| 天天操日日干夜夜撸| 成人18禁高潮啪啪吃奶动态图| 成年人午夜在线观看视频| 久久狼人影院| 美女大奶头黄色视频| 女人精品久久久久毛片| 亚洲精品国产av成人精品| 亚洲成人免费av在线播放| 亚洲精品中文字幕一二三四区 | 中文欧美无线码| 亚洲欧洲日产国产| 国产亚洲av高清不卡| 国产淫语在线视频| 男女之事视频高清在线观看| 亚洲成国产人片在线观看| 日日摸夜夜添夜夜添小说| 亚洲欧美激情在线| 在线观看人妻少妇| 日韩熟女老妇一区二区性免费视频| 亚洲视频免费观看视频| 纵有疾风起免费观看全集完整版| 老司机午夜十八禁免费视频| 亚洲综合色网址| 国精品久久久久久国模美| 黑人操中国人逼视频| 超碰97精品在线观看| 亚洲国产精品999| 如日韩欧美国产精品一区二区三区| 一进一出抽搐动态| 美女扒开内裤让男人捅视频| 国产在视频线精品| videos熟女内射| 99精国产麻豆久久婷婷| 亚洲伊人久久精品综合| av不卡在线播放| 少妇精品久久久久久久| 亚洲七黄色美女视频| 各种免费的搞黄视频| 一本综合久久免费| 狠狠精品人妻久久久久久综合| 亚洲第一av免费看| 黄色视频不卡| 99国产精品一区二区蜜桃av | 成人影院久久| 男女下面插进去视频免费观看| 女警被强在线播放| 91精品国产国语对白视频| 国产在线一区二区三区精| 90打野战视频偷拍视频| 视频区图区小说| 咕卡用的链子| 美女午夜性视频免费| 欧美精品人与动牲交sv欧美| 国产日韩一区二区三区精品不卡| 悠悠久久av| www.av在线官网国产| e午夜精品久久久久久久| 久久精品亚洲熟妇少妇任你| 免费在线观看完整版高清| 99香蕉大伊视频| netflix在线观看网站| 97在线人人人人妻| 最近最新中文字幕大全免费视频| 法律面前人人平等表现在哪些方面 | 国产高清videossex| 少妇 在线观看| 18禁国产床啪视频网站| 亚洲精品久久久久久婷婷小说| av在线播放精品| 999久久久精品免费观看国产| 1024香蕉在线观看| 免费人妻精品一区二区三区视频| 亚洲国产欧美在线一区| 脱女人内裤的视频| 免费在线观看影片大全网站| 黄片小视频在线播放| 天天操日日干夜夜撸| 十分钟在线观看高清视频www| 咕卡用的链子| 丝袜喷水一区| 亚洲九九香蕉| 国产精品秋霞免费鲁丝片| 亚洲精品久久午夜乱码| 国产老妇伦熟女老妇高清| 精品亚洲乱码少妇综合久久| 美女脱内裤让男人舔精品视频| 日韩制服骚丝袜av| 久久精品亚洲熟妇少妇任你| 久久精品国产亚洲av香蕉五月 | 精品久久久精品久久久| 精品国产一区二区三区四区第35| 又大又爽又粗| av在线播放精品| 精品人妻1区二区| 最近中文字幕2019免费版| 18禁裸乳无遮挡动漫免费视频| 亚洲精品av麻豆狂野| 国产亚洲av高清不卡| 日韩电影二区| 麻豆av在线久日| 久久午夜综合久久蜜桃| 久热爱精品视频在线9| 久久av网站| 久久性视频一级片| 亚洲色图综合在线观看| 午夜两性在线视频| 欧美黄色片欧美黄色片| 如日韩欧美国产精品一区二区三区| 多毛熟女@视频| 午夜两性在线视频| 在线观看人妻少妇| 91成年电影在线观看| 爱豆传媒免费全集在线观看| 中文字幕精品免费在线观看视频| 丝袜喷水一区| 成在线人永久免费视频| 国产熟女午夜一区二区三区| 99九九在线精品视频| 精品乱码久久久久久99久播| 国产精品偷伦视频观看了| 欧美少妇被猛烈插入视频| 亚洲熟女毛片儿| 国产一区二区三区在线臀色熟女 | 法律面前人人平等表现在哪些方面 | 久久av网站| 在线av久久热| 在线观看www视频免费| 叶爱在线成人免费视频播放| 法律面前人人平等表现在哪些方面 | 男人添女人高潮全过程视频| 一区二区三区精品91| 亚洲精品国产av成人精品| 欧美黄色片欧美黄色片| 两人在一起打扑克的视频| 国产一级毛片在线| 高清黄色对白视频在线免费看| 男女之事视频高清在线观看| 19禁男女啪啪无遮挡网站| 精品久久久久久电影网| 九色亚洲精品在线播放| 成年女人毛片免费观看观看9 | 日韩中文字幕欧美一区二区| 99精品欧美一区二区三区四区| 久久精品久久久久久噜噜老黄| 真人做人爱边吃奶动态| 在线观看一区二区三区激情| 久久性视频一级片| 9191精品国产免费久久| 在线永久观看黄色视频| 国产精品秋霞免费鲁丝片| 一二三四在线观看免费中文在| 美女高潮喷水抽搐中文字幕| 香蕉国产在线看| 国产成人精品久久二区二区免费| 考比视频在线观看| 亚洲熟女精品中文字幕| a级毛片黄视频| 捣出白浆h1v1| av又黄又爽大尺度在线免费看| 国产精品久久久久久精品古装| 高清黄色对白视频在线免费看| 欧美大码av| 日韩欧美国产一区二区入口| 国产黄色免费在线视频| 久久精品国产综合久久久| 考比视频在线观看| 午夜老司机福利片| 一级黄色大片毛片| 最近中文字幕2019免费版| 国产人伦9x9x在线观看| 国产又爽黄色视频| 亚洲国产av影院在线观看| 在线亚洲精品国产二区图片欧美| 桃红色精品国产亚洲av| 国产精品免费视频内射| 欧美在线一区亚洲| 亚洲成人手机| 国产一区二区激情短视频 | 国产视频一区二区在线看| 视频区欧美日本亚洲| 欧美精品av麻豆av| 国产又爽黄色视频| 国产精品秋霞免费鲁丝片| 亚洲五月色婷婷综合| 日韩欧美免费精品| 国产日韩一区二区三区精品不卡| 精品人妻一区二区三区麻豆| 一区二区三区精品91| 久久精品国产综合久久久| 久久久欧美国产精品| 精品久久久久久电影网| 麻豆乱淫一区二区| 少妇裸体淫交视频免费看高清 | 精品一品国产午夜福利视频| 高清在线国产一区| 国产高清视频在线播放一区 | 高潮久久久久久久久久久不卡| 国产真人三级小视频在线观看| 国产精品成人在线| 不卡av一区二区三区| 久久久水蜜桃国产精品网| 国产亚洲欧美精品永久| 日本91视频免费播放| 日韩制服骚丝袜av| 黑人巨大精品欧美一区二区蜜桃| 欧美日韩视频精品一区| 一级黄色大片毛片| 91精品三级在线观看| av国产精品久久久久影院| 国产精品二区激情视频| 国产成人啪精品午夜网站| www.自偷自拍.com| 成人手机av| 麻豆乱淫一区二区| 国产精品 国内视频| 五月天丁香电影| 男女之事视频高清在线观看| 大片免费播放器 马上看| 女性生殖器流出的白浆| 在线观看免费高清a一片| kizo精华| svipshipincom国产片| 巨乳人妻的诱惑在线观看| 亚洲免费av在线视频| 精品人妻1区二区| 一进一出抽搐动态| 欧美激情久久久久久爽电影 | av福利片在线| 妹子高潮喷水视频| 超色免费av| 麻豆乱淫一区二区| 十八禁人妻一区二区| a 毛片基地| 日本a在线网址| 99九九在线精品视频| 午夜福利在线免费观看网站| 美女福利国产在线| 国产精品九九99| 亚洲精品av麻豆狂野| 成年人免费黄色播放视频| 国产成人av激情在线播放| 国产不卡av网站在线观看| www.精华液| 青春草视频在线免费观看| 欧美另类亚洲清纯唯美| 亚洲第一av免费看| 少妇的丰满在线观看| 亚洲精品一区蜜桃| 精品国产一区二区久久| 欧美日韩黄片免| 香蕉国产在线看| 999精品在线视频| 99热国产这里只有精品6| av在线app专区| 欧美黄色淫秽网站| 搡老岳熟女国产| 国产精品欧美亚洲77777| 少妇被粗大的猛进出69影院| svipshipincom国产片| 国产真人三级小视频在线观看| 精品熟女少妇八av免费久了| 另类精品久久| 国产精品成人在线| 建设人人有责人人尽责人人享有的| 国产精品 国内视频| 老司机影院成人| 欧美xxⅹ黑人| 岛国在线观看网站| 欧美在线一区亚洲| 丝袜脚勾引网站| 视频在线观看一区二区三区| 激情视频va一区二区三区| 91国产中文字幕| 国产精品久久久久久精品古装| 亚洲色图 男人天堂 中文字幕| 亚洲三区欧美一区| 欧美精品一区二区免费开放| 狂野欧美激情性bbbbbb| 伊人亚洲综合成人网| 777久久人妻少妇嫩草av网站| 满18在线观看网站| 成年人午夜在线观看视频| 亚洲欧美一区二区三区黑人| 热99国产精品久久久久久7| 狠狠婷婷综合久久久久久88av| netflix在线观看网站| 91av网站免费观看| 一边摸一边抽搐一进一出视频| 电影成人av| 丝袜在线中文字幕| 黄色片一级片一级黄色片| av天堂在线播放| av片东京热男人的天堂| tube8黄色片| 亚洲精品国产一区二区精华液| 国产精品香港三级国产av潘金莲| √禁漫天堂资源中文www| 亚洲,欧美精品.| 中文字幕av电影在线播放| 色94色欧美一区二区| 精品久久久精品久久久| 日本a在线网址| 亚洲欧美精品综合一区二区三区| 黄片大片在线免费观看| 狠狠精品人妻久久久久久综合|