“海鱘4000”水下滑翔機水動力特性與滑翔性能研究

劉來連, 閔強利, 張光明

“海鱘4000”水下滑翔機水動力特性與滑翔性能研究

劉來連, 閔強利, 張光明

(中國船舶重工集團公司 第710研究所, 湖北 宜昌, 443003)

與一般水下滑翔機相比, 4 000 m水下滑翔機使用深度更深, 體型相對更大。為了克服海水密度變化的影響, 其浮力調(diào)節(jié)系統(tǒng)的體積調(diào)節(jié)量也更大。準(zhǔn)確的動力學(xué)模型以及精確的水動力系數(shù)是實現(xiàn)其控制系統(tǒng)設(shè)計以及精確導(dǎo)航的基礎(chǔ)。文中將在結(jié)構(gòu)參數(shù)已定, 動力學(xué)模型已知的前提下獲得精確的水動力系數(shù), 對水下滑翔機性能進行預(yù)測, 并為控制系統(tǒng)提供控制變量和優(yōu)化方向作為研究的重點。以“海鱘4 000”水下滑翔機為研究對象, 采用理論計算、計算流體動力學(xué)(CFD)仿真與試驗數(shù)據(jù)分析相結(jié)合的方法, 分析了水平翼參數(shù)對其滑翔性能的影響; 獲得了其在縱平面內(nèi)做定常運動的水動力系數(shù), 并對其性能進行預(yù)測; 通過試驗數(shù)據(jù)分析, 獲得典型滑翔狀態(tài)下不同凈浮力對應(yīng)的滑翔速度, 并與理論計算進行對比分析, 驗證理論計算的正確性。研究結(jié)果表明: 試驗實際俯仰角控制與理論計算誤差在5%以內(nèi); 試驗實際凈浮力下的滑翔速度與理論計算誤差在15%以內(nèi), 滿足控制系統(tǒng)相應(yīng)控制量20%的誤差需求。文中研究可為水下滑翔機控制系統(tǒng)的進一步優(yōu)化和完善提供參考。

水下滑翔機; 滑翔性能; 水動力性能; 滑翔速度

0 引言

海洋是一座天然寶藏, 它既是人類實現(xiàn)可持續(xù)發(fā)展的重要物質(zhì)來源, 又具有舉足輕重的軍事戰(zhàn)略意義。在海洋開發(fā)的熱潮及現(xiàn)代化戰(zhàn)爭呈現(xiàn)“信息化局部戰(zhàn)爭”態(tài)勢下, 無人水下航行器(un- manned undersea vehicle, UUV)越來越被各國海軍所重視。水下滑翔機就是近年來誕生的一種令人耳目一新的UUV, 其獨特的設(shè)計顯示出了極大的優(yōu)勢。它不需要任何燃料作為動力, 通過浮力調(diào)節(jié)裝置實現(xiàn)自身浮力改變, 在凈浮力、阻力及升力的作用下, 實現(xiàn)鋸齒形滑翔觀測運動。自1989年Stommel[1]提出水下滑翔機的概念之后, 水下滑翔機發(fā)展迅速, 較著名的有美國Webb研究所的slocum水下滑翔機[2]、華盛頓大學(xué)的Sea- glider水下滑翔機[3]和Scripps海洋研究所的spray水下滑翔機[4]?？傮w來說, 這些水下滑翔機的機體尺寸和質(zhì)量相差不多(體長約2 m, 重約50 kg), 機身形狀均為細(xì)長回轉(zhuǎn)體, 這種設(shè)計尺寸小、阻力低, 同時保證了一定的裝載量和續(xù)航力。

上述機構(gòu)還對水下滑翔機運動及控制特性進行了研究, 并對其運動路徑進行數(shù)值仿真, 驗證了其運動及控制規(guī)律。然而這種僅僅依靠電池驅(qū)動的水下滑翔機需要獲得長航程和高續(xù)航力, 除了需要配備高能電池外, 最重要的是要提高滑翔效率, 而升力和阻力是決定效率的關(guān)鍵因素。我國在水下滑翔機理論和技術(shù)方面, 天津大學(xué)、中國科學(xué)院沈陽自動化研究所、中船重工第710研究所、華中科技大學(xué)、浙江大學(xué)、中船重工第702研究所、中國海洋大學(xué)以及國家海洋技術(shù)中心等單位都相繼開展了研究工作。其中: 天津大學(xué)研制的“海燕”混合驅(qū)動水下滑翔機綜合了典型的水下滑翔機和傳統(tǒng)自主水下航行器(autonomous un- derwater vehicle, AUV)的特點[5]; 中國科學(xué)院沈陽自動化研究所開展了“海翼”水下滑翔機研制, 用于深水海洋環(huán)境監(jiān)測[6]; 浙江大學(xué)在水下滑翔機的運動穩(wěn)定性和控制方面開展研究, 并開發(fā)了試驗樣機[7]。

文中以第710研究所“問海計劃”項目中“海鱘4 000 ”水下滑翔機為研究對象, 采用理論計算、計算流體動力學(xué)(computational fluid dynamics, CFD)仿真與試驗數(shù)據(jù)分析相結(jié)合的方法, 分析了水平翼參數(shù)對水下滑翔機滑翔性能的影響; 獲得了其在縱平面內(nèi)做定常運動的水動力系數(shù), 并對水下滑翔機的性能進行預(yù)測; 通過試驗數(shù)據(jù)分析, 獲得典型滑翔狀態(tài)下不同凈浮力對應(yīng)的滑翔速度, 并與理論計算進行對比分析, 驗證了理論計算的正確性。

1 結(jié)構(gòu)設(shè)計與水動力參數(shù)確定

1.1 “海鱘4000”水下滑翔機結(jié)構(gòu)

“海鱘4000”水下滑翔機主要由滑翔機殼體、浮力調(diào)節(jié)裝置、姿態(tài)調(diào)節(jié)裝置及通信天線等組成,如圖1所示?？v剖面的主要運動方式為滑翔運動, 與傳統(tǒng)滑翔機的滑翔運動類似, 通過浮力調(diào)節(jié)裝置和質(zhì)心調(diào)節(jié)裝置實現(xiàn)縱剖面鋸齒形滑翔運動。

圖1 “海鱘4000”水下滑翔機結(jié)構(gòu)示意圖

1.2 動力學(xué)模型

水下滑翔機的運動形式主要是在縱平面的“之”字形運動。因此, 首先需要建立其在縱平面內(nèi)的動力學(xué)模型。

1) 坐標(biāo)系建立

圖2 “海鱘4000”水下滑翔機坐標(biāo)系

2) 水動力模型

水下滑翔機在水下運動時, 會受到水動力的作用, 而水動力的大小直接決定著滑翔機的運動結(jié)果。水下滑翔機在縱平面運動時, 所受到的水動力主要有升力, 阻力和縱傾力矩。文獻[9]研究表明, 小攻角情況下,和與攻角之間近似為線性關(guān)系、與攻角之間近似為二次函數(shù)關(guān)系, 并可表示為

1.3 水平翼參數(shù)特性分析

水下滑翔機的阻力主要由主體的濕表面積和外形決定, 升力由水平翼提供, 水平翼升阻比的計算公式為

由圖3可知, 展弦比一定的情況下, 水平翼升阻比隨攻角先急劇增加后逐漸減小, 且存在1個最大值; 隨著后掠角增大, 最大升阻比保持不變, 其對應(yīng)的攻角有所增大; 最大升阻比對應(yīng)的攻角一般在6°~10°之間。

圖3 水平翼升阻比隨后掠角、攻角變化關(guān)系曲線

由圖4可知, 后掠角一定, 水平翼升阻比隨攻角先增后減, 同樣存在1個最大值; 隨著展弦比增大, 最大升阻比增大, 其對應(yīng)的攻角有所增大; 最大升阻比對應(yīng)的攻角一般在6°~10°之間。

圖4 水平翼升阻比隨展弦比、攻角變化關(guān)系曲線

由于滑翔機正?；锠顟B(tài)下所對應(yīng)的攻角為小攻角, 當(dāng)攻角一定時(取小攻角), 升阻比隨后掠角、展弦比的關(guān)系如圖5所示。升阻比隨后掠角的增加逐漸減小; 隨著展弦比的增加, 升阻比增加, 但隨著后掠角的增加, 減小的趨勢更快。基本可以認(rèn)為, 在小攻角情況下, 后掠角越小, 相應(yīng)升阻比越大。

圖5 水平翼升阻比隨展弦比、后掠角變化關(guān)系曲線

根據(jù)上述計算分析, 可為“海鱘4000”水下滑翔機確定合適的展弦比和后掠角。

1.4 水動力參數(shù)確定

常用的水動力系數(shù)獲取方法包括理論計算、試驗和CFD方法。理論計算方法主要采用經(jīng)驗公式進行計算, 僅對于細(xì)長體外形及運動過程較為簡單的模型有效[11]。試驗方法主要包括風(fēng)洞試驗和拖曳試驗等, 是目前比較常用且有效的手段, 但是試驗成本高, 且耗時較長。CFD方法以計算機為硬件基礎(chǔ), 通過各類數(shù)學(xué)方法, 能夠解決各種流體力學(xué)難題。

文中針對“海鱘4000”水下滑翔機, 采用CFD方法, 計算出了其水動力參數(shù)。計算狀態(tài)見表1。

表1 水動力參數(shù)計算狀態(tài)表

計算域為長度10、直徑3.5的包圍滑翔機的圓柱體(為滑翔機長度), 如圖6所示。其旋轉(zhuǎn)軸與滑翔機模型的對稱軸重合, 進流邊界面(入口)為圓柱體的前端面, 距模型首端的距離為4, 出流邊界面(出口)為圓柱體的后端面, 距模型尾端的距離為5。整個計算域采用非結(jié)構(gòu)化網(wǎng)格, 網(wǎng)格總數(shù)為120萬, 網(wǎng)格質(zhì)量0.3以上, 以保證求解精度。圖7給出了計算結(jié)果示意圖。

圖6 計算域示意圖

圖7 水下滑翔機表面壓力云圖

圖8 不同速度和攻角下的阻力曲線

2 滑翔性能分析

根據(jù)以上水動力模型, 對“海鱘4000”水下滑翔機的滑翔運動進行受力分析, 得到以下平衡方程

圖10 不同速度和攻角下的俯仰力矩曲線

聯(lián)合式(1)、式(2)和式(13)可得

為使式(14)有解, 需滿足

根據(jù)式(2)、式(10)和式(17), 可得速度與滑翔角的關(guān)系

圖11給出了“海鱘4000”水下滑翔機在有效滑翔角范圍內(nèi)攻角與滑翔角的關(guān)系曲線。由圖可知, 該滑翔機攻角范圍為(0°~ 4.9°)＆(–4.4°~ 0°); 俯仰角范圍為(12.6° ~ 90°)＆(–11.1° ~ –90°)。

圖11 攻角與滑翔角關(guān)系曲線

圖12給出了“海鱘4000”水下滑翔機在不同凈浮力下水平速度與滑翔角的關(guān)系。從圖中可以看出, 水平速度隨著滑翔角的增加先增后減。在確定的凈浮力下, 在有效滑翔角內(nèi)存在1個最大水平速度, 該速度所對應(yīng)的滑翔角約為36o, 對應(yīng)的攻角為–1.45o, 俯仰角為34.55o。

圖12 水平速度與滑翔角、凈浮力的關(guān)系曲線

圖13~圖16給出了“海鱘4000”水下滑翔機在典型滑翔角下, 滑翔速度、水平速度、垂直速度與凈浮力的關(guān)系曲線。

圖13 滑翔角為18o時速度與凈浮力關(guān)系曲線

圖14 滑翔角為27o時速度與凈浮力關(guān)系曲線

圖15 滑翔角為36°時速度與凈浮力關(guān)系曲線

由圖可知, 凈浮力增加可以同時增加滑翔的水平速度和垂直速度, 滑翔角超過36°以后, 水平速度將低于垂直速度。該計算結(jié)果即可與試驗數(shù)據(jù)進行對比, 驗證上述計算過程的正確性。

3 試驗驗證

2019年初對“海鱘4000”水下滑翔機進行了湖上試驗(見圖17), 獲得了不同姿態(tài)下穩(wěn)定滑翔的俯仰角度和滑翔速度數(shù)據(jù)。俯仰角度通過集成在滑翔機機體上的磁羅盤測量獲得,滑翔速度通過機體上安裝的壓力傳感器測量深度間接獲得, 文中選取了2組典型滑翔角下有效試驗數(shù)據(jù)曲線, 分別如圖18~圖20、圖21~圖23所示。需要說明的是, “海鱘4000”水下滑翔機控制程序中滑翔角是設(shè)定變量, 通過如圖11所示的攻角與滑翔角的關(guān)系, 計算出對應(yīng)的俯仰角, 從而對俯仰角進行控制。典型有效試驗數(shù)據(jù)曲線如圖18~圖23所示。

圖16 滑翔角為52°時速度與凈浮力關(guān)系曲線

圖17 “海鱘4000”水下滑翔機湖上試驗

圖18 滑翔角為18°時深度變化曲線

圖19 滑翔角為18°時體積(凈浮力)變化曲線

圖20 滑翔角為18°時俯仰角變化曲線

表2 滑翔角18°時滑翔數(shù)據(jù)表

從試驗實測數(shù)據(jù)與分析計算的數(shù)據(jù)對比來看, 試驗數(shù)據(jù)與分析計算數(shù)據(jù)基本吻合, 誤差均在15%以內(nèi)。存在誤差的原因有: 1) 仿真計算模型過于理想, 與實際樣機存在差異, 造成計算的水動力系數(shù)存在誤差; 2) 實際樣機浮力調(diào)節(jié)系統(tǒng)浮力調(diào)節(jié)量、俯仰角控制精度存在誤差, 控制精度誤差均為20%。因此, 從結(jié)果上看, 理論計算分析完全可以滿足工程設(shè)計的需要, 并對工程設(shè)計具有很好的指導(dǎo)意義。

圖21 滑翔角為27°時深度變化曲線

圖23 滑翔角27°時俯仰角變化曲線

表3 滑翔角為27°時滑翔數(shù)據(jù)表

4 結(jié)束語

文中以“海鱘4 000 ”水下滑翔機為研究對象, 詳細(xì)分析研究了其水動力特性和滑翔性能, 可得如下結(jié)論。

1) 參考相關(guān)文獻對水平翼參數(shù)特性進行了分析研究, 獲得了展弦比、后掠角、攻角與水平翼升阻比的關(guān)系, 最大升阻比對應(yīng)的攻角一般為6°~10°, 可為工程設(shè)計提供參考。

2) 采用CFD方法計算出“海鱘4 000 ”水下滑翔機各個水動力參數(shù), 并以此為基礎(chǔ), 詳細(xì)分析研究了其滑翔性能, 獲得其攻角范圍為(0°~4.9°)＆(–4.4°~0°), 俯仰角范圍為(12.6°~90°)＆(–11.1° ~ –90°), 其最佳抗流能力對應(yīng)的滑翔角為36°。最后詳細(xì)計算了其在不同典型滑翔角下, 滑翔速度、水平速度、垂直速度與凈浮力的關(guān)系曲線, 對其滑翔控制具有很好的指導(dǎo)意義。

3) 對試驗實測數(shù)據(jù)進行了分析研究, 獲得了“海鱘4 000 ”水下滑翔機典型狀態(tài)下的滑翔速度和俯仰角度數(shù)據(jù), 并與分析計算數(shù)據(jù)進行了對比, 試驗數(shù)據(jù)與分析計算數(shù)據(jù)基本吻合, 誤差均在15%以內(nèi), 驗證了分析計算過程的正確性, 滿足控制系統(tǒng)相應(yīng)控制量20%的誤差需求, 并對控制系統(tǒng)的進一步優(yōu)化和完善提供了依據(jù)。

文中的工作可對水下滑翔機的設(shè)計提供參考, 但其理論計算并未考慮海水密度變化及滑翔機受海水壓力產(chǎn)生體積壓縮量的影響, 也并未建立滑翔機改變航向的理論模型, 下一步將對此進行深入研究。

[1] Stommel H. The Slocum Mission[J]. Oceanography, 1989, 2(1): 22-25.

[2] Schofleld O, Kohut J, Aragon D, et a1．Slocum Gliders: Robust and Ready[J]. Journal of Field Robotics, 2010, 24(6): 473-485．

[3] Eriksen C C, Osse T J, Light R D, et al. Seaglider: A Long—Range Autonomous Underwater Vehicle for Oceano Graphic Research[J]．IEEE JournaI of Oceanic Engineering, 2001, 26(4): 424-436．

[4] Rudnick D L, Davis R E, Sherman J T. Spray Underwater Glider Operations[J]．Journal of Atmospheric and Oceanic Technology, 2016, 33(6): 1113-1122．

[5] Liu F, Wang Y H, Wang S X. Development of the Hybrid Underwater Glider PetreI-II[J]. Sea Technology, 2014, 55(4): 51-54

[6] Yu J C, Zhang A Q, Jin W M, et al. Development and Experiments of the Sea-wing Underwater Glider[J]. China Ocean Engineering, 2011, 25(4): 721-736.

[7] Fan S S, Woolsey C. Elements of Underwater Glider Performance and Stability[J]. Marine Technology Society Journal, 2013, 47(3): 81-98.

[8] Bhatta P．Nonlinear Stability and Control of Gliding Vehicles[J]．Recent Developments In Ruminant Nutrition, 2006, 4(10): 325-352．

[9] Leonard N E, Graver J G. Model-based Feedback Control of Autonomous Underwater Gliders[J]. Oceanic Engine- ering, 2001, 26(4): 633-645.

[10] Breshears D D, Whicker J J, Johansen M P, et a1．Wind and Water Erosion and Transport in Semi-aridshrubland, Grassland and Forest Ecosystems: Quantifing Dominance of Horizontal Wind-driven Transport[J]．Earth Surface Pr- ocesses and Landforms, 2003, 28(11): 1189-1209．

[11] Hess J L．Calculation of Potential Flow about Arbitrary Three-dimensional Lifting Bodies[M]. LongBeach: Doug- las Aircraft Company, 1969．

Research on Hydrodynamic Characteristics and Gliding Performance of Underwater Glider “HaiXun 4000”

LIU Lai-lian, MIN Qiang-li, ZHAN Guang-ming

(The 710 Research Institute, China Shipbuilding Industry Corporation, Yichang 443003, China)

Compared with common underwater glider, the 4 000-meter underwater glider needs larger volume adjustment of its buoyancy control system in order to overcome the influence of seawater density change. Accurate dynamic model and hydrodynamic coefficient are the basis for realizing its control system design and accurate navigation. On the premise that its structural parameters are determined and the dynamic model is known, this paper aims to obtain accurate hydrodynamic coefficient and predict its performance, so as to provide control variables and optimization directions for the control system. Taking the 4 000-meter underwater glider “HaiXun 4000” as the research object, theoretical calculation, CFD simulation and experimental data analysis are combined to analyze the influence of the horizontal wing parameters on the gliding performance of the underwater glider, hence the hydrodynamic coefficient of constant motion in longitudinal plane is obtained, and the performance of the glider is predicted. Through the analysis of test data, the gliding speeds corresponding to different net buoyancy under typical gliding conditions are obtained and compared with the theoretical calculation to verify the correctness of the theoretical calculation. The results show that the error between actual pitch angle control and theoretical calculation is less than 5%, and the error between gliding speed and theoretical calculation under the actual net buoyancy of the test is within 15%, meeting the error requirements of 20% of the corresponding control amount for the control system. This research may provide a reference for further optimization and improvement of the control system.

underwater glider; gliding performance; hydrodynamic characteristic; gliding speed

TJ630; U674.941; O352

A

2096-3920(2019)05-0488-08

10.11993/j.issn.2096-3920.2019.05.003

劉來連, 閔強利, 張光明. “海鱘4000”水下滑翔機水動力特性與滑翔性能研究[J]. 水下無人系統(tǒng)學(xué)報, 2019, 27(5): 488-495.

2019-07-09;

2019-08-20.

青島海洋科學(xué)與技術(shù)國家實驗室“問海計劃”項目資助(2017WHZZB0101).

劉來連(1978-), 男, 高級工程師, 主要研究方向為水下機器人控制技術(shù).

(責(zé)任編輯: 楊力軍)