劉洪生，王新海，陳祖瑞，于 俊

(1. 海軍駐431廠軍事代表室，遼寧 葫蘆島 125004；2. 上海船舶設(shè)備研究所，上海 200031)

0 引 言

水下機器人在海洋工程和深海探索中運用廣泛，但水下機器人轉(zhuǎn)向力不足，抗風浪擺蕩能力差，不易保持航行方向和姿態(tài)穩(wěn)定。為水下機器人增加矢量推進裝置，并通過算法控制液壓驅(qū)動矢量推進裝置，能按要求調(diào)整船體姿態(tài)，為了提升整機系統(tǒng)的響應速度、準確度等性能指標，需對液壓系統(tǒng)的參數(shù)進行優(yōu)化[1–2]。由于液壓系統(tǒng)存在效率損失，因此會造成機械系統(tǒng)的運動與預期產(chǎn)生誤差，而且這些損失隨著液壓泵、液壓馬達、液壓缸排量、液壓元件參數(shù)以及液壓系統(tǒng)壓力的不同而不同，不同參數(shù)組合得到的運動控制效果也不盡相同，而遺傳算法在解決這類NP問題、非線性、多目標和多峰值問題方面效果顯著[3–4]?？死锝鸱ㄔ趨?shù)聯(lián)系緊密且數(shù)量確定的情況下可以快速模擬響應面，進而得到優(yōu)化值。因此可采用遺傳算法和克里金法進行液壓系統(tǒng)的參數(shù)優(yōu)化。

1 水下機器人矢量推進液壓液壓驅(qū)動系統(tǒng)

為分析液壓參數(shù)對矢量推進裝置的運動影響，首先要考慮矢量推進裝置的運動過程及其液壓驅(qū)動回路。矢量推進裝置要完成升降、回轉(zhuǎn)和鎖緊3個動作。3個液壓驅(qū)動系統(tǒng)并聯(lián)，故可分3部分分別進行優(yōu)化設(shè)計。若對每個回路的所有液壓參數(shù)進行優(yōu)化設(shè)計，費時費力，且不易發(fā)現(xiàn)關(guān)鍵參數(shù)的作用。故可選取液壓系統(tǒng)中關(guān)鍵元器件的參數(shù)進行優(yōu)化設(shè)計。其關(guān)鍵元器件參數(shù)選取如圖1～圖3所示。

2 理論基礎(chǔ)

在保證完成既定動作的前提下需要協(xié)調(diào)好液壓系統(tǒng)中的各主要參數(shù)，使得系統(tǒng)具有良好的作業(yè)性能、可靠性和高效性。因此需要進行液壓系統(tǒng)參數(shù)優(yōu)化。遺傳算法相比實驗統(tǒng)計方法不用求解系統(tǒng)的近似模型，不必考慮試驗點的選取問題，可以給出各因素較大的取值范圍?？死锝鸱ㄔ趨?shù)數(shù)量確定情況下可以快速靈活對模型進行預測，因此液壓參數(shù)部分采用遺傳算法進行優(yōu)化設(shè)計，液壓控制模塊的PID參數(shù)采用克里金算法。

2.1 遺傳算法

遺傳算法（Genetic Algorithm，GA）是以自然選擇和遺傳理論為基礎(chǔ)，將生物進化過程中適者生存規(guī)則與群體內(nèi)部染色體的隨機信息交換機制相結(jié)合的高效全局尋優(yōu)搜索算法。遺傳算法工具箱采用Sheffield遺傳算法工具箱[5–7]。具體步驟為：

步驟1編碼與解碼，確定待優(yōu)化液壓參數(shù)范圍，將其寫成二進制形式所需的位數(shù)。通過產(chǎn)生隨機數(shù)來生成參數(shù)對應的二進制編碼，并將所有二進制編碼首尾連接成一條“染色體”，即為一個個體，多個個體形成種群。二進制解碼公式如下：

步驟2適應度計算，將隨機數(shù)產(chǎn)生的液壓參數(shù)帶入模型并仿真，將實際位移曲線和期望位移曲面的差值的均方根植最小作為控制目標，即適應度函數(shù)為：

步驟3設(shè)定迭代步、種群值、選擇算子、交叉和變異算子等初始參數(shù)，代入Sheffield工具箱模板進行計算，求得參數(shù)優(yōu)化值。

2.2 克里金法

克里金法（Kriging）是一項實用的空間估計技術(shù)，其由一個參數(shù)模型和一個非參數(shù)隨機過程聯(lián)合構(gòu)成?？死锝鸱ū葐蝹€的參數(shù)化模型具有更強的靈活性和預測能力，同時又克服了非參數(shù)化模型處理高維數(shù)據(jù)的局限性?？死锝鹉Ｐ鸵话惆?部分：多項式和隨機分布[8–9]。具體模型為：

拉丁超立方抽樣，該方法能避免抽樣點在小鄰域內(nèi)重合?；驹砣缦拢喝绻M行n次抽樣，把m個隨機變量都分成等概率的n個區(qū)間，故整個抽樣空間被分成等概率的個小格子。n次抽樣一定分別落在每個小區(qū)間中，故實際得到的抽樣點等概率地分散在整個隨機空間中。具體步驟如下：

步驟1樣本點選取，采用拉丁超立方抽樣，設(shè)計空間為各液壓元件參數(shù)取值的范圍。

步驟2系統(tǒng)響應計算，在選取個樣本點后，將液壓元件設(shè)為樣本點參數(shù)值，進行仿真計算。將機構(gòu)的期望位移與實際位移差值的均方根值作為系統(tǒng)的響應，即

步驟3將樣本點和響應值作為已知信息Ysampling；選擇二次函數(shù)作為擬合函數(shù)和高斯相關(guān)方程；將待優(yōu)化變量取值范圍均分網(wǎng)格，各變量交叉形成節(jié)點。利用克里金插值法預測各節(jié)點響應值，并利用最小值函數(shù)求均方根值最小處變量值，即優(yōu)化值。

3 優(yōu)化設(shè)計實驗

3.1 升降回路液壓參數(shù)優(yōu)化設(shè)計

升降部分的運動為開環(huán)控制，抗干擾能力弱，控制精度較低，無自動糾偏能力。升降運動過程固定，各液壓元件的壓力流量可以根據(jù)實際要求計算得出，但元件參數(shù)較多，計算復雜，可以通過液壓元件參數(shù)的匹配使得升降部分的運動符合預期的運動。

選擇圖1中4個元器件液壓參數(shù)進行參數(shù)優(yōu)化優(yōu)化設(shè)計，如表1所示。優(yōu)化算法采用遺傳算法，期望運動曲線和優(yōu)化曲線如圖4所示。從圖中可以看到升降部分的實際位移曲線與期望曲線吻合較好，兩者的最大誤差小于2 mm，控制在運動誤差允許范圍內(nèi)。

表1 升降運動液壓參數(shù)優(yōu)化范圍和結(jié)果Tab.1 Optimization range and results of hydraulic parameters for lifting motion

3.2 鎖緊回路液壓參數(shù)優(yōu)化設(shè)計

鎖緊回路的運動也是開環(huán)控制，選擇圖2中4個液壓元件，如表2所示。優(yōu)化算法采用遺傳算法，期望運動曲線和優(yōu)化運動曲線如圖5所示。

表2 鎖緊運動液壓參數(shù)優(yōu)化范圍和結(jié)果Tab.2 Optimization range and results of hydraulic parameters for locking motion

從圖中可看到實際位移曲線與期望曲線并沒有完全重合，這是因為優(yōu)化參數(shù)初始范圍給定不合理，最優(yōu)值在初始邊界上獲得，可通過調(diào)整優(yōu)化范圍來獲得更優(yōu)的優(yōu)化值。但是在此回路中鎖緊銷實際能夠在響應時間達到相應的位置，符合機構(gòu)運動的要求。

3.3 回轉(zhuǎn)回路液壓參數(shù)優(yōu)化設(shè)計

回轉(zhuǎn)部分采用PID閉環(huán)控制，液壓參數(shù)選擇圖3中4個液壓元件，如表3所示。優(yōu)化算法采用遺傳算法，期望運動曲線和優(yōu)化曲線如圖6所示。

表3 回轉(zhuǎn)運動液壓參數(shù)優(yōu)化范圍和結(jié)果Tab.3 Optimization range and results of hydraulic parameters for rotary motion

從圖中可以看到回轉(zhuǎn)回路的優(yōu)化運動曲線與期望曲線重合度高，最大誤差角度不超過0.9°。運動控制符合實際要求。

3.4 回轉(zhuǎn)回路PID控制參數(shù)優(yōu)化設(shè)計

回轉(zhuǎn)回路PID參數(shù)數(shù)量較少且聯(lián)系緊密，可采用克里金法進行優(yōu)化。優(yōu)化參數(shù)范圍及結(jié)果如表4所示。

表4 回轉(zhuǎn)回路PID參數(shù)優(yōu)化范圍及結(jié)果Tab.4 Optimization range and results of PID of rotary circuit

其中可以看出I相比P和D范圍變化較小，選擇P和D作克里金曲面如圖7所示。從圖中可看到在D=6.020 4，P=0.142 9處，曲面有一個凹處，即在該處取得最優(yōu)值。

4 結(jié) 語

矢量推進裝置可以增加水下機器人任意工況下的轉(zhuǎn)向力和抗風浪擺蕩能力，有利于保持航行方向穩(wěn)定和姿態(tài)穩(wěn)定。通過對矢量推進裝置液壓系統(tǒng)的參數(shù)進行優(yōu)化，可以提升整機系統(tǒng)的響應速度、準確度等性能指標，并改善部分元器件的工作狀態(tài)。