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

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

（中國運載火箭技術(shù)研究院研究發(fā)展中心，北京，100076）

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

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

（中國運載火箭技術(shù)研究院研究發(fā)展中心，北京，100076）

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

速度增益；顯式制導；中段修正；軌道優(yōu)化

0 引 言

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

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

1 軌道動力學模型

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

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

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

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

2.1 軌道修正策略基本思想

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

圖1 終段控制法

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

2.2 速度增益制導方法

速度增益制導是一種閉路制導。這種閉路制導根據(jù)需要速度對飛行器進行導引和控制，以使其在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 中段修正優(yōu)化軌道方法

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

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

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

表1 不同性能指標的優(yōu)化結(jié)果

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

4 仿真結(jié)果

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

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

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

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

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

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

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

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

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

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

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

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

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

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

表4 不同工況的仿真結(jié)果對比

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

5 結(jié)束語

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

[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]. 中國空間科學技術(shù). 2010, 5: 7-13.

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

[8] 盧山, 段佳佳, 徐世杰. 基于蟻群算法的轉(zhuǎn)移軌道中途修正問題研究[J].系統(tǒng)仿真學報, 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-），男，碩士，主要研究方向為飛行器設計-彈道制導導航設計