阮吉利，李洪軒，邱科悅，郝立麗

（東北林業(yè)大學(xué)，黑龍江 哈爾濱 150040）

“天問一號(hào)”火星探測(cè)器于2021年5月中旬在火星烏托邦平原上順利著陸，我國首次火星探測(cè)任務(wù)取得了巨大的成功。在此之后，“祝融號(hào)”火星車將對(duì)著陸點(diǎn)依次執(zhí)行全球成像、自我檢查、脫離著陸平臺(tái)和巡邏檢測(cè)等工作[1]?！疤靻栆惶?hào)”的成功著陸，是我國航天事業(yè)特別是行星探測(cè)領(lǐng)域發(fā)展的又一個(gè)里程碑，同時(shí)也表明我國成為世界上第二個(gè)著陸火星的國家。本文收集關(guān)于火星探測(cè)器的信息資料，就“天問一號(hào)”的著陸過程建立模型，基于時(shí)間最短、能耗最少的原則，對(duì)操控方案進(jìn)行設(shè)計(jì)，并結(jié)合實(shí)際情況對(duì)操控方案進(jìn)行檢驗(yàn)。

1 “天問一號(hào)”著陸過程耗時(shí)最短控制策略分析

1.1 “天問一號(hào)”落火過程分析

搜集相關(guān)信息后，本文將“天問一號(hào)”著陸過程簡(jiǎn)化為以下5個(gè)階段。

第一階段：兩器分離段。探測(cè)器于火星停泊軌道上運(yùn)動(dòng)，尋覓到合適的著陸時(shí)機(jī)后通過發(fā)動(dòng)機(jī)點(diǎn)火減速降軌。降低到一定高度后，探測(cè)器分離為環(huán)繞器與著陸器，著陸器向火星墜落；環(huán)繞器則在分離半小時(shí)后點(diǎn)火加速升軌，返回停泊軌道。

第二階段：氣動(dòng)減速段。著陸器進(jìn)入火星大氣后，受到空氣阻力的影響，速度大幅減小。

第三階段：傘系減速段。著陸器打開降落傘，由于降落傘與空氣的接觸面積較大，空氣阻力大，能有效減小著陸器的速度。

第四階段：動(dòng)力減速段。探測(cè)器通過發(fā)動(dòng)機(jī)點(diǎn)火的方式進(jìn)一步減小了速度。

第五階段：懸停避障與緩速下降段。在高度降為100 m時(shí)，探測(cè)器的速度基本降到0。探測(cè)器懸停避障，尋找最安全的著陸地點(diǎn)，最后通過4個(gè)著陸腿的緩沖作用在火星表面軟著陸。

1.2 模型建立

為設(shè)計(jì)控制方案使得著陸的時(shí)間最短，主要控制3個(gè)環(huán)節(jié)：環(huán)繞器與著陸器分離；打開阻尼傘的時(shí)機(jī)；發(fā)動(dòng)機(jī)點(diǎn)火的時(shí)機(jī)。3個(gè)環(huán)節(jié)息息相關(guān)，必須綜合決策。查閱資料可得：近火點(diǎn)高度Hc=280 km；遠(yuǎn)火點(diǎn)高度Hf=59 000 km；火星平均半徑R=3 397 km；軌道為橢圓。決定著陸前，探測(cè)器將環(huán)繞火星運(yùn)動(dòng)。

1.2.1 環(huán)繞器與著陸器的兩器分離段

設(shè)探測(cè)器在距離火星h1時(shí)進(jìn)行減速降軌，根據(jù)近火點(diǎn)高度Hc和遠(yuǎn)火點(diǎn)高度Hf可得280 km＜h1＜59 000 km。探測(cè)器在發(fā)現(xiàn)合適的著陸時(shí)間后，通過發(fā)動(dòng)機(jī)噴氣減速。為了簡(jiǎn)化問題的分析，假設(shè)這一減速過程所需時(shí)間可忽略不計(jì)。

根據(jù)萬有引力定律F萬=mg，火星重力加速度為其中m為“天問一號(hào)”在各個(gè)階段的質(zhì)量；M為火星質(zhì)量；G為引力常量。但由于探測(cè)器距離火星的高度相對(duì)于火星半徑不可忽略，則探測(cè)器在距離火星h1處的火星重力加速度為

為方便計(jì)算，本文取g″為降軌時(shí)探測(cè)器所處高度重力加速度與到達(dá)火星大氣時(shí)重力加速度的平均值，即

式中：h2為火星大氣層距離火星地面的高度。

探測(cè)器沿著切線方向具有初速度，且該速度小于火星的第一宇宙速度。而在靠近的過程中，探測(cè)器將分離為環(huán)繞器與著陸器，環(huán)繞器通過加速升軌，不影響著陸器的運(yùn)動(dòng)。綜上，可將探測(cè)器進(jìn)入火星大氣前的運(yùn)動(dòng)視為平拋運(yùn)動(dòng)。

根據(jù)動(dòng)能與勢(shì)能的轉(zhuǎn)換公式，動(dòng)能為

式中：v0為垂直火星方向上探測(cè)器的速度；h=h1-h2。從而推出著陸器進(jìn)入火星大氣時(shí)，探測(cè)器在豎直方向的速度為

探測(cè)器降軌到進(jìn)入火星大氣耗時(shí)t0為

1.2.2 氣動(dòng)減速段和傘系減速段

進(jìn)入火星大氣后，“天問一號(hào)”主要受重力和空氣阻力。著陸器進(jìn)入火星大氣后，高度對(duì)重力加速度的影響可以忽略不計(jì)，因此，取重力加速度g=3.71 m/s2；著陸器所受重力即為G1=mg。

設(shè)著陸器進(jìn)入火星大氣后，經(jīng)過時(shí)間t1打開阻尼傘，傘降時(shí)間t2后進(jìn)入動(dòng)力減速段。為得到著陸器在氣動(dòng)減速段、傘系減速段中，其位移與時(shí)間的關(guān)系，本文做了如下嘗試。下面為便于公式的書寫，令則空氣阻力F=kv2。

在氣動(dòng)減速段，著陸器在下落過程中只受到重力和空氣阻力的作用，且重力的方向與空氣阻力的方向相反。由牛頓第二定律mg-kv2=ma，可得出

通過一定的轉(zhuǎn)化，可將該方程式（6）轉(zhuǎn)化為黎卡提方程[2]，從而解得

由于著陸器長(zhǎng)、寬各約4 m，底盤是圓形[3]，因此迎風(fēng)面積S1約為12.56 m2；火星大氣的空氣密度ρ為0.015 8 kg/m3[4]；空氣阻力系數(shù)C為1.5。因此，上文設(shè)定的

在氣動(dòng)減速段，著陸器進(jìn)入火星時(shí)的初速度即為兩器分離段的末速度。

經(jīng)過時(shí)間t1后，著陸器的速度為

著陸器的運(yùn)動(dòng)路程為

與氣動(dòng)減速段類似，在傘系減速段，著陸器同樣做變加速運(yùn)動(dòng)，主要區(qū)別在于傘系減速段中著陸器展開的降落傘增加了著陸器與空氣的接觸面積，故k值發(fā)生了改變。此時(shí)，降落傘面積S2為300 m2，ρ與C保持不變，計(jì)算得

同理，這一階段（傘系減速段）的初速度為上一階段（氣動(dòng)減速段）的末速度，經(jīng)過時(shí)間t2后，著陸器的速度為

著陸器的運(yùn)動(dòng)路程為

1.2.3 動(dòng)力減速段

考慮到“天問一號(hào)”在動(dòng)力減速段已經(jīng)接近火星表面，且速度比前幾個(gè)階段較低，因此該階段忽略火星大氣阻力對(duì)著陸器的影響[5]。

在動(dòng)力減速段，著陸器扔掉大底和背罩，發(fā)動(dòng)機(jī)點(diǎn)火向地面噴氣，通過反推力減速，著陸器共有27臺(tái)發(fā)動(dòng)機(jī)，包括20臺(tái)250 N、6臺(tái)25 N和1臺(tái)7 500 N的發(fā)動(dòng)機(jī)。故著陸器的最大推力為12 650 N。

達(dá)到動(dòng)力減速段時(shí)，著陸器的速度為v2，其在最大推力的情況下加速度為

由于動(dòng)力減速段著陸器為了懸停，末速度必須為0，因此動(dòng)力減速段所用時(shí)間著陸器在動(dòng)力減速段經(jīng)過的位移x3為

1.3 模型求解

該模型得到的方程組約束條件是決策變量的非線性表達(dá)式，故利用MATLAB軟件編程，先使用蒙特卡羅模型獲得決策變量的初始值，再使用非線性規(guī)劃模型進(jìn)行最優(yōu)化求解[6]。目標(biāo)函數(shù)為

其約束條件為

式中：h=h1-h2，火星半徑R=3 397 km；火星上的重力加速度g=3.71 m/s2；k1=0.149；k2=2.830；探測(cè)器與火星的距離h1為自變量；火星大氣層高度h2=125km。不考慮燃料質(zhì)量，著陸器質(zhì)量m=1 285 kg，模型求解步驟如下。

1）通過蒙特卡羅模型模擬，在給定決策變量范圍內(nèi)生成1億組隨機(jī)數(shù)，從實(shí)驗(yàn)點(diǎn)中找到可行點(diǎn)。

2）通過“打擂臺(tái)”思想，從可行點(diǎn)中找出最小可行點(diǎn)x0，得h1=5 967.365 km；t1=285.77 s；t2=74.68 s；t3=67.43 s。

3）再將得到的x0代入非線性規(guī)劃模型，分別通過內(nèi)點(diǎn)法、序列二次規(guī)劃法、有效集法和信賴域反射算法4種算法，搜索尋求最優(yōu)的決策變量。

4）通過比較4種算法獲得的最優(yōu)決策變量，選取更優(yōu)解。

5）最終選取內(nèi)點(diǎn)法獲得的決策變量。得到h1=6 135.637 km；t0=2 505.3 s；t1=263.67 s；t2=73.23 s；t3=64.65 s；z=2 906.85 s。即“天問一號(hào)”著陸過程的最短總時(shí)長(zhǎng)z為2 906.85 s。

2 “天問一號(hào)”著陸過程耗能最少控制策略分析

2.1 模型建立

該處主要探討在給定的著陸時(shí)間下，著陸器耗能最少的方案，故假定著陸器從降軌到進(jìn)入大氣的持續(xù)時(shí)間t0、氣動(dòng)減速段持續(xù)時(shí)間t1、傘系減速段持續(xù)時(shí)間t2均為已知量。

通過對(duì)探測(cè)器著陸各個(gè)階段的研究，本文認(rèn)為能量消耗的主要階段為降軌減速段以及進(jìn)入火星大氣后的動(dòng)力減速段。

2.1.1 探測(cè)器降軌減速段

當(dāng)探測(cè)器在火星停泊尋找到合適的著陸時(shí)機(jī)后，發(fā)動(dòng)機(jī)將點(diǎn)火制動(dòng)，減小速度以實(shí)現(xiàn)降軌。在該階段中，發(fā)動(dòng)機(jī)點(diǎn)火消耗的能量等價(jià)于探測(cè)器動(dòng)能的變化量。

根據(jù)動(dòng)能定理，該階段著陸器消耗的能量為

為使能量消耗最少，vt應(yīng)該盡量大；但同時(shí)vt又應(yīng)該小于火星的第一宇宙速度，否則著陸器將無法墜向火星，故取

依據(jù)著陸器從降軌到進(jìn)入大氣的持續(xù)時(shí)間t0，可知探測(cè)器降軌時(shí)的停泊軌道距離火星大氣的高度由此可得探測(cè)器距離火星的距離h1=h0+h2=h0+125 m。

由開普勒定律，任何橢圓天體軌道的中心天體一定在橢圓的一個(gè)焦點(diǎn)上。同時(shí)根據(jù)橢圓的定義b2=a2-c2，得出方程組為

聯(lián)系萬有引力定律與牛頓第二定律，可以列出“天問一號(hào)”在近火點(diǎn)和遠(yuǎn)火點(diǎn)的運(yùn)動(dòng)學(xué)方程組為

式中：vc和vf分別為著陸器在近火點(diǎn)的速度和在遠(yuǎn)火點(diǎn)的速度；ρc與ρf分別為近火點(diǎn)和遠(yuǎn)火點(diǎn)的曲率半徑，聯(lián)立以上各式，代入?yún)?shù)，得近火點(diǎn)速度vc=4.691 5×103m/s，遠(yuǎn)火點(diǎn)速度vf=890.538 2 m/s。

根據(jù)開普勒第二定律可得

根據(jù)能量守恒定律可得

由此得著陸器在火星海拔h1上的速度v0為

2.1.2 動(dòng)力減速段

著陸器在進(jìn)入大氣后將先后經(jīng)過氣動(dòng)減速段、傘系減速段以及動(dòng)力減速段。

根據(jù)傘系減速段的運(yùn)動(dòng)規(guī)律，著陸器受到的空氣阻力與其運(yùn)動(dòng)方向相反，且該阻力的值與速度的平方成正比。因此，起初著陸器受到的空氣阻力大于其重力，著陸器做減速運(yùn)動(dòng)。而空氣阻力的值隨著速度的減小而減小，直到與重力達(dá)成平衡，著陸器將近似做勻速運(yùn)動(dòng)[7]。由于t0、t1為已知量，因此根據(jù)式（8）可以算出著陸器做勻速運(yùn)動(dòng)時(shí)的速度，不妨令其為v。

在傘系減速到某一高度后，著陸器將通過發(fā)動(dòng)機(jī)點(diǎn)火來實(shí)現(xiàn)減速。假設(shè)著陸器在距離火星地面高度為h3時(shí)開始進(jìn)行動(dòng)力減速，此時(shí)速度為v，為了進(jìn)行后續(xù)的懸停避障，在距離地面高度為100 m時(shí)著陸器速度應(yīng)減小為0。根據(jù)能量守恒定律，在動(dòng)力減速段中，發(fā)動(dòng)機(jī)所做的功為

設(shè)發(fā)動(dòng)機(jī)推力為F推（0 N＜F推＜12 650 N）。動(dòng)力減速段，豎直方向上著陸器受重力以及發(fā)動(dòng)機(jī)推力。因此，著陸器在動(dòng)力減速段做勻減速運(yùn)動(dòng)，加速度為

為滿足著陸器距離地面100 m時(shí)懸停的條件，可得

2.2 模型求解

2.2.1 降軌減速段求解

在探測(cè)器的降軌減速段，本文通過開普勒定律、牛頓第二定律等對(duì)探測(cè)器在火星停泊軌道上的運(yùn)動(dòng)做了分析，在著陸器從降軌到進(jìn)入大氣的時(shí)間已知的前提下，推導(dǎo)出了探測(cè)器降軌時(shí)的速度。并在基于能量消耗最小化的考量下，得到了探測(cè)器降軌前后的速度，從而利用動(dòng)能定理得到了降軌減速段的能量消耗。

整理相關(guān)公式可得方程組為

此處取t0為1 000 s，解得v0=4.396 3×103m/s；E1=1.678 2×107kJ。

2.2.2 動(dòng)力減速段求解

有關(guān)動(dòng)力減速段所描述的問題，可由一個(gè)有約束條件的線性規(guī)劃問題解釋。因此，本文用遺傳算法對(duì)本問題的模型進(jìn)行求解[8-9]。

目標(biāo)函數(shù)為

其約束條件為

模型求解步驟如下。

1）將決策變量排列為1個(gè)個(gè)體，初始種群由2 000個(gè)隨機(jī)產(chǎn)生的個(gè)體組成。

2）編寫適應(yīng)函數(shù)，來計(jì)算個(gè)體的性能指標(biāo)。

3）選擇出指標(biāo)較好的個(gè)體留存，并遺傳給下一代。

4）隨機(jī)將個(gè)體兩兩配對(duì)，并按照規(guī)定的交叉概率進(jìn)行交叉。

5）若最大進(jìn)化代數(shù)達(dá)到了20 000，則進(jìn)行最優(yōu)解的輸出；否則繼續(xù)進(jìn)行編碼、評(píng)價(jià)、選擇、交叉等操作。

通過約束條件的分析，本問題的模型確定了交叉概率為0.6。對(duì)模型進(jìn)行求解，得到v取90 m/s時(shí)，h的最優(yōu)解為660.21 m。此時(shí)，由于E2=8 351.702 kJ，因此E=E1+E2=1.679×107kJ。

3 “天問一號(hào)”著陸模型與實(shí)際情況的匹配度檢驗(yàn)

探測(cè)器的著陸數(shù)據(jù)如下。

1）兩器分離段結(jié)束，著陸器飛行3 h后進(jìn)入火星大氣（火星大氣層高度為125 km）。

2）氣動(dòng)減速段：著陸器在約290 s內(nèi)，將速度從4.8 km/s減速到460 m/s。

3）傘系減速段：在距離火星地面11 km左右處，探測(cè)器打開降落傘，在約90 s內(nèi)將速度由460 m/s降至95 m/s。

4）動(dòng)力減速段：距離火星地面1.5 km左右處，探測(cè)器通過發(fā)動(dòng)機(jī)進(jìn)一步減速，用80 s左右的時(shí)間將速度減小到3.6 m/s。

5）著陸緩沖段：在距離火星表面約100 m時(shí)，著陸器的速度基本降為0。

前文對(duì)著陸器的著陸過程建立了運(yùn)動(dòng)模型。為驗(yàn)證與評(píng)估該模型的科學(xué)性，本文代入3個(gè)參數(shù)值：進(jìn)入火星大氣時(shí)的速度（4 800 m/s）、氣動(dòng)減速段持續(xù)的時(shí)間（290 s）、傘系減速段持續(xù)的時(shí)間（90 s），則用于求解的方程組為

解得氣動(dòng)減速段的位移為103.120 1 km，傘系減速段的位移為11.642 2 km，動(dòng)力減速段的位移為1.337 7 km。對(duì)比“天問一號(hào)”著陸時(shí)在這3個(gè)階段的實(shí)際位移114 km、9.5 km、1.5 km，利用公式求得相對(duì)誤差分別為9.54%、22.6%、10.82%。其中，Er為相對(duì)誤差；T為真實(shí)值。由此可知，本文建立的模型與實(shí)際的著陸過程存在可忽略不計(jì)的誤差。

隨后，本文取10組著陸器進(jìn)入火星大氣時(shí)速度與氣動(dòng)減速段持續(xù)時(shí)間的數(shù)據(jù)，代入模型運(yùn)算，得到氣動(dòng)減速段的位移，對(duì)取得的數(shù)據(jù)擬合后，得到計(jì)算值與真實(shí)值的擬合圖（見圖1）。

圖1 計(jì)算值與真實(shí)值的擬合圖

由圖1可知，真實(shí)值與擬合所得曲面的計(jì)算值較為接近。綜上，可認(rèn)為本文對(duì)于著陸過程所建立的運(yùn)動(dòng)模型與實(shí)際情況較相符，準(zhǔn)確度較高。

4 結(jié)論

本文通過數(shù)學(xué)和物理方法，對(duì)“天問一號(hào)”探測(cè)器著陸的各個(gè)過程進(jìn)行分析，建立運(yùn)動(dòng)模型方程。首先，以著陸時(shí)間最短為優(yōu)化目標(biāo)，明確自變量、約束條件后，建立了單目標(biāo)優(yōu)化模型，用非線性規(guī)劃模型進(jìn)行最優(yōu)化求解，解得著陸的最短總時(shí)長(zhǎng)為2 906.85 s。其次，在給定著陸時(shí)間的前提下，建立以能量消耗最小為目標(biāo)，以發(fā)動(dòng)機(jī)推力、點(diǎn)火時(shí)間為自變量，以著陸器在火星海拔100 m時(shí)速度為0為約束條件的優(yōu)化模型，并利用遺傳算法進(jìn)行了規(guī)劃求解，解得在確立各個(gè)階段用時(shí)的前提下，著陸最少消耗能量為1.678 2×107kJ。最后，為檢驗(yàn)運(yùn)動(dòng)模型的準(zhǔn)確性，本文搜集了“天問一號(hào)”著陸的真實(shí)值數(shù)據(jù)，代入模型計(jì)算后發(fā)現(xiàn)誤差較小，具有一定可信性。