基于DAWSON法小水線面雙體船模型興波阻力及浮態預報

范井峰李云波

（1.海軍駐上海地區艦艇設計研究軍事代表室 上海 200011；2.哈爾濱工程大學 哈爾濱 150001）

基于DAWSON法小水線面雙體船模型興波阻力及浮態預報

范井峰1李云波2

（1.海軍駐上海地區艦艇設計研究軍事代表室 上海 200011；2.哈爾濱工程大學 哈爾濱 150001）

針對小水線面雙體船型，基于Dawson法，通過船體局部網格快速劃分迭代求解船體模型在不同航速下的興波阻力、縱傾、升沉以及不同航態下的濕表面面積變化，并給出預報結果與模型試驗結果的對比。

DAWSON法；小水線面雙體船；興波阻力；縱傾；升沉

引 言

小水線面雙體船（SWATH）優秀的耐波性能已經得到業內共識。例如，“卡瑪林諾”號和“海鷗”號小水線面雙體船的垂向運動加速度均小于0.1重力加速度；“海鷗”號小水線面雙體船4級海況下橫搖幅值僅為常規單體船的1/4。 SWATH的耐波性能使其使用效率比單體船高得多［1］，且已廣泛用于海洋水文考察船、海上作業試驗船、交通艇等。

SWATH由于水線面積較小，航行時縱傾和升沉變化較大，通常SWATH船型均要安裝首尾穩定鰭來改善航行時的埋首現象。

本文基于Dawson方法［2-4］，通過船體局部網格［5］快速劃分迭代求解船體模型在不同航速下的興波阻力、縱傾、升沉以及不同航態下的濕表面面積變化，并給出預報結果與模型試驗結果的對比［6］。

1 計算原理

1.1 興波阻力基本理論與方法

取隨船坐標系，坐標原點位于水平面上，見下頁圖1。

圖1 坐標系示意圖

假設流體理想，流動定常，無限水深，船舶興波為微幅波。流場速度勢ψ分解為無自由面影響的重疊模速度勢?和考慮自由液面作用的擾動速度勢φ，即

速度勢均滿足無窮遠處趨于零的遠方輻射條件。考慮到?關于靜水面對稱，則

將自由面條件（5）分別在z=ζ和z=0上速度勢展開，并假設????，忽略?的非線性項，可得簡化的自由面邊界條件

式中：l為重疊模流線方向。

船體受到的興波阻力Rw、興波阻力系數Cw、升力Fz、力矩Mz表示為：

1.2 網格處理

為了在較大傅汝德數范圍內獲得較準確和光滑而真實的興波波形，將自由面網格沿縱向從船首向遠前方、船尾向后方呈輻射狀網格長度按比例遞增。橫向網格尺度按比例向外遞增，波形變化較大區域，采用較密集的均勻網格寬度。

2 計算結果與試驗結果對比

本文主要對正常排水量下的船舶阻力及航行姿態進行分析。該船模主尺度見表1。

表1 船舶主尺度

對船模及自由表面進行數值離散，船體表面離散為1 500單元，自由表面單元數為2 300單元。船體表面離散后網格模型如圖2。

圖2 小水線面雙體船網格劃分

計算求得的傅汝德數從0.1到0.45共15個速度點的興波阻力數值計算結果見表2。

表2 小水線面雙體船興波阻力及航態預報結果

圖3 -圖6給出了計算結果與模型試驗結果的對比。

圖3 計算興波阻力系數與試驗剩余阻力系數比較

圖4 計算縱傾值與試驗剩余阻力系數比較

圖5 計算升沉值與試驗數據比較

圖6 計算濕表面面積變化曲線

圖3為計算興波阻力系數與試驗剩余阻力系數數據的比較，在計算傅汝德數范圍內，該小水線面雙體船興波阻力系數出現明顯峰值和谷值，在傅汝德數0.3附近出現明顯峰值，在傅汝德數0.35附近出現明顯谷值。

圖4為計算縱傾值與試驗數據比較，在低于傅汝德數0.42時，小水線面雙體船均為埋首狀態，在0.37附近出現最大埋首狀態，在傅汝德數0.3附近縱傾曲線出現拐點，傅汝德數超過0.42后出現抬首現象。從該曲線可知，小水線面雙體船航行縱傾隨航速變化出現明顯的波動情況。

圖5為計算升沉值與試驗數據比較，從圖中可知，小水線面雙體船航行時船中吃水都會變大，即出現下沉情況。

從圖4、圖5可以看出，數值模擬的升沉、縱傾值與試驗值的變化趨勢相符。縱傾值在傅汝德數0.3左右有數值偏差，升沉預報與試驗結果吻合較好。

圖6為濕表面面積隨航速變化的曲線，在計算的航速范圍內，濕表面面積隨航速增加而增加。

該船型在Fn=0.3時產生的興波如圖7、圖8所示。

圖7 波形圖（Fn=0.3）

圖8 波形云圖（Fn=0.3）

從圖中可見，所形成的興波很好地滿足了遠方輻射條件，同時波峰、波谷的位置也與實際情況相符。

3 結 論

本文基于DAWSON方法對小水線面雙體船興波阻力、縱傾及升沉進行了預報，并給出船模計算結果與模型試驗結果的對比。根據計算結果與模型試驗結果對比可見，數值模擬的升沉、縱傾值與試驗值的變化趨勢相符。縱傾值在傅汝德數0.3左右有數值偏差，升沉預報與試驗結果符合較好。濕表面面積隨航速變化較大，其隨航速增加而增加。

［1］ 黃鼎良.小水線面雙體船性能原理［M］.北京： 國防工業出版社，1993.

［2］ Dawson C W. A practical computer method for solving ship-wave problems［C］ . the Second ICNSH，Beckly：University Extension Publication， 1977：30-38.

［3］ 劉應中.船舶興波阻力理論［M］.北京：國防工業出版社，2001.

［4］ Dawson CW. Calculations with the xyz free surface programfor five ship models［A］.Proceedings of Workshop on Ship Wave-Resistance Computations［C］//David W.Taylor Naval Ship Research and Development Center，Bethesda，MD，1979.

［5］ Hess J L， SmithAMO.Calculation of non-lifting potential flowabout arbitrary three-dimensional bodies［J］. Journal of Ship Research，1964（2）： 22-44.

［6］ 蔡新功.基于Michell理論的多體船阻力計算［J］.船舶，2007（2）：5-9.

Prediction of wave-making resistance and fl oating state of SWATH based on DAWSON method

FAN Jing-feng1LI Yun-bo2
(1. Representative Offi ce of Naval Warship Design & Research, Shanghai 200011, China; 2. Harbin Engineering University, Harbin 150001, China)

A small waterplane area twin hull （SWATH） is investigated based on Dawson method. It iteratively solves the wave-making resistance， trim， sinkage and the wetted surface area under the diff erent fl oating states of the ship model at diff erent speeds through the quick division of the local ship grid， and the prediction results are compared with the model test results.

DAWSON method； SWATH； wave-making resistance； trim； sinkage

U661.31

A

1001-9855（2015）06-0020-04

2015-01-13；

2015-03-18

范井峰（1978-），男，工程師，研究方向：艦船設計審查。

李云波（1968-），女，博士，教授，研究方向：船舶水動力性能。