基于速度增益制導的大橢圓轉移中段制導方法研究

任金磊，李 君，莊學彬，李羅鋼，謝澤兵

（中國運載火箭技術研究院研究發展中心，北京，100076）

基于速度增益制導的大橢圓轉移中段制導方法研究

任金磊，李 君，莊學彬，李羅鋼，謝澤兵

（中國運載火箭技術研究院研究發展中心，北京，100076）

針對大橢圓轉移中初始偏差和長時間飛行中攝動因素導致的方法偏差較大的現象，需進行中段軌道修正。基于顯式制導的思路，采用速度增益制導，研究了考慮項攝動有限推力條件下的中段軌道修正方法，并通過優化修正時間設計了終端偏差小、燃料消耗少的中段修正優化軌道。針對不同初始偏差工況通過數學建模仿真進行了驗證，結果表明該方法可行。

速度增益；顯式制導；中段修正；軌道優化

0 引 言

對于大橢圓轉移軌道、大范圍遠程交會等飛行時間較長的飛行任務，由于推力大小誤差、開關機時間誤差、推力方向誤差等多種誤差會導致一定的初始誤差，這些誤差經過長時間飛行往往造成很大的終端偏差，導致定點入軌等多種飛行任務失敗，因此有必要進行中段軌道修正。如土星探測器卡西尼[1]、火星探測器勇氣號和機遇號為完成飛行任務[2，3]，進行了多次中段修正；中國嫦娥三號的地-月轉移也采用了軌道修正技術[4]。

中制導方法可以分為攝動制導與顯式制導。攝動制導是一種基于標稱軌道的制導方法，簡單方便，但缺乏理論依據。顯式制導是根據目標數據與飛行器當時的運動參數，實時計算并重新設計飛行軌跡的方法，自適應性好，魯棒性強。文獻[5]提供了基于狀態方程線性化的攝動制導方法，對項攝動進行了補償修正。但該方法對標稱軌道依賴性較高，沒有進行優化設計；文獻[6]利用顯式制導法，采用Breakwell間距比策略，給出了基于分段落點預報顯式制導的月-地返回軌道中途修正方案。本文基于顯式制導思想，根據大橢圓轉移軌道終端入軌精度要求，設計了以終端偏差小、燃料消耗少為目標的中段修正優化軌道，制導方法采用速度增益制導，并給出了仿真算例。

1 軌道動力學模型

當考慮地球的中心引力和其他攝動力的作用時，受攝二體問題的運動方程如下：

由于實際地球并非球對稱，具有扁平度、梨形和赤道變形等特性，其產生微小的非球形引力對飛行器軌道會產生攝動作用。由于J2項攝動的量級遠大于其他攝動，僅考慮地球扁率項（J2項）的影響就能滿足要求。因此，設計飛行器的機動軌道時，軌道動力學方程可以簡化為

式中 g為地球引力加速度矢量；aJ2為考慮地球扁率項（J2項）影響的攝動加速度矢量。

2 軌道修正策略基本思想與制導方法

2.1 軌道修正策略基本思想

采用終段控制法對大橢圓轉移軌道進行中段修正，其基本思想為：以理想終端狀態點為控制目標，根據當前機動時刻的實際狀態與理想的終端狀態點，通過蘭伯特（Lambert）理論，重新設計一條新的軌跡來滿足飛行要求。圖1為終端控制修正方法在轉移軌道中的示意圖。

圖1 終段控制法

如圖1所示，根據理想的入軌位置A和目標位置p0，以及轉移時間T,采用Lambert方法，設計1條二體標稱軌道，該軌道為理想軌道。實際飛行中，飛行器入軌時會存在一定的初始誤差，入軌點A1往往會偏離理想入軌位置A；其次飛行中受攝動影響，也會偏離標稱軌道，而且飛行時間越長，終端偏差越大。由初始誤差和攝動影響會導致實際位置p1與目標位置p0存在較大偏差，這種偏差往往無法滿足末端入軌精度要求，會導致末端入軌精度要求較高的飛行任務失敗。此時，通過優化算法以最終位置誤差和修正燃料消耗相關的綜合指標最優為目的，在途中選擇C點進行修正，根據修正點C和目標位置p0，采用Lambert方法，設計1條修正軌道。受攝動等誤差影響，終端入軌點p2仍與目標位置p0存在一定的位置偏差，但是相對于初始入軌點已有p1很大改善。在中途修正時采用速度增益顯式制導進行制導控制。

2.2 速度增益制導方法

速度增益制導是一種閉路制導。這種閉路制導根據需要速度對飛行器進行導引和控制，以使其在tf時刻通過預定入軌點。

設t時刻，飛行器的位置矢量為r，由任務要求可得此刻需要的速度vr，t時刻的實際速度為v已知量，則求得兩者之差：

式中 vg為速度增益。

對式（3）求導，并代入式（2）可得：

式中 a1=g+aJ2；可將dvr/dt分為兩部分：

定義b=dvr/dt-a1=δvr/dt，將式（5）代入式（4）可得：

確定aT的原則是要求dvg/dt能有效地消除vg，這樣aT應滿足如下必要條件：

此時可以通過控制變量γ（γ取值范圍為（0，1）），滿足(γb-aT)與的方向相反。其中，eg為vg的單位矢量。則，aT應滿足：

式（8）中的制導方法即速度增益制導。

3 中段修正優化軌道方法

遺傳算法在求解多目標優化問題時，具有廣泛的應用價值。對于本文研究的問題而言，優化的目標是最終位置誤差和修正燃料消耗相關的綜合指標，優化變量是修正的時間點[7，8]。假設，飛行器最終與目標點的位置誤差為Δd，修正所消耗的速度脈沖為Δv，相應的燃料消耗為Δm，修正的時間點為t1，則根據pareto最優解原理，設終端偏差的權重系數為k1，修正燃料消耗的權重系數為k2，優化目標函數和約束條件為

式中 J(X)為目標函數，X為優化變量，為進入轉移軌道后到修正點所飛行的時間；k1和k2分別為終端偏差和燃料消耗的加權系數，且有k1,k2∈[01],k1+k2=1；ΔSi為地球轉移軌道（GTO）最終目標點的軌道約束，包含6個軌道根數（α為軌道半長軸，e為偏心率，i為軌道傾角，Ω為升交點赤經，ω為近地點幅角，f為真近點角），能保證軌道轉移進入固定目標點。此外，考慮飛行器其它約束條件，飛行器修正時間不得大于12 000 s不得小于3000 s。

當k1=1，k2=2時，式（9）優化問題為終端偏差最優；當k1=0，k2=1時，優化問題為燃料消耗最優；當k1≠0，k2≠0時，優化問題為終端偏差和燃料消耗綜合最優。以轉移軌道飛行中存在初始偏差為例，沿J2000坐標系3個軸方向加入初始誤差，位置誤差取10 km，速度誤差取5 m/s，利用遺傳算法對上述問題進行尋優。不同性能指標的優化結果如表1所示。

表1 不同性能指標的優化結果

通過表1可發現，在燃料消耗和終端偏差加權系數都為0.5時，優化綜合性能較好，既可以達到比較小的位置偏差，燃料消耗也較小。在下文的仿真中本文主要以k1=0.5，k2=0.5的優化函數進行優化。

4 仿真結果

航天器由停泊軌道進入GTO，飛行轉移時間為20 292 s，外推至預定的變軌終端點。此時，標稱軌道的變軌終端點在J2000坐標系下的位置、速度參數列于表2。

表2 標稱軌道GTO高軌入軌點位置速度

在仿真中，沿J2000坐標系3個軸方向加入初始誤差，位置誤差取10 km，速度誤差取5 m/s，不進行修正，變軌終端點存在較大偏差，其偏差如表3所示。

表3 未加修正GTO高軌入軌點位置速度偏差

采用速度增益終端控制法對考慮J2項攝動影響的軌道進行修正，修正推力采用有限推力模型。仿真中，取飛行器進入轉移軌道后的質量為814 kg，發動機比沖為2600 m/s，推力為3000 N。根據齊奧爾科夫斯基公式估算發動機工作時間Δts為32.47 s，然后在變軌點前Δts/2時刻開機。

采用遺傳優化算法通過控制中段修正時間以達到終端偏差和燃料消耗最小，優化后修正點選取在進入轉移軌道后9019 s。選取參數γ=1，對應的中制導過程中需要速度、實際速度及增益速度仿真結果，如圖2～4所示。

由中段制導發動機關機按時間t1外推至預定的變軌終端點其相對位置曲線如圖5所示。

速度增益終端控制中段修正后，飛行器的燃料消耗為21.26 kg，關機點處的增益速度偏差達到0.057 m/s，而以此關機點將GTO積分至變軌終端點時，其位置偏差為6.84 m。

圖2 中段修正需要速度變化曲線

圖3 中段修正實際速度變化曲線

圖4 中段修正增益速度變化曲線

圖5 中段修正后自由飛行段相對位置變化曲線

速度增益終端控制中段修正后，飛行器的燃料消耗為21.26 kg，關機點處的增益速度偏差達到0.057 m/s，而以此關機點將GTO積分至變軌終端點時，其位置偏差為6.84 m。

為進一步驗證該方法的有效性，分別對位置偏差取為5 km，10 km，20 km，速度偏差取為5 m/s，10 m/s，6種不同工況進行仿真，并對最終入軌點偏差（相對標稱軌道終點）、中制導燃料消耗、關機點增益速度,優化修正時間進行了對比。其結果列于表4。

表4 不同工況的仿真結果對比

由表4可以看出，采用速度增益終端控制法進行修正，終端位置偏差可以控制在比較小的范圍，所需的燃料也較少。隨著初始誤差的增大，燃料消耗會相對略有增大。相對于位置偏差，速度偏差比較大，所以可以看出速度增益制導是一種攔截式制導律。另外對于GTO進入GEO需要速度增量在1800 m/s以上來說，修正后的終端速度偏差相對較小，可以在末制導做進一步修正。通過修正時間優化結果可以看出，修正時間不宜過早也不宜過晚，根據不同的初始偏差，在進行適當調整。

5 結束語

本文以大橢圓轉移軌道等長時間中段飛行的航天任務為研究背景，研究了基于速度增益制導的中段修正方法，并通過遺傳優化算法設計了以修正時間點為約束條件，以終端偏差和燃料消耗為優化目標的中段修正優化軌道。在考慮項攝動有限推力的影響下，該方法可以有效地提高終端位置精度。同時本文對不同工況進行了對比分析，進一步驗證了該制導控制方法的可行性。

[1] Burk T. Attitude control performance during Cassini trajectory correction maneuvers[C]. San Francisco: AIAA Guidance，Navigation，and Control Conference and Exhibit, 2005.

[2] D’Amario L A. Mars exploration rovers navigation results[C]. Rhode Island: AIAA/AAS Atrodynamics Specialist Conference and Exhibit, Providence, 2004.

[3] Potts C L, Raofi B, Kangas J A. Mars exploration rovers propulsive maneuver design[C]. Rhode Island: AIAA/AAS Astrodynamics Specialist Conference and Exhibit, Providence, 2004.

[4] Kallemeyn P H, Knocke P C, Burkhart P D, et al. Navigationand guidance for the Mars Surveyor '98 mission[R]. AIAA Paper 98-4562, 1998.

[5] 田野, 張洪波, 吳杰. 基于軌道方程線性化的中段飛行軌道修正方法.飛行力學[J]. 2011, 29(3): 68-71.

[6] 胡軍, 周亮. 基于顯示制導的月地返回軌道中途修正研究[J]. 中國空間科學技術. 2010, 5: 7-13.

[7] 盧山, 徐世杰. 基于遺傳算法的轉移軌道中途修正技術研究[C].全國第十二屆空間及運動體控制技術學術會議論文集, 2006.

[8] 盧山, 段佳佳, 徐世杰. 基于蟻群算法的轉移軌道中途修正問題研究[J].系統仿真學報, 2009, 21(14): 4400-4404.

Research on Midcourse Guidance Method of Large Elliptic Transfer Orbits Based on Velocity Gain Gudance

Ren Jin-lei, Li Jun, Zhuang Xue-bin, Li Luo-gang, Xie Ze-bing
(R &D Center, China Academy of Launch Vehicle Technology, Beijing, 100076)

For the phenomenon of large deviation caused by initial deviation and a long-time flight perturbation factors, mid-course trajectory correction is an essential technique to ensure the completion of elliptic transfer. Based on explicit guidance theory, this paper use velocity gain guidance, considering perturbation on finite-thrust condition in the midcourse orbit correction method, and a fixed optimization orbit with small terminal deviation and less-used fuel is designed through optimizing fixed time. In the end, mathematical modeling is carried out for different working conditions of the initial deviation to validate the result of this method feasible.

Velocity gain guidance; Explicit guidance; Midcourse trajectory correction; Orbit optimization

V448.11

A

1004-7182(2016)05-0096-05

10.7654/j.issn.1004-7182.20160520

2015-10-13；修改日期：2016-02-23

任金磊（1989-），男，碩士，主要研究方向為飛行器設計-彈道制導導航設計