王澤國，孟 斌

(1.北京控制工程研究所，北京100190;2.空間智能控制技術(shù)重點實驗室，北京100190)

再入動力學(xué)的性質(zhì)及其在軌跡優(yōu)化中的應(yīng)用*

王澤國1，2，孟 斌1，2

(1.北京控制工程研究所，北京100190;2.空間智能控制技術(shù)重點實驗室，北京100190)

考慮具有終端約束和過程約束的探月返回飛行器再入軌跡設(shè)計問題，通過將性能指標(biāo)泛函定義為再入終端位置誤差的平方和，再入軌跡設(shè)計問題轉(zhuǎn)化為具有過程約束和狀態(tài)方程約束的優(yōu)化問題.首先僅考慮狀態(tài)方程約束，利用最大值原理，得到該優(yōu)化問題的必要條件，選取間接法中的共軛梯度算法求解最優(yōu)控制量.進(jìn)而針對軌跡約束問題，研究了再入過載和軌道飛行段飛行距離與航跡角以及傾側(cè)角的關(guān)系，在此基礎(chǔ)上，提出了采用調(diào)整初始傾側(cè)角序列的方法實現(xiàn)過程約束.該算法克服了罰函數(shù)方法中需要調(diào)節(jié)參數(shù)較多的問題，并且物理意義明確，實現(xiàn)簡單.最后，給出了Apollo再入軌跡優(yōu)化的數(shù)值仿真算例，驗證了所給出算法的有效性.

探月返回;跳躍式再入;軌跡優(yōu)化;再入動力學(xué)的性質(zhì);初值調(diào)整

近年來，探月飛行器的研究繼Apollo后得到廣 泛關(guān)注.對于探月返回再入軌跡規(guī)劃問題，由于飛行器運(yùn)動方程復(fù)雜，一般采用數(shù)值方法求解.數(shù)值方法可以分為直接法和間接法，兩種方法各有優(yōu)缺點[1].直接法將優(yōu)化問題轉(zhuǎn)化為非線性規(guī)劃問題，采用某種優(yōu)化算法來求解.由于飛行器運(yùn)動精度要求較高，因此需要更多的離散點個數(shù)，并且直接法的初始狀態(tài)難以選擇，這給軌跡設(shè)計帶來不便.本文考慮間接求解方法.在探月返回再入軌跡規(guī)劃的間接法研究中，Istratie等[2]給出了在給定初始條件和最終條件下的最優(yōu)軌跡，但未考慮過載和熱流的限制.南英等[3]主要以軌跡中的過載、熱流和動壓的目標(biāo)給出了最優(yōu)軌跡，并比較不同次數(shù)再入情況下各個目標(biāo)，但尚未考慮終端約束的限制.王樂等[4]給出了采用共軛梯度法設(shè)計臨近空間飛行器再入軌跡的一般方法，給出了終端時刻固定的設(shè)計方法.在間接法中，當(dāng)考慮過載等過程約束問題時，一般采用罰函數(shù)方法[3-4]，這將導(dǎo)致需要調(diào)節(jié)的參數(shù)增加，并且增加了算法收斂的難度.

本文研究探月返回飛行器有約束的跳躍式再入軌跡規(guī)劃問題，軌跡約束包括終端約束、過載約束和熱流約束.通過將性能指標(biāo)泛函定義為再入終端位置誤差的平方和，再入軌跡設(shè)計問題轉(zhuǎn)化為具有過程約束和狀態(tài)方程約束的優(yōu)化問題.首先僅考慮狀態(tài)方程約束，利用最大值原理，得到該優(yōu)化問題的必要條件，選取間接法中的共軛梯度算法求解最優(yōu)控制量.針對過載約束問題，通過分析過載、軌道飛行階段飛行距離和傾側(cè)角的關(guān)系，給出了初始控制量調(diào)節(jié)方法以實現(xiàn)過載約束.該算法克服了罰函數(shù)方法中需要調(diào)節(jié)參數(shù)較多的問題，并且物理意義明確，實現(xiàn)簡單.最后，給出了Apollo再入軌跡優(yōu)化的數(shù)值仿真算例，驗證了所給出算法的有效性.

1 問題描述

忽略地球轉(zhuǎn)動，探月返回飛行器無量綱運(yùn)動方程為

對于具有低升阻比特性的再入飛船，若要求其在執(zhí)行探月返回再入任務(wù)時實現(xiàn)長的飛行縱程，唯一可行的辦法是采用跳躍式再入飛行.圖 1為典型的跳躍式軌跡示意圖.為了論述清晰，將跳躍式軌跡分為5段，如圖 1所示.一般把高度120km作為再入起始點，再入初期空氣阻力幾乎為零，控制量對于飛行器的影響很弱，從120km到=0的再入初始階段作為再入的第1階段=0到=0段為第2階段.假設(shè)飛行器高度第2次到達(dá)80km以上時，忽略大氣的影響，飛行器以橢圓軌道運(yùn)行.把飛行器第2次到達(dá)80km時稱為躍出點，從=0到躍出點為第3階段.把橢圓軌道飛行階段稱為第4階段.第5階段為飛行器從80km第2次再入大氣層的飛行段.分別記每個階段的起始時刻為ti，其中i表示階段數(shù)，如圖1所示.本文僅考慮在前3階段存在一次˙=0和=0的情形，但所給出的結(jié)果可以推廣到一般情形.

本文考慮探月返回飛行器受約束的跳躍式軌跡規(guī)劃問題，軌跡約束包括過程約束、終端約束和控制量約束，過程約束包括過載和熱流約束，終端約束為給定的飛行器的狀態(tài)終值.具體地說，本文研究在給定初始狀態(tài)(r0、θ0、φ0、V0、γ0和 ψ0)下，設(shè)計傾側(cè)角變化曲線，使得飛行器再入過程滿足狀態(tài)方程(1)、過載限制、熱流限制和傾側(cè)角約束，并達(dá)到給定終端狀態(tài)(Vf、θsite和φsite)，即滿足終端約束.針對上述問題，本文研究了共軛梯度算法中軌跡約束過載的實現(xiàn)問題，提出了初始傾側(cè)角序列調(diào)整的方法實現(xiàn)過程約束.下面首先給出傾側(cè)角與軌跡約束過載的關(guān)系，在此基礎(chǔ)上，給出探月返回飛行器再入軌跡優(yōu)化設(shè)計方法.

圖1 跳躍式軌跡Fig.1 Skip trajectory

2 過載和第4階段飛行距離與傾側(cè)角的關(guān)系

本節(jié)考慮軌跡約束過載和第4階段飛行距離與傾側(cè)角的關(guān)系，首先考慮過載問題.由于第二次是以第一宇宙速度再入，與神舟飛船再入問題相同，因此這里主要考慮前3個階段的過載問題.

從方程(1)可以看到，在縱向平面控制量傾側(cè)角σ只是直接影響航跡角γ.因此，本節(jié)通過分析過載和第4階段飛行距離與航跡角 γ的關(guān)系，以及控制量σ與航跡角γ的關(guān)系，給出過載和第4階段飛行距離與控制量σ的關(guān)系.從而可以據(jù)此調(diào)節(jié)最優(yōu)控制算法中的控制量序列.

2.1 過載和航跡角的關(guān)系

本小節(jié)分析第2、3階段過載和航跡角的關(guān)系.過載為

其中動壓q=ρV2/2.易見，由于小，動壓q是過載n的主要影響因素.下面通過對動壓q和航跡角γ關(guān)系的分析找到過載n和航跡角γ的關(guān)系.

引理1.針對系統(tǒng)(1)，˙q=0當(dāng)且僅當(dāng)sinγ=γth，其中，γth=2/(βV2).并且如果sinγ＞γth，那么＜0;反之亦然.

證明.為了表達(dá)簡單，取q中的變量，定義

注意到 V＞0，β＞0，h＞0，以及 γth的定義，引理得證.

注意到閾值 γth是變量，下面證明在第2階段，至少存在一個時刻，使得sinγ=γth.

引理2.針對系統(tǒng)(1)，關(guān)于γth的一些性質(zhì):

1)在第2、3階段，γth＜0.

2)假設(shè)CD為常數(shù)，航跡角γ幅值較小.則在第2階段，γth減小;第3階段γth增大.

3)假設(shè)CD為常數(shù)，航跡角γ較小.那么在第2階段滿足如下關(guān)系:如果sinγ(t2)＜γth(t2)，則在第2階段，存在一個時刻 tm，t2＜tm＜t3，使得sinγ(tm)=γth(tm)，其中t2和t3分別表示第2階段的初始時刻和結(jié)束時刻;進(jìn)而如果在第2階段有 ˙γ＞0，那么sin(γ)與γth有唯一的交點.

4)假設(shè)CD為常數(shù)，航跡角γ較小.則在第3階段sinγ＞γth.

2)當(dāng)航跡角γ幅值較小時，由式(1)和Vs的定義可得，

由上式可知，如果CD為常數(shù)，那么當(dāng) ˙h＞0時，γth增大;反之亦然.亦即在第2階段，γth減小;第3階段γth增大.

3)由定義可知，第2階段末端t3時刻的航跡角為0，由 1)可知，γth＜0，所以如果 sinγ(t2) ＜γth(t2)，則由各狀態(tài)變量的連續(xù)性可知，在第二階段存在一個時刻tm，t2＜tm＜t3，使得sinγ(tm)= γth(tm).進(jìn)而如果第二階段 ˙γ＞0，結(jié)合 γth單調(diào)遞減，則在第2階段sinγ與γth僅有一個交點.

4)由第3階段γ＞0＞γth，顯然得證.

由上述引理，可以得到第2、3階段動壓的性質(zhì).引理3.針對系統(tǒng)(1)，假設(shè)CD為常數(shù)，航跡角γ較小.如果sinγ(t2)＜γth(t2)，那么第2、3階段的動壓在第2階段達(dá)到最大.進(jìn)而如果第2階段有 ˙γ＞0，那么第2階段動壓有唯一的極值點.證明.假設(shè)CD為常數(shù)，航跡角γ較小.由引理1和引理2，第3階段γ＞0＞γth，則˙q＜0，因此，q(t)＜q(t3)，t3＜t＜t4，即，動壓在第2階段達(dá)到極大.進(jìn)而如果第2階段有 ˙γ＞0，由引理1和引理2，sinγ與γth有唯一的交點.也即第2階段動壓有唯一極大值點.

由引理3可知，再入過程第2、3階段的過載最大值出現(xiàn)在第2階段，我們可以將該問題歸納為極小極大問題，給出過載最大值極小的嚴(yán)格的理論結(jié)果.本文僅對該問題進(jìn)行簡單分析.為了減少過載最大值，需要增加或者減小航跡角，但是這可能導(dǎo)致增加高度的同時增加速度，或者減小高度的同時減小速度，因此對動壓產(chǎn)生復(fù)雜的影響.

動壓與高度和速度均相關(guān)，由式(1)可知，航跡角對于速度和高度具有相反的作用，因此對于動壓的影響較復(fù)雜.下面通過比較3組具有不同高度和速度的狀態(tài)(hm，Vm)、(hm(1+a)，Vm)和 (hm，Vm(1-b))的動壓值，給出航跡角對于動壓的影響關(guān)系.實際中軌跡的高度與速度的變化不超過50%，即a∈(0，0.5)，b∈(0，0.5).分別記三組狀態(tài)所對應(yīng)的動壓分別為qm、qm1、qm2.則

假設(shè)動壓最大值出現(xiàn)在高度為60km左右，則e-βhm約為6.85×10-4.如果qm1=qm2，則有b=1-0.0262a.當(dāng)a∈(0，0.5)時，b＞a，即如果動壓改變相同量，則速度減小幅度要大于高度增加幅度，因此高度增加對減小動壓更加有效.所以在第2階段通過增大航跡角，增加高度，可以有效減小動壓極值.

2.2 第4階段飛行距離和航跡角的關(guān)系

按照第1節(jié)的定義，把橢圓軌道飛行階段稱為第4階段，如圖1所示.一般來說第4階段飛行器處于不受控的狀態(tài)，本部分討論該段飛行距離和航跡角的關(guān)系，從而給出控制量傾側(cè)角對于第4階段飛行距離的影響關(guān)系.

引理4.假設(shè)飛行器高度第2次到達(dá)80km以上時，忽略大氣的影響，飛行器以橢圓軌道運(yùn)行，躍出點速常數(shù).則第4階段飛行距離與躍出點航跡角正相關(guān).證明.由假設(shè)，飛行器高度第二次到達(dá)80km以上時，忽略大氣的影響，飛行器以橢圓軌道運(yùn)行.下面由躍出點真近點角與躍出點航跡角的關(guān)系，得出第4階段飛行距離與航跡角的關(guān)系.由文獻(xiàn)[6]可知，躍出點飛行器的真近點角?t(如圖2所示)滿足

下面證明

由式(2)，求得cos?t關(guān)于航跡角γt的導(dǎo)數(shù)

圖2 跳躍式再入軌跡Fig.2 Skip reentry trajectory

由假設(shè)、式(1)及γt較小，忽略阻力D和sinγt，得

由于?t＞90°，所以cos?t＜0.由于Ct＜1＜1/cosγt＜1/cos2γt，可得

和

因為γt為大于0的較小值，我們有

結(jié)合上面的結(jié)果，可得

即cos?t與γt正相關(guān)，?t與γt負(fù)相關(guān).由于第4階段飛行距離與?t負(fù)相關(guān)，所以第4階段飛行距離與γt正相關(guān).

由上述引理可知，第4階段飛行距離與航跡角γt正相關(guān).而第4階段飛行距離一般是確定的，也就是躍出點航跡角γt不能有較大變化.

2.3 過程約束與傾側(cè)角的關(guān)系

2.1 和2.2中分析了過載和第4階段飛行距離和航跡角的關(guān)系，本部分首先分析傾側(cè)角和航跡角的關(guān)系，進(jìn)而給出過程約束與傾側(cè)角的關(guān)系.

由cosγ≈1、1/r≈1，和式(1)可得

由于傾側(cè)角可以在-180°至180°之間變化，cosσ∈[-1，1]，若L＞V2-1，那么傾側(cè)角的選擇會影響航跡角的導(dǎo)數(shù)符號，我們說傾側(cè)角對于航跡角有主要影響.

引理5.假設(shè)升力系數(shù)CL是常數(shù).則當(dāng)V＞1，h＜導(dǎo)數(shù)的符號可以由傾側(cè)角的幅值控制.

時，L＞V2-1.由式(3)可知，當(dāng)L＞V2-1時，航跡角的導(dǎo)數(shù)的符號可以由傾側(cè)角的幅值控制.由于ln(V2)-ln(V2-1)＞0，可知若V＞1，h＜

通過將常數(shù)CL=0.4717及仿真中表1的常數(shù)可估算出當(dāng) V＞1，高度小于68.92km時傾側(cè)角的調(diào)整對于航跡角起主要作用.實際上所忽略的[ln(V2)-ln(V2-1)]/β對于高度估計也有影響，所以在更高的高度就可以使得控制量起主要作用.而且由于在較低高度的阻力很大，使得速度迅速下降，這樣使得以后的控制量一直占主要影響.

綜合上述分析，可以考慮減小第1階段和第2階段的傾側(cè)角，并增加第3階段的傾側(cè)角，以減小動壓極大值并保證軌道飛行階段的飛行距離.

軌跡約束中關(guān)于熱流的限制問題也是需要注意的.但由于熱流最大值一般出現(xiàn)在再入初始階段，因此控制量的選取對于熱流的影響較小[7]，所以不對熱流限制做單獨設(shè)計.

3 最優(yōu)軌跡規(guī)劃

本節(jié)考慮探月返回飛行器受約束的跳躍式軌跡規(guī)劃問題.通過將性能指標(biāo)泛函定義為再入終端位置誤差的平方和，再入軌跡設(shè)計問題轉(zhuǎn)化為具有過程約束和狀態(tài)方程約束的優(yōu)化問題.首先僅考慮狀態(tài)方程約束，利用最大值原理，得到該優(yōu)化問題的必要條件，選取間接法中的共軛梯度算法求解最優(yōu)控制量，利用上節(jié)給出的傾側(cè)角和過程約束的關(guān)系，通過調(diào)整初始控制量序列以滿足過程約束.

選擇目標(biāo)泛函為飛行器的終端經(jīng)度 θf和緯度φf與目標(biāo)點(θsite，φsite)的差的平方和J=k1(θfθsite)2+k2(φf-φsite)2最小[4]，其k1和k2為待定常數(shù).則Hamiltonian函數(shù)為

其中λr，λθ，λφ，λV，λγ，λψ為相應(yīng)的協(xié)態(tài)變量.根據(jù)最大值原理，可得協(xié)態(tài)變量滿足:

協(xié)態(tài)變量的終值為

其他終值為零;最優(yōu)控制量σ滿足

和控制輸入約束.

下面采用共軛梯度法，求解上述無過程約束問題.共軛梯度法的主要思想是先給出一組初始控制量序列，在每次迭代中計算Hamiltonian函數(shù)關(guān)于控制量的偏導(dǎo)，控制量沿著令Hamiltonian函數(shù)減小的方向，最后達(dá)到使得目標(biāo)泛函最小的最優(yōu)解[8-9].

上述算法中需要調(diào)節(jié)的主要是初始控制量σ0.初始控制量序列σ0的選擇根據(jù)第2節(jié)的結(jié)果進(jìn)行調(diào)節(jié)，首先選取某個常數(shù)序列使得飛行器的縱程誤差較小，并且算法收斂，然后通過翻轉(zhuǎn)實現(xiàn)橫程調(diào)節(jié)，進(jìn)一步通過調(diào)整第2階段的傾側(cè)角降低過載，增加精度.

4 數(shù)值仿真

本節(jié)針對探月返回飛行器Apollo，利用初值選取原則和共軛梯度法設(shè)計受限再入軌跡.Apollo的參數(shù)和給定的初始狀態(tài)和終端狀態(tài)見表 1.飛行器在運(yùn)行中需要考慮過載、熱流和控制輸入的限制.在文獻(xiàn)[7]中給出了幾類過載限制，分別是正常乘員逃生時的最大過載，正常乘員再入時的最大過載，和受傷乘員承受的最大過載.這里采用其中最嚴(yán)格的限制，即受傷乘員承受的最大過載.仿真中考慮了熱流的限制[10](見表 1)，控制量的限制為 -180°≤σ≤180°.

下面按照第3節(jié)的算法給出再入軌跡，通過初始控制量序列調(diào)節(jié)方法以滿足過載約束.首先選擇一常數(shù)控制量序列，然后通過嘗試得到飛行器翻轉(zhuǎn)時刻，具體設(shè)計為常值傾側(cè)角-45°，在時刻85s翻轉(zhuǎn).通過優(yōu)化，飛行器剛進(jìn)入大氣層時過載比較大，在10s附近與乘員能夠承受的最大過載重合.

下面通過調(diào)整初始控制量的方式來減少過載.由2.3節(jié)的結(jié)論，選擇再入初始時較小的傾側(cè)角幅值和躍出點較大的傾側(cè)角幅值，第2次再入后選擇一個飛行器翻轉(zhuǎn)時刻以達(dá)到橫向控制.采用這種控制策略，經(jīng)過試驗調(diào)節(jié)，可以得到過載和落點誤差更小的軌跡.圖 3為修改后的初始傾側(cè)角序列，再入初始時傾側(cè)角為零，然后迅速上升至180°以防止飛行器橢圓軌道飛行距離過大，接著按 9°飛行，在700s翻轉(zhuǎn).圖4是優(yōu)化后的傾側(cè)角變化曲線，得到仿真結(jié)果如圖 5～7所示.由圖 5可見過載明顯降低，最大過載由9.8832g降至7.2557g，滿足過載約束.由圖 6可見開傘點與期望點幾乎重合，誤差為8.2050km，與調(diào)整前相比明顯減少.從圖7可見算法收斂較快.熱流的最大值為1.7242×105W/m2，低于表1中要求的熱流最大值.

表1 常數(shù)Tab.1 Constants

圖3 調(diào)整后初始控制量序列Fig.3 Initial control profile after ad justment

圖4 優(yōu)化后控制量序列Fig.4 Optimized control profile

圖5 過載和限制Fig.5 Overload and limit

圖6 飛行器的經(jīng)緯度變化Fig.6 Latitude and longitude of the vehicle

圖7 性能指標(biāo)變化曲線Fig.7 Performance index curve

5 結(jié) 論

本文考慮探月返回飛行器受限軌跡規(guī)劃問題.研究了過載和軌道飛行段飛行距離與傾側(cè)角的關(guān)系，給出了過載極值與航跡角的關(guān)系和軌道飛行段飛行距離與躍出點航跡角的關(guān)系，以及控制量傾側(cè)角與過程限制的關(guān)系.在此基礎(chǔ)上，給出了共軛梯度算法過程限制的初始控制量序列調(diào)節(jié)方法.最后通過數(shù)值仿真給出了滿足約束的Apollo再入最優(yōu)軌跡，表明所給出算法的有效性.

[1] Betts J.Survey of numerical methods for trajectory optimization[J].Journal of Guidance，Control，and Dynamics，1998，21(2):193-207

[2] Istratie V.Optimal skip entry into atmosphere with minimum heat and constraints[C].AIAA Atmospheri Flight Mehanis Conferene，Denver，2000

[3] 南英，陸宇平，龔平.登月返回地球再入軌跡的優(yōu)化設(shè)計[J].宇航學(xué)報，2009，30(5):1842-1847 Nan Y，Lu Y P，Gong P.Optimal reentry trajectory design for mooncraft returning to the Earth[J].Journal of Astronautics，2009，30(5):1842-1847

[4] 王樂，方群，歐岳峰.臨近空間飛行器再入軌跡優(yōu)化設(shè)計[J].科學(xué)技術(shù)與工程，2011，11(3):543-547 Wang L，F(xiàn)ang Q，Ou Y F.The optimization design of the near space vehicles're-entry trajectory[J].Science Technology and Engineering，2011，11(3):543-547

[5] U.S.national oceanic and atmospheric administration，US standard atmosphere[R].Washington，D.C.:U.S.National Oceanic and Atmosphemic Administration，1976

[6] 耿長福.航天器動力學(xué)[M].北京:中國科學(xué)技術(shù)出版社，2006:273-275 Geng C F.Spacecraft dynamics[M].Beijing:China Science and Technology Press，2006:273-275

[7] Bairstow S，Barton G.Orion reentry guidance with extended range capability using PredGuid[C].AIAA Guidance，Navigation and Control Conference and Exhibit，Hilton Head，2007

[8] Bryson A E，Ho Y C.Applied optimal control[M]. Waltham，Mass.:Blaisdell，1969

[9] 宮錫芳.最優(yōu)控制問題的計算方法[M].北京:科學(xué)出版社，1969:97-104 Gong X F.Computation methods on optimal control problems[M].Beijing:Science Press，1969:97-104

[10] Xue S，Lu P.Constrained predictor-corrector entry guidance[J].AIAA Guidance，Control，and Dynamics，2010，33(4):1273-1281

Properties of Reentry Dynam ics and Their App lication in Trajectory Op tim ization

WANG Zeguo1，2，MENG Bin1，2
(1.Beijing Institute of Control Engineering，Beijing 100190，China; 2.Science and Technology on Space Intelligent Control Laboratory，Beijing 100190，China)

The reentry trajectory optim ization problem with term inal constraints and path constraints of a lunar return vehicle is considered in this paper.By defining the performance index function as a square sum of the reentry term inal position errors，the reentry trajectory design problem is transformed into the optimization problem with terminal constraints and state equation constraints.First，in the presence of the state equation constraints only，the necessary conditions of the optim ization problem are obtained by using the maximum principle，and the optimal control is solved using the conjugate gradientmethod.Then，in the presence of path constraints，we have studied the relationship among the overload，the distance of the orbit phase，the flight path angle and the bank angle.Based on the above，amethod of initial bank angle value ad justment is presented to satisfy the path constraints.This approach，with definite physicalmeaning and simple implementation，circumvents the shortcoming of the punishment-functionmethod which involves too many ad justed parameters.Finally，a numerical example for Apollo reentry trajectory optimization is given to illustrate the effectiveness of the algorithm.

lunar return;skip reentry;trajectory optim ization;property of reentry dynam ics;initial value adjustment

V448.235

A

1674-1579(2012)06-0006-07

王澤國(1988—)，男，碩士研究生，研究方向為航天器再入制導(dǎo);孟 斌(1973—)，女，高級工程師，研究方向為航天器控制.

*國家自然科學(xué)基金資助項目(61273153，60736023，60704014).

2012-08-03

DO I:10.3969/j.issn.1674-1579.2012.06.002