變后掠變展長飛行器動力學(xué)建模與動態(tài)響應(yīng)分析

楊貫通,唐勝景,趙林東,郭杰

(北京理工大學(xué)飛行器動力學(xué)與控制教育部重點實驗室,北京 100081）

變后掠變展長飛行器動力學(xué)建模與動態(tài)響應(yīng)分析

楊貫通,唐勝景,趙林東,郭杰

(北京理工大學(xué)飛行器動力學(xué)與控制教育部重點實驗室,北京 100081）

飛行器機翼后掠角和展長動態(tài)變化時,氣動力、轉(zhuǎn)動慣量、壓力中心和質(zhì)心等關(guān)鍵參數(shù)都發(fā)生劇烈變化?；贜ewton-Euler方法建立了包含5個剛體的變后掠變展長飛行器的動力學(xué)模型,推導(dǎo)了由于變形所產(chǎn)生的附加氣動力和氣動力矩。在縱向上解耦簡化了動力學(xué)方程,基于準(zhǔn)定常氣動力模型,分析了不同變形形式和變形速度下的縱向動態(tài)響應(yīng)。研究結(jié)果表明,后掠角和展長的變化對動力學(xué)特性影響較大,但不同的變化速度對結(jié)果影響不大。通過合適的變形方式和變形速度的結(jié)合可以弱化變形過程中的動力學(xué)參數(shù)波動,減小控制系統(tǒng)的負擔(dān)。

應(yīng)用力學(xué);無人機;變后掠變展長;動力學(xué)建模;動態(tài)響應(yīng)

0 引言

傳統(tǒng)的飛行器最優(yōu)功能和性能往往只是體現(xiàn)在某些特定的飛行任務(wù)和飛行狀態(tài)上。變氣動外形飛行器能夠滿足大空域和大速域的飛行要求,可以在多種飛行狀態(tài)下保持最優(yōu)飛行,使得該飛行器能夠執(zhí)行常規(guī)飛行器所不能完成的任務(wù)[1-2]。William提煉出了在變形飛機概念設(shè)計階段最為關(guān)鍵的6個變量:推重比、翼載、翼厚度與弦長比、翼根梢比、翼后掠角、翼展弦比[3]。變后掠變展長無人飛行器綜合了翼后掠角和翼展弦比二種大尺度變形方式,是可變外形飛行器的一個重要研究方向。

對于變外形無人飛行器,氣動外形改變是一個關(guān)鍵的飛行過程,在變形過程中飛行器的質(zhì)量特性、氣動特性都發(fā)生顯著變化,需要對該過程進行詳細的建模和分析,把握關(guān)鍵參數(shù)的變化規(guī)律,為控制系統(tǒng)設(shè)計提供依據(jù)。在變后掠變展長飛行器方面,已有一些氣動特性研究成果[4],動力學(xué)研究還比較少見;在變外形飛行器動力學(xué)建模和動態(tài)響應(yīng)分析方面,已有對變后掠飛機以及折疊飛機等有意義的研究[5-8]。

1 變后掠變展長無人飛行器

在某原型無人機(如圖1所示)的基礎(chǔ)上,設(shè)計了后掠角和展長均可變的無人飛行器。特殊的氣動外形使其能夠?qū)崿F(xiàn)偵察和打擊一機融合,實現(xiàn)自主空中偵察、戰(zhàn)場監(jiān)視、戰(zhàn)場決策和打擊等任務(wù)。其3種典型的飛行狀態(tài)如圖2所示。

圖1 作為基準(zhǔn)模型的小型無人機Fig.1 The referential UAV

圖2 3種典型飛行狀態(tài)Fig.2 Typical flight states

2 多體動力學(xué)建模

變外形飛行器是一個空間運動的多體系統(tǒng),傳統(tǒng)的6自由度飛行器動力學(xué)模型不再適用,需要基于多體動力學(xué)建模方法重新推導(dǎo)動力學(xué)方程。

2.1 坐標(biāo)系的建立

變外形飛行器一般有3個特征點可以選作坐標(biāo)系原點:總質(zhì)心、機體質(zhì)心和飛機頭部。由于總質(zhì)心是變化的,坐標(biāo)系也會在飛行器上運動,增加了建模難度;飛機頭部作為原點,會直接增加力矩項,使建模過程變得復(fù)雜。而機體坐標(biāo)系原點選在機身質(zhì)心上,建模最為直觀和簡潔。在變形之前,即0°后掠角時,整個無人機的總質(zhì)心與機身的質(zhì)心重合。變形時,飛機質(zhì)心位置和坐標(biāo)系原點不再重合,并且隨后掠角和展長的變化在空間運動,將產(chǎn)生附加動力和力矩。為方便研究,將可變后掠翼和伸縮翼的質(zhì)量集中于其質(zhì)心。

如圖3所示,OXYZ為機體坐標(biāo)系,Oxyz為地面坐標(biāo)系。地面坐標(biāo)系到機體坐標(biāo)系的矢徑為ro,gi點為單個可活動機翼的質(zhì)心。機體坐標(biāo)系原點到gi點的矢徑為Si,地面坐標(biāo)系到gi點的矢徑為rgi.用下標(biāo)f表示機身,下標(biāo)1表示左后掠翼,2表示右后掠翼;3表示左伸縮翼,4表示右伸縮翼。對機身和機翼易得速度矢量(1)式。

圖3 坐標(biāo)系示意圖Fig.3 Coordinate system

2.2 平動動力學(xué)

對機身、4片機翼分別列寫動量定理

式中:Faf表示空氣動力對機身的作用力;F1f表示左側(cè)后掠翼對機身的作用力;G1代表左側(cè)后掠翼的重力;P代表發(fā)動機推力;其他類推。整理(2)式并根據(jù)(1)式得到

式中:m是總質(zhì)量;G是總重力?？梢园l(fā)現(xiàn)與傳統(tǒng)飛行器6自由度方程相比,該方程多了等式右側(cè)機翼運動附加力。進一步推導(dǎo)得到

式中:ω為飛行器的旋轉(zhuǎn)角速度,將ω連同S表示成機體坐標(biāo)系內(nèi)的向量形式,帶入(4)式可得到附加力。需要注意,該力是表示在機體坐標(biāo)系下的。

2.3 轉(zhuǎn)動動力學(xué)

對第i塊機翼,根據(jù)動量矩定義[9]知

式中:上標(biāo)*表示機翼上的單位微元。

直接對兩側(cè)求導(dǎo)得到

整理之后得

(10)式等式左側(cè)除Ma為飛行器受到的氣動力矩,其余均是附加力矩項。

3 縱向動力學(xué)模型

3.1 縱向動力學(xué)模型

兩側(cè)對稱變形時,由于左右兩側(cè)具有運動的一致性和對稱性,易得

此時附加力和附加力矩可以按后掠翼和伸縮翼分為2組。根據(jù)(11)式對力和力矩方程進行簡化,得到由后掠翼引起的機體坐標(biāo)系OX軸、OY軸上的附加力

由后掠翼引起的關(guān)于Z軸的附加力矩為

vX、vY為速度在機體坐標(biāo)系上的投影。最終得到變后掠變展長飛行器縱向動力學(xué)方程

3.2 幾何關(guān)系

飛行器的幾何關(guān)系如圖4所示。

由圖4易知以下幾何關(guān)系

式中:a、b是與外形有關(guān)的定值;后掠角χ和伸縮展長l均是控制量。(12)式中所需的對時間一階和二階導(dǎo)數(shù)均易導(dǎo)出。

圖4 幾何關(guān)系示意圖Fig.4 The geometric relationship

4 動態(tài)響應(yīng)仿真與分析

對變形過程進行動力學(xué)仿真,變形之前飛行器處于定直平飛狀態(tài),主要參數(shù)列于表1.變形過程中采用準(zhǔn)定常氣動力計算。

表1 平飛主要參數(shù)Tab.1 Parameters of straight flight

4.1 不同后掠角變化時的動態(tài)響應(yīng)

保持飛行器的展長不變,后掠角分別以4°/s、8°/s、10°/s的速度變化,得到的后掠角、高度、速度、攻角和俯仰角速度等運動參數(shù)的動力學(xué)響應(yīng)如圖5所示。

從圖5可以看出,由于后掠角的增加,導(dǎo)致氣動焦點、壓心均劇烈后移,使得氣動力產(chǎn)生的低頭力矩迅速增加,滾轉(zhuǎn)角速度增加,俯仰角減小,但彈道傾角變化相對緩慢,進而導(dǎo)致飛行器的攻角迅速降低。隨著后掠角的增加,展弦比減小,升力系數(shù)顯著降低,飛行器開始迅速掉高。注意到附加力矩是抬頭力矩,可知變形過程實際是氣動特性與質(zhì)量特性的博弈,但氣動參數(shù)的變化更加劇烈,因此飛行器低頭趨勢占據(jù)了主導(dǎo),同時需注意在設(shè)計時考慮放寬靜穩(wěn)定性。從圖中還可以看出,后掠角變化速率快時,使得各參數(shù)的變化趨勢提前發(fā)生,但總體影響不顯著。

圖5 不同后掠角變化速率下的動態(tài)響應(yīng)Fig.5 The dynamicre sponses at different morphingspeeds of sweep angle

4.2 不同展長變化時的動態(tài)響應(yīng)

保持飛行器的后掠角不變,展長分別以0.2 m/s、0.3 m/s、0.5 m/s的速度變化,得到的伸縮翼展長、高度、速度、攻角和俯仰角速度等運動參數(shù)的動力學(xué)響應(yīng)如圖6所示。

圖6 不同展長變化速率下的動態(tài)響應(yīng)Fig.6 The dynamic responses at different morphingspeeds of span

從圖6可以看出,由于展長的減小,飛行器呈現(xiàn)先增高后掉高的過程。主要因為在0°后掠角狀態(tài)展長變化時,質(zhì)心位置不發(fā)生改變,氣動作用中心變化也不顯著,但此時由于零升力矩系數(shù)的減小,導(dǎo)致飛行器呈現(xiàn)緩慢的抬頭趨勢,俯仰角速度緩慢增加,進而導(dǎo)致攻角增加。由于推力不夠,速度持續(xù)降低,升力減小,隨后飛行器開始掉高,飛行軌跡進入長周期的高度方向震蕩衰減過程。同時從圖中看出,相較于后掠角,動力學(xué)參數(shù)對展長的變化速率相對敏感。但增加展長的變化速率不會帶來惡劣影響,這對于使用展長的變化來控制飛行器姿態(tài)而言很有意義。

4.3 展長和后掠角同時變化時的響應(yīng)

后掠角和展長同時改變,變化速度為4組,分別為:6°/s,0.3 m/s;6°/s,0.5 m/s;8°/s, 0.3 m/s;8°/s,0.5 m/s.得到的高度、速度、攻角和俯仰角速度等各運動參數(shù)的動力學(xué)響應(yīng)如圖7所示。

圖7 不同后掠角和展長變化速率下的動態(tài)響應(yīng)Fig.7 The dynamicre sponses at different morphing speeds of sweep and span

從圖7可以看出,當(dāng)后掠角和展長同時變化時,后掠角對動態(tài)特性的影響占據(jù)主導(dǎo),但是展長的變化對結(jié)果起到了平衡和緩和的作用,進而可知通過合適的變形方式和變形速度的結(jié)合可以弱化變形過程中的動力學(xué)參數(shù)波動。

5 結(jié)論

本文基于Newton-Euler方法建立了變后掠變展長無人飛行器的多體動力學(xué)模型,推導(dǎo)了由變展長和變后掠角所產(chǎn)生的附加力和力矩項。在縱向上解耦簡化了動力學(xué)方程,基于準(zhǔn)定常氣動力模型,分析了不同變形形式和變形速率下的飛行器縱向動態(tài)響應(yīng)。仿真結(jié)果表明,后掠角和展長的變化對動力學(xué)特性影響較大,但變化速度對結(jié)果影響不大。通過合適的變形方式和變形速度的結(jié)合可以弱化變形過程中的動力學(xué)參數(shù)波動,進而減小控制系統(tǒng)的負擔(dān)。

References)

[1] Vasista S,Tong L Y,Wong K C.Realization of morphing wings: a multidisciplinary challenge[J].Journal of Aircraft,2012, 49(1):11-28.

[2] Smith K,Butt J,Von Spakovsky M R.A study of the benefits of using morphing wing technology in fighter aircraft systems[C]∥Proceedings of the 39th AIAA Thermophysics Conference.Miami, Florilia,US:AIAA,2007:25-28.

[3] William A.Crossley.Morphing aircraft sizing using design optimization[J].Journal of Aircraft.2011,48(2):612-622.

[4] 陳錢,尹維龍,白鵬,等.變后掠變展長翼身組合體系統(tǒng)設(shè)計與特性分析[J].航空學(xué)報,2010,31(3):506-513.

CHEN Qian,YIN Wei-long,BAI Peng,et al.System design and characteristics analysis of a variable-sweep and variable-span wingbody[J].Acta Aeronautica ET Astronautica Sinica,2010, 31(3):506-513.(in Chinese)

[5] Mestrinho J,Gamboa P,Santos P.Design optimization of a variable-span morphing wing for a small UAV[C]∥52nd AIAA/ ASME/AHS/ASC Structures,Structural Dynamics and Materials Conference.Denver,Colorado:AIAA,ASME,ASCE,AHS, ASC,2011:4-7.

[6] An J G,Yan M,Zhou W B,et al.Aircraft dynamic response to variable wing sweep geometry[J].Journal of Aircraft,1988, 25(3):216-221.

[7] 樂挺,王立新,艾俊強.Z型翼變體飛機的縱向多體動力學(xué)特性[J].航空學(xué)報,2010,31(4):679-686.

YUE Ting,WANG Li-xin,AI Jun-qiang.Longitudinal multibody dynamic characteristics of Z-wing morphing Aircraft[J]Acta Aeronautica ET Astronautica Sinica,2010,31(4):679-686. (in Chinese)

[8] Henry J J,Pines D J.A mathematical model for roll dynamics by use of a morphing-Span Wing[C]∥48th AIAA/ASME/AHS/ ASC Structures,Soructural Dynamics,and Materials Conference. Honolulu,Hawaii:AIAA,ASME,ASCE,AHS,ASC,2007:23-26.

[9] Fossen T I.Guidance and control of ocean vehicles[M].US: John Wiley&Sons Inc,1994:25-28.

Dynamic Modeling and Response of a Morphing UAV with Variable Sweep and Variable Span

YANG Guan-tong,TANG Sheng-jing,ZHAO Lin-dong,GUO Jie
(Key Laboratory of Dynamic and Control of Flight Vehicle,Ministry of Education,Beijing Institute of Technology,Beijing 100081,China)

For morphing flight vehicles,the changes of sweep angle and span have great influence on key parameters such as aerodynamic forces,inertial moments,pressure center and gravity center.A dynamic model of morphing vehicle with variable sweep and span is established based on Newton-Euler method. The additional force and moment generalized by morphing are deduced.The dynamic response aroused by morphing is studied using decoupled longitudinal equations and quasi-steady aerodynamic model.The results indicate that the changes of sweep and span influence the dynamic characteristics severely;however,the effect of different morphing speed is not remarkable.Appropriate strategy of morphing could weaken the fluctuation of dynamic parameters and relieve the burden of control system.

applied mechanics;UAV;variable sweep and variable span;dynamic modeling;dynamic response

V212.1

:A

1000-1093(2014)01-0102-06

10.3969/j.issn.1000-1093.2014.01.015

2013-04-22

航空科學(xué)基金項目(2012ZA72002)

楊貫通(1986—),男,博士研究生。E-mail:ygt@bit.edu.cn;

唐勝景(1959—),男,教授,博士生導(dǎo)師。E-mail:tangsj@bit.edu.cn