升力式飛行器最短返回時間離軌點設(shè)計方法

丁 寧 李兆亭 張洪波

國防科技大學(xué)空天科學(xué)學(xué)院，長沙 410073

0 引言

在載人空間飛行、天地往返運輸?shù)瓤臻g任務(wù)結(jié)束后，需要空天飛行器返回預(yù)定著陸場，離軌制動是返回的第一步。離軌制動是指飛行器通過施加推力矢量改變原來的運行軌道，從而進入地球大氣層的軌道改變過程。離軌點設(shè)計是根據(jù)飛行器運行軌道和任務(wù)參數(shù)，確定星下點軌跡和可達域范圍，通過選擇合適的軌道位置進行離軌制動，使飛行器最終滿足再入點約束條件。

針對再入飛行器返回過程中的再入段，國內(nèi)外學(xué)者進行了諸多的研究，包括再入軌跡優(yōu)化方法[1-2]再入制導(dǎo)方法[3-4]、再入姿態(tài)控制方法[5-6]等，而對再入飛行器離軌段的軌道設(shè)計問題研究相對較少。早期研究主要將終端狀態(tài)約束設(shè)定為范圍約束，因此不同目標(biāo)函數(shù)下最優(yōu)解對應(yīng)的大氣進入狀態(tài)存在差異[7]。后期的研究大多將約束設(shè)定為固定的再入點位置和再入?yún)?shù)。其中，史樹峰等[8]通過規(guī)劃離軌制動時機和離軌控制參數(shù)，實現(xiàn)指定再入點位置的飛行器離軌任務(wù)；龔宇蓮等[9]設(shè)計了一種航天器自主離軌制動控制算法，可根據(jù)再入點狀態(tài)約束確定離軌過渡段的平均軌道根數(shù)，實時預(yù)報再入點航跡傾角；郭付明等[10]綜合考慮優(yōu)化精度和運行時間，提出了一種離軌點尋優(yōu)方法，可以降低軌道機動過程中對不確定因素的敏感性。然而上述各成果研究對象多為傳統(tǒng)彈道式飛行器，在設(shè)計過程中往往忽略橫向機動能力對離軌點的影響。而升力式飛行器升阻比大，再入過程中橫向機動能力強[11]，當(dāng)再入橫程確定時，相比于彈道式飛行器，縱程的取值范圍較大，對應(yīng)再入點的位置范圍較大。基于此，綜合考慮再入點位置變化，有必要對約束條件進行設(shè)置，以優(yōu)化離軌段時間，拓展離軌窗口范圍。

本文考慮升力式飛行器的橫向機動能力，提出了一種最短時間返回的離軌點設(shè)計策略。首先，根據(jù)星下點軌跡與著陸場的位置關(guān)系及各階段航程角的幾何關(guān)系，分別建立了位置約束和時間約束模型；然后根據(jù)給定的再入角和離軌能量約束，利用二體軌道理論推導(dǎo)了離軌點和再入點速度傾角、地心距及航程角之間的關(guān)系，并提出了一種基于牛頓迭代公式的離軌航程角尋優(yōu)方法；最后計算返回段總時間，采用非線性優(yōu)化方法求解最短時間返回軌道，確定離軌點。

1 最短返回時間離軌點設(shè)計流程

再入飛行器從運行軌道到達指定著陸場，依次經(jīng)歷等待段、離軌制動段、過渡段、再入段和著陸段。將離軌制動與過渡段合稱離軌段。

如圖1所示，Oe為地心，O，K，E和P分別為飛行器當(dāng)前位置、離軌點、再入點和著陸點，rk和re分別為離軌點和再入點地心距，Δtw，Δtk和Δte分別為等待段、離軌段和再入段時間，Δβw，Δβk和Δβe分別為等待段、離軌段和再入段航程角，Θk和Θe分別為離軌制動結(jié)束后的當(dāng)?shù)厮俣葍A角和再入角。離軌制動段為有限推力制動，但因時間較短，可簡化，以脈沖推力代替。在此基礎(chǔ)上，本文設(shè)計滿足位置、時間及過程約束的離軌軌道。計算步驟可歸納為：

圖1 再入飛行器返回過程示意圖

1)根據(jù)飛行器當(dāng)前軌道根數(shù)[a,e,i,ω,Ω,M]，確定一個周期內(nèi)的星下點軌跡；

2)計算著陸點與星下點軌跡的最小橫向距離dmin、著陸點與飛行器當(dāng)前位置的縱向航程角Δβ；

3)根據(jù)飛行器總體參數(shù)計算再入可達域，判斷是否滿足橫向位置約束，若滿足，則可得到不同橫、縱程對應(yīng)的再入時間Δte；若不滿足，將初始時刻后推一周期重新計算；

4)由軌道根數(shù)計算得到rk，并基于給定的re，Θe和能量約束Δvmax，計算不同增益速度Δv對應(yīng)的離軌時間Δtk，尋優(yōu)得到最小離軌航程角Δβkmin，進而計算由當(dāng)前時刻至著陸場總時間Δt；

5)對Δt尋優(yōu)，得到最短時間Δtmin，計算并輸出最優(yōu)離軌點軌道根數(shù)。

其計算流程圖如圖2所示，下面幾節(jié)分別敘述其計算方法。

圖2 最短返回時間離軌點設(shè)計流程

2 基于再入可達域的時空約束建模與分析

2.1 再入可達域條件設(shè)計

本文以最小返回時間為離軌點設(shè)計指標(biāo)，求解過程中需根據(jù)再入段時間、航程角等參數(shù)確定位置和時間約束，這些參數(shù)可由再入可達域提供。本節(jié)對再入可達域要提供的參數(shù)形式進行說明，以滿足后續(xù)設(shè)計需求。

再入可達域可描述不同縱程條件下的最大橫程，通常為橢圓形或扇形區(qū)域，常在換極坐標(biāo)系中表示。如圖3所示，λpf和φpf分別為落點的經(jīng)度和緯度，Z為縱程，L為橫程。

圖3 可達域示意圖

再入可達域的計算有很多經(jīng)典算法[12-14]，本文不作討論。在進行離軌點設(shè)計時，首先根據(jù)飛行器總體參數(shù)和再入約束對飛行可達域進行計算，扣除再入制導(dǎo)預(yù)留的射程機動余量，得到最大橫程Lmax。將最大縱程等間隔劃分為m個區(qū)域，每區(qū)域橫向n等分，得到m×n個網(wǎng)格，如圖4所示。取網(wǎng)格點為樣本點，以數(shù)表的形式給出每個樣本點的再入縱程Z、橫程L及時間Δte，作為離軌點設(shè)計的輸入量給到系統(tǒng)中，如表1所示。

圖4 可達域網(wǎng)格劃分

表1 可達域輸入?yún)?shù)表

在輸入樣本點再入時間參數(shù)后，對于橫、縱程分別為L和Z的著陸點，可通過拉格朗日插值多項式計算再入時間擬合值

Δte(L,Z)=

(1)

以及該落點對應(yīng)的再入航程角

(2)

式中:Re為地球平均半徑。

2.2 橫向位置約束

圖5 著陸點與星下點軌跡的最小橫向距離示意圖

(3)

實際上，由于地球自轉(zhuǎn)，在協(xié)議地心慣性坐標(biāo)系中，E點方位角是不停變化的。按下述公式進行迭代

(4)

則最小橫向距離

dmin=fd(λE(n),φE)

(5)

2.3 時間約束

由于燃料有限、再入橫程確定，離軌段和再入段時間不可能無限小，因此需要考慮最短時間帶來的時間約束。以當(dāng)前時刻作為起始時刻t0，設(shè)從該時刻至著陸點的飛行時間為Δt，則滿足

Δt=Δtw+Δtk+Δte

(6)

時間約束可表示為

Δtw≥0

(7)

即Δt≥Δtkmin+Δtemin。式中Δtkmin為離軌段最小時間，其計算方法將在3.2節(jié)給出，Δtemin=Δte(dmin,Zmin)為再入段最小時間，可由式求得。當(dāng)Δt<Δtkmin+Δtemin時，不滿足時間約束，一個周期內(nèi)無論離軌點在任何位置，飛行器都不可能到達著陸點，需在軌等待一個周期后重新進行判斷，直到滿足時間約束條件，方可進行后續(xù)計算。

3 能量約束建模及最短返回時間離軌點求解

飛行器由運行軌道進入過渡軌道、飛行到再入點，需要由制動火箭發(fā)動機提供制動速度增量。由于攜帶燃料有限，發(fā)動機所能提供的速度增量也是有限的。本節(jié)首先研究制動能量約束下滿足再入角要求的離軌點位置范圍，再綜合考慮等待段、離軌段及再入段時間，對離軌點位置進行優(yōu)化求解。

3.1 LVLH坐標(biāo)系下軌道方程

假設(shè)地球為均質(zhì)圓球，則制動發(fā)動機關(guān)機后飛行器將沿二體軌道運行。在LVLH坐標(biāo)系[17]中，將飛行器二體軌道運動方程進行投影，可得

(8)

式中:z為側(cè)向位移大小，vr，vβ和vz分別為速度矢量v在r軸、β軸和z軸上的分量，μ為地球引力常數(shù)。

假設(shè)在初始時刻，航天器位于軌道上的O點，當(dāng)前的運動狀態(tài)rk和vk可以由導(dǎo)航系統(tǒng)提供，再入點E的地心距re已知，則需要設(shè)計離軌軌道，使得飛行器到達再入點時，滿足再入角Θe和速度沖量Δvmax的要求。下一節(jié)將根據(jù)二體軌道邊值理論求解該問題。

3.2 基于再入角的離軌速度增量求解

軌道方程建立后，在給定再入角條件下，為快速確定能量約束下的離軌點位置范圍，需要推導(dǎo)離軌速度增量與離軌航程角的解析關(guān)系式。根據(jù)二體軌道理論，有

(9)

式中:f為真近點角，p為半通徑，e為偏心率。注意到Δf=Δβ，整理得

(10)

由上式可知(vrk+vre)和(vβk-vβe)成比例。故有

(11)

將式(11)代入式(10)，得到Θk與Θe的關(guān)系式

(12)

根據(jù)上式，就可以由rk，r，Θe和Δβk確定出在當(dāng)前點處的期望速度傾角

(13)

進而確定期望離軌軌道的半通徑及動量矩

(14)

(15)

式中:c為弦長，pm為最小能量橢圓的半通徑。進而可得到當(dāng)前時刻的增益速度和轉(zhuǎn)移軌道的半長軸

Δv=vR-vk

(16)

(17)

式中:vR為當(dāng)前位置的需要速度，其徑向和周向分量分別為

(18)

當(dāng)轉(zhuǎn)移軌道是橢圓時，轉(zhuǎn)移時間可由Lagrange方程

(19)

求得，其中，α和β表示Lagrange參數(shù)，滿足

(20)

此外，當(dāng)位置向量和速度向量均求得之后，離軌點和再入點的所有6個軌道要素均可以得到，故也可由開普勒方程求得轉(zhuǎn)移時間。

(21)

依據(jù)上述解析關(guān)系式，可以求得能量約束下離軌航程角及離軌時間的取值范圍。應(yīng)用序列二次規(guī)劃算法(SQP)、信賴域反射算法等都可以實現(xiàn)對該問題的求解，此處不做贅述。此外，可采用牛頓迭代公式

(22)

求解最短時間。為提高計算速度，f′(Δβk)可用Richardson外推法近似得到。設(shè)誤差階數(shù)為n，外推次數(shù)為m，則根據(jù)外推公式

(23)

可知，當(dāng)n和m足夠大時，D(n,m)可充分接近f′(Δβk)。通過求解，得到離軌航程角范圍Δβk∈[Δβkmin,Δβkmax]。

3.3 最短返回時間離軌點設(shè)計

在得到了滿足橫向位置約束、時間約束和能量約束的航程角范圍后，可將飛行器等待段航程角表示為：

(24)

式中:

(25)

根據(jù)等待段航程角及初始軌道根數(shù)，可求得離軌點軌道根數(shù)，得到離軌窗口。以返回段離軌時間為設(shè)計指標(biāo)，采用非線性優(yōu)化方法對離軌點位置進行尋優(yōu)，即可得到最短時間離軌點。

4 仿真校驗

本節(jié)將通過數(shù)值仿真驗證所提算法的有效性。設(shè)初始時刻t0為2022年3月29日9:00:00，飛行器位于近地圓軌道上，仿真參數(shù)見表2。

表2 仿真參數(shù)

參照升力體飛行器的性能指標(biāo)，給定再入段可達域。圖6給出了可達域的實際經(jīng)緯度范圍，由于考慮了地球自轉(zhuǎn)，可達域不是南北對稱的，最大橫程Lmax=2036.4km。將可達域縱程平均設(shè)置18個節(jié)點，橫程平均設(shè)置6個節(jié)點，每個節(jié)點處的再入段最小時間如表3所示。

圖6 仿真條件下可達域范圍

表3 各節(jié)點最小再入時間(單位：s)

根據(jù)假設(shè)，計算得到飛行器當(dāng)前點至著陸點的航程角Δβ=1.92rad。利用式(5)可求得目標(biāo)距星下點軌跡的最小橫向距離dmin=556.51km。由于dmin≤Lmax，故滿足位置約束。

由仿真條件，確定再入段時間取值范圍。當(dāng)橫程相同時，隨著縱程的增加，再入時間增加。當(dāng)橫程L=dmin時，計算得縱程取值范圍Rz∈[1083.2,4819.1]km。通過插值計算得到不同縱程對應(yīng)的再入時間和再入航程角，其中最小和最大再入時間分別為Δtemin=810.17s、Δtemax=3098.60s，對應(yīng)再入航程角Δβemin=0.17rad、Δβemax=0.7556rad。

下面計算離軌時間取值范圍。根據(jù)式，由離軌點及再入點高度、再入角計算得到不同離軌航程角Δβk對應(yīng)的離軌時間Δtk。圖7給出了二者的關(guān)系圖，可以看出二者呈正相關(guān)，標(biāo)注的紅色點表示速度沖量Δvr=Δvmax時航程角和離軌時間大小，當(dāng)曲線上點在該點下方時，滿足能量約束。設(shè)離軌航程角求解精度為10-6rad，選取滿足約束的點進行迭代，分別計算離軌段最短飛行時間Δtk，迭代結(jié)果如圖8所示。

圖7 離軌航程角與離軌時間關(guān)系圖

圖8 離軌段最短飛行時間迭代結(jié)果

可以看出，在合適初值條件下，兩種方法均求得Δtkmin=877.26s，對應(yīng)Δβkmin=1.02rad，結(jié)果較為準(zhǔn)確，過程較為平滑。但SQP算法迭代12次，用時0.705965s，而牛頓迭代法迭代9次，用時0.099555s，相比于SQP算法，牛頓迭代法迭代次數(shù)少、求解時間快，可以更好地實現(xiàn)最短離軌時間的求解。類似地，求得離軌段最大飛行時間Δtkmax=2069.9s，對應(yīng)的航程角Δβkmax=2.3284rad。

在求得離軌段及再入段航程角范圍后，可進一步求得等待段航程角Δβw∈[0,0.73]rad，進而確定離軌窗口。如圖9所示，標(biāo)注段表示離軌窗口。

圖9 離軌窗口

由式(6)和(25)計算出不同離軌和再入航程角下的等待時間及返回總時間，在所求離軌窗口的基礎(chǔ)上，應(yīng)用非線性優(yōu)化設(shè)計方法fmincon對返回總時間進行尋優(yōu)，結(jié)果如圖10所示，其中圖線①、②、③、④分別表示Δtw、Δtk、Δte和Δt的迭代過程?？梢钥闯?，在本例中，當(dāng)離軌段和再入段時間均取最小值時，返回段總時間最短，為2379.2s，離軌點軌道根數(shù)如表4所示。

表4 最短返回時間離軌點軌道根數(shù)

為進一步驗證算法的適用性，分別改變飛行器初始時刻位置和著陸點位置進行仿真。飛行器初始時刻位置改變，離軌窗口也相應(yīng)改變。圖11為一個回歸周期內(nèi)飛行器滿足約束的離軌點集合，其中標(biāo)粗部分為各圈離軌窗口?？梢钥闯觯谏κ斤w行器強橫向機動優(yōu)勢下，離軌窗口范圍較大，機動性和靈活性較強。

圖11 一個回歸周期內(nèi)飛行器離軌窗口

當(dāng)飛行器需應(yīng)急返回時，往往需要改變著陸點位置。假設(shè)除了主著陸場1外，還有7個副著陸場，表5給出了各著陸場位置及最短離軌時間仿真結(jié)果。可以看出，對于不同的著陸點位置，本文所提算法均有較好的適應(yīng)性。仿真結(jié)果表明，飛行器在8個著陸場均可正常著陸，其中著陸場7返回段飛行時間最短，因此可優(yōu)先選擇該著陸場進行應(yīng)急著陸。

表5 各著陸場位置及仿真結(jié)果

5 結(jié)論

為了實現(xiàn)升力式飛行器離軌點的快速設(shè)計，本文充分考慮其橫向機動能力強的特性，建立了基于再入可達域的時空約束模型，推導(dǎo)了燃耗約束下的離軌航程角解析計算方法，給出了確定最優(yōu)制動點位置的迭代算法。通過仿真分析可以得到以下結(jié)論：

1) 多約束下的最短時間離軌點確定策略能夠滿足升力式飛行器的離軌要求。經(jīng)過較少次數(shù)迭代后，離軌段航程角誤差收斂至較高精度,速度較快，結(jié)果較為準(zhǔn)確；

2) 對飛行器不同飛行時刻的離軌窗口求解表明，本算法對于不同起始位置均有較好的適應(yīng)性，升力式飛行器的橫向機動能力拓寬了其離軌窗口，增加了返回的機動性和靈活性；

3) 對不同位置著陸點的最優(yōu)離軌點求解表明，當(dāng)飛行器需要應(yīng)急返回時，該算法能夠基于當(dāng)前時刻狀態(tài)進行快速計算，通過對比確定最優(yōu)著陸點，具有較強的可行性和實時性。