龐愛平　何朕

摘 要：針對水下高速超空泡航行體尾部與空泡壁周期性碰撞的現(xiàn)象，提出了采用擾動觀測器來補償和消除重力的影響。推導(dǎo)了擾動觀測器方程，給出了補償器設(shè)計方案。仿真驗證了擾動補償能夠避免超空泡航行體前進中周期性碰撞空泡壁的現(xiàn)象，并且控制面偏轉(zhuǎn)角滿足正常工作范圍的要求。接入擾動估計和補償控制后，航行體由原來的周期性拍打空泡壁能很快恢復(fù)到平穩(wěn)工作狀態(tài)。結(jié)論是擾動觀測補償?shù)脑O(shè)計能夠避免超空泡航行體產(chǎn)生滑行，使運動更為平穩(wěn)，并且減少了前進的阻力。所提的擾動估計和補償?shù)目刂扑枷耄苊饬酥苯友a償時執(zhí)行機構(gòu)行程過大的問題，可用于一般的控制系統(tǒng)設(shè)計。

關(guān)鍵詞：超空泡；超空泡航行體；擾動觀測器；補償；滑行

中圖分類號：TP 273

文獻標(biāo)志碼：A

文章編號：1007-449X（2018）01-0107-07

0 引 言

水下航行體因受到水的阻力作用而極大限制了其速度。超空泡技術(shù)的應(yīng)用突破了水下航行體速度的瓶頸，提高了水下武器的戰(zhàn)斗力，具有重大的科研意義和軍事應(yīng)用前景[1]。已成為國內(nèi)外學(xué)者關(guān)注的熱點[2-4]。

超空泡航行體的流體動力特性極其復(fù)雜，俄羅斯和烏克蘭首先對超空泡航行體的建模與控制開展研究。烏克蘭學(xué)者Savchenko[5]首先研究了航行體與空泡的作用關(guān)系，仿真分析超空泡航行體的開環(huán)運動，在開環(huán)狀態(tài)下，航行體無法維持運動，必須采用主動控制來保證航行體的運動穩(wěn)定。Kirschner等[6]提出可以采用尾翼和空化器作為控制面來聯(lián)合控制航行體的運動穩(wěn)定性，設(shè)計了LQR前饋反饋控制，可以維持航行體的運動穩(wěn)定性。俄羅斯學(xué)者Dzielski等[7]在研究中提出了一種超空泡航行體的標(biāo)準(zhǔn)設(shè)計，建立了超空泡航行體的基準(zhǔn)模型并給定相應(yīng)參數(shù)，該模型既相對簡化，又保留了超空泡航行體的基本特性，被之后的研究者所廣泛采用[7-9]。Dzielski在文獻[7]中還給出了一個線性狀態(tài)反饋控制律，在這個控制律的作用下，航行體尾部周期性地拍打著空泡下壁高速前行。在Dzielski的研究基礎(chǔ)上，Lin等在文獻[9]中討論了反饋增益對航行體極限環(huán)振蕩的影響。

由于航行體被空泡包裹，失去水的浮力作用，在重力的作用下自然下沉，當(dāng)尾部接觸并刺穿空泡壁時產(chǎn)生反彈的滑行力，所產(chǎn)生的滑行力又將航行體推回空泡內(nèi)，航行體受到重力作用再次下沉而與空泡壁再次碰撞，如此反復(fù)形成周期性的滑水現(xiàn)象。雖然這種周期性滑水的極限環(huán)穩(wěn)態(tài)運行模式也是超空泡航行體一種典型的工作模式，但由于空泡的內(nèi)徑有限，當(dāng)滑行力這個擾動外力過大時，會將航行體推向空泡的另一側(cè)，來回沖撞導(dǎo)致航行體最終失去穩(wěn)定，這就要求控制設(shè)計具有良好的抗擾動能力。此外航行體與空泡壁頻繁的碰撞也產(chǎn)生阻力，消耗能量，不但影響前進速度，還有可能破壞空泡的穩(wěn)定性，故良好的工作狀態(tài)應(yīng)該是避免滑行力的產(chǎn)生，消除這種周期性滑水現(xiàn)象。為了消除重力的影響，文獻[7]和文獻[9]中都采用反饋線性化或直接對消的方法設(shè)計了避免滑行力產(chǎn)生的控制方案，但是其設(shè)計的結(jié)果所要求的控制面偏轉(zhuǎn)角過大，達到了1.5 rad，顯然超出了正常的工作范圍，因而無法得到實際應(yīng)用。此外反饋線性化還要求精確的數(shù)學(xué)模型，這在實際中也是難于實現(xiàn)的。

這里需要指出的是，對于超空泡航行體來說，需要補償?shù)牟⒉粌H僅是簡單的一個重力，還包含其軌跡控制中的航跡變化問題，需要補償?shù)氖且粋€未知的常值擾動。本文分析了超空泡航行體的常值擾動，采用擾動觀測器的方法對這個常值擾動進行估計和補償[10]。文中推導(dǎo)了擾動觀測器算法，設(shè)計了基于擾動觀測和補償?shù)某张莺叫畜w反饋控制方案， 避免了超空泡航行體的周期性滑水現(xiàn)象，無需加大執(zhí)行機構(gòu)的行程，且具有一定的魯棒性。

1 超空泡航行體動力學(xué)模型

本文采用文獻[7]中給出的超空泡航行體的基準(zhǔn)模型和參數(shù)。超空泡航行體的示意圖如圖1所示，其整體為長度比為1：2的梭形，總長度為L0，半徑為R，與水的密度比為m。超空泡航行體的頭部與空化器相連，空化器除了產(chǎn)生和維持空泡，還能偏轉(zhuǎn)一定的角度而產(chǎn)生流體動力，空化器和尾翼共同控制航行體的運動穩(wěn)定。

3 仿真驗證

為了說明反饋控制器（20）在加上擾動觀測器后的補償效果，下面分3種情況進行驗證說明。第1種情況是在航行體初始下沉速度w（0）=1，這種情況是文獻[7]和文獻[9]為了避免產(chǎn)生滑行而給出的典型算例，最終都以δc超出實際可能值而放棄了這個方案，而在本文設(shè)計的擾動補償下，δc完全可以在正常工作范圍內(nèi)消除滑行。第2種情況是航行體在反饋控制作用下典型的拍打著滑水前進的情況，在接上擾動觀測器并進行補償后，航行體很快就平穩(wěn)了下來。第3種情況則是航行體按照航跡要求航行的例子，可以看到，估計器隨航跡要求的xeq給出不同的擾動估計值。這3種情況都表明補償后航行體能平穩(wěn)運行，不再產(chǎn)生滑水的拍打現(xiàn)象。

第1種情況：初始速度w（0）=1，航跡輸入xeq=0。各個狀態(tài)變量的響應(yīng)曲線如圖5所示， 對應(yīng)的控制輸入及滑行力如圖6所示，觀測器輸出如圖7所示。從圖中可見，航行體僅與空泡壁碰撞1次，經(jīng)過1個波的振蕩后就逐漸收斂到平衡點，相應(yīng)的控制偏轉(zhuǎn)角均小于0.2 rad。約0.4 s后觀測器輸出趨于穩(wěn)定值，航行體恢復(fù)穩(wěn)定運行狀態(tài)。

第2種情況：初始速度w（0）=0，航跡輸入xeq=0，并遲后一秒接入觀測器補償。各個狀態(tài)變量的響應(yīng)曲線如圖8所示，對應(yīng)的控制輸入及滑行力如圖9所示。從圖中可見，在未接入擾動觀測器前，w在1 m/s和2 m/s之間波動，即在無滑行力和有滑行力之間波動（參見圖2），而俯仰角θ的平均值為0.02 rad，即航行體是微微的抬著頭，尾部周期性的拍打著空泡下壁高速前進，這就是當(dāng)前的航行體在重力影響下典型工作模式。接入擾動補償后，航行體很快平穩(wěn)下來，體現(xiàn)了補償?shù)男Ч?。在此過程中控制面的偏轉(zhuǎn)角均小于0.2 rad，滿足正常工作的范圍。

第3種情況：初始速度w（0）=0，軌跡輸入xeq變化曲線如圖10所示，即航行體以0.1 m/s的速度爬坡 （注：w、z向下為正方向），爬升至1 m后，保持5 s后，再以0.1 m/s的速度下降回到原來的深度。各個狀態(tài)變量的響應(yīng)曲線如圖11所示， 從圖中可見，除初始的過渡過程外，航行體按照預(yù)定的深度航跡平穩(wěn)運行，不與空泡壁碰撞。對應(yīng)的控制輸入及滑行力如圖12所示，圖中顯示，在此過程中控制面的偏轉(zhuǎn)角均小于0.2 rad，滿足實際工作要求。觀測器輸出如圖13所示，從圖中可見，x^s曲線在10 s和15 s上均有一個變化，這是由于航跡變化引起的。圖11和圖12表明，在整個航跡變化過程中已沒有拍打現(xiàn)象，體現(xiàn)了補償?shù)男Ч?/p>

4 結(jié) 論

超空泡航行體典型的工作模式是周期性拍打空泡壁高速前行的。這種周期性滑水現(xiàn)象產(chǎn)生的根源是重力的作用，如果補償了重力，即可以避免航行體下沉與空泡壁碰撞產(chǎn)生滑行力。對于超空泡航行體來說，常值擾動還包含航跡控制中的平衡點變化問題，所以這個擾動是一個未知的變化著的擾動。采用擾動觀測估計來補償系統(tǒng)的常值擾動，可以避免滑行力的產(chǎn)生，從而消除周期性滑水的運行模式。本文所設(shè)計的控制器能夠保證航行體與空泡壁不會連續(xù)碰撞，在空泡內(nèi)平穩(wěn)運行，且能夠保證控制面的偏轉(zhuǎn)角在正常工作范圍內(nèi)。

參 考 文 獻：

[1] VASIN A D. The principle of independence of the cavity sections expansion （Logvinovich's principle） as the Basis for investigation on cavitation flows[C]// National Tax，AssociationTax，America. 2001：161-162.

[2] SANABRIA D E， BALAS G J， ARNDT R E A.Modeling， control， and experimental validation of a highspeed supercavitating vehicle[J]. IEEE Journal of Oceanic Engineering， 2015，40（2）：362.

[3] SANABRIA D E， BALAS G J， ARNDT R E A. Planing avoidance control for supercavitating vehicles[C]// American Control Conference， 2014： 4979-498zz.

[4] MAO Xiaofeng， WANG Q. Delaydependent control design for a timedelay supercavitating vehicle model[J]. Journal of Vibration & Control， 2011，17（3）：431.

[5] SAVCHENKO Y N. Control of supercavitation flow and stability of supercavitating motion of bodies[C]// Vki Lecture Series Supercavitating Flows， 2001：11.

[6] KIRSCHNER I N， KRING D C， STOKES A W， et al. Control strategies for supercavitating vehicles[J]. Journal of Vibration and Control， 2002， 8：219.

[7] DZIELSKI J， KURDILA A. A benchmark control problem for supercavitating vehicles and an initial investigation of solutions[J]. Journal of Vibration and Control， 2003：791.

[8] DZIELSKI J. Longitudinal stability of a supercavitating vehicle[J]. IEEE Journal of Oceanic Engineering， 2011：562.

[9] GUO Jianlin， BANACHANDRAN B， ABED E H. Dynamics and control of supercavitating vehicles[J]. Journal of Dynamic Systems， Measurement， and Control， 2008， 130（1）：1.

[10] 王廣雄，何朕. 控制系統(tǒng)設(shè)計[M]. 北京：清華大學(xué)出版社，2005.

[11] LOGVINOVICH G V. Hydrodynamics of flows with free boundaries[M]. New York： Halsted Press， 1972.

（編輯：劉琳琳）