薛源, 徐浩軍, 張建林, 葉蔚

(1.空軍工程大學(xué) 航空航天工程學(xué)院, 陜西 西安 710038;2.中國人民解放軍94270部隊, 山東 濟南 250114)

0 引言

徐浩軍(1965-)，男，浙江余姚人，教授，博士生導(dǎo)師，研究方向為飛行安全與作戰(zhàn)效能評估。

在無人機系統(tǒng)越來越強調(diào)體系協(xié)同能力的背景下,無人機編隊飛行仿真技術(shù)成為一個重要的研究方向,擁有廣闊的發(fā)展前景[1-2]。通過無人機編隊飛行仿真研究可以拓寬無人機的應(yīng)用范圍、提高執(zhí)行任務(wù)效能、增加系統(tǒng)冗余度,使其較好地適應(yīng)復(fù)雜環(huán)境,完成單機無法自主完成的任務(wù)。

前機尾流是研究無人機編隊問題需要考慮的一個重要影響因素。國外學(xué)者較早認(rèn)識到編隊情形下尾流場對后機氣動性能及安全性的影響,并進(jìn)行了相關(guān)研究[3-6]。國內(nèi)研究人員對這一問題的認(rèn)識還不夠,定量研究不足,對無人機編隊情形下長機尾流場對僚機氣動影響的研究較少。研究工作大多是分析尾流對后機的氣動影響,并沒有涉及編隊飛行控制算法。針對編隊飛行控制算法的研究大多是基于理想無尾流擾動狀況[7-10],很少涉及前機尾流對后機的氣動影響或者研究不深,更少有研究編隊視景仿真技術(shù)的工作;因此有必要將前機尾流引入到飛行仿真與跟蹤控制的迭代過程中,并開展視景仿真技術(shù)研究。

本文基于分布式仿真框架[11]構(gòu)建視景仿真系統(tǒng),通過提供通用的、相對獨立的支撐軟件將仿真應(yīng)用與底層通訊功能分離,使各部分運用各自的技術(shù),最大限度地提高系統(tǒng)的通用性以及各分系統(tǒng)的互操作性和重用性。本文無人機編隊視景及態(tài)勢采用基于OPENGL的圖形建模技術(shù)進(jìn)行開發(fā),并采用衛(wèi)星圖片與航空照片、貼圖與三維場景模型相結(jié)合的地景構(gòu)建方式,視景仿真中的大地景精度較高,真實感較強。

1 尾流影響下的無人機編隊飛行建模

1.1 長機尾流模型

目前,已有多種尾流場中的滾轉(zhuǎn)效應(yīng)模型,這些模型對尾流場中后機的飛行動力學(xué)仿真是可以借鑒的。本文采用Proctor提出的非奇尾渦模型[12],有側(cè)風(fēng)情況下長機在后方兩個展長處的尾流速度場模型如圖1所示。

圖1 尾流場內(nèi)的速度分布Fig.1 The velocity distribution of wake vortex

1.2 編隊飛行動力學(xué)建模

長機尾流對僚機有效風(fēng)分量及有效風(fēng)梯度的計算采用構(gòu)型積分算法,忽略長機尾流對尾翼、襟副翼、設(shè)備艙、掛載等結(jié)構(gòu)外形的氣動影響。飛機本體六自由度微分方程解算的參數(shù)為:

x=[V,α,β,q0,q1,q2,q3,p,q,r,x,y,z] (1)

所得參數(shù)經(jīng)相應(yīng)的坐標(biāo)變換、單位轉(zhuǎn)換、微積分變換后得到所需的其他飛行參數(shù),[q0,q1,q2,q3]T即為姿態(tài)角轉(zhuǎn)換而來的四元數(shù)。

使用Real Time Workshop(RTW)將SIMULINK搭建的運動方程及控制律轉(zhuǎn)化為實時系統(tǒng)VxWorks支持的C代碼,并將其下載到實時仿真機,仿真步長為20 ms。

1.3 基于改進(jìn)線性二次型的編隊保持算法

2 構(gòu)建分布式視景仿真系統(tǒng)

在對尾流情形下的無人機編隊飛行建模的基礎(chǔ)上,根據(jù)分布式仿真系統(tǒng)建模路線[11],構(gòu)建無人機編隊飛行分布式仿真框架,如圖2所示。仿真計算延遲小于20 ms,幀速率為50 Hz,仿真單元間的同步誤差不高于10 ms。該系統(tǒng)具有開放式分布仿真體系結(jié)構(gòu),可以根據(jù)需要任意擴展或重載具有分布式體系結(jié)構(gòu)的仿真單元,仿真單元間可以互操作,提高了仿真組件的重用性。

各單元的特點及功能如下:

初始化單元:通過Graphical User Interface (GUI)設(shè)置初始仿真參數(shù),并對初始編隊狀態(tài)進(jìn)行配平,通過仿真總線向各單元發(fā)送初始化命令及成員初始化信息;然后對仿真中的無人機任務(wù)進(jìn)行規(guī)劃并將仿真所有權(quán)交予編隊飛行實時計算單元。

編隊飛行計算單元:編隊飛行過程中,飛行動力學(xué)仿真單元根據(jù)規(guī)劃的無人機指令及地面站操縱信號解算飛機運動方程,僚機通過改進(jìn)的動態(tài)線性二次型算法跟蹤長機的飛行狀態(tài)。

硬件測控單元:響應(yīng)無人機任務(wù)狀態(tài),利用仿真總線中的飛參數(shù)據(jù)對模擬器轉(zhuǎn)臺、舵面等硬件設(shè)備進(jìn)行驅(qū)動,最終表現(xiàn)為舵面和油門信號的變化。測控系統(tǒng)將從硬件系統(tǒng)采集的舵面和油門信號以數(shù)組的格式送入Run Time Infrastructure (RTI) 仿真總線。

虛擬儀表單元:從仿真總線接收仿真計算的數(shù)據(jù),轉(zhuǎn)發(fā)給各虛擬儀表,驅(qū)動儀表參數(shù)的變化,從而更加直觀、形象地觀察無人機的姿態(tài)、位置等參數(shù)的變化。

視景單元:視景單元主要包括三維視景、二維態(tài)勢、三維態(tài)勢三個子單元。當(dāng)視景單元從總線接到仿真開始交互類后開始等待接收數(shù)據(jù)。一旦接收到仿真所有權(quán),視景單元根據(jù)接收的飛參數(shù)組對二維和三維場景中的對象進(jìn)行控制,并將仿真場景效果實時渲染到最終的顯示設(shè)備上。

三維視景為我國某地500 km×500 km范圍的真實地景;整個飛行場景區(qū)域紋理精度為5 m,高程精度為15 m;精細(xì)地面目標(biāo)包括12個機場模型;機場起飛著陸區(qū)域小于15 km×15 km,紋理精度為1 m。

二維態(tài)勢單元接收仿真總線飛行參數(shù)的飛機位置和姿態(tài)信息,并對真實情況下的二維地圖進(jìn)行實時刷新,不斷更新顯示態(tài)勢信息,勾畫出飛機的飛行航跡。

三維態(tài)勢單元對仿真區(qū)域內(nèi)的真實環(huán)境進(jìn)行模擬,并在三維地圖上顯示飛機航跡和態(tài)勢的實時變化情況。

數(shù)據(jù)庫單元:當(dāng)數(shù)據(jù)庫從總線接到仿真開始交互類后開始等待接收數(shù)據(jù),一旦接收到仿真所有權(quán)后即負(fù)責(zé)所有仿真數(shù)據(jù)的記錄及存儲。

數(shù)據(jù)回放單元:對整個仿真過程進(jìn)行二維曲線或三維視景的回放。

評估反饋單元:當(dāng)從仿真總線中接到仿真結(jié)束的交互類后,評估整個編隊飛行過程,分析飛參數(shù)據(jù)。

RTI是分布式接口規(guī)范的具體實現(xiàn),在系統(tǒng)運行過程中為各單元的“同步”和“數(shù)據(jù)交換”提供公共接口。

圖2 無人機編隊飛行分布式仿真系統(tǒng)Fig.2 Distributed simulation system for UAV formation flight

3 YOYO機動仿真及分析

本文對典型的空中戰(zhàn)術(shù)動作YOYO機動進(jìn)行模擬仿真。YOYO機動對于無人機規(guī)避敵方攻擊、撤離任務(wù)區(qū)域、追蹤移動目標(biāo)等行為具有重要的戰(zhàn)術(shù)意義,其實質(zhì)是短時間內(nèi)勢能和動能的轉(zhuǎn)化。在無人機執(zhí)行任務(wù)的過程中,轉(zhuǎn)彎半徑隨著飛機速度的增大而增大;因此,在短時間內(nèi)獲得較大的轉(zhuǎn)彎角度并建立無人機的指向優(yōu)勢,最有效的手段是拉升一定的高度,使勢能增大、動能減小,然后向轉(zhuǎn)彎方向滾轉(zhuǎn)并俯沖下降。

仿真試驗地區(qū)設(shè)置在我國東南某地。試驗起始點在某機場上空4 000 m,長僚機編隊以速度100 m/s保持配平狀態(tài)下的勻速飛行,試驗開始時的經(jīng)緯度為(118.583 128°,24.787 605°)。編隊狀態(tài)為典型的V字形編隊,僚機位置相對于長機位置在地軸系x,y,z軸上的增量為[-10 m,-8 m,0 m]。

仿真開始時,雙機在4 000 m高度以100 m/s的速度保持平飛狀態(tài)。YOYO機動開始時,長機首先下降到高度3 959.52 m,此時速度為107.75 m/s;右滾轉(zhuǎn)的同時增大飛機的俯仰角使飛機爬升,當(dāng)仿真時間為99.00 s時速度為67.08 m/s,高度達(dá)到4 475.32 m;之后調(diào)整升降舵使飛機迅速低頭并增大滾轉(zhuǎn)角,飛機開始快速轉(zhuǎn)彎并急劇掉高度,當(dāng)仿真時間為136.00 s時速度為146.88 m/s,高度為3 724.49 m,偏航角達(dá)到180°,完成了180°轉(zhuǎn)向;此后,長機繼續(xù)右滾轉(zhuǎn),當(dāng)仿真時間達(dá)到151.40 s時,速度達(dá)到本次仿真的最大值164.93 m/s,高度為3 333.21 m。在長機執(zhí)行YOYO機動的過程中,僚機跟蹤長機狀態(tài)進(jìn)行編隊伴隨飛行。在155.81 s時仿真結(jié)束,此時經(jīng)緯度坐標(biāo)為(118.613 039°,4.814 918°),仿真過程及結(jié)果如圖3～圖5所示。

圖3為仿真過程的飛行軌跡,可以看出雙機編隊經(jīng)YOYO機動后偏航角增大了180°以上,獲得了較大的指向優(yōu)勢。無人機編隊通過YOYO機動使偏航角變化180°所需的時間為136 s,比定常盤旋變化180°偏航角所需的200～300 s約快一倍。圖4為本次仿真過程中無人機編隊在不同時刻的視景截圖。圖5為仿真過程中長僚機主要參數(shù)的變化情況,僚機跟蹤長機各參數(shù)的均方差不超過5%。

圖3 無人機編隊進(jìn)行YOYO機動的飛行軌跡Fig.3 The flight path during YOYO maneuver

圖4 無人機編隊進(jìn)行YOYO機動的3D仿真視景Fig.4 3D visual simulation of UAV formation flightduring YOYO maneuver

圖5 仿真過程中的長僚機參數(shù)變化Fig.5 Flight parameters of UAVs during the simulation

由圖5可以看出,在初始的配平狀態(tài),僚機的副翼偏角為-1.362°。這是由于長機尾流對僚機產(chǎn)生了滾轉(zhuǎn)力矩增量,使得僚機有向左滾轉(zhuǎn)的趨勢,故僚機副翼偏轉(zhuǎn)以配平尾流引起的滾轉(zhuǎn)力矩。初始配平狀態(tài)時長機升降舵偏角為-1.739°、迎角為4.493°、俯仰角為4.493°,而僚機升降舵偏角為-0.820°、迎角為2.770°、俯仰角為2.770°,這是由于長機尾流對僚機的升力系數(shù)產(chǎn)生了增量,使得僚機有向上的升力增量,故僚機不必保持和長機相同的迎角來獲得配平升力。綜上所述,如果編隊飛行的配平狀態(tài)不考慮尾流影響,將會使得計算所得僚機操縱面偏角與有尾流影響的真實情況存在較大誤差。因此,在工程應(yīng)用中應(yīng)考慮長機尾流誘導(dǎo)的升力增量和滾轉(zhuǎn)力矩作用,通過飛控動態(tài)調(diào)整舵面以保持編隊飛行狀態(tài)。

通過對圖5中飛行參數(shù)的分析,可看出在YOYO機動的過程中,長機尾流對僚機的影響不是一個持續(xù)的過程,其影響程度主要取決于長僚機相對姿態(tài)角的大小。在YOYO機動的初始階段和結(jié)束階段長僚機的相對姿態(tài)較小,僚機受到尾流的影響較大。例如在YOYO機動開始之后,長機升降舵偏轉(zhuǎn)使得俯仰角從初始的4.493°開始減小,在此過程中,長僚機的相對姿態(tài)角減小,長機尾流對僚機的作用增強,僚機的升力系數(shù)顯著增加。而僚機副翼偏角此時由于尾流作用的加強而增大偏轉(zhuǎn)角度。在67 s后,隨著滾轉(zhuǎn)角的增大,長機尾流軌跡的位置偏離僚機機翼表面,使得尾流對僚機的影響逐漸減弱,觀察圖5亦可發(fā)現(xiàn),67 s后僚機和長機副翼、升降舵的響應(yīng)較為接近,姿態(tài)角的變化情況也較為一致。

YOYO機動的仿真結(jié)果表明了所設(shè)計的控制系統(tǒng)具有較強的魯棒性和自適應(yīng)跟蹤性能。視景單元在YOYO機動過程中表現(xiàn)出了較好的穩(wěn)定性和可靠性,對編隊飛行的模擬效果較為形象直觀。分布式系統(tǒng)中飛參數(shù)據(jù)由計算單元向視景單元傳輸?shù)耐叫院蛯崟r性亦得到驗證。

4 結(jié)束語

本文構(gòu)建了編隊飛行時長僚機動力學(xué)模型和尾流場模型,設(shè)計了僚機跟蹤控制算法。將尾流效應(yīng)考慮到無人機編隊飛行的仿真計算當(dāng)中,使其更加符合真實的編隊飛行狀況;建立了真實地理信息下高保真度的視景單元,其沉浸感和參與感較強;定量分析了YOYO機動過程中僚機受到的尾流影響和對長機狀態(tài)的跟蹤情況。分布式無人機編隊視景仿真方法具有較高的保真度,可以實現(xiàn)無人機編隊飛行研究中的宏觀與微觀、定性分析與定量分析相結(jié)合;可作為預(yù)研性的仿真試驗工具,同時也為其他類似的飛行仿真試驗系統(tǒng)的設(shè)計提供了一定的參考經(jīng)驗。

參考文獻(xiàn)：

[1] 樊瓊劍,楊忠,方挺,等.多無人機協(xié)同編隊飛行控制的研究現(xiàn)狀[J].航空學(xué)報,2009,30(4):683-691.

[2] 李文皓,張珩.無人機編隊飛行技術(shù)的研究現(xiàn)狀與展望[J].飛行力學(xué),2007,25(1):9-11.

[3] Ning S A,Flanzer T C,Kroo I M.Aerodynamic performance of extended formation flight[J].Journal of Aircraft,2011,48(3):855-865.

[4] Dogan A,Venkataramanan S.Modeling of aerodynamic coupling between aircraft in close proximity[J].Journal of Aircraft,2005,42(4):941-955.

[5] Blake W,Gingras D.Comparison of predicted and measured formation flight interference effects[J].Journal of Aircraft,2004,41(2):201-207.

[6] Jake V M,Ray R J,Ennix K,et al.F/A-18 performance benefits measured during the autonomous formation flight project[R].NASA/TM-2003-210732,2003.

[7] 朱杰斌,秦世引.無人機編隊飛行的分布式控制策略與控制器設(shè)計[J].智能系統(tǒng)學(xué)報,2010,5(5):392-399.

[8] Binetti P,Ariyur K B,Krstic M,et al.Formation flight optimization using extremism seeking feedback[J].Jouranal of Guidance,Control,and Dynamics,2003,26(1):132-142.

[9] 肖亞輝,王新民,王曉燕,等.無人機三維編隊飛行模糊PID控制器設(shè)計[J].西北工業(yè)大學(xué)學(xué)報,2011,29(6):834-838.

[10] 宋敏,魏瑞軒,沈東,等.基于非線性動態(tài)逆的無人機編隊協(xié)同控制[J].控制與決策,2011,26(3):448-452.

[11] 韓守鵬.分布式仿真系統(tǒng)動態(tài)重構(gòu)技術(shù)研究[D].長沙:國防科學(xué)技術(shù)大學(xué),2007.

[12] Gerz T,Holz?pfel F,Darracq D.Commercial aircraft wake vortices[J].Progress in Aerospace Sciences,2002,38(3):181-208.