非對稱雙體船片體構型及阻力性能的數值研究

謝云平，袁雙雙，徐曉森，胡冬芳

(江蘇科技大學 船舶與海洋工程學院，江蘇 鎮江 212003)

非對稱雙體船片體構型及阻力性能的數值研究

謝云平，袁雙雙，徐曉森，胡冬芳

(江蘇科技大學 船舶與海洋工程學院，江蘇 鎮江 212003)

對于海上風電雙體運維船而言，興波阻力是其總阻力的重要組成部分。本文基于CFD技術對非對稱雙體船進行初步研究，通過改變雙體船左右片體角度以及改變單側片體肥瘦程度來實現內外片體排水體積差的目的2個方面來研究非對稱雙體船的阻力性能，比較各方案的阻力，綜合考慮最終得出最優方案，為今后雙體風電運維船的船型設計與優化提供一個新的研究方向。

海上風電；運維船；非對稱雙體船；CFD

0 引 言

近年來，隨著海上風電場的迅速發展，風電維護船的需求日益增大。雙體船被廣泛用于風電維護船船型[1]，因為其具有阻力性能好、甲板面積大、舒適性好、穩性好、操縱性優良等優點，但是雙體船片體間存在興波干擾，一般來講會增加附加干擾阻力。本文受三體船不對稱型片體構型方案的啟發[2]來研究雙體船不對稱時阻力的性能，從通過改變雙體船左右片體角度以及單個片體肥瘦程度來實現片體內外排水體積差的目的2個方面來研究非對稱雙體船的阻力性能，為今后雙體風電運維船的發展提供一個新的研究方向。

1 運維船主尺度及不對稱方案

1.1 運維船主尺度

本研究工作主要針對某風電運維船展開進行，其模型如圖1所示，其主要要素見表1。

1.2 運維船不對稱方案

1.2.1 左右片體呈一定角度的不對稱方案

根據海上風電場維護船航行的環境和使用特點，本方案以某雙體風電維護船為基本船型，通過改變2個片體之間的角度，并且保持船中處的型寬不變，從而生成左右片體呈一定角度的不對稱方案，如表2所示，并在此基礎上研究各方案在波浪中的阻力性能。

1.2.2 單側片體內外不對稱

本方案以某雙體風電維護船為基本船型，在保證船體總排水量、浮心位置以及水線寬度不變的條件下，運用NAPA軟件分別對船體的型線進行修改以此來達到單個片體內外排水比之差的目的（片體內外排水量之差相對船體總排水量的比值），如圖3所示，利用Gambit進行模型處理，并在此基礎上進一步研究4種方案在波浪中的阻力性能。

表 1 某風電運維船主要要素Tab. 1 Parameters of the maintenance ship

表 2 某風電運維船片體旋轉角度Tab. 2 The rotation angle of the maintenance ship

表 3 內外片體排水量之差Tab. 3 The difference between internal and external displacement

2 阻力數值方法研究

本文選取荷蘭代爾夫特大學Delft372雙體船船模進規則波中的計算驗證，模型如圖4所示，船模主尺度如表4所示。

表 4 Delft船模主尺度Tab. 4 Main scale of Delft372

2.1 計算域大小的設置

相關理論的研究和實踐經驗表明，計算域的大小應根據船模的尺度從而控制在一個合理的范圍內[4]，一方面，計算域的大小在船長和吃水方向仍以船長L的倍數為量度；另一方面，根據雙體船船寬船長比一般單體船大的特點，計算域大小在船寬方向初步選取船寬的倍數為量度。一般計算區域的長度應為8～10倍船長L，寬度則應保證在2.5～3.0倍船長L[5]。為了節省計算的時間成本，根據船型和流場的對稱性，本文初步采用半模計算，模型縮尺比為1:10，

2.2 網格劃分

本文是用全六面體非結構網格生成器HEXPRESS來生成網格，其單元基本按笛卡爾坐標方向排列，在物面附近網格被細化并投影到物面上，形成貼體網格，準確描述物體外形。在計算域中的絕大部分區域，網格單元都接近于長方體，網格質量非常高，能達到很高的計算精度，并能很好地模擬粘性邊界層。本文計算雙體船在規則波中的阻力，考慮流場的流動特性，因此各方案網格的劃分方法相同。

本文運用Fine-Marine軟件對各方案的模型進行網格劃分，以片體呈內八2°為例，計算域和網格如圖5和圖6所示。

根據試驗資料[3]，船模在一系列傅氏數下的實驗數據和在CFD下的仿真結果如表5和圖7所示。

表 5 各傅氏數下仿真值和實驗值的大小Tab. 5 The simulation and experimental values with different Froude number

比較圖7中的2種方案，并且隨著傅氏數的變化仿真值與試驗值具有相同的變化趨勢。從阻力角度來看，仿真值與試驗值相比差別很小，阻力值的最大誤差均在7%以內。因此本文可以沿用該方法進行雙體船的阻力仿真計算研究。

3 左右片體呈角度對阻力性能的研究

根據前面所選定的方案，本文選取當Fr=0.546和Fr=0.729時分別對表2中各工況進行阻力仿真計算，得出的結果如表6所示。

表 6 不同傅氏數情況下各方案阻力的大小Tab. 6 The resistance of each scheme in different Froude number

如圖8所示，當Fr=0.546的時候，表6中方案1中片體向外旋轉（呈外八狀）各個角度，根據圖8所示，船體的阻力呈現增大趨勢，并且隨著角度的增大，船體興波阻力增加的幅度也越大。而方案2中片體向內旋轉（呈內八狀）各個角度后，根據圖表所示，船體在旋轉2°和4°的時候的興波阻力小于對稱型片體的興波阻力尤以旋轉角度為2°的時候，船體阻力有明顯的減小，當片體旋轉角度繼續增加，船體的興波阻力開始大于對稱型船體即0°時的興波阻力，并且由圖8總體趨勢看出，片體呈“內八狀”時的興波阻力比呈“外八狀”時的阻力小，這是由于向內旋轉的時候片體間互相的興波干擾變小的緣故，當Fr=0.729時，由圖8所知，船體的興波阻力始終呈增大趨勢，并且隨著角度的增加，增大的趨勢也變大，這是因為由于航速的增加，片體間的興波干擾進一步變大，興波阻力也隨之增加。

以上分析表明，片體的不對稱構型以及船體的航速對阻力的影響規律極其復雜，隨著航速的不同，影響的趨勢也不同，較為傾向性的結論是當Fr=0.546的時候，片體旋轉2°，船體興波阻力最小并且在高航速階段時，船體的阻力呈增長趨勢。

4 單個片體內外不對稱對阻力性能的研究

選取Fr=0.546，Fr=0.729分別對表3里各工況進行阻力仿真計算。得出的結果如表7所示。

表 7 不同傅氏數情況下各方案阻力的大小Tab. 7 The resistance of each scheme in different Froude number

當Fr=0.546時，表7中方案1中即內瘦外肥的片體由圖9所知，在9%處的阻力大于原對稱型船體阻力，當排水比為3%，5%，7%的時候，船體的阻力小于原對稱型阻力，并當排水比之差為5%（見圖10）時，阻力有明顯的減小，而方案2中曲線所示，船體阻力較原對稱型片體呈明顯的上升勢頭，這是由于內肥外瘦的片體在水中相互干擾的作用變大，導致船體興波阻力增大的緣故。而當Fr=0.729時，根據圖9所示，船體阻力的變化規律與Fr=0.546時幾乎一樣。

以上結果表明內瘦外肥型片體的阻力小于內肥外瘦型片體，并且在片體內瘦外肥型片體的排水比之差為5%的時候，船體阻力明顯小于原對稱型船體。

5 角度及片體肥瘦綜合因素對阻力的影響

由上述2個方案可知，當雙體船在片體向內旋轉2°呈內八時以及單側片體左右排水比之差（內瘦外肥）為5%的時候，為上述2個方案的最佳方案，阻力有明顯減小，本章節結合上述2個方案綜合考慮，選取Fr=0.546，Fr=0.729分別對單側片體的排水比之差為5%（內瘦外肥）并且向內旋轉2°呈內八形狀時的情況下進行阻力仿真計算，得出的結果如圖10所示，當Fr=0.546時阻力為12.96 N，小于上述兩方案的最優解，而當Fr=0.729時阻力為17.97 N，參考表6和表7所知，阻力大小介于上述兩方案最優解之間。

以上分析表明，由于船體特殊的構型，在不同航速下對阻力的影響規律極其復雜，隨著航速的不同，影響的趨勢也不同，這是由于在不同航速下，片體間相互的興波干擾也會隨著航速發生不規律的變化。由于運維船在工作情況下，航速會隨時發生變化，所以即使本方案在Fr=0.546時，阻力為各方案的最優解，但是由于航速的提高，阻力并非是該航速情況下的最優解。較之片體為內瘦外肥時的情況，并非是最佳方案。

6 結 語

本文基于CFD技術，采用FINE/Marine軟件探究了非對稱雙體船不同片體構型情況下阻力研究，得出以下結論：

1）在片體角度對阻力的影響方案中，在Fr=0.546時，片體角度為向內旋轉2°時的阻力性能最優，但在各個航速下，該方案并非是最優解。

2）單側片體內外不對稱對阻力的影響的方案中，方案1中（片體為內瘦外肥）當片體內外排水比之差為5%的時候，船體的阻力較之對稱型船體有明顯的減小，可視為不對稱雙體船片體構型的最佳方案。

3）綜合上述2個方案考慮研究的情況下，雖然在Fr=0.546及Fr=0.729時，船體的阻力均小于對稱型船體時的阻力，且當Fr=0.546時為各個方案的最優解，不過隨著航速增加，在Fr=0.729時阻力大于內瘦外肥型片體，并非是最優，結合運維船的工作情況及實際建造情況，本方案并非為最佳方案。

本文提出的非對稱雙體船阻力優化對解決工程實際問題有一定的借鑒作用；對雙體船片體構型的研究為風電維護船的船型優化提供了一定的參考價值。同時，本文僅對非對稱雙體船的阻力進行了初步研究，對于片體在波浪中的耐波性以及水動力響應方面還有待進一步研究。

張秀萍. 海上風電場維護船船型及相關性能研究[D]. 鎮江: 江

[1]蘇科技大學, 2014.

[2]高雷, 常亮. 三體船不對稱型側片體構型及阻力性能研究[C]//第十九屆中國國際高性能船學術報告會暨中國國際游艇設計建造技術論壇, 2014. GAO Lei, CHANG Liang. Trimaran type asymmetric side body configuration and resistance performance study[C]//The 19th China International High Performance Ship Academic Report and China International Yacht Design and Construction Technology Forum, 2014.

[3]TERESA C, FREDERICK S, SERGIO B, et al. Numerical investigation of the seakeeping behavior of a catamaran advancing in regular head waves[J]. Ocean Engineering, 2011, (38): 1806–1822.

[4]STERN F, WILSON R V, COLEMAN H W, et al. Comprehensive approach to verification and validation of CFD simulations. 1: Methodology and procedures[J]. Journal of Fluids Engineering, 2001, 123(4): 793–802.

[5]謝云平, 徐殷茵. 海上長球首球尾雙體風電維護船波浪增阻研究[J]. 艦船科學技術. 2013, 38(2): 74–79 XIE Yun-ping, XU Yin-yin. The research of wave added resistance of maintenance catamaran in wind farm with big bulbous head and stern in regular waves[J]. Ship Science and Technology. 2013, 38(2): 74–79.

Asymmetric catamaran hull configuration and the relevant research of resistance performance

XIE Yun-ping, YUAN Shuang-shuang, XU Xiao-sen, HU Dong-fang
(School of Naval Architecture and Ocean Engineering, Jiangsu University of Science and Technology, Zhenjiang 212003, China)

For offshore wind power catamaran maintenance ship, wave resistance is an important part of the total resistance. In this paper, based on the CFD technology to a preliminary study of asymmetric catamaran, by changing the angle about catamaran hull and change one piece body lean degree to achieve the purpose of the drainage volume difference inside and outside of block two aspects to study the resistance, comparing the resistance of the programs, finally it is concluded that the optimal solution,which provide a new research direction for the future offshore wind power catamaran maintenance ship design and optimization.

offshore wind；maintenance of the ship；asymmetric catamaran；CFD

U661.1

A

1672 – 7649(2017)08 – 0028 – 04

10.3404/j.issn.1672 – 7649.2017.08.006

2016 – 09 – 07

謝云平(1964 – )，男，研究員，研究方向為船舶設計。