黃衛(wèi)剛 姜治芳 邱遼原

中國(guó)艦船研究設(shè)計(jì)中心，湖北武漢 430064

1 引 言

球鼻艏是船舶設(shè)計(jì)的重要方面，設(shè)計(jì)優(yōu)良的球鼻艏具有較好的減阻效果，能有效改善船舶運(yùn)營(yíng)的經(jīng)濟(jì)性。在船舶設(shè)計(jì)領(lǐng)域，常見(jiàn)的球鼻艏形式主要包括水滴形球鼻艏和上翹球鼻艏［1］等。在國(guó)內(nèi)外現(xiàn)有的一些大型艦船上，出于布置需求，安裝有水滴形球鼻艏，但這類球鼻艏的減阻效果不佳，有時(shí)反而會(huì)使船舶阻力增加［2-3］。高速船加裝新型撞角節(jié)能球鼻艏是近年來(lái)國(guó)內(nèi)外的前沿研究課題。這類球鼻艏具有近水面、大前伸的外形特征，既能滿足類似于水滴形球鼻艏的設(shè)備布置需求，又具有良好的減阻效果，目前已得到廣泛應(yīng)用，例如，美海軍的LHA-6兩棲攻擊艦和韓國(guó)海軍的“獨(dú)島”號(hào)兩棲攻擊艦等均是采用撞角球鼻艏。

為確定撞角球鼻艏形式及參數(shù)變化對(duì)船舶阻力的影響，研究了球鼻艏形式和參數(shù)變化對(duì)阻力的影響。首先針對(duì)3個(gè)形式各異的球鼻艏型線方案開(kāi)展了CFD計(jì)算和模型試驗(yàn)，結(jié)果表明，撞角球鼻艏具有較好的減阻效果。然后設(shè)計(jì)了8個(gè)主要參數(shù)各異的撞角球鼻艏型線方案，采用CFD技術(shù)分析了撞角球鼻艏的參數(shù)變化對(duì)興波阻力的影響，并提出了一些可用于指導(dǎo)球鼻艏設(shè)計(jì)的結(jié)論。

2 球鼻艏形式對(duì)阻力影響的分析

2.1 球鼻艏形式分析

本文主要分析水滴形球鼻艏、上翹球鼻艏和撞角球鼻艏這3種球鼻艏形式。

水滴形球鼻艏主要是為了滿足布置需求，其減阻效果并不明顯。水滴形球鼻艏的外形類似于水滴（圖1），其首端與首垂線之間的距離較小，垂向中心的高度與設(shè)計(jì)吃水間的比值也較小。

上翹球鼻艏是一種能取得較好減阻效果的球鼻艏（圖2），其特點(diǎn)為：

1）球鼻艏首端與首垂線之間的距離大于2%設(shè)計(jì)水線長(zhǎng)；

2）球鼻艏首端存在一定的上翹，首端離水面距離較近；

3）球鼻艏各剖面的垂向中心由首垂線起逐漸升高。

上翹球鼻艏是通過(guò)球鼻艏興波與主船體興波之間的相互干擾來(lái)減小興波阻力。這種球鼻艏在傅汝德數(shù)Fn較高時(shí)的減阻效果比較明顯。

撞角球鼻艏是一種正逐漸得到應(yīng)用的球鼻艏形式（圖3），不僅能滿足布置需求，還具有明顯的減阻效果。其特點(diǎn)為：

1）球鼻艏首端與首垂線之間的距離大于2%設(shè)計(jì)水線長(zhǎng)；

2）球鼻艏首端的縱剖線和橫剖線均類似于橢圓；

3）球鼻艏垂向中心線平行于基線。

撞角球鼻艏的減阻原理類似于上翹球鼻艏。當(dāng)Fn較低時(shí)，撞角球鼻艏的減阻效果要優(yōu)于上翹球鼻艏；當(dāng)Fn較高時(shí)，撞角球鼻艏的減阻效果則稍劣于上翹球鼻艏。

2.2 球鼻艏形式對(duì)阻力影響的分析

2.2.1 不同形式球鼻艏型線方案設(shè)計(jì)

針對(duì)方形系數(shù)接近0.6、設(shè)計(jì)Fn接近0.3的單體船方案，設(shè)計(jì)出3套球鼻艏型線方案，如表1所示。

表1 不同形式球鼻艏方案的主要尺度Tab.1 Main dimensions of different bulbous bow schemes

2.2.2 興波阻力計(jì)算及結(jié)果

興波阻力計(jì)算軟件適用于船舶各種水動(dòng)力特性研究，并給出了波浪模式、壓力分布等船體特征參數(shù)。因這一領(lǐng)域的技術(shù)比較成熟［4-7］，故本文將直接采用興波阻力軟件進(jìn)行數(shù)值計(jì)算分析。計(jì)算興波阻力主要采用勢(shì)流計(jì)算法，即采用基于Dawson方法的二階面源法計(jì)算勢(shì)流，通過(guò)迭代法，能很好地滿足自由液面的邊界條件。Dawson方法其實(shí)就是在疊模繞流的基礎(chǔ)上，采用慢船條件下的線性化理論。該方法比疊模線性法還要深入一步，在獲得線性解后并沒(méi)有結(jié)束計(jì)算，而是通過(guò)迭代不斷重復(fù)該過(guò)程，直至迭代收斂，就可以得到完全非線性解。

本文采用興波阻力計(jì)算軟件計(jì)算興波阻力系數(shù)Cw，為使各方案之間的計(jì)算結(jié)果具備可比性，采取下列措施控制網(wǎng)格劃分對(duì)計(jì)算產(chǎn)生的影響：

1）流體域計(jì)算范圍以及網(wǎng)格設(shè)置參數(shù)需保持一致；

2）球鼻艏影響區(qū)域的網(wǎng)格劃分參數(shù)需保持一致。

球鼻艏引起的濕表面積差異小于1%，此時(shí)，各方案濕表面積差異對(duì)阻力的影響可忽略不計(jì)。因此，本文主要比較各方案的興波阻力系數(shù)。4個(gè)球鼻艏方案的阻力計(jì)算結(jié)果如表2所示。

表2 球鼻艏多方案興波阻力計(jì)算結(jié)果Tab.2 Wave coefficient calculation results of different bulbous bow schemes

在計(jì)算中，確定了各方案沿船體表面的興波分布情況，如圖4～圖7所示。

根據(jù)表2確定的球鼻艏形式對(duì)興波阻力系數(shù)的影響趨勢(shì)如圖8所示。

由圖8 可知，當(dāng) Fn＝0.27～0.29 時(shí)，采用水滴形球鼻艏方案的興波阻力最大，但隨著Fn的增加，水滴形球鼻艏的興波阻力系數(shù)逐漸接近于無(wú)球鼻艏方案的興波阻力系數(shù)，這表明當(dāng)Fn足夠大時(shí)，水滴形球鼻艏會(huì)具備一定的減阻效果；當(dāng)Fn＝0.27～0.29 時(shí)，無(wú)球鼻艏方案的興波阻力系數(shù)介于水滴形球鼻艏、撞角球鼻艏與上翹球鼻艏的興波阻力系數(shù)之間；當(dāng) Fn=0.27～0.29 時(shí)，采用撞角球鼻艏的方案A和采用上翹球鼻艏的方案B均具有良好的減阻效果。當(dāng)Fn=0.29時(shí)，與無(wú)球鼻艏方案相比，撞角球鼻艏的興波阻力系數(shù)約小7%。撞角球鼻艏與上翹球鼻艏的興波阻力系數(shù)的變化趨勢(shì)基本一致。

2.3 球鼻艏形式對(duì)阻力的驗(yàn)證性模型試驗(yàn)

為進(jìn)一步驗(yàn)證CFD計(jì)算結(jié)果的可信度，并確定撞角球鼻艏的減阻效果，開(kāi)展了球鼻艏形式對(duì)阻力的驗(yàn)證性模型試驗(yàn)（圖9）。

由圖9可知，當(dāng)Fn=0.28時(shí)，撞角球鼻艏的減阻效果最好，上翹球鼻艏的減阻效果次之，兩者之間的差異不大；當(dāng)Fn=0.28時(shí)，撞角球鼻艏和上翹球鼻艏的剩余阻力系數(shù)比水滴形球鼻艏的剩余阻力系數(shù)約小9%，比無(wú)球鼻艏的剩余阻力系數(shù)約小20%。模型試驗(yàn)結(jié)果表明，撞角球鼻艏具有良好的減阻效果。

由圖8、圖9可知，CFD流場(chǎng)計(jì)算確定的無(wú)球鼻艏方案的興波阻力比水滴形球鼻艏方案的小，模型試驗(yàn)確定的無(wú)球鼻艏方案的剩余阻力比水滴球鼻艏方案的大。出現(xiàn)這種情況的原因是：剩余阻力包括興波阻力和粘壓阻力兩部分，無(wú)球鼻艏方案因首尾差異較大而導(dǎo)致粘壓阻力較大，水滴形球鼻艏方案的粘壓阻力又較小，因此在當(dāng)前Fn（Fn=0.27～0.29）下，粘壓阻力在剩余阻力系數(shù)中所占的比例便較大，此時(shí)就會(huì)出現(xiàn)興波阻力較小而剩余阻力較大的情況。

3 撞角球鼻艏參數(shù)變化對(duì)阻力影響的分析

鑒于撞角球鼻艏具有良好的減阻效果，為進(jìn)一步掌握球鼻艏參數(shù)變化對(duì)阻力的影響，開(kāi)展了撞角球鼻艏的多方案設(shè)計(jì)，并用CFD方法進(jìn)行了計(jì)算，初步確定了撞角球鼻艏參數(shù)變化對(duì)阻力的影響。

3.1 球鼻艏參數(shù)分析

球鼻艏的幾何特征形狀可以通過(guò)以下參數(shù)進(jìn)行描述：

1）縱向長(zhǎng)度比：即球鼻艏最前端至首柱的距離與設(shè)計(jì)水線長(zhǎng)之間的比值；

2）相對(duì)浸深：球鼻艏中心或球鼻艏最前點(diǎn)或者最大寬度處距靜水面的距離與吃水之間的比值；

3）寬度比：首柱處球鼻艏最大橫剖面寬度與設(shè)計(jì)水線寬之間的比值；

4）面積比：首柱處球鼻艏橫剖面面積與最大橫剖面面積的比值；

5）體積比：球鼻艏增加的排水體積與設(shè)計(jì)排水體積之間的比值。

撞角球鼻艏的縱剖面和橫剖面均類似于橢圓，此時(shí)，設(shè)計(jì)吃水、縱向長(zhǎng)度比和寬度比一旦確定，其相對(duì)浸深、面積比和體積比也就得以確定。本文主要研究球鼻艏的縱向長(zhǎng)度比與寬度比對(duì)阻力的影響。

3.2 球鼻艏寬度比對(duì)阻力的影響分析

以撞角球鼻艏方案A為基礎(chǔ)，設(shè)計(jì)4個(gè)球鼻艏寬度比變化的方案，各方案的具體參數(shù)如表3所示。

表3 球鼻艏寬度比方案Tab.3 Different ratio of bulbous bow breath and waterline breath for bulbous bow scheme

球鼻艏寬度比的范圍為 5.0%～10.2%，制定該范圍的原因是：

1）球鼻艏的寬度比增加，會(huì)出現(xiàn)在拋錨時(shí)錨會(huì)撞到球鼻艏的情況。因此，從工程實(shí)現(xiàn)的可行性角度考慮，球鼻艏寬度比的上限應(yīng)取10.2%；

2）球鼻艏寬度比的下限取5%，這是因?yàn)榇藭r(shí)球鼻艏的體積已經(jīng)很小，球鼻艏的減阻效果會(huì)較差。

針對(duì)這5個(gè)球鼻艏型線方案，開(kāi)展了CFD計(jì)算，其結(jié)果如表4、圖10所示。

表4 球鼻艏多方案CFD阻力計(jì)算結(jié)果Tab.4 The resistance calculation results of different bulbous bow schemes

由表4和圖10可知，球鼻艏寬度比對(duì)興波阻力系數(shù)存在著一定的影響，當(dāng)寬度比為7.4%時(shí)，其興波阻力系數(shù)最??；當(dāng)寬度比為8.4%時(shí)，興波阻力系數(shù)最大。當(dāng)球鼻艏的寬度比變化范圍為5%～10.2%，F(xiàn)n＝0.27～0.29 時(shí)，球鼻艏寬度比的變化將導(dǎo)致興波阻力系數(shù)的變化范圍約小于5%。球鼻艏寬度比對(duì)興波阻力系數(shù)的影響趨勢(shì)不明顯，出現(xiàn)這種情況的原因是，球鼻艏減阻效果在寬度方向存在非線性情況。

3.3 球鼻艏縱向長(zhǎng)度比對(duì)阻力的影響分析

以撞角球鼻艏方案E為基礎(chǔ)，設(shè)計(jì)了4套球鼻艏縱向長(zhǎng)度比變化的型線方案，各方案的具體參數(shù)如表5所示。

球鼻艏縱向長(zhǎng)度比的范圍為2%～5%，確定該范圍的原因?yàn)椋焊鶕?jù)設(shè)計(jì)經(jīng)驗(yàn)，當(dāng)球鼻艏縱向長(zhǎng)度比為2%時(shí)，基本沒(méi)有減阻效果，因此球鼻艏縱向長(zhǎng)度比的下限取2%的設(shè)計(jì)水線長(zhǎng)。球鼻艏若太長(zhǎng)，會(huì)導(dǎo)致建造成本增加，影響船舶的經(jīng)濟(jì)性，考慮到工程的可行性，球鼻艏縱向長(zhǎng)度比的上限取5%的設(shè)計(jì)水線長(zhǎng)。

表5 球鼻艏長(zhǎng)度比方案Tab.5 Different ratio of bulbous bow length to waterline length for bulbous bow schemes

針對(duì)這5個(gè)球鼻艏型線方案，開(kāi)展了CFD計(jì)算，其結(jié)果如表6、圖11所示。

表6 球鼻艏多方案CFD阻力計(jì)算結(jié)果Tab.6 The resistance calculation results of different bulbous bow schemes

由表6及圖11可知，球鼻艏縱向長(zhǎng)度比對(duì)阻力的影響較明顯，本文研究的主船體方案具有可取得較好減阻效果的最佳縱向長(zhǎng)度比。當(dāng)前方案中減阻效果最好時(shí)的縱向長(zhǎng)度比為3%。在Fn=0.27～0.29，球鼻艏縱向長(zhǎng)度比變化范圍介于 2%～5%之間時(shí)，當(dāng)縱向長(zhǎng)度比小于3%時(shí)，球鼻艏縱向長(zhǎng)度比變化會(huì)導(dǎo)致興波阻力系數(shù)的變化范圍約小于7%；當(dāng)縱向長(zhǎng)度比大于3%時(shí)，球鼻艏縱向長(zhǎng)度比變化會(huì)導(dǎo)致興波阻力系數(shù)的變化范圍約小于3%。

4 結(jié) 論

本文開(kāi)展了球鼻艏形式和參數(shù)變化的多方案球鼻艏研究，得到了以下結(jié)論：

1）對(duì)于方形系數(shù)約為 0.6、 設(shè)計(jì)Fn接近0.3的常規(guī)單體船方案而言，撞角球鼻艏和上翹形球鼻艏均具有明顯的減阻效果。當(dāng)Fn=0.28時(shí)，撞角球鼻艏的剩余阻力系數(shù)比水滴形球鼻艏的剩余阻力系數(shù)約小9%，比無(wú)球鼻艏的剩余阻力系數(shù)約小20%。

2）CFD計(jì)算結(jié)果表明，撞角球鼻艏的寬度比和縱向長(zhǎng)度比均對(duì)興波阻力系數(shù)具有一定的影響，并且縱向長(zhǎng)度比對(duì)興波阻力系數(shù)的影響更明顯。

［1］葉明，吳正廉.大伸出上翹型球鼻在一艘小方形系數(shù)船上的應(yīng)用［J］.船舶，2003（4）：5-7.YE M，WU Z L.Application of large-protruded upwardbending bulbous bow in a small block coefficient ship［J］.Ship and Boat，2003（4）：5-7.

［2］吳曉光，徐青，吳宏敏，等.國(guó)外水面艦船50年［M］.北京：兵器工業(yè)出版社，2006.

［3］張廣思.水面艦艇球鼻艏設(shè)計(jì)中的幾個(gè)問(wèn)題［J］.艦船科學(xué)技術(shù)，1983（12）：1-11.ZHANG G S.Some problems in bulbous bow design for surface warship［J］.Ship Science and Technology，1983（12）：1-11.

［4］王中，盧曉平.基于細(xì)分曲面的興波阻力計(jì)算方法［J］.海軍工程大學(xué)學(xué)報(bào)，2006，18（增刊）：20-25.WANG Z，LU X P.A calculation method for ship wave making resistance based on subdivision surface theory［J］.Journal of Naval University of Engineering，2006，18（Supplement）：20-25.

［5］王中，盧曉平.水面艦船加裝減阻節(jié)能球鼻艏研究［J］.水動(dòng)力學(xué)研究與進(jìn)展（A 輯），2006（6）：789-795.WANG Z，LU X P.Research on high speed surface warship with bulbous bow［J］.Journal of Hydrodynamics（Ser.A），2006（6）：789-795.

［6］DATLA R，KIM H Y，STEBE J R.Evaluation of a CFD program AEGIR （trademark）for bare hull resistance and seakeeping prediction capability［R］.USA Maryland：Naval Surface Warfare Center，2009.

［7］BOHYUN Y，DAVID W.A ray theory for nonlinear ship waves and wave resistance［M］.USA Maryland：David Taylor Naval Ship Research and Development Center，1981.