(南通中遠(yuǎn)海運(yùn)川崎船舶工程有限公司，江蘇 南通 226005)

目前，船舶阻力預(yù)報(bào)最精確的方法主要是船模試驗(yàn)，船舶設(shè)計(jì)初期，在沒(méi)有型線的前提下，無(wú)法應(yīng)用船模試驗(yàn)和CFD仿真來(lái)確定阻力，一般采用近似的方法估算船舶阻力。近似估算阻力的方法有很多，這些估算方法雖然有著耗時(shí)少，成本低等優(yōu)點(diǎn)，但是存在估算精度低，適用范圍有限等不足。改進(jìn)的方法[1-5]雖探討了不同阻力估算公式，但適用范圍有限。為此，考慮從船舶設(shè)計(jì)者角度，尋求一種快速的、精度高的，以及適用范圍廣的預(yù)報(bào)方法。搜集30多艘三大主力船型的資料，形成船型數(shù)據(jù)庫(kù)。采用Holtrop方法分別對(duì)船型數(shù)據(jù)庫(kù)中的散貨船、集裝箱船、油船的阻力進(jìn)行估算，通過(guò)數(shù)學(xué)統(tǒng)計(jì)分析，改進(jìn)相關(guān)參數(shù)的經(jīng)驗(yàn)公式，提出一種輸入?yún)?shù)少，適用范圍廣以及估算精度高的阻力估算的半經(jīng)驗(yàn)公式。

1 Holtrop阻力估算公式

Holtrop法[6-7]其估算公式如下。

Rt=Rf(1+k)+Rapp+Rw+Rb+Rtr+Ra(1)

式中：Rf為摩擦阻力，根據(jù)ITTC1957公式得到；1+k為船舶形狀因子；Rapp為船舶附體阻力；Rw為興波阻力；Rb為球鼻艏引起的附加阻力；Rtr為方艉引起的附加阻力；Ra為船模與實(shí)船相關(guān)的修正因子。

Holtrop給出了式(1)中每個(gè)參數(shù)的回歸公式，綜合這些回歸公式，可以總結(jié)出其總阻力的一般函數(shù)表達(dá)式。

Rt=f(L,B,T,Δ,Cb，Cp，Cm，Lcb，ABT，hB，AT，SApp)

(2)

式中：L、B、T分別為船長(zhǎng)、船寬、吃水；Δ為船舶的排水量；Cb、Cp、Cm分別為方形系數(shù)、棱形系數(shù)以及中橫剖面系數(shù)；Lcb為船舶浮心縱向位置；ABT、hB分別為球鼻艏橫向面積和球鼻艏橫截面中心高度；AT為艉封板浸水面積；SApp為附體濕表面積。

采用Holtrop公式對(duì)船型數(shù)據(jù)庫(kù)中的船舶進(jìn)行阻力估算，并與試驗(yàn)值進(jìn)行對(duì)比，誤差分布見(jiàn)圖1。

圖1 三大主力船型Holtrop法阻力估算誤差分布

由圖1可見(jiàn)，Holtrop法估算船舶阻力的誤差較大，大部分在10%～30%之內(nèi)。其中散貨船和肥大型油船的阻力誤差相對(duì)較小，部分誤差在10%以下；集裝箱船的阻力誤差較大，基本在15%以上。從航速上來(lái)看，速度較低時(shí)，阻力誤差相對(duì)較小，當(dāng)速度增加時(shí)，阻力誤差隨之增大。

2 改進(jìn)Holtrop阻力估算公式

Holtrop公式采用三因次法估算船舶的阻力，除了黏性阻力和興波阻力外，還考慮了球鼻艏、方尾、附體引起的附加阻力以及船模與實(shí)船的尺度效應(yīng)，輸入變量較多，對(duì)于不同船型，其阻力估算的誤差有一定的差別。在船舶設(shè)計(jì)初期，已知的主要參數(shù)有限，尤其是球鼻艏、方尾以及附體的參數(shù)無(wú)從得知，很難利用公式進(jìn)行估算。鑒于此，分別對(duì)球鼻艏、方尾以及附體引起的附加阻力進(jìn)行討論，在保證精度的前提下，尋求一種可靠的替代方案。

1)球鼻艏引起的附加阻力Rb。采用球鼻艏對(duì)于某些船型可以有效降低船舶的興波阻力，這部分的阻力變化在興波阻力中已經(jīng)考慮，而由球鼻艏本身引起的附加阻力較球鼻艏減阻效果為小量，可放到船模與實(shí)船的修正因子Ca中考慮。另外，Holtrop公式中也對(duì)球艏重心高度作了限制，不能超出0.6Tf(艏吃水)，而球艏的重心普遍較高，限制了Rb的計(jì)算。

2)方艉引起的附加阻力Rtr。對(duì)于低速船，方艉處會(huì)形成漩渦，對(duì)阻力不利，這方面的影響可在船舶形狀因子K中考慮；對(duì)于中高速船，方艉會(huì)增加船舶的虛長(zhǎng)度，一定程度上將降低船舶的阻力，這方面的影響可在船模與實(shí)船的修正因子Ca中考慮。

3)附體阻力Rapp。附體對(duì)阻力的影響為小量，主要影響船舶的推進(jìn)性能，因此，附體對(duì)阻力的影響可在船模與實(shí)船的修正因子Ca中考慮。

綜合上述分析，對(duì)Holtrop公式進(jìn)行一定的簡(jiǎn)化，將球鼻艏、方艉,以及附體阻力的影響考慮到船模與實(shí)船的修正因子中，得到船舶總阻力系數(shù)為

Ct=Cf(1+K)+Cw+Ca

(3)

式中：Cf為摩擦阻力系數(shù)，采用ITTC1957公式；Cw為興波阻力系數(shù)；K為船舶形狀因子；Ca為船模與實(shí)船的修正因子。

主要輸入?yún)?shù)見(jiàn)表1。

表1 主要的船型參數(shù)

注：在船舶設(shè)計(jì)初期，不考慮有縱傾的工況；海水密度取常數(shù)，即1.025 kg/m3。

式(2)中的方形系數(shù)、棱形系數(shù)和中橫剖面系數(shù)，采用表1中已知的參數(shù)計(jì)算得到。其中中橫剖面系數(shù)Cm，可以考慮扣除中橫剖面處舭部圓弧(一般圓弧半徑為1.8 m)的面積近似估算得到。

Holtrop公式中沒(méi)有涉及對(duì)浮心縱向位置Lcb的估算，Lcb的值對(duì)摩擦阻力的影響較小，但是其值會(huì)影響船體的艏部和艉部形狀，進(jìn)而影響船舶的興波阻力和黏壓阻力。根據(jù)已有船型數(shù)據(jù)庫(kù)，總結(jié)Lcb與Cb的關(guān)系，見(jiàn)圖2。

圖2 浮心縱向位置隨方形系數(shù)變化規(guī)律

根據(jù)圖2，擬合Lcb的估算公式如下。

2 563.3Cb-387.63

(4)

Holtrop公式中給出的濕表面積公式如下。

0.286 2Cm-0.003 467B/T+

0.369 6Cwp)+2.38ABT/Cb

(5)

對(duì)數(shù)據(jù)庫(kù)中的數(shù)據(jù)進(jìn)行擬合，在式(5)的基礎(chǔ)上進(jìn)行修改，得到濕表面積估算公式如下。

Smod=L(2T+B)(0.623 5+0.333 2Cb-

0.002 611B/T)+BL(0.264 1Cb-0.261 3)

(6)

運(yùn)用式(5)和式(6)估算數(shù)據(jù)庫(kù)中三種主力船型的濕表面積，并與實(shí)際值進(jìn)行對(duì)比，見(jiàn)圖3。

圖3 兩種方法估算不同船型濕表面積誤差對(duì)比

由圖3可見(jiàn)，由式(6)估算得到的濕表面積的誤差明顯降低。

對(duì)于形狀因子，Holtrop方法中給出了1+K的估算公式。

1+K=0.93+0.487 118c14(B/L)1.068 06·

(T/L)0.461 06(L/LR)0.121 563·

(L3/▽)0.364 86(1-Cp)-0.604 247

(7)

式中：LR為Holtrop公式中定義的船型因子，與L、Cp、Lcb有關(guān)。運(yùn)用式(7)對(duì)現(xiàn)有船型的形狀因子進(jìn)行估算，并與數(shù)據(jù)庫(kù)中的數(shù)據(jù)進(jìn)行對(duì)比。通過(guò)分析對(duì)比，修改式(7)，得到改正后的形狀因子公式如下。

1+Kmod=1.046 6+0.312(B/L)1.068 06·

(T/L)0.461 06(L/LR)0.121 563·

(L3/▽)0.364 86(1-Cp)-0.604 247

(8)

興波阻力估算仍采用Holtrop法中給出的公式。式(3)中并未單獨(dú)考慮球艏、方尾、附體等的影響，而是將這些影響放到了船模與實(shí)船的修正因子Ca中。通過(guò)對(duì)比改進(jìn)后的經(jīng)驗(yàn)公式和船型數(shù)據(jù)庫(kù)中船舶試驗(yàn)結(jié)果，發(fā)現(xiàn)Ca的值不僅與船長(zhǎng)和船速有關(guān)，還與船型有關(guān)。根據(jù)水池試驗(yàn)的結(jié)果和式(3)，回歸分析得到不同船型修正因子Ca的關(guān)系。

1)集裝箱船(R2=0.83)。

Ca×103=-0.434 2-13.61Fr2+

1.858Fr+125.6/L

(9)

2)散貨船(R2=0.71)。

Ca×103=0.256 1-26.06Fr2+

0.800 6Fr+47.95/L

(10)

3)油船(R2=0.69)。

Ca×103=-1.482-65.07Fr2+

18.11Fr+109.9/L

(11)

式中：Fr為弗勞德數(shù)；R2為回歸系數(shù)，表示模型擬合度的好壞，其值越接近1，說(shuō)明擬合效果越好，從回歸系數(shù)可以看出，回歸公式的擬合度較好。

3 改進(jìn)的船舶阻力公式驗(yàn)證

為進(jìn)一步驗(yàn)證上述公式的精度及適用性，以不屬于船型數(shù)據(jù)庫(kù)中的3艘主力船為例，采用上述改進(jìn)的半經(jīng)驗(yàn)公式進(jìn)行估算船舶的阻力。3艘船的主要參數(shù)見(jiàn)表2。估算結(jié)果見(jiàn)表3～5。

表2 三種主力船型基本參數(shù)

表3 集裝箱船阻力估算值與試驗(yàn)值對(duì)比

表4 散貨船阻力估算值與試驗(yàn)值對(duì)比

表5 油船阻力估算值與試驗(yàn)值對(duì)比

由表3～5可見(jiàn)，改進(jìn)后的Holtrop法估算結(jié)果對(duì)3種主力船型都有很好的適用性，并且誤差也得到改善，均在7.0%以內(nèi)，達(dá)到一定的修正效果。

4 結(jié)論

通過(guò)對(duì)比現(xiàn)有新造船的船型數(shù)據(jù)庫(kù)，基于三大主力船型的試驗(yàn)數(shù)據(jù)，對(duì)Holtrop法中部分公式進(jìn)行改進(jìn)，修正船模與實(shí)船相關(guān)因子，并減少該方法的輸入?yún)?shù)。經(jīng)過(guò)驗(yàn)證，采用改進(jìn)的方法，3種船型的阻力估算值的精度明顯得到提高，具有一定的工程應(yīng)用價(jià)值。

將已有的估算公式與已建造船的數(shù)據(jù)結(jié)合起來(lái)，通過(guò)總結(jié)和歸納的方法，尋求適合新造船阻力估算的方法，在船舶設(shè)計(jì)初期不失為一種有效的阻力估算方法。鑒于本文的船型數(shù)據(jù)庫(kù)樣本容量有限，統(tǒng)計(jì)公式仍需進(jìn)一步完善。另外，本文僅分析了阻力的估算方法，對(duì)推進(jìn)方面的歸納還需進(jìn)一步探討。