于志軍，陶 楊，孫吉良，孫偉奇

（1.91206 部隊，山東 青島 266108；2.海軍航空工程學(xué)院 研究生管理大隊，山東 煙臺 264001）

在地面附近，由于空氣粘性的作用，風(fēng)速往往是不均勻的，會形成圖1所示的風(fēng)剖面，這種現(xiàn)象稱為風(fēng)切變。尤其是低空風(fēng)切變，會嚴重威脅到飛機在空中的飛行姿態(tài)，能使飛機連續(xù)不斷地顛簸，甚至造成機毀人亡的后果[1-3]。所以，加強對風(fēng)切變影響下飛機運動特征的研究具有重要意義。如未作特殊說明，本文接下來的研究內(nèi)容均為鉛垂平面內(nèi)飛機的運動。

圖1 風(fēng)切變剖面圖

1 飛機縱向運動方程

1.1 未考慮風(fēng)切變時的運動方程

在未考慮風(fēng)切變對飛機縱向運動的影響時，風(fēng)速w=0，飛機的地速v 與空速u相等。此時，在航跡—機體坐標(biāo)系下的運動方程為：

式中：m為飛機質(zhì)量；P為發(fā)動機推力；α為幾何迎角；Q為飛機所受的空氣阻力；θ為航跡傾斜角；Y為飛機所受的升力；zω為繞機體z軸（z軸垂直于飛機對稱平面，指向右為正）的轉(zhuǎn)動角速度；Mz為繞機體z軸的力矩；eP為發(fā)動機推力線與飛機質(zhì)心之間的偏心距離。

此外，還存在如下關(guān)系：

式中，?為俯仰角。

1.2 考慮風(fēng)切變時的運動方程[4]

在存在風(fēng)切變的情況下，飛機的地速v，空速u與風(fēng)速w的關(guān)系式為：

把式(4)投影到氣流坐標(biāo)系中，有

式(6)展開后得到

顯然，式(1)的第三個方程可以拋棄。

式(8)與式(1)的第三個方程共同組成風(fēng)切變下的縱向運動方程。

2 縱向運動方程的小擾動線性化方程

2.1 方程線性化的一般方法

飛行器的任何一個運動都可以寫成如下形式：

式中：x,y,…是運動參數(shù)或者它們的時間導(dǎo)數(shù)。那么運動參數(shù)可以表示成基準(zhǔn)量 a*與偏量 ?a 之和x=x*+?x,y=y*+?y,…，則式(9)可寫成：

在基準(zhǔn)點(x*,y*,…)處展開成Taylor級數(shù)，并只保留一階小量[5-7]，則式(10)可表示為：

因為基準(zhǔn)運動也應(yīng)滿足式(9)，故式(11)簡化為

該結(jié)果即為運動方程線性化的一般形式。

在運動方程線性化時，將出現(xiàn)力和力矩的偏量。一般來說，若某個力或力矩是若干變量的函數(shù)A=A(a,b,…)，則A的偏量應(yīng)為：

式中，Aa,Ab,…表示在基準(zhǔn)點的偏導(dǎo)數(shù)，即

2.2 未考慮風(fēng)切變時縱向運動方程的線性化

若不考慮飛行高度的變化對力和力矩的影響，則按照2.1節(jié)所述方法，對式(1)進行線性化，得

同時，對式(2)、(3)分別進行線性化，得到：

將式(15)、(16)代入式(14)，消去變量 ?θ 和 ?ωz；把?v、? ?、?ωz、?α看做狀態(tài)變量，?Pc看做輸入變量，得到如下的耦合方程組：

2.3 考慮風(fēng)切變時縱向運動方程的線性化

依據(jù)2.2節(jié)，同理可得存在風(fēng)切變下的線性狀態(tài)方程組，如下所示

3 平飛情況下縱向模態(tài)算例

第2節(jié)已給出了在不考慮飛行高度變化時有無風(fēng)切變飛機的縱向運動?，F(xiàn)給出在基準(zhǔn)運動為水平飛行，即 θ*=0時飛機縱向模態(tài)算例。以一定的初始條件[8]分別代入式(18)、(19)，求出對應(yīng)的特征值，結(jié)果如表1、表2所示。

表1 未考慮風(fēng)切變的縱向模態(tài)參數(shù)

表2 考慮風(fēng)切變的縱向模態(tài)參數(shù)

由上述結(jié)果可以看出，在受到風(fēng)切變擾動后，飛機長周期和短周期模態(tài)的半衰期T12、振蕩頻率ω 均有所縮短，而周期T 延長。

4 結(jié)論

第3節(jié)出現(xiàn)的結(jié)果，主要因為飛機在受到風(fēng)切變影響時，相當(dāng)于間接減小了升阻比，導(dǎo)致了飛機半衰期的縮短和周期的延長，同時在長周期模態(tài)下，飛機的飛行品質(zhì)也有一定幅度的下降，而短周期模態(tài)下則下降不大。在近似處理中，由于略去了轉(zhuǎn)動能，且阻力與推力始終平衡時，長周期模態(tài)體現(xiàn)了飛機在保守力場中的等幅振蕩運動，位能增加時動能減少，位能減少時動能增加?？紤]風(fēng)切變時，高度增加，氣流速度也增加，則對飛機來說，就從大氣流動中獲得了一部分使動能增加的能量，這部分能力增加的值越大，飛機振蕩頻率就越低。最后，當(dāng)增加的能力超過某一值后，就會出現(xiàn)非周期發(fā)散現(xiàn)象。

[1]葉汝杰.低空風(fēng)切變對飛行試驗的影響[J].電子對抗試驗,2001,34(1):68-70.

[2]徐可,唐長紅,李洲圣.風(fēng)切變與飛機響應(yīng)研究[J].理論與應(yīng)用力學(xué)學(xué)報,1997,24(3):35-39.

[3]劉志華,顧志濤.低空風(fēng)切變對飛行的危害及對策研究[J].飛行試驗,1998,14(1):35-37.

[4]肖業(yè)倫.飛行器運動方程[M].北京:航空工業(yè)出版社,1987:95-99.

[5]JOE IMAE,SATORU NAKATANI,JUNYA TAKAHASH.GP based flight control in the windshear[C]//IEEE SNNC’99 Conference Proceedings.1999:650-653.

[6]GEORGE LEITMANN,SANDEEP PANDEY.Adaptive control of aircraft in windeshear[C]//Proceedings of 30th Conference on Decision and Control.1991:2157-2161.

[7]陳揚鑒,李翰芝,呂玉虎.飛機在風(fēng)切變下進場的模擬試驗研究[J].飛行力學(xué),1995,13(1):75-83.

[8]熊海泉,劉昶,鄭本武.飛機飛行動力學(xué)[M].北京:航空工業(yè)出版社,1990:182.