盧海洋,查 智
(中國船舶重工第七一〇研究所,湖北 宜昌 443003)
21世紀是海洋的世紀,在廣闊無垠的海洋空間,蘊藏著豐富的海洋生物資源、海洋礦物資源,是人類社會可持續(xù)發(fā)展的重要財富。水下機器人在海洋資源勘探和開發(fā)中的應用越來越廣泛,也發(fā)揮著重要的作用,具有廣泛的應用前景。對于大型或重型UUV運載平臺,其在無人干預遠程自主航行過程中,由于所處外部環(huán)境海水溫度、鹽度、深度的變化會引起海水密度的變化,導致UUV的浮力發(fā)生變化,對UUV的航行控制產(chǎn)生不利的影響,需要進行浮力補償。如果不對這種浮力變化作出補償則會使UUV可能沉入海底或浮于水面,影響作業(yè)任務的完成[1]。
圖1 油介質(zhì)型原理圖
通常,物體在水中實現(xiàn)沉浮動作通常有三種途徑:1)改變物體的體積而不改變質(zhì)量;2)改變物體的質(zhì)量而不改變體積;3)增加或減少所施加的外力,第一種和第二種方案的典型工作原理分別如圖1、圖2所示。第一種方案主要通過向皮囊中注/排油的方式來改變其體積大小,從而實現(xiàn)浮力調(diào)節(jié);第二種方案是通過向水艙中注/排水的方式來改變航行體質(zhì)量來實現(xiàn)浮力調(diào)節(jié);第三種方案是通過航行體自身的動力裝置來施加一個外力,從而實現(xiàn)浮力平衡作用[2-3]。
圖2 水介質(zhì)型原理圖
綜上所述,第一種方案,需要在UUV內(nèi)部存貯相當數(shù)量的液壓油,由于UUV內(nèi)部空間有限,采用此種方法將占用大量的空間,不利于UUV的整體設計。第三種方案,需要UUV自身攜帶的動力裝置來實現(xiàn)浮力的平衡,將消耗過多的能量,而UUV在水下航行時,所攜帶的能量是有限的。根據(jù)本文的研究背景,故采用第二種方案實現(xiàn)UUV的浮力調(diào)節(jié)功能[4-5]。
隨著水下機器人技術不斷發(fā)展和成熟,其在民用領域和軍用領域都將有越來越多的應用。近年來,以UUV為代表的各種無人作戰(zhàn)平臺在現(xiàn)代戰(zhàn)爭中扮演著重要的角色,成為兵力的倍增器,這一趨勢大大促進了人們對UUV技術的研究和應用。根據(jù)美國海軍發(fā)布的無人潛航器(UUV)總體規(guī)劃,待發(fā)展的UUV根據(jù)尺寸和重量分為便攜式航行器、輕型航行器、重型航行器、大型航行器四類,本文研究的浮力補償系統(tǒng),主要針對重型或大型航行器,解決其在遠程航行過程中的浮力補償問題[6]。
下面以中國南海11月份海水密度分布圖為例,分析UUV航行過程中浮力變化情況[7]。
圖3 11月份南海表層密度分布圖(多年平均)
圖4 11月份南海100 m層密度分布圖(多年平均)
11月份南海表層海水密度變化范圍基本在1.020~1.022 g/cm3之間,局部地區(qū)大于1.022 g/cm3但是小于1.0225 g/cm3,如圖3所示(圖中標示的是現(xiàn)場條件密度)。11月份南海100 m層海水密度變化范圍基本在1.022~1.024 g/cm3之間,局部地區(qū)大于1.024 g/cm3但是小于1.0245 g/cm3,如圖4所示(圖4中標示的是現(xiàn)場條件密度)。則在此海域內(nèi),0~100 m深度范圍內(nèi),海水密度最大變化量約為0.0045 g/cm3。
如果不考慮碼頭附近由于淡水的注入密度遠遠小于外海,按UUV排水量從小到大,依次計算其浮力變化值如表1所示。
表1 浮力變化值
由表1可知,若UUV排水量為10 t,海水密度變化為0.002 g/cm3時,UUV需要調(diào)節(jié)的浮力值為196 N。在南海航行時,UUV需要調(diào)節(jié)的浮力值為441 N,不管其為正浮力或負浮力,都將給UUV航行和安全帶來災難性的影響,所以必須對UUV的浮力變化采取相應的補償措施。
浮力補償水艙、調(diào)節(jié)裝置和控制器在UUV中的位置如圖5所示,在UUV的艏、艉段分別布置一個水艙,保證進行浮力補償時力矩平衡的作用。調(diào)節(jié)裝置和控制器布置在水艙旁的空艙內(nèi),浮力補償控制器采用主從方式,主控制器從UUV中央控制器接收指令,然后再分別向兩個從控制器發(fā)送控制指令,從而實現(xiàn)UUV的浮力補償[8]。
圖5 水艙位置圖
UUV浮力補償原理框圖如圖6所示,UUV在航行過程中,需要進行浮力補償時,浮力補償控制器根據(jù)當前海水密度和初始海水密度進行比較,計算其浮力變化值,然后對AUV前后水艙進行浮力分配,通過浮力調(diào)節(jié)裝置注/排水的方式實現(xiàn)UUV的浮力補償。
圖6 浮力補償原理框圖
浮力補償系統(tǒng)主要由浮力補償控制器、液位傳感器、浮力調(diào)節(jié)裝置和CTD傳感器等幾部分組成。浮力補償控制器根據(jù)CTD傳感器提供的海水溫度、鹽度、深度信息,通過海水狀態(tài)方程公式進行計算,得出當前UUV所處環(huán)境的海水密度。浮力補償控制器運用適當?shù)乃惴ǎㄟ^浮力調(diào)節(jié)裝置向UUV艙內(nèi)注/排水的方式來改變UUV自身的質(zhì)量,從而達到浮力補償?shù)哪康?,系統(tǒng)原理框圖如圖7所示[9]。
圖7 浮力補償系統(tǒng)原理框圖
由阿基米德原理可知,當一個物體全部浸在水中且處于平衡狀態(tài)時,流體作用在此物體上力的大小等于被此物體排開的水的重力,方向與重力相反,由于本文是研究UUV在水下航行的過程,所以其受到的浮力為[10]:
對于UUV在水下航行過程中,由于UUV所處外部環(huán)境海水密度變化所導致的浮力變化計算公式推導如下:
將式(2)與式(3)相除,得:
所以,UUV航行過程中由于海水密度變化導致浮力變化計算公式為:
式(5)中:ΔF表示浮力變化值;
Δρ表示海水密度的變化值;
ρ0表示淡水的密度;
m0表示UUV質(zhì)量。
海水密度是研究海洋學的一個重要參數(shù),海水密度不僅與海水的溫度、壓力和鹽度有關,而且還和純水的密度有關。本文根據(jù)80國際海水狀態(tài)方程(IESS-80)計算UUV所處環(huán)境海水密度,具體流程如圖8所示[11]。
浮力補償算法軟件流程圖如圖9所示。具體執(zhí)行過程如下:浮力補償算法根據(jù)CTD采集到的數(shù)據(jù)計算海水密度,并根據(jù)當前密度值與初始值之間的變化量Δρ與e1進行比較。當Δρ≥e1時,需要進行浮力調(diào)節(jié),并計算需要補償?shù)母×χ?,驅動浮力調(diào)節(jié)裝置工作,以實現(xiàn)浮力自適應之功能,浮力調(diào)節(jié)完成后將初始密度值置為當前密度值,即 ρ0=ρ1。當Δρ<e1時,不需要進行浮力調(diào)節(jié)。
圖8 海水密度計算軟件流程
圖9 浮力補償算法軟件流程
本文對浮力補償算法進行半實物仿真分析,CTD采集的海水溫度、鹽度、深度數(shù)據(jù)用計算機模擬輸入,計算出當前海水密度ρ1,與初始海水密度 ρ0進行比較,當 Δρ=|ρ1-ρ0|≥e1時,進行浮力補償, Δρ=|ρ1-ρ0|<e1時,不進行浮力補償。根據(jù)實際使用情況,本文取e1=1 kg/m3,UUV排水量V排=15 t,海水泵的流量Q=20 L/min。
根據(jù)UUV進行浮力補償時的實際應用情況,下面分別從Δρ<e1、Δρ≥e1兩種情況進行仿真分析。因為在UUV航行過程中,海水鹽度、溫度變化,對海水密度影響較大,所以在本文仿真試驗中主要改變海水鹽度和溫度。UUV在航行過程中,假設溫度、鹽度、深度變化情況如表2所示,UUV航行深度用壓強表示。
表2 海水密度計算值
從表2可知,當溫度變化Δt=4℃或鹽度變化ΔS=1時,對應表中序號分別為1和2、1和3,則海水密度變化分別為 Δρ=0.84 kg/m3<e1、Δρ=0.77 kg/m3<e1,水艙中水的體積仿真變化曲線如圖10、圖11所示。
當溫度變化Δt=5℃或鹽度變化ΔS=2時,對應表中序號分別為1和4、1和5,則海水密度變化分別為Δρ=1.059 kg/m3>e1、Δρ=1.542 kg/m3>e1。為了便于對比分析,仿真時將海水密度先由 ρ0增大 ρ1,再由 ρ1減小為 ρ0,看注入水量和排出水量是否相同,水艙中水的體積仿真變化曲線如圖12、圖13所示。
從圖10~圖13可知,當海水密度變化Δρ<e1=1 kg/m3時,浮力自適應控制系統(tǒng)沒有進行浮力補償工作;當海水密度變化Δρ≥e1=1 kg/m3時,海水密度先由 ρ0增大到 ρ1,再由 ρ1減小為 ρ0,浮力自適應控制系統(tǒng)進行浮力補償工作,并完成相應注排水工作,注入和排出水量相同,實現(xiàn)浮力補償功能。
圖10 溫度變化Δt=4℃
圖11 鹽度變化ΔS=1
圖12 溫度變化Δt=5℃
圖13 鹽度變化ΔS=2
本文針對重型、大型UUV在遠程航行過程的浮力補償問題進行研究,設計了一種基于CTD的浮力自動補償系統(tǒng),有效地解決了UUV在航行過程中浮力補償問題,為今后重型、大型UUV運載體的應用奠定了關鍵的技術基礎。
[1]陳強,張林根.美國軍用UUV現(xiàn)狀及發(fā)展趨勢分析[J].艦船科學技術,2010,32(7):129-134.
[2]錢東,孟慶國,薛蒙,等.美國海軍UUV的任務與能力需求[J].魚雷技術,2005,13(4):7-13.
[3]蘇紀蘭,袁業(yè)立.中國近海水文[M].北京:海洋出版社,2005.
[4]鄭榮,常海龍.浮力調(diào)節(jié)系統(tǒng)在作業(yè)型AUV上的應用研 究 [J].微 計 算 機 信 息 ,2006, 22 (9-2):207-208.
[5]楊鋼,郭晨冰,李寶仁.浮力調(diào)節(jié)裝置實驗研究[J].機床與液壓,2008,36(10):52-53.
[6]嚴安慶,方學紅,楊邦清.淺談潛水器浮力調(diào)節(jié)系統(tǒng)的研究現(xiàn)狀[J].水雷戰(zhàn)與艦船防護,2009,17(2):55-59.
[7]邱中梁,冷建興,陳建國,等.深海載人潛水器可調(diào)壓載系統(tǒng)研究[J].液壓與氣動,2003(11):9-11.
[8]盧海洋,林平.基于模糊PID控制的UUV定深控制系統(tǒng)設計[J].重慶理工大學學報,2011,25(214):58-62.
[9]趙文德,李建朋.基于浮力調(diào)節(jié)的AUV升沉運動控制技術[J].南京航空航天大學學報,2010,42(4):411-417.
[10]王恩華.阿基米德原理教學案例[J].物理教學討論,2007,25(295):17-18.
[11]姚勇,衛(wèi)薇,陳華秋.90國際溫標對80國際海水狀態(tài)方程影響的研究[J].海洋技術,2001,20(3):74-78.