• <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
    国产麻豆69| 亚洲成国产人片在线观看| 黄色 视频免费看| 久久99精品国语久久久| 日韩一区二区三区影片| 久久久久久久久久久久大奶| 日本欧美国产在线视频| 狠狠婷婷综合久久久久久88av| 色播在线永久视频| 国产精品.久久久| 脱女人内裤的视频| 看免费成人av毛片| 啦啦啦在线免费观看视频4| 成年人免费黄色播放视频| 脱女人内裤的视频| 五月开心婷婷网| 午夜福利视频精品| 亚洲人成77777在线视频| 亚洲色图综合在线观看| 一二三四社区在线视频社区8| 黄色怎么调成土黄色| 国产欧美日韩精品亚洲av| 久久 成人 亚洲| 夫妻性生交免费视频一级片| 精品国产国语对白av| 欧美xxⅹ黑人| 下体分泌物呈黄色| 日韩大码丰满熟妇| 久久久久久免费高清国产稀缺| 99久久精品国产亚洲精品| 日韩制服骚丝袜av| av网站免费在线观看视频| 亚洲国产日韩一区二区| 无遮挡黄片免费观看| 999久久久国产精品视频| 欧美日本中文国产一区发布| 日日摸夜夜添夜夜爱| 黄频高清免费视频| 一二三四在线观看免费中文在| 成人手机av| 亚洲图色成人| 亚洲国产欧美在线一区| 婷婷丁香在线五月| 成人影院久久| av在线播放精品| svipshipincom国产片| av一本久久久久| 一级毛片电影观看| 日本wwww免费看| 久久久亚洲精品成人影院| 亚洲中文av在线| 日韩 欧美 亚洲 中文字幕| 高清不卡的av网站| 国产成人精品久久二区二区91| 妹子高潮喷水视频| av片东京热男人的天堂| 两个人免费观看高清视频| 日本av手机在线免费观看| 成人亚洲欧美一区二区av| 日韩 亚洲 欧美在线| 青草久久国产| 丝瓜视频免费看黄片| 亚洲av欧美aⅴ国产| 日本av手机在线免费观看| 香蕉丝袜av| 亚洲欧洲日产国产| 亚洲视频免费观看视频| 久久 成人 亚洲| 亚洲av欧美aⅴ国产| 久久久久久久大尺度免费视频| 久久久久久亚洲精品国产蜜桃av| 久久精品aⅴ一区二区三区四区| xxxhd国产人妻xxx| 51午夜福利影视在线观看| 日韩一区二区三区影片| e午夜精品久久久久久久| 80岁老熟妇乱子伦牲交| 久久国产亚洲av麻豆专区| 满18在线观看网站| 晚上一个人看的免费电影| 日韩中文字幕视频在线看片| 午夜福利,免费看| 精品久久久精品久久久| 国产成人影院久久av| 十八禁高潮呻吟视频| 丰满人妻熟妇乱又伦精品不卡| 成人免费观看视频高清| 老司机在亚洲福利影院| 国产国语露脸激情在线看| 久久精品aⅴ一区二区三区四区| 久久国产亚洲av麻豆专区| 一级,二级,三级黄色视频| 国产爽快片一区二区三区| 美女扒开内裤让男人捅视频| 日韩,欧美,国产一区二区三区| 免费在线观看视频国产中文字幕亚洲 | 婷婷色综合大香蕉| 90打野战视频偷拍视频| 成年人午夜在线观看视频| 欧美激情极品国产一区二区三区| 2021少妇久久久久久久久久久| 久久精品亚洲av国产电影网| 国产成人欧美| 国产精品三级大全| 欧美日韩综合久久久久久| 久久精品久久久久久久性| 九色亚洲精品在线播放| 丝袜脚勾引网站| 我要看黄色一级片免费的| 精品少妇一区二区三区视频日本电影| 欧美性长视频在线观看| 好男人视频免费观看在线| 黄色a级毛片大全视频| 黑人猛操日本美女一级片| 一本色道久久久久久精品综合| 看十八女毛片水多多多| 国产成人欧美| 亚洲情色 制服丝袜| 国产精品 国内视频| 丝袜喷水一区| 啦啦啦中文免费视频观看日本| 免费看十八禁软件| 在线观看一区二区三区激情| 国产1区2区3区精品| 国产精品一国产av| av视频免费观看在线观看| 超碰成人久久| 两人在一起打扑克的视频| 亚洲av日韩精品久久久久久密 | 国产一级毛片在线| 国产人伦9x9x在线观看| 永久免费av网站大全| 丁香六月天网| 成人18禁高潮啪啪吃奶动态图| 亚洲七黄色美女视频| av在线app专区| 91精品国产国语对白视频| 中文欧美无线码| 性色av一级| 欧美精品av麻豆av| 大陆偷拍与自拍| 亚洲欧美精品综合一区二区三区| 女人精品久久久久毛片| 女人久久www免费人成看片| 妹子高潮喷水视频| 两个人看的免费小视频| 男女边摸边吃奶| 青春草视频在线免费观看| 超碰成人久久| 大码成人一级视频| 欧美xxⅹ黑人| 搡老岳熟女国产| 一本综合久久免费| 亚洲av男天堂| 国产有黄有色有爽视频| 十八禁高潮呻吟视频| 咕卡用的链子| 桃花免费在线播放| 国产高清不卡午夜福利| 午夜福利免费观看在线| 国产一区二区三区av在线| 亚洲欧洲日产国产| 97人妻天天添夜夜摸| 脱女人内裤的视频| 丰满少妇做爰视频| 在现免费观看毛片| 久久久久网色| a 毛片基地| 一级a爱视频在线免费观看| 久热爱精品视频在线9| 国产日韩一区二区三区精品不卡| 国产精品成人在线| 国产精品二区激情视频| 丝袜在线中文字幕| 亚洲精品久久午夜乱码| 久久精品久久久久久噜噜老黄| 国产精品一区二区在线观看99| 亚洲av在线观看美女高潮| 一本综合久久免费| 亚洲国产欧美在线一区| 亚洲国产欧美一区二区综合| 十分钟在线观看高清视频www| 好男人视频免费观看在线| 国产成人欧美| 国产一区二区三区综合在线观看| 成年美女黄网站色视频大全免费| 国产熟女午夜一区二区三区| 在线观看国产h片| 少妇 在线观看| 亚洲视频免费观看视频| 电影成人av| 别揉我奶头~嗯~啊~动态视频 | 一边亲一边摸免费视频| 夜夜骑夜夜射夜夜干| 新久久久久国产一级毛片| 亚洲中文av在线| 国产极品粉嫩免费观看在线| av国产精品久久久久影院| 成年av动漫网址| 精品人妻在线不人妻| 叶爱在线成人免费视频播放| 美女主播在线视频| 在线观看国产h片| 在线 av 中文字幕| 久久久久久久大尺度免费视频| 久久精品国产亚洲av涩爱| 国产精品偷伦视频观看了| 在线av久久热| 高清不卡的av网站| 国产精品免费视频内射| 涩涩av久久男人的天堂| 亚洲精品美女久久久久99蜜臀 | 中文精品一卡2卡3卡4更新| 国语对白做爰xxxⅹ性视频网站| 丝袜喷水一区| 中文字幕另类日韩欧美亚洲嫩草| 久久99精品国语久久久| 大码成人一级视频| 不卡av一区二区三区| 国产又爽黄色视频| 国产成人精品久久二区二区91| 亚洲国产av影院在线观看| 国产精品免费大片| 国产片特级美女逼逼视频| videos熟女内射| 视频区欧美日本亚洲| 少妇人妻 视频| 高清不卡的av网站| 国产亚洲av片在线观看秒播厂| 国产成人精品久久久久久| 午夜激情久久久久久久| 日韩中文字幕欧美一区二区 | 国产精品三级大全| 久久久久国产精品人妻一区二区| 另类亚洲欧美激情| 1024视频免费在线观看| 欧美精品av麻豆av| 亚洲精品中文字幕在线视频| 老汉色av国产亚洲站长工具| 美女主播在线视频| 啦啦啦视频在线资源免费观看| 在线亚洲精品国产二区图片欧美| 黄片小视频在线播放| 老司机亚洲免费影院| 少妇精品久久久久久久| 亚洲欧美一区二区三区黑人| 最新在线观看一区二区三区 | 热re99久久精品国产66热6| av欧美777| 国产伦人伦偷精品视频| 国产97色在线日韩免费| 亚洲,一卡二卡三卡| 91字幕亚洲| 久久久国产精品麻豆| 国产熟女午夜一区二区三区| 在线观看免费高清a一片| 免费在线观看影片大全网站 | 欧美成狂野欧美在线观看| 欧美少妇被猛烈插入视频| 水蜜桃什么品种好| 精品亚洲成国产av| 亚洲九九香蕉| 国产一区亚洲一区在线观看| 欧美精品啪啪一区二区三区 | 国产精品一区二区精品视频观看| 人人妻,人人澡人人爽秒播 | 色综合欧美亚洲国产小说| 高清欧美精品videossex| 国产精品麻豆人妻色哟哟久久| 久久中文字幕一级| 久久热在线av| 手机成人av网站| 女性生殖器流出的白浆| 天堂8中文在线网| 国产精品香港三级国产av潘金莲 | 97在线人人人人妻| 国产黄频视频在线观看| 亚洲欧洲精品一区二区精品久久久| 黄色 视频免费看| 一级片免费观看大全| 亚洲精品国产av成人精品| 国产一区二区三区综合在线观看| 青青草视频在线视频观看| 成年动漫av网址| 日本vs欧美在线观看视频| 另类精品久久| 91国产中文字幕| 欧美精品亚洲一区二区| 男女下面插进去视频免费观看| 午夜日韩欧美国产| 女人久久www免费人成看片| 两个人看的免费小视频| 18禁国产床啪视频网站| 欧美精品av麻豆av| 亚洲人成网站在线观看播放| 一级毛片电影观看| av视频免费观看在线观看| 大码成人一级视频| 国产在线观看jvid| 国产麻豆69| 欧美 日韩 精品 国产| 亚洲人成电影观看| 精品少妇一区二区三区视频日本电影| 日韩av不卡免费在线播放| 超碰成人久久| 51午夜福利影视在线观看| 王馨瑶露胸无遮挡在线观看| 亚洲成色77777| 久久青草综合色| 久久精品熟女亚洲av麻豆精品| 欧美黄色淫秽网站| 男女午夜视频在线观看| 日本vs欧美在线观看视频| 另类精品久久| 久久人妻熟女aⅴ| www.熟女人妻精品国产| 自拍欧美九色日韩亚洲蝌蚪91| 午夜激情av网站| 久久国产精品大桥未久av| 久久精品国产亚洲av高清一级| 男女高潮啪啪啪动态图| 好男人电影高清在线观看| 欧美精品人与动牲交sv欧美| 亚洲,一卡二卡三卡| 人妻 亚洲 视频| 国产欧美日韩综合在线一区二区| 亚洲欧美中文字幕日韩二区| 女人久久www免费人成看片| 日本一区二区免费在线视频| 美女国产高潮福利片在线看| 亚洲中文av在线| 2018国产大陆天天弄谢| 国产精品国产三级国产专区5o| 国产成人精品在线电影| 国产精品欧美亚洲77777| 欧美+亚洲+日韩+国产| 亚洲黑人精品在线| 国产av国产精品国产| 极品少妇高潮喷水抽搐| 色精品久久人妻99蜜桃| 少妇猛男粗大的猛烈进出视频| 久久国产精品男人的天堂亚洲| 午夜av观看不卡| 久久人人爽人人片av| 久久久亚洲精品成人影院| 9热在线视频观看99| 欧美xxⅹ黑人| 亚洲国产精品成人久久小说| 伊人亚洲综合成人网| 中文乱码字字幕精品一区二区三区| 欧美日韩亚洲综合一区二区三区_| 男男h啪啪无遮挡| 啦啦啦 在线观看视频| 我的亚洲天堂| 国产精品九九99| 高清不卡的av网站| 午夜福利免费观看在线| 丰满少妇做爰视频| 欧美xxⅹ黑人| 精品福利观看| 十八禁人妻一区二区| 亚洲五月色婷婷综合| 亚洲精品自拍成人| 热re99久久精品国产66热6| 国产精品国产三级专区第一集| 国产亚洲av高清不卡| 亚洲精品自拍成人| 久久精品久久久久久久性| 免费人妻精品一区二区三区视频| 超碰97精品在线观看| 后天国语完整版免费观看| 精品国产超薄肉色丝袜足j| 大香蕉久久成人网| 悠悠久久av| 69精品国产乱码久久久| 亚洲av成人不卡在线观看播放网 | 一级毛片电影观看| 亚洲男人天堂网一区| 久久久久久久久免费视频了| 精品人妻熟女毛片av久久网站| 青青草视频在线视频观看| 国产精品二区激情视频| 国产1区2区3区精品| 亚洲av日韩精品久久久久久密 | 人人妻人人爽人人添夜夜欢视频| av在线播放精品| 五月开心婷婷网| 精品少妇内射三级| 久久中文字幕一级| 国产精品久久久久久精品电影小说| 男女下面插进去视频免费观看| 曰老女人黄片| 亚洲人成电影观看| 女人被躁到高潮嗷嗷叫费观| 亚洲欧美成人综合另类久久久| 一二三四在线观看免费中文在| 99国产综合亚洲精品| 国产在视频线精品| 美女脱内裤让男人舔精品视频| 国产精品一区二区免费欧美 | 青草久久国产| 久久久久久久国产电影| 亚洲激情五月婷婷啪啪| 蜜桃国产av成人99| 久久久精品国产亚洲av高清涩受| 少妇精品久久久久久久| √禁漫天堂资源中文www| 久久天堂一区二区三区四区| 精品少妇久久久久久888优播| 九草在线视频观看| 丰满少妇做爰视频| 人成视频在线观看免费观看| 国产精品一二三区在线看| 欧美精品一区二区免费开放| 久久久久久久国产电影| 日韩大片免费观看网站| 人妻 亚洲 视频| 精品第一国产精品| 大香蕉久久网| 欧美乱码精品一区二区三区| 大话2 男鬼变身卡| 精品少妇一区二区三区视频日本电影| 麻豆国产av国片精品| 伊人久久大香线蕉亚洲五| 精品久久久久久电影网| 国产xxxxx性猛交| 一级毛片电影观看| 精品人妻熟女毛片av久久网站| 在线观看免费高清a一片| 美女脱内裤让男人舔精品视频| 国产精品熟女久久久久浪| a级毛片在线看网站| 日韩中文字幕欧美一区二区 | 国产真人三级小视频在线观看| 多毛熟女@视频| 亚洲国产精品999| av在线播放精品| 欧美日韩综合久久久久久| 女性被躁到高潮视频| 久久天躁狠狠躁夜夜2o2o | 国产97色在线日韩免费| 国产主播在线观看一区二区 | 国产av精品麻豆| 免费日韩欧美在线观看| 91九色精品人成在线观看| 五月开心婷婷网| 最近手机中文字幕大全| 久久天躁狠狠躁夜夜2o2o | 各种免费的搞黄视频| 一区福利在线观看| 日本av手机在线免费观看| 人人妻人人澡人人爽人人夜夜| 高清欧美精品videossex| 久久免费观看电影| 婷婷色av中文字幕| 亚洲伊人久久精品综合| 久久精品aⅴ一区二区三区四区| 国产一区二区 视频在线| 国产精品九九99| 夜夜骑夜夜射夜夜干| 久久免费观看电影| 纯流量卡能插随身wifi吗| 桃花免费在线播放| 久久久欧美国产精品| 看免费av毛片| 黄频高清免费视频| 亚洲三区欧美一区| 久久久精品国产亚洲av高清涩受| 日韩精品免费视频一区二区三区| 人人妻人人添人人爽欧美一区卜| 黄色视频不卡| 超碰成人久久| 成人手机av| 亚洲国产精品国产精品| 午夜久久久在线观看| 久久av网站| 啦啦啦在线免费观看视频4| 久久中文字幕一级| 一边摸一边做爽爽视频免费| 国产欧美日韩综合在线一区二区| 美女午夜性视频免费| 搡老乐熟女国产| 五月开心婷婷网| 免费日韩欧美在线观看| 欧美成狂野欧美在线观看| 亚洲少妇的诱惑av| 丰满人妻熟妇乱又伦精品不卡| 久久鲁丝午夜福利片| 午夜福利影视在线免费观看| 母亲3免费完整高清在线观看| 免费看av在线观看网站| 女人精品久久久久毛片| 欧美老熟妇乱子伦牲交| 午夜91福利影院| 国产精品 国内视频| 久久久久久人人人人人| 色综合欧美亚洲国产小说| 国产成人一区二区在线| 亚洲欧美清纯卡通| 国产成人影院久久av| 韩国精品一区二区三区| 夫妻午夜视频| 久久久国产精品麻豆| 曰老女人黄片| 操出白浆在线播放| 男人添女人高潮全过程视频| 精品少妇久久久久久888优播| 少妇裸体淫交视频免费看高清 | 男女午夜视频在线观看| 中文字幕最新亚洲高清| 黑丝袜美女国产一区| 97精品久久久久久久久久精品| 国产国语露脸激情在线看| 亚洲国产看品久久| 啦啦啦在线免费观看视频4| 最近最新中文字幕大全免费视频 | 国产精品亚洲av一区麻豆| 免费看不卡的av| 精品高清国产在线一区| 久久天躁狠狠躁夜夜2o2o | 男女边吃奶边做爰视频| 成人免费观看视频高清| 国产精品久久久久久精品电影小说| 免费观看a级毛片全部| 少妇精品久久久久久久| 人人妻人人添人人爽欧美一区卜| 亚洲欧美精品自产自拍| 国产精品久久久久久精品古装| 少妇的丰满在线观看| 午夜福利,免费看| 无限看片的www在线观看| 嫁个100分男人电影在线观看 | 搡老乐熟女国产| 国产高清不卡午夜福利| 亚洲中文字幕日韩| 十分钟在线观看高清视频www| 亚洲,欧美精品.| 日韩免费高清中文字幕av| 亚洲人成电影观看| 免费黄频网站在线观看国产| 天堂8中文在线网| 天堂中文最新版在线下载| 性色av一级| 国产精品人妻久久久影院| 91精品三级在线观看| 在现免费观看毛片| 日韩伦理黄色片| 久久久精品免费免费高清| 亚洲av片天天在线观看| 欧美亚洲 丝袜 人妻 在线| 精品亚洲成a人片在线观看| 久久国产精品人妻蜜桃| 欧美xxⅹ黑人| 在线亚洲精品国产二区图片欧美| 多毛熟女@视频| 精品久久久久久久毛片微露脸 | 日本五十路高清| 亚洲国产精品一区三区| 成年美女黄网站色视频大全免费| 欧美 日韩 精品 国产| 国产激情久久老熟女| 日韩熟女老妇一区二区性免费视频| 亚洲专区国产一区二区| 国产成人免费观看mmmm| 少妇被粗大的猛进出69影院| 国语对白做爰xxxⅹ性视频网站| 两性夫妻黄色片| 国产日韩欧美在线精品| 黄色毛片三级朝国网站| 这个男人来自地球电影免费观看| 久久久国产欧美日韩av| 91老司机精品| 国产成人精品久久二区二区91| 男人舔女人的私密视频| 国产高清国产精品国产三级| 国产亚洲精品久久久久5区| 欧美 日韩 精品 国产| 国产精品欧美亚洲77777| 成人亚洲欧美一区二区av| 国产精品免费视频内射| 七月丁香在线播放| 精品人妻熟女毛片av久久网站| 丝袜脚勾引网站| 久久久精品国产亚洲av高清涩受| 黄片小视频在线播放| 在线 av 中文字幕| 国产精品秋霞免费鲁丝片| 亚洲国产精品999| 国产老妇伦熟女老妇高清| 少妇人妻 视频| 亚洲国产精品999| 久久久久网色| av福利片在线| 国产欧美日韩一区二区三 | 国产无遮挡羞羞视频在线观看| 亚洲综合色网址| 亚洲av日韩精品久久久久久密 | 电影成人av| 91成人精品电影| 高清不卡的av网站| 免费高清在线观看视频在线观看| 日本欧美视频一区| 在线观看免费高清a一片| 一区二区三区精品91| 欧美精品高潮呻吟av久久| 亚洲,一卡二卡三卡| 亚洲黑人精品在线| 免费在线观看影片大全网站 | 精品久久蜜臀av无| 久久久久精品人妻al黑| 一级a爱视频在线免费观看| 日韩大片免费观看网站| 国产一区二区激情短视频 | 十八禁网站网址无遮挡| 亚洲欧美一区二区三区国产|