• <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天美| 97超视频在线观看视频| 乱人视频在线观看| 国产熟女欧美一区二区| 亚洲成a人片在线一区二区| 香蕉av资源在线| 久久精品国产亚洲av涩爱 | 女生性感内裤真人,穿戴方法视频| 国产亚洲精品久久久久久毛片| 亚洲欧美激情综合另类| 国产精品日韩av在线免费观看| 日韩在线高清观看一区二区三区 | 亚洲人成网站在线播| АⅤ资源中文在线天堂| 国产女主播在线喷水免费视频网站 | 91精品国产九色| 美女大奶头视频| 一区二区三区激情视频| 国产白丝娇喘喷水9色精品| 窝窝影院91人妻| 亚洲四区av| 99热只有精品国产| 日本成人三级电影网站| 国产综合懂色| 久久天躁狠狠躁夜夜2o2o| 国产成年人精品一区二区| 欧美区成人在线视频| 国产真实乱freesex| 欧美日韩中文字幕国产精品一区二区三区| 亚洲av二区三区四区| 国内精品久久久久久久电影| 亚洲精品色激情综合| 亚洲自拍偷在线| 久久热精品热| 国产一级毛片七仙女欲春2| 一个人看视频在线观看www免费| 久久久久久久久久久丰满 | 99热这里只有精品一区| 亚洲欧美精品综合久久99| 99九九线精品视频在线观看视频| 长腿黑丝高跟| 国产私拍福利视频在线观看| 久久精品夜夜夜夜夜久久蜜豆| 国产精品av视频在线免费观看| 99久久精品热视频| 亚洲av美国av| 亚洲av中文av极速乱 | 少妇人妻一区二区三区视频| 欧美日韩瑟瑟在线播放| 亚洲欧美精品综合久久99| 夜夜爽天天搞| 国产91精品成人一区二区三区| av.在线天堂| 欧美潮喷喷水| 91久久精品国产一区二区成人| 久久久成人免费电影| 波野结衣二区三区在线| 又爽又黄无遮挡网站| 国产精品日韩av在线免费观看| 午夜免费男女啪啪视频观看 | 亚洲无线在线观看| 中出人妻视频一区二区| 热99在线观看视频| 久久久国产成人免费| 亚洲成人中文字幕在线播放| 身体一侧抽搐| 欧美又色又爽又黄视频| 国产免费男女视频| 亚洲成a人片在线一区二区| 日韩,欧美,国产一区二区三区 | 日本一本二区三区精品| 18禁裸乳无遮挡免费网站照片| 五月玫瑰六月丁香| 久久精品91蜜桃| 亚洲精华国产精华液的使用体验 | 午夜福利高清视频| 成人永久免费在线观看视频| 欧美高清性xxxxhd video| 久久久久国内视频| 精品久久国产蜜桃| 精品人妻视频免费看| 一级黄色大片毛片| 啦啦啦观看免费观看视频高清| 欧美色视频一区免费| 搡女人真爽免费视频火全软件 | 国产精品嫩草影院av在线观看 | 国产中年淑女户外野战色| 一级黄片播放器| 男人的好看免费观看在线视频| 少妇被粗大猛烈的视频| 亚洲第一区二区三区不卡| 国产伦精品一区二区三区四那| av专区在线播放| 国产淫片久久久久久久久| 国产私拍福利视频在线观看| 91午夜精品亚洲一区二区三区 | 久久中文看片网| 久久久午夜欧美精品| av专区在线播放| aaaaa片日本免费| 久久天躁狠狠躁夜夜2o2o| 亚洲国产欧美人成| 搡老岳熟女国产| 国内少妇人妻偷人精品xxx网站| 欧美日韩国产亚洲二区| 真人做人爱边吃奶动态| 亚洲av中文字字幕乱码综合| 免费一级毛片在线播放高清视频| 国内精品久久久久精免费| bbb黄色大片| 日韩强制内射视频| 亚洲内射少妇av| 搡老岳熟女国产| 亚洲av二区三区四区| 赤兔流量卡办理| 欧美另类亚洲清纯唯美| 久99久视频精品免费| 日日干狠狠操夜夜爽| 国产成年人精品一区二区| 日韩欧美精品v在线| 欧美区成人在线视频| 国产亚洲精品综合一区在线观看| 成年免费大片在线观看| 日韩欧美精品免费久久| 亚洲最大成人av| 搡女人真爽免费视频火全软件 | 1024手机看黄色片| av在线天堂中文字幕| 久久精品国产亚洲av香蕉五月| 成人毛片a级毛片在线播放| 不卡视频在线观看欧美| 精品一区二区三区视频在线| 中亚洲国语对白在线视频| 一区福利在线观看| 久久精品91蜜桃| 中文字幕高清在线视频| 国产aⅴ精品一区二区三区波| 日韩欧美在线二视频| 一级黄片播放器| 亚洲av第一区精品v没综合| 国产精品久久久久久精品电影| 亚洲精品一区av在线观看| 免费av不卡在线播放| 亚洲av中文字字幕乱码综合| 亚洲内射少妇av| 亚洲av一区综合| a级毛片免费高清观看在线播放| 午夜福利18| 成人三级黄色视频| avwww免费| 国产精品三级大全| 色精品久久人妻99蜜桃| 中国美女看黄片| 少妇猛男粗大的猛烈进出视频 | 哪里可以看免费的av片| 午夜福利欧美成人| 99视频精品全部免费 在线| 99热6这里只有精品| 亚洲专区中文字幕在线| 麻豆成人av在线观看| 亚洲欧美日韩东京热| 国产激情偷乱视频一区二区| 亚洲图色成人| 欧美中文日本在线观看视频| 亚洲第一电影网av| 久久久国产成人免费| 日韩欧美精品免费久久| 欧美成人性av电影在线观看| 99国产精品一区二区蜜桃av| 国产精品三级大全| 丰满乱子伦码专区| 麻豆国产97在线/欧美| 噜噜噜噜噜久久久久久91| 男女之事视频高清在线观看| 啪啪无遮挡十八禁网站| 亚洲无线在线观看| 亚洲av成人av| 日韩强制内射视频| 亚洲狠狠婷婷综合久久图片| 91久久精品国产一区二区成人| 国产一区二区三区在线臀色熟女| 午夜福利欧美成人| 午夜福利18| 国产精品久久久久久亚洲av鲁大| 国产精品久久久久久久电影| 黄色配什么色好看| 亚洲国产日韩欧美精品在线观看| 久久精品国产亚洲av天美| 久久久久性生活片| 亚洲av免费高清在线观看| 天美传媒精品一区二区| 久久久久久大精品| 国产国拍精品亚洲av在线观看| 他把我摸到了高潮在线观看| 少妇丰满av| 黄片wwwwww| 村上凉子中文字幕在线| 嫩草影视91久久| 一进一出抽搐gif免费好疼| 国产午夜福利久久久久久| 亚洲中文字幕一区二区三区有码在线看| 国产亚洲精品av在线| 黄色一级大片看看| 亚洲成人久久性| 日本色播在线视频| 色综合亚洲欧美另类图片| 有码 亚洲区| 成人毛片a级毛片在线播放| 又粗又爽又猛毛片免费看| 欧美激情久久久久久爽电影| 又紧又爽又黄一区二区| 国内精品美女久久久久久| 波野结衣二区三区在线| 看十八女毛片水多多多| 日韩国内少妇激情av| 亚洲精华国产精华液的使用体验 | 亚洲经典国产精华液单| 久久久精品大字幕| 成年女人毛片免费观看观看9| 日本爱情动作片www.在线观看 | 日本精品一区二区三区蜜桃| 99久久中文字幕三级久久日本| 久久亚洲精品不卡| 男女下面进入的视频免费午夜| 午夜福利欧美成人| 久久精品国产自在天天线| 亚洲欧美日韩无卡精品| 精品人妻视频免费看| 美女大奶头视频| 精品久久久久久,| 男人舔女人下体高潮全视频| 一进一出抽搐gif免费好疼| 亚洲成人中文字幕在线播放| АⅤ资源中文在线天堂| 日本a在线网址| 51国产日韩欧美| 日韩欧美国产在线观看| av在线亚洲专区| 日韩一本色道免费dvd| 日本免费一区二区三区高清不卡| 亚洲成人精品中文字幕电影| 日韩欧美国产一区二区入口| 97人妻精品一区二区三区麻豆| 亚洲一级一片aⅴ在线观看| 国内毛片毛片毛片毛片毛片| 我要看日韩黄色一级片| 我要看日韩黄色一级片| 亚洲av一区综合| 欧美zozozo另类| 欧美在线一区亚洲| 一边摸一边抽搐一进一小说| 亚洲欧美清纯卡通| 麻豆国产av国片精品| 久久6这里有精品| 啪啪无遮挡十八禁网站| 乱系列少妇在线播放| 亚洲avbb在线观看| 国产欧美日韩精品一区二区| 一级黄色大片毛片| 91狼人影院| 天堂影院成人在线观看| 午夜激情福利司机影院| 国产中年淑女户外野战色| 日韩欧美在线二视频| 真人做人爱边吃奶动态| 一个人看的www免费观看视频| 亚洲av中文av极速乱 | 嫩草影院入口| 午夜精品久久久久久毛片777| 极品教师在线视频| 亚洲久久久久久中文字幕| 欧美潮喷喷水| 欧美精品啪啪一区二区三区| 国产精品一区二区性色av| eeuss影院久久| 日韩欧美一区二区三区在线观看| 国产高清激情床上av| 男女下面进入的视频免费午夜| 黄色女人牲交| 欧美黑人欧美精品刺激| 国产精品日韩av在线免费观看| 免费大片18禁| 日韩av在线大香蕉| 久久久久国内视频| 两人在一起打扑克的视频| 97热精品久久久久久| 男人和女人高潮做爰伦理| 少妇人妻精品综合一区二区 | 成年版毛片免费区| 亚洲电影在线观看av| 国产极品精品免费视频能看的| av在线蜜桃| 日本 av在线| www.色视频.com| 亚洲熟妇熟女久久| 在线观看美女被高潮喷水网站| 嫩草影院入口| 简卡轻食公司| 91麻豆av在线| 国产又黄又爽又无遮挡在线| 免费一级毛片在线播放高清视频| 欧美人与善性xxx| 一级黄色大片毛片| 亚洲专区中文字幕在线| 热99在线观看视频| 国产男人的电影天堂91| 免费观看人在逋| 久99久视频精品免费| 亚洲天堂国产精品一区在线| 日本与韩国留学比较| 91久久精品电影网| 亚洲一级一片aⅴ在线观看| 欧美zozozo另类| 97超级碰碰碰精品色视频在线观看| 精品久久久久久久久久久久久| 波多野结衣巨乳人妻| 永久网站在线| 一进一出好大好爽视频| 免费搜索国产男女视频| 老司机福利观看| 亚洲在线观看片| 欧美日本视频| 天天一区二区日本电影三级| 动漫黄色视频在线观看| 人人妻,人人澡人人爽秒播| 亚洲av美国av| 最后的刺客免费高清国语| 国内久久婷婷六月综合欲色啪| 波野结衣二区三区在线| 免费av毛片视频| 国产亚洲精品久久久com| 噜噜噜噜噜久久久久久91| 精品欧美国产一区二区三| 色av中文字幕| 欧美+日韩+精品| 美女高潮喷水抽搐中文字幕| 一夜夜www| 精品人妻熟女av久视频| 在线观看舔阴道视频| 久久精品久久久久久噜噜老黄 | 亚洲欧美清纯卡通| 精品一区二区三区av网在线观看| 国产午夜福利久久久久久| 国产av不卡久久| 在线看三级毛片| 欧美日本视频| 成人三级黄色视频| 日本一二三区视频观看| 亚洲av电影不卡..在线观看| 国产免费男女视频| 国产一区二区三区av在线 | 搞女人的毛片| 最新中文字幕久久久久| 免费人成视频x8x8入口观看| 中文字幕免费在线视频6| 亚洲最大成人中文| 麻豆国产97在线/欧美| 麻豆av噜噜一区二区三区| 亚洲av中文字字幕乱码综合| av天堂在线播放| 免费高清视频大片| 极品教师在线视频| 国产精品亚洲一级av第二区| 国产精品,欧美在线| 嫁个100分男人电影在线观看| av专区在线播放| 男女下面进入的视频免费午夜| 一边摸一边抽搐一进一小说| 午夜精品在线福利| 国产在线精品亚洲第一网站| 中文字幕精品亚洲无线码一区| 国产精品嫩草影院av在线观看 | 99精品在免费线老司机午夜| 日本黄大片高清| 国产一区二区激情短视频| 国产v大片淫在线免费观看| 又紧又爽又黄一区二区| 丰满人妻一区二区三区视频av| 一本久久中文字幕| 欧美xxxx性猛交bbbb| 中文字幕免费在线视频6| 国内毛片毛片毛片毛片毛片| 日本与韩国留学比较| 日日干狠狠操夜夜爽| 欧美另类亚洲清纯唯美| 国产精品亚洲美女久久久| 国产爱豆传媒在线观看| 欧美成人性av电影在线观看| av黄色大香蕉| 久久精品91蜜桃| 欧美精品国产亚洲| 国模一区二区三区四区视频| 国产精品伦人一区二区| 一卡2卡三卡四卡精品乱码亚洲| 我要看日韩黄色一级片| 久久热精品热| 十八禁网站免费在线| 真人一进一出gif抽搐免费| 别揉我奶头 嗯啊视频| 国产精品乱码一区二三区的特点| 日韩欧美国产一区二区入口| 天堂√8在线中文| 99riav亚洲国产免费| 欧美极品一区二区三区四区| av在线亚洲专区| 波多野结衣高清无吗| 亚洲av成人av| 99久久精品国产国产毛片| 国产久久久一区二区三区| 熟妇人妻久久中文字幕3abv| 日韩欧美国产一区二区入口| 国产成人福利小说| 国产精品嫩草影院av在线观看 | 天堂影院成人在线观看| 一个人看的www免费观看视频| 一级a爱片免费观看的视频| 噜噜噜噜噜久久久久久91| 51国产日韩欧美| 少妇熟女aⅴ在线视频| 国产精品,欧美在线| 一进一出抽搐动态| 国产探花在线观看一区二区| 国产精品乱码一区二三区的特点| 日韩精品青青久久久久久| 淫秽高清视频在线观看| 一级av片app| 免费在线观看日本一区| 全区人妻精品视频| 3wmmmm亚洲av在线观看| 色综合站精品国产| 干丝袜人妻中文字幕| 搡老熟女国产l中国老女人| 久久热精品热| 亚洲狠狠婷婷综合久久图片| 国产亚洲欧美98| 伦理电影大哥的女人| 亚洲av日韩精品久久久久久密| 国产中年淑女户外野战色| 国产精品亚洲美女久久久| 亚洲综合色惰| 成人av在线播放网站| 欧美人与善性xxx| 最新中文字幕久久久久| 久久国内精品自在自线图片| 日韩大尺度精品在线看网址| 国产三级在线视频| 九九热线精品视视频播放| 欧美色欧美亚洲另类二区| 国产精品亚洲美女久久久| 久久这里只有精品中国| 国产麻豆成人av免费视频| 啦啦啦韩国在线观看视频| 日本黄大片高清| 亚洲最大成人av| 欧美极品一区二区三区四区| 国产午夜福利久久久久久| 亚洲精品国产成人久久av| 91精品国产九色| 在线观看免费视频日本深夜| 国产免费av片在线观看野外av| 成人午夜高清在线视频| 国产v大片淫在线免费观看| 少妇熟女aⅴ在线视频| 色噜噜av男人的天堂激情| 啪啪无遮挡十八禁网站| 日本三级黄在线观看| 麻豆成人av在线观看| 国产精品国产三级国产av玫瑰| 搡老熟女国产l中国老女人| 亚洲男人的天堂狠狠| 国语自产精品视频在线第100页| 亚洲av.av天堂| 三级毛片av免费| 久久香蕉精品热| 观看美女的网站| 啦啦啦观看免费观看视频高清| av国产免费在线观看| 最近最新免费中文字幕在线| 国产在线精品亚洲第一网站| 国产免费一级a男人的天堂| 成年女人永久免费观看视频| 一级黄片播放器| 91狼人影院| 3wmmmm亚洲av在线观看| 精品久久久久久久末码| 久久久久久久亚洲中文字幕| 国产免费一级a男人的天堂| a级毛片免费高清观看在线播放| 国产视频内射| 亚洲va日本ⅴa欧美va伊人久久| 日韩大尺度精品在线看网址| 亚洲国产精品成人综合色| 亚洲不卡免费看| 女同久久另类99精品国产91| 色吧在线观看| 免费观看的影片在线观看| 国产激情偷乱视频一区二区| av天堂在线播放| 丝袜美腿在线中文| 十八禁国产超污无遮挡网站| 久久久久久九九精品二区国产| 亚洲性夜色夜夜综合| 老司机福利观看| 亚洲电影在线观看av| 欧美日韩中文字幕国产精品一区二区三区| 亚洲avbb在线观看| 男人舔奶头视频| 国产亚洲精品久久久com| 国内精品一区二区在线观看| 亚洲午夜理论影院| 1000部很黄的大片| 女人十人毛片免费观看3o分钟| 日韩一区二区视频免费看| 88av欧美| 亚洲人与动物交配视频| 99视频精品全部免费 在线| 久久久久久大精品| 蜜桃久久精品国产亚洲av| 成人av一区二区三区在线看| or卡值多少钱| 麻豆国产97在线/欧美| 校园春色视频在线观看| 国产成人一区二区在线| 亚洲,欧美,日韩| 色综合亚洲欧美另类图片| 51国产日韩欧美| 日韩欧美国产一区二区入口| 嫩草影院精品99| 老司机午夜福利在线观看视频| 国产免费一级a男人的天堂| 88av欧美| 中文字幕免费在线视频6| 日本熟妇午夜| 九九热线精品视视频播放| 男女那种视频在线观看| 欧美日韩精品成人综合77777| 麻豆av噜噜一区二区三区| 亚洲专区中文字幕在线| 欧美xxxx黑人xx丫x性爽| 国产免费男女视频| 男女边吃奶边做爰视频| 日本一本二区三区精品| 高清日韩中文字幕在线| 老司机午夜福利在线观看视频| 男女啪啪激烈高潮av片| 国产精华一区二区三区| 国产精品亚洲美女久久久| 少妇裸体淫交视频免费看高清| 国产亚洲精品综合一区在线观看| 久久久久久久久久黄片| 欧美另类亚洲清纯唯美| 村上凉子中文字幕在线| 午夜福利视频1000在线观看| 午夜福利在线观看免费完整高清在 | 国产午夜精品久久久久久一区二区三区 | 国产aⅴ精品一区二区三区波| 日韩欧美国产一区二区入口| 国产欧美日韩一区二区精品| 国产精品98久久久久久宅男小说| 日韩欧美免费精品| 可以在线观看毛片的网站| 十八禁网站免费在线| 国产av麻豆久久久久久久| 亚洲乱码一区二区免费版| 国产麻豆成人av免费视频| 国产一级毛片七仙女欲春2| 亚洲性夜色夜夜综合| 亚洲中文日韩欧美视频| 哪里可以看免费的av片| 一本精品99久久精品77| 久久精品久久久久久噜噜老黄 | 中文字幕人妻熟人妻熟丝袜美| 日日撸夜夜添| 午夜福利成人在线免费观看| 美女大奶头视频| 在线免费观看不下载黄p国产 | 美女高潮的动态| 此物有八面人人有两片| 99精品久久久久人妻精品| 亚洲中文字幕日韩| 国产亚洲精品久久久com| 老司机深夜福利视频在线观看| 丰满乱子伦码专区| 真人做人爱边吃奶动态| 深夜精品福利| 国产精品一区二区免费欧美| 国产真实伦视频高清在线观看 | 日本与韩国留学比较| 中文字幕熟女人妻在线| 精华霜和精华液先用哪个| 成人性生交大片免费视频hd| 麻豆久久精品国产亚洲av| 午夜激情欧美在线| av福利片在线观看| 国产精品久久久久久av不卡| av福利片在线观看| 日本免费a在线| 九九热线精品视视频播放| 久久久久久久午夜电影| 久久国产精品人妻蜜桃| 免费av毛片视频| 极品教师在线视频| 色综合亚洲欧美另类图片| 日本 欧美在线| 特级一级黄色大片| a级毛片a级免费在线| 啦啦啦观看免费观看视频高清| 久久久久久大精品| 中文资源天堂在线| 亚洲性夜色夜夜综合| 欧美人与善性xxx| 极品教师在线免费播放| 亚洲成av人片在线播放无| 亚洲无线在线观看| 日韩国内少妇激情av|