朱軍，朱韜，王智宇，夏齊強(qiáng)，黃昆侖，葛義軍

1 海軍工程大學(xué)艦船與海洋學(xué)院，湖北武漢430033

2 海軍研究院特種勤務(wù)研究所，上海200235

0 引 言

半個(gè)多世紀(jì)前，研究人員就注意到，船舶在沒(méi)有橫向作用力的縱向波浪中會(huì)發(fā)生大幅度的橫搖運(yùn)動(dòng)，即船舶參數(shù)橫搖現(xiàn)象。Kerwin［1］采用一階諧波形式表達(dá)了回復(fù)力矩的變化量，運(yùn)用馬蒂厄方程構(gòu)建了橫搖運(yùn)動(dòng)模型；Paulling 等［2］認(rèn)為，垂蕩和縱搖的非線性耦合至橫向會(huì)引發(fā)橫搖運(yùn)動(dòng)，其在垂蕩或縱搖為諧振運(yùn)動(dòng)模式的假定下，導(dǎo)出了一階諧波回復(fù)力矩形式的馬蒂厄方程，經(jīng)模型試驗(yàn)，測(cè)定了靜水中由垂蕩引發(fā)的參數(shù)橫搖運(yùn)動(dòng)。

1998 年，巴拿馬型C11 大型集裝箱船因受極端氣象影響而發(fā)生了參數(shù)橫搖運(yùn)動(dòng)，導(dǎo)致甲板集裝箱遭受巨大損失，再次引起了人們對(duì)船舶參數(shù)橫搖運(yùn)動(dòng)的關(guān)注。有關(guān)船舶參數(shù)橫搖運(yùn)動(dòng)計(jì)算，可以歸納為2 類：一是采用馬蒂厄方程構(gòu)建橫搖運(yùn)動(dòng)模型，其中穩(wěn)性參數(shù)采用遭遇頻率的一階諧波模式；二是采用垂蕩、縱搖和橫搖的耦合數(shù)值計(jì)算模型，即直接的搖蕩耦合運(yùn)動(dòng)數(shù)值模擬。

France 等［3］對(duì)C11 大型集裝箱船遭遇的海況進(jìn)行了調(diào)查，確認(rèn)了參數(shù)橫搖發(fā)生的海況條件。Shin 等［4-6］針對(duì)船舶發(fā)生參數(shù)橫搖運(yùn)動(dòng)的現(xiàn)象給出了相同的物理解釋，即波浪通過(guò)船體時(shí)水線面面積會(huì)發(fā)生周期性的改變，并用遭遇頻率的一階諧波模式描述變化的穩(wěn)性參數(shù)，以垂蕩和縱搖的靜態(tài)或準(zhǔn)靜態(tài)平衡方式估計(jì)穩(wěn)性參數(shù)的改變量，以馬蒂厄標(biāo)準(zhǔn)方程的數(shù)值圖解作為穩(wěn)定性的判據(jù)。Belenky 等［5］采用Spyrou［7］給出的參數(shù)橫搖判別的衡準(zhǔn)取決于穩(wěn)性參數(shù)改變量的大小。

在數(shù)值計(jì)算方面，F(xiàn)rance 等［3］采用多個(gè)耐波性專業(yè)軟件對(duì)參數(shù)橫搖運(yùn)動(dòng)進(jìn)行了計(jì)算，這些軟件盡管不是完全針對(duì)參數(shù)橫搖運(yùn)動(dòng)開(kāi)發(fā)的，但其計(jì)算結(jié)果與試驗(yàn)結(jié)果在趨勢(shì)上相同。傅汝德-克雷洛夫波浪力被認(rèn)為是船舶發(fā)生參數(shù)橫搖運(yùn)動(dòng)的關(guān)鍵。Ma 等［8-10］采用瞬時(shí)三維壓力分布的濕表面積積分計(jì)算波浪力，計(jì)算研究了船舶的參數(shù)橫搖運(yùn)動(dòng)。魯江和卜淑霞等［11-12］運(yùn)用單自由度（類似馬蒂厄方程）和三自由度數(shù)學(xué)模型模擬了參數(shù)橫搖運(yùn)動(dòng)，并在三自由度模型中采用STF 方法處理了波浪力，結(jié)果顯示，與單自由度模型相比，三自由度模型的效果并非更好。

上述研究均未指明參數(shù)橫搖運(yùn)動(dòng)的能量來(lái)源與機(jī)制（力學(xué)機(jī)理），Shin 和Belenky 等［4-5］只是泛泛地認(rèn)為穩(wěn)性參數(shù)的變化是橫搖運(yùn)動(dòng)的能量來(lái)源，這個(gè)機(jī)理的認(rèn)識(shí)是數(shù)學(xué)性的。另外，在國(guó)際海事組織（IMO）正在制定的船舶第2 代完整穩(wěn)性標(biāo)準(zhǔn)中，參數(shù)橫搖的薄弱性衡準(zhǔn)是采用穩(wěn)性變化量的大小來(lái)衡量［6］，然而，不同類型的船舶其計(jì)算結(jié)果離散較大，這表明以穩(wěn)性變化量大小作為參數(shù)橫搖衡準(zhǔn)只是權(quán)宜辦法［5］。

本文擬采用所提出的搖蕩耦合切片計(jì)算方法［13］，克服大幅運(yùn)動(dòng)時(shí)參數(shù)含義的模糊性，以及會(huì)誘發(fā)更多自由度運(yùn)動(dòng)耦合的缺陷，數(shù)值計(jì)算船舶的參數(shù)橫搖運(yùn)動(dòng)，由此分析發(fā)生參數(shù)橫搖運(yùn)動(dòng)的力學(xué)機(jī)理，并基于能量原理提出船舶參數(shù)橫搖的衡準(zhǔn)。

1 坐標(biāo)系及廣義參數(shù)

1.1 坐標(biāo)系

3 個(gè)坐標(biāo)系分別為慣性坐標(biāo)系、船體運(yùn)動(dòng)坐標(biāo)系和船體固定坐標(biāo)系，如圖1 所示。

圖1 坐標(biāo)系Fig.1 Coordinate system

慣性坐標(biāo)系E-ξηζ的原點(diǎn)E固定于靜水面，Eξ軸指向船舶前進(jìn)方向，Eζ軸垂直向上，Eη軸指向船體右舷。

船體運(yùn)動(dòng)坐標(biāo)系M-xyz在初始時(shí)刻與慣性坐標(biāo)系重合，坐標(biāo)系沿慣性坐標(biāo)系的Eξ軸方向以速度U做直線運(yùn)動(dòng)。

船體固定坐標(biāo)系B-xByBzB固定于船體，原點(diǎn)位于船體龍骨基線左右對(duì)稱面的船舯處，BzB軸線垂直向上，ByB軸線指向船體右舷，BxB軸線指向船艏。

1.2 廣義參數(shù)

繞船體軸線轉(zhuǎn)動(dòng)定義的常規(guī)橫傾角和縱傾角只適合小角度情況，它存在轉(zhuǎn)動(dòng)順序問(wèn)題，即轉(zhuǎn)動(dòng)順序不同，船體姿態(tài)亦不同。在極端情況，例如當(dāng)橫傾角和縱傾角均為90°時(shí)，先后不同地轉(zhuǎn)動(dòng)橫傾角和縱傾角，船體姿態(tài)會(huì)有很大的差異。為此，本文定義橫傾角繞船體縱向坐標(biāo)軸轉(zhuǎn)動(dòng)，縱傾角繞靜水面橫軸線（My 軸）轉(zhuǎn)動(dòng)，具體定義如下。

假定船舶靜水漂浮時(shí)吃水為T，船體繞Mx軸線轉(zhuǎn)動(dòng)橫搖角φ，再沿Mz軸平移吃水增量ΔT，最后再繞My軸轉(zhuǎn)動(dòng)縱傾角θ（圖2）。這里，稱ΔT，θ分別為廣義吃水增量和廣義縱傾角。根據(jù)轉(zhuǎn)動(dòng)關(guān)系，不難得到船體運(yùn)動(dòng)坐標(biāo)系和船體固定坐標(biāo)系的轉(zhuǎn)換關(guān)系為：

圖2 坐標(biāo)系轉(zhuǎn)動(dòng)與平移Fig.2 Coordinate systems transformation by rotation and translation

2 船體搖蕩耦合運(yùn)動(dòng)方程組

習(xí)慣上，將運(yùn)動(dòng)方程構(gòu)建在船體固定坐標(biāo)系中，其優(yōu)點(diǎn)是繞坐標(biāo)軸的轉(zhuǎn)動(dòng)慣量不隨運(yùn)動(dòng)變化，其為常數(shù)，缺點(diǎn)是大幅度運(yùn)動(dòng)時(shí)易誘發(fā)更多自由度的運(yùn)動(dòng)耦合。例如，頂浪航行時(shí)，垂向波浪力在轉(zhuǎn)換到船體固定坐標(biāo)系后會(huì)有沿船體橫向的分量，這個(gè)力會(huì)導(dǎo)致船體固定坐標(biāo)下的橫蕩和搖艏運(yùn)動(dòng)。這就意味著原本在靜水面觀測(cè)到的3 個(gè)自由度運(yùn)動(dòng)，在船體固定坐標(biāo)下會(huì)演變?yōu)? 個(gè)自由度的運(yùn)動(dòng)，這將導(dǎo)致交叉耦合項(xiàng)的水動(dòng)力出現(xiàn)，使耦合運(yùn)動(dòng)的計(jì)算變得更加復(fù)雜。

因此，本文將搖蕩耦合運(yùn)動(dòng)方程構(gòu)建在船體運(yùn)動(dòng)坐標(biāo)的慣性系中。在該慣性坐標(biāo)系下，隨時(shí)間不斷變化的是船體水下形狀，在垂直面內(nèi)，船體重力和入射波浪力作用會(huì)產(chǎn)生垂蕩運(yùn)動(dòng)，重力和入射波浪力對(duì)船舯的力矩會(huì)產(chǎn)生縱搖運(yùn)動(dòng)。采用普通切片法，將切片的位移、速度和加速度的作用力沿船體縱向積分得到作用于船體的力/力矩，從而不難得到垂蕩和縱搖運(yùn)動(dòng)方程式為：

式中：ZΔT?，ZΔT?，Zθ?，Zθ?，ZUθ?，ZUθ為船體垂蕩的水動(dòng)力系數(shù)；Mθ?，Mθ?，MΔT?，MΔT?，MUθ?，MUθ為船體縱搖的水動(dòng)力系數(shù)；Fζ和Mζ，F(xiàn)ζ?和Mζ?，F(xiàn)ζ?和Mζ?分別為入射波浪相對(duì)船體位移、速度和加速度的垂向力及縱搖力矩；m為船體質(zhì)量；Jθ為船體縱向轉(zhuǎn)動(dòng)慣量；xG為船體重心縱向坐標(biāo)位置；ρg?為船體在靜水中的排水量，其中ρ為水的密度，g為重力加速度，?為靜水中的排水體積。

橫搖運(yùn)動(dòng)方程將忽略入射波浪對(duì)船體速度和加速度產(chǎn)生的偏心作用，只考慮入射波浪力的位移作用。因此，繞船體重心的橫搖運(yùn)動(dòng)方程式可寫(xiě)為

式中：ΔIφ和Kφ?分別為橫搖附加慣性矩和阻尼系數(shù)；Iφ為船體橫向轉(zhuǎn)動(dòng)慣量；yζ為作用力的橫向坐標(biāo)位置；yG船體重心橫向坐標(biāo)位置。

式（2）和式（3）即構(gòu)成了船體垂蕩、縱搖和橫搖的耦合運(yùn)動(dòng)方程組。值得注意的是，F(xiàn)ζ是將船體姿態(tài)（ΔT，θ，φ）和瞬時(shí)波浪位置ζ作為一個(gè)整體來(lái)計(jì)算的，方程組通過(guò)力Fζ=f(ΔT，θ，φ，ζ)耦合了船體的垂蕩、縱搖和橫搖運(yùn)動(dòng)，具體計(jì)算公式將在后文給出。

3 水動(dòng)力計(jì)算式

3.1 瞬時(shí)波面下壓力分布

由勢(shì)流理論，可知船體運(yùn)動(dòng)坐標(biāo)系下小振幅規(guī)則波升ζ可表達(dá)為

水下壓力p的分布式為e 指數(shù)形式

式中：k為波數(shù)；χ為浪向角；t 為時(shí)間；A 為波幅；ωe為遭遇頻率；Θ為計(jì)入了時(shí)間和位置的波浪相位。為方便起見(jiàn)，令Θ=ωet+k(xcosχ+ysinχ) 。式（4）并不適合直接做壓力積分計(jì)算，為此，通常采用靜水壓力和拉伸修正的近似方法［14］，或不做修正。本文采用線性近似，為此，做一個(gè)變量代換：

z=z1+AcosΘ

該代換的含義是，z坐標(biāo)是以靜水面為零點(diǎn)，而z1坐標(biāo)則是以瞬時(shí)波面為零點(diǎn)。將其代入式（3），忽略高階項(xiàng)后，可得壓力的近似表達(dá)式為

p=-ρ′gz1（5）

式中，ρ′=ρ(1-α0cosΘ) ，為等效水的密度，其中α0=Ak，為最大波傾角。式（4）滿足瞬時(shí)波面（z1=0）處壓力p=0，這就意味著瞬時(shí)波面下的壓力用等效靜水壓力近似，等效靜水的密度ρ′在波峰區(qū)域（cosΘ＞0）較實(shí)際水的密度ρ小，在波谷區(qū)域（cosΘ＜0）較實(shí)際水的密度ρ大。

3.2 入射波浪的位移力

利用高斯定理，將瞬時(shí)波面下船體濕表面的壓力積分轉(zhuǎn)換成體積分：

式中，F(xiàn)和M分別為作用于船體的力和繞坐標(biāo)原點(diǎn)的力矩；s為船體濕表面的面積；V為波浪下船體濕表面的體積；n為船體表面的單位外法向矢量；r為微分濕表面ds到坐標(biāo)原點(diǎn)的位置矢量；i，j，k分別為船體運(yùn)動(dòng)坐標(biāo)系的3 個(gè)單位坐標(biāo)矢量。將壓力的表達(dá)式（4）代入到上述力的計(jì)算公式中，不難導(dǎo)出瞬時(shí)波面下船體入射波浪3 個(gè)方向的作用力為：

式中：L 為船長(zhǎng)；Sx為瞬時(shí)波面下船體切片的面積（圖3 中陰影部分面積）。式（7）表明，瞬時(shí)波面下船體入射波浪的作用力只有垂向力和縱向力，橫向力為0。顯然，縱向力Fx與波傾角α0成正比，因波傾角α0為小量，因此縱向力Fx是α0級(jí)小量，波浪的主要作用力是垂向力Fz。

圖3 瞬時(shí)波面下船體切片示意圖Fig.3 Diagram of hull strip under transient wave surface

同樣，可以導(dǎo)出繞坐標(biāo)軸的橫搖力矩K 和縱搖力矩M：

式中，xs，ys和zs分別為瞬時(shí)波面下船體切片面積中心在船體運(yùn)動(dòng)坐標(biāo)系的坐標(biāo)值。其中，縱搖力矩M由縱向力和垂向力2 部分構(gòu)成，也即式（2）中的Mζ。由此，由式（7）和式（8）不難確定出垂向力的橫向坐標(biāo)位置yζ=K/Fz。推導(dǎo)中，利用了頂浪航行條件sinχ=0。

3.3 船體附加質(zhì)量和阻尼系數(shù)估算

式（2）需要計(jì)算船體的附加質(zhì)量和阻尼系數(shù)等，本文則將按經(jīng)驗(yàn)公式來(lái)估算。本節(jié)中的經(jīng)驗(yàn)公式和近似圖譜來(lái)自于文獻(xiàn)［15］。

1）船體橫搖附加慣性矩和阻尼系數(shù)。

船體橫搖的慣性矩、附加慣性矩和阻尼系數(shù)按照常規(guī)單自由度運(yùn)動(dòng)方式，以船體平均位置（靜水平衡狀態(tài)）近似估算。橫搖附加慣性矩和阻尼系數(shù)由經(jīng)驗(yàn)公式估算：

式中：ρφ為橫搖慣性半徑，ρφ=CB，其中B為船寬，經(jīng)驗(yàn)系數(shù)C的取值范圍為0.33～0.45；2μφ為無(wú)量綱橫搖衰減系數(shù)，其取值范圍為0.11～0.14；GM0為初穩(wěn)性高；Fn 為傅汝德數(shù)。

2）切片垂蕩附加質(zhì)量和阻尼。

用二維浮體的寬度、寬吃水比和面積系數(shù)這3 個(gè)參數(shù)來(lái)估算船體切片的垂向附加質(zhì)量和阻尼系數(shù)，沿船體縱向積分得出船體的垂蕩附加質(zhì)量、縱搖附加慣性矩和對(duì)應(yīng)的阻尼系數(shù)。

二維浮體單位長(zhǎng)度的附加質(zhì)量Δm和阻尼系數(shù)ΔNμ估算公式為

式中：Bn為二維浮體水面處的寬度；C1為附加質(zhì)量系數(shù)；Aˉ為衰減系數(shù)，是輻射波幅值與垂蕩幅值之比。附加質(zhì)量系數(shù)C1和衰減系數(shù)Aˉ的大小與遭遇頻率ωe、寬吃水比和水下面積系數(shù)有關(guān)，本文采用文獻(xiàn)［15］的圖譜插值近似。

3）船體附加質(zhì)量和阻尼系數(shù)。

對(duì)式（2）等號(hào)左邊的船體垂蕩水動(dòng)力系數(shù)和縱搖水動(dòng)力系數(shù)，按照切片在瞬時(shí)波面下的船體姿態(tài)（ΔT，θ，φ），由式（10）來(lái)估算切片附加質(zhì)量Δm和阻尼系數(shù)ΔNμ，然后再沿船體縱向積分得到船體垂蕩和縱搖的附加質(zhì)量及附加慣性矩等。

入射波浪相對(duì)于船體速度和加速度的垂向力及縱搖力矩也可按照同樣的方法積分計(jì)算。

4 參數(shù)橫搖數(shù)值計(jì)算及力學(xué)機(jī)理分析

本文采用上述動(dòng)力學(xué)模型和基于圖形面域技術(shù)開(kāi)發(fā)的軟件計(jì)算了一艘船舶在頂浪規(guī)則波中的搖蕩運(yùn)動(dòng)，以此分析船舶發(fā)生參數(shù)橫搖運(yùn)動(dòng)的力學(xué)機(jī)理。

4.1 參數(shù)橫搖運(yùn)動(dòng)計(jì)算

計(jì)算船舶的船長(zhǎng)L=142 m，航速U=14 kn，波陡范圍Hw/λ=0.015～0.050。

圖4 所示為波陡Hw/λ=0.025 時(shí)的垂蕩、縱搖和橫搖計(jì)算時(shí)歷曲線。計(jì)算結(jié)果顯示：在該狀態(tài)下呈現(xiàn)出了參數(shù)橫搖運(yùn)動(dòng)，橫搖幅值不斷增加；同時(shí)還表明，橫搖頻率為縱搖頻率的一半，且位于橫搖共振頻率范圍內(nèi)，這與發(fā)生參數(shù)橫搖的頻率條件是吻合的。

圖4 頂浪中的搖蕩運(yùn)動(dòng)計(jì)算時(shí)歷曲線（Hw/λ=0.025）Fig.4 The calculated time history curves of oscillating motions in head wave（Hw/λ=0.025）

表1 參數(shù)橫搖運(yùn)動(dòng)的預(yù)報(bào)結(jié)果Table 1 Prediction results of parametric rolling motion

表1 所示為不同波陡和波長(zhǎng)條件下的參數(shù)橫搖運(yùn)動(dòng)預(yù)報(bào)結(jié)果。表中：符號(hào)“●”表示橫搖運(yùn)動(dòng)振蕩增加，呈現(xiàn)參數(shù)橫搖運(yùn)動(dòng)；符號(hào)“×”表示橫搖運(yùn)動(dòng)衰減，無(wú)參數(shù)橫搖運(yùn)動(dòng)。預(yù)報(bào)結(jié)果顯示，除較小的波陡不發(fā)生參數(shù)橫搖運(yùn)動(dòng)外，隨著波陡的增加，發(fā)生參數(shù)橫搖的波段向高頻區(qū)偏移，整體上波長(zhǎng)/船長(zhǎng)（λ/L）在0.9～1.4 范圍。

4.2 參數(shù)橫搖力學(xué)機(jī)理分析

將式（3）第3 項(xiàng)的回復(fù)力矩寫(xiě)成通常的形式：

式中，參數(shù)Kφ為回復(fù)力矩系數(shù)，相當(dāng)于船舶橫搖系統(tǒng)剛度，其作用是在每半個(gè)橫搖周期內(nèi)，在橫搖角增大的過(guò)程中吸收能量，在橫搖角減小的過(guò)程釋放所吸收的能量。

圖5 回復(fù)力矩系數(shù)和橫搖角計(jì)算曲線（λ/L=1.2）Fig.5 calculated curves of recovery moment coefficient and rolling angle（λ/L=1.2）

圖6 回復(fù)力矩系數(shù)和橫搖角計(jì)算曲線（λ/L=1.5）Fig.6 calculated curve of recovery moment coefficient and rolling angle（λ/L=1.5）

圖5 和圖6 所示為Hw/λ=0.025 時(shí)回復(fù)力矩系數(shù)Kφ與橫搖角φ隨時(shí)間變化的計(jì)算曲線（其中Kφ的單位為N·m，φ的單位為（°），圖5、圖6 的縱坐標(biāo)為它們各自乘以不同系數(shù)后的尺度）。圖5呈現(xiàn)的是參數(shù)橫搖運(yùn)動(dòng)，其λ L=1.2；而圖6 呈現(xiàn)的則是橫搖運(yùn)動(dòng)振蕩衰減，其λ L=1.5。從圖5 所示的計(jì)算曲線可以看出，在橫搖角增加過(guò)程中，參數(shù)Kφ低于平均值，處于較低水平；而在橫搖角減小的過(guò)程中，參數(shù)Kφ高于平均值，處于較高水平。這就意味著在每半個(gè)橫搖周期內(nèi)，在橫搖角增大過(guò)程中所吸收的能量小于橫搖角減小過(guò)程中釋放的能量；圖6 的情況與之相反，因而不會(huì)出現(xiàn)參數(shù)橫搖運(yùn)動(dòng)。

因此，參數(shù)橫搖的力學(xué)機(jī)理是：回復(fù)力矩系數(shù)Kφ隨著波浪的作用呈現(xiàn)出圍繞均值的波動(dòng)，在橫搖角減小過(guò)程中，其釋放的能量要大于橫搖角增大過(guò)程中回復(fù)力矩系數(shù)吸收的能量。該能量差值也即參數(shù)橫搖運(yùn)動(dòng)持續(xù)放大的能量來(lái)源。

4.3 參數(shù)橫搖衡準(zhǔn)

圖5 和圖6 顯示，在不同波長(zhǎng)/船長(zhǎng)、相同波陡情況下，回復(fù)力矩系數(shù)Kφ的均值與振蕩幅值沒(méi)有明顯差異。如上述分析，是否發(fā)生參數(shù)橫搖運(yùn)動(dòng)取決于系統(tǒng)吸收和釋放能量的平衡與否，系數(shù)Kφ與橫搖角φ隨時(shí)間變化的相對(duì)位置，即相位關(guān)系決定了吸收與釋放能量的相對(duì)大小。現(xiàn)用遭遇頻率ωe的諧波函數(shù)擬合系數(shù)Kφ：

任意一個(gè)橫搖半周期內(nèi)，假定橫搖角為正弦函數(shù)，其頻率為用遭遇頻率ωe的一半，即

式中，φ0為半周期內(nèi)的橫搖角幅度，為常數(shù)。因此不難得到，橫搖角半周期內(nèi)回復(fù)力矩系數(shù)Kφ吸收和釋放的能量總和為

式（13）表明，只有系數(shù)Kφ的一階諧波分量的積分不為0，常數(shù)項(xiàng)和高階諧波項(xiàng)的積分均等于0。發(fā)生參數(shù)橫搖的條件是能量J＜0。因此，滿足此條件的相位角ε1可由式（12）得到

這就是參數(shù)橫搖的判別衡準(zhǔn)。

如果將式（11）只保留常數(shù)項(xiàng)和一階諧波項(xiàng)，則與France 等［3-5］研究的參數(shù)橫搖所采用的馬蒂厄方程相類似。但是，判別參數(shù)橫搖的衡準(zhǔn)不同，是否發(fā)生參數(shù)橫搖不取決于一階諧波項(xiàng)系數(shù)的大小，而是一階諧波項(xiàng)的相位角ε1。相位角ε1的變化取決于船波相對(duì)運(yùn)動(dòng)，其與波浪中船舶的垂蕩、縱搖和由入射波浪引起的水線面形狀變化以及水下排水體積變化相關(guān)。

5 結(jié) 論

本文建立了慣性坐標(biāo)下的垂蕩和縱搖耦合運(yùn)動(dòng)方程，以及船體坐標(biāo)下的橫搖運(yùn)動(dòng)方程，由此構(gòu)成垂蕩、縱搖和橫搖的混合動(dòng)力學(xué)模型，并采用普通切片方法計(jì)算了入射波浪力、附加質(zhì)量和阻尼系數(shù)。數(shù)值計(jì)算了一艘船舶的垂蕩、縱搖和橫搖耦合運(yùn)動(dòng)，預(yù)報(bào)了參數(shù)橫搖運(yùn)動(dòng)，通過(guò)對(duì)數(shù)值計(jì)算結(jié)果的分析，得出了參數(shù)橫搖運(yùn)動(dòng)的力學(xué)機(jī)理和參數(shù)橫搖判別衡準(zhǔn)，研究認(rèn)為：

1）參數(shù)橫搖的力學(xué)機(jī)理是：由于回復(fù)力矩系數(shù)Kφ是隨時(shí)間變化的，導(dǎo)致在橫搖過(guò)程中回復(fù)力矩系數(shù)吸收和釋放的能量不等，釋放能量大于吸收能量是發(fā)生參數(shù)橫搖的力學(xué)機(jī)理，也是發(fā)生參數(shù)橫搖后運(yùn)動(dòng)持續(xù)放大的能量來(lái)源。

2）判別參數(shù)橫搖的衡準(zhǔn)是，回復(fù)力矩系數(shù)Kφ的一階諧波分量具有超前的相位角ε1?[0，π]。

船舶參數(shù)橫搖力學(xué)機(jī)理的明確有助于對(duì)參數(shù)橫搖運(yùn)動(dòng)的認(rèn)識(shí)，而如何建立相位角ε1與船體參數(shù)的關(guān)系尚需深入研究。