秋路,高正紅,劉艷
(西北工業(yè)大學(xué)航空學(xué)院,陜西西安710072)
斜置翼布局飛機(jī)可通過任意改變機(jī)翼的斜置角來滿足不同飛行狀態(tài)的要求,例如,進(jìn)行亞聲速巡航時(shí)選擇斜置角為零的大展弦比構(gòu)型,而進(jìn)行超聲速飛行時(shí)則選擇非零斜置角的小展弦比斜掠構(gòu)型。較之于變后掠翼布局,斜置翼布局能更有效地降低超聲速波阻,且沒有前者焦點(diǎn)后移、機(jī)翼轉(zhuǎn)軸結(jié)構(gòu)復(fù)雜等缺點(diǎn)[1]。然而,機(jī)翼斜置后非對(duì)稱的氣動(dòng)外形導(dǎo)致該布局飛機(jī)具有強(qiáng)烈的氣動(dòng)耦合及慣性耦合特性,因此,在研究其飛行動(dòng)力學(xué)特性或進(jìn)行控制律設(shè)計(jì)時(shí),不能采用常規(guī)的縱橫分離小擾動(dòng)模型,必須針對(duì)斜置翼布局飛機(jī)建立更為精準(zhǔn)的飛行動(dòng)力學(xué)模型。
目前,國(guó)外已從理論分析、氣動(dòng)布局設(shè)計(jì)、風(fēng)洞試驗(yàn)、驗(yàn)證機(jī)試制及試飛等方面對(duì)斜置翼技術(shù)的可行性進(jìn)行了大量的研究,但相關(guān)文獻(xiàn)并沒有完整地給出此類布局飛機(jī)動(dòng)力學(xué)模型的建模方法,尤其是線化模型,部分論文雖提供了數(shù)值形式的線化模型,但多是從非線性飛行模擬中提取而來的,沒有理論研究?jī)r(jià)值及指導(dǎo)意義;國(guó)內(nèi)文獻(xiàn)僅側(cè)重于其氣動(dòng)特性分析。因此,本文將在詳細(xì)分析斜置翼布局飛機(jī)耦合特性的基礎(chǔ)上,對(duì)其進(jìn)行六自由度建模,并基于小擾動(dòng)理論推導(dǎo)出縱橫耦合的線化動(dòng)力學(xué)模型。
斜置翼布局的特點(diǎn)在于可任意改變機(jī)翼的斜置角,可設(shè)計(jì)為翼身組合體形式或飛翼布局。圖1為斜置翼布局飛機(jī)F-8 OWRA(Oblique Wing Research Aircraft)的三視圖[2],此類布局通過單根轉(zhuǎn)軸實(shí)現(xiàn)斜掠,而飛翼布局則通過引擎的旋轉(zhuǎn)或控制面的操縱來實(shí)現(xiàn)。
機(jī)翼斜置角非零時(shí),斜置翼布局飛機(jī)具有不同于常規(guī)布局的耦合特性,包括氣動(dòng)耦合與慣性耦合。
氣動(dòng)耦合體現(xiàn)為縱向運(yùn)動(dòng)參數(shù)的變化會(huì)引起橫航向氣動(dòng)力與力矩的變化,橫航向運(yùn)動(dòng)參數(shù)的變化亦會(huì)引起縱向氣動(dòng)力與力矩的變化。
圖1 F-8 OWRA三視圖Fig.1 Three-views of F-8 OWRA
機(jī)翼斜置時(shí),從結(jié)構(gòu)上分為前掠的前翼和后掠的后翼,如圖2所示。首先,展向流動(dòng)使機(jī)翼后緣下洗渦沿展向呈非對(duì)稱變化,導(dǎo)致前后翼的附加升力分布不同,造成后翼加載,從而產(chǎn)生滾轉(zhuǎn)力矩,同時(shí)非對(duì)稱的誘導(dǎo)阻力產(chǎn)生偏航力矩[3];其次,超聲速飛行時(shí)展向流動(dòng)導(dǎo)致前后翼壓力分布不同,前翼處于壓縮流場(chǎng),后翼處于加速流場(chǎng),從而產(chǎn)生非對(duì)稱氣動(dòng)力[3];此外,偏轉(zhuǎn)副翼進(jìn)行滾轉(zhuǎn)操縱時(shí),前后副翼上的氣動(dòng)力會(huì)產(chǎn)生俯仰力矩。
圖2 斜置機(jī)翼示意圖Fig.2 Schematic diagram of an oblique wing
以斜置翼布局無人機(jī)Oblique Wing RPV(Remotely Piloted Vehicle)[4]為例,表 1 對(duì)斜置角為 0°,45°時(shí)的氣動(dòng)導(dǎo)數(shù)[4]作了對(duì)比。
表1 氣動(dòng)導(dǎo)數(shù)參考值Table 1 Reference values of aerodynamic derivatives
續(xù)表1
機(jī)翼斜置后,飛機(jī)無縱向?qū)ΨQ面,外形、質(zhì)量分布不對(duì)稱,故慣性積Ixy,Iyz均不為零,即出現(xiàn)縱向-橫航向的慣性耦合。以F-8 OWRA為例,表2對(duì)機(jī)翼斜置角為 0°,45°時(shí)的慣性矩及慣性積[2]進(jìn)行了對(duì)比。
表2 慣性量參考值Table 2 Reference values of inertial characteristics
特殊的耦合特性使斜置翼布局飛機(jī)不能直接采用常規(guī)的縱橫分離小擾動(dòng)模型,下面將建立適用于該布局的飛行動(dòng)力學(xué)模型。
在一般簡(jiǎn)化假設(shè)[5]的基礎(chǔ)上,對(duì)于斜置翼飛機(jī),還需假設(shè)機(jī)身外形及質(zhì)量對(duì)稱,而機(jī)翼外形及質(zhì)量不對(duì)稱。
氣流坐標(biāo)系Oxayaza定義為:原點(diǎn)O在飛機(jī)質(zhì)心處,坐標(biāo)系與飛機(jī)固連;Oxa軸與飛行速度矢量重合一致;Oza軸在機(jī)身對(duì)稱面內(nèi),并垂直于Oxa軸,指向機(jī)腹下方;Oya軸垂直于Oxaza平面,指向機(jī)身右方。
機(jī)體坐標(biāo)系Oxbybzb定義為:原點(diǎn)O在飛機(jī)質(zhì)心處;Oxb軸沿機(jī)身縱軸方向,指向前;Oyb軸垂直于機(jī)身對(duì)稱面,指向右;Ozb軸在機(jī)身對(duì)稱面內(nèi),垂直于Oxb軸,指向下。
若為飛翼布局,上述定義中機(jī)身對(duì)稱面改為機(jī)翼中心位置剖面。示意簡(jiǎn)圖如圖3所示(未畫出機(jī)身及Oz軸)。
圖3 斜置翼布局坐標(biāo)系定義Fig.3 Coordinate system definition of oblique wing aircrafts
定義在 Oxayaza的力方程見式(1),定義在Oxbybzb的力矩方程見式(2),運(yùn)動(dòng)學(xué)方程見式(3)。
根據(jù)小擾動(dòng)原理以及基于小擾動(dòng)原理的線性化方法[6],取基準(zhǔn)運(yùn)動(dòng)為零滾轉(zhuǎn)角的等速直線側(cè)滑飛行,非零參數(shù)為 V0,α0,β0,θ0,δeL0,δeR0及 δr0,對(duì)式(1)~式(3)進(jìn)行處理,可得到形如下式的斜置翼布局飛機(jī)線化飛行動(dòng)力學(xué)模型:
其中:
式中:γ0為基準(zhǔn)狀態(tài)的航跡傾斜角,γ0= θ0- α0。
選取與上文中解析表達(dá)式對(duì)應(yīng)的F-8 OWRA線化數(shù)值模型[2]進(jìn)行算例分析。斜置翼飛機(jī)的氣動(dòng)耦合與慣性耦合最終表現(xiàn)為運(yùn)動(dòng)耦合,可通過飛機(jī)對(duì)不同舵面的操縱響應(yīng)來分析。取計(jì)算飛行狀態(tài)Λ=45°,Ma=0.8,H=6 096 m,相應(yīng)的基準(zhǔn)運(yùn)動(dòng)參數(shù)為 V0=253 m/s,α0=1.6°,β0= -0.23°,δeL0=-6.3°,δeR0=2.9°,δr0= -3.4°。分別取駕駛員操縱輸入為升降舵對(duì)稱偏轉(zhuǎn)1°、副翼反對(duì)稱偏轉(zhuǎn)1°、方向舵偏轉(zhuǎn)1°的階躍信號(hào),得到飛機(jī)各角度、角速度的時(shí)間歷程曲線分別如圖4~圖6所示。
圖4 升降舵對(duì)稱偏轉(zhuǎn)1°時(shí)的響應(yīng)曲線Fig.4 Response to 1°symmetric elevator deflection
圖5 副翼反對(duì)稱偏轉(zhuǎn)1°時(shí)的響應(yīng)曲線Fig.5 Response to 1°asymmetric aileron deflection
圖6 方向舵偏轉(zhuǎn)1°時(shí)的響應(yīng)曲線Fig.6 Response to 1°rudder deflection
在以上三種操縱反應(yīng)中,p及φ變化均很大;此外,升降舵對(duì)稱偏轉(zhuǎn)1°后引起p及φ的變化量大于副翼偏轉(zhuǎn)1°后相應(yīng)的變化量,其原因一是因?yàn)榭v橫耦合現(xiàn)象嚴(yán)重,二是因?yàn)闄C(jī)翼斜置后副翼的操縱效率降低;相對(duì)于升降舵偏轉(zhuǎn)引起的橫航向運(yùn)動(dòng)參數(shù)變化量,副翼及方向舵偏轉(zhuǎn)時(shí)引起的縱向運(yùn)動(dòng)參數(shù)變化量較小,但方向舵偏轉(zhuǎn)1°時(shí)的橫向耦合較為嚴(yán)重。
考慮到斜置翼布局飛機(jī)非對(duì)稱的氣動(dòng)特性,選擇駕駛員操縱方式分別為升降舵同向差動(dòng)偏轉(zhuǎn)進(jìn)行縱向操縱、副翼反向差動(dòng)偏轉(zhuǎn)進(jìn)行橫向操縱,可以改善縱橫耦合現(xiàn)象;選擇升降舵反向差動(dòng)偏轉(zhuǎn)進(jìn)行橫向操縱可以在機(jī)翼斜置角較大時(shí)提高橫向操縱效率。
本文建立的斜置翼布局飛機(jī)線化小擾動(dòng)模型可用來研究其縱橫耦合的運(yùn)動(dòng)特性,但僅限于運(yùn)動(dòng)參數(shù)變化量較小的情況,此時(shí)主要?dú)鈩?dòng)參數(shù)的變化與擾動(dòng)量呈線性關(guān)系;更精確的動(dòng)力學(xué)模型還需考慮非線性因素及氣動(dòng)彈性的影響;此外,該模型可用于飛行控制律的設(shè)計(jì)。
采用本文所建立的線化模型進(jìn)行相關(guān)計(jì)算及分析前應(yīng)通過試飛數(shù)據(jù)或模擬飛行進(jìn)行數(shù)值模型驗(yàn)證,由于條件不足,本文未進(jìn)行驗(yàn)證。
為了使斜置翼布局飛機(jī)的飛行品質(zhì)滿足要求,可通過以下兩種途徑改善其縱橫耦合特性:(1)改進(jìn)機(jī)翼的氣動(dòng)外形使其產(chǎn)生對(duì)稱的升力分布,如彎扭設(shè)計(jì)、翼尖修形、氣動(dòng)彈性設(shè)計(jì)等;(2)設(shè)計(jì)解耦控制律。
[1] Michael J Hirschberg,David M Hart.A summary of a half-century of oblique wing research[C]//AIAA Proceedings of the 45th Aerospace Science Meeting and Exhibit.Reno Nevada:AIAA,2007:1-35.
[2] Robert W Kempel,Walter E McNeil.A piloted evaluation of an oblique-wing research aircraft motion simulation with decoupling control laws[R].NASA TP-2874,1988.
[3] Desktop Aeronautics Inc.Oblique flying wing:an introduction and white paper[R].USA:Desktop Aeronautics Inc,2005.
[4] Travassos R H,Iliff K W.Determination of an oblique wing aircraft’s aerodynamic characteristics[R].USA:AIAA,1980.
[5] 胡兆豐,何植岱,高浩.飛行動(dòng)力學(xué)[M].第1版.西安:西北工業(yè)大學(xué)出版社,2007:21.
[6] 吳森堂,費(fèi)玉華.飛行控制系統(tǒng)[M].第1版.北京:北京航空航天大學(xué)出版社,2005:44,64-66.