王 超，何 苗，王 偉，熊 鷹，黃 勝

（1.海軍工程大學(xué)船舶與動力學(xué)院，湖北武漢430033；2.哈爾濱工程大學(xué)船舶工程學(xué)院，黑龍江哈爾濱150001）

艏艉線型優(yōu)化對船舶阻力性能的影響＊

王 超1，2，何 苗2，王 偉2，熊 鷹1，黃 勝2

（1.海軍工程大學(xué)船舶與動力學(xué)院，湖北武漢430033；2.哈爾濱工程大學(xué)船舶工程學(xué)院，黑龍江哈爾濱150001）

為提高母型船阻力性能，以船體阻力性能為優(yōu)化對象，基于改造母型船法，研究船舶球鼻艏以及船尾線型的改變對船舶阻力性能的影響。采用高度集成化的Tribon系統(tǒng)、可視化繪圖軟件Auto CAD及CFD（Computational Fluid Dynamics）通用前處理軟件ICEM聯(lián)合建模的方法來建立船體模型。通過模擬計算結(jié)果與實驗值的對比分析，驗證CFD技術(shù)在船舶阻力性能預(yù)報中的合理性和有效性。通過對比3種不同球鼻艏時的船體阻力得知：從阻力性能方面考慮，對于低速豐滿型船舶選用普通型球鼻艏以及中高速船舶采用上翹型球鼻艏均可以獲取較好地減阻效果。同時比較不同航速下尾部線型對船體總阻力的影響表明，選優(yōu)后的方形尾在相同的航速下阻力低、消耗的功率小、形狀效應(yīng)小、黏壓阻力和摩擦阻力也相對較小。

球鼻艏；阻力性能；傅汝德數(shù)；CFD技術(shù)；數(shù)值模擬

隨著當(dāng)今世界能源緊張，燃油成本在船舶營運總成本中的比重也不斷提升，開發(fā)性能優(yōu)良的船型，節(jié)能降耗已成航運界和造船界共同追求的目標(biāo)。在船舶性能研究領(lǐng)域，降低油耗，提升船舶的快速性指標(biāo)，已成為新船建造必做的研發(fā)內(nèi)容［1］。船舶線型的優(yōu)劣是衡量船舶快速性技術(shù)水平的標(biāo)志性因素之一。在現(xiàn)代艦船設(shè)計中，球鼻艏被廣泛地應(yīng)用在各種軍事和民用船舶中。對于民用船舶來說，球艏產(chǎn)生的波與船首前波相互疊加可以形成有利干擾，能夠有效地減小興波阻力，有利于提高航速以達到更好的經(jīng)濟性能；對于軍用艦船，球鼻艏內(nèi)還可以放置聲納等設(shè)備，實現(xiàn)特殊的功能［2－5］。球鼻艏有常規(guī)型、上翹型、埋首型等多種形式，不同形式的球鼻艏對于船舶性能的影響也各不相同［6－7］。從有利于阻力性能的角度來進行球鼻艏的選型，以獲得更加優(yōu)良的局部型線，為船舶節(jié)能提供一種行之有效的方法。

隨著計算機技術(shù)的快速發(fā)展，為船型優(yōu)化設(shè)計提供了高效快捷的手段。船型優(yōu)化包括船體型線優(yōu)化，船體型線優(yōu)化通?？梢詰?yīng)用數(shù)學(xué)方法對型線進行光順，電子計算機輔助產(chǎn)生型線更是得到了廣泛的應(yīng)用［8］。近年來CFD技術(shù)的發(fā)展迅猛，使得船舶黏性繞流場的數(shù)值模擬也有了長足的進步，對于單體船流場，不僅可以得到船體表面波形以及尾部伴流特征，在阻力預(yù)報上也與試驗吻合較好，為船體線型的優(yōu)化和評價提供了良好手段。

以船體阻力性能為型線優(yōu)化對象時，只能在某些約束條件下完成。本文在保證母型船主要參數(shù)不變的情況下，基于改造母型船法，借助Tribon－Lines模塊進行船舶線型的生成和光順。采用CFD技術(shù)對船舶阻力較為敏感的球艏及尾部線型進行分析研究，通過船舶快速性及流場特性預(yù)報，比較不同球鼻艏及尾部線型方案的阻力性能優(yōu)劣。同時文中還對不同傅汝德數(shù)下船體周圍自由液面三維波形圖和波高等值線進行了比較分析。

1 艏艉的選型

1.1 球鼻艏線型

選取3條與母船型相似的帶有不同球鼻艏的船型進行黏性繞流場的數(shù)值模擬，討論不同形式的艏部線型對船體阻力的影響。為了便于說明，3種新船型記為方案1、2和3。方案1具有與母型船同樣的埋首型球鼻艏，只是尺度有所不同（具體差異見表1）；方案2為普通型球鼻艏；方案3為上翹型球鼻艏，球艏具體形狀見圖1。為了使計算結(jié)果具有可比性，將3條新船主尺度經(jīng)過仿似變換變成與母船型相同的尺寸。

表1 不同船型的球艏幾何參數(shù)Table 1 The parameters of different bulbous bows

上表中球艏的主要幾何參數(shù)［7］定義如下：

突出量（長度系數(shù)）：CLpr＝Lpr／Lpp（Lpr—球艏伸出首垂線的長度，Lpp—垂線間長）；

高度系數(shù)：CZb＝Zb／Tft（Zb—球艏最前緣點距基線的高度，Tft—首吃水）；

橫剖面面積系數(shù)：CS＝S／SM（S—首垂線處橫剖面面積，SM—船中橫剖面面積）；

縱剖面面積系數(shù)：CSL＝SL／SM（SL球艏縱剖面面積；SM—船中橫剖面面積）；

圖1 各方案船型球鼻艏橫剖面形狀Fig.1 The bulbous bows section shapes of different projects

1.2 尾部線型

母船型尾部為巡洋艦艉，方案船型為兩形狀不同的方形艉。為討論方便，現(xiàn)將母型船及改造后的船型分別編號為RCSM－Ⅰ、RCSM－Ⅱ、RCSM－Ⅲ。在進行尾部線型改進時，假定參考船型與母型船的浮心縱向位置Xb完全相同，將參考船型的總尺寸縮尺比為母型船的主尺度，然后取其后體型值進行方案船型的組合。在進行船艉型線生成時特別要注意過渡區(qū)的型線光順，其過渡質(zhì)量用Tribon－Lines的3向光順模塊來確定。圖2給出了母型船的橫剖線圖以及改型后的方案船型的橫剖線圖。RCSM系列船型的濕表面積變化、排水體積變化、方形系數(shù)Cb以及棱形系數(shù)Cp變化見表2。

圖2 艉部型線Fig.2 The section lines of stern

表2 RCSM系列船型參數(shù)Table 2 The parameters of RCSM ship hulls

2 湍流模型的選取

SST（Shear Stress Transport）湍流模型考慮了湍流剪應(yīng)力的輸運特性，能夠準(zhǔn)確地預(yù)報由于逆壓梯度導(dǎo)致的流動分離點和分離區(qū)域，因而在帶分離的復(fù)雜湍流的預(yù)報研究中能發(fā)揮很大的優(yōu)勢［9］。由于該湍流模型是在結(jié)合模化后的k－ω湍流模型和k－ω湍流模型的基礎(chǔ)上建立起來的，并且在湍流黏性系數(shù)中考慮了輸運效應(yīng)，因而比k－ω二方程湍流模型更加先進［7］。

湍動能k的輸運方程：

湍流耗散ω的輸運方程：

湍流黏性系數(shù)：

混合函數(shù)（Blending function）F1和F2有如下定義：

3 計算模型的建立

3.1 船體與計算域生成

采用高度集成化的Tribon系統(tǒng)、可視化繪圖軟件Auto CAD和CFD通用前處理軟件ICEM聯(lián)合建模的方法來建立船體模型。其中采用Tribon系統(tǒng)的Lines模塊來完成型值細化的內(nèi)容，界面見圖3。在Lines模塊中導(dǎo)出.dxf格式文件（三維視圖），在可視化繪圖軟件AutoCAD里做一下簡單處理：為了避免導(dǎo)進ICEM后圖形失真，可在折角線處打斷；并且為了有利于在ICEM里建模方便，在盡量滿足真實船體曲面的前提下，可以將一些不必要的線條刪掉。最后就可以在通用CFD前處理軟件ICEM里生成船體曲面了，最終的實體模型見圖4。

圖3 Lines模塊系統(tǒng)運行界面Fig.3 The run interface of module system

圖4 船體曲面Fig.4 The ship hull

計算域選取為1個橢圓形［10］，縱向從船艏向前延伸1個船長，船尾向后延伸3個船長，橫向從中縱剖面向兩側(cè)各延伸1.5個船長，靜水面向下延伸1個船長。由于流動對稱，計算域取整個流域的一半，計算域的示意圖見圖5。其中入口處給定均勻來流速度V及湍流參數(shù)，出口采用開放式邊界，對稱面使用對稱邊界條件，計算域頂部及遠方邊界為滑移壁面，船體表面為不可滑移壁面。

圖5 計算域分塊示意圖Fig.5 The three dimension chart of calculative region

3.2 網(wǎng)格劃分

高質(zhì)量的網(wǎng)格是實現(xiàn)數(shù)值模擬成功的首要條件。結(jié)構(gòu)化網(wǎng)格的特點是每一節(jié)點與其相鄰點之間的聯(lián)接關(guān)系固定不變且隱含在所生成的網(wǎng)格中，因而不必專門設(shè)置數(shù)據(jù)去確認節(jié)點和相鄰點之間的聯(lián)系，可以減少一些常用差分格式的數(shù)值擴散誤差；同時在計算過程中收斂性快，穩(wěn)定性好，可以較準(zhǔn)確地體現(xiàn)固壁對流場的影響。非結(jié)構(gòu)化網(wǎng)格適應(yīng)不規(guī)則區(qū)域性強，局部和全局分解、加密和整合性好，能較好地處理邊界和控制網(wǎng)格單元的大小、形狀及網(wǎng)格的位置［9］。本文網(wǎng)格劃分是借助CFD通用前處理軟件ICEM來實現(xiàn)的，采用分塊結(jié)構(gòu)化網(wǎng)格，艏艉部均采取C型網(wǎng)格，其他部分采取H型網(wǎng)格。為了更好的捕捉流動特征，網(wǎng)格在船艏、船尾和靜水面附近進行了加密。整個計算域的網(wǎng)格分布見圖6。

圖6計算域網(wǎng)格劃分Fig.6 The mesh distribution of ship calculation domain

4 數(shù)值計算結(jié)果分析

本文在考慮黏性理論的范疇內(nèi)，將目前較流行的CFD技術(shù)應(yīng)用于船體的阻力性能研究中。在考慮自由液面效應(yīng)的前提下，模擬在不同航速下船體的黏性流場，獲取了不同球艏及尾部線型時船體的阻力以及自由液面的變化情況等。由于研究對象屬于排水船型，運動過程中不可避免地會受到興波阻力的作用，所以需要考慮到自由液面問題［11］。計算時采用基于VOF（Volume of Fluid）法的均勻多相流模型來模擬自由液面。

4.1 阻力計算結(jié)果驗證

圖7給出了算例船在不同傅汝德數(shù)下船體阻力變化結(jié)果，通過與實驗值［12］對比可以發(fā)現(xiàn)二者存在一定的誤差，并且隨著航速的增加差距也在加大。誤差的產(chǎn)生主要應(yīng)由試驗誤差和數(shù)值模擬誤差2部分組成，但從圖7中看出總體來說數(shù)值計算和試驗結(jié)果變化趨勢比較吻合，可以滿足工程運用需要。進而認為采用CFD技術(shù)對考慮自由液面的船體阻力進行預(yù)報可得比較合理的結(jié)果，CFD技術(shù)將成為船舶阻力性能預(yù)報中的1種快捷且有效的方法。

圖7 阻力計算值與試驗值比較Fig.7 Comparison of resistance between calculation and experimentation

4.2 不同球鼻艏時的船體阻力

圖8給出了在不同傅汝德數(shù)（Fr）時具有不同類型球鼻艏船體的阻力變化曲線，從中可以看出不同航速情況下方案1（埋首型球鼻艏）產(chǎn)生的船體阻力比另外方案產(chǎn)生的阻力都大；在低速段（Fr＜0.23）時，方案3（上翹型大球艏）船體阻力值介于方案1和2（普通型小球艏）之間；而當(dāng)Fr＞0.23時，方案3阻力值最小，方案1值最大。從而以阻力性能方面考慮，對于低速豐滿型（Fr＜0.22）船舶選用方案2以及中高速（Fr＞0.22）船舶采用方案3均可以獲取較好地減阻效果。

圖8 船體阻力曲線Fig.8 The curve of ship resistance

4.3 不同尾部線型時的船體阻力

3個船型的數(shù)值模型均參照船舶水池試驗方案數(shù)值模型的相關(guān)設(shè)定來建立。數(shù)值計算模擬了這3個方案分別在6個航速下的黏性自由面流動。圖9給出了母船型阻力計算值與試驗值比較以及方案船型的阻力計算值。由圖9可以看出：（1）計算所得的RCSM－Ⅰ船型靜水阻力曲線和試驗所得曲線發(fā)展趨勢基本上一致，大體上能較好地反映單體船模在靜水中航行時的阻力特性。（2）計算得到的總阻力與模型阻力試驗得到的總阻力在航速較低時比較接近，但隨著航速的增加，總阻力系數(shù)計算值與試驗值產(chǎn)生偏差；從整個模擬結(jié)果看，模擬計算值要比試驗值小。（3）從計算所得的3種方案船型阻力值可以看出，RCSM－Ⅱ的方形尾的阻力值明顯低于RCSM－Ⅰ和RCSM－Ⅲ船型阻力值。

圖9 RCSM船型阻力曲線Fig.9 Resistance curves of RCSM ship hulls

4.4 球鼻艏對自由液面的影響

圖10 自由表面三維波形圖和波高等值線（Fr＝0.173 6）Fig.10 The mesh distribution of ship calculation domain

圖10給出了船體周圍自由液面三維波形圖和波高等值線。通過興波凱爾文角大小及橫波寬度可以看出，在相同航速下，低速時：上翹型大球艏比埋首型和普通型小球艏興起的波浪稍大些。同時計算過程中發(fā)現(xiàn)同一條船橫向比較時，隨著航速的增加艏艉部的興波越來越明顯，并且有船舶首肩波波峰后移的現(xiàn)象。

4.5 尾流場分析

取Fr＝0.231 5時船舶的尾流場為分析對象，數(shù)值模擬船舶的自由表面波形和尾流場波形見圖11和12所示。

圖11 船尾波高分布圖Fig.11 Distribution of RCSM－Ⅱwake wave Contours

圖12 尾流場分布Fig.12 Distribution of wake flow

由圖11可以看出RCSM－Ⅲ船尾產(chǎn)生的興波明顯比RCSM－Ⅰ、RCSM－Ⅱ產(chǎn)生的興波大；從圖12可以看出RCSM－Ⅱ和RCSM－Ⅲ的方形艉部出現(xiàn)了明顯的“虛長度”［13］，而且計算結(jié)果顯示隨著航速的增加“虛長度”加長，這就相當(dāng)于在船體濕表面積不變的情況下增加了設(shè)計水線長，對阻力而言在摩擦阻力不增大的情況下減小了剩余阻力。對于中低速船舶尾部興波阻力的貢獻并不是很大，黏壓阻力和摩擦阻力占主要地位，RCSM－Ⅱ船型的尾部濕表面積比RCSM－Ⅰ、RCSM－Ⅲ濕面積小，形狀較規(guī)則，所以形狀效應(yīng)小，因而RCSM－Ⅱ的黏壓阻力和摩擦阻力也相對較小。

5 結(jié)論

該文采用CFD技術(shù)分析了球鼻艏及尾部線型對船體阻力性能的影響，通過對計算結(jié)果的分析，得出以下結(jié)論：

（1）CFD方法在計算船體阻力時考慮了黏性影響，由計算結(jié)果與實驗值的比較認為CFD技術(shù)對考慮自由液面的船體阻力進行預(yù)報可得出比較合理的結(jié)果。

（2）從而以阻力性能方面考慮，對于低速豐滿型（Fr＜0.22）船舶選用普通型球鼻艏以及中高速（Fr＞0.22）船舶采用上翹型球鼻艏均可以獲取較好地減阻效果。同時對同一條船橫向比較時，隨著航速的增加艏艉部的興波越來越明顯，并且有船舶艏肩波波峰后移的現(xiàn)象。

（3）從計算所得的3種不同尾部船型阻力值比較可以看出，模擬計算值要比試驗值小，且RCSM－Ⅱ的方形尾的阻力值明顯低于RCSM－Ⅰ和RCSM－Ⅲ船型阻力值。同時可以看出RCSM－Ⅱ和RCSM－Ⅲ的方形艉部出現(xiàn)了明顯的“虛長度”，而且計算結(jié)果顯示隨著航速的增加“虛長度”加長。

由本文的研究可以看出采用CFD技術(shù)輔助船型設(shè)計和性能優(yōu)化，是一種提高研究開發(fā)效率的實用途徑，有著重要的工程實用價值。但作者論文中對球鼻艏和尾部線型是分開研究的，接下來將系統(tǒng)地進行二者組合的船型設(shè)計和性能優(yōu)化的研究。

［1］ 陳霞萍，陳昌運.基于CFD與EFD的船體線型優(yōu)化設(shè)計研究［C］∥.《船舶力學(xué)》創(chuàng)刊十周年紀(jì)念學(xué)術(shù)會議論文集，無錫：［s.l.］，2007：243－250.

［2］ 王中，盧小平.水面艦船加裝減阻節(jié)能球鼻艏研究［J］.水動力學(xué)研究與進展，2006，21（6）：789－795.

［3］ 李為.護衛(wèi)艦球鼻首型線設(shè)計的探討［J］.艦船工程研究，1997（3）：19－23.

［4］ 蔡榮泉，張彥明.艦船流場計算程序系統(tǒng)及其在淺吃水肥大船球艏選型中的應(yīng)用［J］.艦船科學(xué)技術(shù)，1997，2：16－20.

［5］ 劉文璽.對高速艦船雙球首構(gòu)形優(yōu)化的初步研究［D］.哈爾濱：哈爾濱工程大學(xué)，2004.

［6］ 盛振邦，劉應(yīng)中.船舶原理［M］.上海：上海交通大學(xué)出版社，2003.

［7］ 張廣思.水面艦艇球鼻艏設(shè)計中的幾個問題［J］.艦船科學(xué)技術(shù)，1983，12：1－10.

［8］ 王言英.基于阻力性能船體型線精細優(yōu)化的CFD方法［J］.大連理工大學(xué)學(xué)報，2002，42（2）：127－131.

［9］ 王福軍.計算流體動力學(xué)分析—CFD軟件原理與應(yīng)用［M］.北京：清華大學(xué)出版社，2004：113－142.

［10］ 孫榮，吳曉光，姜治芳.三體船阻力數(shù)值計算及方案優(yōu)選［J］.中國艦船研究，2008，3（5）：21－24.

［11］ LI Tingqiu.Computation of turbulent free－surface flows around modern ships［J］.Int J Numer Meth Fluids，2003，43：407－430.

［12］ 王偉.大型艦船流場的數(shù)值計算［D］.哈爾濱：哈爾濱工程大學(xué)，2011.

［13］ 李云波.船舶阻力［M］.哈爾濱：哈爾濱工程大學(xué)出版社，2004：103－106.

Influence on the Ship Resistance Performance of the Lines Optimization of the Bow and Stern

WANG Chao1，2，HE Miao2，WANG Wei2，XIONG Ying1，HUANG Sheng2
（1.College of Naval Architecture and Marine Power，Naval University of Engineering，Wuhan 430033，China；2.College of Shipbuilding Engineering，Harbin Engineering University，Harbin 150001，China）

In order to improve the performance of the original hull resistance，with the hull resistance as optimal object，the effect of varying bulbous stem and mold lines of stern to hull resistance was studied based on the modifying parent ship method.The highly integrated Tribon system，the visual mapping software Auto CAD and the current CFD pre－treatment software ICEM were used in the establishment of the hull model.The rationality and validity of the CFD technology used in the research of resistance performance was validated through the comparison of the results of simulation and experiment.Considered from the resistance performance aspect，the low speed plentiful ships selecting the popular type bulbous bow and the medium and high speed ships selecting the goes up bulbous bow may obtain the good effect of drag reduction through the comparison of the resistance of the ships with different kind of bulbous bow.Meanwhile，by comparing hull resistance of various mold lines stern under different speed show that the optimal ship hull featuring lower resistance and the minimum power consumption at the same speed，whose shape effect is smaller，the pressure resistance and friction resistance are relatively smaller than others through the effect of stern line on the hull resistance under different velocities.

bulbous stem；resistance performance；froude number；computational fluid dynamics（CFD）；numerical simulation

U661

A

1672－5174（2012）04－075－07

國家自然科學(xué)基金項目（10702016）；中央高校基本科研業(yè)務(wù)費專項（HEUCFT1001）資助

2011－04－28；

2011－10－15

王 超（1982－），男，講師。E－mail：zhitao0213＠sina.com

責(zé)任編輯 陳呈超