李偉杰, 沈作軍

(北京航空航天大學(xué) 航空科學(xué)與工程學(xué)院, 北京 100191)

?

基于μ分析的高超聲速飛行器再入軌跡評估

李偉杰, 沈作軍

(北京航空航天大學(xué) 航空科學(xué)與工程學(xué)院, 北京 100191)

為解決高超聲速飛行器再入過程中存在的不確定性問題,將μ分析理論應(yīng)用到標(biāo)稱再入軌跡的評估中。采用線性穩(wěn)定儲備準(zhǔn)則,在密度、升力系數(shù)和阻力系數(shù)存在不確定性時,對標(biāo)稱軌跡在制導(dǎo)律下的制導(dǎo)穩(wěn)定性進(jìn)行了分析;根據(jù)所獲得的最壞參數(shù)組合,采用系統(tǒng)的開環(huán)尼克爾斯曲線驗證了μ分析的結(jié)果;分析了系統(tǒng)所能容納的最大不確定性參數(shù),并研究了系統(tǒng)對不同不確定性參數(shù)的敏感性。結(jié)果表明:在密度、升力系數(shù)和阻力系數(shù)不確定性中,阻力系數(shù)對該標(biāo)稱軌跡的魯棒穩(wěn)定性影響最為明顯。

μ分析; 不確定性; 再入; 線性分式變換

0　引言

高超聲速飛行器是指馬赫數(shù)不低于5的飛行器[1]。高超聲速飛行器再入過程中,飛行包線大、飛行環(huán)境復(fù)雜、飛行高度及馬赫數(shù)跨度大,氣動特性變化劇烈;同時,由于缺乏足夠的飛行試驗,導(dǎo)致高超聲速飛行器再入過程中存在很大的不確定性。因此,必須在各種不確定性存在的條件下,對設(shè)計的標(biāo)稱再入軌跡及制導(dǎo)律進(jìn)行魯棒穩(wěn)定性評估。

在傳統(tǒng)的制導(dǎo)穩(wěn)定性評估中,常采用拉偏氣動升阻和環(huán)境參數(shù)的方法,通過再入制導(dǎo)飛行仿真來確定再入軌跡及其制導(dǎo)律對主要不確定參數(shù)的適應(yīng)性。工程上常采用的網(wǎng)格法將不確定性參數(shù)表達(dá)為參數(shù)不確定性并進(jìn)行網(wǎng)格化,在網(wǎng)格點針對評估準(zhǔn)則逐一評估。由于網(wǎng)格法僅在網(wǎng)格化的離散點進(jìn)行評估,因此可靠性得不到保證;此外,其工作量會隨不確定性參數(shù)的增加呈指數(shù)增加[2]?；诰€性分式變換(Linear Fractional Transformation, LFT)的μ分析方法克服了網(wǎng)格法的不足,能在連續(xù)的頻域范圍內(nèi)計算帶有不確定性系統(tǒng)的結(jié)構(gòu)奇異值,且其計算量是一個關(guān)于不確定參數(shù)個數(shù)的多項式函數(shù)[3]。

在采用μ分析方法進(jìn)行魯棒穩(wěn)定性評估時,其可靠性很大程度上依賴于線性分式變換的質(zhì)量?；谖锢砟Ｐ偷慕７椒?圖操作方法)能充分考慮不確定參數(shù)之間的關(guān)聯(lián),避免了數(shù)值建模方法導(dǎo)致的保守性[4-5]。

1　模型建立

1.1三自由度縱向再入模型

在再入問題中,常應(yīng)用瞬時平衡假設(shè)忽略姿態(tài)運(yùn)動及地球自轉(zhuǎn),得到北天東坐標(biāo)系中的高超聲速飛行器三自由度縱向動力學(xué)方程[6],如式(1)所示:

(1)

1.2基于物理模型的LFT建模方法

基于物理模型的LFT建模方法直接將不確定性參數(shù)引入飛行器非線性SIMULINK模型。以密度為例,假設(shè)密度可在其標(biāo)稱值附近±20%范圍內(nèi)變化,則密度可以表示為ρ=ρnom(1+WρΔρ)。其中ρnom為密度的標(biāo)稱值,Wρ值為0.2,Δρ∈[-1,1]。于是,包含不確定性參數(shù)的密度在模塊中可以表示為如圖1的形式,在密度不確定性處引入了一個虛擬輸出zρ和虛擬輸入wρ。這樣,便直接將不確定參數(shù)引入到飛行器的非線性模型中,且由于考慮到了不確定參數(shù)之間可能存在的關(guān)聯(lián),不會引起數(shù)值建模方法帶來的保守性問題[4,7]。

圖1　基于物理模型的建模方法Fig.1　Modeling method based on physical models

1.3線性穩(wěn)定儲備準(zhǔn)則

圖2(a)所示的尼克爾斯圖所對應(yīng)的區(qū)域A和B分別對應(yīng)±6 dB,±36.87°和±4.5 dB,±28.44°的儲備要求。若系統(tǒng)的開環(huán)尼克爾斯曲線不穿越橢圓區(qū)域,則系統(tǒng)滿足相應(yīng)的幅值和相位儲備要求[7]。由文獻(xiàn)[8]可知,圖2(b)中區(qū)域A和B分別對應(yīng)奈奎斯特圖中以-1.25為圓心、0.75為半徑和以-1.14為圓心、0.54為半徑的圓。系統(tǒng)的開環(huán)尼克爾斯曲線不穿越橢圓區(qū)域B,等價于系統(tǒng)的奈奎斯特曲線不穿越圓區(qū)域B。

圖2　尼克爾斯橢圓區(qū)域和奈奎斯特圓區(qū)域Fig.2　Elliptical regions in Nichols plane and circular regions in Nyquist plane

在分析制導(dǎo)穩(wěn)定性時,通過構(gòu)造一個虛擬的復(fù)不確定塊,將準(zhǔn)則作為飛機(jī)對象虛擬的乘積型不確定性考慮,加入閉環(huán)系統(tǒng)。本文采用圖2中區(qū)域B對應(yīng)的儲備準(zhǔn)則。由文獻(xiàn)[8-9]可知,所有包含評估準(zhǔn)則的不確定對象的集合為:

(2)

(3)

圖3　引入幅值和相位裕度的反饋閉環(huán)系統(tǒng)Fig.3　Feedback closed-loop system with gain and phase margins introduced

2　結(jié)構(gòu)奇異值理論

設(shè)結(jié)構(gòu)Δ∈Cn×n為復(fù)數(shù)對角塊,可描述為:

(4)

式中:S和F為重復(fù)標(biāo)量子塊和子滿塊的個數(shù)。

(5)

當(dāng)不存在這樣的攝動使det(I-MΔ)=0時,則定義μΔ(M)=0。

定理1:假設(shè)M∈RH∞,對于所有滿足‖Δ‖∞≤1的Δ,圖4所示系統(tǒng)魯棒穩(wěn)定的充要條件為:

(6)

圖4　線性分式變換反饋框圖Fig.4　Feedback block diagram of LFT

3　仿真計算及結(jié)果分析

3.1仿真計算

假設(shè)再入過程中密度、升力系數(shù)和阻力系數(shù)存在±20%不確定參數(shù)。按照基于物理模型建模原則搭建非線性動力學(xué)模型,并考慮±4.5 dB,±28.44°的線性穩(wěn)定儲備,利用Linmod函數(shù)可方便地得到含不確定性系統(tǒng)的狀態(tài)空間實現(xiàn)。至此,便得到了關(guān)于密度、升力系數(shù)、阻力系數(shù)不確定參數(shù)的線性分式變換,且將線性穩(wěn)定儲備引入了原系統(tǒng)中。同時,復(fù)不確定參數(shù)的引入有效地解決了MATLAB無法計算實不確定參數(shù)的結(jié)構(gòu)奇異值的上、下界的問題[11]。計算出的結(jié)構(gòu)奇異值的上界可用于判斷魯棒穩(wěn)定性,下界可用于確定“最壞情況”。

已知X-33的一條標(biāo)稱軌跡,相關(guān)的三自由度縱向狀態(tài)信息和標(biāo)稱控制量如圖5所示。

圖5　X-33再入軌跡標(biāo)稱狀態(tài)量和控制量Fig.5　Nominal statement profiles and control profiles of X-33 reentry trajectory

從σ剖面可以看出,再入過程中飛行器存在兩次明顯的翻轉(zhuǎn)。翻轉(zhuǎn)點也是再入過程中較為危險的點。本文以第一次翻轉(zhuǎn)前后的狀態(tài)點為例進(jìn)行分析。

3.2魯棒穩(wěn)定性分析

假設(shè)ρ,CL和CD不確定性取值范圍為[-20%,+20%],分別在r,V,γ和S2go通道加入線性穩(wěn)定儲備,得到第一次翻轉(zhuǎn)前后狀態(tài)點的結(jié)構(gòu)奇異值頻域響應(yīng)如圖6所示。由定理1可知,標(biāo)稱軌跡在給定不確定性參數(shù)攝動下滿足魯棒穩(wěn)定性要求。

圖6　第一次翻轉(zhuǎn)前后結(jié)構(gòu)奇異值Fig.6　μ values before and after the first bank-reversal

為了驗證整條標(biāo)稱軌跡在制導(dǎo)律閉環(huán)后存在不確定參數(shù)時的魯棒穩(wěn)定性,可以更大范圍地選取軌跡上的點進(jìn)行驗證。以r通道為例,其結(jié)果見圖7。

圖7　第一次翻轉(zhuǎn)附近區(qū)域r通道結(jié)構(gòu)奇異值Fig.7　μ values of r loop around the first bank-reversal point

結(jié)果表明,r通道的結(jié)構(gòu)奇異值均小于1,說明標(biāo)稱軌跡r通道滿足魯棒穩(wěn)定性要求。

3.3最壞情況驗證

為了驗證μ分析的結(jié)果,現(xiàn)以第一次翻轉(zhuǎn)前狀態(tài)點的r通道為例加以驗證,其結(jié)果如圖8所示。

由圖8可知,在所求得的最壞情況下,系統(tǒng)的開環(huán)尼克爾斯曲線不穿越圖2(b)中區(qū)域B所對應(yīng)的穩(wěn)定儲備區(qū)域,即系統(tǒng)在密度、升力系數(shù)、阻力系數(shù)具有[-20%,+20%]不確定參數(shù)變化范圍內(nèi)的最壞情況仍滿足魯棒穩(wěn)定性,從而驗證了μ分析的結(jié)論。

3.4穩(wěn)定性邊界參數(shù)確定及“敏感性”分析

在對標(biāo)稱軌跡進(jìn)行制導(dǎo)穩(wěn)定性分析時,總希望求得其所能容納的最大不確定性參數(shù)范圍?，F(xiàn)以第一次翻滾前r通道為例,在原有的20%基礎(chǔ)上,按一定比例(在此選4%)同時增大3個不確定參數(shù),其結(jié)構(gòu)奇異值上界如圖9所示?？梢钥闯?隨著不確定參數(shù)的增加,閉環(huán)系統(tǒng)的結(jié)構(gòu)奇異值逐漸增加。當(dāng)不確定參數(shù)達(dá)到約35%時,閉環(huán)系統(tǒng)臨界穩(wěn)定。

圖9　不同不確定參數(shù)對應(yīng)的結(jié)構(gòu)奇異值上界Fig.9　The upper boundard of μ corresponding to different uncertain parameters

為了分析密度、升力系數(shù)、阻力系數(shù)三個不確定參數(shù)對閉環(huán)系統(tǒng)的標(biāo)稱軌跡穩(wěn)定性的影響程度,在原有的20%不確定參數(shù)的基礎(chǔ)上,固定其中兩個不確定參數(shù),增加另外一個不確定參數(shù)。通過計算系統(tǒng)的結(jié)構(gòu)奇異值,來檢驗閉環(huán)系統(tǒng)的穩(wěn)定性,從而分析不確定參數(shù)的“敏感性”,計算結(jié)果如圖10所示。

圖10　不確定參數(shù)對系統(tǒng)穩(wěn)定性的影響Fig.10　Influence of uncertain parameters on the system stability

由圖10可知:當(dāng)密度和升力系數(shù)不確定性分別達(dá)到80%和92%時,閉環(huán)系統(tǒng)臨界穩(wěn)定;而當(dāng)阻力系數(shù)不確定性參數(shù)增加到40%時,閉環(huán)系統(tǒng)便失穩(wěn)。因此可知,閉環(huán)系統(tǒng)的魯棒穩(wěn)定性對阻力系數(shù)不確定性參數(shù)更為敏感,因此在標(biāo)稱軌跡設(shè)計時需要著重考慮。

4　結(jié)束語

本文將結(jié)構(gòu)奇異值理論應(yīng)用到高超聲速飛行器再入過程中存在參數(shù)不確定性時標(biāo)稱軌跡的魯棒穩(wěn)定性分析當(dāng)中。應(yīng)用該方法可對用不同方法設(shè)計的標(biāo)稱再入軌跡在存在參數(shù)不確定性時進(jìn)行魯棒穩(wěn)定性分析,從而為高超聲速飛行器再入軌跡的設(shè)計及制導(dǎo)律評估提供一定的參考。

[1]吳宏鑫,孟斌.高超聲速飛行器控制研究綜述[J].力學(xué)進(jìn)展,2009,39(6):756-765.

[2]劉林.現(xiàn)代飛行控制系統(tǒng)的評估與確認(rèn)方法[M].北京:國防工業(yè)出版社,2010:76-100.

[3]Juliana S.Re-entry flight clearance[M].Delft:Delft University of Technology,2006:15-18.

[4]Kureemun R,Bates D G,Hayes M J.On the generation of LFT-based uncertainty models for flight control law robustness analysis[R].AIAA-2001-4396,2001.

[5]Zerz E.LFT representations of parametrized polynomial systems[J].IEEE Transactions on Circuits and Systems,1999,46(3):410-416.

[6]Vinh N X,Busemann A,Culp R D.Hypersonic and planetary entry flight mechanics[M].MI, Ann Arbor:University of Michigan Press,1980:26-27.

[7]Bates D G,Kureemun R,Mannchen T.Improved clearance of a flight control law using-analysis techniques[J].Journal of Guidance,Control,and Dynamics,2003,26(6):869-884.

[8]Deodhare G,Patel V V.A "modern" look at gain and phase margins:anH∞/μapproach[R].AIAA-98-4134,1998.

[9]Kureemun R,Bates D G,Postlethwaite I.Quantifying the robustness of flight control systems using Nichols exclusion regions and the structured singular value [J].Proceedings of the Institution of Mechanical Engineers Part Ⅰ Journal of Systems and Control Engineering,2001,215(6):625-638.

[10]Kemin Zhou.Essentials of robust control[M].Louisiana:Louisiana State University,1999:165-220.

[11]Balas G J,Doyle J G,Glover K.MATLABμ-analysis and synthesis toolbox user’s guide[M].Version 3.USA:The Math Works Inc,1998:61-66.

(編輯：李怡)

Evaluation of re-entry trajectory of hypersonic vehicles based onμ-analysis

LI Wei-jie, SHEN Zuo-jun

(School of Aeronautic Science and Engineering, BUAA, Beijing 100191, China)

Theμ-analysis theory was applied to the evaluation of the nominal re-entry trajectory to solve the uncertainties in the course of re-entry of hypersonic vehicles. Based on the linear stability reserve criteria, the guidance stability of the nominal trajectory and guidance law was analyzed when the density, lift coefficient and drag coefficient are uncertainty. According to the worst-case combination of parameters, the Nichols curve of the open-loop system was used to verify the conclusion ofμ-analysis results. The maximum uncertainty parameters which the system could accommodate were analyzed, and the sensitivity of the system to different uncertain parameters was studied. The results indicate that the drag coefficient has the most obvious effect on the robust stability of the nominal trajectory in the uncertainties of the density, lift coefficient and drag coefficient.

μ-analysis; uncertainties; re-entry; linear fractional transformation

2015-09-14;

2016-01-06; 網(wǎng)絡(luò)出版時間:2016-02-29 16:37

李偉杰(1991-),男,四川簡陽人，碩士研究生,研究方向為高超聲速飛行器制導(dǎo)與控制。

V412.4

A

1002-0853(2016)04-0050-04