陳佳倫,谷海濤,林揚(yáng),孫原,孟令帥
(1. 中國(guó)科學(xué)院沈陽(yáng)自動(dòng)化研究所機(jī)器人學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,遼寧 沈陽(yáng) 110016;2. 中國(guó)科學(xué)院機(jī)器人與智能制造創(chuàng)新研究院,遼寧 沈陽(yáng) 110169;3. 中國(guó)科學(xué)院大學(xué),北京 100049)
海洋探索的主要方式是母船攜帶AUV至指定位置,布放其進(jìn)行水下作業(yè),AUV完成任務(wù)后,人工回收至母船。母船式海洋作業(yè)流程經(jīng)過數(shù)十年的經(jīng)驗(yàn)累積,安全性、可靠性較高,但效率較低,人類已探索的海底只有5%,還有95%未知海底,海洋探索方式進(jìn)入一個(gè)瓶頸期。
基于目前情況,美國(guó)、日本提出下一步設(shè)計(jì)規(guī)劃,以及國(guó)際組織XPRIZE舉辦Shell Ocean Discovery(以下簡(jiǎn)稱海洋探測(cè)大賽),旨在提高海洋探索的自主性、規(guī)模、速度、深度和分辨率,突破目前瓶頸。伍茲霍爾海洋研究所提出了“未來海洋作業(yè)三步走發(fā)展規(guī)劃”[1],提出ASV(Autonomous Surface Vehicles,ASV)充當(dāng)母船通信的中繼信標(biāo)協(xié)調(diào)AUV完成大范圍、超視距作業(yè),逐步實(shí)現(xiàn)AUV長(zhǎng)時(shí)海洋探測(cè)工作。日本針對(duì)其周圍海域提出了“自盤古”計(jì)劃[2],使用ASV與水面艇協(xié)助AUV作業(yè),并提出可在海床移動(dòng)、采集樣品的ROV(Remote Operated Vehicle,ROV),完成研究礦床成因等科研任務(wù)。海洋探測(cè)大賽是一項(xiàng)為期3年、耗資700萬(wàn)美元的全球競(jìng)賽,日本GEBCO-NF團(tuán)隊(duì)采用無(wú)人艇與AUV協(xié)同作業(yè)方案,高效完成在有限條件下對(duì)規(guī)定區(qū)域探測(cè)并贏得比賽冠軍[3]。

圖 1 “自盤古”計(jì)劃示意圖Fig. 1Schematic illustration of the Zipangu project
本文針對(duì)未來海洋探測(cè),基于目前海洋技術(shù)水平,結(jié)合USV與AUV各自優(yōu)勢(shì),彌補(bǔ)AUV攜帶能源有限、USV探測(cè)深度有限等缺陷,提出USV與AUV一體化系統(tǒng)概念設(shè)計(jì),并針對(duì)一體化系統(tǒng)的需求提出一種三體船型USV。回收方案設(shè)計(jì)為USV水面動(dòng)態(tài)回收AUV,并根據(jù)水面回收流程,以某便攜式AUV為模型,研究AUV靜水面航行的水動(dòng)力特性分析,通過CFD軟件分析結(jié)果、水動(dòng)力簡(jiǎn)化模型理論推導(dǎo)以及外場(chǎng)水域試驗(yàn)數(shù)據(jù)對(duì)比,證明該CFD軟件分析的可靠性和準(zhǔn)確性。針對(duì)回收過程中USV尾流對(duì)AUV水面航行的影響進(jìn)行了15種工況下的對(duì)比仿真分析,得出AUV回收最佳工況。
USV與AUV一體化系統(tǒng)包括一對(duì)AUV與三體船型USV組合,USV具備運(yùn)輸、收放、補(bǔ)給、跟蹤、通信、定位AUV的能力,在水下一定深度、一定面積進(jìn)行USV與AUV協(xié)同作業(yè)。一體化系統(tǒng)作業(yè)模式包括定位作業(yè)流程、收放作業(yè)流程、通信作業(yè)流程、運(yùn)載作業(yè)流程、艇載數(shù)據(jù)處理作業(yè)流程、補(bǔ)給作業(yè)流程。一體化系統(tǒng)適應(yīng)于深遠(yuǎn)海的自主作業(yè),作業(yè)期間無(wú)需母船直接協(xié)助,降低作業(yè)成本、提高作業(yè)效率。

圖 2 三體船示意圖Fig. 2Schematic illustration of the trimaran
收放作業(yè)流程分為AUV回收作業(yè)流程和AUV布放作業(yè)流程。
1)回收作業(yè)流程
AUV在超短基線定位系統(tǒng)導(dǎo)航下歸航至進(jìn)入無(wú)線通信引導(dǎo),AUV在無(wú)線通信下與USV同向航行,通過速度差進(jìn)入U(xiǎn)SV主艙體;當(dāng)AUV靠近USV尾部,導(dǎo)向板微調(diào)AUV前行方向,AUV首部觸碰傳送帶,傳送帶上表面滾動(dòng)方向朝向船首,AUV緊貼傳送帶進(jìn)入主艙體,AUV首部碰觸限位裝置限位,限位信號(hào)經(jīng)由艇載微型處理器處理,下達(dá)指令1傳送帶停止工作、下達(dá)指令2經(jīng)由無(wú)線通信至AUV停止推進(jìn)器工作、下達(dá)指令3夾緊裝置固定AUV,完成收放作業(yè)流程回收作業(yè)。
2)AUV布放作業(yè)流程
夾緊裝置停止工作,傳送帶上表面滾動(dòng)方向朝向三體船船尾,AUV緊貼傳送帶進(jìn)入水域,完成收放作業(yè)流程布放作業(yè)。

圖 3 水面回收?qǐng)鼍皥DFig. 3Schematic illustration of the recovery scenario on the water
針對(duì)一體化系統(tǒng)AUV回收采取水面回收,相對(duì)于水下回收方式,其通信環(huán)境較惡劣、控制較復(fù)雜、回收裝置限定條件較多[4-7]。水面回收優(yōu)勢(shì)是將水下三維對(duì)點(diǎn)轉(zhuǎn)換為水面二維回收,并采取無(wú)線導(dǎo)引,增大回收成功概率。針對(duì)一體化作業(yè)要求多次收放AUV,水面回收對(duì)AUV外形無(wú)需改動(dòng),減少AUV額外功耗,USV尾部借鑒登陸艦尾部斜坡設(shè)計(jì),便于運(yùn)放AUV。
動(dòng)態(tài)回收考慮USV在水面無(wú)法固定于一點(diǎn)、易受擾動(dòng),降低回收成功概率,故USV水面直線航行保證方向穩(wěn)定性,AUV通過無(wú)線導(dǎo)引方式以較快于USV的速度同方向追趕,通過速度差實(shí)現(xiàn)進(jìn)入U(xiǎn)SV后艙,提高回收成功概率。兩者對(duì)接相對(duì)速度較小,同時(shí)減小碰撞力。水面回收的缺點(diǎn)是海浪對(duì)其影響較大,故便攜式AUV回收應(yīng)在低海況下進(jìn)行,仿真環(huán)境設(shè)置為靜水面環(huán)境。
定位作業(yè)流程包括中國(guó)北斗衛(wèi)星導(dǎo)航系統(tǒng)定位USV以及USV超短基線定位AUV。通信作業(yè)流程包括水下通信、水面通信及衛(wèi)星通信3部分,其中水下通信使用水聲通信機(jī)完成信息傳遞,水面通信利用無(wú)線通信設(shè)備形成通信網(wǎng)絡(luò),衛(wèi)星通信使用北斗通信系統(tǒng),三者形成通信鏈路。艇載數(shù)據(jù)處理作業(yè)流程用于對(duì)USV接收數(shù)據(jù)進(jìn)行相關(guān)分析得出下一步運(yùn)行指令;運(yùn)載系統(tǒng)是指AUV固定于USV主艙體,USV運(yùn)輸AUV至目標(biāo)區(qū)域;補(bǔ)給作業(yè)流程對(duì)位于USV主艙體內(nèi)的AUV進(jìn)行無(wú)線充電。

圖 4 通信作業(yè)和定位作業(yè)場(chǎng)景圖Fig. 4Schematic illustration of the communication and location work processes
一體化系統(tǒng)作業(yè)模式作業(yè)流程為:AUV固定于USV主艙體,USV離開基站或母船,USV與基站或母船時(shí)刻通過無(wú)線通信或衛(wèi)星通信保持?jǐn)?shù)據(jù)、指令傳輸,USV自身位置通過定位作業(yè)流程北斗衛(wèi)星定位;運(yùn)輸AUV至指定作業(yè)區(qū)域,USV收放作業(yè)流程布放AUV,AUV開始位于一定深度水域作業(yè),無(wú)人水面艇超短基線定位作業(yè)流程定位AUV,USV通過水下通信相互傳輸信息至AUV,完成水下一定深度作業(yè)后,AUV上浮水面,USV收放裝置回收AUV,AUV固定于USV主艙體,USV補(bǔ)給作業(yè)流程補(bǔ)給AUV,AUV通過無(wú)線數(shù)據(jù)傳輸系統(tǒng)傳輸原始數(shù)據(jù),艇載數(shù)據(jù)處理作業(yè)流程處理原始數(shù)據(jù),USV通過無(wú)線數(shù)據(jù)傳輸系統(tǒng)傳輸指令至AUV,完成一次作業(yè)流程,數(shù)次循環(huán)直至完成任務(wù)要求,USV返回基站或母船。
針對(duì)USV水面回收AUV步驟,第1步AUV需要穩(wěn)定在水面航行,以某便攜式AUV為模型,通過STARCCM+水動(dòng)力軟件構(gòu)建仿真模型。AUV質(zhì)量75 kg,長(zhǎng)度1776 mm,直徑250 mm,靜浮力0.5 kg,浮心距首部830 mm。
AUV尾舵采用NACA0012翼型,弦長(zhǎng)115 mm,展長(zhǎng)110 mm,展弦比0.96。

圖 6 AUV垂直面動(dòng)網(wǎng)格Fig. 6AUV moving grid on vertical plane
AUV仿真流域設(shè)定為45 m×4 m×5 m,本文在計(jì)算域內(nèi)設(shè)置雙重網(wǎng)格加密區(qū),以提高體網(wǎng)格質(zhì)量,在保證計(jì)算精度前提下,提高仿真環(huán)境模型計(jì)算效率。仿真流域網(wǎng)格基礎(chǔ)值0.20 m;第1重自由液面加密區(qū)尺度為45 m×4 m×3 m,基礎(chǔ)值0.06 m;第2重AUV計(jì)算加密區(qū)為3 m×1 m×1 m,基礎(chǔ)值0.03 m。
根據(jù)之前試驗(yàn)情況以及以往工作經(jīng)驗(yàn),AUV穩(wěn)定在水面航行必須保持上浮舵狀態(tài),否則AUV首部將呈現(xiàn)大幅震蕩狀態(tài),不滿足實(shí)際動(dòng)態(tài)水面回收方案對(duì)AUV的航行要求,故仿真過程中需對(duì)不同舵角進(jìn)行仿真計(jì)算,尋求最佳航行狀態(tài)。選取舵角10°、推進(jìn)力20 N為例。
仿真結(jié)果看,AUV仿真計(jì)算收斂后,最后速度和縱傾角穩(wěn)定在一定數(shù)值范圍內(nèi),水平前行運(yùn)動(dòng),AUV的縱傾角與沖角相等,速度水平向前,符合實(shí)際航行現(xiàn)象。

圖 7 AUV縱傾角、速度變化曲線Fig. 7Variation curves of pitch and velocity of AUV

圖 8 AUV水面穩(wěn)定航行Fig. 8AUV surface stability navigation
選取舵角 5°,10°,15°,20°,25°五種工況分析,對(duì)AUV尾部處軸線方向施加5~100 N推進(jìn)力,多種工況進(jìn)行仿真計(jì)算,結(jié)果如圖9和圖10。其中A,B,C,D 分別表示舵角 10°,15°,20°,25°。

圖 9 AUV舵角與縱傾角關(guān)系曲線Fig. 9Curve of relationship between rudder angle and pitch angle of AUV

圖 10 AUV推進(jìn)力與穩(wěn)定速度關(guān)系曲線Fig. 10Curve of relationship between rudder angle and constant velocity of AUV
在相等舵角工況下,縱傾角隨推進(jìn)力增加先增大至平穩(wěn)再減小,舵角越大縱傾角越大,大致呈現(xiàn)線形關(guān)系,AUV極限速度也隨著舵角增大而增大,平均速度1.665 m/s。
本節(jié)主要根據(jù)仿真數(shù)據(jù)相關(guān)結(jié)論,構(gòu)造AUV垂直面水動(dòng)力模型,并進(jìn)行水動(dòng)力計(jì)算推導(dǎo),并與仿真結(jié)果進(jìn)行對(duì)比,驗(yàn)證仿真結(jié)果可靠性。
動(dòng)系原點(diǎn)與AUV重心重合,AUV重心位于浮心正下2 mm,如圖11所示。

圖 11 AUV運(yùn)動(dòng)坐標(biāo)系Fig. 11The moving coordinate system of AUV
AUV一般運(yùn)動(dòng)方程建模,問題只考慮垂直面運(yùn)動(dòng)求解。


式中:X0為AUV直線航行阻力;Z0為AUV上下不對(duì)稱導(dǎo)致z軸方向力;M0為AUV前后不對(duì)稱導(dǎo)致y軸力矩;Xu,Zw,Mw為速度系數(shù);Zq,Mq為角速度系數(shù);,為舵角系數(shù)。

由靜力學(xué)可知,當(dāng)AUV俯仰運(yùn)動(dòng)時(shí)導(dǎo)致重力與浮力作用不在同一鉛垂線上構(gòu)成浮力扶正力矩,其中h為AUV穩(wěn)心高,g為重力加速度。



設(shè)定AUV弱機(jī)動(dòng)運(yùn)動(dòng),故忽略X方向等式,并代入水動(dòng)力式,計(jì)及剩余靜載、扶正力矩和槳力矩,并考慮AUV前后不對(duì)稱時(shí)的垂直面操縱運(yùn)動(dòng)方程式:


上式水動(dòng)力導(dǎo)數(shù)和運(yùn)動(dòng)參數(shù)都是有因次項(xiàng),不便于比較,對(duì)上式進(jìn)行無(wú)因次化,力遍除力矩遍除
上式相應(yīng)無(wú)因次形式為:

通過仿真不同潛深的拘束船模試驗(yàn)測(cè)定,考慮自由水面對(duì)水動(dòng)力系數(shù)的影響[8],校正后:

AUV水面俯仰定常運(yùn)動(dòng)時(shí),可看作帶縱傾角定深航行,定深深度近似為0,,故速度矢方向在水平方向上,由公式可知。

代入相應(yīng)數(shù)值,可得:

對(duì)比仿真與計(jì)算數(shù)值,以及后續(xù)的外場(chǎng)水域試驗(yàn)數(shù)據(jù),得出下表。
通過水動(dòng)力計(jì)算縱傾角與仿真縱傾角對(duì)比,仿真穩(wěn)定數(shù)值與水動(dòng)力計(jì)算數(shù)值相差較小,仿真結(jié)果AUV穩(wěn)定航行時(shí)縱傾角與舵角呈線形關(guān)系,縱傾角與AUV外形結(jié)構(gòu)以及靜浮力相關(guān)。
求解AUV不同舵角下的極限速度,由公式化簡(jiǎn)可得:

將其中與速度相關(guān)水動(dòng)力系數(shù)代入u計(jì)算:


化簡(jiǎn):

代入公式:
在外場(chǎng)水域進(jìn)行實(shí)際試驗(yàn),AUV上浮下潛舵機(jī)轉(zhuǎn)動(dòng)角度范圍在-20°~+20°;AUV內(nèi)部安裝有TCM5電子羅盤,時(shí)刻記錄AUV縱傾角數(shù)值。AUV底部安裝有DVL多普勒測(cè)速儀,時(shí)刻測(cè)量AUV速度。縱傾角、速度數(shù)值通過無(wú)線網(wǎng)線發(fā)送至岸基,試驗(yàn)人員可時(shí)刻觀察AUV航行情況,保證試驗(yàn)安全進(jìn)行。
考慮離岸水域水深較淺,為避免AUV下潛可能性,AUV從岸基出發(fā)時(shí)速度控制在1 m/s、上浮舵20°,監(jiān)控AUV實(shí)際縱傾角數(shù)值,將AUV速度控制在平均速度1 m/s、均方差0.10 m/s、航行距離約50 m范圍內(nèi),縱傾角穩(wěn)定平均數(shù)值為11.0°,按照此步驟以此完成上浮舵角15°,10°,5°試驗(yàn)航行,平均縱傾角為7.4°,4.2°,在上浮舵為5°航行時(shí),AUV航行狀況不穩(wěn)定。
將實(shí)際航行數(shù)據(jù)與計(jì)算縱傾角、仿真縱傾角列于表1中,對(duì)比三者10°,15°,20°三種工況數(shù)據(jù),證實(shí)了仿真模型以及水動(dòng)力簡(jiǎn)化模型的可靠性與準(zhǔn)確性,為第二步AUV追逐USV的仿真環(huán)境搭建提供基礎(chǔ)。

表 1 AUV縱傾角對(duì)比表Tab. 1Contrast of AUV pitch angles
針對(duì)一體化系統(tǒng)對(duì)USV外形需求,以某便攜式AUV為基準(zhǔn),選取三體船外形作為一體化系統(tǒng)USV。
三體船相比單體船提高了船舶的穩(wěn)定性和耐波性,可將推進(jìn)器安裝至2個(gè)側(cè)體上,提高三體船的機(jī)動(dòng)性,并考慮在回收AUV時(shí),減小螺旋槳尾流對(duì)AUV航行影響,以及當(dāng)回收失敗時(shí)避免AUV撞擊螺旋槳。相比雙體船,三體船的主船體可提供運(yùn)輸空間,在船體尾部設(shè)計(jì)AUV回收裝置。
USV的推進(jìn)器布置在2個(gè)側(cè)體上且低速航行,故數(shù)值仿真時(shí)可忽略螺旋槳尾流對(duì)AUV航行特性的影響,但航行過程中,不能忽略USV航行時(shí)尾流對(duì)AUV航行特性的影響。通過CFD仿真,對(duì)其影響進(jìn)行判斷。
追逐仿真流域設(shè)定為40 m×16 m×8 m,并在計(jì)算域內(nèi)針對(duì)AUV和USV設(shè)置雙重網(wǎng)格加密區(qū)。仿真流域網(wǎng)格基礎(chǔ)值0.20 m;第一重自由液面加密區(qū)尺度為40 m×16 m×1 m,基礎(chǔ)值0.06 m;第二重AUV加密區(qū)為3 m×1 m×1 m,基礎(chǔ)值0.03 m,USV加密區(qū)10 m×5 m×1 m,基礎(chǔ)值0.03 m。
USV仿真過程中保持速度約為0.5,0.75,1.00,1.25,1.50 m/s五種工況,AUV距離USV尾部20 m時(shí)兩者同時(shí)前向運(yùn)動(dòng),AUV工況選取舵角25°,仿真收斂速度1.75 m/s、縱傾角7.8°,AUV追逐USV過程中速度、縱傾角變化如圖15所示,仿真停止條件是AUV碰觸USV。

圖 12 AUV外場(chǎng)航行數(shù)據(jù)Fig. 12AUV outfield navigation data

圖 13 三體船型USV模型Fig. 13Model of trimaran USV

表 2 三體船模型參數(shù)表Tab. 2Parameter table of trimaran USV model
通過不同工況下仿真結(jié)果看,AUV在逐漸靠近USV過程中,AUV速度和縱傾角均發(fā)生變化。
1)靠近過程中AUV速度降低、縱傾角增大,速度差越近越明顯;
2)AUV即將碰觸USV尾部時(shí),速度和縱傾角會(huì)發(fā)生突增,速度差越近越明顯。
由于USV速度越快產(chǎn)生尾流對(duì)AUV影響越大,尾流擾動(dòng)使AUV縱傾角增大速度降低。靠近USV時(shí)AUV速度突增,原因是USV后艙內(nèi)壁產(chǎn)生正伴流現(xiàn)象,即其帶動(dòng)附近水域有向前運(yùn)動(dòng)趨勢(shì),導(dǎo)致流速相對(duì)于AUV速度差減小,AUV粘性水阻力隨之減小,USV速度越快越明顯,并且速度差較小時(shí),AUV處在尾部區(qū)域相對(duì)時(shí)間也增加,速度突增幅度越大,2個(gè)原因?qū)е滤俣炔钤叫UV靠近USV速度突增越明顯。

圖 14 水面回收仿真場(chǎng)景計(jì)算域Fig. 14Calculation domain of docking simulation on the water

圖 15 回收過程AUV速度1.75 m/s仿真航行數(shù)據(jù)Fig. 15AUV navigation data in docking simulation process at speed of 1.75 m/s
由于縱傾角增加至約16°,回收過程中AUV底部有較大可能性率先碰觸USV斜坡,由于AUV底部布置傳感器等探測(cè)裝置,對(duì)其將會(huì)產(chǎn)生損害,且考慮回收時(shí)AUV速度突增情況,故針對(duì)便攜式AUV不建議使用AUV速度1.75 m/s工況回收。

圖 16 回收過程仿真水面場(chǎng)景Fig. 16Surface water scene of recovery simulation process

圖 17 水面速度矢量圖Fig. 17Velocity vector map on the water
USV仿真過程中保持速度約為0.5,0.75,1.00,1.25 m/s四種工況,AUV距離USV尾部15 m時(shí)兩者同時(shí)前進(jìn)方向運(yùn)動(dòng),AUV工況選取舵角20°,仿真收斂速度1.50 m/s,縱傾角8.7°,AUV追逐USV過程中速度、縱傾角變化如圖18所示,仿真停止條件是AUV碰觸USV。
分析結(jié)果與AUV1.75 m/s大致相同,AUV碰觸USV時(shí)縱傾角較理想;USV在0.75 m/s,0.5 m/s速度突增較小,在此基礎(chǔ)上USV速度越大,方向穩(wěn)定性越強(qiáng),USV工況在0.75 m/s時(shí)仿真結(jié)果針對(duì)一體化水面回收方案較理想。
鑒于上述仿真結(jié)果,相對(duì)速度差過小,AUV速度、縱傾角均有較大影響,對(duì)于舵角25°縱傾角過大,不適用于水面回收狀態(tài)。針對(duì)AUV 1.25 m/s采用舵角10°,15°,20°三種舵角,USV 速度 0.5,0.75 m/s六種工況進(jìn)行仿真對(duì)比,如圖19所示。
USV同速情況下,AUV回收時(shí)速度、縱傾角變化趨勢(shì)、幅度基本相同,綜合考慮前2種工況,優(yōu)先選擇縱傾角較小的20°舵角;USV速度0.75,0.5 m/s影響差別不大,故選取0.75 m/s。

圖 18 回收過程AUV速度1.50 m/s仿真航行數(shù)據(jù)Fig. 18AUV navigation data in docking simulation process at speed of 1.50 m/s
本文提出一種基于AUV與USV的一體化作業(yè)系統(tǒng),面向該系統(tǒng)設(shè)計(jì)USV水面回收AUV作業(yè)流程,對(duì)回收過程中AUV的水面穩(wěn)定航行特性進(jìn)行了深入分析,USV航行尾流對(duì)其后方AUV航行在垂直面上三自由度的影響進(jìn)行仿真分析。通過仿真計(jì)算及外場(chǎng)試驗(yàn)得出如下結(jié)論:
1)以某便攜式AUV為模型進(jìn)行水面垂直面三自由度航行水動(dòng)力仿真,通過使用上浮舵角使AUV可在水面穩(wěn)定水平前向航行,AUV的縱傾角與沖角相等,符合實(shí)際航行現(xiàn)象。
2)建立AUV水面航行水動(dòng)力模型,求解出AUV舵角與AUV水面航行縱傾角呈線形關(guān)系,經(jīng)過外場(chǎng)試驗(yàn)驗(yàn)證了水動(dòng)力模型和仿真結(jié)果的準(zhǔn)確性和可靠性。
3)以三體船型USV為仿真模型,進(jìn)行動(dòng)回收回收AUV水動(dòng)力仿真分析,通過15種工況對(duì)比分析USV航行尾流對(duì)較高速AUV航行的垂直面三自由度影響,得出USV速度0.75 m/s航行、AUV舵角20°速度1.25~1.5 m/s回收可靠性最佳。

圖 19 回收過程AUV速度1.25 m/s仿真航行數(shù)據(jù)Fig. 19AUV navigation data in docking simulation process at speed of 1.25 m/s