• <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网站免费在线观看视频| 亚洲精品日本国产第一区| 亚洲精品日本国产第一区| 乱码一卡2卡4卡精品| 亚洲国产日韩一区二区| 国产一区二区三区av在线| 国产亚洲欧美精品永久| a 毛片基地| 精品亚洲成国产av| 成人亚洲欧美一区二区av| av福利片在线| 免费观看性生交大片5| 亚洲人与动物交配视频| 日韩欧美一区视频在线观看| a级毛片在线看网站| 亚洲国产欧美在线一区| 一本久久精品| 一个人免费看片子| 寂寞人妻少妇视频99o| 99热网站在线观看| 午夜91福利影院| 亚洲精品456在线播放app| a级片在线免费高清观看视频| 伊人亚洲综合成人网| 国产一区亚洲一区在线观看| 中文乱码字字幕精品一区二区三区| 高清在线视频一区二区三区| 99国产综合亚洲精品| 精品久久久久久电影网| 日韩欧美精品免费久久| 国产成人精品一,二区| 熟女电影av网| 18禁裸乳无遮挡动漫免费视频| 啦啦啦中文免费视频观看日本| 在现免费观看毛片| 97精品久久久久久久久久精品| 曰老女人黄片| 在线 av 中文字幕| 亚洲欧美中文字幕日韩二区| 丝袜在线中文字幕| 国产精品.久久久| 美女主播在线视频| 午夜91福利影院| 精品人妻一区二区三区麻豆| 黑丝袜美女国产一区| 日韩不卡一区二区三区视频在线| 一区二区三区免费毛片| av黄色大香蕉| 午夜福利在线观看免费完整高清在| 久久久久视频综合| 欧美日韩av久久| 一级毛片aaaaaa免费看小| 国产免费一级a男人的天堂| 国产欧美另类精品又又久久亚洲欧美| 简卡轻食公司| 人人妻人人爽人人添夜夜欢视频| 亚洲美女黄色视频免费看| 黑人欧美特级aaaaaa片| 91精品国产九色| 成年人免费黄色播放视频| 亚洲精品日韩av片在线观看| 亚洲国产最新在线播放| 黄色欧美视频在线观看| 汤姆久久久久久久影院中文字幕| 亚洲精品国产色婷婷电影| 在线免费观看不下载黄p国产| 久久久精品免费免费高清| 伊人亚洲综合成人网| 欧美少妇被猛烈插入视频| 欧美xxxx性猛交bbbb| 久热这里只有精品99| 亚洲av欧美aⅴ国产| 有码 亚洲区| 水蜜桃什么品种好| 91精品三级在线观看| 精品久久久久久久久亚洲| 久久精品久久精品一区二区三区| 国产乱人偷精品视频| 老司机影院毛片| 久久久久久伊人网av| 能在线免费看毛片的网站| 久久久国产欧美日韩av| 伊人亚洲综合成人网| 亚洲精品国产色婷婷电影| 青春草国产在线视频| 人妻一区二区av| 亚洲熟女精品中文字幕| 午夜91福利影院| 欧美精品高潮呻吟av久久| 人妻夜夜爽99麻豆av| 日本猛色少妇xxxxx猛交久久| 丝袜脚勾引网站| 亚洲人成网站在线观看播放| 日本午夜av视频| 久久 成人 亚洲| 欧美成人午夜免费资源| 午夜福利网站1000一区二区三区| 美女xxoo啪啪120秒动态图| 国产精品久久久久久精品古装| 免费大片黄手机在线观看| 日韩伦理黄色片| 精品少妇内射三级| 乱码一卡2卡4卡精品| 国产成人精品久久久久久| 国产亚洲欧美精品永久| 一边摸一边做爽爽视频免费| 国产在线免费精品| 少妇人妻 视频| 午夜老司机福利剧场| 午夜影院在线不卡| 免费人妻精品一区二区三区视频| 黄片播放在线免费| 涩涩av久久男人的天堂| 亚洲精品久久久久久婷婷小说| 9色porny在线观看| 精品亚洲乱码少妇综合久久| 22中文网久久字幕| 国产高清有码在线观看视频| 免费少妇av软件| 日韩av在线免费看完整版不卡| 美女xxoo啪啪120秒动态图| 亚洲国产成人一精品久久久| 免费观看av网站的网址| av.在线天堂| 国产精品免费大片| 香蕉精品网在线| 黑丝袜美女国产一区| 欧美精品国产亚洲| 国产成人aa在线观看| 欧美老熟妇乱子伦牲交| 日韩制服骚丝袜av| 少妇 在线观看| 欧美人与善性xxx| 26uuu在线亚洲综合色| videossex国产| 久久久久久久亚洲中文字幕| 汤姆久久久久久久影院中文字幕| 少妇猛男粗大的猛烈进出视频| 一区二区三区免费毛片| 少妇被粗大的猛进出69影院 | 少妇精品久久久久久久| 免费高清在线观看视频在线观看| 性色avwww在线观看| 亚洲三级黄色毛片| 欧美日韩精品成人综合77777| 好男人视频免费观看在线| 天天影视国产精品| 哪个播放器可以免费观看大片| 99九九在线精品视频| 一级毛片aaaaaa免费看小| 只有这里有精品99| 精品卡一卡二卡四卡免费| 成人综合一区亚洲| av福利片在线| 51国产日韩欧美| 久久精品国产亚洲av涩爱| 国产伦精品一区二区三区视频9| 国产亚洲精品久久久com| 国产不卡av网站在线观看| av福利片在线| 最近的中文字幕免费完整| 精品少妇黑人巨大在线播放| 精品午夜福利在线看| 99热这里只有是精品在线观看| 亚洲中文av在线| 女人精品久久久久毛片| 韩国高清视频一区二区三区| 欧美xxⅹ黑人| 日日爽夜夜爽网站| 两个人免费观看高清视频| 看非洲黑人一级黄片| 亚洲综合精品二区| 国产探花极品一区二区| 2018国产大陆天天弄谢| 国产精品蜜桃在线观看| 国产精品 国内视频| 国产黄片视频在线免费观看| 久久狼人影院| 亚洲国产av影院在线观看| 老司机影院成人| 国产伦理片在线播放av一区| 国产男人的电影天堂91| 久久精品国产亚洲av天美| 亚洲av日韩在线播放| 国产片内射在线| 国产成人aa在线观看| 国产精品麻豆人妻色哟哟久久| 十八禁网站网址无遮挡| 我的老师免费观看完整版| 久久久久精品久久久久真实原创| 在线观看免费高清a一片| 亚洲成人一二三区av| 欧美老熟妇乱子伦牲交| 汤姆久久久久久久影院中文字幕| 18禁裸乳无遮挡动漫免费视频| 日本欧美国产在线视频| 一个人看视频在线观看www免费| 成人国语在线视频| 狠狠精品人妻久久久久久综合| 亚洲第一区二区三区不卡| 青春草亚洲视频在线观看| 黑人巨大精品欧美一区二区蜜桃 | 自拍欧美九色日韩亚洲蝌蚪91| 韩国高清视频一区二区三区| 两个人的视频大全免费| 高清不卡的av网站| 国产午夜精品一二区理论片| 一本—道久久a久久精品蜜桃钙片| 青春草国产在线视频| 精品国产露脸久久av麻豆| 青春草亚洲视频在线观看| av在线app专区| 亚洲av中文av极速乱| 久久午夜福利片| 18+在线观看网站| 草草在线视频免费看| 水蜜桃什么品种好| 日本-黄色视频高清免费观看| 黄片播放在线免费| 国产成人精品一,二区| 日本色播在线视频| 99热全是精品| 人成视频在线观看免费观看| 亚洲国产最新在线播放| 黄色一级大片看看| 成人漫画全彩无遮挡| 91成人精品电影| 免费人妻精品一区二区三区视频| 亚洲四区av| 精品人妻熟女毛片av久久网站| 黑人高潮一二区| 嘟嘟电影网在线观看| 两个人的视频大全免费| 婷婷色综合www| 91精品三级在线观看| 亚洲综合色网址| 中文字幕亚洲精品专区| 男的添女的下面高潮视频| 亚洲成人手机| 免费看光身美女| 亚洲内射少妇av| 国产爽快片一区二区三区| 亚洲精品自拍成人| 曰老女人黄片| 成年美女黄网站色视频大全免费 | 在线天堂最新版资源| 人人澡人人妻人| 在线观看一区二区三区激情| 国产精品秋霞免费鲁丝片| 少妇 在线观看| 国产精品久久久久久久电影| 卡戴珊不雅视频在线播放| 欧美变态另类bdsm刘玥| 国产亚洲欧美精品永久| 国产精品麻豆人妻色哟哟久久| 欧美日韩精品成人综合77777| 欧美精品亚洲一区二区| 我的老师免费观看完整版| 热re99久久国产66热| 日韩精品免费视频一区二区三区 | 精品久久国产蜜桃| 久久久久久久久久久丰满| 蜜臀久久99精品久久宅男| 免费观看的影片在线观看| 国产在线视频一区二区| 18禁在线播放成人免费| 午夜视频国产福利| av又黄又爽大尺度在线免费看| 精品99又大又爽又粗少妇毛片| 亚洲精品视频女| 国产在视频线精品| 国产高清不卡午夜福利| 国产成人免费无遮挡视频| 97超碰精品成人国产| 日韩中文字幕视频在线看片| 亚洲国产精品国产精品| 看非洲黑人一级黄片| 99热全是精品| 欧美精品一区二区大全| 午夜福利影视在线免费观看| 国产男女超爽视频在线观看| 黄色毛片三级朝国网站| 伊人久久精品亚洲午夜| 激情五月婷婷亚洲| 男女无遮挡免费网站观看| 色婷婷av一区二区三区视频| 狠狠精品人妻久久久久久综合| 超色免费av| 曰老女人黄片| 国产女主播在线喷水免费视频网站| 丰满少妇做爰视频| 色94色欧美一区二区| 综合色丁香网| 亚洲性久久影院| 国国产精品蜜臀av免费| 午夜免费观看性视频| 亚洲精品色激情综合| 99久久人妻综合| 成人黄色视频免费在线看| 国产黄色视频一区二区在线观看| 国产熟女欧美一区二区| 成人无遮挡网站| 国产精品一区二区三区四区免费观看| 一级二级三级毛片免费看| 国产淫语在线视频| 国产亚洲欧美精品永久| 久久狼人影院| 日韩av免费高清视频| 久久久a久久爽久久v久久| 精品国产一区二区久久| 久久精品国产亚洲av天美| 五月伊人婷婷丁香| 婷婷成人精品国产| 日本色播在线视频| 欧美成人精品欧美一级黄| 18在线观看网站| 在线 av 中文字幕| 午夜福利网站1000一区二区三区| a级毛片黄视频| 国产精品久久久久久久电影| 婷婷色综合www| 波野结衣二区三区在线| av卡一久久| 不卡视频在线观看欧美| av国产精品久久久久影院| 亚洲,欧美,日韩| 91久久精品国产一区二区三区| 一级毛片aaaaaa免费看小| 免费不卡的大黄色大毛片视频在线观看| 熟妇人妻不卡中文字幕| 亚洲怡红院男人天堂| 欧美 日韩 精品 国产| 欧美日韩成人在线一区二区| 韩国av在线不卡| 男女啪啪激烈高潮av片| 97精品久久久久久久久久精品| 人成视频在线观看免费观看| 国产 精品1| 国产成人精品久久久久久| av电影中文网址| 色网站视频免费| 18禁在线无遮挡免费观看视频| 九草在线视频观看| 午夜福利视频精品| 久久国产精品大桥未久av| 美女福利国产在线| 国产高清不卡午夜福利| 亚洲精品国产av成人精品| 欧美日韩综合久久久久久| 大片免费播放器 马上看| 黄色配什么色好看| 美女福利国产在线| 成人亚洲欧美一区二区av| 99久久综合免费| 精品酒店卫生间| 精品少妇久久久久久888优播| 欧美激情国产日韩精品一区| 观看美女的网站| 国产毛片在线视频| 久久99热这里只频精品6学生| 中国美白少妇内射xxxbb| 狂野欧美激情性xxxx在线观看| 少妇被粗大的猛进出69影院 | 国产成人免费无遮挡视频| 国产欧美亚洲国产| 久久国产亚洲av麻豆专区| 日韩,欧美,国产一区二区三区| xxxhd国产人妻xxx| 亚洲人成网站在线观看播放| 国产伦精品一区二区三区视频9| 国产片特级美女逼逼视频| 亚洲人成网站在线观看播放| 街头女战士在线观看网站| 夜夜骑夜夜射夜夜干| 黄色配什么色好看| 在线观看三级黄色| 亚洲精品国产av蜜桃| 这个男人来自地球电影免费观看 | 男人爽女人下面视频在线观看| 国产精品久久久久久久电影| 国产老妇伦熟女老妇高清| 欧美精品一区二区大全| 韩国av在线不卡| 免费人妻精品一区二区三区视频| 久热久热在线精品观看| 精品国产一区二区久久| 中文字幕人妻熟人妻熟丝袜美| 草草在线视频免费看| 国产精品一二三区在线看| 精品人妻熟女毛片av久久网站| 大码成人一级视频| 国产精品久久久久久久久免| 3wmmmm亚洲av在线观看| 伊人久久国产一区二区| 爱豆传媒免费全集在线观看| av专区在线播放| 男女国产视频网站| 在线精品无人区一区二区三| 91久久精品国产一区二区成人| 人人妻人人澡人人爽人人夜夜| 欧美日韩国产mv在线观看视频| 九九在线视频观看精品| 亚洲国产色片| 午夜日本视频在线| 蜜桃国产av成人99| 九九久久精品国产亚洲av麻豆| 高清在线视频一区二区三区| 欧美精品国产亚洲| 国产精品人妻久久久久久| 黑人欧美特级aaaaaa片| 国产精品.久久久| 丰满少妇做爰视频| 亚洲精品aⅴ在线观看| 最近的中文字幕免费完整| av在线观看视频网站免费| 婷婷成人精品国产| 国产一区二区在线观看av| 91成人精品电影| 国产男人的电影天堂91| kizo精华| 国产av国产精品国产| 女性生殖器流出的白浆| 欧美 亚洲 国产 日韩一| 男的添女的下面高潮视频| 久久人人爽人人片av| 国产精品成人在线| 中文字幕免费在线视频6| 草草在线视频免费看| 一个人免费看片子| 男女高潮啪啪啪动态图| 久久久久久久久久成人| 91精品国产国语对白视频| 精品国产一区二区三区久久久樱花| 国内精品宾馆在线| 18在线观看网站| 91精品国产九色| 亚洲国产av影院在线观看| 国产精品麻豆人妻色哟哟久久| 日本猛色少妇xxxxx猛交久久| 亚洲欧美成人精品一区二区| av在线播放精品| 黑丝袜美女国产一区| 国产精品熟女久久久久浪| 男女啪啪激烈高潮av片| 嘟嘟电影网在线观看| 国产一区二区在线观看av| 国产成人精品一,二区| 免费av中文字幕在线| 亚洲国产精品999| 久久精品久久精品一区二区三区| 人体艺术视频欧美日本| 亚洲欧美一区二区三区国产| 国产精品成人在线| 久久久久久久久久久丰满| 亚洲国产欧美在线一区| 欧美激情 高清一区二区三区| 日日爽夜夜爽网站| 天堂中文最新版在线下载| 亚洲人与动物交配视频| 精品久久久久久久久亚洲| videos熟女内射| 精品久久国产蜜桃| 麻豆成人av视频| 久久国产亚洲av麻豆专区| 欧美日韩一区二区视频在线观看视频在线| 久久久久国产精品人妻一区二区| 日韩伦理黄色片| 亚洲国产精品国产精品| 高清黄色对白视频在线免费看| 午夜激情av网站| 国产 一区精品| 老熟女久久久| 国产亚洲午夜精品一区二区久久| 亚洲久久久国产精品| 日本爱情动作片www.在线观看| 国产高清不卡午夜福利| 熟女人妻精品中文字幕| 亚洲精品日韩av片在线观看| 一个人免费看片子| 青春草亚洲视频在线观看| 国产白丝娇喘喷水9色精品| 高清欧美精品videossex| 这个男人来自地球电影免费观看 | 亚洲综合色网址| 少妇被粗大的猛进出69影院 | 日本欧美国产在线视频| 亚洲精品,欧美精品| 老熟女久久久| 免费av中文字幕在线| 久久鲁丝午夜福利片| 一区二区三区乱码不卡18| 久久久久久久久久久久大奶| 在线亚洲精品国产二区图片欧美 | 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 2022亚洲国产成人精品| 亚洲不卡免费看| 国产男女内射视频| 99re6热这里在线精品视频| 一级毛片我不卡| 视频在线观看一区二区三区| 18禁观看日本| 精品久久国产蜜桃| av有码第一页| 精品一区二区三区视频在线| 欧美人与性动交α欧美精品济南到 | 一边亲一边摸免费视频| 欧美日本中文国产一区发布| 日韩亚洲欧美综合| 国产午夜精品久久久久久一区二区三区| 国产乱人偷精品视频| 多毛熟女@视频| 亚洲中文av在线| 亚洲国产色片| 欧美bdsm另类| 大又大粗又爽又黄少妇毛片口| 国产精品女同一区二区软件| 国产精品久久久久久久电影| 久久久精品区二区三区| 中国三级夫妇交换| 国产毛片在线视频| 免费观看的影片在线观看| 大香蕉久久成人网| 国产日韩欧美亚洲二区| 成人漫画全彩无遮挡| 人体艺术视频欧美日本| 青春草国产在线视频| 色婷婷av一区二区三区视频| 97超碰精品成人国产| 国模一区二区三区四区视频| 国产精品国产三级专区第一集| 国产乱人偷精品视频| 国产高清国产精品国产三级| 九九在线视频观看精品| 丰满饥渴人妻一区二区三| 桃花免费在线播放| 人妻一区二区av| 中文字幕精品免费在线观看视频 | 欧美日韩视频高清一区二区三区二| 精品国产露脸久久av麻豆| a级毛色黄片| 亚洲综合色惰| 国产精品人妻久久久久久| 丰满少妇做爰视频| 国产黄频视频在线观看| 国产一区有黄有色的免费视频| 观看美女的网站| 大话2 男鬼变身卡| 在线亚洲精品国产二区图片欧美 | 看免费成人av毛片| 国产精品一区www在线观看| 亚洲欧洲国产日韩| 国产精品人妻久久久久久| 丰满少妇做爰视频| 街头女战士在线观看网站| 一级爰片在线观看| 亚洲,一卡二卡三卡| 高清黄色对白视频在线免费看| 欧美日韩亚洲高清精品| 男女边吃奶边做爰视频| 亚洲人成77777在线视频| 国产日韩欧美视频二区| 我的女老师完整版在线观看| 色5月婷婷丁香| 国产精品熟女久久久久浪| 看非洲黑人一级黄片| 80岁老熟妇乱子伦牲交| 王馨瑶露胸无遮挡在线观看| 在线观看免费日韩欧美大片 | 一级a做视频免费观看| 91精品一卡2卡3卡4卡| 免费观看无遮挡的男女| 在线观看一区二区三区激情| 汤姆久久久久久久影院中文字幕| 夫妻午夜视频| 不卡视频在线观看欧美| 99久久精品国产国产毛片| 日产精品乱码卡一卡2卡三| 久久久久久人妻| 亚洲欧美一区二区三区黑人 | 日韩av不卡免费在线播放| 久热久热在线精品观看| 一级毛片 在线播放| 国产又色又爽无遮挡免| 日日摸夜夜添夜夜爱| 日本-黄色视频高清免费观看| 亚洲精品第二区| 国产精品无大码| 极品少妇高潮喷水抽搐| a级毛片免费高清观看在线播放| 国产在线免费精品| 国内精品宾馆在线| 亚洲精品久久久久久婷婷小说| 久久久久久久国产电影| 亚洲av福利一区| 国产毛片在线视频| 成人无遮挡网站| 男女边摸边吃奶| 久久久久久久国产电影| 亚洲综合色网址| 欧美激情国产日韩精品一区| 最后的刺客免费高清国语| 日韩一本色道免费dvd| 久久久久久久大尺度免费视频| 一本—道久久a久久精品蜜桃钙片| 80岁老熟妇乱子伦牲交| 亚洲精品av麻豆狂野| av.在线天堂| 日韩一区二区三区影片| 日韩精品有码人妻一区| 青春草视频在线免费观看| 视频中文字幕在线观看| 国产有黄有色有爽视频| 51国产日韩欧美| 精品国产一区二区久久| h视频一区二区三区|