淺水中浮體在波浪作用下運(yùn)動(dòng)的三維時(shí)域計(jì)算模型

張俊生，滕 斌，叢培文

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

淺水中浮體在波浪作用下運(yùn)動(dòng)的三維時(shí)域計(jì)算模型

張俊生，滕 斌，叢培文

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

港口中系泊船在波浪作用下運(yùn)動(dòng)問題的本質(zhì)是淺水波浪與浮體的相互作用。與深水情況不同，淺水問題應(yīng)當(dāng)考慮水底、水域邊界的影響及淺水波浪自身的特性，單一模型很難實(shí)現(xiàn)該模擬過程。為此，建立了Boussinesq方程計(jì)算入射波和Laplace方程計(jì)算散射波的全時(shí)域組合計(jì)算模型。有限元法求解的Boussinesq方程能使入射波充分考慮到水底、水域邊界的影響和淺水波浪的特性；散射波被線性化，采用邊界元法求解，并以浮體運(yùn)動(dòng)時(shí)的物面條件為入射波和散射波求解的匹配條件。該方法為完全的時(shí)域方法，計(jì)算網(wǎng)格不隨時(shí)間變動(dòng)，計(jì)算過程較為方便。通過與實(shí)驗(yàn)及其他數(shù)值方法的結(jié)果進(jìn)行比較，驗(yàn)證了本模型對(duì)非線性波面、浮體的運(yùn)動(dòng)都有比較理想的計(jì)算結(jié)果，顯示了本模型對(duì)非線性問題具有較好的計(jì)算能力。

淺水中浮體運(yùn)動(dòng)；非線性波浪作用；Boussinesq方程；Laplace方程；三維時(shí)域數(shù)值模型

Abstract:The calculation of the wave-induced motion of a moored ship in a harbor is essentially a problem of calculating the wave-induced motion of a floating body in shallow water.In this case,the influence of the sea bed and boundaries as well as waves must be considered,as it is different from the problem in deep water.It is difficult to solve the problem by using a single model.Therefore,a combination of an improved Boussinesq equation,to calculate incident waves,and the Laplace equation,to calculate scattered waves,is used to establish a numerical model.The Boussinesq equation solved by the finite element method can embody the characteristics in shallow water while simulating incident waves,and scattered waves are solved by the boundary element method under the linear approximation.The two parts are matched through the boundary condition on the surface of the floating body.The present model is a full time-domain one,but computational grids do not change with the time,which makes the calculation convenient.While calculating nonlinear wave elevations and the motion of the floating body induced by nonlinear waves,the comparison of present results with experimental data and results of other numerical methods shows a good agreement and proves the combined model feasible to deal with nonlinear problems.The model is worthy of being applied to calculate the motion of the moored ship in engineering.

Keywords:motion of the floating body in shallow water; action of nonlinear waves; Boussinesq equation; Laplace equation; 3D time-domain numerical model

計(jì)算港口內(nèi)系泊船在波浪作用下的運(yùn)動(dòng)是一個(gè)傳統(tǒng)而遠(yuǎn)未解決的問題。準(zhǔn)確計(jì)算需要考慮多方面問題：波浪的非線性、港口邊界和地形、浮體的運(yùn)動(dòng)響應(yīng)、系泊系統(tǒng)等，整個(gè)計(jì)算是一個(gè)復(fù)雜的耦合系統(tǒng)。其中，核心問題就是淺水中波浪作用下浮體運(yùn)動(dòng)響應(yīng)的計(jì)算。最早比較完整地給出該問題計(jì)算的是Sawaragi & Kubo[1]。其計(jì)算是一種頻域解析的方法，整個(gè)計(jì)算域被分為外域、船體與港池邊界間的內(nèi)域、船底下方與水底間計(jì)算域等三個(gè)部分，各部分之間通過速度勢(shì)和壓力連續(xù)相匹配，進(jìn)行求解。計(jì)算只能假定水底為平地，港池也只能為規(guī)則的矩形，而計(jì)算條件完全采用線性假設(shè)。Weiler等[2]提出了一個(gè)計(jì)算岸壁前駐波對(duì)船體作用的頻域計(jì)算方法，計(jì)算中，駐波被分離為兩個(gè)反向的行進(jìn)波進(jìn)行考慮。鑒于頻域方法對(duì)入射波特性、計(jì)算域形狀等的限制，為使求解更能適用于真實(shí)、復(fù)雜的近岸港口水域，于是，一些非線性時(shí)域數(shù)值方法被不斷地提出。

最有代表性的是Bingham[3]提出的一種組合型的計(jì)算方法，其通過差分法求解Boussinesq方程(以下簡(jiǎn)稱“B方程”)計(jì)算入射波，通過源匯法在頻域內(nèi)求解線性的散射波激振力、附加質(zhì)量和輻射阻尼，再通過傅里葉變換轉(zhuǎn)化為時(shí)域的散射波作用力、附加質(zhì)量和遲滯函數(shù)，最后把入射波和散射波的激振力相加，進(jìn)而求解含有卷積積分的時(shí)域運(yùn)動(dòng)方程。B方程的特點(diǎn)是可以將三維問題簡(jiǎn)化為二維計(jì)算，使計(jì)算量大為減少，從而能夠在大面積水底地形變化的淺水水域中模擬非線性波浪。因此，在系泊船運(yùn)動(dòng)的計(jì)算模型中應(yīng)用B方程在此后的研究中被不斷地采用，以使計(jì)算模型能在一定程度上考慮水域地形、邊界的影響和波浪的非線性。Van Der Molen[4]也提出了類似的計(jì)算模型，但散射勢(shì)的求解采用了卷積積分求解時(shí)域線性解的方法。不過，歸根到底，該方法還是一種利用脈沖響應(yīng)函數(shù)的算法，不是完全意義上的時(shí)域方法。王大國(guó)等[5]、Wang等[6]及王海龍等[7]給出了幾個(gè)耦合時(shí)域模型，盡管形式不同，但本質(zhì)上皆是內(nèi)外域耦合計(jì)算的方法，外域采用B方程求解，內(nèi)域采用Euler方程或Laplace方程(以下簡(jiǎn)稱“L方程”)計(jì)算。該計(jì)算思路雖是非線性耦合，也是完全時(shí)域的，但最大的不足是這些模型只能適用于物體不動(dòng)的情況，即模型不能考慮物體的運(yùn)動(dòng)，只能計(jì)算波面。故在其進(jìn)行的實(shí)驗(yàn)中，浮體也固定不動(dòng)，從而僅驗(yàn)證了波面計(jì)算的準(zhǔn)確性，但該模型無法滿足實(shí)際工程的需求。

這里也采用“組合”計(jì)算的方法，以B方程模擬入射波浪，以L方程計(jì)算散射波，通過運(yùn)動(dòng)浮體的物面條件連接兩部分的計(jì)算并求解運(yùn)動(dòng)方程。與Bingham[3]和Van Der Molen[4]的模型相比，在浮體運(yùn)動(dòng)響應(yīng)和散射波的計(jì)算上，采用沿時(shí)間步遞進(jìn)求解的方法，從而構(gòu)建一個(gè)嚴(yán)格意義上的完全時(shí)域計(jì)算模型。在數(shù)值方法上，則采用有限元法求解B方程，以滿足曲邊界的計(jì)算，避免差分法處理邊界不靈活的弊端；采用邊界元法計(jì)算浮體的運(yùn)動(dòng)響應(yīng)和散射波，從而既能滿足浮體形狀多樣性的需求，又能簡(jiǎn)化三維計(jì)算。本文對(duì)散射波也做線性化處理。正如Van Der Molen[4]指出的：入射波非線性的影響是重要的，但船在港口內(nèi)受波浪作用的運(yùn)動(dòng)幅度比較小，輻射波和繞射波比入射波小很多，因此采用線性方法計(jì)算波浪對(duì)船體的作用是有效的。也正因此，計(jì)算散射波時(shí)，地形的影響可以忽略，視水底為平地。

通過與Wang等[6]實(shí)驗(yàn)結(jié)果的對(duì)比，驗(yàn)證了本模型對(duì)非線性波面的計(jì)算能力；通過與Bai等[8]的浮體運(yùn)動(dòng)時(shí)域計(jì)算模型的線性波作用結(jié)果和Teng等[9]的二階Stokes波作用下浮體運(yùn)動(dòng)頻域計(jì)算模型的結(jié)果相比，驗(yàn)證了本模型的組合計(jì)算方法和銜接條件的正確性，考察了本模型對(duì)非線性波浪作用下的浮體運(yùn)動(dòng)的計(jì)算能力。

1 模型的建立

1.1入射波的計(jì)算

入射波浪的計(jì)算無需考慮浮體的存在。鑒于等參單元有限元法無法計(jì)算四階及其以上導(dǎo)數(shù)，故采用Beji等[10]的改進(jìn)型B方程計(jì)算入射波浪，方程形式為：

空間離散可以采用四邊形或三角形單元；入射邊界條件則為Fenton[12]的高精度非線性波解析解；固邊壁上的邊界條件按照Engelman等[13]方法處理，可對(duì)曲邊界準(zhǔn)確求解；根據(jù)計(jì)算的具體情況和需要，阻尼層和松弛層可分別布置在計(jì)算域出口和入射邊界前，以吸收計(jì)算域內(nèi)傳出的波浪。時(shí)間積分方面，前四步，采用Runge-Kutta法求解，此后，采用四階Adams-Bashforth-Moulton預(yù)報(bào)-校正方法進(jìn)行計(jì)算。

1.2散射波的計(jì)算

采用L方程求解散射波，其形式為

其中，ηs為散射波波高。物面取浮體平均濕表面，物面條件為

式中：U為浮體上某一點(diǎn)運(yùn)動(dòng)速度矢量，n=(n1,n2,n3)為該點(diǎn)處指向物體內(nèi)的物面單位法向向量，φI為入射勢(shì)。浮體被看作剛體，其位移和轉(zhuǎn)角矢量分別為ξ=(ξ1,ξ2,ξ3)、χ=(χ1,χ2,χ3)，在小轉(zhuǎn)角假定下有

其中，x(x,y,z)為所求點(diǎn)坐標(biāo)，xr(xr,yr,zr)為轉(zhuǎn)動(dòng)中心坐標(biāo)。將式(7)代入式(6)，并進(jìn)一步轉(zhuǎn)化，可得

令(ξ4,ξ5,ξ6)=χ，(n4,n5,n6)=(x-xr)n，式(8)可簡(jiǎn)寫為

其中，N=(n1,n2,n3,n4,n5,n6)為指向物體內(nèi)的廣義物面單位法向向量，ξ也擴(kuò)為6維向量(ξ1,ξ2,ξ3,ξ4,ξ5,ξ6)。φI/n=φIn=uIn是散射波和入射波相銜接的條件。uI為入射波的質(zhì)點(diǎn)速度，可由求得，其水平分量uI、vI近似到二階的沿水深分布表達(dá)式為

對(duì)應(yīng)的垂向分量wI的表達(dá)式為

wI=-

以邊界元法求解散射波時(shí)，水底被看作不透水的平地。采用Rankine源及其關(guān)于海底的像作為格林函數(shù)

計(jì)算中，網(wǎng)格固定不變。通過式(4)、(5)可求得每一時(shí)間步自由水面上的ηs、φs；通過式(9)，可求得物面上的φs/n。通過這些邊界條件，利用邊界元法，可求得自由水面上的φs/n(線性情況下，為z=0面上的φs/z)和物面上的φs，前者結(jié)合ηs通過式(4)、(5)可進(jìn)一步求得下一時(shí)間步自由水面上的ηs、φs，后者可用于計(jì)算浮體受力，利用運(yùn)動(dòng)方程求得下一時(shí)間步的浮體運(yùn)動(dòng)位移，再通過式(9)計(jì)算物面上的φs/n，以此往復(fù)，在時(shí)間域上持續(xù)求解。對(duì)于計(jì)算散射波時(shí)的時(shí)間積分，這里采用四階Runge-Kutta法。

1.3浮體運(yùn)動(dòng)的計(jì)算

對(duì)于波浪與浮體作用的問題，浮體運(yùn)動(dòng)以及其產(chǎn)生的輻射波必須予以考慮。不同于Bingham和Van Der Molen在計(jì)算中使用含有遲滯函數(shù)和附加質(zhì)量的運(yùn)動(dòng)方程，本模型在時(shí)域內(nèi)直接計(jì)算輻射波及其作用力，浮體的運(yùn)動(dòng)方程采用下式求解：

其中，M、B、C分別為浮體質(zhì)量陣、黏性阻尼系數(shù)陣、恢復(fù)力矩陣，均為66的方陣，M和C需要根據(jù)物體質(zhì)心和質(zhì)量分布計(jì)算得出，具體求解過程參見李玉成和滕斌[14]給出的方法；F(t)為六個(gè)分量的廣義波浪激振力，包括入射波和散射波的共同作用；G(t)為系泊系統(tǒng)等其它外部作用力和力矩。該方程包含了6個(gè)廣義方向的運(yùn)動(dòng)：縱蕩(ξ1沿x軸)、橫蕩(ξ2沿y軸)、垂蕩(ξ3沿z軸)、橫搖(ξ4繞x軸)、縱搖(ξ5繞y軸)、回轉(zhuǎn)(ξ6繞z軸)。系泊系統(tǒng)等外部作用力G(t)，如纜繩、護(hù)舷等，在具體工程估算中，可以由計(jì)算中求出的纜繩、護(hù)舷變形結(jié)合產(chǎn)品力學(xué)性能曲線獲得；理論計(jì)算中，可由纜繩、護(hù)舷理論模型求得。本文不對(duì)G(t)進(jìn)行具體分析求解，只集中考慮波浪激振力F(t)對(duì)浮體的作用，計(jì)算后文算例時(shí)，僅在式(13)中加入剛度陣K(式左側(cè)加入Kξ(t)項(xiàng))，代替系泊系統(tǒng)，以保證浮體能圍繞某一平衡位置運(yùn)動(dòng)。

通過流體動(dòng)壓強(qiáng)p在物面的平均濕表面Ωb上積分求解波浪激振力：

同樣，p可分為入射波和散射波產(chǎn)生的動(dòng)壓強(qiáng)pI、ps。根據(jù)B方程理論，pI沿水深分布的二階表達(dá)式為

ps的線性結(jié)果為

運(yùn)動(dòng)方程的時(shí)間積分也采用Runge-Kutta法，與散射波時(shí)間積分同步求解。如果把式(13)抽象表述為

則浮體運(yùn)動(dòng)位移和速度的時(shí)間積分表達(dá)式為

2 模型驗(yàn)證

2.1波面計(jì)算的驗(yàn)證

Wang等[6]開展了一個(gè)波浪對(duì)箱型船繞射的實(shí)驗(yàn)，實(shí)驗(yàn)中箱型船被固定住，并通過浪高儀測(cè)定了一些位置的波高。現(xiàn)利用該實(shí)驗(yàn)結(jié)果檢驗(yàn)本模型對(duì)非線性問題的處理能力。

實(shí)驗(yàn)水池長(zhǎng)31 m、寬14 m，實(shí)驗(yàn)水深為0.3 m，實(shí)驗(yàn)所用箱型船長(zhǎng)l=2.0 m、寬B=0.6 m、高h(yuǎn)=0.45 m、吃水dr=0.24 m。具體船與浪高儀的布置如圖1所示。在靜水面上，以箱形船的幾何中心為原點(diǎn)，建立方向如圖1所示的坐標(biāo)系，表1給出了24個(gè)測(cè)點(diǎn)的坐標(biāo)。實(shí)驗(yàn)中，入射波浪周期為T=3.0 s，波高為H0=0.03 m。

由實(shí)驗(yàn)布置可以看出，水池的上下兩個(gè)邊壁與船模型的距離并不遠(yuǎn)，對(duì)散射波浪場(chǎng)亦會(huì)產(chǎn)生影響，故計(jì)算散射波時(shí)需要對(duì)兩側(cè)邊壁作些處理。如1.2節(jié)所述，在格林函數(shù)中，加入源點(diǎn)關(guān)于一側(cè)邊壁的像，從而去掉該側(cè)邊壁的面積積分，簡(jiǎn)化計(jì)算。另一側(cè)邊壁，則需進(jìn)行網(wǎng)格劃分，參與計(jì)算，其上的邊界條件為φs/n=0。計(jì)算中為保證船不動(dòng)，需取K足夠大，文中取其對(duì)角線元素為1020。

圖1 Wang等實(shí)驗(yàn)布置圖 Fig.1 Layout of Wang et al.’s experiment

表1 測(cè)點(diǎn)位置Tab.1 Positions of the test points

Wang等給出了波浪場(chǎng)穩(wěn)定后，兩個(gè)周期內(nèi)各測(cè)點(diǎn)上波高變化曲線。圖2展示了本文模型計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果的對(duì)比。盡管本模型對(duì)散射波做了線性化處理，致使在某些測(cè)點(diǎn)上存在次波峰的差異，但從波高變化曲線的整體上觀察，本計(jì)算與實(shí)驗(yàn)結(jié)果相當(dāng)一致。可見，在散射波線性化處理的情況下，本模型對(duì)波浪非線性還是有較強(qiáng)的模擬能力，這從一方面體現(xiàn)了本模型的有效性，也證明了Van Der Molen的說法：散射波相比于入射波很小，可線性化處理。

圖2 各測(cè)點(diǎn)上計(jì)算波高與實(shí)驗(yàn)結(jié)果比較Fig.2 Comparison of wave elevations at the test points

2.2浮體運(yùn)動(dòng)計(jì)算的驗(yàn)證

Bai等[8]和Teng等[9]分別給出了開敞水域中浮體在二階Stokes波作用下運(yùn)動(dòng)的時(shí)域和頻域計(jì)算模型，并對(duì)模型的有效性進(jìn)行了充分驗(yàn)證。本文計(jì)算浮體運(yùn)動(dòng)的結(jié)果將分別與該兩種模型的結(jié)果進(jìn)行比較，以驗(yàn)證本文對(duì)浮體運(yùn)動(dòng)計(jì)算的能力。

圖3 被測(cè)試浮體的示意Fig.3 Sizes of the floating body tested

如圖3所示，取開敞的計(jì)算水域水深為d=0.5 m，漂浮于水面上的箱形船長(zhǎng)l=4.2 m、寬B=0.6 m、高h(yuǎn)=0.5 m、吃水dr=0.3 m，計(jì)算中取箱形船為均質(zhì)，質(zhì)心oc為其形心。以質(zhì)心垂直投影到靜水面所得之點(diǎn)為原點(diǎn)o，建立坐標(biāo)系，x軸沿船長(zhǎng)方向，y軸沿船寬方向，z軸沿垂直方向。計(jì)算時(shí)，取原點(diǎn)o為轉(zhuǎn)動(dòng)中心。

測(cè)試波浪為4組，周期皆為T=1.8 s，波高由低到高分別取H0=0.012 m、0.036 m、0.072 m和0.108 m，相應(yīng)的波陡分別為H0/L=0.34%、1%、2%、3%，均沿y方向入射。根據(jù)《The Shore Protection Manual (1984)》，可以大體判斷出四組波浪的非線性和適用理論，如圖4所示。很明顯，Case 1(H0=0.012 m)可以被看作線性波，適用于線性理論求解；Case 2(H0=0.036 m)、Case 3(H0=0.072 m)、Case 4(H0=0.108 m)則皆適用于二階Stokes波理論求解。

為便于計(jì)算結(jié)果的比較，需要保證浮體在波浪激振力作用下于某一平衡位置做簡(jiǎn)諧振動(dòng)，故在計(jì)算中需要給出剛度陣K和黏性阻尼系數(shù)陣B。剛度陣中，取k22=450，k44=24，其余元素取0；黏性阻尼系數(shù)陣中，取b22=470，b44=25，其余元素取0。

圖4 測(cè)試波浪的適用理論Fig.4 Appropriate theories for the tested waves

圖5給出了本模型和Bai等的浮體運(yùn)動(dòng)時(shí)域計(jì)算模型對(duì)Case 1的計(jì)算結(jié)果。縱蕩ξ1、縱搖ξ5、回轉(zhuǎn)ξ6始終為0，故只給出橫蕩ξ2、垂蕩ξ3、橫搖ξ4的歷時(shí)曲線。可以明顯地看到，除因初始狀態(tài)不同而致使前面若干個(gè)周期有差異外，兩者的計(jì)算結(jié)果在浮體運(yùn)動(dòng)穩(wěn)定之后完全一致。這證明了本模型作為一種組合模型，計(jì)算入射波的B方程和計(jì)算散射波的L方程相互間的銜接條件是合理的，匹配過程是準(zhǔn)確的。圖6給出了本模型和Teng等的二階Stokes波作用下浮體運(yùn)動(dòng)頻域計(jì)算模型對(duì)Case 2、3、4的計(jì)算結(jié)果，這里只計(jì)算浮體的周期運(yùn)動(dòng)，未考慮浮體因平均漂移力而產(chǎn)生的定常位移。同樣只給出ξ2、ξ3、ξ4的歷時(shí)曲線，其余方向的運(yùn)動(dòng)始終為0。可以直觀地看到，整體上，本模型與二階Stokes理論的頻域計(jì)算結(jié)果吻合較好，顯示出：盡管對(duì)散射波做了線性化處理，但本模型仍能在一定程度上對(duì)浮體在非線性波浪作用下的運(yùn)動(dòng)進(jìn)行模擬和計(jì)算。進(jìn)一步觀察可以發(fā)現(xiàn)，在波陡從1%增加到3%的過程中，文中計(jì)算的運(yùn)動(dòng)幅值會(huì)逐漸小于二階Stokes理論的頻域結(jié)果。波浪非線性越強(qiáng)，偏小的程度就越略大一些，這是散射波線性化的必然結(jié)果。由圖4可知，Case 4的波浪非線性已經(jīng)比較可觀，但這里計(jì)算結(jié)果與二階Stokes理論結(jié)果仍比較接近；而另一方面，一般認(rèn)為在港池中波浪的非線性不會(huì)太強(qiáng)，所以，本模型就系泊船的運(yùn)動(dòng)而言，對(duì)非線性波浪的計(jì)算能力是可以接受的。

圖5 線性波作用下浮體運(yùn)動(dòng)的計(jì)算結(jié)果比較(T=1.8 s,H0=0.012 m)Fig.5 Comparison of calculated results for the motion of the floating body induced by a linear wave (T=1.8 s,H0=0.012 m)

圖6 二階Stokes波作用下浮體運(yùn)動(dòng)的計(jì)算結(jié)果比較Fig.6 Comparison of calculated results for the motion of the floating body induced by 2nd order Stokes waves

3 結(jié) 語

港口中波浪作用下系泊船運(yùn)動(dòng)問題的核心是計(jì)算淺水中浮體在波浪作用下的運(yùn)動(dòng)。而相比于深水問題，海底地形、水域邊界的影響和淺水波浪特性成為應(yīng)該考慮的必要因素，因此，建立含有B方程的組合模型成為一種有效的方法。相比于已有的模型，文中建立的B方程和L方程的組合模型是一個(gè)能有效計(jì)算浮體運(yùn)動(dòng)的完全時(shí)域模型。采用有限元法求解B方程模擬入射波，在體現(xiàn)淺水波浪特性和水底影響之外，與差分法相比又能較好地處理曲邊界的影響。散射波的線性化使計(jì)算得到簡(jiǎn)化，浮體的物面積分也只需在平均濕表面上進(jìn)行，在時(shí)域求解過程中網(wǎng)格也始終保持不變，從而使整個(gè)求解過程相對(duì)簡(jiǎn)單、方便。

通過與Wang等實(shí)驗(yàn)結(jié)果的對(duì)比，體現(xiàn)了本模型對(duì)非線性波浪的計(jì)算能力；在計(jì)算線性波作用下浮體的運(yùn)動(dòng)時(shí)，本文結(jié)果與Bai等的時(shí)域模型結(jié)果完全一致，說明了本模型作為組合模型，各計(jì)算部分的匹配和計(jì)算過程是準(zhǔn)確的。在計(jì)算二階Stokes波作用下浮體的運(yùn)動(dòng)時(shí)，本文結(jié)果與Teng等的二階Stokes理論頻域模型的結(jié)果也基本一致，表現(xiàn)出本模型對(duì)非線性問題有較理想的計(jì)算能力。

綜上可見，本模型是一個(gè)計(jì)算過程較簡(jiǎn)單、能在較大程度上計(jì)算出淺水中浮體在非線性波浪作用下運(yùn)動(dòng)響應(yīng)的組合模型。可以預(yù)期，該模型在不斷完善之后，能夠成為工程中估算系泊船運(yùn)動(dòng)的有效工具。

[1] SAWARAGI T,KUBO M.The motion of a ship in a harbor basin[C]//Proc.18th ICCE.1982,3：2743-2762.

[2] WEILER O M,DEKKER J.Mooring container ships exposed to long waves[C]//Proc.13th Int.Harbor Congress.2003.

[3] BINGHAM H B.A hybrid Boussinesq-panel method for predicting the motion of a moored ship[J].Coastal Engineering,2000,40：21-38.

[4] VAN DER MOLEN W.Behavior of moored ships in harbors[D].Delft：The Water Research Center of The Delft University of Technology,2006：55-68.

[5] 王大國(guó)，鄒志利，唐春安.波浪對(duì)箱形船作用的三維耦合計(jì)算模型[J].船舶力學(xué)，2007，11(4)：533-544.(WANG Daguo,ZOU Zhili,TANG Chunan.Time stepping solutions of nonlinear wave forces on a three-dimensional box-shaped ship in a harbor[J].Journal of Ship Mechanics,2007,11(4)：533-544.(in Chinese))

[6] WANG D G,ZOU Z L,TANG C A.A three-dimensional coupled numerical model of nonlinear waves in a harbor[J].Science in China Series E,Technological Sciences,2008,51(12)：2195-2206.

[7] 王海龍，鄒志利，周亞龍，等.波浪對(duì)圓柱作用的三維耦合計(jì)算模型[J].中國(guó)海洋平臺(tái)，2010，25(5)：38-48.(WANG Hailong,ZOU Zhili,ZHOU Yalong,et al.The three-dimensional coupled numerical model of non-linear wave acting on a circular cylinder[J].China Offshore Platform,2010,25(5)：38-48.(in Chinese))

[8] BAI W,TENG B.Second-order wave diffraction around 3-D bodies by a time-domain method[J].China Ocean Engineering,2001,15(1)：73-85.

[9] TENG B,EATOCK R.A BEM method for the second-order wave action on a floating body in monochromatic waves[C]//Proc.1st Int.Conf.on Hydrodynamics.1994：274-281.

[10] BEJI S,NADAOKA K.A formal derivation and numerical modeling of the improved Boussinesq equations for varying depth[J].Ocean Engineering,1996,23(8)：691-704.

[11] LI Y S,LIU S X,YU Y X,et al.Numerical modeling of Boussinesq equations by finite element method[J].Coastal Engineering,1999,37：97-122.

[12] FENTON J D.The numerical solution of steady water wave problems[J].Computers & Geosciences,1988,14：357-368.

[13] ENGELMAN M S,SANI R L.The implementation of normal and/or tangential boundary conditions in finite element codes for incompressible fluid flow[J].International Journal for Numerical Methods in Fluids,1982,2：225-238.

[14] 李玉成，滕斌.波浪對(duì)海上建筑物的作用(第二版)[M].北京：海洋出版社，2002：70-73.(LI Y C,TENG B.Wave action on maritime structures (2nd edition)[M].Beijing:China Ocean Press,2002：70-73.(in Chinese))

A 3D time-domain numerical model for floating body motion
induced by shallow water waves

ZHANG Junsheng,TENG Bin,CONG Peiwen

(State Key Laboratory of Coastal and Offshore Engineering,Dalian University of Technology,Dalian 116024,China)

TV139.2

A

10.16483/j.issn.1005-9865.2016.01.001

1005-9865(2016)01-0001-09

2014-11-23

國(guó)家自然科學(xué)基金面上項(xiàng)目(51379032)；國(guó)家自然科學(xué)基金重大項(xiàng)目(51490672)

張俊生(1980-)，男，遼寧大連人，助理研究員，博士生，從事波浪對(duì)海上結(jié)構(gòu)物作用的研究。E-mail：yzjz009@126.com

滕 斌。E-mail:bteng@dlut.edu.cn