楊永輝，劉 睿，李 華

(1.西北核技術(shù)研究所，陜西 西安 710024；2.西安航空學院 理學院，陜西 西安 710077)

嫦娥三號軟著陸軌道設計與控制策略

楊永輝1，劉 睿2，李 華2

(1.西北核技術(shù)研究所，陜西 西安 710024；2.西安航空學院 理學院，陜西 西安 710077)

針對嫦娥三號軟著陸的軌道設計和控制策略進行研究。首先，采用開普勒第一、二定律求出嫦娥三號在近月點與遠月點的相對速度與著陸點的相對位置。其次,將嫦娥三號軟著陸過程分為六個階段,利用動力學微分方程，針對每個階段，建立消耗燃料最少為目標的優(yōu)化模型。由開關(guān)控制理論與極小值原理,得到六個階段的點火時間和推力。最后，通過穩(wěn)定性分析，可知本文所設計的控制策略軌道具備誤差小、敏感性低的優(yōu)點。

開關(guān)控制；平滑度矩陣；搜索算法

0 引言

開發(fā)月球資源，建立月球基地已成為世界航天活動的必然趨勢與競爭熱點。我國已邁出航天深空探測第一步，開展月球探測工作。而在月球表面實現(xiàn)軟著陸則是月球勘探的重要前提。嫦娥三號于2013年12月2日1時30分成功發(fā)射，12月6日抵達月球軌道。嫦娥三號在著陸準備軌道上的運行質(zhì)量為2.4t，其安裝在下部的主減速發(fā)動機能夠產(chǎn)生1500N到7500N的可調(diào)節(jié)推力，其比沖(即單位質(zhì)量的推進劑產(chǎn)生的推力)為2940m/s，可以滿足調(diào)整速度的控制要求。在四周安裝有姿態(tài)調(diào)整發(fā)動機，在給定主減速發(fā)動機的推力方向后，能夠自動通過多個發(fā)動機的脈沖組合實現(xiàn)各種姿態(tài)的調(diào)整控制。嫦娥三號的預定著陸點為19.51W，44.12N，海拔為-2641m。嫦娥三號已于北京時間12月14號在月球表面實施軟著陸，其如何實現(xiàn)軟著陸已成為全世界關(guān)注的焦點。

Apollo系列則是通過宇航員觀測著陸區(qū)并操縱人控系統(tǒng)實現(xiàn)了避障和安全著陸。鳳凰號探測器成功實現(xiàn)了火星軟著陸,其采用了事先篩選高概率安全著陸區(qū)來避免大障礙的危害[1]。NEAR探測器首次實現(xiàn)了小行星著陸任務,但由于姿控誤差較大等原因,探測器未能實現(xiàn)避障就繼續(xù)向小行星表面下降了[2-8]。

嫦娥三號的核心任務是實施高可靠高安全的月面軟著陸,要求著陸器必需具備自主障礙識別與規(guī)避能力。本文針對嫦娥三號軟著陸的軌道設計和控制策略進行研究。

首先，采用開普勒第一、二定律求出嫦娥三號在近月點與遠月點的相對速度與著陸點的相對位置。

其次,將嫦娥三號軟著陸過程分為六個階段,利用動力學微分方程，分別針對每個階段，建立消耗燃料最少為目標的優(yōu)化模型。由開關(guān)控制的理論與極小值原理,得到六個階段的點火時間和推力。

最后，通過穩(wěn)定性分析，可知本文所設計的控制策略軌道具備誤差小、敏感性低的優(yōu)點。

1 嫦娥三號軟著陸過程的六個階段

嫦娥三號軟著陸過程分為六個階段，如圖1所示。

圖1 嫦娥三號軟著陸過程示意圖

2 嫦娥三號在近月點與遠月點的相對速度與著陸點的相對位置

為了確定著陸準備軌道中近月點和遠月點的相對速度與相對位置，本文將其轉(zhuǎn)化為求解主減速階段與快速調(diào)整階段中探月器產(chǎn)生水平位移之和的問題。

圖2 近月點與遠月點示意圖

2.1 嫦娥三號在近月點和遠月點的速度

由開普勒第一定律知，嫦娥三號圍繞月球的運動軌道為橢圓軌道,在由近月點A到遠月點B的過程中，根據(jù)機械能守恒定律得：

(1)

由開普勒第二定律得到：

(2)

其中，m為嫦娥三號的質(zhì)量，M為月球質(zhì)量，RA為近月點到月心的距離，RB為遠月點到月心的距離，vA為嫦娥三號在近月點的速度，vB為嫦娥三號在遠月點的速度。由式(2)可得：

vA=1.69km/s，方向沿橢圓軌道的切線方向指向著陸點；vB=1.14km/s，方向與近月點速度方向相反。

2.2 近月點及遠月點的位置

嫦娥三號在近月點、遠月點的位置取決于主減速階段、調(diào)整階段的嫦娥三號水平位移。

2.2.1 嫦娥三號在主減速階段水平位移

對嫦娥三號在主減速階段降落的分析得圖3：

圖3 主減速階段示意圖

此階段嫦娥三號距月球表面的海拔變化為15000m～3000m。t為嫦娥三號從近月點到距月球表面某一高度的時間間隔，u為嫦娥三號發(fā)動機的推力，v為嫦娥三號的瞬時飛行速度，g為月球表面的重力加速度，v0為嫦娥三號在近月點的速度，即：v0=vA，v1為嫦娥三號到達距離月球表面3000m時的水平速度，v2為嫦娥三號到達距離月球表面3000m時豎直速度，v12為嫦娥三號到達距離月球表面3000m的速度，ve是以米/秒為單位的比沖，x1為此階段嫦娥三號的水平位移。則：

(3)

利用Matlab可求得非線性方程的解為：

v1=56.0425m/s，v2=11.0699m/s，t1= 137.2116s，x1=214667.5m

2.2.2 嫦娥三號在快速調(diào)整階段降落過程中的水平位移

此階段嫦娥三號距月球表面的海拔變化為3000m～2400m。假設此階段的推力不變，且嫦娥三號到距月球表面2400m時，水平速度為零。為保證降落過程中燃料消耗最少，此階段的推力方向平行于月球表面。

設u1為此階段發(fā)動機的推力，x2為此階段嫦娥三號的水平位移。則

(4)

利用Matlab中的牛頓迭代法求非線性方程的解：

t2=30.6822m/s，u1=4055.8N，x2=859.8m

近月點到著陸點之間的水平距離為：x=x1+x2=215530m

由3.2.1與3.2.2可知，近月點距著陸點的水平距離為215.5km，豎直方向距離15km。遠月點距著陸點的水平距離為215.5km，豎直方向距離100km。

3 嫦娥三號軟著陸軌道設計與控制策略

根據(jù)圖1所示的六個階段，逐段進行軟著陸軌道與控制策略的設計。

3.1 主減速階段

主減速階段的狀態(tài)方程如下：

(5)

(6)

終端的橫截條件為：

(λ1(tf),λ2(tf),λ3(tf))=(γ1,γ2,-1)

(7)

得最優(yōu)控制u*(t)為：

圖4 水平速度分量的v-t圖像

圖5 豎直速度分量的v-t圖像

(8)

其中M′為嫦娥三號無燃料時的質(zhì)量。

m(t)=-ka(t-τ)+m(τ)

(9)

為了直觀形象地研究探月器的運動，在相平面(v(t),h(t))上繼續(xù)討論探月器的運動軌跡。

在主減速階段的后半段，u*(t)=a，且當t=tf時h(tf)=3000，v(tf)=57，tf-τ=td，代入式(9)中得到

(10)

(11)

在主減速階段前半段u=0，由(9)得

兩條運動軌跡作圖如下：

圖6 運動軌跡

兩條運動軌跡的交點為最優(yōu)控制的切換點，由圖6可以看出交點的橫坐標無限接近于坐標原點，可認為從探月器進入主減速階段的時刻開始，開關(guān)函數(shù)控制主發(fā)動機提供最大推力。

主減速階段嫦娥三號從近月點開始點火，主發(fā)動機以最大推力來提供減速動力，主減速階段結(jié)束時，水平速度v1=56.0425m/s，豎直速度v2=11.0699m/s，歷時t1= 137.2116s，嫦娥三號的質(zhì)量：m1=2262.8kg。

3.2 快速調(diào)整階段

此階段嫦娥三號距月球表面的海拔變化為3000m～2400m。根據(jù)主減速階段的最優(yōu)控制策略，我們將探月器在快速調(diào)整階段主發(fā)動機提供的減速動力也進行開關(guān)控制，使發(fā)動機在切換時刻τ之前提供恒定的水平最大推力，使之水平速度減小到零，而在切換時刻τ之后不提供推力，以此來減小燃料消耗。

以月球表面3000米的時刻為零時刻，建立動力學方程：

求解得到：x3=464.93m，t3=16.6s，h1=358.23m

3.3 粗避障階段

此階段嫦娥三號距月球表面的海拔變化為2400m～100m。分為兩個部分，一個部分是控制啟動主發(fā)動機控制豎直推力的時刻，使燃料消耗最少；另一部分是啟動姿態(tài)調(diào)整發(fā)動機控制水平推力，使探月器粗步避開大隕石坑。

3.3.1 控制豎直推力

根據(jù)最優(yōu)控制策略，減小燃料消耗，將探月器在粗避障階段主發(fā)動機提供的減速動力也進行開關(guān)控制，使發(fā)動機在切換時刻τ之前不提供推力，而在切換時刻τ之后，提供恒定的豎直向上的最大推力。

根據(jù)上述減速動力控制原則，嫦娥三號在豎直方向上先進行自由落體運動，再進行加速度不斷變大的減速運動。當?shù)竭_距月球表面100m時豎直方向上的速度減小到零，并保持懸停。

設t5表示自由落體的時間，制動時間為t6，h2表示自由落體過程中下降的高度。

由Matlab解得：t5=30.78s，h2=2041m，t6=16.76s。到達距月球表面距離為100米時，嫦娥三號的質(zhì)量m3=2169.1kg。

3.3.2 控制水平推力

將2400m處高程圖(圖7)導入到2300×2300的矩陣A中，

并計算矩陣中各點的平滑度η：

得到平滑度矩陣

以2400米處高程圖的圖像中心點為中心輻射，搜索平滑度矩陣中4×4范圍內(nèi)平滑度之和小于篩選精度μ1的點，令μ1=0.5，求得圖像中心左上方點的位置坐標。

以此類推，找到圖像全部范圍內(nèi)符合篩選精度的點，從而找到距離圖像中心點最近的點作為粗避障的目標點，坐標為(1150,1145)，距預定著陸點的水平距離為5m。

圖7 2400m處高程圖

3.3 精避障階段

將100m處的高程圖的圖像(圖8)導入到1000×1000的矩陣B中，

計算矩陣中各點的平滑度η(同粗避障階段)，得到平滑度矩陣，以圖像中心點輻射，搜索平滑度矩陣中10×10范圍內(nèi)平滑度之和小于篩選精度μ的點，此時我們設定μ=10，求得圖像中心左上方所有點的位置坐標。

圖8 100m處高程圖

以此類推，找到圖像全部范圍內(nèi)符合篩選精度的點，從而找到距離圖像中心點最近的點作為精避障的目標點，即最終的實際著陸點，坐標為(445,500)，距預定著陸點的垂直距離為5.5m。

3.4 緩速下降階段

此階段嫦娥三號距月球表面的海拔變化為100m～4m，根據(jù)燃料消耗最小原則，將探月器在上述兩個階段主發(fā)動機提供的減速動力進行開關(guān)控制，使發(fā)動機在切換時刻τ之前不提供推力，而在切換時刻τ之后，提供恒定的豎直向上的最大推力。

根據(jù)上述減速動力控制原則，嫦娥三號在豎直方向上先進行自由落體運動，再進行加速度不斷變大的減速運動。當?shù)竭_距月球表面100m時豎直方向上的速度減小到零，并保持懸停。

此過程中的自由落體的時間t7=10.4s，制動時間t8=5.2s，此過程中自由落體過程中下降的高度：h3=70.3m。

3.5 嫦娥三號軟著陸軌道設計與控制策略

六個階段控制策略如表1所示：

表1 嫦娥三號軟著陸控制策略表

4 誤差與靈敏度分析

4.1 誤差分析

本文用實際著陸點和預定著陸點的直線距離來刻畫偏離預定著陸點位置的程度。

在快速調(diào)整階段，距月球表面2400米高度時水平方向產(chǎn)生的偏差：

x1=859.8-464.93=394.87m

在避障階段(包括粗避障和精避障)，實際著陸點與預定著陸點之間的水平距離x2=5m，垂直距離x3=5.5m

4.2 敏感度分析

在粗避障與精避障階段，改變μ1,μ2的值，可得不同的篩選數(shù)量，見表2：

表2 粗避障與精避障篩選精度與篩選數(shù)量關(guān)系表

參數(shù)μ1的敏感度為0。

[1]BonfiglioEP,AdamsD,CraigL,etal.LandingsitedispersionanalysisandstatisticalassessmentfortheMarsPhoenixLander[C]//AstrodynamicsSpecialistConferenceandExhibit.Hawaii:AIAA,2008:7348.

[2]KubotaT,HashimotoT,KawaguchiJ,etal.GuidanceandnavigationofHayabusaspacecraftforasteroidexplorationandsamplereturnmission[C]//InternationalJointConference.Busan:SICE-ICASE,2006:2793-2796.

[3] 劉少創(chuàng)，賈陽，馬友青,等.嫦娥三號月面巡視探測器高精度定位[J].科學通報,2015(4):372-378.

[4] 趙育善，師鵬.航天器飛行動力學建模理論與方法[M].北京：北京航空航天大學出版社，2012:85-115.

[5] 吳滄浦，夏元清，楊毅.最優(yōu)控制的理論與方法[M].北京：國防工業(yè)出版社，2013:105-134.

[6] 胡壽松，王執(zhí)銓，胡維禮.最優(yōu)控制的理論與系統(tǒng)[M].北京：科學出版社，2005:46-76.

[7] 馬明亮，王超，施潤和，等.嫦娥一號IIM高光譜數(shù)據(jù)和月球軌道器LOLADEM數(shù)據(jù)的配準與月表地形校正及評價[J].地球信息科學學報,2015(1):118-125.

[8] 李培佳，黃勇，昌勝騏,等.基于地基觀測的嫦娥三號著陸器與巡視器高精度定位[J].科學通報,2014(32):3162-3173.

[責任編輯、校對：周 千]

Trajectory-design and Control-strategy of Chang'E-3 Soft-landing

YANGYong-hui1,LIURui2,LIHua2

(1.Northwest Institute of Nuclear Technology, Xi'an 710024,China;2.Faculty of Science, Xi'an Aeronautical University, Xi'an 710077, China)

In this paper, the trajectory-design and control-strategy of the Chang'E-3 soft-landing is studied.At first, the relative velocity and position at the pericynthion and apocynthion are given by Kepler's first and second laws.Then, the Chang'E-3 soft-landing process is divided into six phases.For every phase, the optimization model whose object is minimum fuel consumption is found by the dynamic differential equation.We get the ignition time and thrust of six phases by switch control theory and minimal value principle.At last, it is concluded through stability analysis that the trajectory-design and control-strategy of the Chang'E-3 soft-landing in this paper has advantages such as smaller error and lower sensitivity.

switch control theory; smoothness matrix; search algorithm

2015-03-26

楊永輝(1980-)，男，山西聞喜人，工程師，從事空間技術(shù)方面的研究。

V

A

1008-9233(2015)03-0014-06