李裕龍，朱仁傳，繆國(guó)平，范 菊

（上海交通大學(xué) 船舶海洋與建筑工程學(xué)院，海洋工程國(guó)家重點(diǎn)試驗(yàn)室，上海 200240）

基于全時(shí)域勢(shì)流理論的船舶與液艙晃蕩耦合運(yùn)動(dòng)的數(shù)值計(jì)算

李裕龍，朱仁傳，繆國(guó)平，范 菊

（上海交通大學(xué) 船舶海洋與建筑工程學(xué)院，海洋工程國(guó)家重點(diǎn)試驗(yàn)室，上海 200240）

基于三維線(xiàn)性有航速時(shí)域勢(shì)流理論計(jì)算船體時(shí)域運(yùn)動(dòng)外域波浪力，同時(shí)采用三維全非線(xiàn)性時(shí)域勢(shì)流理論來(lái)計(jì)算艙內(nèi)液體的非線(xiàn)性晃蕩所誘導(dǎo)力與力矩，進(jìn)而建立了波浪中載液船舶耦合運(yùn)動(dòng)方程。該方法能夠完整地考慮波浪、船體和液艙晃蕩之間的實(shí)時(shí)耦合作用。研究結(jié)果表明：通過(guò)模型實(shí)驗(yàn)和數(shù)值模擬計(jì)算的對(duì)比，數(shù)值模擬計(jì)算能夠清晰顯現(xiàn)出液艙晃蕩對(duì)船體全局運(yùn)動(dòng)影響，無(wú)航速船體運(yùn)動(dòng)RAO與模型實(shí)驗(yàn)結(jié)果吻合良好，有航速運(yùn)動(dòng)計(jì)算合乎預(yù)期。

液艙晃蕩；耦合運(yùn)動(dòng)；非線(xiàn)性時(shí)域；有航速理論

0 引 言

船舶在波浪上的運(yùn)動(dòng)激勵(lì)液艙內(nèi)的液體產(chǎn)生晃蕩，液艙晃蕩誘導(dǎo)的沖擊力同時(shí)作用在艙壁上，進(jìn)而影響到船舶的運(yùn)動(dòng)姿態(tài)。船舶運(yùn)動(dòng)與液艙晃蕩之間的耦合效應(yīng)研究一直是船舶與海洋工程界一直關(guān)注的問(wèn)題，合適且高效的理論與數(shù)值計(jì)算方法對(duì)探討分析船舶運(yùn)動(dòng)與液艙晃蕩之間的耦合效應(yīng)是十分重要的。目前國(guó)內(nèi)外許多學(xué)者已對(duì)船舶運(yùn)動(dòng)與液艙晃蕩之間耦合效應(yīng)進(jìn)行了一些分析與研究工作[1-3]，綜合來(lái)說(shuō)，基于勢(shì)流理論的耦合效應(yīng)的相關(guān)研究主要可分為兩種途徑：相關(guān)領(lǐng)域的模型試驗(yàn)研究[4]，以及用來(lái)輔助實(shí)驗(yàn)研究且專(zhuān)門(mén)針對(duì)本問(wèn)題的數(shù)值計(jì)算工作[5-6]。

船舶在波浪上的運(yùn)動(dòng)計(jì)算是船舶與海洋工程領(lǐng)域的經(jīng)典問(wèn)題。如果考慮到晃蕩流體的非線(xiàn)性作用，則需要在時(shí)域下來(lái)分析船舶運(yùn)動(dòng)。時(shí)域三維Rankine源邊界元方法相對(duì)于復(fù)雜時(shí)域格林函數(shù)方法，雖然其需要加以粘性處理輻射面，或者是借助于匹配面的方法處理艦船時(shí)域遠(yuǎn)場(chǎng)波浪，但其編程相對(duì)簡(jiǎn)單以及其對(duì)船型幾何外形適應(yīng)性更強(qiáng)的優(yōu)點(diǎn)仍然在學(xué)術(shù)界引起足夠重視。本文使用時(shí)域三維Rankine源邊界元方法以及脈沖響應(yīng)理論模擬時(shí)域艦船有航速運(yùn)動(dòng)。

液艙流體晃蕩問(wèn)題在學(xué)術(shù)界與工程界中同樣是經(jīng)典的水動(dòng)力學(xué)問(wèn)題。自從液化天然氣船等載液高附加值船舶的迅速發(fā)展，液艙流體晃蕩業(yè)已成為水動(dòng)力學(xué)研究工作的熱點(diǎn)問(wèn)題。已有很多學(xué)者分別采用了頻域勢(shì)流理論、時(shí)域非線(xiàn)性勢(shì)流理論、粘性流理論方法來(lái)研究模擬液艙晃蕩問(wèn)題及其誘導(dǎo)的水動(dòng)力砰擊載荷。本文采用時(shí)域非線(xiàn)性勢(shì)流理論，藉此配合艦船時(shí)域線(xiàn)性有航速勢(shì)流理論，構(gòu)建了完全時(shí)域下的艦船耦合液艙晃蕩時(shí)域運(yùn)動(dòng)計(jì)算方法，為基于勢(shì)流理論研究艦船耦合液艙晃蕩現(xiàn)象提供了新的途徑。

本文將船舶在波浪上的運(yùn)動(dòng)問(wèn)題和液艙內(nèi)液體晃蕩問(wèn)題采用不同的時(shí)域勢(shì)流理論求解，船舶的運(yùn)動(dòng)采用基于脈沖響應(yīng)函數(shù)理論的三維線(xiàn)性有航速時(shí)域勢(shì)流理論，計(jì)算首先需要求得有航速水動(dòng)力系數(shù)以及脈沖響應(yīng)函數(shù)，以及時(shí)域下的繞射問(wèn)題，進(jìn)而求解時(shí)域有航速艦船時(shí)域運(yùn)動(dòng)方程。本文同時(shí)應(yīng)用三維非線(xiàn)性時(shí)域勢(shì)流理論模擬非線(xiàn)性液艙晃蕩問(wèn)題，將實(shí)時(shí)誘導(dǎo)力與力矩在同一時(shí)刻添加入船舶有航速時(shí)域運(yùn)動(dòng)方程，從而建立了載液船舶在波浪上運(yùn)動(dòng)的有航速時(shí)域運(yùn)動(dòng)方程。這種方法在零航速時(shí)刻能夠模擬常見(jiàn)的無(wú)航速計(jì)算工況與載液浮式結(jié)構(gòu)物的運(yùn)動(dòng)響應(yīng)，同時(shí)能夠進(jìn)行有航速的載液艦船時(shí)域運(yùn)動(dòng)模擬，具有很好的適用范圍。通過(guò)對(duì)加載方形液艙的S175船模進(jìn)行耦合運(yùn)動(dòng)的數(shù)值模擬及與實(shí)驗(yàn)結(jié)果的對(duì)比，以及有航速艦船的數(shù)值計(jì)算結(jié)果，驗(yàn)證了本文數(shù)值計(jì)算方法的準(zhǔn)確性，表明使用本文提出的方法進(jìn)行載液有航速艦船時(shí)域運(yùn)動(dòng)是完全可行的?；诖吮疚耐ㄟ^(guò)對(duì)一艘載液15000GT船橫搖運(yùn)動(dòng)進(jìn)行了計(jì)算分析，更進(jìn)一步的拓展了計(jì)算代碼的實(shí)用性與可靠性。

1 載液船舶耦合運(yùn)動(dòng)相關(guān)基本理論

圖1 船體運(yùn)動(dòng)的縱坐標(biāo)系示意圖Fig.1 The coordinate system of ship motion

1.1 載液船舶在波浪上的運(yùn)動(dòng)

本文采用勢(shì)流理論求解船舶在波浪上的運(yùn)動(dòng)。這里給出加載液艙的船舶在波浪上運(yùn)動(dòng)滿(mǎn)足的數(shù)學(xué)模型。線(xiàn)性入射波速度勢(shì)有以下形式：

式中：ζI，ω，k分別為入射波的波幅，自然頻率與波數(shù)，ωe為遭遇頻率。定義隨船體運(yùn)動(dòng)的笛卡爾坐標(biāo)系坐標(biāo)系示意圖見(jiàn)圖1。在規(guī)則波激勵(lì)下船舶保持固定的航速作六自由度運(yùn)動(dòng)。

本文采用三維有航速時(shí)域勢(shì)流理論求解外域流場(chǎng)的速度勢(shì)函數(shù)。船舶運(yùn)動(dòng)輻射問(wèn)題速度勢(shì)所滿(mǎn)足的流場(chǎng)邊界條件如下：

根據(jù)時(shí)域分析的脈沖響應(yīng)方法，能夠?qū)⑤椛鋭?shì)進(jìn)一步為瞬時(shí)項(xiàng)和記憶項(xiàng)。在線(xiàn)性系統(tǒng)下，時(shí)域輻射系統(tǒng)的輸入為浮體的k模態(tài)的運(yùn)動(dòng)速度。假定浮體的初始位移則輻射問(wèn)題的物面條件能夠有如下的表達(dá)式[7]：

輻射勢(shì)的脈沖響應(yīng)函數(shù)根據(jù)物面條件分解如下：

進(jìn)而便可獲得艦船輻射運(yùn)動(dòng)的輻射力：

其中：附加質(zhì)量μjk取決于船體的外表面幾何形狀，bjk和cjk取決于船的幾何形狀和船速。較為復(fù)雜的記憶效應(yīng)體現(xiàn)在時(shí)延函數(shù)中，其取決于船體外表面的的幾何外形與航行速度。

同樣，將波浪的運(yùn)動(dòng)視為簡(jiǎn)諧運(yùn)動(dòng)，且將入射波波高ζI作為繞射系統(tǒng)的輸入，同樣使用脈沖響應(yīng)理論求解繞射問(wèn)題，便可求得艦船運(yùn)動(dòng)的外部波浪力。

1.2 基于三維全非線(xiàn)性勢(shì)流理論的液艙晃蕩問(wèn)題

本文選取動(dòng)坐標(biāo)系作為計(jì)算坐標(biāo)系，其優(yōu)勢(shì)在于壁面相對(duì)于動(dòng)坐標(biāo)系保持靜止，這樣便給勢(shì)流計(jì)算帶來(lái)不少便利。液艙動(dòng)坐標(biāo)系如圖2所示。

圖2 所選液艙坐標(biāo)系示意圖Fig.2 The coordinate system of the tank

速度勢(shì)函數(shù)滿(mǎn)足的控制即拉普拉斯方程如下：

在流域內(nèi)

自由表面上動(dòng)力學(xué)邊界條件：

自由表面上運(yùn)動(dòng)學(xué)邊界條件：

艙壁上的邊界條件：

為了要計(jì)算實(shí)時(shí)自由面的變化，自然必須計(jì)及自由液面處的速度勢(shì)和波高隨時(shí)間的變化。不同時(shí)刻的自由面升高與速度勢(shì)，需要采用時(shí)間的差分方式來(lái)計(jì)算。得到邊界上每點(diǎn)的速度勢(shì)后，根據(jù)非線(xiàn)性拉格朗日方程：

便能夠方便計(jì)算出邊界處該時(shí)刻的壓力。通過(guò)將壓力沿著液艙所有壁面的積分就能夠計(jì)算當(dāng)前時(shí)刻液艙晃蕩對(duì)壁面誘導(dǎo)的晃蕩力與力矩。在計(jì)算中同時(shí)需要實(shí)時(shí)對(duì)自由液面進(jìn)行光順化處理。

2 載液船舶在波浪上時(shí)域運(yùn)動(dòng)的數(shù)值算法

在液艙晃蕩問(wèn)題的非定常計(jì)算求解過(guò)程中，在某一計(jì)算時(shí)刻通過(guò)對(duì)液艙壁面的壓力積分可以得到當(dāng)前時(shí)刻下液艙晃蕩在艙壁上誘導(dǎo)的水動(dòng)力幅值?；问幷T導(dǎo)的水動(dòng)力進(jìn)而會(huì)影響到船舶的運(yùn)動(dòng)姿態(tài)，因此載液船舶時(shí)域運(yùn)動(dòng)方程需要在當(dāng)前時(shí)間步下建立。耦合液艙晃蕩作用下的時(shí)域船舶全局運(yùn)動(dòng)方程有如下形式：

圖3 S175船型線(xiàn)Fig.3 Main parameters of LNG ship

3 數(shù)值計(jì)算

3.1 S175船模主尺度與有航速時(shí)域計(jì)算

數(shù)值計(jì)算工作采用的是高速集裝箱船型S175，其船型線(xiàn)見(jiàn)圖3所示。其無(wú)航速船舶耦合液艙晃蕩的模型試驗(yàn)是在中國(guó)船舶科學(xué)研究中心耐波性水池中進(jìn)行[8]，試驗(yàn)?zāi)Ｐ褪且粭l加載了一個(gè)方形液艙的S175船模，模型與實(shí)船的縮尺比為1:55，其具體的主尺度參數(shù)見(jiàn)表1。

表1 LNG主尺度Tab.1 Principal dimensions of LNG

所加載的液艙為方形液艙，其尺寸和安裝位置見(jiàn)圖4。液艙長(zhǎng)度為600 mm，寬度為300 mm，高度為250 mm，艙內(nèi)液體的深度為125 mm，液艙重心位置與船模重心位置重合。

圖5和圖6給出了S175船的時(shí)域運(yùn)動(dòng)響應(yīng)幅值的數(shù)值計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果的對(duì)比。數(shù)值結(jié)果與實(shí)驗(yàn)結(jié)果兩者吻合良好。無(wú)論是水動(dòng)力系數(shù)還是幅值響應(yīng)，計(jì)算程序都能夠給出較為令人滿(mǎn)意的計(jì)算結(jié)果。這驗(yàn)證了有航速時(shí)域計(jì)算的有效性。

圖4 加載液艙的尺寸與位置Fig.4 Size and location of tank

圖5 迎浪工況下數(shù)值與實(shí)驗(yàn)所得的S175船模垂蕩RAOFig.5 Comparison of heave RAO of ship in head sea by experiment and calculation

圖6 迎浪工況下數(shù)值與實(shí)驗(yàn)所得的S175船模縱搖RAOFig.6 Comparison of Pitch RAO of ship in head sea by experiment and calculation

3.2 無(wú)航速耦合液艙晃蕩船舶時(shí)域運(yùn)動(dòng)的數(shù)值計(jì)算

本節(jié)給出S175船舶在無(wú)航速情況下的計(jì)算結(jié)果與模型實(shí)驗(yàn)結(jié)果作為耦合計(jì)算的驗(yàn)證。這里將有航速耦合計(jì)算程序中的航速設(shè)為零。為了驗(yàn)證程序的有效性，這里給出了S175船的無(wú)航速的RAO曲線(xiàn)與實(shí)驗(yàn)結(jié)果的對(duì)比?？梢灶A(yù)計(jì)，在無(wú)航速下的結(jié)果如果能夠符合實(shí)驗(yàn)結(jié)果，便可以首先驗(yàn)證時(shí)域運(yùn)動(dòng)方程與液艙晃蕩問(wèn)題計(jì)算無(wú)誤。這里就加載了方形液艙的S175船模，在迎浪和橫浪兩個(gè)工況下進(jìn)行了無(wú)航速耦合運(yùn)動(dòng)的時(shí)域數(shù)值計(jì)算。

圖7 迎浪工況下數(shù)值與實(shí)驗(yàn)所得的S175船模縱搖RAOFig.7 Comparison of pitch RAO of ship in head sea by experiment and calculation

圖8 橫浪工況下數(shù)值與實(shí)驗(yàn)所得的S175船模橫搖RAOFig.8 Comparison of roll RAO of ship in head sea by experiment and calculation

從圖中可以看出數(shù)值計(jì)算與模型實(shí)驗(yàn)的結(jié)果同樣吻合較好，在此階段驗(yàn)證了耦合計(jì)算程序中液艙晃蕩模塊與耦合計(jì)算模塊的有效性，并且經(jīng)由有航速船舶的計(jì)算結(jié)果的準(zhǔn)確性，便能夠預(yù)測(cè)出本文提出的有航速船舶耦合液艙晃蕩時(shí)域計(jì)算方法的可靠與有效性。

3.3 有航速耦合液艙晃蕩船舶時(shí)域運(yùn)動(dòng)的數(shù)值計(jì)算

在上一節(jié)無(wú)航速計(jì)算結(jié)果準(zhǔn)確的基礎(chǔ)上，以及有航速時(shí)域數(shù)值計(jì)算能夠保證計(jì)算精度的前提下，便可以方便地進(jìn)行有航速船舶耦合液艙晃蕩的時(shí)域數(shù)值計(jì)算工作。圖9至圖12給出了不同航速下船舶耦合液艙晃蕩時(shí)域運(yùn)動(dòng)的縱搖以及橫搖模態(tài)的響應(yīng)曲線(xiàn)。

圖9 迎浪工況下加載與不加載液艙所得的S175船?？v搖RAO（Fn=0.10）Fig.9 Comparison of pitch RAO of ship in head sea of with tank and without tank(Fn=0.10)

圖10 橫浪工況下加載與不加載液艙所得的S175船模橫搖RAO（Fn=0.10）Fig.10 Comparison of roll RAO of ship in head sea of with tank and without tank(Fn=0.10)

在圖9至圖12中給出了加載液艙有航速S175船橫搖與縱搖的運(yùn)動(dòng)響應(yīng)曲線(xiàn)，其中為了對(duì)比同時(shí)給出了不加載液艙的S175船的時(shí)域運(yùn)動(dòng)曲線(xiàn)。在RAO中的橫坐標(biāo)的ω為波浪自然頻率。從計(jì)算結(jié)果可以看出，液艙晃蕩對(duì)于不同航速下船舶的縱向運(yùn)動(dòng)影響不是很大，并且在計(jì)算中觀(guān)察到垂蕩運(yùn)動(dòng)所受液艙晃蕩的影響也并不明顯，但是隨著船舶的縱向艙室加多或者幾何形狀的改變對(duì)縱搖的影響需要在實(shí)際計(jì)算中加以確定。而在船舶的有航速橫浪運(yùn)動(dòng)時(shí)，其運(yùn)動(dòng)規(guī)律與無(wú)航速時(shí)趨勢(shì)大致相同。加載液艙船舶在低于船舶自身的橫搖自振頻率下運(yùn)動(dòng)響應(yīng)高于不加載液艙的運(yùn)動(dòng)幅值，其峰值也大于無(wú)載液船舶橫搖運(yùn)動(dòng)峰值，其原因是此時(shí)液艙所誘導(dǎo)的力矩與船舶的波浪誘導(dǎo)力矩幅值疊加所造成的。并且從上一節(jié)可以看到，本文的無(wú)航速計(jì)算結(jié)果同實(shí)驗(yàn)結(jié)果對(duì)比良好，即這種數(shù)值方法可以直接對(duì)無(wú)航速載液船舶的時(shí)域運(yùn)動(dòng)響應(yīng)進(jìn)行預(yù)報(bào)。

圖11 迎浪工況下加載與不加載液艙所得的S175船?？v搖RAO（Fn=0.30）Fig.11 Comparison of pitch RAO of ship in head sea of with tank and without tank(Fn=0.30)

圖12 橫浪工況下加載與不加載液艙所得的S175船模橫搖RAO（Fn=0.30）Fig.12 Comparison of roll RAO of ship in head sea of with tank and without tank(Fn=0.30)

3.4 載液15000GT集裝箱船橫搖運(yùn)動(dòng)時(shí)域數(shù)值計(jì)算

為了拓展計(jì)算代碼的實(shí)用性，本文選取了一艘具有實(shí)驗(yàn)數(shù)據(jù)的集裝箱船15000GT作為計(jì)算算例[9]，針對(duì)不同加載深度的液艙與15000GT船模的耦合橫搖運(yùn)動(dòng)做了數(shù)值計(jì)算研究。圖13是15000GT的邊界元計(jì)算網(wǎng)格。

模型試驗(yàn)中對(duì)船舶在無(wú)航速規(guī)則波作用下的橫向運(yùn)動(dòng)響應(yīng)做了測(cè)量與分析。液艙尺寸與幾何位置見(jiàn)表2和圖14。模型試驗(yàn)重點(diǎn)在于測(cè)量不同深度裝載的液艙對(duì)橫搖運(yùn)動(dòng)的影響，并且對(duì)于不加載液艙的船模在同樣來(lái)波工況中的運(yùn)動(dòng)響應(yīng)做了測(cè)量。加載不同深度的液艙的模型工況見(jiàn)表3。

圖13 15000GT船邊界元計(jì)算網(wǎng)格Fig.13 The boundary mesh of 15000GT container ship

圖14 液艙位置簡(jiǎn)圖Fig.14 Schematic profile of the rectangle tank

表2 液艙的幾何尺寸與位置Tabl.2 Main dimension of rectangle tank

表3 模型試驗(yàn)工況Tab.3 Experiment parameters

圖15 水深8 cm液艙15000GT集裝箱船橫搖RAOFig.15 Roll RAO of 15000GT container ship model with 8 cm depth liquid tank

圖15-17為在液艙不同裝載深度下的船模橫搖運(yùn)動(dòng)RAO曲線(xiàn)與實(shí)驗(yàn)值的對(duì)比結(jié)果，如圖所示數(shù)值計(jì)算結(jié)果趨勢(shì)與實(shí)驗(yàn)結(jié)果吻合良好，驗(yàn)證了本文所提出的數(shù)值計(jì)算方法對(duì)于不同工況橫搖運(yùn)動(dòng)計(jì)算的可行性，拓展并延伸了計(jì)算代碼的實(shí)用性。

從圖15-17可以看出，不同裝載下的液艙對(duì)橫搖運(yùn)動(dòng)有著很大的影響，且程度隨裝載的深度的增加而增大。首先是幅值，隨著裝載深度的增加，其幅值逐漸減小，其原因是液艙誘導(dǎo)的內(nèi)部力隨著裝載深度的增加而增大，且由于液艙幾何外形所致，在當(dāng)前激勵(lì)頻率下其相位與波浪所激勵(lì)的外部力的相位一直有著180度左右的偏差，其合成作用便是減小船舶受到的外部力，因此船舶的橫搖幅值逐漸減小。并且隨著裝載深度的增加，液艙內(nèi)部的液體的固有頻率也會(huì)逐漸變化，并且與當(dāng)前運(yùn)動(dòng)模態(tài)的自由液面尺度相關(guān)，并且晃蕩隨著激勵(lì)頻率接近共振頻率時(shí)，其液面高度隨時(shí)間的變化會(huì)頗為劇烈，當(dāng)前裝載下所誘導(dǎo)的最大波面升高所包含的非線(xiàn)性因素也會(huì)影響船舶的橫搖共振頻率。并且由于液艙內(nèi)部液體的形狀隨著運(yùn)動(dòng)的幅值會(huì)逐漸變化，從簡(jiǎn)單的角度來(lái)說(shuō)液艙內(nèi)部液體的重心由于液體形狀變化同樣產(chǎn)生偏移，這樣便致使船舶的穩(wěn)心同樣會(huì)發(fā)生變化，因此橫搖運(yùn)動(dòng)自然會(huì)受到不同程度的影響。

圖16 水深10 cm液艙15000GT集裝箱船橫搖RAOFig.16 Roll RAO of 15000GT container ship model with 10 cm depth liquid tank

圖17 水深12 cm液艙15000GT集裝箱船橫搖RAOFig.17 Roll RAO of 15000GT container ship model with 12 cm depth liquid tank

4 結(jié) 論

本文在成功地將船舶在波浪上運(yùn)動(dòng)和液艙內(nèi)液體晃蕩問(wèn)題分析同時(shí)采用時(shí)域勢(shì)流理論求解。時(shí)域船舶運(yùn)動(dòng)問(wèn)題采用基于勢(shì)流理論的脈沖響應(yīng)函數(shù)及Rankine源邊界元法求解，非定常液艙晃蕩問(wèn)題采用非線(xiàn)性時(shí)域勢(shì)流理論和Rankine源邊界元方法求解。本文基于目前的研究結(jié)果，能夠得出以下結(jié)論：

液艙晃蕩效應(yīng)對(duì)載液船舶在波浪上的運(yùn)動(dòng)有著較為重要的影響，特別是對(duì)于橫浪工況下的橫搖運(yùn)動(dòng)。船舶運(yùn)動(dòng)的遭遇頻率在船舶響應(yīng)頻率范圍附近時(shí)，液艙晃蕩誘導(dǎo)的橫搖力矩同波浪誘導(dǎo)的橫搖力矩之間的相位差在一定程度上減小了船舶橫搖運(yùn)動(dòng)的幅值。在較低頻率的入射波激勵(lì)下，液艙晃蕩誘導(dǎo)的橫搖力矩同波浪誘導(dǎo)的橫搖力矩之間的相位差較小，兩種橫搖力矩的疊加會(huì)增大船舶的橫搖運(yùn)動(dòng)幅值。在橫浪工況下，液艙晃蕩對(duì)船舶的垂蕩運(yùn)動(dòng)影響較小。在迎浪工況下，隨著航速提高，本文所取算例來(lái)看液艙晃蕩所誘導(dǎo)的外力與外力矩對(duì)船舶的垂蕩與縱搖等縱向運(yùn)動(dòng)的幅值影響均比較小。

本文采用的數(shù)值計(jì)算方法適用于計(jì)算分析載液船舶在波浪上的時(shí)域運(yùn)動(dòng)，不僅適用于船舶的無(wú)航速時(shí)域運(yùn)動(dòng)，同樣還可以適用和拓展到有航速的船舶時(shí)域運(yùn)動(dòng)。本研究所用的計(jì)算液艙的幾何形狀規(guī)則，而其他各異的幾何形狀的液艙，與不同類(lèi)型的液艙布置，均會(huì)對(duì)載液船舶運(yùn)動(dòng)造成不同的影響，其運(yùn)動(dòng)規(guī)律取決于船舶的幾何外形與航行工況。因此本文所得出的結(jié)論證明了開(kāi)發(fā)的計(jì)算程序的準(zhǔn)確性，但是對(duì)于不同類(lèi)型的載液船舶由于其幾何因素與設(shè)計(jì)理念差異頗大，本文的計(jì)算結(jié)論對(duì)其他各異的載液船型應(yīng)該不具有完全的推廣性，針對(duì)不同的船舶需要更進(jìn)一步的詳細(xì)計(jì)算。

本文方法由于基于脈沖響應(yīng)函數(shù)方法，具有較高的計(jì)算效率，數(shù)值計(jì)算與試驗(yàn)結(jié)果吻合良好，為設(shè)計(jì)載液船舶或減搖水艙等的前期設(shè)計(jì)提供了快速有效的分析方法和技術(shù)手段。這種方法的局限性在于無(wú)法處理自由液面破碎等強(qiáng)非線(xiàn)性現(xiàn)象，但可以較好地預(yù)報(bào)船體與液艙晃蕩耦合運(yùn)動(dòng)的共振頻率區(qū)間與基本幅值。

[1]Kim Kyung sung,Kim Moo Hyun.Sloshing effects on Multi-Vessel motions by using moving particle simulation[C]//Proceedings of the 21th International Offshore and Polar Engineering Conference,June 19-24,2011.Maui,Hawaii,USA, 2011:39-45.

[2]Huang Z J,Danaczko M A,Esenkov O E,Martin C B,et al.Coupled tank sloshing and LNG carrier motions[C]//Proceedings of the 19th International Offshore and Polar Engineering Conference,June 21-26,2009.Osaka,Japan,2009: 92-99.

[3]Wang X,Arai Makoto.A study on coupling effect between sea-keeping and sloshing for membrane-type LNG carrier[C]// Proceedings of the 20th International Offshore and Polar Engineering Conference,June 20-25,2010.Beijing,China, 2010:245-251.

[4]Nam Bo-Woo,Kim Yonghwan,Kim Dae-Woong,Kim Yong-soo.Experimental and numerical studies on ship motion responses coupled with sloshing in waves[J].Journal of Ship Research,2009,53(2):68-82.

[5]Tabri Kristjan,Matusiak Jerzy,Varsta Petri.Sloshing interaction in ship collisions-An experimental and numerical study [J].Ocean Engineering,2009,36:1366-1376.

[6]Kaminski Miroslaw Lech,Bogaert Hannes.Full scale sloshing impact tests[C]//Proceedings of the 19th International Offshore and Polar Engineering Conference,June 21-26,2009.Osaka,Japan,2009:125-134.

[7]戴遺山,段文洋.船舶在波浪中運(yùn)動(dòng)的時(shí)域理論[M].北京:國(guó)防工業(yè)出版社,2008.

[8]鄒 康.耦合液艙晃蕩的船舶運(yùn)動(dòng)性能研究[D].鎮(zhèn)江:江蘇科技大學(xué),2009.

[9]洪 亮,朱仁傳,繆國(guó)平,范 菊.波浪中船體與液艙晃蕩耦合運(yùn)動(dòng)的時(shí)域數(shù)值計(jì)算[J].哈爾濱工程大學(xué)學(xué)報(bào),2012, 33(5):635-641.

Numerical method of ship motions coupled with tank sloshing based on fully time domain potential flow theory

LI Yu-long,ZHU Ren-chuan,MIAO Guo-ping,FAN Ju
(The State Key Laboratory of Ocean Engineering,School of Naval Architecture,Ocean and Civil Engineering, Shanghai Jiao Tong University,Shanghai 200240,China)

Based on three dimensional time domain potential theory,ship motion is solved by using an impulse response function method and boundary element method.Using three dimensional fully nonlinear time domain potential theory to simulate the nonlinear tank sloshing.Then the ship motion coupled with tank sloshing is established that the interactions of wave,ship body and tank sloshing are completely taken into considerations.Numerical simulation and experimental studies indicate that the numerical results of ship motion coupled with tank sloshing can clearly show the coupling effect of tank sloshing on the ship global motion.The results of the ship motion RAO of both computational and experiment ones are in good agreement, the numerical result of ship which has forward speed with liquid tank is also agreed with expectation.

tank sloshing;coupling motion;nonlinear time domain;ship speed

O35

A

10.3969/j.issn.1007-7294.2016.07.003

1007-7294（2016）11-1369-12

2016-01-09

工業(yè)和信息化部高技術(shù)船舶科研項(xiàng)目資助（2007K24161）

李裕龍（1985-），男，博士研究生；朱仁傳（1969-），男，教授，博士生導(dǎo)師，通訊作者，E-mail：renchuan@sjtu.edu.cn。