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

    水流對(duì)復(fù)合筒型基礎(chǔ)波浪爬升影響機(jī)制的數(shù)值模擬研究

    2022-10-31 07:05:18于通順宋昊陽(yáng)劉梅梅章譽(yù)天
    關(guān)鍵詞:波流波面水槽

    于通順,宋昊陽(yáng),劉梅梅,章譽(yù)天

    水流對(duì)復(fù)合筒型基礎(chǔ)波浪爬升影響機(jī)制的數(shù)值模擬研究

    于通順1,宋昊陽(yáng)1,劉梅梅2,章譽(yù)天3

    (1. 中國(guó)海洋大學(xué)工程學(xué)院,青島 266100;2. 中交四航工程研究院有限公司,廣州 510230;3.重慶市墊江縣水利局,重慶 408300)

    復(fù)合筒型基礎(chǔ)(CBF)是一種新型的海上風(fēng)電基礎(chǔ),當(dāng)波浪海流與其大尺寸反弧段作用時(shí),會(huì)產(chǎn)生復(fù)雜的波浪爬升現(xiàn)象,嚴(yán)重威脅基礎(chǔ)外平臺(tái)的安全,因此,有必要對(duì)波流作用下復(fù)合筒型基礎(chǔ)的波浪爬升現(xiàn)象開(kāi)展研究.首先,采用RNG-模型建立了波流與復(fù)合筒型基礎(chǔ)相互作用的數(shù)值水槽.其次,基于試驗(yàn)數(shù)據(jù)驗(yàn)證了數(shù)值模型的準(zhǔn)確性.最后,探究了波流作用下復(fù)合筒型基礎(chǔ)周?chē)牟ɡ伺郎?guī)律,揭示了不同流速比cw下水流對(duì)波浪爬升的貢獻(xiàn).研究發(fā)現(xiàn):①波浪單獨(dú)作用下復(fù)合筒型基礎(chǔ)前方水體回落到最低點(diǎn)時(shí),水質(zhì)點(diǎn)的運(yùn)動(dòng)方向均為斜向下方,而在波流聯(lián)合作用時(shí),復(fù)合筒型基礎(chǔ)前端形成漩渦,背面的水質(zhì)點(diǎn)將朝著不同的方向進(jìn)行擴(kuò)散;②波浪單獨(dú)作用時(shí),復(fù)合筒型基礎(chǔ)背浪面水體呈“公雞尾”狀,而波流聯(lián)合作用時(shí),復(fù)合筒型基礎(chǔ)背浪面水體呈“雙髻鯊”狀;③當(dāng)波陡參數(shù)為0.025時(shí),流速比由0.46增大到0.59時(shí),波浪爬升高度增大約0.55倍.

    復(fù)合筒型基礎(chǔ);波浪爬升;波流相互作用;數(shù)值模擬

    2020年,我國(guó)提出力爭(zhēng)于2030年前二氧化碳排放達(dá)到峰值,2060年前實(shí)現(xiàn)碳中和.我國(guó)海上風(fēng)能資源豐富,可開(kāi)發(fā)量達(dá)7.5億千瓦[1],在“碳中和”的背景下海上風(fēng)電迎來(lái)歷史性的發(fā)展機(jī)遇.海上風(fēng)電基礎(chǔ)處于復(fù)雜的海洋環(huán)境中,波流作用于風(fēng)電基礎(chǔ)時(shí)會(huì)產(chǎn)生強(qiáng)烈的波浪爬升現(xiàn)象(如圖1所示),引發(fā)靠船設(shè)備及工作平臺(tái)的沖擊破壞,海上風(fēng)電基礎(chǔ)波浪爬升的研究成為近年來(lái)海洋工程領(lǐng)域的研究熱點(diǎn).

    圖1?海上風(fēng)電基礎(chǔ)波浪爬升

    復(fù)合筒型基礎(chǔ)[2-5](如圖2所示)作為一種新型的海上風(fēng)電基礎(chǔ),因其抗傾覆能力和對(duì)軟土地基的適應(yīng)能力強(qiáng)等優(yōu)勢(shì)受到越來(lái)越多的關(guān)注.復(fù)合筒型基礎(chǔ)大尺寸的反弧段結(jié)構(gòu)及寬淺式的筒裙為其包括水動(dòng)力特征、承載特性、傳力規(guī)律等工作性能的明確帶來(lái)了挑戰(zhàn),有學(xué)者對(duì)復(fù)合筒型基礎(chǔ)承載性能[6]、應(yīng)力響應(yīng)及傳力特性[7]、地基沖刷[8-9]等方面開(kāi)展了大量研究,然而不同波流環(huán)境下復(fù)合筒型基礎(chǔ)波浪爬升特征尚缺乏系統(tǒng)性分析,這也在一定程度上影響了復(fù)合筒型基礎(chǔ)適用水深的確定,增加了靠船設(shè)備及平臺(tái)的運(yùn)行風(fēng)險(xiǎn),限制了其大規(guī)模推廣.因此,研究復(fù)雜波流環(huán)境中復(fù)合筒型基礎(chǔ)的波浪爬升特性很有必要.

    圖2?海上風(fēng)電復(fù)合筒型基礎(chǔ)

    關(guān)于波流作用下海上風(fēng)電基礎(chǔ)波浪爬升的研究已經(jīng)展開(kāi),這些研究的對(duì)象主要集中于單樁基礎(chǔ)[10-13]、多樁基礎(chǔ)[14]及重力式基礎(chǔ)[15-16],針對(duì)復(fù)合筒型基礎(chǔ)這種新型結(jié)構(gòu)波浪爬升的研究相對(duì)較少.通過(guò)單樁基礎(chǔ)波浪爬升的研究結(jié)果[10-11],可以對(duì)復(fù)合筒型基礎(chǔ)上波浪爬升特征有初步的了解,但是考慮到復(fù)合筒型基礎(chǔ)大尺寸反弧段結(jié)構(gòu)這種復(fù)雜的邊界條件,復(fù)合筒型基礎(chǔ)周?chē)ɡ伺郎^(guò)程中流場(chǎng)特征會(huì)更加復(fù)雜. Yu等[17]基于物理模型試驗(yàn)的方法探究了波浪單獨(dú)作用下復(fù)合筒型基礎(chǔ)周?chē)呐郎卣?,獲取了不同工況下的最大爬升高度.然而,復(fù)合筒型基礎(chǔ)波浪爬升過(guò)程中的水動(dòng)力機(jī)制尚需深入探究,水流對(duì)波浪爬升影響機(jī)理的研究也需進(jìn)一步開(kāi)展.因此,為深入探討波浪、水流作用過(guò)程中復(fù)合筒型基礎(chǔ)波浪爬升機(jī)制,本文借鑒其他風(fēng)電基礎(chǔ)波浪爬升的數(shù)值模擬研究方?法[13-15],建立了波流與復(fù)合筒型基礎(chǔ)相互作用的數(shù)值模型,研究了海流對(duì)復(fù)合筒型基礎(chǔ)爬升特性的影響,闡述了海流和不同的波浪參數(shù)對(duì)復(fù)合筒型基礎(chǔ)附近波面擾動(dòng)的影響.

    1?控制方程及求解方法

    本文中涉及的流體為不可壓縮黏性牛頓流體,選用連續(xù)性方程和不可壓縮運(yùn)動(dòng)流體的N-S動(dòng)量方程作為流體運(yùn)動(dòng)的控制方程,其具體表達(dá)式如下:

    為了準(zhǔn)確模擬波浪、水流和結(jié)構(gòu)物相互作用時(shí)流體產(chǎn)生的變形破碎,本文采用重組化RNG-模型,其表達(dá)式如下:

    表1?RNG-模型系數(shù)取值

    Tab.1?Coefficient value of RNGk-ε model

    2?三維數(shù)值波流水槽

    2.1?三維數(shù)值波流水槽的建立

    建立計(jì)算區(qū)域大小(長(zhǎng)×寬×高)為10.2m×3m× 0.37m的數(shù)值波浪水槽,數(shù)值波浪水槽最左側(cè)是波流邊界(wave-current),最右側(cè)是出流邊界(outflow),上、下兩側(cè)邊界條件為對(duì)稱(chēng)邊界(symmetry),底部邊界為壁面(wall),頂部邊界為壓力(specified pressure),如圖3所示.為了減少波浪的反射,在最右側(cè)出流邊界前放置了大小為2m×3m×0.37m的消波塊,左側(cè)的阻尼系數(shù)定為1,右側(cè)的阻尼系數(shù)定為4.2.

    2.2?網(wǎng)格劃分及網(wǎng)格獨(dú)立性驗(yàn)證

    數(shù)值模擬的準(zhǔn)確性與網(wǎng)格劃分結(jié)果息息相關(guān),因此,為了更好地分析復(fù)合筒型基礎(chǔ)周?chē)牟ɡ伺郎?yīng),需要對(duì)波高部分的網(wǎng)格進(jìn)行加密,本文選用3種算例,在波高部分加密的網(wǎng)格形式為細(xì)網(wǎng)格、中網(wǎng)格、粗網(wǎng)格,如表2所示.圖3所示水槽內(nèi)先不放置復(fù)合筒型基礎(chǔ),提取波高()0.06m、周期()1.25s、水深()0.25m的波浪作用下基礎(chǔ)位置處的波面時(shí)程曲線(xiàn),如圖4所示.比較發(fā)現(xiàn)中等網(wǎng)格和細(xì)網(wǎng)格的計(jì)算結(jié)果相差不大,粗網(wǎng)格的計(jì)算結(jié)果與其他兩種相差較大.所以,綜合考慮計(jì)算精度和計(jì)算速度,選用中網(wǎng)格對(duì)波高部分進(jìn)行加密.

    圖3?三維數(shù)值波浪水槽

    表2?網(wǎng)格數(shù)量和尺寸

    Tab.2?Meshnumber and size

    圖4?3種網(wǎng)格尺寸的自由面高程歷時(shí)曲線(xiàn)

    2.3?三維數(shù)值波流水槽驗(yàn)證

    提取中等網(wǎng)格狀態(tài)下波高0.06m、周期1.25s、水深0.25m的入射波浪作用下距離水槽入口邊界5m處波面時(shí)程曲線(xiàn)與理論值進(jìn)行對(duì)比如圖5所示,可以看出計(jì)算值與理論值吻合較好,表明此網(wǎng)格狀態(tài)下水槽準(zhǔn)確性較高.

    提取波高0.06m、周期1.4s、水深0.25m時(shí)不同流速下距造波板5m處波高的計(jì)算值與理論值[18]對(duì)比如圖6所示.

    圖5?自由波面計(jì)算值與理論值的時(shí)程曲線(xiàn)

    圖6?波流作用下波高的計(jì)算值與理論值的對(duì)比

    為了驗(yàn)證數(shù)值水槽中基礎(chǔ)周?chē)ɡ伺郎臏?zhǔn)確性,將復(fù)合筒型基礎(chǔ)模型放置于距離造波板5m處,將波高0.06m、周期1.4s、水深0.25m的波浪與3種流速(0.12m/s、0.16m/s、0.20m/s)的水流共同作用下基礎(chǔ)最大波浪爬升高度與已有試驗(yàn)數(shù)據(jù)[17]對(duì)比,如圖7所示.可以看出,計(jì)算值與試驗(yàn)值的波浪爬升高度誤差在5%以?xún)?nèi),說(shuō)明此數(shù)值水槽能夠較好地模擬復(fù)合筒型基礎(chǔ)上的波浪爬升特征.

    圖7 復(fù)合筒型基礎(chǔ)A1處最大波浪爬升高度計(jì)算值和試驗(yàn)值的對(duì)比

    2.4?虛擬觀測(cè)點(diǎn)的位置

    為了深入研究復(fù)合筒型基礎(chǔ)周?chē)牟ɡ伺郎匦裕谄渲車(chē)O(shè)置6組浪高儀,相鄰兩組之間的角度為36°,觀測(cè)點(diǎn)的具體布置如圖8所示.浪高儀至復(fù)合筒型基礎(chǔ)中心的水平距離如表3所示.

    圖8?沿徑向分布的虛擬浪高儀

    表3?虛擬浪高儀的位置

    Tab.3?Positions of virtual wave probes

    注:為復(fù)合筒型基礎(chǔ)在其相應(yīng)水深下對(duì)應(yīng)的水面交界處筒的半徑,即0.25m水深下=8.07cm.

    2.5?工況設(shè)置

    波浪爬升的數(shù)值模擬采用如表4所示工況,分波浪單獨(dú)作用和波流(同向傳播)共同作用兩組.波浪單獨(dú)作用時(shí)選取了兩組波陡(0.025、0.035),波流共同作用時(shí)選取了2種海流速度(0.12m/s、0.20m/s),波陡變化范圍為0.025~0.035.

    表4?波浪爬升工況

    Tab.4?Wave conditions for wave run-up

    注:流速比cw=c/(c+wm),wm為波浪誘導(dǎo)水質(zhì)點(diǎn)流速幅值,c為水流流速;cw=0代表單獨(dú)波浪情況,cw=1代表單獨(dú)水流情況.

    3?結(jié)果與討論

    3.1?波流作用下復(fù)合筒型基礎(chǔ)周?chē)鲌?chǎng)特性

    3.1.1?基礎(chǔ)周?chē)|(zhì)點(diǎn)速度分布

    波浪單獨(dú)作用及波流共同作用下復(fù)合筒型基礎(chǔ)周?chē)Q向剖面流場(chǎng)分布如圖9、圖10所示.從圖9中可以看出:①由于復(fù)合筒型基礎(chǔ)的阻礙作用,水體的動(dòng)能轉(zhuǎn)化為勢(shì)能從而迅速向上爬升;②當(dāng)其前端波浪爬升至最高點(diǎn)時(shí),峰值附近的水體速度接近于零(區(qū)域Q),雖然仍有少量質(zhì)點(diǎn)的向上爬升,但大部分水質(zhì)點(diǎn)的矢量方向?yàn)樾毕蛳?,說(shuō)明水體正在回落,這與單鐵兵等[11]針對(duì)立柱所得到的規(guī)律一致.

    (a)Ucw=0(工況2) ?(b)Ucw=0.54(工況7)

    從圖10可以看出:①波流共同作用下,當(dāng)前端波浪回落到較低位置時(shí),復(fù)合筒型基礎(chǔ)前端流體質(zhì)點(diǎn)緊靠反弧段產(chǎn)生順時(shí)針?lè)较蛐D(zhuǎn)的漩渦(如區(qū)域S1),出現(xiàn)這種情況的原因可能是由于入射波浪遇到復(fù)合筒型基礎(chǔ)反射后的水質(zhì)點(diǎn)、沿復(fù)合筒型基礎(chǔ)兩側(cè)回流的水體,這兩種水體的水質(zhì)點(diǎn)矢量方向不同,發(fā)生碰撞,再加上海流的作用,產(chǎn)生強(qiáng)烈的非線(xiàn)性作用.而僅有波浪作用時(shí),水質(zhì)點(diǎn)的運(yùn)動(dòng)比較規(guī)律,都斜向左下方運(yùn)動(dòng),沒(méi)有出現(xiàn)明顯的漩渦現(xiàn)象;②當(dāng)復(fù)合筒型基礎(chǔ)前端波浪回落到最低點(diǎn)時(shí),復(fù)合筒型基礎(chǔ)背浪面水位上升,在波流共同作用下,復(fù)合筒型基礎(chǔ)背面產(chǎn)生漩渦,圖(b)中的漩渦靠近壁面,水質(zhì)點(diǎn)沿著不同方向迅速擴(kuò)散,如圖(b)中S2區(qū)域所示,而在圖(a)中水質(zhì)點(diǎn)的運(yùn)動(dòng)比較規(guī)律,矢量方向都為斜向右下方.因此可以看出,海流的存在將在一定程度上增加復(fù)合筒型基礎(chǔ)周?chē)娴臄_動(dòng).

    3.1.2?基礎(chǔ)周?chē)鲌?chǎng)分布

    取初始時(shí)刻0=40s、周期=1.4s.波流共同作用及波浪單獨(dú)作用下復(fù)合筒型基礎(chǔ)周?chē)鲌?chǎng)分布如圖11所示.從波流共同作用時(shí)復(fù)合筒型基礎(chǔ)周?chē)牧鲌?chǎng)分布可以看出:在0+/5時(shí),(見(jiàn)圖(a)、(b))波峰到達(dá)復(fù)合筒型基礎(chǔ)前沿,使復(fù)合筒型基礎(chǔ)兩側(cè)出現(xiàn)正向傳播的邊緣波,并且對(duì)稱(chēng)分布.當(dāng)波峰穿過(guò)復(fù)合筒型基礎(chǔ)時(shí),由于基礎(chǔ)的阻礙作用,復(fù)合筒型基礎(chǔ)前方水體壅高,動(dòng)能轉(zhuǎn)化為爬升勢(shì)能,前方的自由表面以對(duì)稱(chēng)形狀向外擾動(dòng)和輻射;0+2/5(見(jiàn)圖(c)、(d))時(shí)波浪沿著復(fù)合筒型基礎(chǔ)壁面繞射,兩股波浪沿著復(fù)合筒型基礎(chǔ)的兩側(cè)向后面?zhèn)鞑?,在?fù)合筒型基礎(chǔ)背面進(jìn)行疊加,波面迅速放大;0+3/5和0+4/5(見(jiàn)圖(e)~(h))時(shí),當(dāng)復(fù)合筒型基礎(chǔ)背浪面波速接近于零時(shí),水體開(kāi)始回落,波浪沿著復(fù)合筒型基礎(chǔ)兩側(cè)進(jìn)行回流,這時(shí),由于基礎(chǔ)前端水體的速度比較大,并且水質(zhì)點(diǎn)的運(yùn)動(dòng)方向與波的傳播方向相反,波浪會(huì)向著復(fù)合筒型基礎(chǔ)迎浪面迅速移動(dòng),形成邊波;0+(見(jiàn)圖(i)、(j))時(shí),接近入射波谷到達(dá)的相位,擾動(dòng)被掃離基礎(chǔ)以產(chǎn)生一對(duì)對(duì)稱(chēng)但非同心的波面,與Sun等[19]在數(shù)值水槽中所觀察到的Type-2散射波場(chǎng)類(lèi)似.反方向流的邊波將與入射波浪或海流相互作用引起波面的非線(xiàn)性擾動(dòng)和大變形,并且波形不斷向遠(yuǎn)方傳播;最后,由于水體能量的不斷耗散,水質(zhì)點(diǎn)的速度也會(huì)不斷減小,波面擾動(dòng)效應(yīng)的影響逐漸變?。ɡ藛为?dú)作用時(shí),復(fù)合筒型基礎(chǔ)周?chē)牧鲌?chǎng)分布,與波流聯(lián)合作用時(shí)的情況進(jìn)行比較得出:①兩者水體回蕩過(guò)程相似,比較(c)、(d)圖可以看出,波浪單獨(dú)作用時(shí),復(fù)合筒型基礎(chǔ)背浪面水體呈“公雞尾”狀,而波流聯(lián)合作用時(shí),復(fù)合筒型基礎(chǔ)背浪面水體呈“雙髻鯊”狀;②復(fù)合筒型基礎(chǔ)背面的水體回落,水質(zhì)點(diǎn)的運(yùn)動(dòng)方向與波的傳播方向相反,波浪會(huì)向著復(fù)合筒型基礎(chǔ)迎浪面迅速移動(dòng),將會(huì)形成邊波,波浪單獨(dú)作用時(shí),水體自由表面分布比較光滑和均勻,水體沒(méi)有出現(xiàn)較大的變形,波流聯(lián)合作用時(shí)的波面速度變化大,擾動(dòng)劇烈,擾動(dòng)范圍約為2,說(shuō)明海流在一定程度上增加波面的非線(xiàn)性特征.

    圖10 波浪回落到最低點(diǎn)時(shí)復(fù)合筒型基礎(chǔ)周?chē)Q向剖面流場(chǎng)分布

    3.2?水流對(duì)復(fù)合筒型基礎(chǔ)波浪爬升高度的影響

    本文中對(duì)波浪爬升值Ru用入射波幅/2進(jìn)行了無(wú)量綱化處理.圖12為波浪沿著復(fù)合筒型基礎(chǔ)反弧段在0°、36°、72°、108°、144°和180°時(shí)的爬升分布規(guī)律.從圖中可以看出:①波浪爬升倍數(shù)Ru/(/2)的最大值都出現(xiàn)在復(fù)合筒型基礎(chǔ)迎浪面正中心即0°附近,并且在同一速度條件下,隨著波陡增大,爬升倍數(shù)也相應(yīng)增大.波浪沿著復(fù)合筒型基礎(chǔ)兩側(cè)傳播過(guò)程中爬升倍數(shù)逐漸降低,但在背浪面108°附近又開(kāi)始增大,出現(xiàn)這種現(xiàn)象的原因是波浪遇到復(fù)合筒型基礎(chǔ)阻擋后,會(huì)有兩股波浪沿著復(fù)合筒型基礎(chǔ)側(cè)面繼續(xù)傳播,在復(fù)合筒型基礎(chǔ)背面進(jìn)行疊加,使得波面較高;②波浪爬升倍數(shù)的最小值約在108°處,隨著水流流速的增大,波浪爬升倍數(shù)的最小值也相應(yīng)增大,這是因?yàn)楹A骱筒ɡ说南嗷プ饔眉哟罅怂w的能量,從而使得轉(zhuǎn)化的勢(shì)能更大導(dǎo)致了爬升倍數(shù)增大.這與Lin等[16]得出單樁支撐結(jié)構(gòu)的最大爬升高度出現(xiàn)在0°的位置,而相對(duì)最小的爬升高度分別位于112.5°附近基本一致;③波流聯(lián)合作用時(shí),144°到180°附近的爬升倍數(shù)變化較為平緩,而波浪單獨(dú)作用時(shí),爬升值變化幅度較大,與圖11(c)、(d)相對(duì)應(yīng);④復(fù)合筒型基礎(chǔ)迎浪側(cè)的相對(duì)爬升倍數(shù)隨波陡增大而增大,而在復(fù)合筒型基礎(chǔ)背浪一側(cè),對(duì)于更高的波陡相對(duì)爬升倍數(shù)更小.這與Lin等[15]所得到的規(guī)律:高波陡情況下爬升厚度比低波陡要薄,較薄波舌的能量傳遞對(duì)前側(cè)爬升高度增加的貢獻(xiàn)很大,但會(huì)導(dǎo)致后側(cè)的上升高度降低基本一致.

    圖12 不同波流組合下復(fù)合筒型基礎(chǔ)反弧段不同周向上波浪爬升的分布規(guī)律

    圖13~16分別給出了A1、C1、E1、F1處波浪爬升時(shí)程和相應(yīng)的頻譜圖.當(dāng)波浪和水流沿著復(fù)合筒型基礎(chǔ)繞射時(shí)將會(huì)出現(xiàn)以下特征:①對(duì)比圖13~16可以看出,隨著流速比的增大,波浪爬升高度也相應(yīng)增加,功率譜也相應(yīng)增大,波浪爬升在復(fù)合筒型基礎(chǔ)迎浪面正中央(浪高儀A1)處達(dá)到最大,此后,由于水體黏性和復(fù)合筒型基礎(chǔ)的阻礙作用,導(dǎo)致波能不斷耗散,波幅呈現(xiàn)出減小的趨勢(shì)(浪高儀C1至E1的變化)繼而因?yàn)閮蓚?cè)邊波在復(fù)合筒型基礎(chǔ)背浪面(浪高儀F1處)疊加導(dǎo)致波浪爬升迅速放大;②對(duì)比波浪單獨(dú)作用和波流共同作用時(shí),隨著流速比的增大,波浪爬升時(shí)程曲線(xiàn)的波峰形狀逐漸變尖瘦,從相應(yīng)的頻譜分析中可以發(fā)現(xiàn),流速比越大,所占的非線(xiàn)性成分越多;③浪高儀C1、E1處波面時(shí)程曲線(xiàn)發(fā)現(xiàn)了對(duì)波浪爬升高度的顯著高階貢獻(xiàn),出現(xiàn)二次波峰(如圖14、15所示),并且隨著流速比增大二次波峰的幅值也相應(yīng)增大,與圖11中的(g)、(i)出現(xiàn)的邊波和波面的變形相對(duì)應(yīng).出現(xiàn)這種現(xiàn)象的原因,可以在 Swan和Sheikh[20]所做的物理試驗(yàn)中找到一些解釋?zhuān)渲胁ㄉ⑸涞钠鹪磁c圍繞基礎(chǔ)的流體循環(huán)有關(guān),并且波流共同作用下,二次波峰的出現(xiàn)要早于波浪單獨(dú)作用.

    圖13 浪高儀A1處波浪爬升時(shí)程和相應(yīng)的頻譜圖

    圖14?浪高儀C1處波浪爬升時(shí)程和相應(yīng)的頻譜圖

    圖15?浪高儀E1處波浪爬升時(shí)程和相應(yīng)的頻譜圖

    圖17為不同波流組合下復(fù)合筒型基礎(chǔ)周?chē)鼽c(diǎn)的爬升分布.其中,為浪高儀至復(fù)合筒中心的水平距離,m;為復(fù)合筒型基礎(chǔ)在其相應(yīng)水深下對(duì)應(yīng)的水面交界處筒的半徑,0.25m水深下=8.07cm.從圖中可以得出:①在0°和36°(即A組和B組)方向上各監(jiān)測(cè)點(diǎn)的波浪爬升值隨著與基礎(chǔ)距離的增大而減小,在D組方向上各監(jiān)測(cè)點(diǎn)的波浪爬升值反而隨著與基礎(chǔ)距離的增大而增大,這種現(xiàn)象可以用伯努利方程來(lái)解釋?zhuān)拷A(chǔ)的水質(zhì)點(diǎn)速度較大,通過(guò)自由液面的伯努利方程可知波面高程隨著與基礎(chǔ)距離的增大而增大.C組、E組和F組方向上波浪爬升值隨著與復(fù)合筒型基礎(chǔ)距離的增大變化不明顯.②流速比增大0.28倍,波浪爬升高度增大約0.55倍,這表明增大海流的速度會(huì)引起爬升高度的顯著增加.

    圖16?浪高儀F1處波浪爬升時(shí)程和相應(yīng)的頻譜圖

    圖17?H/L=0.025時(shí)不同波流組合下各個(gè)浪高儀組的波浪爬升分布

    4?結(jié)論

    (1) 當(dāng)復(fù)合筒型基礎(chǔ)前端波浪回落到波谷位置,波流共同作用時(shí),其前端流體質(zhì)點(diǎn)朝著多個(gè)方向運(yùn)動(dòng),較為紊亂,產(chǎn)生順時(shí)針?lè)较蛐D(zhuǎn)的漩渦;背面水質(zhì)點(diǎn)沿著不同方向傳播,迅速擴(kuò)散,僅有波浪作用時(shí)水質(zhì)點(diǎn)的方向比較一致.

    (2) 波浪單獨(dú)作用和波流聯(lián)合作用時(shí),水體回蕩過(guò)程相似,波浪單獨(dú)作用時(shí),復(fù)合筒型基礎(chǔ)背浪面水體呈“公雞尾”狀,而波流聯(lián)合作用時(shí),復(fù)合筒型基礎(chǔ)背浪面水體呈“雙髻鯊”狀;波浪爬升的最大值都出現(xiàn)在復(fù)合筒型基礎(chǔ)迎浪面正中心即0°(A1)附近,波浪爬升的最小值約在108°(D1).

    (3) 對(duì)比波浪單獨(dú)作用和波流共同作用時(shí)發(fā)現(xiàn),隨著流速比的增大,波浪爬升時(shí)程曲線(xiàn)的波峰形狀逐漸變尖瘦,從相應(yīng)的頻譜分析中可以發(fā)現(xiàn),隨著流速比的增大,非線(xiàn)性成分所占的比例也隨之增大.

    (4) 復(fù)合筒型基礎(chǔ)在A組(0°)和B組(36°)時(shí)的波浪爬升隨著與復(fù)合筒型基礎(chǔ)距離的變化與在D組(108°)時(shí)的波浪爬升變化趨勢(shì)呈現(xiàn)出相反的規(guī)律.而在C組(72°)、E組(144°)和F組(180°)的爬升情況隨著與復(fù)合筒型基礎(chǔ)距離的增大變化并不明顯.流速比由0.46增大到0.59時(shí),波浪爬升高度增大約0.55倍.

    [1] 王月普. 風(fēng)力發(fā)電現(xiàn)狀與發(fā)展趨勢(shì)分析[J]. 電力設(shè)備管理,2020(11):21-22.

    Wang Yuepu. Analysis of current situation and develop-ment trend of wind power generation[J]. Power Equip-ment Management,2020(11):21-22(in Chinese).

    [2] 曲?罡. 海上風(fēng)電筒型基礎(chǔ)結(jié)構(gòu)設(shè)計(jì)研究[D]. 天津:天津大學(xué),2010.

    Qu Gang. The Structural Design Study of Bucket Foundation in Offshore Wind Power[D]. Tianjin:Tianjin University,2010(in Chinese).

    [3] Lian Jijian,Ding Hongyan,Zhang Puyang,et al. Design of large-scale prestressing bucket foundation for offshore wind turbines[J]. Transactions of Tianjin University,2012,18(2):79-84.

    [4] Liu Meimei,Yang Min,Wang Haijun. Bearing behavior of wide shallow bucket foundation for offshore wind turbines in drained silty sand[J]. Ocean Engineering,2014,82:169-179.

    [5] Ding Hongyan,Liu Yonggang,Zhang Puyang,et al. Model tests on the bearing capacity of wide-shallow composite bucket foundations for offshore wind turbines in clay[J]. Ocean Engineering,2015,103:114-122.

    [6] Ding Hongyan,Hu Ruiqi,Zhang Puyang,et al. Load bearing behaviors of composite bucket foundations for offshore wind turbines on layered soil under combined loading[J]. Ocean Engineering,2020,198:106997.

    [7] 張浦陽(yáng),張?志,丁紅巖,等. 筒型基礎(chǔ)預(yù)應(yīng)力過(guò)渡段結(jié)構(gòu)傳力特性研究[J]. 太陽(yáng)能學(xué)報(bào),2017,38(11):3144-3150.

    Zhang Puyang,Zhang Zhi,Ding Hongyang,et al. Research on force transmission chatacteristics of prestressed transition structure of bucket foundation[J]. Journal of Solar Energy,2017,38(11):3144-3150(in Chinese).

    [8] 于通順,練繼建,齊?越,等. 復(fù)合筒型風(fēng)電基礎(chǔ)單向流局部沖刷試驗(yàn)研究[J]. 巖土力學(xué),2015,36(4):1015-1020.

    Yu Tongshun,Lian Jijian,Qi Yue,et al. Experimental study of the local scour around the composite bucket foundations of wind turbines under unidirectional current[J]. Rock and Soil Mechanics,2015,36(4):1015-1020(in Chinese).

    [9] 于通順,張舒博,章譽(yù)天,等. 復(fù)合筒基細(xì)砂地基波流沖刷過(guò)程及最大沖深研究[J]. 太陽(yáng)能學(xué)報(bào),2021,42(8):434-439.

    Yu Tongshun,Zhang Shubo,Zhang Yutian,et al. Study on scour process and maximum scour depth of fine sand composite bucket foundation induced by wave and current action[J]. Journal of Solar Energy,2021,42(8):434-439(in Chinese).

    [10] Li Jinxuan,Wang Zhanhang,Liu Shuxue. Experimen-tal study of interactions between multi-directional focused wave and vertical circular cylinder,Part I:Wave run-up[J]. Coastal Engineering,2012,64(6):151-160.

    [11]單鐵兵,楊建民,李?欣,等. 水流對(duì)立柱周?chē)ɡ伺郎匦缘挠绊慬J]. 上海交通大學(xué)學(xué)報(bào),2014,48(1):116-124,131.

    Shan Tiebing,Yang Jianmin,Li Xin,et al. Current efforts on wave run-up characteristics around column[J]. Journal of Shanghai Jiao Tong University,2014,48(1):116-124,131(in Chinese).

    [12] Bonakdar L,Oumeraci H,Etemad-Shahidi A,et al. Run-up on vertical piles due to regular waves:Small-scale model tests and prediction formulae[J]. Coastal Engineering,2016,118:1-11.

    [13] Cheng M H,Hwang R R,Hsieh C M. Numerical study on the transformation of an internal solitary wave propagating across a vertical cylinder[J]. Applied Ocean Research,2020,95:102016.

    [14] Cao Hongjian,Wan Decheng. Benchmark computations of wave run-up on single cylinder and four cylinders by naoe-FOAM-SJTU solver[J]. Applied Ocean Research,2017,65:327-337.

    [15] Lin Yuhsien,Chen Jingfu,Lu Poying. A CFD model for simulating wave run-ups and wave loads in case of different wind turbine foundations influenced by nonlinear waves[J]. Ocean Engineering,2017,129:428-440.

    [16] Lin Y H,Lu P Y,Lin C W,et al. Numerical simulation of maximum wave loads and run-up heights on offshore wind turbine foundations influenced by the instability of bichromatic wave groups[J]. Marine Structures,2019,67:102648.

    [17] Yu Tongshun,Zhao Zishuai,Shi Zhongyu,et al. Experimental investigation of wave load and run-up on the composite bucket foundation influenced by regular waves[J]. Journal of Ocean University of China,2021,20(2):271-284.

    [18]李玉成. 波流共同作用下的波浪要素[J]. 大連工學(xué)院學(xué)報(bào),1982(4):81-89.

    Li Yucheng. Wave elements under the combined action of wave and current[J]. Journal of the Dalian University of Technology,1982(4):81-89(in Chinese).

    [19] Sun L,Zang J,Chen L F,et al. Regular waves onto a truncated circular column:A comparison of experiments and simulations[J]. Applied Ocean Research,2016,59(1):650-662.

    [20] Swan C,Sheikh R. The interaction between steep waves and a surface-piercing column[J]. Philos Trans A Math Phys Eng,2015,373(2033):41251-41257.

    Numerical Simulation of the Effect of Current on Wave Run-Up Around Composite Bucket Foundation

    Yu Tongshun1,Song Haoyang1,Liu Meimei2,Zhang Yutian3

    (1. College of Engineering,Ocean University of China,Qingdao 266100,China;2. CCCC Fourth Harbor Engineering Institute Co.,Ltd.,Guangzhou 510230,China;3. Dianjiang CountyWater Conservancy Bureau of Chongqing City,Chongqing 408300,China)

    Composite bucket foundation(CBF)is a new type of offshore wind turbine foundation. When wave and current interact with its large-scale arc transition structure,complex wave run-up occurs,which heavily affects the safety of the platforms. Thus,it is necessary to study the wave run-up effect on CBF. First,a numerical wave-current tank was established using the RNG-model. Secondly,the accuracy of the numerical model was verified based on the experimental data. Finally,the law of the wave run-up around the CBF under the action of wave and current was explored,and the contribution of current to wave run-up under different velocity ratioscwwas revealed. The results show that:①Under the action of wave alone,when the water in front of the CBF falls back to the lowest point,the motion trajectory of the water particles is relatively regular,and the vector directions are inclined downward. However,under the combined action of wave and current,a vortex is formed at the front end of the CBF,and the water particles on the back will diffuse in different directions. ②When the wave acts alone,the water on the back of the CBF is a“cock tail”shape,while when the wave and current act together,the water on the back of the CBF is a “hammerhead shark”shape. ③When the wave steepness parameter is 0.025 and the velocity ratio increases from 0.46 to 0.59,the wave run-up height increases about 0.55 times.

    composite bucket foundation(CBF);wave run-up;wave-current interaction;numerical simulation

    10.11784/tdxbz202110024

    TV92

    A

    0493-2137(2022)12-1309-09

    2021-10-26;

    2021-12-09.

    于通順(1986—??),男,博士,副教授.

    于通順,tshyu707@ouc.edu.cn.

    泰山學(xué)者工程專(zhuān)項(xiàng)經(jīng)費(fèi)資助項(xiàng)目(ts20190914);國(guó)家自然科學(xué)基金資助項(xiàng)目(52071304,51509230);山東省重點(diǎn)科研開(kāi)發(fā)項(xiàng)目(2019GHY112044).

    Supported by the Special Funds for the Taishan Scholars Project(No.ts20190914),the National Natural Science Foundation of China (No.52071304,No.51509230),the Key Research and Development Program of Shandong Province(No.2019GHY112044).

    (責(zé)任編輯:樊素英)

    猜你喜歡
    波流波面水槽
    波流耦合下樁周珊瑚砂沖刷機(jī)理研究
    可升降折疊的飲水機(jī)水槽
    波流聯(lián)合作用下海上輸油漂浮軟管動(dòng)力響應(yīng)分析
    基于恒定陡度聚焦波模型的分析與討論
    水道港口(2020年6期)2020-02-22 11:33:50
    偶感
    多普勒效應(yīng)中觀察者接收頻率的計(jì)算
    為什么水槽管要做成彎曲狀
    淺談光的干涉和衍射的區(qū)別和聯(lián)系
    中文信息(2018年2期)2018-05-30 11:45:10
    槽道內(nèi)渦波流場(chǎng)展向渦的分布特征
    要挑好水槽,就看這里了!
    幸福(2016年6期)2016-12-01 03:08:13
    熟妇人妻久久中文字幕3abv| 免费搜索国产男女视频| 中文字幕人成人乱码亚洲影| 少妇熟女aⅴ在线视频| 成人亚洲精品av一区二区| 女警被强在线播放| 国产一区二区激情短视频| 99国产极品粉嫩在线观看| 一二三四在线观看免费中文在| xxx96com| 色av中文字幕| www国产在线视频色| 亚洲五月天丁香| 嫁个100分男人电影在线观看| 亚洲伊人色综图| 免费在线观看影片大全网站| 丝袜美足系列| 正在播放国产对白刺激| 久久久精品国产亚洲av高清涩受| 日本免费a在线| 美女高潮到喷水免费观看| 精品国产美女av久久久久小说| 美女午夜性视频免费| 波多野结衣av一区二区av| 久久久久亚洲av毛片大全| av中文乱码字幕在线| 人人妻,人人澡人人爽秒播| av天堂在线播放| 亚洲第一av免费看| 亚洲电影在线观看av| 少妇熟女aⅴ在线视频| 91精品国产国语对白视频| 国产成人欧美在线观看| 性欧美人与动物交配| aaaaa片日本免费| 一级作爱视频免费观看| 亚洲aⅴ乱码一区二区在线播放 | 亚洲精品美女久久av网站| 日韩成人在线观看一区二区三区| 不卡一级毛片| 亚洲avbb在线观看| 一区二区三区精品91| 视频区欧美日本亚洲| 欧美精品啪啪一区二区三区| www.精华液| 1024视频免费在线观看| 黄色视频,在线免费观看| 宅男免费午夜| 精品一品国产午夜福利视频| 久久久久亚洲av毛片大全| 熟妇人妻久久中文字幕3abv| 高清在线国产一区| 精品久久蜜臀av无| 中亚洲国语对白在线视频| 亚洲男人天堂网一区| 成人国产一区最新在线观看| 久久精品亚洲精品国产色婷小说| 精品午夜福利视频在线观看一区| 成人手机av| 久久中文字幕一级| 国产精品久久久久久精品电影 | 精品久久蜜臀av无| 男人的好看免费观看在线视频 | 久久人人爽av亚洲精品天堂| 制服诱惑二区| 免费女性裸体啪啪无遮挡网站| 亚洲视频免费观看视频| 一进一出抽搐gif免费好疼| 国产精品电影一区二区三区| 很黄的视频免费| 中亚洲国语对白在线视频| videosex国产| 1024香蕉在线观看| 巨乳人妻的诱惑在线观看| 亚洲五月色婷婷综合| www日本在线高清视频| 国产精品野战在线观看| 大型av网站在线播放| 欧美成人午夜精品| 国产精品乱码一区二三区的特点 | 丝袜在线中文字幕| 亚洲一区中文字幕在线| 天天一区二区日本电影三级 | 日韩精品免费视频一区二区三区| 超碰成人久久| 亚洲 欧美一区二区三区| 国产区一区二久久| 变态另类丝袜制服| 亚洲精品中文字幕在线视频| 免费在线观看影片大全网站| 精品不卡国产一区二区三区| 国产精品国产高清国产av| 久久久久久久久中文| 黄色丝袜av网址大全| 国产区一区二久久| 99riav亚洲国产免费| 两个人视频免费观看高清| 丰满的人妻完整版| 天天躁夜夜躁狠狠躁躁| 午夜两性在线视频| 午夜激情av网站| 日本在线视频免费播放| 欧美日韩福利视频一区二区| 制服丝袜大香蕉在线| 黄片小视频在线播放| 国产男靠女视频免费网站| 此物有八面人人有两片| 亚洲午夜理论影院| 亚洲国产精品久久男人天堂| 国产精品亚洲av一区麻豆| 中文字幕人妻熟女乱码| or卡值多少钱| 亚洲精品美女久久久久99蜜臀| 精品国产一区二区三区四区第35| 日韩欧美国产一区二区入口| 又黄又爽又免费观看的视频| 久久国产亚洲av麻豆专区| 最新在线观看一区二区三区| 国产精品亚洲美女久久久| 天天添夜夜摸| 91麻豆精品激情在线观看国产| 日韩成人在线观看一区二区三区| 老司机靠b影院| 国产精品一区二区精品视频观看| 欧美性长视频在线观看| 99国产精品一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 精品久久久久久成人av| 国内毛片毛片毛片毛片毛片| 色综合亚洲欧美另类图片| 免费在线观看影片大全网站| 国产精品亚洲美女久久久| 国内精品久久久久精免费| 人人澡人人妻人| 国产精品 欧美亚洲| 免费不卡黄色视频| e午夜精品久久久久久久| 久99久视频精品免费| 最近最新中文字幕大全免费视频| 99精品欧美一区二区三区四区| 国产亚洲精品第一综合不卡| 亚洲三区欧美一区| 色老头精品视频在线观看| www日本在线高清视频| 国产av一区在线观看免费| 亚洲欧美精品综合久久99| 97超级碰碰碰精品色视频在线观看| 久久九九热精品免费| 午夜福利在线观看吧| 亚洲国产精品成人综合色| 日韩av在线大香蕉| or卡值多少钱| 久久久国产精品麻豆| 亚洲第一av免费看| 久久欧美精品欧美久久欧美| 999精品在线视频| 久热爱精品视频在线9| 美女国产高潮福利片在线看| 在线十欧美十亚洲十日本专区| 国产99久久九九免费精品| 免费久久久久久久精品成人欧美视频| 大香蕉久久成人网| 母亲3免费完整高清在线观看| 免费人成视频x8x8入口观看| 男女之事视频高清在线观看| 国产一区二区三区在线臀色熟女| 深夜精品福利| 88av欧美| 19禁男女啪啪无遮挡网站| 日韩欧美免费精品| 香蕉丝袜av| 大型av网站在线播放| 一个人观看的视频www高清免费观看 | 国产又爽黄色视频| 99re在线观看精品视频| 亚洲视频免费观看视频| 国产亚洲欧美在线一区二区| 亚洲精华国产精华精| 国产麻豆69| 亚洲avbb在线观看| 精品第一国产精品| 麻豆一二三区av精品| 在线观看舔阴道视频| 欧美 亚洲 国产 日韩一| 久久国产精品影院| 亚洲午夜精品一区,二区,三区| 欧美绝顶高潮抽搐喷水| 亚洲专区国产一区二区| 夜夜看夜夜爽夜夜摸| 此物有八面人人有两片| 亚洲av成人不卡在线观看播放网| 午夜福利一区二区在线看| 国产高清videossex| 十八禁人妻一区二区| 18禁黄网站禁片午夜丰满| 久久青草综合色| 亚洲第一电影网av| 国产单亲对白刺激| 99国产精品99久久久久| 少妇熟女aⅴ在线视频| 久久国产乱子伦精品免费另类| 久久香蕉国产精品| 操出白浆在线播放| 99精品久久久久人妻精品| 国产精品一区二区在线不卡| 日韩一卡2卡3卡4卡2021年| 色av中文字幕| 欧洲精品卡2卡3卡4卡5卡区| 欧美绝顶高潮抽搐喷水| 国产高清视频在线播放一区| 亚洲中文av在线| 日本五十路高清| 亚洲精品粉嫩美女一区| 国产成人精品无人区| 两性夫妻黄色片| 国产人伦9x9x在线观看| 91九色精品人成在线观看| 免费在线观看视频国产中文字幕亚洲| 十八禁网站免费在线| 久久久久久久午夜电影| 免费女性裸体啪啪无遮挡网站| 亚洲欧美激情综合另类| 国产熟女xx| 亚洲人成网站在线播放欧美日韩| 亚洲色图av天堂| 午夜福利一区二区在线看| 亚洲伊人色综图| 国产一区二区三区在线臀色熟女| 欧美成狂野欧美在线观看| 久久青草综合色| 大型黄色视频在线免费观看| 免费在线观看完整版高清| videosex国产| videosex国产| 亚洲一区高清亚洲精品| 欧美日韩福利视频一区二区| 女人被狂操c到高潮| 中文字幕精品免费在线观看视频| 性欧美人与动物交配| 人妻丰满熟妇av一区二区三区| 午夜激情av网站| 亚洲欧美日韩无卡精品| 日韩av在线大香蕉| 国产成人av教育| 在线观看舔阴道视频| 91成人精品电影| 岛国视频午夜一区免费看| 日韩av在线大香蕉| 熟女少妇亚洲综合色aaa.| 国语自产精品视频在线第100页| 成人三级黄色视频| 久久人妻福利社区极品人妻图片| 免费在线观看亚洲国产| 国产精品乱码一区二三区的特点 | 制服丝袜大香蕉在线| 亚洲国产中文字幕在线视频| 亚洲欧美精品综合一区二区三区| 日韩成人在线观看一区二区三区| 国产成人精品在线电影| 女人精品久久久久毛片| 国产区一区二久久| 91九色精品人成在线观看| 国产午夜福利久久久久久| 老熟妇仑乱视频hdxx| 欧美日本中文国产一区发布| 在线观看免费日韩欧美大片| 极品人妻少妇av视频| 久久热在线av| 国产成人欧美在线观看| 熟妇人妻久久中文字幕3abv| 欧美+亚洲+日韩+国产| 搞女人的毛片| 天天躁夜夜躁狠狠躁躁| 久久精品aⅴ一区二区三区四区| 黄片小视频在线播放| 黄片大片在线免费观看| 免费看十八禁软件| 91老司机精品| 欧美国产精品va在线观看不卡| 日韩欧美在线二视频| 欧美日韩黄片免| 国产免费av片在线观看野外av| 国产精品99久久99久久久不卡| 亚洲成人国产一区在线观看| 中亚洲国语对白在线视频| 日本三级黄在线观看| 日韩精品中文字幕看吧| 久久久久精品国产欧美久久久| 不卡一级毛片| 宅男免费午夜| 99香蕉大伊视频| 久久久久精品国产欧美久久久| 亚洲一区二区三区色噜噜| 色在线成人网| 国产精品久久视频播放| 欧美国产日韩亚洲一区| 国产区一区二久久| 午夜激情av网站| 亚洲久久久国产精品| 亚洲一卡2卡3卡4卡5卡精品中文| 色精品久久人妻99蜜桃| 亚洲 欧美一区二区三区| 久久久久亚洲av毛片大全| 美女大奶头视频| 亚洲一码二码三码区别大吗| 国产aⅴ精品一区二区三区波| 一区二区三区激情视频| or卡值多少钱| 成人三级做爰电影| 一二三四在线观看免费中文在| 老汉色av国产亚洲站长工具| 国产麻豆69| 色老头精品视频在线观看| 一进一出好大好爽视频| 国产麻豆成人av免费视频| 免费无遮挡裸体视频| 免费高清视频大片| 国产成人一区二区三区免费视频网站| 美女国产高潮福利片在线看| 国产成年人精品一区二区| 99久久精品国产亚洲精品| 九色亚洲精品在线播放| 亚洲久久久国产精品| 国产97色在线日韩免费| 日韩中文字幕欧美一区二区| 午夜福利免费观看在线| svipshipincom国产片| 黄色 视频免费看| 免费女性裸体啪啪无遮挡网站| 国产aⅴ精品一区二区三区波| 此物有八面人人有两片| 夜夜爽天天搞| 精品国内亚洲2022精品成人| 亚洲精品在线观看二区| 老鸭窝网址在线观看| 国产日韩一区二区三区精品不卡| 午夜精品国产一区二区电影| 亚洲av第一区精品v没综合| 成人18禁高潮啪啪吃奶动态图| 精品人妻1区二区| 1024视频免费在线观看| 真人做人爱边吃奶动态| 99riav亚洲国产免费| 亚洲精品av麻豆狂野| 精品久久久久久,| 给我免费播放毛片高清在线观看| 熟妇人妻久久中文字幕3abv| 国产精品,欧美在线| 日本vs欧美在线观看视频| 又紧又爽又黄一区二区| 久久久国产精品麻豆| 欧美中文综合在线视频| 曰老女人黄片| 女生性感内裤真人,穿戴方法视频| 成年女人毛片免费观看观看9| 又紧又爽又黄一区二区| 午夜福利免费观看在线| 欧美人与性动交α欧美精品济南到| 老司机深夜福利视频在线观看| 国产私拍福利视频在线观看| 日本vs欧美在线观看视频| 亚洲中文字幕日韩| 国产精品av久久久久免费| 久久午夜综合久久蜜桃| 国产精品影院久久| 黄片大片在线免费观看| 女人爽到高潮嗷嗷叫在线视频| 欧美激情 高清一区二区三区| 国产精品秋霞免费鲁丝片| 国产精品98久久久久久宅男小说| 视频区欧美日本亚洲| 亚洲精品国产色婷婷电影| 日韩精品青青久久久久久| 国产精品 国内视频| 一级毛片精品| av欧美777| x7x7x7水蜜桃| 亚洲人成网站在线播放欧美日韩| 国内毛片毛片毛片毛片毛片| 人妻丰满熟妇av一区二区三区| 国产又色又爽无遮挡免费看| 日本一区二区免费在线视频| 国产精品一区二区精品视频观看| 亚洲人成77777在线视频| av天堂久久9| 国产野战对白在线观看| 大香蕉久久成人网| 免费在线观看视频国产中文字幕亚洲| 亚洲av成人av| 国产单亲对白刺激| 一区二区三区高清视频在线| 麻豆av在线久日| 亚洲国产日韩欧美精品在线观看 | 九色亚洲精品在线播放| 国产精品爽爽va在线观看网站 | 操美女的视频在线观看| 日韩中文字幕欧美一区二区| 制服丝袜大香蕉在线| 日日夜夜操网爽| 国产av一区二区精品久久| 97超级碰碰碰精品色视频在线观看| 动漫黄色视频在线观看| 高清黄色对白视频在线免费看| 可以在线观看毛片的网站| 国产精品久久久人人做人人爽| 性少妇av在线| 动漫黄色视频在线观看| 最近最新中文字幕大全免费视频| 亚洲中文日韩欧美视频| 久久中文看片网| 成人三级做爰电影| 乱人伦中国视频| tocl精华| 两性午夜刺激爽爽歪歪视频在线观看 | 一区二区三区精品91| 精品一区二区三区av网在线观看| 丁香欧美五月| 三级毛片av免费| 国产精品98久久久久久宅男小说| 久久人妻福利社区极品人妻图片| 女性生殖器流出的白浆| 国产精品 欧美亚洲| 亚洲自拍偷在线| 国产av又大| 亚洲精品一卡2卡三卡4卡5卡| 亚洲av第一区精品v没综合| 色播在线永久视频| 深夜精品福利| 日本在线视频免费播放| 午夜影院日韩av| 国产国语露脸激情在线看| 精品国产超薄肉色丝袜足j| 色播亚洲综合网| 亚洲国产欧美网| 亚洲国产精品sss在线观看| 久久精品国产亚洲av高清一级| 看片在线看免费视频| 国产私拍福利视频在线观看| 熟女少妇亚洲综合色aaa.| 91成人精品电影| 国产麻豆69| 亚洲国产高清在线一区二区三 | 黄色视频不卡| 淫妇啪啪啪对白视频| 久久中文看片网| 国产成人啪精品午夜网站| 一级a爱片免费观看的视频| e午夜精品久久久久久久| 午夜精品国产一区二区电影| 制服诱惑二区| 女同久久另类99精品国产91| 国产高清有码在线观看视频 | 黄色视频,在线免费观看| 久久国产亚洲av麻豆专区| 亚洲片人在线观看| 欧美一级毛片孕妇| 久久婷婷成人综合色麻豆| 午夜福利视频1000在线观看 | 最近最新免费中文字幕在线| 老熟妇乱子伦视频在线观看| 国产亚洲精品久久久久5区| 黄色a级毛片大全视频| 18禁国产床啪视频网站| 亚洲国产中文字幕在线视频| 久久国产精品男人的天堂亚洲| 性欧美人与动物交配| 波多野结衣高清无吗| bbb黄色大片| 精品电影一区二区在线| 在线观看免费日韩欧美大片| 精品电影一区二区在线| 久久久久九九精品影院| 亚洲精品一卡2卡三卡4卡5卡| 亚洲黑人精品在线| 国产精品国产高清国产av| 日本a在线网址| 成人永久免费在线观看视频| 久久久精品欧美日韩精品| 日日爽夜夜爽网站| 久久精品影院6| av天堂久久9| 国产私拍福利视频在线观看| 成人18禁在线播放| 国产精品香港三级国产av潘金莲| 可以免费在线观看a视频的电影网站| av天堂在线播放| 免费在线观看黄色视频的| 国产亚洲av高清不卡| 大香蕉久久成人网| 一区二区三区激情视频| 久久精品91无色码中文字幕| 老司机在亚洲福利影院| 一级片免费观看大全| 欧美性长视频在线观看| 精品免费久久久久久久清纯| 熟妇人妻久久中文字幕3abv| 免费女性裸体啪啪无遮挡网站| 老司机深夜福利视频在线观看| 久久中文字幕人妻熟女| 搡老熟女国产l中国老女人| 丝袜美腿诱惑在线| 国产成人精品在线电影| 精品久久久久久久人妻蜜臀av | 最新在线观看一区二区三区| 欧美日韩黄片免| 真人一进一出gif抽搐免费| 成人国产综合亚洲| 亚洲五月天丁香| 成人国产综合亚洲| 少妇裸体淫交视频免费看高清 | 少妇裸体淫交视频免费看高清 | 高清黄色对白视频在线免费看| 亚洲一区中文字幕在线| 色播亚洲综合网| 又紧又爽又黄一区二区| 又大又爽又粗| 中文亚洲av片在线观看爽| 青草久久国产| 99久久久亚洲精品蜜臀av| 别揉我奶头~嗯~啊~动态视频| 亚洲电影在线观看av| 欧美成狂野欧美在线观看| 一进一出抽搐gif免费好疼| 女警被强在线播放| 乱人伦中国视频| 免费不卡黄色视频| 午夜影院日韩av| 人妻久久中文字幕网| 男女下面进入的视频免费午夜 | 两人在一起打扑克的视频| 国产野战对白在线观看| 亚洲欧美精品综合久久99| 午夜a级毛片| 黄色丝袜av网址大全| 亚洲久久久国产精品| 国产精品野战在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 国产激情欧美一区二区| 亚洲av日韩精品久久久久久密| 91在线观看av| 最近最新中文字幕大全免费视频| 欧美日韩精品网址| 久久久久九九精品影院| 妹子高潮喷水视频| 欧美日韩福利视频一区二区| 99国产极品粉嫩在线观看| 久久久久久亚洲精品国产蜜桃av| 一级a爱片免费观看的视频| 午夜久久久在线观看| 国产91精品成人一区二区三区| 黄色a级毛片大全视频| 美女高潮到喷水免费观看| 亚洲av片天天在线观看| 国内毛片毛片毛片毛片毛片| 激情在线观看视频在线高清| 这个男人来自地球电影免费观看| 精品国产超薄肉色丝袜足j| 无人区码免费观看不卡| 亚洲国产毛片av蜜桃av| 美女扒开内裤让男人捅视频| 国产私拍福利视频在线观看| 成人欧美大片| 欧美日韩亚洲综合一区二区三区_| 少妇熟女aⅴ在线视频| www.www免费av| 欧美一区二区精品小视频在线| 桃色一区二区三区在线观看| 午夜a级毛片| 久久久久久久久免费视频了| 久久精品人人爽人人爽视色| 欧美性长视频在线观看| 亚洲五月天丁香| 免费在线观看影片大全网站| av片东京热男人的天堂| 中文字幕人成人乱码亚洲影| 国产av一区在线观看免费| 天堂√8在线中文| 亚洲国产欧美网| 亚洲黑人精品在线| 精品一区二区三区视频在线观看免费| 日本一区二区免费在线视频| 亚洲国产精品999在线| 亚洲国产精品久久男人天堂| 老汉色av国产亚洲站长工具| 美国免费a级毛片| 电影成人av| 免费在线观看影片大全网站| 国产97色在线日韩免费| 久久精品国产亚洲av香蕉五月| 久久影院123| 亚洲熟妇熟女久久| 久久精品亚洲精品国产色婷小说| 日本 欧美在线| 黄色 视频免费看| 国产亚洲欧美98| 婷婷六月久久综合丁香| 精品国内亚洲2022精品成人| 叶爱在线成人免费视频播放| 国产精品国产高清国产av| 后天国语完整版免费观看| 亚洲伊人色综图| 色尼玛亚洲综合影院| 亚洲精品国产精品久久久不卡| 99国产精品一区二区三区| 丝袜人妻中文字幕| 97碰自拍视频| 亚洲精品国产一区二区精华液| 国产精品二区激情视频| 国产亚洲精品一区二区www| 午夜福利在线观看吧| 老熟妇仑乱视频hdxx| 成人亚洲精品av一区二区| 国产欧美日韩精品亚洲av| 国产精品野战在线观看| 男女下面插进去视频免费观看|