匡培欽，柳存根，汪學(xué)鋒，周翀劍

基于數(shù)值模擬的水翼五體船耐波性分析

匡培欽，柳存根，汪學(xué)鋒，周翀劍

為改善五體船在小搖蕩時的耐波性，提出在側(cè)體和主船體之間加裝水翼構(gòu)造水翼五體船，借助HydroStar水動力學(xué)軟件，依托Wigley片體五體船進(jìn)行水翼五體船耐波性計算對比。結(jié)果表明，由于水翼的減搖減震及升力作用，在側(cè)體和主船體之間加裝水翼可以明顯降低五體船的縱搖幅值、橫搖幅值及垂蕩幅值，提高五體船的耐波性。

五體船；水翼；縱搖穩(wěn)性；垂蕩；橫搖穩(wěn)性

隨著高性能船舶設(shè)計研發(fā)熱度的增加，五體船以其在快速性、耐波性及節(jié)能減排方面的出色表現(xiàn)，逐漸引起國內(nèi)外多方學(xué)者的關(guān)注[1-4]。經(jīng)過多年的船型優(yōu)化，五體船在側(cè)體分布方面已經(jīng)相對比較成熟，2個后側(cè)體保證了其快速前進(jìn)時的橫搖穩(wěn)定性，而在發(fā)生縱搖或垂蕩運動時，2個前側(cè)體入水，大大減小其縱搖和垂蕩幅值[5]。

在空載或小搖蕩時由于2個前側(cè)體在水線之上，則五體船各性能與三體船無異[6]，抗風(fēng)浪性能相應(yīng)削弱。1981年美國華盛頓大學(xué)的Calkins首次提出多體水翼船的概念[7-8]，多體水翼船是水翼船與多體船的融合船型，兼具多體船和水翼船的優(yōu)點。基于此，提出水翼五體船的船型，以Wigley片體五體船為依托，通過在側(cè)體與主船體之間加裝水翼構(gòu)造水翼五體船；并借助耐波性通用軟件HydroStar計算在波浪中的縱搖、垂蕩、橫搖等運動，以此研究加裝水翼對五體船耐波性的改善情況。

1 數(shù)值模擬基本理論

1.1 橫搖運動方程

(1)

由am=am0sin (ωt)可將式(1)改寫為

(2)

得到橫搖運動方程：

(3)

1.2 縱搖、垂蕩耦合運動方程

根據(jù)切片法思想，某處垂向相對位置為

(4)

式中：X為縱向位置；θ為縱搖角。

等效波面方程為

則可得切片所受流體動力如下。

流體靜力：

(6)

興波阻力：

(7)

附加慣性力：

(8)

式中：NH為阻尼系數(shù)；MH為附加質(zhì)量；V為船速。

沿船長積分可得垂蕩力及縱搖力矩：

(9)

根據(jù)牛頓第二定律可得船體垂蕩和縱搖運動方程為

(10)

式中：D為排水量；IYY為縱向慣性矩。

2 模型參數(shù)及計算設(shè)置

2.1 模型參數(shù)

主要對比研究加裝水翼對五體船耐波性的改善問題，采用Wigley船型作為主船體及側(cè)體船型進(jìn)行研究。Wigley船型瘦長，符合線性理論小擾動假定[9]。船型參數(shù)見表1。

表1 船型參數(shù)

模型水翼采用NACA 4位數(shù)字標(biāo)準(zhǔn)翼型。NACA 4位數(shù)字翼型是美國NACA最早建立的一個低速翼型系列，有較高的最大升力系數(shù)和較低的阻力系數(shù)，適合用于快速艇水翼[10]。當(dāng)相對低速時，NACA4415翼型的失速攻角為14°，且在-4°～14°范圍內(nèi)升力系數(shù)隨攻角增大而增大[11]。保守起見，采用水翼攻角10°。水翼參數(shù)見表2。

表2 翼型NACA4415水翼參數(shù)

2.2 計算設(shè)置

五體船的側(cè)體布局優(yōu)化基本朝同一趨勢發(fā)展，前側(cè)體縱向位置基本穩(wěn)定在主船體船舯附近，后側(cè)體縱向位置則主要集中于主船體艉部，若假定前后側(cè)體的橫向坐標(biāo)一致，則其大致分布于距船舯縱剖面1.5倍半船寬位置[12]。據(jù)此將前側(cè)體置于主船體船舯稍前位置，后側(cè)體置于船艉，前后側(cè)體橫向位置一致，位于1.5倍船寬處。見圖1。

水翼均置于側(cè)體中間部分，用于連接4個側(cè)體和主船體。后水翼橫向水平，前水翼從前側(cè)體向下斜插入水中。

不同海域海況不同，以北大西洋為例，波高在5 m以下出現(xiàn)的概率為90.42%，其中2～3 m波高出現(xiàn)頻率最高為30.10%；波浪周期在11 s以下的概率為96.17%，其中波浪周期在7～9 s范圍內(nèi)出現(xiàn)頻率最高，可達(dá)35.87%[13]。選取周期為1～10 s有義波高為3 m的非規(guī)則波作為計算海況，分別計算迎浪、斜浪和橫浪3種情況下有無水翼2種五體船型耐波性指標(biāo)，進(jìn)行對比研究。

3 計算結(jié)果及分析

由于后水翼橫向水平，前水翼向中間橫向斜插入水面以下，因此前后水翼主要對縱搖、垂蕩，以及橫搖產(chǎn)生較大影響。

3.1 縱搖分析

為了便于不同波幅下對比分析，引入標(biāo)準(zhǔn)化變量縱搖響應(yīng)因子HP。

(11)

式中：H為縱搖角度幅值，(°)；A為有義波高，m。

縱搖主要受迎浪或斜浪的縱向分量影響，因此對于縱搖這里僅對斜浪和迎浪情況進(jìn)行對比分析。在斜浪和迎浪2種浪向下，有無水翼兩種五體船縱搖狀態(tài)見圖2。

由圖2可以看出，在低頻波段和高頻波段加裝水翼的水翼五體船縱搖響應(yīng)因子整體要高于無水翼的五體船，但在中頻波段，水翼五體船的縱搖響應(yīng)因子要明顯低于五體船。且無論是斜浪還是迎浪浪向，水翼五體船縱搖響應(yīng)因子的最大值總要小于五體船。按縱搖響應(yīng)因子最大值分析，斜浪下五體船加裝水翼使得減縱搖效果高達(dá)21.8%，而迎浪下減縱搖效果可達(dá)26.6%。水翼五體船后水翼橫向水平，在縱搖運動時除提供升力外還起到了減搖鰭的作用；而前水翼由于半出水，在艏揚時出水增加，所提供升力減小，在艏沉?xí)r，出水減少，所提供升力增加，從而增加了縱搖運動時的回復(fù)力矩，使得水翼五體船的縱搖穩(wěn)定性要優(yōu)于無水翼的五體船。

3.2 垂蕩分析

引入標(biāo)準(zhǔn)化變量垂蕩響應(yīng)因子PP。

(12)

式中：P為垂蕩，m；A為有義波高，m。

在斜浪和迎浪2種浪向下，有無水翼2種五體船垂蕩狀態(tài)見圖3。

由圖3可見，無論斜浪還是迎浪情況下，在低、中頻波段，加裝水翼的水翼五體船和無水翼的五體船垂蕩響應(yīng)因子隨波頻的增加急劇下降，且水翼五體船的垂蕩響應(yīng)因子一直低于五體船。在高頻波段，2種船型的垂蕩響應(yīng)因子均趨于平穩(wěn)，但水翼五體船的垂蕩響應(yīng)因子要略高于五體船。在垂蕩運動時，同縱搖運動相似，水翼五體船后水翼橫向水平，增加其垂向運動的運動阻尼；前水翼除了可以增加運動阻尼意外，通過升沉?xí)r增減出水部分還使得所提供升力改變從而增加升沉回復(fù)力。因而總體上水翼五體船垂蕩性能總體要優(yōu)于五體船。

3.3 橫搖分析

引入一個標(biāo)準(zhǔn)化變量，橫搖響應(yīng)因子RP。

(13)

式中：R為橫搖角度幅值，(°)；A為有義波高，m。

橫搖主要受橫浪和斜浪的橫向分量影響，因此對于橫搖這里僅對斜浪和橫浪2種情況進(jìn)行對比分析。在斜浪和橫浪兩種浪向下，有無水翼兩種五體船橫搖狀態(tài)見圖4。

根據(jù)圖4可得，在斜浪和橫浪兩種浪向下，低頻波段加裝水翼的水翼五體船橫搖響應(yīng)因子要高于無水翼的五體船，但中、高頻波段水翼五體船橫搖響應(yīng)因子要明顯低于五體船。另外，在2種浪向下，水翼五體船橫搖響應(yīng)因子的最大值要遠(yuǎn)低于五體船，按橫搖響應(yīng)因子的最大值分析，斜浪情況下水翼五體船橫搖穩(wěn)定性要優(yōu)于五體船36.3%，而橫浪情況下水翼五體船橫搖穩(wěn)定性要優(yōu)于五體船50.3%。在橫搖運動時，水翼五體船前后水翼相當(dāng)于4個長減搖鰭，大大增加了橫搖回復(fù)力矩，同時前水翼左右兩側(cè)進(jìn)出水配合改變左右兩側(cè)的升力，進(jìn)一步增大橫搖回復(fù)力矩。

4 結(jié)論

通過在五體船主船體與側(cè)體間加裝水翼，創(chuàng)新性提出水翼五體船的概念，并借助數(shù)值模擬對水翼五體船耐波性進(jìn)行比較分析。在小搖蕩情況下，主船體與側(cè)體間水翼通過姿態(tài)變化，可以增大搖蕩回復(fù)力矩，彌補(bǔ)常規(guī)五體船前側(cè)體出水帶來的耐波性損失。水翼五體船借助其水翼升力作用及搖蕩阻尼作用，縱搖穩(wěn)定性、垂蕩穩(wěn)定性及橫搖穩(wěn)定性較常規(guī)五體船得到顯著提高。

另外，分析發(fā)現(xiàn)水翼五體船的縱搖響應(yīng)因子及橫搖響應(yīng)因子隨波頻變化異于常規(guī)五體船，出現(xiàn)2次明顯波峰，這可能由于水翼諧振頻率異于船體固有頻率所致，對此將展開進(jìn)一步研究。

[1] 王澤鴻,柳存根,周翀劍,等.應(yīng)用五體船型降低集裝箱船EEDI[J].船舶工程,2014(4):3.

[2] 熊文海,毛筱菲,李毓江.船舶耐波性衡準(zhǔn)及其評價方法淺析[J].船海工程,2007,36(4):42-45.

[3] GEE N, DUDSON E, STEIGER H. The pentamaran-A new hull concept for fast freight and car ferry applications[C].The 13th Fast Ferry International Conference. Singapore,1997.

[4] DUDSON E, GEE N. Optimisation of the seakeeping and performance of a 40 knot Pentamaran container vessel[C].Proceedings of the 6th International Conference on Fast Sea Transportation FAST2001, I.2001,225:233.

[5] 賀俊松.五體船型的性能與載荷研究[D].上海:上海交通大學(xué),2007.

[6] 徐敏,張世聯(lián).三體船/五體船的運動及波浪載荷[J].大連海事大學(xué)學(xué)報(自然科學(xué)版),2011,37(3):1-4.

[7] CALKINS D E. HYCAT: hybrid hydrofoil catamaran concept[J]. Ocean engineering,1984,11(1):1-21.

[8] CALKINS D E. Development of HYCAT[J]. AIAA paper,1986:86-2381.

[9] LIU S K. Numerical simulation of large amplitude ship motions and applications to ship design and safe operation[D]. Ph. D Thesis, National Technical University of Athens School,2002.

[10] KATAM V, LEBEAU R P, JACOB J D. Experimental and computational investigation of a modified NACA 4415 in low-Re flows[C].22nd Applied Aerodynamics Conference and Exhibit,2004:4972.

[11] 曲立群,汪建文,朱德臣,等.NACA4415翼型失速特性的二維數(shù)值研究[J].內(nèi)蒙古大學(xué)學(xué)報(自然科學(xué)版),2009,28(1):48-51.

[12] 羅少澤,馬寧,顧解忡,等.五體船阻力與耐波性協(xié)同優(yōu)化方法研究[C].2013年船舶水動力學(xué)學(xué)術(shù)會議論文集.西安：中國造船工程學(xué)會船舶力學(xué)學(xué)術(shù)委員會,2013.

[13] 毛筱菲,唐衛(wèi)軍,熊云峰,等.艦船耐波性設(shè)計及其評價指標(biāo)的計算[J].中國艦船研究,2006(1):41-45.

(上海交通大學(xué) a.高新船舶與深海開發(fā)裝備協(xié)同創(chuàng)新中心；b.船舶與海洋工程國家重點實驗室，上海 200240)

Numerical Research on Seakeeping of Hydrofoil-Pentamaran

KUANG Pei-qin, LIU Cun-gen, WANG Xue-feng, ZHOU Chong-jian

(a.Collaborative Innovation Center for Advanced Ship and Deep-Sea Exploration;b.State Key Laboratory of Ocean Engineering, Shanghai Jiaotong University, Shanghai 200240, China)

To improve the seakeeping performance of the pentamaran in small waves, a new pentamaran with four hydrofoil between main hull and accessory hulls was put forward. Its seakeeping was analyzed by HydroStar software and contrasted with normal pentamaran. It was found from the analysis that the amplitude of pitch, roll and heave of the pentamaran can be reduced by the hydrofoil between main hull and accessory hull.

pentamaran; hydrofoil; pitch; heave; roll stability

10.3963/j.issn.1671-7953.2017.01.009

2016-05-07

國家重點基礎(chǔ)研究發(fā)展計劃(2013CB036103)

匡培欽(1992—)，男，碩士生研究方向:五體船性能優(yōu)化

U661.32

A

1671-7953(2017)01-0037-04

修回日期：2016-06-13