楊雨，陸宇平，戴正升

(南京航空航天大學(xué)自動化學(xué)院，江蘇南京210016)

0 引言

大型運(yùn)輸機(jī)一般擁有非常大的翼展，且機(jī)翼有一定的下反角，使得飛機(jī)執(zhí)行超低空空投任務(wù)時(shí)機(jī)翼相對于飛機(jī)質(zhì)心距離地面更近。由于太陽照射等因素，低空氣流相對于高空更加復(fù)雜，飛機(jī)易受到強(qiáng)側(cè)風(fēng)影響［1-3］，一旦控制律不能對飛機(jī)的軌跡實(shí)施良好的控制，會使飛機(jī)偏離空投區(qū)域，錯(cuò)過空投機(jī)會，且姿態(tài)不穩(wěn)定可能導(dǎo)致翼尖直接觸地，致使飛機(jī)墜毀。理想飛行情況下，飛機(jī)執(zhí)行超低空空投任務(wù)時(shí)橫側(cè)向應(yīng)滿足:(1)飛機(jī)地速與空投區(qū)域中心線方向一致;(2)飛機(jī)相對目標(biāo)空投區(qū)域的側(cè)向偏移為0;(3)機(jī)體軸與目標(biāo)空投區(qū)域中心線方向一致;(4)飛機(jī)的滾轉(zhuǎn)角為0°，以避免大翼展飛機(jī)翼尖觸地。但多數(shù)情況下，飛機(jī)幾乎不能在理想狀況下飛行，因而在橫側(cè)向上對飛機(jī)姿態(tài)的要求也會改變，側(cè)向偏離和航跡偏角應(yīng)根據(jù)目標(biāo)空投區(qū)域具體情況適當(dāng)修改，在飛機(jī)進(jìn)入空投段時(shí)可適當(dāng)放寬偏航角來進(jìn)行抗側(cè)風(fēng)控制，此外，考慮到空投區(qū)域地面可能會有起伏，且飛機(jī)高度并不穩(wěn)定，所以需在橫側(cè)向留有一定余度，這將進(jìn)一步減小飛機(jī)允許的最大滾轉(zhuǎn)角。

飛機(jī)在正常的低空飛行和著陸的抗側(cè)風(fēng)控制律設(shè)計(jì)上已比較成熟，主要根據(jù)飛機(jī)的橫側(cè)向控制目的設(shè)計(jì)基于側(cè)滑法和側(cè)航法的橫側(cè)向控制系統(tǒng)，同時(shí)，針對可能提供直接側(cè)力的控制面設(shè)計(jì)的直接側(cè)力法，都取得了較好的控制效果［4-5］。

本文以超低空空投橫側(cè)向非線性模型及紊流模型為基礎(chǔ)設(shè)計(jì)了PID控制器和橫側(cè)向抗側(cè)風(fēng)控制器，仿真結(jié)果表明，具有更理想的控制效果。

1 數(shù)學(xué)模型

分析飛機(jī)的動力學(xué)與運(yùn)動學(xué)方程［6］，結(jié)合飛機(jī)空投時(shí)的受力情況，可建立飛機(jī)橫側(cè)向數(shù)學(xué)模型如下:

當(dāng)貨物移動時(shí)，式(1)中:

式中，m1為飛機(jī)自身質(zhì)量;m2為貨物質(zhì)量;f2為貨物對飛機(jī)的作用力在機(jī)體坐標(biāo)系y軸的分量;r1，r3分別為貨物質(zhì)心相對飛機(jī)質(zhì)心的位移在機(jī)體坐標(biāo)系x軸和z軸的分量;Tpy為牽引傘拉力在機(jī)體坐標(biāo)系y軸的分量。其它參數(shù)含義見文獻(xiàn)［6］。

2 基于側(cè)航法的PID控制器設(shè)計(jì)及仿真

2.1 控制器設(shè)計(jì)

橫側(cè)向飛行模態(tài)控制律設(shè)計(jì)均以飛機(jī)增穩(wěn)系統(tǒng)為內(nèi)回路，內(nèi)回路結(jié)合側(cè)滑法控制律，在副翼通道采用側(cè)滑角反饋，方向舵通道采用側(cè)滑角和偏航角速率的綜合反饋，目的是改善橫側(cè)向飛行品質(zhì)，以滿足飛行品質(zhì)中關(guān)于荷蘭滾阻尼和頻率、滾轉(zhuǎn)模態(tài)時(shí)間常數(shù)和螺旋模態(tài)時(shí)間常數(shù)的要求，控制結(jié)構(gòu)如圖1所示。

圖1 內(nèi)回路控制結(jié)構(gòu)圖

圖中，Gδ(s)為舵機(jī)傳遞函數(shù);GT(s)為清洗網(wǎng)絡(luò)。此時(shí)控制律可以用下式表示:

以圖1為內(nèi)回路，利用滾轉(zhuǎn)角反饋的PID控制，實(shí)現(xiàn)滾轉(zhuǎn)角保持，控制結(jié)構(gòu)如圖2所示。

圖2 滾轉(zhuǎn)角保持控制結(jié)構(gòu)圖

最后側(cè)向偏離的控制通過把滾轉(zhuǎn)角保持回路作為內(nèi)回路，通過PD控制來實(shí)現(xiàn)，為了保持空投過程的滾轉(zhuǎn)角，在滾轉(zhuǎn)角命令上加入飽和環(huán)節(jié)進(jìn)行限制，其控制律結(jié)構(gòu)如圖3所示。

圖3 外回路控制結(jié)構(gòu)圖

在方向舵通道不做任何改變，仍保持增穩(wěn)回路。主要通過副翼通道控制飛機(jī)的側(cè)向偏離，同時(shí)對滾轉(zhuǎn)角指令進(jìn)行限制，也限制了飛機(jī)的滾轉(zhuǎn)角。

2.2 數(shù)值仿真

假設(shè)飛機(jī)重120 t，貨物重20 t，飛行速度Ma=0.23。根據(jù)飛行品質(zhì)設(shè)置橫側(cè)向增穩(wěn)系統(tǒng)參數(shù)，令Gδ(s)=1/(10s+1)，GT(s)=s/(s+1)，Kaβ=1.117，。滾轉(zhuǎn)角保持回路的PID控制參數(shù)設(shè)置為:Kpφ=15.7，Kiφ=0.35，Kdφ=16。在滾轉(zhuǎn)角保持回路基礎(chǔ)上設(shè)計(jì)側(cè)向偏離回路，PD控制器參數(shù)選取為:Kpy=1.1，Kdy=4.2。

在模型中加入大氣擾動來對橫側(cè)向控制系統(tǒng)抗側(cè)風(fēng)進(jìn)行仿真。令側(cè)向偏離指令為零，滾轉(zhuǎn)角限制在±5°。飛機(jī)側(cè)向偏離、滾轉(zhuǎn)角、偏航角、側(cè)滑角仿真曲線如圖4所示。

圖4 仿真曲線

從圖4中可以看出，飛機(jī)機(jī)頭偏轉(zhuǎn)了一個(gè)偏航角來抵抗側(cè)風(fēng)，滾轉(zhuǎn)角和側(cè)滑角都穩(wěn)定到0°，滾轉(zhuǎn)角最大偏轉(zhuǎn)已經(jīng)接近3°，側(cè)向偏離也穩(wěn)定到0 m。對比C-130對于超低空空投側(cè)風(fēng)的要求，6 m/s的側(cè)風(fēng)在超低空空投中屬于較強(qiáng)的側(cè)風(fēng)，所以滾轉(zhuǎn)角已經(jīng)偏大。

系統(tǒng)的副翼和方向舵的輸入曲線如圖5所示。

圖5 舵面的輸入曲線

由以上仿真結(jié)果可知，基于側(cè)航法設(shè)計(jì)的控制律使飛機(jī)在橫側(cè)向紊流和常值風(fēng)擾動下偏轉(zhuǎn)一定的偏航角來對抗側(cè)風(fēng)，飛機(jī)的側(cè)向位移和滾轉(zhuǎn)角都能回到初值0，使飛機(jī)對準(zhǔn)目標(biāo)空投區(qū)域并保持飛機(jī)水平。需說明的是，以上仿真中飛機(jī)帶貨物并進(jìn)行空投，所以完全可反映空投時(shí)飛機(jī)的真實(shí)飛行狀態(tài)。從仿真曲線可以看出，基于側(cè)航法對抗側(cè)風(fēng)時(shí)，副翼是主要的控制舵面，方向舵只是予以輔助，轉(zhuǎn)動幅度很小。

3 基于動態(tài)逆的抗側(cè)風(fēng)控制器設(shè)計(jì)及仿真

基于側(cè)航法對抗側(cè)風(fēng)時(shí)，主要利用了副翼通道，舵面偏轉(zhuǎn)很大，且滾轉(zhuǎn)角瞬時(shí)值偏大，下面設(shè)計(jì)一種基于最優(yōu)調(diào)節(jié)器的動態(tài)逆控制器，在利用副翼的基礎(chǔ)上，同時(shí)利用方向舵通道以減小滾轉(zhuǎn)角的瞬時(shí)值。

3.1 橫側(cè)向模型的精確線性化

由于非線性系統(tǒng)不能像線性系統(tǒng)那樣得到狀態(tài)的解析表達(dá)式，在一定條件下通過狀態(tài)和輸入變量的適當(dāng)變化實(shí)現(xiàn)所有或部分狀態(tài)的線性化，從而將復(fù)雜的非線性控制系統(tǒng)設(shè)計(jì)問題轉(zhuǎn)化為簡單的線性系統(tǒng)設(shè)計(jì)問題。

為實(shí)現(xiàn)模型精確線性化，需利用非線性系統(tǒng)分析中的重要概念和運(yùn)算符號:李導(dǎo)數(shù)和李括號，這將在后續(xù)的分析中大量應(yīng)用，詳細(xì)定義見文獻(xiàn)［7］。

考慮輸入輸出維數(shù)相等的系統(tǒng)如下式:

式中，x∈Rn，u∈Rm及y∈Rm分別為系統(tǒng)的狀態(tài)向量、輸入向量和輸出向量;f(x)為n維充分光滑的向量場;G(x)=(g1(x)，…，gm(x))，gi(i=1，…，m)為n維充分光滑的向量場;H(x)=(h1(x)，…，hm(x))，hi(i=1，…，m)為充分光滑的標(biāo)量函數(shù)。

對于非線性系統(tǒng)式(4)，系統(tǒng)中第j個(gè)輸出yj對時(shí)間的導(dǎo)數(shù)為:

在上式中，如果所有Lgihj≡0，則導(dǎo)數(shù)與控制無關(guān)。若γj使導(dǎo)數(shù) y(γj)j至少依賴一個(gè)輸入的最小正數(shù)，即:

若如上定義的相對階向量有定義，則式(5)可表示為:

由于矩陣A(x)可逆，可取狀態(tài)反饋控制律為:

可以得到輸入/輸出的動態(tài)方程為:

那么，通過輸入/輸出之間的精確線性化后所得的反饋表達(dá)式(9)就是原非線性系統(tǒng)式(4)的逆系統(tǒng)，如果將式(8)的輸出看成式(7)的輸入，可得到式(9)，從而實(shí)現(xiàn)對非線性系統(tǒng)式(4)的動態(tài)逆控制［8］。同時(shí)，通過逆控制作用，還實(shí)現(xiàn)了對非線性系統(tǒng)的輸入/輸出解耦，因此矩陣A(x)也稱為解耦矩陣。

從式(2)中可以看出，r·1，Tpy對于橫側(cè)向的影響有限，令r·1=Tpy=0。式(1)中，側(cè)力 C、滾轉(zhuǎn)力矩L、偏航力矩N可以通過相關(guān)的側(cè)力導(dǎo)數(shù)、滾轉(zhuǎn)力矩導(dǎo)數(shù)和偏航力矩導(dǎo)數(shù)計(jì)算。

把飛機(jī)橫側(cè)向模型化作式(4)的形式，可以得到:將各導(dǎo)數(shù)代入模型并展開可得到上式的各個(gè)量。

系統(tǒng)的輸入只有升降舵和副翼，考慮到超低空空投時(shí)橫側(cè)向的控制目標(biāo)，飛機(jī)滾轉(zhuǎn)角和側(cè)向偏離是決定空投能否成功的關(guān)鍵變量，所以把滾轉(zhuǎn)角和側(cè)向偏離作為輸出進(jìn)行輸入/輸出線性化。經(jīng)計(jì)算:Lg1Lfh2≠0，Lg2Lfh2≠0，那么 γ2=1。同理可得 γ6=2。

代入式(7)得:

取狀態(tài)反饋控制律為:

可以得到輸入/輸出的動態(tài)方程為:

3.2 動態(tài)逆控制器設(shè)計(jì)

對橫側(cè)向定義向量 x=［x1，x2］T，其中:

式中，φd為給定的滾轉(zhuǎn)角值，控制目標(biāo)為0。

對式(14)作微分運(yùn)算，可得到:

式中，u1=

那么，如果能夠得到u1，u2的表達(dá)式，就可以得到動態(tài)逆控制器，考慮到便于計(jì)算和方便應(yīng)用，使用最優(yōu)狀態(tài)調(diào)節(jié)器來求解，其表達(dá)式如下:

式中，矩陣P1，P2滿足以下黎卡提矩陣方程;R1，R2為選定的正定對稱矩陣。

那么控制器的輸入可以表達(dá)為:

最后得到橫側(cè)向抗側(cè)風(fēng)動態(tài)逆控制律為:

控制律結(jié)構(gòu)如圖6所示。

圖6 動態(tài)逆控制器結(jié)構(gòu)圖

需要說明的是，系統(tǒng)需要貨物的位置作為輸入，這里把貨物的位置設(shè)定為初始位置，如果貨物的位置和速度都可以測量，則可以直接代入方程，如果不能測量，那么可以代入貨物的初始位置，不做改變。

3.3 仿真結(jié)果

經(jīng)過對飛機(jī)風(fēng)洞數(shù)據(jù)分析和擬合，得到在h=5 m，Ma=0.23的狀態(tài)飛行時(shí)，飛機(jī)橫側(cè)向氣動導(dǎo)數(shù)如下:CYβ=－1.575 2，CYˉp=CYˉr=0，CYδr=－0.008，Clβ=－0.201，Clˉp=－0.182 3，Clˉr=－0.073 6，Clδr=－0.049，Clδa=0.001 1，Cnβ=－0.201，Cnˉp=0，Cnˉr= － 0.003 2，Cnδr= － 0.090 2，Cnδa=0。

飛機(jī)的轉(zhuǎn)動慣量為:Ix=7 300 000 kg·m·s2，Iz=15 800 000 kg·m·s2，α0=θ0=1.022 6°。給定的滾轉(zhuǎn)角和側(cè)向偏離均為零。最優(yōu)調(diào)節(jié)器的參數(shù)選取如下:

在大氣擾動有紊流和6 m/s側(cè)風(fēng)的環(huán)境下進(jìn)行仿真。仿真曲線如圖7所示。

從圖中可以看出，在相同的側(cè)風(fēng)干擾下，滾轉(zhuǎn)角峰值幾乎是側(cè)航法的一半，原因在于系統(tǒng)的輸入上，方向舵在有側(cè)風(fēng)干擾時(shí)也做出較大偏轉(zhuǎn)來減小側(cè)風(fēng)的影響，同時(shí)系統(tǒng)能夠很快達(dá)到穩(wěn)定，而且對于紊流的抑制也表現(xiàn)出非常好的穩(wěn)定性，有紊流干擾時(shí)，并不需要副翼來回偏轉(zhuǎn)用以對抗紊流帶來的影響。

圖7 仿真曲線

4 結(jié)束語

本文設(shè)計(jì)的兩種控制器均能消除橫側(cè)向紊流及常值風(fēng)影響，避免了因滾轉(zhuǎn)角過大使飛機(jī)翼尖觸地，確保貨物順利出艙，并使貨物空投側(cè)向偏離較小。而基于最優(yōu)調(diào)節(jié)器的動態(tài)逆控制器比基于側(cè)航法的PID控制器更能減小滾轉(zhuǎn)角的幅值，且側(cè)向偏離和偏航角能很快達(dá)到穩(wěn)定，同時(shí)副翼不需要來回較大偏轉(zhuǎn)，具有更好的控制效果。

［1］呂新波，張欣，劉振欽．農(nóng)林飛機(jī)超低空飛行大氣紊流影響研究［J］．飛行力學(xué)，2008，26(2):64-68．

［2］張婧，陳瀾，李曉曦，等．針對大氣紊流改進(jìn)的飛控系統(tǒng)設(shè)計(jì)及仿真研究［J］．計(jì)算機(jī)測量與控制，2011，19(4):860-862．

［3］Blick M E．Effectof wind turbulence and shear on landing performance of jet transports［R］．AIAA-78-20771，1978．

［4］王鵬，陳懷民，吳成富，等．飛翼無人機(jī)著陸過程中的抗側(cè)風(fēng)控制研究［J］．計(jì)算機(jī)仿真，2009，26(12):78-81．

［5］嵇鼎毅，陸宇平．飛翼布局無人機(jī)抗側(cè)風(fēng)自動著陸控制［J］．飛機(jī)設(shè)計(jì)，2007，27(2):25-33．

［6］吳森堂，費(fèi)玉華．飛行控制系統(tǒng)［M］．第1版．北京:北京航空航天大學(xué)出版社，2005:9-15．

［7］Lu Qiang，Sun Yuanzhang，Mei Shengwei．Nonlinear control systems and power system dynamics［M］．Netherlands:Kluwer Academic Publishers，2001:38-45．

［8］Binh Dang Vu．Nonlinear dynamic inversion control［M］．Berlin Heidelberg:Springer，1997:102-111．