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

    非周期波浪與直墻作用的非線性數(shù)值研究1)

    2017-11-11 01:54:41張崇偉寧德志
    力學(xué)學(xué)報(bào) 2017年5期
    關(guān)鍵詞:直墻入射波波浪

    李 翔 張崇偉 寧德志 蘇 朋

    (大連理工大學(xué)海岸和近海國(guó)家重點(diǎn)實(shí)驗(yàn)室,大連116024)

    非周期波浪與直墻作用的非線性數(shù)值研究1)

    李 翔 張崇偉 寧德志2)蘇 朋

    (大連理工大學(xué)海岸和近海國(guó)家重點(diǎn)實(shí)驗(yàn)室,大連116024)

    基于時(shí)域高階邊界元方法,建立了完全非線性二維數(shù)值波浪水槽,對(duì)非周期波浪與直墻的相互作用問(wèn)題進(jìn)行了模擬和研究.自由表面滿足完全非線性自由水面運(yùn)動(dòng)學(xué)和動(dòng)力學(xué)邊界條件,采用混合歐拉--拉格朗日方法追蹤瞬時(shí)自由面流體質(zhì)點(diǎn),采用四階Runge-Kutta法對(duì)下一時(shí)間步的波面和自由面速度勢(shì)進(jìn)行更新.采用加速度式法求解直墻表面速度勢(shì)的時(shí)間導(dǎo)數(shù),對(duì)瞬時(shí)物體濕表面上的水動(dòng)力壓強(qiáng)積分,得到作用在物體上的瞬時(shí)波浪力.首先,將全非線性與Serre-Green-Naghdi(SGN)模型的結(jié)果進(jìn)行了對(duì)比分析,發(fā)現(xiàn)對(duì)于大幅值雙入射波問(wèn)題,僅滿足弱色散關(guān)系的SGN模型大大低估了最大波浪爬高;其次,研究了雙入射波與直墻的非線性作用問(wèn)題,發(fā)現(xiàn)線性預(yù)報(bào)對(duì)波浪最大爬高有較大低估,而波浪的非線性成分不只導(dǎo)致了自由面爬高的異常增大,也引起了局部自由面的高頻振蕩,該物理過(guò)程中,直墻所受的波浪載荷,也展示出了與波浪爬高相似的非線性特性;最后,對(duì)波浪爬升和波浪力的時(shí)間歷程進(jìn)行了頻譜分析,發(fā)現(xiàn)入射主頻波的部分能量傳遞給了更高頻的波浪成分,反映出該問(wèn)題具有典型的非線性特性.

    非線性波浪,波浪爬高,時(shí)域模擬,高階邊界元,GN理論

    引言

    近年來(lái),由于大型海洋結(jié)構(gòu)的廣泛應(yīng)用,波浪與海洋結(jié)構(gòu)相互作用的問(wèn)題越來(lái)越被人們所重視,并在近年來(lái)得到了廣泛的研究[13].其中,直墻式結(jié)構(gòu)是一種常見(jiàn)的海岸工程結(jié)構(gòu),其安全性和造價(jià)很大程度上取決于作用在其上的極限波浪載荷.工程實(shí)際中,人們多通過(guò)研究波浪沿直墻的爬高,來(lái)預(yù)報(bào)該極限波浪載荷.

    針對(duì)波浪爬高問(wèn)題,目前人們圍繞連續(xù)波浪、孤立波與直墻的相互作用已經(jīng)做了大量研究.例如,對(duì)于連續(xù)波,Isaacson和Cheung[4]提出了二階時(shí)域模擬的理論模型,研究了波浪在直壁物體上的爬高問(wèn)題.Grilli等[5]基于全非線性勢(shì)流理論,建立了三維數(shù)值水槽,利用邊界元方法來(lái)模擬波浪爬高.Jamois等[6]使用基于高階Boussinesq方程的有限差分模型,對(duì)三維波浪在豎直板的爬高問(wèn)題進(jìn)行了研究.Fuhrman和Madsen[7]在高階Boussinesq模型中引入了一種插值技術(shù)來(lái)捕捉非線性自由面沿壁面的爬高.Pelinovsky等[8]求解了直墻附近非線性淺水方程,并給出了直墻附近出現(xiàn)畸形波的概率分布.孫英偉[9]以Naiver-Stokes方程為基本控制方程,采用流體體積和大渦模擬方法,研究了規(guī)則波作用于直墻過(guò)程中的流場(chǎng)形態(tài)和波浪力.Chatjigeorgiou和Molin[10]采用了改進(jìn)型Boussinesq方程來(lái)研究波浪與直墻的相互作用問(wèn)題,發(fā)現(xiàn)波浪爬升與入射波和反射波的三階非線性作用密切相關(guān).

    對(duì)于孤立波,Su和 Mirie[11]通過(guò)攝動(dòng)展開(kāi)法推導(dǎo)了孤立波沿直墻爬高的三階解析解.Fenton和Rienecker[12]采用譜方法研究了孤立波在直墻上的碰撞問(wèn)題,發(fā)現(xiàn)該過(guò)程對(duì)反射波的波高和波速都產(chǎn)生了影響,而考慮三階非線性效應(yīng)的模擬結(jié)果可準(zhǔn)確預(yù)測(cè)最大波浪爬高,但無(wú)法準(zhǔn)確預(yù)測(cè)反射波相位的改變.Kim等[13]最早使用二維邊界元方法,模擬了考慮完全非線性自由面條件的孤立波沿直墻的爬高問(wèn)題.Power和Chwang[14]對(duì)孤立波的直墻反射問(wèn)題,給出了Boussinesq方程的解析解和數(shù)值解,也發(fā)現(xiàn)了反射孤立波的相位延遲現(xiàn)象.Cooker等[15]使用二維時(shí)域全非線性邊界元方法研究了大幅孤立波與直墻相互作用,發(fā)現(xiàn)孤立波碰撞直墻后能量耗散給一群色散波,降低了反射波波高和波速;同時(shí),直墻上波浪下降要比其爬升過(guò)程消耗更多時(shí)間.Maiti和Sen[16]采用時(shí)域全非線性邊界元方法分析了孤立波與直墻非線性作用問(wèn)題,他們發(fā)現(xiàn),比起傾斜墻壁,垂直墻壁上波浪作用的非線性效應(yīng)更顯著.Liu和Al-Banaa[17]對(duì)孤立波在下方有開(kāi)孔的直墻上的爬高進(jìn)行了物理模型實(shí)驗(yàn).Jian等[18]采用光滑粒子水動(dòng)力學(xué)方法,對(duì)包含波浪破碎過(guò)程的孤立波爬高問(wèn)題進(jìn)行了模擬.Park和Cox[19]建立了近岸波浪和風(fēng)暴潮爬高的經(jīng)驗(yàn)公式.Chambarel等[20]采用考慮全非線性自由表面的二維邊界元方法模擬了兩相同孤立波的對(duì)撞問(wèn)題(等效于孤立波與直墻的對(duì)撞問(wèn)題),觀測(cè)到了射流現(xiàn)象.宣瑞韜等[21]通過(guò)物理試驗(yàn)研究了雙孤立波的直墻爬高問(wèn)題,著重分析了兩波相對(duì)波高和波距對(duì)波浪爬高的影響.

    近期Carbone等[22]指出,某些特殊波浪(如非周期波浪)也值得關(guān)注.非周期波浪與直墻作用時(shí)可能會(huì)產(chǎn)生比連續(xù)波浪更大的波浪爬高,導(dǎo)致更大的波浪載荷,從而威脅結(jié)構(gòu)的安全性.為了提高計(jì)算效率,他們?cè)跀?shù)值模擬中采用了基于淺水假設(shè)且僅滿足弱色散關(guān)系的Serre-Green-Naghdi(SGN)[2326]波浪模型.而Cooker等[12]提出波浪爬高問(wèn)題會(huì)伴隨一群色散波.因此,本文將基于全非線性勢(shì)流模型,進(jìn)一步對(duì)文獻(xiàn)[19]中涉及到的非周期波浪與直墻相互作用的問(wèn)題進(jìn)行研究,通過(guò)與SGN模型結(jié)果的對(duì)比,研究SGN模型在該問(wèn)題中的適用性.另外,本文將重點(diǎn)研究直墻結(jié)構(gòu)上雙入射波爬升和載荷的非線性特性.通過(guò)頻譜分析,研究波浪爬升過(guò)程中各頻率組分的能量分布,進(jìn)而分析波與波相互作用的非線性物理特征.

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

    本文將建立如圖1所示的二維數(shù)值波浪水槽.水槽左側(cè)為造波推板,右側(cè)為直墻.水槽長(zhǎng)度為L(zhǎng),平均水深為d,造波板、自由水面、直墻和水底邊界分別用Γw,Γf,Γb,Γd來(lái)表示,流域用V來(lái)表示.

    圖1 波浪水槽示意圖Fig.1 Sketch of the wave fl ume

    建立笛卡爾坐標(biāo)系oxz,其中o位于造波板與自由水面平均位置的交點(diǎn)處,ox指向右側(cè)水平方向,oz豎直向上.假定流體為理想流體且其運(yùn)動(dòng)無(wú)旋,則可采用勢(shì)流理論模型來(lái)描述水槽內(nèi)水體運(yùn)動(dòng),并得到下述初邊值問(wèn)題

    其中,速度勢(shì)?(x,y,z)為標(biāo)量,其空間梯度為流體的速度,η表示自由面波高,g為重力加速度,??/?n表示速度勢(shì)在物面法線方向的偏導(dǎo)數(shù).造波板的速度u(t)設(shè)為

    式中,a和ω分別為波浪的幅值和頻率,T為波浪周期,H(t)為Heaviside函數(shù)

    作用在直墻上的水平波浪力可通過(guò)對(duì)Bernoulli方程得到的壓力進(jìn)行積分求得

    式中?t為速度勢(shì)對(duì)時(shí)間的偏導(dǎo)數(shù).可通過(guò)求解下列邊值問(wèn)題得到

    2 數(shù)值方法

    根據(jù)格林第二定理,方程(1)可轉(zhuǎn)換為下述邊界積分方程[27]

    式中,p={x0,z0}為配置點(diǎn)坐標(biāo),q={x,z}為邊界上一點(diǎn)的坐標(biāo),C為固角系數(shù)[27],Γ為全部流域邊界,G為簡(jiǎn)單格林函數(shù).對(duì)于該二維問(wèn)題,G(p,q)可定義為

    進(jìn)一步采用邊界單元來(lái)離散流域邊界.為了更準(zhǔn)確地描述流域邊界,這里采用曲線單元,即用二次曲線來(lái)逼近真實(shí)的幾何邊界,它比直線單元具有更強(qiáng)的適應(yīng)性.單元上任一點(diǎn)的坐標(biāo)和速度勢(shì)等物理量可通過(guò)如下插值得到

    式中,M為單元上節(jié)點(diǎn)個(gè)數(shù),S用來(lái)表示x,z,φ和?φ/?n等物理量,Sn為單元上第n個(gè)節(jié)點(diǎn)處物理量的值,ξ∈[?1.0,1.0]為固有坐標(biāo),hn表示第n個(gè)形狀函數(shù)

    這里采用三節(jié)點(diǎn)二次高階單元來(lái)離散流域邊界,即M=3.該方法比常數(shù)單元和線性單元具有更高的精度.

    式中,d(e,m)表示第e個(gè)單元上第 m個(gè)節(jié)點(diǎn)的編號(hào),Ne為邊界單元的數(shù)量.J(ξ)是雅可比(Jacobian)行列式,可寫(xiě)為如下形式

    積分方程經(jīng)過(guò)離散后,可通過(guò)求解線性方程組得到未知量.

    在每一個(gè)時(shí)間步,自由面上的速度勢(shì)和其他物面上的速度勢(shì)法向?qū)?shù)是已知的,通過(guò)求解線性方程組(20),可以得到物面上速度勢(shì)和自由面上速度勢(shì)法向?qū)?shù)的值.通過(guò)差分可以得到邊界上速度勢(shì)的切向?qū)?shù),并通過(guò)下列關(guān)系求得速度勢(shì)的空間導(dǎo)數(shù)值

    其中,n={nx,nz}為單位法向量.

    根據(jù)自由面邊界條件,可得到速度勢(shì)和波高的時(shí)間偏導(dǎo)數(shù),可應(yīng)用數(shù)值積分的方法,進(jìn)行時(shí)間步進(jìn),計(jì)算得到下一時(shí)刻自由面上的速度勢(shì)和波高.本文選用四階Runge-Kutta法對(duì)自由面條件進(jìn)行時(shí)間積分.為了敘述方便,將自由面邊界條件寫(xiě)成如下常微分方程的形式

    則t+?t時(shí)刻波面高度η(t+?t)和速度勢(shì)?(t+?t)可寫(xiě)為

    其中

    3 計(jì)算結(jié)果及分析

    下面,我們將利用前面建立的數(shù)值波浪水槽,研究不連續(xù)波對(duì)直墻的作用.在本文中,以d,g和ρ作為無(wú)量綱基準(zhǔn),其他參數(shù)均為無(wú)量綱數(shù).為防止反射波在造波板處發(fā)生二次反射,對(duì)計(jì)算域長(zhǎng)度和模擬時(shí)間進(jìn)行如下控制

    其中,Nw為非周期波浪數(shù),λ為入射波波長(zhǎng).

    在本文的數(shù)值模擬中,設(shè)置自由面上單元的長(zhǎng)度為?x=λ/600,水深方向布置10個(gè)單元,時(shí)間步長(zhǎng)取?t=T/200.該模型的收斂性和精度等已在此前的工作[28]中得到驗(yàn)證,故此處不再贅述.

    3.1 單入射波作用的情況

    首先考慮單個(gè)入射波(Nw=1)的情況.為與文獻(xiàn)[19]中的SGN模型結(jié)果進(jìn)行對(duì)比,設(shè)定入射波的頻率為ω0=0.145,幅值為a=0.05.圖2給出了直墻處自由面爬升的時(shí)間歷程.本文數(shù)值模擬結(jié)果與文獻(xiàn)[19]中的SGN模型的結(jié)果非常吻合,驗(yàn)證了本模型計(jì)算結(jié)果的可靠性.圖3顯示了自由面沿直墻爬高達(dá)到最大值時(shí)(t=1.67T)的自由表面形狀,本文的模擬結(jié)果同樣與SGN模型的結(jié)果相吻合.

    圖4進(jìn)一步測(cè)試了不同入射頻率下單入射波沿直墻的最大波浪爬高R,本模型結(jié)果和SGN模型結(jié)果吻合得依然很好.從圖中可以看出,低頻入射波的最大波浪爬高約是入射波波幅的2倍,接近線性波浪理論的預(yù)測(cè)結(jié)果.而頻率在ω0附近的單入射波可產(chǎn)生達(dá)到2.3倍入射波波高的自由面爬高,表現(xiàn)出了明顯的非線性特征.

    圖2 單入射波在直墻上爬高的時(shí)間歷程Fig.2 History of wave run-up along the vertical wall for the single incident wave

    圖3 單入射波在t=1.67T時(shí)刻的波面形狀Fig.3 Free surface pro fi le at t=1.67T for the single incident wave

    圖4 不同頻率的單入射波沿直墻的最大爬高Fig.4 Maximum wave run-up along the vertical wall for single incident wave of di ff erent frequencies

    3.2 雙入射波作用的情況

    接下來(lái)考慮雙入射波(Nw=2)的情況.入射波頻率設(shè)為ω0=0.02.圖5給出了不同波幅(a=0.048,a=0.049和 a=0.05)的雙入射波撞擊直墻后,自由面沿直墻爬升的時(shí)間歷程.與單入射波的情況相比,雙入射波的兩個(gè)波峰撞擊直墻時(shí),自由面爬高均顯示出高頻振蕩的特征.該高頻振蕩在直墻反射波與第二個(gè)入射波相遇時(shí)尤為明顯.通過(guò)與SGN模型結(jié)果的對(duì)比,可以發(fā)現(xiàn)入射波波幅較小時(shí)(即a=0.048),本模型結(jié)果與SGN模型的結(jié)果有很好的吻合.隨著入射波波幅的增大,本文與SGN模型的結(jié)果,在第一個(gè)入射波爬高、回落和第二個(gè)入射波回落階段,比較吻合;但在第二個(gè)入射波的爬高階段,二者的差別隨入射波波高的增大而增大.對(duì)于幅值為a=0.05的入射波,本文模型得到的最大波浪爬高可達(dá)到4.4倍入射波波高,而SGN模型預(yù)報(bào)的最大爬高為入射波波高的3.8倍.同時(shí),本文波浪爬高歷程的震蕩也更加劇烈.

    圖5 不同波幅下雙入射波直墻上爬高的時(shí)間歷程Fig.5 History of wave run-up along the vertical wall for double incident waves with di ff erent amplitudes

    本文進(jìn)一步考慮了其他頻率的情況.圖6給出了不同入射波頻率下雙入射波沿直墻的最大爬高.從圖中可以看出,隨著入射波波幅的增大,本文模型得到的最大波浪爬高越來(lái)越高于SGN模型的預(yù)報(bào)值.對(duì)于幅值為a=0.05的入射波,兩個(gè)模型結(jié)果的差別隨頻率的增大而逐漸減小,并最終趨于穩(wěn)定.本文的非線性數(shù)值結(jié)果顯示,波浪沿直壁的最大爬高可達(dá)到5倍入射波幅,大于SGN模型預(yù)報(bào)的3.8倍入射波幅.非線性結(jié)果中激發(fā)最大波浪爬高的入射波的頻率,則略小于SGN結(jié)果中最大爬高對(duì)應(yīng)的入射波頻率.在圖6中可以發(fā)現(xiàn),入射波頻率越小,兩種模型的差別越大,說(shuō)明波浪頻率越低,產(chǎn)生的非線性越強(qiáng).

    圖6 不同波幅下不同頻率的雙入射波沿直墻的最大爬高Fig.6 Maximum wave run-up along the vertical wall for double incident waves of di ff erent frequencies and amplitudes

    由于SGN模型是基于淺水假設(shè)且僅滿足弱色散關(guān)系,而本文所采用的完全非線性模型并不受SGN模型假設(shè)的限制,因此有理由相信本文模型能更準(zhǔn)確地反映本文所涉及的非線性問(wèn)題.進(jìn)一步聯(lián)系到3.1節(jié)單入射波的情況,可以得知SGN模型可對(duì)單個(gè)非線性入射波進(jìn)行準(zhǔn)確模擬,而對(duì)于雙入射波情況下的非線性波與波相互作用,SGN模型的模擬精度是有限的,其結(jié)果低估了可能的最大波浪爬高.

    3.3 雙入射波作用的非線性分析

    下面將基于本文模型,對(duì)幅值為a=0.05的雙入射波與直墻相互作用的非線性特性做更進(jìn)一步研究.圖7給出了不同頻率(ω=0.02,ω=0.1和ω=0.2)的雙入射波撞擊直墻后,自由面沿直墻爬升的時(shí)間歷程.其中,線性結(jié)果可通過(guò)將非線性自由面條件(方程(2)和方程(3))線性化得到的.圖中各非線性結(jié)果明顯不同于相應(yīng)的線性結(jié)果.一方面,波浪的非線性成分導(dǎo)致了自由面爬高的異常增大,使得線性預(yù)報(bào)明顯低估了波浪的最大爬高.另一方面,在第二個(gè)波浪的爬升階段,非線性結(jié)果顯示出大幅的高頻振蕩,該振蕩的頻率隨著入射波頻率的增加而減小.這反映出,隨著入射波頻率的增大,第二個(gè)波與前一個(gè)波的反射波發(fā)生相互作用的非線性效應(yīng)在減弱,同時(shí)也導(dǎo)致其最大爬高減小.

    圖7 不同頻率的雙入射波在直墻上爬高的時(shí)間歷程Fig.7 History of wave run-up along the vertical wall for double incident waves of di ff erent wave frequencies

    圖8 雙入射波波形的時(shí)間歷程Fig.8 Evolution of wave pro fi les for double incident waves

    圖8給出了數(shù)值水槽內(nèi)頻率為ω=0.02的雙入射波的波形隨時(shí)間的變化.從圖中可以看出,波浪與直墻的作用過(guò)程可分為撞擊前、撞擊中、撞擊后三個(gè)階段.第二個(gè)波與第一個(gè)波的反射波相遇時(shí),匯合成一個(gè)波峰,其幅值在極短的時(shí)間內(nèi)迅速增大至約兩倍于入射波波高.隨后兩波分離,繼續(xù)沿各自來(lái)向傳播.在該過(guò)程中,兩波之間的水面形成了一系列的振蕩波.隨后,與第一個(gè)波的反射波作用之后的第二個(gè)入射波碰撞直墻,開(kāi)始了爬高過(guò)程,其波高的高頻振蕩也在直墻上檢測(cè)到,此時(shí)在直墻上波浪的最大爬高遠(yuǎn)大于入射波波幅的二倍.

    圖9進(jìn)一步給出了各頻率的雙入射波作用于直墻的波浪力的時(shí)間歷程.其展示出來(lái)的非線性特性與波浪爬高過(guò)程基本相同.

    圖9 不同頻率的雙入射波對(duì)直墻作用力的時(shí)間歷程Fig.9 History of wave force acting on the vertical wall for double incident waves of di ff erent wave frequencies

    圖10給出了不同頻率的雙入射波,沿直墻的最大爬高和對(duì)直墻的最大波浪力,其中設(shè)ω0=0.02.從圖中可以看出,最大波浪爬高和最大波浪力隨入射波頻率的變化趨勢(shì)基本相同,其峰值都出現(xiàn)在略低于ω0的頻率處.峰值處最大波浪爬高可高達(dá)約為5倍的入射波高,而最大波浪力可達(dá)到2.5倍線性預(yù)報(bào)值.對(duì)于ω>ω0,非線性波浪爬高和波浪力的最大值逐漸減小,但始終大于線性預(yù)報(bào)值.

    圖10 不同頻率的雙入射波沿直墻的最大爬高和對(duì)直墻的最大波浪力Fig.10 Maximum wave run-up along the vertical wall and maximum wave force acting on the vertical wall for double incident waves of di ff erent wave frequencies

    圖11進(jìn)一步對(duì)圖 7和圖9中ω0=0.02時(shí)的波浪爬升和波浪力時(shí)間歷程進(jìn)行了頻譜分析.頻率為ω=ω0的組分在線性和非線性結(jié)果中均占主導(dǎo)地位.而相較于線性結(jié)果,非線性結(jié)果在ω=2ω0,ω=3ω0,ω=4ω0等入射頻率整數(shù)倍的高頻處也均出現(xiàn)不可忽略的峰值.這反映了,在入射波與反射波相互作用的過(guò)程中,由于波浪的非線性特性,不同頻率的波浪發(fā)生了能量交換,入射主頻波的能量傳遞給了更高頻的波浪成分.該現(xiàn)象表明雙入射波沿直墻爬高問(wèn)題有典型的非線性特性.

    圖11 雙入射波沿直墻爬高和對(duì)直墻作用力的頻譜Fig.11 Spectra of wave run-up and force on vertical wall for double incident waves

    4 結(jié)論

    本文基于時(shí)域高階邊界元方法,建立了完全非線性二維數(shù)值波浪水槽模型,對(duì)非周期波浪與直墻的相互作用問(wèn)題進(jìn)行了模擬和研究.

    本文首先將全非線性的結(jié)果與SGN模型的結(jié)果進(jìn)行了對(duì)比分析.對(duì)于單入射波撞擊直墻的情況,全非線性結(jié)果與SGN模型的結(jié)果非常吻合.而對(duì)于非線性效應(yīng)更高的大幅值雙入射波問(wèn)題,僅滿足弱色散關(guān)系的SGN模型的模擬結(jié)果大大低估了最大波浪爬高.尤其在第二個(gè)入射波撞擊直墻時(shí),全非線性方法模擬的波面爬高顯示出了更劇烈的振蕩特征.

    本文更進(jìn)一步研究了大幅值雙入射波與直墻的作用問(wèn)題,確認(rèn)了該問(wèn)題的強(qiáng)非線性特性.發(fā)現(xiàn)線性預(yù)報(bào)對(duì)波浪最大爬高有較大低估,而波浪的非線性成分不只導(dǎo)致了自由面爬高的異常增大,也引起了局部自由面的高頻振蕩,而該振蕩的頻率隨著入射波頻率的增加而減小.該物理過(guò)程中,雙入射波作用于直墻的波浪載荷,也展示出了與波浪爬高相似的非線性特性.

    最后,本文對(duì)波浪爬升和波浪力的時(shí)間歷程進(jìn)行了頻譜分析.結(jié)果顯示,由于入射波與反射波的非線性作用,入射主頻波的部分能量傳遞給了更高頻的波浪成分,體現(xiàn)出了典型的非線性特性.

    1張新曙,胡曉峰,尤云祥等.深海多立柱浮式平臺(tái)渦激運(yùn)動(dòng)特性研究.力學(xué)學(xué)報(bào),2016,48(3):593-598(Zhang Xinshu,Hu Xiaofeng,You Yunxiang,et al.Investigation on the characteristics of vortex induced motion of a deep sea muti-column fl oating platform.Chinese Journal of Theoretical and Applied Mechanics,2016,48(3):593-598(in Chinese))

    2趙海峰,石俊,王斌斌.海上浮油定位浮標(biāo)的力學(xué)分析與尺寸優(yōu)化.力學(xué)學(xué)報(bào),2016,48(1):235-242(Zhao Haifeng,Shi Jun,Wang Binbin.Mechanical analysis and sizing optimization of o ff shore oil slick tracking buoy.Chinese Journal of Theoretical and Applied Mechanics,2016,48(1):235-242(in Chinese))

    3徐道臨,盧超,張海成.海上浮動(dòng)機(jī)場(chǎng)動(dòng)力學(xué)建模及非線性動(dòng)力響應(yīng)特性.力學(xué)學(xué)報(bào),2015,47(2):289-300(Xu Daolin,Lu Chao,Zhang Haicheng.Dynamic modeling and nonlinear characteristics of fl oating airport.Chinese Journal of Theoretical and Applied Mechanics,2015,47(2):289-300(in Chinese))

    4 Isaacson M,Cheung KF.Second order wave di ff raction around twodimensional bodies by time-domain method.Applied Ocean Research,1991,13(4):175-186

    5 Grilli ST,Guyenne P,Dias F.A fully non-linear model for threedimensional overturning waves over an arbitrary bottom.International Journal for Numerical Methods in Fluids,2001,35(7):829-867

    6 Jamois E,Fuhrman DR,Bingham HB,et al.A numerical study of nonlinear wave run-up on a vertical plate.Coastal Engineering,2006,53(11):929-945

    7 FuhrmanDR,MadsenPA.Simulationofnonlinearwaverun-upwith a high-order Boussinesq model.Coastal Engineering,2008,55(2):139-154

    8 Pelinovsky E,Kharif C,Talipova T.Large-amplitude long wave interaction with a vertical wall.European Journal of Mechanics-B/Fluids,2008,27(4):409-418

    9孫英偉.直立堤前的破碎波流場(chǎng)特征及波浪力研究.[博士論文].大連:大連理工大學(xué),2012(Sun Yingwei.A study of breaking wave fl ow fi eld in front of the vertical breakwater and wave force.[PhD Thesis].Dalian:Dalian University of Technology,2012(in Chinese))

    10 Chatjigeorgiou IK,Molin B.Third-order interactions,wave run-up and hydrodynamic loading on a vertical plate in an in fi nite wave fi eld.Applied Ocean Research,2013,41:57-64

    11 SuCH,MirieRM.Onhead-oncollisionbetweentwosolitarywaves.Journal of Fluid Mechanics,1980,98:509-525

    12 Fenton JD,Rienecker MM.A Fourier method for solving nonlinear water-wave problems:Application to solitary-wave interactions.Journal of Fluid Mechanics,1982,118:411-443

    13 Kim SK,Liu PLF,Liggett JA.Boundary integral equation solutions for solitary wave generation,propagation and run-up.Coastal Engineering,1983,7(4):299-317

    14 Power H,Chwang AT.On re fl ection of a planar solitary wave at a vertical wall.Wave Motion,1984,6(2):183-195

    15 Cooker MJ,Weidman PD,Bale DS.Re fl ection of a high-amplitude solitary wave at a vertical wall.Journal of Fluid Mechanics,1997,342(1):141-158

    16 Maiti S,Sen D.Computation of solitary waves during propagation and runup on a slope.Ocean Engineering,1999,26(11):1063-1083

    17 Liu PLF,Al-Banaa K.Solitary wave runup and force on a vertical barrier.Journal of Fluid Mechanics,2004,505:225-233

    18 Jian W,Sim SY,Huang Z,et al.Modelling of solitary wave runup on an onshore coastal cli ffby smoothed particle hydrodynamics method.Procedia Engineering,2015,116:88-96

    19 Park H,Cox DT.Empirical wave run-up formula for wave,storm surge and berm width.Coastal Engineering,2015,115:67-78

    20 Chambarel J,Kharif C,Touboul J.Head-on collision of two solitary waves and residual falling jet formation.Nonlinear Processes in Geophysics,2009,16:111-122

    21宣瑞韜,吳衛(wèi),劉樺等.雙孤立波直墻爬高的實(shí)驗(yàn)研究.水動(dòng)力學(xué)研究與進(jìn)展A輯,2013,03:241-251(Xuan Ruitao,Wu Wei,Liu Ye,et al.An experimental study on run-up of double solitary waves against a vertical wall.Chinese Journal of Hydrodynamics A,2013,03:241-251(in Chinese))

    22 Carbone F,Dutykh D,Dudley JM,et al.Extreme wave runup on a vertical cli ff.Geophysical Research Letters,2013,40(12):3138-3143

    23 Green AE,Laws N,Naghdi PM.On the theory of water waves//Proceedings of the Royal Society of London A:Mathematical,Physical and Engineering Sciences.The Royal Society,1974,338:43-55

    24 Green AE,Naghdi PM.A derivation of equations for wave propagation in water of variable depth.Journal of Fluid Mechanics,1976,78:237-246

    25 Zhao BB,Ertekin RC,Duan WY.A comparative study of di ff raction of shallow-water waves by high-level IGN and GN Equations.Journal of Computational Physics,2015,283:129-147

    26 Zhao BB,Duan WY,Ertekin RC.Application of higher-level GN theory to some wave transformation problems.Coastal Engineering,2014,83:177-189

    27 Ning D,Teng B,Zhao H,et al.A comparison of two methods for calculating solid angle coefficients in a BIEM numerical wave tank.Engineering Analysis with Boundary Elements,2010,34(1):92-96

    28 Ning D.Zhuo X,Teng B.Focused waves interaction with a vertical wall//Twenty-third(2013)International O ff shore and Polar Engineering,Alaska,USA,June 30–July 5,2013

    NONLINEAR NUMERICAL STUDY OF NON-PERIODIC WAVES ACTING ON A VERTICAL CLIFF1)

    Li Xiang Zhang Chongwei Ning Dezhi2)Su Peng
    (State Key Laboratory of Coastal and O ff shore Engineering,Dalian University of Technology,Dalian 116024,China)

    In this study,a 2D fully-nonlinear numerical wave tank is developed based on the time-domain higher-order boundary element method.Non-periodic waves acting on a vertical cli ffare investigated.The fully nonlinear kinematic and dynamic boundary conditions are satis fi ed on the instantaneous free surface.The mixed Eulerian-Lagrangian method is adopted to track the transient water particle on the free surface and the fourth order Runge-Kutta method is used to predict the velocity potential and wave elevation on the free surface.Then the acceleration potential technique is adopted to calculate the temporal derivative of the potential on the vertical wall surface,and transient wave loads are obtained by integrating the Bernoulli equation along the wetted wall surface.The obtained nonlinear results are fi rstly compared with solutions of the Serre-Green-Naghdi(SGN)theory.It is observed that,for the highly nonlinear case of double-incidentwaves,the SGN model which only satis fi es the weak dispersion relationship greatly underestimates the maximum wave run-up(MWR).Then,the nonlinear interaction between double-incident-waves and a vertical cli ffis further studied.It is found that the linear prediction also underestimates the MWR.The nonlinearity not only leads to an evident increase of the MWR,but also results in a high-frequency oscillation of the free surface.During this process,nonlinear properties of wave loads are similar to those of the wave run-up.Finally,spectral analysis is performed on histories of wave run-up and wave loads.The dominant frequency wave component is found to transfer its energy to higher frequency components,as a typical nonlinear wave-wave interaction phenomenon.

    nonlinear wave,wave run-up,time-domain simulation,higher-order BEM,GN theory

    O353.2

    A

    10.6052/0459-1879-16-337

    2016–11–20收稿,2017–07–21 錄用,2017–07–21 網(wǎng)絡(luò)版發(fā)表.

    1)國(guó)家自然科學(xué)基金(51679036,51490672),水文水資源與水利工程科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室開(kāi)放基金(2016490111)和中央高校基本科研業(yè)務(wù)費(fèi)專項(xiàng)資金(DUT17GJ202,DUT16RC(3)113)資助項(xiàng)目.

    2)寧德志,教授,主要研究方向:海岸和海洋工程中水動(dòng)力學(xué).E-mail:dzning@dlut.edu.cn

    李翔,張崇偉,寧德志,蘇朋.非周期波浪與直墻作用的非線性數(shù)值研究.力學(xué)學(xué)報(bào),2017,49(5):1042-1049

    Li Xiang,Zhang Chongwei,Ning Dezhi,Su Peng.Nonlinear numerical study of non-periodic waves acting on a vertical cli ff.Chinese Journal of Theoretical and Applied Mechanics,2017,49(5):1042-1049

    猜你喜歡
    直墻入射波波浪
    “直墻平底”之于碑刻·篆刻·書(shū)法的藝術(shù)價(jià)值
    波浪谷和波浪巖
    SHPB入射波相似律與整形技術(shù)的試驗(yàn)與數(shù)值研究
    直墻半圓拱巷道圍巖應(yīng)力分布解析
    深埋直墻拱形隧道穩(wěn)定性研究
    波浪谷隨想
    去看神奇波浪谷
    直墻半圓拱可縮性U型鋼支架卡纜臨界約束力理論研究
    瞬態(tài)激勵(lì)狀態(tài)下樁身速度以及樁身內(nèi)力計(jì)算
    對(duì)機(jī)械波半波損失現(xiàn)象的物理解釋
    電子科技(2015年11期)2015-03-06 01:32:24
    七月丁香在线播放| 亚洲四区av| 免费看av在线观看网站| 男女无遮挡免费网站观看| 宅男免费午夜| 自线自在国产av| 欧美精品一区二区免费开放| 精品国产乱码久久久久久小说| 亚洲激情五月婷婷啪啪| 国产免费福利视频在线观看| www.av在线官网国产| 日韩成人伦理影院| 免费看不卡的av| 欧美国产精品va在线观看不卡| 97在线视频观看| av免费在线看不卡| 1024视频免费在线观看| 亚洲丝袜综合中文字幕| 午夜免费鲁丝| 亚洲综合色网址| 一级毛片黄色毛片免费观看视频| 在线观看一区二区三区激情| videossex国产| 我的女老师完整版在线观看| 久久久久久久久久人人人人人人| 欧美+日韩+精品| 亚洲av在线观看美女高潮| 免费看av在线观看网站| 一二三四中文在线观看免费高清| 亚洲精品色激情综合| 香蕉国产在线看| 久久久久久久久久久久大奶| 国产av码专区亚洲av| 欧美精品高潮呻吟av久久| 天堂8中文在线网| 亚洲美女视频黄频| 精品熟女少妇av免费看| 国产成人精品久久久久久| 插逼视频在线观看| 亚洲av日韩在线播放| 国产极品天堂在线| 日韩不卡一区二区三区视频在线| 大码成人一级视频| 久久青草综合色| 国产在视频线精品| 亚洲,一卡二卡三卡| 搡老乐熟女国产| 日本与韩国留学比较| 成年人免费黄色播放视频| 中文字幕另类日韩欧美亚洲嫩草| 男女边摸边吃奶| 精品一区二区免费观看| 国产成人免费观看mmmm| 综合色丁香网| a 毛片基地| 97在线人人人人妻| 麻豆乱淫一区二区| 在现免费观看毛片| 人妻少妇偷人精品九色| 婷婷色综合www| 亚洲五月色婷婷综合| 丝袜在线中文字幕| 校园人妻丝袜中文字幕| 伊人亚洲综合成人网| 午夜福利在线观看免费完整高清在| 午夜免费观看性视频| 超碰97精品在线观看| 国产爽快片一区二区三区| 成年女人在线观看亚洲视频| 午夜免费男女啪啪视频观看| av有码第一页| 精品酒店卫生间| 18禁观看日本| 亚洲欧洲日产国产| 亚洲精品aⅴ在线观看| av在线观看视频网站免费| 欧美精品亚洲一区二区| 少妇人妻精品综合一区二区| 色视频在线一区二区三区| 99久久中文字幕三级久久日本| 日韩大片免费观看网站| 男女边吃奶边做爰视频| www.熟女人妻精品国产 | 波野结衣二区三区在线| 热re99久久精品国产66热6| 一级毛片电影观看| 免费人妻精品一区二区三区视频| 人人妻人人爽人人添夜夜欢视频| 亚洲图色成人| 国产精品久久久久久av不卡| 亚洲av中文av极速乱| 久久久亚洲精品成人影院| 免费人成在线观看视频色| 欧美日本中文国产一区发布| 91午夜精品亚洲一区二区三区| 日韩av在线免费看完整版不卡| 成年美女黄网站色视频大全免费| 久久99蜜桃精品久久| 亚洲av.av天堂| 青春草亚洲视频在线观看| av电影中文网址| 另类精品久久| 巨乳人妻的诱惑在线观看| 国产探花极品一区二区| 边亲边吃奶的免费视频| 国产一区二区在线观看av| 激情五月婷婷亚洲| 性色avwww在线观看| 久久久久国产精品人妻一区二区| 伦精品一区二区三区| 一区在线观看完整版| 亚洲成人手机| 天天影视国产精品| 久久免费观看电影| 啦啦啦中文免费视频观看日本| 18禁观看日本| 男女下面插进去视频免费观看 | 久久av网站| 69精品国产乱码久久久| 嫩草影院入口| 久久狼人影院| 久久久久精品久久久久真实原创| 亚洲综合精品二区| 熟女av电影| 观看av在线不卡| 丝袜美足系列| 两个人免费观看高清视频| 一区二区三区精品91| 看免费av毛片| 国产色婷婷99| 亚洲熟女精品中文字幕| 国产精品熟女久久久久浪| 日韩一区二区视频免费看| 成人手机av| 一二三四中文在线观看免费高清| 免费少妇av软件| 视频在线观看一区二区三区| av国产精品久久久久影院| 乱人伦中国视频| 寂寞人妻少妇视频99o| 夫妻午夜视频| 女的被弄到高潮叫床怎么办| 人体艺术视频欧美日本| 久久午夜综合久久蜜桃| 一个人免费看片子| 两性夫妻黄色片 | 亚洲国产精品999| 观看av在线不卡| 一区二区三区乱码不卡18| 精品一区二区三卡| 国产精品一区二区在线观看99| 精品福利永久在线观看| 蜜臀久久99精品久久宅男| 国产一区二区在线观看av| 狂野欧美激情性xxxx在线观看| 好男人视频免费观看在线| 最近的中文字幕免费完整| 18禁动态无遮挡网站| 久久久精品免费免费高清| 天堂中文最新版在线下载| 午夜日本视频在线| 久久影院123| 国产黄色视频一区二区在线观看| 极品少妇高潮喷水抽搐| 亚洲在久久综合| 国产永久视频网站| xxx大片免费视频| 一级毛片我不卡| 免费看av在线观看网站| 国产精品人妻久久久久久| 成年美女黄网站色视频大全免费| 又粗又硬又长又爽又黄的视频| 啦啦啦视频在线资源免费观看| 婷婷色av中文字幕| 赤兔流量卡办理| 欧美激情极品国产一区二区三区 | 精品少妇久久久久久888优播| 国产无遮挡羞羞视频在线观看| 午夜免费男女啪啪视频观看| 一级片免费观看大全| www.色视频.com| av电影中文网址| 新久久久久国产一级毛片| 九草在线视频观看| 免费观看a级毛片全部| 中国国产av一级| 国产亚洲一区二区精品| 日韩制服丝袜自拍偷拍| 欧美老熟妇乱子伦牲交| 男女免费视频国产| 亚洲av.av天堂| 9191精品国产免费久久| 三上悠亚av全集在线观看| 久久韩国三级中文字幕| 久久国内精品自在自线图片| 麻豆精品久久久久久蜜桃| 国产高清三级在线| 最近中文字幕2019免费版| 99香蕉大伊视频| 天美传媒精品一区二区| 建设人人有责人人尽责人人享有的| 热99久久久久精品小说推荐| 久久久久精品人妻al黑| 满18在线观看网站| 最黄视频免费看| 777米奇影视久久| 满18在线观看网站| 免费高清在线观看视频在线观看| www.熟女人妻精品国产 | 久久精品人人爽人人爽视色| 视频区图区小说| 国产在线视频一区二区| 波野结衣二区三区在线| 日本午夜av视频| 免费久久久久久久精品成人欧美视频 | 成人国产av品久久久| 亚洲精品国产av成人精品| 国产欧美日韩一区二区三区在线| 中文字幕另类日韩欧美亚洲嫩草| 在线观看三级黄色| 99久国产av精品国产电影| 一区二区三区四区激情视频| 大陆偷拍与自拍| 中文字幕另类日韩欧美亚洲嫩草| 欧美人与性动交α欧美精品济南到 | 国产精品麻豆人妻色哟哟久久| 国产成人午夜福利电影在线观看| 满18在线观看网站| 五月伊人婷婷丁香| 色网站视频免费| www.色视频.com| 久久99热这里只频精品6学生| 国产精品99久久99久久久不卡 | 美女视频免费永久观看网站| 国产午夜精品一二区理论片| 国产黄色免费在线视频| 成人漫画全彩无遮挡| 青春草国产在线视频| 亚洲欧美成人精品一区二区| 两性夫妻黄色片 | 久久人人97超碰香蕉20202| 性高湖久久久久久久久免费观看| 久久精品久久久久久噜噜老黄| 99九九在线精品视频| 丝袜喷水一区| 国产成人精品无人区| 亚洲内射少妇av| 久久久久久人人人人人| av在线观看视频网站免费| 日本欧美国产在线视频| 免费高清在线观看日韩| 亚洲欧美日韩卡通动漫| 少妇被粗大的猛进出69影院 | 丝袜在线中文字幕| 久久99精品国语久久久| 欧美日韩av久久| 国产免费现黄频在线看| 啦啦啦视频在线资源免费观看| 少妇被粗大的猛进出69影院 | 自线自在国产av| 日产精品乱码卡一卡2卡三| 热99国产精品久久久久久7| 日韩,欧美,国产一区二区三区| 欧美另类一区| 亚洲三级黄色毛片| 少妇人妻 视频| 欧美成人午夜精品| 亚洲激情五月婷婷啪啪| 亚洲四区av| 亚洲成色77777| 在线观看国产h片| 精品一区二区三卡| 日本vs欧美在线观看视频| 久久久精品区二区三区| 日韩av在线免费看完整版不卡| 亚洲欧美清纯卡通| 欧美最新免费一区二区三区| 亚洲 欧美一区二区三区| 18+在线观看网站| 欧美精品国产亚洲| 巨乳人妻的诱惑在线观看| 精品国产一区二区久久| 久久人人97超碰香蕉20202| 成人手机av| 欧美xxxx性猛交bbbb| 国产av码专区亚洲av| 久久99热这里只频精品6学生| 亚洲人成网站在线观看播放| 免费av中文字幕在线| 在线观看一区二区三区激情| 少妇人妻久久综合中文| 亚洲美女搞黄在线观看| 黄片无遮挡物在线观看| 欧美成人午夜精品| 国产极品天堂在线| 国产精品国产三级国产专区5o| 男的添女的下面高潮视频| 又黄又粗又硬又大视频| 免费观看av网站的网址| av电影中文网址| 久久久国产精品麻豆| 啦啦啦中文免费视频观看日本| 亚洲精华国产精华液的使用体验| 少妇的逼好多水| av黄色大香蕉| av免费在线看不卡| 久久久久久人妻| 观看av在线不卡| 另类精品久久| 男人操女人黄网站| 成人二区视频| 2022亚洲国产成人精品| 97精品久久久久久久久久精品| 肉色欧美久久久久久久蜜桃| 人妻人人澡人人爽人人| 你懂的网址亚洲精品在线观看| 99久久综合免费| 十八禁网站网址无遮挡| 伊人久久国产一区二区| 国产熟女午夜一区二区三区| 女的被弄到高潮叫床怎么办| 26uuu在线亚洲综合色| 99热这里只有是精品在线观看| 日韩人妻精品一区2区三区| 五月天丁香电影| 肉色欧美久久久久久久蜜桃| 午夜老司机福利剧场| 成年人免费黄色播放视频| 人人澡人人妻人| 热99久久久久精品小说推荐| av天堂久久9| 秋霞伦理黄片| 国产精品三级大全| 国产午夜精品一二区理论片| 久久久久人妻精品一区果冻| 国产成人午夜福利电影在线观看| 亚洲av电影在线观看一区二区三区| 另类亚洲欧美激情| 亚洲国产精品一区三区| 亚洲,欧美精品.| 国产精品一区二区在线不卡| 亚洲精品久久成人aⅴ小说| 一本色道久久久久久精品综合| 日日啪夜夜爽| 久久久久久伊人网av| 久久精品国产鲁丝片午夜精品| 精品国产一区二区久久| 亚洲少妇的诱惑av| 我的女老师完整版在线观看| 午夜影院在线不卡| 在线观看三级黄色| 久久精品aⅴ一区二区三区四区 | 一二三四中文在线观看免费高清| 亚洲成色77777| 精品少妇内射三级| 国产白丝娇喘喷水9色精品| 日本av手机在线免费观看| 欧美+日韩+精品| 国产成人av激情在线播放| 宅男免费午夜| 最近中文字幕高清免费大全6| 深夜精品福利| 9色porny在线观看| 亚洲久久久国产精品| 亚洲av电影在线观看一区二区三区| kizo精华| 久久99热这里只频精品6学生| 不卡视频在线观看欧美| 日日摸夜夜添夜夜爱| 精品人妻熟女毛片av久久网站| 久久综合国产亚洲精品| 人人妻人人爽人人添夜夜欢视频| av在线播放精品| 少妇高潮的动态图| 好男人视频免费观看在线| 久久久久久久久久人人人人人人| 欧美人与性动交α欧美软件 | 国产精品久久久久久精品古装| 黄片无遮挡物在线观看| 午夜精品国产一区二区电影| 国产成人91sexporn| 涩涩av久久男人的天堂| 啦啦啦在线观看免费高清www| 国产精品麻豆人妻色哟哟久久| xxx大片免费视频| 国产淫语在线视频| 激情五月婷婷亚洲| 欧美激情国产日韩精品一区| 久久韩国三级中文字幕| 国产精品久久久久久久电影| 人妻人人澡人人爽人人| 永久免费av网站大全| 国产精品99久久99久久久不卡 | 高清黄色对白视频在线免费看| 亚洲精品日本国产第一区| av网站免费在线观看视频| 97人妻天天添夜夜摸| 国产黄色免费在线视频| 飞空精品影院首页| 色婷婷av一区二区三区视频| 亚洲国产毛片av蜜桃av| 免费久久久久久久精品成人欧美视频 | 你懂的网址亚洲精品在线观看| 久久国内精品自在自线图片| 久久久亚洲精品成人影院| 一边摸一边做爽爽视频免费| 香蕉国产在线看| 亚洲国产成人一精品久久久| 国产成人免费观看mmmm| 久久精品熟女亚洲av麻豆精品| 免费久久久久久久精品成人欧美视频 | 22中文网久久字幕| 激情视频va一区二区三区| 五月伊人婷婷丁香| 日韩制服丝袜自拍偷拍| xxx大片免费视频| 久热这里只有精品99| 制服人妻中文乱码| 欧美日本中文国产一区发布| 欧美日韩国产mv在线观看视频| 亚洲第一区二区三区不卡| 女性被躁到高潮视频| 亚洲精品视频女| 99热这里只有是精品在线观看| 国产熟女午夜一区二区三区| 国产日韩欧美亚洲二区| 久久人人爽人人爽人人片va| 另类亚洲欧美激情| 亚洲欧美日韩卡通动漫| 免费观看av网站的网址| 99热6这里只有精品| 久久99蜜桃精品久久| 午夜精品国产一区二区电影| 少妇的丰满在线观看| 精品国产一区二区三区四区第35| 国产片内射在线| 久久精品国产亚洲av涩爱| 高清不卡的av网站| 熟女电影av网| 国产精品一区www在线观看| 欧美日韩成人在线一区二区| 成年av动漫网址| 久久久久视频综合| 考比视频在线观看| 午夜免费鲁丝| 黑丝袜美女国产一区| 成年人午夜在线观看视频| 亚洲精品久久久久久婷婷小说| 欧美成人精品欧美一级黄| 少妇人妻 视频| 在线观看一区二区三区激情| 中国国产av一级| 亚洲中文av在线| 一级毛片我不卡| av网站免费在线观看视频| 天堂俺去俺来也www色官网| 日本午夜av视频| 国产日韩欧美亚洲二区| 久久久精品94久久精品| 大香蕉久久网| 女性被躁到高潮视频| 超色免费av| 久久人人97超碰香蕉20202| 国产一区亚洲一区在线观看| 黄色 视频免费看| 亚洲色图 男人天堂 中文字幕 | 婷婷色综合大香蕉| av在线老鸭窝| 男女高潮啪啪啪动态图| 日韩,欧美,国产一区二区三区| 午夜老司机福利剧场| 美女大奶头黄色视频| 亚洲精品456在线播放app| 99热网站在线观看| 午夜免费男女啪啪视频观看| 少妇猛男粗大的猛烈进出视频| 在线天堂最新版资源| 久久久久精品性色| 香蕉国产在线看| 国产日韩欧美亚洲二区| 伦理电影免费视频| 久久久久久久久久久免费av| 老司机影院成人| 日产精品乱码卡一卡2卡三| 亚洲美女搞黄在线观看| 色哟哟·www| 亚洲久久久国产精品| 丰满少妇做爰视频| 亚洲一区二区三区欧美精品| 日韩制服丝袜自拍偷拍| 爱豆传媒免费全集在线观看| 在线观看人妻少妇| 国产欧美另类精品又又久久亚洲欧美| 久久久a久久爽久久v久久| 亚洲一级一片aⅴ在线观看| 高清毛片免费看| 国产永久视频网站| 久久99精品国语久久久| 亚洲五月色婷婷综合| 一级,二级,三级黄色视频| 亚洲精品日本国产第一区| 你懂的网址亚洲精品在线观看| 午夜福利影视在线免费观看| 国产精品.久久久| 亚洲精品久久午夜乱码| 久久av网站| 中国国产av一级| 亚洲欧美成人精品一区二区| 男女国产视频网站| 国产又色又爽无遮挡免| 欧美少妇被猛烈插入视频| 秋霞伦理黄片| 最近的中文字幕免费完整| 久久久久久久久久久久大奶| 九色亚洲精品在线播放| 男女啪啪激烈高潮av片| 狂野欧美激情性bbbbbb| 午夜日本视频在线| 日韩一区二区三区影片| 日本av免费视频播放| 欧美日韩亚洲高清精品| av免费观看日本| 国产精品久久久久久精品古装| 国产 一区精品| 亚洲av综合色区一区| 国产成人91sexporn| 亚洲人成网站在线观看播放| 9191精品国产免费久久| 2022亚洲国产成人精品| 国产亚洲精品第一综合不卡 | 日韩精品免费视频一区二区三区 | 熟女av电影| 国语对白做爰xxxⅹ性视频网站| 日产精品乱码卡一卡2卡三| 极品少妇高潮喷水抽搐| 久久热在线av| 在线天堂最新版资源| 丰满饥渴人妻一区二区三| 美女国产视频在线观看| 日韩一区二区视频免费看| 久久国产亚洲av麻豆专区| av网站免费在线观看视频| 免费黄网站久久成人精品| 97在线视频观看| 亚洲人成77777在线视频| 狂野欧美激情性bbbbbb| 国产成人91sexporn| 中国美白少妇内射xxxbb| 亚洲国产最新在线播放| 国产精品偷伦视频观看了| av在线app专区| 亚洲情色 制服丝袜| 天堂中文最新版在线下载| 飞空精品影院首页| 男女边吃奶边做爰视频| 国产男女超爽视频在线观看| 韩国高清视频一区二区三区| 黑人高潮一二区| 欧美人与善性xxx| 久久国内精品自在自线图片| 婷婷色综合大香蕉| 成人午夜精彩视频在线观看| 卡戴珊不雅视频在线播放| 成人毛片a级毛片在线播放| 内地一区二区视频在线| 国产精品一区二区在线不卡| 免费观看在线日韩| 另类精品久久| 国产精品免费大片| 国产白丝娇喘喷水9色精品| 一级黄片播放器| av片东京热男人的天堂| 亚洲成av片中文字幕在线观看 | 超碰97精品在线观看| 中文精品一卡2卡3卡4更新| 伦理电影大哥的女人| 免费黄色在线免费观看| 亚洲精品乱久久久久久| 国产欧美亚洲国产| 侵犯人妻中文字幕一二三四区| 我的女老师完整版在线观看| 免费看光身美女| 亚洲人与动物交配视频| 十分钟在线观看高清视频www| 不卡视频在线观看欧美| 亚洲欧美一区二区三区国产| 极品人妻少妇av视频| 高清视频免费观看一区二区| 曰老女人黄片| 久久久久久伊人网av| 亚洲精品美女久久久久99蜜臀 | 99热国产这里只有精品6| 精品人妻在线不人妻| 国产免费一区二区三区四区乱码| 街头女战士在线观看网站| 成年女人在线观看亚洲视频| 亚洲国产最新在线播放| av免费观看日本| 热99久久久久精品小说推荐| 亚洲四区av| 搡女人真爽免费视频火全软件| 国产伦理片在线播放av一区| 美女国产高潮福利片在线看| 一区二区av电影网| 亚洲,欧美,日韩| 九色亚洲精品在线播放| 国产精品国产三级国产av玫瑰| 岛国毛片在线播放| 国产乱来视频区| 五月玫瑰六月丁香| av有码第一页| 久久精品久久久久久久性| 有码 亚洲区| 一级毛片黄色毛片免费观看视频| 中文字幕制服av|