寧德志，蘇曉杰，滕斌

（1.大連理工大學(xué)海岸和近海工程國(guó)家重點(diǎn)實(shí)驗(yàn)室，遼寧大連 116024）

波浪與帶窄縫多箱體作用共振現(xiàn)象的模擬研究

寧德志1，蘇曉杰1，滕斌1

（1.大連理工大學(xué)海岸和近海工程國(guó)家重點(diǎn)實(shí)驗(yàn)室，遼寧大連 116024）

針對(duì)波浪與帶有窄縫多箱體結(jié)構(gòu)作用產(chǎn)生的流體共振問(wèn)題，建立了基于域內(nèi)源造波技術(shù)的二維非線性時(shí)域數(shù)值波浪水槽模型，其中自由水面滿(mǎn)足完全非線性運(yùn)動(dòng)學(xué)和動(dòng)力學(xué)邊界條件，窄縫內(nèi)流體引入人工阻尼來(lái)等效由于渦旋運(yùn)動(dòng)和流動(dòng)分離引起的黏性耗散，計(jì)算域邊界采用高階邊界元進(jìn)行離散。通過(guò)模擬三箱體間兩窄縫內(nèi)相對(duì)波高變化，并與已發(fā)表的數(shù)值與實(shí)驗(yàn)結(jié)果對(duì)比，驗(yàn)證了本模型的準(zhǔn)確性。同時(shí)通過(guò)大量的數(shù)值計(jì)算，分析了箱體數(shù)量對(duì)窄縫內(nèi)水體共振頻率、共振波高以及對(duì)結(jié)構(gòu)反射波高和透射波高的影響。

窄縫；流體共振；域內(nèi)源造波技術(shù)；高階邊界元；非線性數(shù)值波浪水槽

1 引言

為了充分利用海洋空間資源，由多模塊組成的超大型浮體已在海洋工程領(lǐng)域得到利用［1］，譬如用作海上機(jī)場(chǎng)、軍事基地、儲(chǔ)存器和海難救助點(diǎn)等，國(guó)內(nèi)外學(xué)者已經(jīng)開(kāi)展了很多有關(guān)超大型浮體水動(dòng)力特性的研究，如波浪沖擊［2］，水彈性運(yùn)動(dòng)響應(yīng)［3］，非穩(wěn)定外荷載作用下的動(dòng)力響應(yīng)［4］等。由于組成超大型浮體的各個(gè)模塊之間并不是無(wú)縫連接，大都存在相對(duì)模塊特征長(zhǎng)度很小的窄縫，縫隙內(nèi)的水體在某些頻率波浪作用下會(huì)發(fā)生共振現(xiàn)象，誘發(fā)很大的波浪爬高和荷載，進(jìn)而對(duì)海洋結(jié)構(gòu)的作業(yè)安全帶來(lái)很大的影響。

有關(guān)海洋結(jié)構(gòu)物間窄縫內(nèi)波浪水動(dòng)力特性問(wèn)題，國(guó)內(nèi)外已經(jīng)開(kāi)展了許多相關(guān)理論、模型試驗(yàn)和數(shù)值模擬研究，但大多數(shù)研究目前還只是局限于兩個(gè)結(jié)構(gòu)或3個(gè)結(jié)構(gòu)間縫隙水體共振問(wèn)題。譬如，Miao等［5］采用漸近匹配法研究了帶狹縫二維雙箱的共振現(xiàn)象，給出了狹縫很小時(shí)雙箱的理論共振頻率。Saitoh等［6］對(duì)不同入射波浪作用下兩個(gè)方箱間窄縫的波高變化進(jìn)行了試驗(yàn)研究，發(fā)現(xiàn)窄縫內(nèi)最大共振波高可以達(dá)到入射波高的5倍。Kristiansen和Faltinsen［7—8］對(duì)波浪與具有窄縫的二維固定及浮式箱體和固定式岸壁結(jié)構(gòu)之間的相互作用分別開(kāi)展了數(shù)值和試驗(yàn)?zāi)M，發(fā)現(xiàn)線性和非線性勢(shì)流數(shù)值結(jié)果都比共振條件下的窄縫內(nèi)試驗(yàn)波高要大，尤其是線性結(jié)果更大；并研究了共振條件下浮體3個(gè)自由度的運(yùn)動(dòng)響應(yīng)與水動(dòng)力系數(shù)。在三體結(jié)構(gòu)雙縫隙方面，Iwata等［9］分別對(duì)不同入射波浪作用下3個(gè)方箱之間窄縫的波高變化進(jìn)行了試驗(yàn)研究，發(fā)現(xiàn)共振頻率與箱體的吃水深度、窄縫寬度和箱體個(gè)數(shù)成一定的函數(shù)關(guān)系。何廣華等［10］采用比例邊界有限元方法分別研究了波浪與三箱結(jié)構(gòu)作用下窄縫內(nèi)水體共振現(xiàn)象。Lu等［11—12］采用黏性流模型和改進(jìn)的線性勢(shì)流模型對(duì)三體結(jié)構(gòu)間窄縫內(nèi)水體共振引起的波浪爬高和波浪荷載進(jìn)行了研究，并發(fā)現(xiàn)在窄縫內(nèi)自由水面上布置一定的人工阻尼時(shí)，勢(shì)流模型也可以得到與試驗(yàn)和黏性流模型相一致的結(jié)果。通過(guò)上述研究發(fā)現(xiàn)，三箱體間窄縫內(nèi)水體共振頻率明顯不同于雙箱浮體情況，甚至出現(xiàn)多個(gè)自振頻率等復(fù)雜情況，也進(jìn)一步說(shuō)明箱體數(shù)量對(duì)窄縫內(nèi)水體運(yùn)動(dòng)特性具有重要影響，而目前對(duì)于3個(gè)及以上浮體間窄縫內(nèi)水體共振問(wèn)題的研究還很少。

本文將在前人研究基礎(chǔ)上，建立自由水面滿(mǎn)足完全非線性邊界條件的二維時(shí)域數(shù)值波浪水槽模型，采用域內(nèi)源造波方法產(chǎn)生入射波浪并布置前置阻尼層消除二次反射影響，實(shí)現(xiàn)在較小計(jì)算域內(nèi)進(jìn)行長(zhǎng)時(shí)間模擬；參照Lu等［11］的方法在窄縫水體自由水面上引入常人工阻尼系數(shù)來(lái)等效共振條件下黏性耗散。通過(guò)模擬波浪與3箱體間雙窄縫內(nèi)水體運(yùn)動(dòng)問(wèn)題，并與Iwata等［9］實(shí)驗(yàn)結(jié)果及Lu等［12］的CFD數(shù)值結(jié)果進(jìn)行對(duì)比，證明了本模型的準(zhǔn)確性。進(jìn)而通過(guò)大量數(shù)值計(jì)算研究了箱體數(shù)量（最多至6個(gè)）對(duì)窄縫內(nèi)水體共振頻率、共振波高以及對(duì)結(jié)構(gòu)迎浪側(cè)反射波高和背浪側(cè)透射波高的影響，總結(jié)一般性規(guī)律，為如超大型浮體這樣由多模塊組成的具有窄縫的海洋結(jié)構(gòu)水動(dòng)力分析提供參考。

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

考慮單向規(guī)則波浪與具有窄縫的多固定箱體相互作用問(wèn)題，其布置如圖1所示。建立二維笛卡爾坐標(biāo)系Oxz，z＝0位于靜水面上，且z軸向上為正，x軸右方向?yàn)檎Ｓ?jì)算域包含自由水面Гf和固體邊界ГN（包括水底Гd和箱體Гb）。波浪由控制垂直源造波面（Гs）的流量密度產(chǎn)生。圖中h為水槽靜水深，W為箱體寬度，D為箱體吃水深度，Wg為兩箱體間窄縫的寬度，N為箱體個(gè)數(shù)，本研究中取值為2、3、4、5、6?？紤]問(wèn)題一般性，本文中假定各個(gè)箱體的寬度、吃水深度和縫隙寬度均一致，圖1中迎浪側(cè)e點(diǎn)、背浪側(cè)f點(diǎn)和各個(gè)箱體件窄縫（gap）內(nèi)波高變化規(guī)律是本文研究的重點(diǎn)。在流體無(wú)黏、不可壓縮和流動(dòng)無(wú)旋的假定下，整個(gè)流域的速度可用速度勢(shì)的梯度來(lái)描述，上述問(wèn)題的控制方程為由速度勢(shì)滿(mǎn)足的泊松（Poisson）方程［13］，即

式中，q*（xs，z，t）＝2Vδ（x－xs）為造波源強(qiáng)度；造波位置x＝xs（本文均取xs＝0）；V為流體質(zhì)點(diǎn)水平速度，本文給定二階Stokes速度解析解。

在自由水面上，滿(mǎn)足完全非線性動(dòng)力學(xué)和運(yùn)動(dòng)學(xué)邊界條件，本文中采用混合歐拉-拉格朗日方法更新自由水面，利用物質(zhì)導(dǎo)數(shù)，并在計(jì)算域的上游和下游區(qū)域的自由水面分別布置人工阻尼層來(lái)吸收從結(jié)構(gòu)物反射回來(lái)的波浪與出流波浪［14］，在窄縫內(nèi)布置一常參數(shù)人工阻尼來(lái)近似由于渦和分流引起的黏性耗散［15］。自由水面邊界條件可以寫(xiě)成以下形式：

式中，η代表自由水面的鉛垂位移，g是重力加速度，X0＝（x0，0）是指水質(zhì)點(diǎn)初始靜止時(shí)的位置，阻尼系數(shù)

用于計(jì)算域邊界的兩個(gè)阻尼層，x1和x2分別是左右兩側(cè)阻尼層的起點(diǎn)位置；L為阻尼層長(zhǎng)度，本文取1.5倍波長(zhǎng)（即1.5λ）；ω是波浪角頻率。阻尼系數(shù)μ2用于窄縫內(nèi)自由水面，其數(shù)值根據(jù)試驗(yàn)黏性耗散來(lái)確定；k是波數(shù)，滿(mǎn)足如下線性色散方程關(guān)系

在固定的結(jié)構(gòu)表面和底面邊界上，流體法向速度為0，滿(mǎn)足固壁不可滲透邊界條件。由于本研究在時(shí)域內(nèi)進(jìn)行，自由水面滿(mǎn)足靜初始條件，即起始時(shí)刻速度勢(shì)和波面均為0。

對(duì)于滿(mǎn)足上述控制方程和邊界條件的定解問(wèn)題，在整個(gè)流域內(nèi)對(duì)速度勢(shì)應(yīng)用格林第二定理，可轉(zhuǎn)換如下邊界積分方程［13］：

式中，p＝（x0，z0）為源點(diǎn)；q＝（x，z）為場(chǎng)點(diǎn)；C為固角系數(shù)；Ω代表整個(gè)流域；G是簡(jiǎn)單格林函數(shù)，考慮到水底鏡像，可以表示為如下形式：

式中，r1為p和q兩點(diǎn)距離，

r2為p和q關(guān)于水底鏡像之間距離，

本文用三節(jié)點(diǎn)高階邊界元離散計(jì)算域成一些曲線單元，單元內(nèi)任一點(diǎn)的幾何坐標(biāo)和速度勢(shì)等物理量可以用如下二次形狀函數(shù)hi（ξ）插值得到，

式中，ξ代表固有坐標(biāo)，取值范圍（－1.0，1.0）。

這樣在任一邊界單元內(nèi)，物理量和幾何量都可以通過(guò)形狀函數(shù)插值得到，也即，邊界單元都是等參的。積分方程（4）經(jīng)高階邊界元離散后，可以表示成如下形式：

式中，Ne1、Ne2、Ne3分別為自由水面、物面和出入流邊界及造波源面上劃分的單元個(gè)數(shù)，J（ξ）是聯(lián)系大地坐標(biāo)和固有坐標(biāo)的雅可比行列式。

由于本文采用了高階邊界元方法，自由水面邊界條件式（2）中用到的速度勢(shì)的空間導(dǎo)數(shù)項(xiàng)進(jìn)而可以表示成如下形式：

式中，n＝（nx，nz）為單位法向量，指出水體為正。

最后把未知量都移到方程（9）的左側(cè)，該方程組就可以寫(xiě)出如下矩陣乘積形式：

式中，X是未知的速度勢(shì)和速度勢(shì)法向?qū)?shù)，A為空間系數(shù)矩陣，F(xiàn)為由已知的速度勢(shì)和速度在邊界上積分得到的列向量。

計(jì)算中認(rèn)為當(dāng)前時(shí)刻物面上的速度勢(shì)法向?qū)?shù)和自由水面上的速度勢(shì)是已知的，根據(jù)積分方程計(jì)算當(dāng)前時(shí)刻物面上的速度勢(shì)和自由水面上的速度勢(shì)法向?qū)?shù)，然后應(yīng)用四階Runga-Kutta法，根據(jù)自由水面條件式（2）計(jì)算下一時(shí)刻的水質(zhì)點(diǎn)位置和自由水面上的速度勢(shì)，再用二次形狀函數(shù)在舊單元上插值求得新節(jié)點(diǎn)上的物理量來(lái)對(duì)自由水面網(wǎng)格重新劃分，重新應(yīng)用積分方程計(jì)算下一時(shí)刻物面上的速度勢(shì)和自由水面上的速度勢(shì)法向?qū)?shù)。這樣計(jì)算周而復(fù)始，直到計(jì)算結(jié)束［16—17］。

3 數(shù)值計(jì)算及討論

3.1 模型準(zhǔn)確性和穩(wěn)定性

作為算例，本文以Iwata等［9］的實(shí)驗(yàn)來(lái)進(jìn)行數(shù)值模擬波浪與具有窄縫的三箱體的相互作用的研究，驗(yàn)證本文數(shù)學(xué)模型的準(zhǔn)確性。這里選用的實(shí)驗(yàn)參數(shù)為水槽靜水深h＝0.5 m，箱體的寬度W＝0.5 m，吃水深度D＝0.252 m，入射波高H0＝0.024 m，箱體間窄縫寬度0.05 m。在數(shù)值模型中，計(jì)算域長(zhǎng)度取7.5倍波長(zhǎng)，在水槽的左右兩端各布置1.5倍波長(zhǎng)的阻尼層，造波源位于x＝0，箱體1左側(cè)面邊界位于距離造波源2.5倍波長(zhǎng)的位置，然后依次按Wg調(diào)整箱體2和3的位置。通過(guò)開(kāi)展數(shù)值收斂性實(shí)驗(yàn)，自由水面上每個(gè)波長(zhǎng)布置15個(gè)單元，窄縫內(nèi)布置2個(gè)單元，計(jì)算域垂向邊界和造波源面分別布置10個(gè)單元，箱體側(cè)面邊界上均布置6個(gè)單元，底面邊界上均布置12個(gè)單元；時(shí)間步長(zhǎng)Δt＝T／60 s，每個(gè)算例模擬30個(gè)周期。

圖2給出了窄縫寬度Wg為0.05 m情況下兩個(gè)窄縫中心位置無(wú)量綱波高Hg／H0隨入射波波數(shù)kh的變化關(guān)系，及本文數(shù)值結(jié)果與實(shí)驗(yàn)數(shù)據(jù)［9］、黏性流模型數(shù)值結(jié)果［12］的比較。由于窄縫內(nèi)水體黏性與波浪條件無(wú)關(guān)，故在某一波浪頻率下通過(guò)取一系列不同人工黏性系數(shù)進(jìn)行數(shù)值模擬，并與實(shí)驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比，確定窄縫中間水面的人工黏性系數(shù)取為μ2＝0.03；數(shù)值模擬中，波高Hg以測(cè)試點(diǎn)的穩(wěn)定時(shí)間段內(nèi)波面時(shí)間序列中波峰值與波谷絕對(duì)值和的平均來(lái)計(jì)算。從圖中可以看出，箱體1與箱體2之間窄縫內(nèi)數(shù)值模擬波面高在kh＝1.35時(shí)達(dá)到最大，也即窄縫內(nèi)流體發(fā)生共振，波高為入射波高的4.6倍。箱體2與箱體3之間窄縫內(nèi)數(shù)值模擬波面高在kh＝1.4時(shí)達(dá)到最大，波高為入射波高的4.4倍。同時(shí)發(fā)現(xiàn)窄縫1內(nèi)出現(xiàn)2個(gè)共振頻率，這與單縫隙內(nèi)的共振規(guī)律是不同的，也進(jìn)一步說(shuō)明浮體數(shù)量對(duì)窄縫內(nèi)水體運(yùn)動(dòng)規(guī)律的重要影響。整體上兩種數(shù)值結(jié)果與實(shí)驗(yàn)數(shù)據(jù)均符合的很好，在個(gè)別位置處甚至本文結(jié)果比黏性流模型結(jié)果與實(shí)驗(yàn)數(shù)據(jù)吻合的更好，說(shuō)明所建立模型在取得合適的人工黏性系數(shù)情況下可以準(zhǔn)確模擬多箱體窄縫內(nèi)流體共振問(wèn)題。

圖3給出了波數(shù)kh＝1.40，窄縫寬度Wg＝0.05 m，入射波高H0＝0.024 m情況下t＝26T和30T時(shí)的整個(gè)計(jì)算域波面分布，圖中虛線分界處分別對(duì)應(yīng)3個(gè)箱體所在的位置。從圖中可以看出，兩個(gè)時(shí)刻的波面曲線已經(jīng)完全重合，包括箱體1前的反射波，箱體3后的透射波和3個(gè)箱體之間的窄縫內(nèi)波面；并且兩端的波面基本趨于0，說(shuō)明兩端阻尼層吸收波浪的效果很理想；箱體1前的波浪形成了穩(wěn)定的立波，說(shuō)明其反射回去的波浪透過(guò)造波源被前端阻尼層完全吸收，而對(duì)入射波浪沒(méi)有產(chǎn)生影響。以上現(xiàn)象說(shuō)明本模型的模擬結(jié)果已達(dá)到穩(wěn)定。

圖2 兩窄縫中無(wú)因次波高Hg／H0與波數(shù)kh間的關(guān)系Fig.2 Distribution of non-dimensional wave height versus incident wave number at two narrow gaps

圖3t＝26T和30T兩個(gè)時(shí)刻的水槽波面分布Fig.3 Snapshot of wave elevation along the wave flume att＝26Tand 30T

3.2 數(shù)值結(jié)果

下面仍以上述工況為例，保持各參數(shù)不變，只是改變箱體數(shù)量，進(jìn)一步分析箱體數(shù)量對(duì)窄縫內(nèi)水體共振頻率、波浪爬高以及結(jié)構(gòu)前反射波高和結(jié)構(gòu)后透射波高的影響。

圖4和圖5是箱體個(gè)數(shù)N分別為4和5時(shí)各個(gè)窄縫內(nèi)無(wú)因次波高Hg／H0隨入射波波數(shù)kh的變化關(guān)系。各窄縫的主共振頻率相仿，其所對(duì)應(yīng)的共振波高在處于中間位置的窄縫內(nèi)達(dá)到最大。在大于主共振頻率的某一頻率處，還會(huì)有次共振現(xiàn)象發(fā)生，且處于兩端位置的窄縫內(nèi)較中間位置的窄縫內(nèi)次共振現(xiàn)象更為明顯，尤為突出的是迎浪側(cè)的窄縫內(nèi)，在圖5的窄縫1中高共振頻率所對(duì)應(yīng)的共振波高已大于其低共振頻率所對(duì)應(yīng)的波高。其原因可能是各窄縫內(nèi)自振頻率下的波浪會(huì)向兩側(cè)傳播并與迎浪側(cè)透射的波浪相互作用，處于中間位置的窄縫由于對(duì)稱(chēng)關(guān)系，從窄縫兩側(cè)反射來(lái)的波浪疊加達(dá)到最大，而處于邊界的窄縫由于其相對(duì)位置關(guān)系則會(huì)發(fā)生復(fù)雜的波浪干涉現(xiàn)象。

圖4 4箱體時(shí)各窄縫內(nèi)波高隨波數(shù)的變化關(guān)系Fig.4 Dimensionless wave height againstkhat various gaps for four boxes

圖6和圖7分別為不同數(shù)量箱體時(shí)窄縫1（即迎浪側(cè)第一個(gè)窄縫）和窄縫N－1（背浪側(cè)第一個(gè)窄縫）內(nèi)無(wú)因次波高隨波數(shù)的分布情況。可以看出隨著箱體數(shù)量的增加，主頻共振發(fā)生時(shí)，窄縫內(nèi)波高減小，對(duì)應(yīng)的共振頻率向低頻偏移。并且箱體數(shù)量較多時(shí)，窄縫會(huì)在多個(gè)頻率發(fā)生共振，且高共振頻率處的波高隨箱體個(gè)數(shù)的增加而增大，窄縫1內(nèi)在箱體數(shù)N為5和6時(shí)已大于低共振頻率處的波高，窄縫5內(nèi)高頻共振波高也大于低頻共振波高。各個(gè)工況下窄縫N－1內(nèi)主共振頻率都不小于窄縫1內(nèi)主共振頻率，而窄縫1和N－1內(nèi)相鄰箱體數(shù)對(duì)應(yīng)的次頻共振頻率的差值近似常數(shù)（Δkh≈0.1）。

圖5 5箱體時(shí)各窄縫內(nèi)波高隨波數(shù)的變化關(guān)系Fig.5 Dimensionless wave height againstkhat various gaps for five boxes

圖6 窄縫1內(nèi)無(wú)因次波高隨波數(shù)的變化關(guān)系Fig.6 Dimensionless wave height againstkhat gap 1

圖7 窄縫N內(nèi)無(wú)因次波高隨波數(shù)的變化關(guān)系Fig.7 Dimensionless wave height againstkhat gapN

表1和表2列出了不同箱體數(shù)情況下各個(gè)窄縫內(nèi)對(duì)應(yīng)共振頻率和共振波高分布情況。可以看出，除了得到與圖4、圖5、圖6和圖7相同的規(guī)律外，還會(huì)發(fā)現(xiàn)對(duì)于同一箱體個(gè)數(shù)N，各窄縫內(nèi)的低共振頻率會(huì)從迎浪側(cè)向背浪側(cè)逐漸增大，對(duì)應(yīng)的共振波高則是中間窄縫內(nèi)最大，向迎浪側(cè)和背浪側(cè)兩側(cè)減小；而高共振頻率則是處于中間位置窄縫內(nèi)最大，向迎浪側(cè)和背浪側(cè)兩側(cè)減小，對(duì)應(yīng)的共振波高則是中間位置窄縫內(nèi)最小，向兩側(cè)增大。

表1 各窄縫對(duì)應(yīng)的低共振波數(shù)和共振波高（kh，Hg／H0）Tab.1 Low resonant wave number and wave height at various gaps（kh，Hg／H0）

表2 各窄縫對(duì)應(yīng)的高共振波數(shù)和共振波高（kh，Hg／H0）Tab.2 High resonant wave number and wave height at various gaps（kh，Hg／H0）

圖8 迎浪側(cè)e點(diǎn)波高隨波數(shù)的變化關(guān)系Fig.8 Wave height at the weather-side pointeagainstkh

圖8和圖9分別給出了迎浪側(cè)e點(diǎn)和背浪側(cè)f點(diǎn)無(wú)因次波高隨波數(shù)的變化關(guān)系。從圖8中可以看出，在窄縫1內(nèi)主頻共振頻率處，迎浪側(cè)e點(diǎn)波高出現(xiàn)最小值，且箱體數(shù)N越大，對(duì)應(yīng)的極值越?。∟＝6時(shí)，最小的Hg／H0＝1.4），在低頻處無(wú)因次波高接近于2，相當(dāng)于波浪線性化并被全反射，而在高頻處無(wú)因次波高會(huì)大于2，這是由于波浪的非線性增強(qiáng)，高階諧波貢獻(xiàn)所致。從圖9可以看出，在窄縫N內(nèi)主頻共振頻率處，背浪側(cè)f點(diǎn)波高也出現(xiàn)明顯的峰值，但其變化規(guī)律與圖8相反，隨著箱體數(shù)N的增大，透射波高的峰值也增大（N＝6時(shí)，最大的Hg／H0＝0.78），在高頻處，透射波高接近于0。

圖9 背浪側(cè)f點(diǎn)波高隨波數(shù)的變化關(guān)系Fig.9 Wave height at the lee-side pointfagainstkh

4 結(jié)論

本文基于域內(nèi)源造波的時(shí)域高階邊界元方法建立波浪與具有窄縫的多箱體結(jié)構(gòu)相作用的完全非線性數(shù)值水槽模型，對(duì)不同個(gè)數(shù)箱體時(shí)窄縫內(nèi)流體的共振頻率、共振波高及結(jié)構(gòu)迎浪側(cè)波高和背浪側(cè)透射波高等進(jìn)行了模擬研究。通過(guò)與已發(fā)表實(shí)驗(yàn)數(shù)據(jù)和數(shù)值結(jié)果進(jìn)行對(duì)比驗(yàn)證，表明本文所建立數(shù)學(xué)模型可以準(zhǔn)確模擬波浪與具有窄縫的多箱體相互作用過(guò)程，且在較小計(jì)算域內(nèi)可長(zhǎng)時(shí)間模擬得到穩(wěn)定的結(jié)果，沒(méi)有在入射邊界發(fā)生二次反射現(xiàn)象。研究發(fā)現(xiàn)：箱體寬度、吃水、窄縫寬度一定時(shí)，不同數(shù)量箱體，窄縫內(nèi)水體發(fā)生共振的頻率是不同的且箱體個(gè)數(shù)較多時(shí)會(huì)有多個(gè)共振頻率。隨著箱體數(shù)量的增加，主頻共振發(fā)生時(shí)，窄縫內(nèi)波高減小，主頻共振頻率值減?。淮晤l共振頻率處的波高隨箱體個(gè)數(shù)的增加而增大。對(duì)于一給定箱體數(shù)N的情況，各窄縫內(nèi)的主頻共振頻率會(huì)從迎浪側(cè)向背浪側(cè)逐漸增大，對(duì)應(yīng)的共振波高則是中間窄縫內(nèi)最大，向迎浪側(cè)和背浪側(cè)兩側(cè)減?。欢喂舱耦l率則從處于中間位置窄縫內(nèi)最大，向迎浪側(cè)和背浪側(cè)兩側(cè)減小，對(duì)應(yīng)的共振波高則是中間位置窄縫內(nèi)最小，向兩側(cè)增大。迎浪側(cè)波高在窄縫1內(nèi)共振發(fā)生時(shí)最小，背浪側(cè)波高則在窄縫N－1內(nèi)共振發(fā)生時(shí)出現(xiàn)峰值；隨著箱體個(gè)數(shù)的增加，共振發(fā)生時(shí)迎浪側(cè)波高減小，背浪側(cè)波高增大。

［1］ 繆國(guó)平，劉應(yīng)中.征服海洋之夢(mèng)——超大型海洋浮式結(jié)構(gòu)物［J］.自然雜志，1996，18（1）：26－30.

Miao Guoping，Liu Yingzhong.A dream to conquer the ocean——super large floating ocean structures［J］.Ziran Zazhi，1996，18（1）：26－30.

［2］ Yashimoto H，Ohmatsu S，Hoshino K，et al.Slamming load on a very large floating body with a shallow drift［J］.Journal of Marine Science and Technology，1997，2：163－172.

［3］ Kagemoto H，F(xiàn)ujino M，Murai M.Theoretical and experimental predictions of the hydro-elastic response of a very large floating structure in waves［J］.Applied Ocean Research，1998，20：135－144.

［4］ Qiu L，Liu H.Three-dimensional time-domain analysis of very large floating structures subjected to unsteady external loading［J］.Journal of Offshore Mechanics and Arctic Engineering，2006，129（1）：21－28.

［5］ Miao G P，Ishida H，Saitoh T.Influence of gaps between multiple floating bodies onwaveforces［J］.China Ocean Engineering，2000，14（4）：407－422.

［6］ Saitoh T，Miao G P，Ishida H.Theoretical analysis on appearancecondition of fluid resonance in a narrow gap between two modulesof very large floating structure［C］／／Proceedings of the Third Asia－Pacific Workshop on Marine Hydrodynamics，Shanghai，China，2006：170－175.

［7］ Kristiansen T，F(xiàn)altinsen O M.Studies on resonant water motion between a ship and a fixed terminal in shallow water［J］.Journal of Offshore Mechanics and Arctic Engineering，2009，131：021102.

［8］ Kristiansen T，F(xiàn)altinsen O M.A two-dimensional numerical andexperimental study of resonant coupled ship and piston-modemotion［J］.Applied O-cean Research，2010，32：158－176.

［9］ Iwata H，Saitoh T，Miao G P.Fluid resonance in narrow gaps ofvery large floating structure composed of rectangular modules［C］／／Proceedings of the Fourth International Conference on Asian and Pacific Coasts，Nanjing，China，2007：815－826.

［10］ 何廣華，滕斌，李博寧，等.應(yīng)用比例邊界有限元法研究波浪與帶狹縫三箱作用的共振現(xiàn)象［J］.水動(dòng)力學(xué)研究與進(jìn)展，2006，21（3）：418－424.

He Guanghua，Teng Bin，Li Boning，et al.Research on the hydrodynamic influence from the gaps between threeidentical boxes by a scaled boundary finite element method［J］.Journal of Hydrodynamics，2006，21（3）：418－424.

［11］ Lu L，Cheng L，Teng B，et al.Numerical investigation of fluid resonance in two narrow gaps of three identical rectangular structures［J］.Applied Ocean Research，2010，32（2）：177－190.

［12］ Lu L，Teng B，Cheng L，et al.Modelling of multi-bodies in close proximity under water waves——fluid resonance in narrow gaps［J］.Science China Physics，Mechanics＆Astronomy，2011，54（1）：16－25.

［13］ Ning D Z，Teng B，Eatoack Taylor R，et al.Numerical simulation of non-linear regular and focused waves in an infinite water-depth［J］.Ocean Engineering，2008，35（8／9）：887－899.

［14］ Tanizawa K.Long time fully nonlinear simulation of floating body motions with artificial damping zone［J］.The Society of Naval Architects of Japan，1996，180：311－319.

［15］ Kim Y W.Artificial damping in water wave problems：Ⅰ.Constant damping［J］.International Journal of Offshore and Polar Engineering，2003，13（2）：88－93.

［16］ 周斌珍，寧德志，滕斌.造波板運(yùn)動(dòng)造波實(shí)時(shí)模擬［J］.水動(dòng)力學(xué)研究與進(jìn)展，2009，24（4）：1－12.

Zhou Binzhen，Ning Dezhi，Teng Bin.Real-time simulation of waves generated by a wave maker［J］.Journal of Hydrodynamics，2009，24（4）：1－12.

［17］ 陳麗芬，寧德志，滕斌，等.潛堤上波流傳播的完全非線性數(shù)值模擬［J］.力學(xué)學(xué)報(bào)，2011，43（5）：834－843.

Chen Lifen，Ning Dezhi，Teng Bin，et al.Full nonlinear numerical simulation for wave-current propagation over a submerged bar［J］.Chinese Journal of Theoretical and Applied Mechanics，2011，43（5）：834－843.

Numerical study of fluid resonance induced by wave action on multi-boxes with narrow gaps

Ning Dezhi1，Su Xiaojie1，Teng Bin1

（1.State Key Laboratory of Coastal and Offshore Engineering，Dalian University of Technology，Dalian 116024，China）

Based on wave generation technique by a inner-domain source，a two-dimensional nonlinear numerical wave flume over time is developed to investigate the fluid resonance induced by the interaction between wave and multi-objects with narrow gaps.In the numerical model，the fully nonlinear kinematic and dynamic boundary conditions are set for the instantaneous free surface；the artificial damping is introduced into the gap to simulate the viscous dissipation due to vortex motion and flow separation；the computational domain is discretized using higher-order boundary elements.The proposed model is validated by the published experimental and numerical data of the relative wave height at two narrow gaps of three boxes.Numerical experiments are performed to study the following：the effects of the number of the boxes on the resonant frequency，wave height at various gaps and the reflected and transmitted wave heights of the objects.

narrow gap；fluid resonance；inner-domain source generation technique；higher-order boundary element；nonlinear numerical wave flume

O353.2

A

0253-4193（2015）03-0126-08

寧德志，蘇曉杰，滕斌.波浪與帶窄縫多箱體作用共振現(xiàn)象的模擬研究［J］.海洋學(xué)報(bào)，2015，37（3）：126—133，

10.3969／j.issn.0253-4193.2015.03.013

Ning Dezhi，Su Xiaojie，Teng Bin.Numerical study of fluid resonance induced by wave action on multi-boxes with narrow gaps［J］.Haiyang Xuebao，2015，37（3）：126—133，doi：10.3969／j.issn.0253-4193.2015.03.013

2014-01-22；

2014-04-21。

國(guó)家自然科學(xué)基金項(xiàng)目（51179028，51222902，51221961）；教育部新世紀(jì)優(yōu)秀人才支持計(jì)劃（NCET－13－0076）；中央高?；究蒲袠I(yè)務(wù)費(fèi)專(zhuān)項(xiàng)資金（DUT13YQ104）。

寧德志（1975—），男，黑龍江省五常市人，教授，博士生導(dǎo)師，主要從事非線性波浪及其與結(jié)構(gòu)物相互作用的研究。E-mail：dzning＠dlut.edu.cn