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

    群遮效應(yīng)對(duì)海上結(jié)構(gòu)物波漂移力的低減作用

    2017-12-06 06:45:18何廣華張志剛張子豪張世軍
    關(guān)鍵詞:浮體入射波波數(shù)

    何廣華,張志剛,張子豪,張世軍

    (哈爾濱工業(yè)大學(xué)(威海) 船舶與海洋工程學(xué)院,山東 威海 264209)

    群遮效應(yīng)對(duì)海上結(jié)構(gòu)物波漂移力的低減作用

    何廣華,張志剛,張子豪,張世軍

    (哈爾濱工業(yè)大學(xué)(威海) 船舶與海洋工程學(xué)院,山東 威海 264209)

    針對(duì)波浪作用下海上浮式建筑物的波漂移力問(wèn)題,考慮在海上建筑物周圍合理地布設(shè)小型浮體群以降低海上建筑物的波漂移力。本文基于波浪交互理論和高階邊界元方法,建立了分析三維多浮體間群遮效應(yīng)的數(shù)值模型。計(jì)算結(jié)果表明:通過(guò)優(yōu)化外圍浮體的尺寸可以降低海上建筑物的波漂移力;本數(shù)值模型計(jì)入了浮體間的波浪相互干涉作用,對(duì)群遮效應(yīng)的預(yù)報(bào)更加準(zhǔn)確和接近工程實(shí)際問(wèn)題。隨著入射波頻率的不同,發(fā)生群遮效應(yīng)所對(duì)應(yīng)的外圍浮體尺寸不同;群遮效應(yīng)不僅可以降低所需防護(hù)海上建筑物的波漂移力,也可以降低整個(gè)浮體群的波漂移力。

    群遮效應(yīng); 多浮體; 波浪交互理論; 高階邊界元方法; 波浪干涉; 波漂移力

    在海洋資源和空間利用過(guò)程中,波浪荷載是海上浮式建筑物所受的主要荷載之一[1]。由于波漂移力等波浪載荷的作用,海上浮式建筑物經(jīng)常需要進(jìn)行錨泊定位,波漂移力等波浪荷載直接影響著錨鏈等系纜的安全性和經(jīng)濟(jì)性[2]。目前人們主要關(guān)注:通過(guò)改變海上浮式建筑物的水動(dòng)力特性,提高結(jié)構(gòu)強(qiáng)度和安全系數(shù)等方式來(lái)確保建筑物的安全。但是,建筑物仍置于惡劣海洋環(huán)境中這一事實(shí)并沒(méi)有改變。本文基于一種全新的理念,通過(guò)在海上浮式建筑物周圍布設(shè)一種海洋工程裝備——浮體群(防護(hù)衣),利用浮體群之間的群遮效應(yīng),從根本上降低波浪對(duì)海上浮式建筑物的作用。

    cloaking的相關(guān)研究是近年非常新穎的國(guó)際前沿研究。cloaking最早是由Pendry等[3]在電磁波研究領(lǐng)域發(fā)現(xiàn)。cloaking是指沒(méi)有散射波向外輻射傳播的現(xiàn)象。隨后,這種現(xiàn)象被應(yīng)用到其他領(lǐng)域中,Zigoneanu等在實(shí)驗(yàn)的基礎(chǔ)上分別研究了聲波和地震波的cloaking現(xiàn)象[4-5]。水波領(lǐng)域cloaking現(xiàn)象的研究正處于初步發(fā)展階段,Alam通過(guò)在浮體外圍的水底安裝Ripples,將表面重力波激發(fā)為內(nèi)波,從而實(shí)現(xiàn)了對(duì)水面浮體的保護(hù)[6]。Porter理論上驗(yàn)證了改變海底地形,實(shí)現(xiàn)單圓柱浮體cloaking現(xiàn)象的可能性[7]。Newman通過(guò)在單圓柱浮體周圍布設(shè)一定數(shù)量的圓柱浮體研究了群遮效應(yīng)下的cloaking現(xiàn)象[8]。計(jì)算結(jié)果表明:通過(guò)優(yōu)化外圍圓柱的尺寸參數(shù),中心圓柱的散射波能量可以減小到幾乎為0。

    波浪交互理論最先由Ohkusu[9]應(yīng)用于水波計(jì)算問(wèn)題中,隨后Kagemoto等[10-11]進(jìn)一步拓展了波浪交互理論在多浮體水動(dòng)力計(jì)算中的應(yīng)用。

    本文中浮體群的布置方式與Newman[8]的相同。為了提高計(jì)算精度,基于Kashiwagi[12]的高階邊界元方法進(jìn)行開(kāi)發(fā)和模擬。浮體波漂移力采用Maruo[13]的遠(yuǎn)場(chǎng)方法計(jì)算得到。基于波浪交互理論和高階邊界元法,本文計(jì)算了在不同外圍浮體數(shù)時(shí)的群遮效應(yīng)下,中心浮體、外圍浮體以及整個(gè)浮體群的波漂移力。計(jì)算結(jié)果表明,群遮效應(yīng)不僅可以降低中心浮體的波漂移力,而且可以降低整個(gè)浮體群的波漂移力。

    1 數(shù)學(xué)模型及數(shù)值方法

    1.1速度勢(shì)

    數(shù)值模型以拉普拉斯方程為基本控制方程。假設(shè)流體不可壓縮并且無(wú)粘、無(wú)旋,這樣流場(chǎng)可以用速度勢(shì)表達(dá)。速度勢(shì)Φ(x,y,z;t)滿足線性邊界條件。

    根據(jù)線性疊加定理,速度勢(shì)可以表達(dá)為

    其中

    φ(P)=φI(P)+ψ(P)

    式中:ζI和ω表示入射波波幅和圓頻率,P(x,y,z)表示流場(chǎng)中的點(diǎn),φ(P)為入射勢(shì)φI(P)與擾動(dòng)勢(shì)ψ(P)之和。

    入射勢(shì)可以在柱坐標(biāo)系中表示為

    (1)

    其中

    式中:β表示入射波與x軸的夾角,h表示有限水深。

    Kagemoto等[10]詳細(xì)地介紹了波浪交互理論。這一理論,可充分考慮浮體群之間的波浪相互干涉作用。當(dāng)外圍浮體數(shù)為N時(shí),ψ(P)可以表示為

    式中 Hm(2)表示m階的第二類漢克爾函數(shù)。

    由Maruo[13]的遠(yuǎn)場(chǎng)方法可知,浮體群外一定距離處速度勢(shì)的遠(yuǎn)場(chǎng)表達(dá)式為

    為進(jìn)一步得到第n個(gè)浮體的速度勢(shì),在第n個(gè)浮體上引入局部柱坐標(biāo)系。此時(shí)第n個(gè)浮體的入射勢(shì)由式(1)給出的全局入射勢(shì)和其他浮體產(chǎn)生的擾動(dòng)勢(shì)兩部分組成,因此總的速度勢(shì)可表示為

    Z0(z)e-imθn

    (2)

    其中

    1.2波漂移力

    Kashiwagi等采用朗斯基矩陣推導(dǎo)出了波漂移力的表達(dá)式[14]。本文采用同樣的方法,利用式(2)可以得到第n個(gè)浮體的波漂移力:

    (3)

    2 波漂移力數(shù)值模擬結(jié)果與分析

    為了準(zhǔn)確計(jì)及浮體間相互作用情況下的二階漂移力,本文采用了兩種坐標(biāo)系,如圖1所示。一種是放置在中心浮體上的全局坐標(biāo)系o-xyz,一種是放置在第n個(gè)浮體上的局部坐標(biāo)系on-xnynzn。除了笛卡爾坐標(biāo)系外,本文還采用了柱坐標(biāo)系;z軸豎直向上,每一個(gè)坐標(biāo)系的原點(diǎn)位于未擾動(dòng)的自由表面(z=0)。

    圖1 全局坐標(biāo)系與局部坐標(biāo)系的關(guān)系Fig.1 Global and local coordinate systems

    2.1浮體群模型

    Newman[8]研究了入射波數(shù)為K=1.0時(shí),多浮體群遮效應(yīng)下的散射波能量。

    計(jì)算結(jié)果表明:通過(guò)優(yōu)化外圍浮體的吃水、直徑以及外圍浮體與中心浮體的距離等參數(shù),可以使浮體群的散射波能量減小到幾乎為0,即cloaking現(xiàn)象產(chǎn)生。對(duì)應(yīng)的外圍浮體尺寸參數(shù)與散射波能量,如表1所示。

    本文采用同Newman[8]相同的布置方式,入射波方向?yàn)閤軸負(fù)方向,如圖2所示。多個(gè)圓柱(N=4,8,16)規(guī)律地環(huán)繞在中心圓柱的周圍,并構(gòu)成一個(gè)同心圓。外圍圓柱尺寸相同且等間距分布,中心圓柱的半徑r=0.5,吃水d=1.0,如圖3所示。本文中波漂移力進(jìn)行了無(wú)量綱化處理:F=2Fx/ρgζI2d,其中Fx為x方向的波漂移力,ρ為流體密度,g為重力加速度,ζI為入射波波幅,d為浮體吃水。表1中同時(shí)給出了本文關(guān)于所需保護(hù)浮體波漂移力的計(jì)算結(jié)果。

    表1 外圍浮體尺寸參數(shù)的優(yōu)化值

    圖2 入射波方向與浮體編號(hào)(俯視圖)Fig.2 The direction of incident waves and notations (top view)

    圖3 浮體形狀尺寸(浮體吃水與半徑)Fig.3 Geometrical dimensions of floating bodies(draft and radius of floating body)

    圖4為表1中衰減系數(shù)與外圍浮體數(shù)N的關(guān)系。從圖4可以看出,隨著外圍浮體數(shù)的增多,散射波能量逐漸較小并趨近于0。本文針對(duì)需保護(hù)浮體的波漂移力計(jì)算結(jié)果(圖中F/F0)與散射波能量(圖中E/E0)的變化趨勢(shì)相同,驗(yàn)證了群遮效應(yīng)對(duì)波漂移力的低減作用。需要指出的是:在外圍浮體數(shù)N>8后,群遮效應(yīng)對(duì)波漂移力的低減效果增加并不明顯,建議工程應(yīng)用時(shí),外圍浮體數(shù)不大于8即可。

    圖4 衰減系數(shù)與外圍浮體數(shù)N的關(guān)系Fig.4 The relationship between attenuation coefficient and the number of outer floating bodies

    2.2波漂移力與波數(shù)的關(guān)系

    本節(jié)采用表1中的外圍浮體尺寸參數(shù),計(jì)算了波漂移力與入射波數(shù)K之間的關(guān)系。圖5給出了不同外圍浮體數(shù)(N=4, 8, 16)的中心浮體、浮體群以及外圍浮體的漂移力(合力)。

    在圖5中可以看出孤立圓柱(N=0)的波漂移力隨著波數(shù)的增加而增加。圖5(a)是中心圓柱的波漂移力,可以發(fā)現(xiàn)在外圍布設(shè)圓柱浮體后,中心圓柱的波漂移力相比于單個(gè)圓柱發(fā)生了明顯的變化,這是因?yàn)槎喔◇w之間產(chǎn)生了復(fù)雜的波浪相互干涉。在波數(shù)K=1.0時(shí),中心圓柱的波漂移力相對(duì)于單個(gè)圓柱浮體的波漂移力明顯較小,顯示了群遮效應(yīng)對(duì)波漂移力的低減作用;與Newman[8]關(guān)于散射波能量的研究結(jié)論一致,并且外圍浮體的數(shù)量越多,群遮效果越明顯。圖5(b)、(c)分別給出了對(duì)應(yīng)的浮體群和外圍浮體的漂移力??梢钥吹?,波漂移力在波數(shù)K=1.0時(shí),幾乎同樣減小到0,這表明群遮效應(yīng)不僅可以降低中心海上建筑物的波漂移力,也可以降低整個(gè)浮體群的波漂移力。

    圖6給出了每個(gè)外圍浮體的波漂移力,由于浮體布設(shè)的對(duì)稱性,只需給出其中一半浮體結(jié)果即可(浮體編號(hào)如圖2所示)。從圖6(a)中可以看出,兩圓柱的波漂移力在K=1.0時(shí),方向相反、大小近似;圖5(b)中1號(hào)圓柱與4號(hào)圓柱的的波漂移力方向相反、大小近似,2號(hào)和3號(hào)圓柱的波漂移力幾乎為0;圖5(c)中1號(hào)圓柱與8號(hào)圓柱、2號(hào)圓柱與7號(hào)圓柱、3號(hào)圓柱與6號(hào)圓柱的波漂移力在K=1.0時(shí)方向相反、大小近似,而4號(hào)圓柱與5號(hào)圓柱的波漂移力幾乎為0。這表明正是由于結(jié)構(gòu)關(guān)于y軸的對(duì)稱性,才使得在群遮效應(yīng)發(fā)生時(shí)整個(gè)浮體群的漂移力合力也減小到幾乎為0。還需要指出的是,迎浪側(cè)圓柱的波漂移力絕對(duì)值變化較大,特別是1號(hào)圓柱的變化最大;并且,由于波浪復(fù)雜相互干涉作用,在某些波數(shù)范圍內(nèi)迎浪側(cè)的圓柱波漂移力變?yōu)樨?fù)值(與入射波方向相反)。

    圖5 漂移力與入射波數(shù)的關(guān)系Fig.5 The relationship between wave drift force and wavenumber

    圖6 外圍單個(gè)圓柱的波漂移力Fig.6 The wave drift force on every outer cylinder

    2.3波漂移力與外圍浮體數(shù)N的關(guān)系

    本文中的外圍圓柱浮體的尺寸參數(shù)采用的是Newman[8]在波數(shù)K=1.0時(shí)以散射波能量為目標(biāo)函數(shù)的優(yōu)化結(jié)果,如表1所示。本節(jié)在計(jì)算波數(shù)K=1.0情況下不同外圍浮體數(shù)(N=4, 8, 16)的波漂移力之外,采用同樣的外圍圓柱尺寸還給出了其他波數(shù)(K=0.1, 0.5, 1.5, 2.0)下的波漂移力,計(jì)算結(jié)果如圖7所示。

    從圖7中可以看出,在波數(shù)K=1.0時(shí),中心圓柱的波漂移力與浮體群的波漂移力合力都隨著外圍浮體數(shù)的增加而減小,而在其他波數(shù)則不存在這樣的規(guī)律。說(shuō)明針對(duì)不同的入射波頻率,特別是大波數(shù)情況,需要分別優(yōu)化外圍浮體的尺寸參數(shù)才能實(shí)現(xiàn)多浮體的群遮效應(yīng)。

    圖7 波漂移力與外圍浮體數(shù)N的關(guān)系Fig.7 The relationship between wave drift force and the number of outer cylinders

    2.4本文方法與Fang等提出的方法的對(duì)比

    對(duì)于多浮體波漂移力的計(jì)算,F(xiàn)ang等[15]提出了一種忽略浮體間相互干涉效應(yīng)的簡(jiǎn)便計(jì)算方法。針對(duì)中心浮體波漂移力,本節(jié)給出了本文方法與Fang等[15]方法的對(duì)比結(jié)果,如圖8所示。

    圖8(a)為孤立圓柱時(shí)的波漂移力,本文計(jì)算結(jié)果與Fang等[15]計(jì)算結(jié)果一致。圖8(b)、(c)中兩種方法的計(jì)算結(jié)果出現(xiàn)了差異,特別是在波數(shù)K=1.0附近。群遮效應(yīng)是由于波浪復(fù)雜

    相互干涉作用即產(chǎn)生的,F(xiàn)ang等[15]的方法因忽略了波浪的相互干涉,而在波數(shù)K=1.0(即群遮效應(yīng)最顯著點(diǎn)),計(jì)算結(jié)果誤差較大。此外,在小波數(shù)(長(zhǎng)波)范圍內(nèi),F(xiàn)ang等[15]的方法與本方法的計(jì)算結(jié)果吻合較好,這是因?yàn)樾〔〝?shù)情況下,波浪的相互干涉作用并不明顯,可以忽略。但是在大波數(shù)情況下,浮體間的波浪干涉作用更加復(fù)雜,此時(shí)不可以忽略波浪干涉對(duì)浮體波漂移力的影響。

    圖8 本文計(jì)算方法與Fang等[15]的計(jì)算方法的對(duì)比(中心圓柱波漂移力)Fig.8 Comparison of wave drift forces between Fang′s results[15] and the present solution

    3 結(jié)論

    1)群遮效應(yīng)不僅可以降低中心海上建筑物的波漂移力,也可以降低整個(gè)浮體群的波漂移力。

    2)在目標(biāo)波數(shù)下(K=1.0),中心浮體的波漂移力隨著外圍浮體數(shù)的增多而單調(diào)遞減并逐漸趨于0。

    3)不同入射波頻率,發(fā)生群遮效應(yīng)所對(duì)應(yīng)的外圍浮體尺寸不同。

    4)本文建立的數(shù)值模型計(jì)及了浮體間的波浪相互干涉作用,對(duì)群遮效應(yīng)的預(yù)報(bào)更加準(zhǔn)確。

    [1] 桂琴琴, 鄒志利, 王大國(guó). 箱型船體波浪和波浪壓力的非線性時(shí)域模型[J]. 工程力學(xué), 2011, 28(2): 239-245.

    GUI Qinqin, ZOU Zhili, WANG Daguo. Time-domian model of nonlinear wave forces on a box-shaped ship [J]. Engineering mechanics, 2011, 28(2): 239-245.

    [2] 繆泉明, 顧民, 楊占明, 等. 極限海況下浮標(biāo)運(yùn)動(dòng)及錨鏈?zhǔn)芰浪鉡J]. 船舶力學(xué), 2003, 7(5): 21-27.

    MIU Quanming, GU Min, YANG Zhanming, et al. Estimation of buoy motions and chain forces in extreme sea state [J]. Journal of ship mechanics, 2003, 7(5): 21-27.

    [3] PENDRY J B, SCHURIG D, SMITH D R. Controlling electromagnetic filed[J]. Science, 2006, 312: 1780-1782.

    [4] ZIGONEANU L, POPA B, CUMMER S A. Three-dimensional broadband omnidirectional acoustic ground cloak[J]. Nature materials, 2014, 13(4): 352-355.

    [5] BRULE S, JAVELAUD E H, ENOCH S, et al. Experiments on seismic metamaterials: molding surface waves[J]. Physical review letters, 2014, 112(13): 133901.

    [6] ALAM M. Broadband cloaking in stratified seas[J]. Physical review letters, 2011, 108: 1-5.

    [7] PORTER R. Cloaking of a cylinder in waves[C]//26thInternational Workshop on Water Waves and Floating Bodies, Greece, 2011, 36: 1-4.

    [8] NEWMAN J N. Cloaking a circular cylinder in deep water[C]//28thInternational Workshop on Water Waves and Floating Bodies, L′isle surla Sorgue, France, 2013: 157-160.

    [9] OHKUSU M. Hydrodynamic forces on multiple cylinders in waves[C]//The Process of International Symposium Dynamics of Marine Vehicles and Structures in Waves, London, 1974: 261-268.

    [10] KAGENOTO H, YUE D K P. Three-dimensional bodies in water waves: an exact algebraic method[J]. Journal of fluid mechanics, 1986, 166: 189-209.

    [11] LINTON C M, EVANS D V. The interaction of waves with arrays of cylinders[J]. Journal of fluid mechanics, 1990, 215: 549-569.

    [12] KASHIWAGI M. 3-D calculation for multiple floating bodies in proximity using wave interaction theory[J]. Journal of offshore and polar engineering, 2008, 18: 1-7.

    [13] MARUO H. The drift of a body floating on waves[J]. Journal of ship research, 1960, 4: 1-10.

    [14] KASHIWAGI M, YOSHIDA S. Wave drift force and moment on a VLFS supported by a great number of floating columns[J]. International journal of offshore and polar engineering, 2001, 11: 176-183.

    [15] FANG M C, CHEN G R. On three-dimension solution of drift forces and moments between two ships in waves[J]. Journal of ship research, 2002, 46(4): 280-28.

    本文引用格式:

    何廣華,張志剛,張子豪,等. 群遮效應(yīng)對(duì)海上結(jié)構(gòu)物波漂移力的低減作用[J]. 哈爾濱工程大學(xué)學(xué)報(bào), 2017, 38(11): 1676-1681.

    HE Guanghua, ZHANG Zhigang, ZHANG Zihao, et al. Reduction in wave drift force on marine structures by cloaking phenomenon[J]. Journal of Harbin Engineering University, 2017, 38(11): 1676-1681.

    Reductioninwavedriftforceonmarinestructures
    bycloakingphenomenon

    HE Guanghua, ZHANG Zhigang, ZHANG Zihao, ZHANG Shijun

    (School of Naval Architecture and Ocean Engineering, Harbin Institute of Technology, Weihai 264209, China)

    This study aims to examine the reduction in the wave drift force acting on a floating structure, such as an offshore platform, by appropriate positioning of small floating bodies around the floating structure. To investigate the cloaking phenomenon on 3D multiple floating bodies, a numerical model using a combination of wave interaction theory and the higher-order boundary element method is introduced. Computed results show that reduction in the wave drift force on the offshore structure is realized using optimized parameters of the surrounding floating bodies. Numerical computations were highly accurate and benefited from considering the interactions between multi-bodies, and the optimized parameters of the surrounding bodies were different in relation to the various wave frequencies. Results show that the cloaking phenomenon can reduce the wave drift force on the central structure and on the entire array of floating bodies.

    cloaking phenomenon; multiple floating bodies; wave interaction theory; higher-order boundary element method; wave interactions; wave drift force

    10.11990/jheu.201607009

    http://www.cnki.net/kcms/detail/23.1390.u.20171016.1546.004.html

    U661.1

    A

    1006-7043(2017)11-1676-06

    2016-07-04.

    網(wǎng)絡(luò)出版日期:2017-10-16.

    國(guó)家自然科學(xué)基金項(xiàng)目(11502059); 大連理工大學(xué)海岸和近海工程國(guó)家重點(diǎn)實(shí)驗(yàn)室開(kāi)放課題基金項(xiàng)目(LP1513); 哈爾濱工業(yè)大學(xué)科研創(chuàng)新項(xiàng)目(HIT.NSRIF.201726); 水沙科學(xué)與水災(zāi)害防治湖南省重點(diǎn)實(shí)驗(yàn)室基金(2015SS02); 衛(wèi)星海洋環(huán)境動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室(國(guó)家海洋局第二海洋研究所)項(xiàng)目(SOED1514).

    何廣華(1980-), 男, 教授, 博士生導(dǎo)師;

    張志剛(1992-), 男, 碩士研究生.

    何廣華, E-mail: ghhe@hitwh.edu.cn.

    猜你喜歡
    浮體入射波波數(shù)
    聲場(chǎng)波數(shù)積分截?cái)嗖〝?shù)自適應(yīng)選取方法
    浮體結(jié)構(gòu)沉浮過(guò)程周圍水流特性研究
    一種基于SOM神經(jīng)網(wǎng)絡(luò)中藥材分類識(shí)別系統(tǒng)
    SHPB入射波相似律與整形技術(shù)的試驗(yàn)與數(shù)值研究
    物探船硬浮體陣列自擴(kuò)變量分析與應(yīng)用
    超大型浮體結(jié)構(gòu)碰撞損傷研究
    有限流動(dòng)水域浮體受力及側(cè)傾研究
    瞬態(tài)激勵(lì)狀態(tài)下樁身速度以及樁身內(nèi)力計(jì)算
    對(duì)機(jī)械波半波損失現(xiàn)象的物理解釋
    電子科技(2015年11期)2015-03-06 01:32:24
    重磁異常解釋的歸一化局部波數(shù)法
    国产精品久久久久久久电影| 亚洲 国产 在线| 成年免费大片在线观看| 亚洲国产精品成人综合色| 久久精品影院6| 婷婷精品国产亚洲av在线| 亚洲欧美日韩高清在线视频| 毛片女人毛片| 搡老妇女老女人老熟妇| avwww免费| 国产精品人妻久久久久久| 国产淫片久久久久久久久 | 国产色婷婷99| 久99久视频精品免费| 亚洲欧美日韩卡通动漫| 精品午夜福利在线看| 久久精品人妻少妇| 99久久99久久久精品蜜桃| 欧美日韩亚洲国产一区二区在线观看| 性色avwww在线观看| 九九热线精品视视频播放| 国产精品av视频在线免费观看| ponron亚洲| 国产精品人妻久久久久久| 免费观看的影片在线观看| 亚洲精品色激情综合| 又紧又爽又黄一区二区| 亚洲中文日韩欧美视频| 亚洲久久久久久中文字幕| 午夜精品久久久久久毛片777| 丝袜美腿在线中文| 我要看日韩黄色一级片| 亚洲熟妇熟女久久| 香蕉av资源在线| 亚洲精品日韩av片在线观看| av天堂在线播放| 在线天堂最新版资源| 亚洲一区二区三区不卡视频| 亚洲天堂国产精品一区在线| 久99久视频精品免费| 亚洲黑人精品在线| 女生性感内裤真人,穿戴方法视频| 91麻豆精品激情在线观看国产| 国产三级中文精品| 中文亚洲av片在线观看爽| 成人永久免费在线观看视频| 可以在线观看毛片的网站| 蜜桃亚洲精品一区二区三区| 好男人在线观看高清免费视频| 舔av片在线| 久久久久久国产a免费观看| 日本免费一区二区三区高清不卡| 欧美潮喷喷水| 国产精品人妻久久久久久| 日日摸夜夜添夜夜添av毛片 | a级一级毛片免费在线观看| 小蜜桃在线观看免费完整版高清| 国产精品乱码一区二三区的特点| 亚洲色图av天堂| 一区二区三区激情视频| 国内少妇人妻偷人精品xxx网站| 亚洲欧美日韩无卡精品| 亚洲美女视频黄频| 亚洲精品粉嫩美女一区| 免费搜索国产男女视频| 午夜激情福利司机影院| 真人做人爱边吃奶动态| 好男人电影高清在线观看| 欧美日韩中文字幕国产精品一区二区三区| 欧美色欧美亚洲另类二区| 亚洲国产日韩欧美精品在线观看| 日韩成人在线观看一区二区三区| 我要搜黄色片| 国产单亲对白刺激| 好看av亚洲va欧美ⅴa在| 午夜免费成人在线视频| 少妇裸体淫交视频免费看高清| 97超级碰碰碰精品色视频在线观看| 麻豆一二三区av精品| av在线蜜桃| 精品乱码久久久久久99久播| 性色avwww在线观看| 中文字幕人妻熟人妻熟丝袜美| 亚洲色图av天堂| av欧美777| 成人午夜高清在线视频| 国产精品伦人一区二区| av福利片在线观看| 色吧在线观看| 好男人在线观看高清免费视频| 精品国内亚洲2022精品成人| 神马国产精品三级电影在线观看| 极品教师在线免费播放| 校园春色视频在线观看| 午夜两性在线视频| 少妇丰满av| 国产私拍福利视频在线观看| 亚洲国产精品合色在线| 99久久久亚洲精品蜜臀av| 精品人妻偷拍中文字幕| 最新中文字幕久久久久| 欧美性猛交╳xxx乱大交人| 一级毛片久久久久久久久女| a在线观看视频网站| 国产人妻一区二区三区在| 首页视频小说图片口味搜索| 真实男女啪啪啪动态图| 久久精品国产亚洲av涩爱 | 亚洲真实伦在线观看| 日韩欧美一区二区三区在线观看| 最近视频中文字幕2019在线8| a在线观看视频网站| 国产69精品久久久久777片| 精品久久久久久久久久免费视频| 韩国av一区二区三区四区| 亚洲专区国产一区二区| 精品人妻视频免费看| 床上黄色一级片| 一级黄片播放器| 欧美最新免费一区二区三区 | 在线观看一区二区三区| 午夜福利视频1000在线观看| 成熟少妇高潮喷水视频| 国产精品乱码一区二三区的特点| 美女大奶头视频| 女人十人毛片免费观看3o分钟| 亚洲精品456在线播放app | 国产激情偷乱视频一区二区| 国产免费一级a男人的天堂| 久久精品影院6| 国产精品美女特级片免费视频播放器| 国产精品野战在线观看| 久久精品人妻少妇| 乱人视频在线观看| aaaaa片日本免费| 人妻丰满熟妇av一区二区三区| 成人国产综合亚洲| 日韩欧美国产一区二区入口| 国产免费av片在线观看野外av| 听说在线观看完整版免费高清| 给我免费播放毛片高清在线观看| 特大巨黑吊av在线直播| 久久人妻av系列| 搞女人的毛片| 一进一出抽搐gif免费好疼| 他把我摸到了高潮在线观看| 观看美女的网站| 日韩亚洲欧美综合| 中文字幕免费在线视频6| 91在线观看av| 久久热精品热| 久久中文看片网| 91久久精品国产一区二区成人| 亚洲一区高清亚洲精品| 国产精品久久久久久久电影| 制服丝袜大香蕉在线| 久久精品国产亚洲av天美| 成年女人看的毛片在线观看| 天美传媒精品一区二区| 成人美女网站在线观看视频| 色在线成人网| 欧美性猛交╳xxx乱大交人| 亚洲综合色惰| 日本免费一区二区三区高清不卡| 熟女电影av网| 国产主播在线观看一区二区| 欧美极品一区二区三区四区| 18+在线观看网站| 精品一区二区三区人妻视频| 亚洲精品成人久久久久久| 国产精品三级大全| 亚洲精华国产精华精| 色视频www国产| 真实男女啪啪啪动态图| 婷婷亚洲欧美| av在线观看视频网站免费| 久久婷婷人人爽人人干人人爱| 免费高清视频大片| 免费搜索国产男女视频| 欧美黄色淫秽网站| 中文字幕久久专区| 特级一级黄色大片| 国产乱人视频| 床上黄色一级片| 男人和女人高潮做爰伦理| 免费看日本二区| 欧美激情国产日韩精品一区| 午夜福利18| 国产国拍精品亚洲av在线观看| 久久午夜福利片| 男女那种视频在线观看| 国产一区二区在线观看日韩| 一级黄片播放器| 黄色配什么色好看| 久久人人爽人人爽人人片va | 99久久精品一区二区三区| 日韩精品青青久久久久久| 国内毛片毛片毛片毛片毛片| 免费人成在线观看视频色| 18禁黄网站禁片免费观看直播| 一边摸一边抽搐一进一小说| 91狼人影院| 久久99热这里只有精品18| 亚洲成人久久爱视频| 国产精品一区二区三区四区免费观看 | 两人在一起打扑克的视频| 老熟妇乱子伦视频在线观看| 国产精品三级大全| 国产色爽女视频免费观看| 亚洲专区中文字幕在线| 青草久久国产| 一个人看的www免费观看视频| 国产黄a三级三级三级人| 亚洲av不卡在线观看| 国产高潮美女av| 岛国在线免费视频观看| 亚洲在线观看片| 亚洲经典国产精华液单 | 在线十欧美十亚洲十日本专区| 中文亚洲av片在线观看爽| 精品久久久久久,| 99国产综合亚洲精品| 精品熟女少妇八av免费久了| 国产aⅴ精品一区二区三区波| 91狼人影院| 日韩欧美国产一区二区入口| 国产精品精品国产色婷婷| bbb黄色大片| 真人做人爱边吃奶动态| 一区二区三区高清视频在线| 美女高潮喷水抽搐中文字幕| 亚洲av第一区精品v没综合| 性欧美人与动物交配| 国产精品久久久久久久电影| 精品免费久久久久久久清纯| 九色国产91popny在线| 午夜影院日韩av| 成人无遮挡网站| 国产欧美日韩精品一区二区| 午夜福利欧美成人| 亚洲av免费高清在线观看| 9191精品国产免费久久| 国产av不卡久久| 亚洲美女搞黄在线观看 | av福利片在线观看| av视频在线观看入口| 天堂av国产一区二区熟女人妻| 久久99热这里只有精品18| 久久久久性生活片| 熟女电影av网| 日本a在线网址| 免费观看精品视频网站| 亚洲av美国av| 欧美国产日韩亚洲一区| 一级a爱片免费观看的视频| 在线观看午夜福利视频| 国产精品一区二区三区四区免费观看 | 国产免费男女视频| 欧美性猛交黑人性爽| 级片在线观看| 久久精品夜夜夜夜夜久久蜜豆| 国产成人aa在线观看| 亚洲专区国产一区二区| 欧美成人性av电影在线观看| 麻豆一二三区av精品| 国产精品,欧美在线| 久久精品国产99精品国产亚洲性色| 国产极品精品免费视频能看的| 99国产精品一区二区蜜桃av| 国产精品99久久久久久久久| 99在线人妻在线中文字幕| 国产精品永久免费网站| 亚洲成人精品中文字幕电影| 麻豆国产av国片精品| h日本视频在线播放| 嫩草影院精品99| 久久人妻av系列| 久久精品久久久久久噜噜老黄 | 9191精品国产免费久久| a级毛片免费高清观看在线播放| 亚洲天堂国产精品一区在线| 日韩欧美一区二区三区在线观看| 亚洲美女搞黄在线观看 | 嫩草影视91久久| 久久性视频一级片| 国产精品亚洲一级av第二区| 中亚洲国语对白在线视频| 国产熟女xx| 变态另类丝袜制服| 亚洲在线观看片| 久久这里只有精品中国| 亚洲第一欧美日韩一区二区三区| 久久国产乱子伦精品免费另类| 成人av在线播放网站| 国语自产精品视频在线第100页| 偷拍熟女少妇极品色| 99热精品在线国产| 久久精品影院6| 蜜桃久久精品国产亚洲av| 亚洲成人精品中文字幕电影| 国产美女午夜福利| 国产精品亚洲一级av第二区| 国产黄a三级三级三级人| 在线免费观看不下载黄p国产 | 亚洲中文字幕一区二区三区有码在线看| 乱人视频在线观看| 午夜免费男女啪啪视频观看 | 男女下面进入的视频免费午夜| 色哟哟·www| 欧美xxxx性猛交bbbb| 人人妻人人看人人澡| 国产免费av片在线观看野外av| 别揉我奶头~嗯~啊~动态视频| 特级一级黄色大片| 国产成年人精品一区二区| 国产美女午夜福利| 人妻制服诱惑在线中文字幕| 免费av毛片视频| 99热6这里只有精品| 午夜视频国产福利| 村上凉子中文字幕在线| 在线观看美女被高潮喷水网站 | 搡老熟女国产l中国老女人| 色在线成人网| 在线观看午夜福利视频| 精品国内亚洲2022精品成人| 午夜福利在线在线| 最近中文字幕高清免费大全6 | 亚洲成人久久爱视频| 色5月婷婷丁香| 国产伦精品一区二区三区视频9| 国产亚洲精品久久久com| 中文资源天堂在线| avwww免费| 国产成人aa在线观看| 不卡一级毛片| 男人和女人高潮做爰伦理| 亚洲人成网站在线播放欧美日韩| 男人的好看免费观看在线视频| 国语自产精品视频在线第100页| 一区二区三区四区激情视频 | 欧美乱妇无乱码| 亚洲中文字幕日韩| 成年人黄色毛片网站| 全区人妻精品视频| 久久午夜亚洲精品久久| 九九在线视频观看精品| 青草久久国产| 999久久久精品免费观看国产| 美女cb高潮喷水在线观看| 国产亚洲精品久久久久久毛片| 99久久无色码亚洲精品果冻| 免费av不卡在线播放| 女人十人毛片免费观看3o分钟| 亚洲va日本ⅴa欧美va伊人久久| 在线观看美女被高潮喷水网站 | 免费看a级黄色片| 亚洲精品一区av在线观看| 五月伊人婷婷丁香| 亚洲欧美日韩卡通动漫| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 精品一区二区三区人妻视频| 嫩草影院精品99| 首页视频小说图片口味搜索| 波野结衣二区三区在线| 99久久精品热视频| 免费在线观看成人毛片| 久久国产精品影院| 国产精品爽爽va在线观看网站| 亚洲第一欧美日韩一区二区三区| 亚洲成人久久性| 亚洲av熟女| 亚洲av中文字字幕乱码综合| 一a级毛片在线观看| 亚洲成人久久爱视频| 亚洲成人免费电影在线观看| 亚洲五月婷婷丁香| 国产一区二区三区视频了| 真人一进一出gif抽搐免费| 久久精品91蜜桃| 国产亚洲av嫩草精品影院| 一区二区三区激情视频| 简卡轻食公司| 在线观看av片永久免费下载| 国产精品美女特级片免费视频播放器| 婷婷精品国产亚洲av在线| 97超级碰碰碰精品色视频在线观看| 最好的美女福利视频网| 国产单亲对白刺激| 精品国产三级普通话版| 国产又黄又爽又无遮挡在线| 变态另类丝袜制服| 亚洲国产日韩欧美精品在线观看| 亚洲七黄色美女视频| 99久久九九国产精品国产免费| 免费人成视频x8x8入口观看| 又紧又爽又黄一区二区| 好男人在线观看高清免费视频| 成人av一区二区三区在线看| 在线免费观看不下载黄p国产 | 18禁在线播放成人免费| 又黄又爽又免费观看的视频| www.www免费av| 午夜福利在线观看免费完整高清在 | 激情在线观看视频在线高清| 一边摸一边抽搐一进一小说| 悠悠久久av| 狂野欧美白嫩少妇大欣赏| av中文乱码字幕在线| 性欧美人与动物交配| 3wmmmm亚洲av在线观看| 中文字幕人妻熟人妻熟丝袜美| 88av欧美| 午夜免费男女啪啪视频观看 | a级毛片免费高清观看在线播放| 9191精品国产免费久久| 女人十人毛片免费观看3o分钟| h日本视频在线播放| 高清毛片免费观看视频网站| 久久久久久九九精品二区国产| 国产美女午夜福利| 深夜a级毛片| 精品久久久久久成人av| 美女高潮喷水抽搐中文字幕| 亚洲中文日韩欧美视频| а√天堂www在线а√下载| 国产精品一区二区三区四区免费观看 | 国产精品av视频在线免费观看| 中文在线观看免费www的网站| 欧美黄色片欧美黄色片| 色综合婷婷激情| 日韩欧美免费精品| 亚洲专区国产一区二区| 一级作爱视频免费观看| 欧美日本亚洲视频在线播放| 欧美另类亚洲清纯唯美| 国产成年人精品一区二区| 99久久精品一区二区三区| 欧美性猛交╳xxx乱大交人| 757午夜福利合集在线观看| 午夜福利在线观看吧| 亚洲 欧美 日韩 在线 免费| 亚洲av二区三区四区| 欧美成人a在线观看| 欧美三级亚洲精品| 国内精品久久久久久久电影| 欧美bdsm另类| 最近最新中文字幕大全电影3| 少妇的逼好多水| 赤兔流量卡办理| 少妇裸体淫交视频免费看高清| 免费无遮挡裸体视频| 禁无遮挡网站| 波多野结衣高清无吗| 欧美中文日本在线观看视频| 欧美午夜高清在线| 99久久精品一区二区三区| 超碰av人人做人人爽久久| 一个人看的www免费观看视频| 久久九九热精品免费| 久久精品人妻少妇| 亚洲av电影不卡..在线观看| 99热只有精品国产| 91在线观看av| 久久天躁狠狠躁夜夜2o2o| 精品一区二区三区视频在线观看免费| 免费在线观看亚洲国产| 亚洲美女视频黄频| 成人欧美大片| 日日干狠狠操夜夜爽| 国产精品国产高清国产av| 嫁个100分男人电影在线观看| 在线观看av片永久免费下载| 亚洲av不卡在线观看| 亚洲,欧美精品.| 亚洲中文日韩欧美视频| 免费人成视频x8x8入口观看| 人妻夜夜爽99麻豆av| 成年免费大片在线观看| 一进一出好大好爽视频| 少妇的逼水好多| 天美传媒精品一区二区| 成人特级黄色片久久久久久久| 成人av在线播放网站| 少妇人妻一区二区三区视频| 国产av一区在线观看免费| 精品久久久久久久久av| netflix在线观看网站| 亚洲国产精品成人综合色| 国产一区二区三区视频了| 亚洲成a人片在线一区二区| 国产伦精品一区二区三区视频9| 国内精品久久久久精免费| 国产精品女同一区二区软件 | 欧美午夜高清在线| 午夜日韩欧美国产| 老鸭窝网址在线观看| 12—13女人毛片做爰片一| 欧美中文日本在线观看视频| 欧美黄色片欧美黄色片| 一本综合久久免费| 美女cb高潮喷水在线观看| 99国产极品粉嫩在线观看| 精品99又大又爽又粗少妇毛片 | 亚洲成人久久爱视频| 亚洲精品成人久久久久久| 国产精品一区二区免费欧美| 亚洲性夜色夜夜综合| 变态另类成人亚洲欧美熟女| 可以在线观看的亚洲视频| 一a级毛片在线观看| 怎么达到女性高潮| 91在线精品国自产拍蜜月| 午夜福利在线观看免费完整高清在 | 又紧又爽又黄一区二区| 很黄的视频免费| 久久这里只有精品中国| 麻豆一二三区av精品| 国产精品久久久久久人妻精品电影| 每晚都被弄得嗷嗷叫到高潮| 亚洲第一欧美日韩一区二区三区| www.色视频.com| 中文字幕高清在线视频| 一区二区三区免费毛片| 神马国产精品三级电影在线观看| 久久久久精品国产欧美久久久| 丰满乱子伦码专区| .国产精品久久| 精品人妻熟女av久视频| 精品国产三级普通话版| 赤兔流量卡办理| 亚洲精品影视一区二区三区av| 久久热精品热| 日本在线视频免费播放| 国产精品爽爽va在线观看网站| av福利片在线观看| 亚洲欧美激情综合另类| 久久中文看片网| 91狼人影院| 久久性视频一级片| 97人妻精品一区二区三区麻豆| 国产av一区在线观看免费| 88av欧美| 赤兔流量卡办理| 国产精品女同一区二区软件 | 内射极品少妇av片p| 免费看美女性在线毛片视频| 亚洲无线观看免费| 观看美女的网站| 婷婷亚洲欧美| 亚洲五月天丁香| 夜夜看夜夜爽夜夜摸| 国产真实乱freesex| 我的老师免费观看完整版| 精品人妻熟女av久视频| 婷婷精品国产亚洲av在线| 欧美3d第一页| avwww免费| 国产 一区 欧美 日韩| 哪里可以看免费的av片| 色综合婷婷激情| 免费大片18禁| 成年版毛片免费区| 亚洲五月婷婷丁香| 国产一区二区三区视频了| 成人性生交大片免费视频hd| 午夜免费成人在线视频| 99久国产av精品| 他把我摸到了高潮在线观看| 高潮久久久久久久久久久不卡| 国产精品三级大全| 国产精品不卡视频一区二区 | 免费av不卡在线播放| 在线看三级毛片| 深夜精品福利| 成人永久免费在线观看视频| 亚洲欧美日韩高清专用| 国产主播在线观看一区二区| netflix在线观看网站| 给我免费播放毛片高清在线观看| 少妇被粗大猛烈的视频| 看黄色毛片网站| 男人的好看免费观看在线视频| 美女大奶头视频| 欧美一级a爱片免费观看看| 女人十人毛片免费观看3o分钟| 国产乱人视频| 国产aⅴ精品一区二区三区波| 搡老妇女老女人老熟妇| 欧美一区二区精品小视频在线| 欧美高清性xxxxhd video| 九色成人免费人妻av| 久久国产精品影院| 国产成人欧美在线观看| 麻豆国产97在线/欧美| 变态另类成人亚洲欧美熟女| 国产精品亚洲av一区麻豆| 99在线人妻在线中文字幕| av福利片在线观看| 国产精品亚洲av一区麻豆| 床上黄色一级片| 欧美日韩综合久久久久久 | 国产成人a区在线观看| 久久6这里有精品| 91午夜精品亚洲一区二区三区 | 我的老师免费观看完整版| 久久伊人香网站| 欧美日韩乱码在线| 91久久精品电影网| 亚洲天堂国产精品一区在线| а√天堂www在线а√下载| 又粗又爽又猛毛片免费看| 国产综合懂色| 99国产精品一区二区蜜桃av| 亚洲真实伦在线观看|