張連洪，張宇航，楊紹瓊，馬?偉，牛文棟，侯鵬榮

一種基于低速、小振幅周期運動CFD數(shù)值模擬的水下滑翔機附加質量求解方法

張連洪1, 2, 3，張宇航1，楊紹瓊1, 2, 3，馬?偉1, 2, 3，牛文棟1, 2, 3，侯鵬榮4

(1. 天津大學機構理論與裝備設計教育部重點實驗室，天津 300350；2. 天津大學青島海洋技術研究院，青島 266237；3. 青島海洋科學與技術試點國家實驗室海洋觀測與探測聯(lián)合實驗室，青島 266237；4. 大連理工大學機械工程學院，大連 116024)

附加質量；水下滑翔機；計算流體動力學；水池試驗

水下滑翔機作為一種新型水下航行器，可通過浮力驅動實現(xiàn)多剖面連續(xù)滑翔運動，由于其強大的續(xù)航能力和較高的隱蔽性，被廣泛應用于海洋調查和軍事探測等領域[1]．隨著海洋科學的不斷發(fā)展，人們對水下滑翔機的滑翔性能以及操縱性能要求越來越高．而附加質量是水下滑翔機水動力重要影響因素之一，因此，人們對附加質量計算精度提出了更高的要求．

Lee等[2]分別采用垂直平面運動機構(vertical plan motion mechanism，VPMM)測試與計算流體動力學(computational fluid dynamics，CFD)方法對水下無人航行器的附加質量進行求解，驗證了VPMM設備在水下航行器水池試驗中的適用性；周景軍等[3]基于相對運動法，利用CFD分別模擬水下航行體進行非定常運動與定常運動，得到流體動力的差值來求解附加質量；楊磊等[4]采用基于勢流理論的H-S面元法和基于PMM試驗的CFD技術求得水下滑翔機附加質量，并將結果進行對比，相對小于11%；莫慧黠?等[5]基于動參考系思想，避免了相對運動法添加動量源項以及CFD模擬時動網格等問題，附加質量求解誤差在10%以內；宋麗華[6]基于有限體積法，求解基本速度勢的拉普拉斯方程計算附加值量；劉智麗等[7]將物體在靜止來流中的運動問題轉換為物體靜止而來流時變的計算問題，避免了在CFD中采用動網格技術，同時計算了在黏性流體中橢球體的附加質量；此外，還可以采用細長體假設或切片理論求解附加質量，但只適用于簡單三維物體，如細長圓柱等[8-9]；對于復雜外形水下航行器可使用拖曳水池試驗的方法求解附加質量，但該方法存在費用高、實驗周期長等問題．目前，最常用的求解形狀復雜的水下航行器附加質量的方法是基于勢流理論的H-S面元法和數(shù)值模擬平面運動機構PMM試驗的CFD方法．其中，PMM試驗法主要包括物體的純升沉運動、純橫蕩運動、純俯仰運動、純搖艏運動，其難點在于模擬過程涉及到流體流動分離、航行器本體和附體之間的干擾作用等．

水下滑翔機依靠機翼提供的升力實現(xiàn)鋸齒形剖面運動，依靠尾舵實現(xiàn)其滑翔的滑翔穩(wěn)定性[10-11]．在進行PMM試驗時，機翼會產生較大的流動分離現(xiàn)象．因此本文提出了一種采用CFD技術模擬水下滑翔機低速、小振幅周期運動的試驗方法，以減小機翼與流體的相互作用．同時，本文采用基于勢流理論的H-S面元法編寫了附加質量計算程序，用以驗證CFD方法的準確性．此外，通過對水下滑翔機等比例模型進行水池試驗，將得到的結果與CFD數(shù)值模擬的結果進行對比，以驗證這種方法的可行性．

1?附加質量求解原理

1.1?研究對象

本文研究對象為天津大學自主研發(fā)的海燕-L長航程水下滑翔機，如圖1所示．海燕-L主要由耐壓艙、機翼、前后導流罩和天線構成．耐壓艙內部包括浮力調節(jié)系統(tǒng)、姿態(tài)調節(jié)系統(tǒng)、控制系統(tǒng)、通信和導航系統(tǒng)等[13]．表1為海燕-L的主要特征參數(shù)．

圖1?海燕-L水下滑翔機

表1?海燕-L主要特征參數(shù)

Tab.1?Main characteristic parameters of Petrel-L

圖2?海燕-L體坐標系

1.2?求解原理

采用空間固定坐標系，平面為豎直面，軸根據(jù)右手定則垂直于豎直面方向．設定海燕-L在流體域中沿、、軸同時進行低速度、小振幅的周期性升沉和橫蕩運動．圖3為水下滑翔機運動軌跡放大后的示意圖，實際模擬過程中振幅很?。?/p>

圖3?海燕-L運動軌跡示意

海燕-L運動方程為

則加速度為

將式(1)、式(2)代入式(3)，可得

在相同流體域中，分別使水下滑翔機繞體坐標系的、、軸做純搖艏運動．設定水下滑翔機運動方程為

將式(6)代入式(7)，可得

2?利用CFD計算水下滑翔機附加質量

2.1?CFD模型建立

本文采用CFD商用軟件STAR-CCM＋對水下滑翔機低速、小振幅周期運動進行模擬．STAR-CCM+具有強大的網格支持能力，支持動網格、滑移網格以及重疊網格等．其中，重疊網格綜合了滑移網格與動網格的優(yōu)點，既能夠很好地描述部件運動，又能夠保證物體運動過程中的網格質量，可以對任意不規(guī)則運動的物體進行求解．因此本文采用重疊網格技術對水下滑翔機運動進行模擬．

本算例需要建立兩個流體域，如圖4所示，分別對這兩個域劃分網格，即組件網格和背景網格．其中，組件網格為任意長方體與其內部水下滑翔機外表面之間區(qū)域的網格(該網格可由長方體域與內部水下滑翔機進行布爾減運算得到，長方體尺寸略大于水下滑翔機整體尺寸)，背景網格為流體域網格，長度為3(2.6m)，寬、高分別為1.5．本文將流體域劃分為切割體網格單元．

圖4?邊界條件

根據(jù)第1.2節(jié)所述方法，模擬水下滑翔機進行低速小振幅運動，監(jiān)控水下滑翔機所受水動力和力矩，即可得到水下滑翔機的8個附加質量．

2.2?網格數(shù)量與計算步長設置

圖5?網格數(shù)量設置

可以看出，不同網格數(shù)量對計算結果的影響不大，附加質量上下波動可能與計算步長的選取有關．為使附加質量結果達到一定精度，并節(jié)省計算資源，本文進行仿真時采用85×104的網格數(shù)量．

圖6?計算步長設置

可以看出當計算步長減小為0.02s，即一個周期計算200次時，計算結果已經不再隨著計算步長的減小而發(fā)生明顯的變化．為節(jié)省計算資源，提高附加質量計算效率，本文采用計算步長0.02s計算水下滑翔機附加質量，此時能夠達到很高的精度．

2.3?速度、振幅影響分析

圖7?速度幅值與附加質量的關系

圖8?角速度幅值與附加質量的關系

2.4?計算結果

本文最終使用網格數(shù)為85×104的模型，計算步長為0.02s，對水下滑翔機低速、小振幅周期運動進行數(shù)值模擬，速度幅值設置為0.01m/s，角速度幅值設置為0.005rad/s，得到水下滑翔機8個附加質量值如表2所示．

表2?基于CFD的附加質量計算值

Tab.2?Value of additionalmassbased on CFD

3?結果驗證及分析

3.1?基于H-S面元法的水下滑翔機附加質量計算

H-S面元法是基于勢流理論，采用分布奇點法來求解物體附加質量的，所針對的是物體在沒有黏性、不可壓縮的無旋流體中的運動問題[13]．

假設物體在靜止的理想流體中以速度作非定常直線運動，則速度勢滿足控制方程為

邊界條件為

由勢流理論可知，流場內任意點的速度勢為

將式(13)兩端取法相導數(shù)，帶入式(12)可以得到

則物體的附加質量為

本文基于上述方法使用Matlab編寫了水下滑翔機附加質量計算程序．與CFD數(shù)值模擬法不同的是，CFD方法是對流體域進行體網格劃分，而H-S面元法是對水下滑翔機進行表面網格劃分．考慮到劃分網格時邊長相同的四邊形面元較三角形面元會節(jié)省大量計算資源，而三角形面元相比四邊形面元更能貼合物體表面，因此本文采用四邊形占優(yōu)來劃分網格．本文編寫的程序能兼容三角形和四邊形兩種形狀面元，從而達到既貼合物體表面又能節(jié)省計算資源的目的．

水下滑翔機三維模型的坐標原點位置應與浮心重合．使用workbench的mesh模塊對水下滑翔機表面進行網格劃分，如圖9所示，最小網格尺寸為15mm，表面網格數(shù)量為17293．

圖9?Petrel-L表面網格

將軟件輸出的網格文件導入程序，即可得到海燕-L水下滑翔機的附加質量．

3.2?水池試驗

本文針對海燕-L水下滑翔機進行了水池試驗，以驗證上述CFD數(shù)值模擬方法的可行性．

試驗地點為哈爾濱工程大學，試驗設備采用水平型循環(huán)水槽，它與普通水池有所不同，被測試的模型沒有前進速度，而周圍流場是持續(xù)運動的，因此試驗時間不受場地限制，可實現(xiàn)長時間的自動化測量，極大縮減了試驗周期和成本．

長航程水下滑翔機的水池試驗主要通過力天平測力傳感器，在循環(huán)水槽的均勻流場中，使用平面運動機構(PMM)驅動水下滑翔機做指定方式的運動. 試驗過程包括水平面定長變漂角六分力測量、水平面純搖艏簡諧運動六分力測量、水平面純橫蕩簡諧運動六分力測量、垂直面定常變攻角六分力測量、重直面純俯仰簡諧運動六分力測量以及垂直面純升沉簡諧運動六分力測量．同時，采集水下滑翔機運動過程中所受的力和力矩，繪制并擬合受力曲線，求取水下滑翔機附加質量．圖10為水下滑翔機模型的加工與裝配圖．圖11為試驗過程．

圖10?海燕-L模型加工與裝配

圖11?試驗過程

3.3?結果與分析

表3?試驗結果

Tab.3?Experimental results

(1) 如圖11(a)所示，天平支撐桿與滑翔機殼體連接處存在較大縫隙，在進行水池試驗時會產生較大的內流場；

(2) 受到航道尺度的限制，水下滑翔機純升沉簡諧運動受到邊界的干擾較為嚴重；

(3) 水下滑翔機機翼材質為碳纖維，在進行純升沉簡諧運動時機翼會發(fā)生較大的抖動．

4?結?語

本文以海燕-L水下滑翔機為研究對象，提出了一種采用CFD數(shù)值模擬法求解水下滑翔機附加質量的新方法，該方法只需要建模一次，在同一個流體域中模擬水下滑翔機的4種運動方式，監(jiān)控其所受的水動力和力矩，即可得到海燕-L的8個附加質量，相比文獻[4]中所述的基于PMM試驗的CFD方法，極大提高了附加質量求解效率．同時，本文采用基于勢流理論的H-S面元法，對海燕-L水下滑翔機的附加質量進行求解，得到兩種方法計算結果較為一致，相對誤差基本可以控制在5%以內．最后，本文通過水池試驗對等比例水下滑翔機模型進行試驗驗證，相對誤差基本可以控制在6%以內，可以看出本文提出的CFD方法具有較高的可行性和準確性．由于水池試驗對試驗設備的要求較高，且費用高、周期長，因此，該CFD數(shù)值模擬方法可以在一定程度上代替水池試驗對水下滑翔機附加質量進行求解，且該方法能適用于任意尺寸、形狀的水下滑翔機，對不同型號水下滑翔機水動力外形設計、動力學建模都具有較大的參考價值．

［1］ 何柏巖，杜金輝，楊紹瓊，等. 基于VMD-LSSVM的水下滑翔機深平均流預測[J]. 天津大學學報(自然科學與工程技術版)，2021，54(4)：388-396.

He Baiyan，Du Jinhui，Yang Shaoqiong，et al. Prediction of underwater glider depth-averaged current velocity based on VMD-LSSVM[J]. Journal of Tianjin University(Science and Technology)，2021，54(4)：388-396(in Chinese).

［2］ Lee S K，Joung T H，Cheon S J，et al. Evaluation of the added mass for a spheroid-type unmanned underwater vehicle by vertical planar motion mechanism test[J]. International Journal of Naval Architecture and Ocean Engineering，2011，3(3)：174-180.

［3］ 周景軍，李育英，項慶睿. 一種水下航行體附加質量數(shù)值計算方法[J]. 魚雷技術，2013，21(4)：246-249.

Zhou Jingjun，Li Yuying，Xiang Qingrui. A numerical computation method of additional mass for underwater vehicle[J]. Torpedo Technology，2013，21(4)：246-249(in Chinese).

［4］ 楊?磊，曹軍軍，姚寶恒，等. 水下滑翔機附加質量數(shù)值計算[J]. 艦船科學技術，2016，38(23)：116-120.

Yang Lei，Cao Junjun，Yao Baoheng，et al. Numerical calculation of additional mass of underwater glider[J]. Ship Science and Technology，2016，38(23)：116-120(in Chinese).

［5］ 莫慧黠，黨建軍，羅?凱，等. 水下航行器附加質量數(shù)值計算方法[J]. 水下無人系統(tǒng)學報，2017，25(4)：250-255.

Mo Huixia，Dang Jianjun，Luo Kai，et al. Numerical calculation method of additional mass of underwater vehicle[J]. Journal of Underwater Unmanned Systems，2017，25(4)：250-255(in Chinese).

［6］ 宋麗華. 水下航行體附加質量的計算方法[J]. 水雷戰(zhàn)與艦船防護，2017，25(4)：58-61，54.

Song Lihua. Calculation method of additional mass of underwater vehicle[J]. Mine Warfare and Ship Protection，2017，25(4)：58-61，54(in Chinese).

［7］ 劉智麗，林獻武. 有黏流場中旋成橢球體的軸向附加質量系數(shù)計算[J]. 廈門大學學報(自然科學版)，2019，58(3)：408-414.

Liu Zhili，Lin Xianwu. Calculation of axial additional mass coefficient of rotating ellipsoid in viscous flow field[J]. Journal of Xiamen University(Natural Science Edition)，2019，58(3)：408-414(in Chinese).

［8］ Wang K，Zhang X，Zhang Z Q，et al. Numercal analysis of added mass and damping of floating production，storage and offloading system[J]. Acta Mechanica Sinica，2012，28(3)：870-876.

［9］ Jensen R L，Watts P B，Lawrence J E，et al. Vertical hand force and forearm emg during a high-step rock-on climbing move with and without added mass[C]// Proceedings of 23rd International Symposium on Biomechanics in Sports. Beijing，China，2005：466-469.

［10］沈新蕊，王延輝，楊紹瓊，等. 水下滑翔機技術發(fā)展現(xiàn)狀與展望[J]. 水下無人系統(tǒng)學報，2018，26(2)：89-106.

Shen Xinrui，Wang Yanhui，Yang Shaoqiong，et al. Development of underwater gliders：An overview and prospect [J]. Journal of Unmanned Undersea Systems，2018，26(2)：89-106(in Chinese).

［11］張?亮，李云波. 流體力學[M]. 哈爾濱：哈爾濱工程大學出版社，2001.

Zhang Liang，Li Yunbo. Fluid Mechanics[M]. Harbin：Harbin Engineering University Press，2001(in Chinese).

［12］ Yang Ming，Wang Yanhui，Wang Shuxin，et al. Motion parameter optimization for gliding strategy analysis of underwater gliders[J]. Ocean Engineering，2019，191：106502-1-106502-12.

［13］戴遺山，段文洋. 船舶在波浪中運動的勢流理論[M]. 北京：國防工業(yè)出版社，2008.

Dai Yishan，Duan Wenyang. The Potential Flow Theory of Ship Moving in Waves[M]. Beijing：National Defense Industry Press，2008(in Chinese).

Method for Solving Added Mass of Underwater Glider Based on CFD of Low-Speed and Small-Amplitude Periodic Motion

Zhang Lianhong1, 2, 3，Zhang Yuhang1，Yang Shaoqiong1, 2, 3，Ma Wei1, 2, 3，Niu Wendong1, 2, 3，Hou Pengrong4

(1. Key Laboratory of Mechanism Theory and Equipment Design of Ministry of Education，Tianjin University，Tianjin 300350，China；2. Qingdao Institute for Ocean Engineering of Tianjin University，Qingdao 266237，China；3. The Joint Laboratory of Ocean Observing and Detection，Pilot National Laboratory for Marine Science and Technology(Qingdao)，Qingdao 266237，China；4. School of Mechanical Engineering，Dalian University of Technology，Dalian 116024，China)

added mass；underwater glider；computational fluid dynamics(CFD)；pool test

U674.941

A

0493-2137(2022)03-0239-08

10.11784/tdxbz202010039

2020-10-21；

2020-11-07.

張連洪（1963—??），男，博士，教授，zhanglh@tju.edu.cn.

馬?偉，wei.ma@tju.edu.cn.

國家重點研發(fā)計劃資助項目(2019YFC0311803，2019YFC0311701)；國家自然科學基金資助項目(11902219)；天津市自然科學基金資助項目(18JCJQJC46400).

Supported by the National Key Research and Development Program of China(No.2019YFC0311803，No.2019YFC0311701)，the National Natural Science Foundation of China(No.11902219)，the Natural Science Foundation of Tianjin，China (No.18JCJQJC46400).

(責任編輯：王曉燕)