火星EDL過程動力學仿真平臺

唐明章 王奇 黃偉 錢凱

火星EDL過程動力學仿真平臺

唐明章1王奇1黃偉1錢凱2

（1 北京空間機電研究所，北京 100094）（2 北京東方創達軟件有限公司，北京 100089)

“進入、減速、著陸（Entry，Descent，Landing，EDL）”火星的過程是實現整個火星探測任務最為重要的階段之一，著陸成功與否直接決定探測任務的成敗。文章對火星探測器自身氣動減速階段、降落傘彈射拉直階段、降落傘開傘充氣階段、降落傘全張滿減速階段等EDL關鍵過程分別建立探測器進入段六自由度動力學模型、降落傘理想拉直過程動力學模型、降落傘開傘充氣過程動力學模型、物傘及彈性吊掛組合系統著陸過程動力學模型。并在此基礎上研制相應的動力學仿真軟件模塊，將各模塊按照EDL過程進行集成，實現模塊間數據傳遞，形成火星EDL全過程動力學仿真系統。最后對“火星探路者”EDL過程進行了動力學仿真分析，校驗了系統的有效性和可行性。該平臺可用于火星EDL全過程動力學性能的分析與預測。

火星探測 進入、減速、著陸動力學 仿真分析 軟件設計

0 引言

2020年7月23日，我國首顆火星探測器“天問一號”發射升空。經歷10個月太空之旅后，于2021年5月15日成功著陸火星烏托邦平原南部預選著陸區。

探測器從進入火星大氣層至著陸主要經歷三個階段：大氣進入段（Entry）、降落傘減速段（Descent）和著陸段（Landing），這三個階段就是火星的“進入、減速、著陸（EDL）”過程。火星EDL過程是整個火星探測任務最為關鍵的環節[1-2]，著陸成功與否直接決定探測任務的成敗。

整個火星EDL過程持續6～8min。此間探測器將經受嚴酷的峰值過載、峰值動壓以及惡劣的氣動加熱環境的考驗。另外，探測器進入初始狀態參數和進入過程中火星復雜多變的環境參數直接影響探測器的著陸精度。因此對探測器和降落傘組成的物傘系統，在火星EDL過程中的動力學性能進行預測與分析有著非常重要的工程意義。文獻[3-4]以“火星探路者”為研究對象，建立了物傘系統工作過程的動力學模型，系統地研究了艙傘系統減速下降過程的運動學特性和動力學特性。對于EDL全過程動力學的仿真研究，美國宇航局噴氣推進實驗室已開發專業仿真軟件DSENDS[5-6]。國內，文獻[7]對傘降過程中物傘系統建立了六自由度模型，并在此基礎上對火星EDL過程進行了虛擬環境下的展示。文獻[8]在降落傘減速階段采用九自由度動力學模型，建立了火星EDL動力學集成仿真框架，實現了EDL過程的三維顯示。

本文對物傘組合體著陸過程建立了精度更高的十二自由度動力學模型，同時對火星EDL其它關鍵過程分別建立了探測器進入段六自由度動力學模型、降落傘理想拉直過程動力學模型、降落傘開傘充氣過程動力學模型。并開發相應的動力學仿真軟件模塊，將各模塊集成實現EDL全過程動力學仿真集成軟件。通過對“火星探路者”EDL過程的仿真分析驗證了軟件的有效性和可行性。

1 EDL過程動力學建模

1.1 進入段六自由度動力學模型

進入段從探測器進入火星大氣層開始，到降落傘開傘時結束。進入段探測器受氣動力和火星引力作用，探測器質心運動方程

式中為探測器質量；為探測器速度；為時間；為火星地心至探測器質心的矢徑；為火星自轉角速度；為探測器氣動力；為探測器所受火星引力。將式（1）中各矢量投影至航跡坐標系下即可建立進入段航跡坐標系下的質心動力學方程[9]。

在探測器坐標系下建立繞質心的轉動方程[10]：

1.2 降落傘開傘過程動力學模型

降落傘拉直過程是由引導傘將折疊于傘包內的降落傘按序拉出的過程，也是降落傘開傘過程的第一步。本文采用Wolf[11]提出的理想拉直動力學模型：即在拉直過程中，傘衣和傘繩從傘包按序連續拉出，并且認為在拉直過程中回收物和引導傘遵循同一彈道傾角。

降落傘充氣過程是指從傘系統全長拉直到傘衣第一次充滿為止的整個工作過程。充氣過程異常復雜，在很短時間內降落傘外形、質量分布等急劇變化。同時降落傘各部分力學特性也發生劇烈而復雜變化。目前尚無完整的理論模型對充氣過程進行精確的計算和分析，通常采用經驗模型來計算[12]。

1.3 物傘組合體著陸過程動力學模型

物傘系統是一個多體系統，動力學建模時可對物、傘分別建立剛體六自由度模型，形成十二自由度的組合系統。物傘之間由彈性吊帶連接，通過平衡法可求得吊帶系統在每個瞬時的約束力。

1.3.1 降落傘及探測器動力學模型

物傘組合體著陸過程中對降落傘模型做如下假設：

1）降落傘處于充滿狀態，且傘衣保持軸對稱形狀不變；

2）傘衣的壓心與幾何中心重合；

3）用附加質量表示降落傘非定常運動產生的流體附加力和附加力矩；

4）忽略尾流對降落傘的影響。

在降落傘坐標系建立六自由度運動方程如下所示

1.3.2 吊掛系統模型

降落傘與探測器之間通過彈性吊掛系統相連。圖1為常用的兩種吊掛連接方式，為傘繩交匯點，1、2、…、為吊掛系統與探測器連接點，下標為自然數。圖1（b）吊掛中間有一個中間自由點，與傘繩交匯點由吊帶連接。

圖1 吊掛結構示意

吊掛系統約束力的求解通常有小質量點法和“平衡法”兩種方法。小質量點法將吊掛系統中間節點視為一個節點，由于中間節點質量很小，容易使運動方程出現剛性，嚴重限制了積分時間步長。故本文采用“平衡法”求解吊掛系統中傘繩、吊帶、連接帶的約束力。

1）對于無中間自由點的吊掛連接方式，在傘繩交匯點建立力的平衡方程

2）對于有中間自由點的吊掛結構，分別在和點建立平衡方程[15-16]

2 火星EDL動力學仿真系統設計

2.1 系統架構

系統總體架構采用三層架構設計，即界面層、業務邏輯層、數據訪問層。

界面層指與用戶交互的界面，用于接收用戶輸入的數據和顯示處理后用戶需要的數據。它包括進入段彈道計算參數輸入界面、降落傘拉直及開傘過程計算界面、物傘著陸過程計算界面、動力學及運動學數值結果和時程曲線顯示界面等等。

業務邏輯層是界面層和數據訪問層之間的橋梁，實現業務邏輯。它包括進入段六自由度彈道分析、降落傘拉直及開傘充氣過程動力學分析、物傘組合體著陸過程動力學分析及后處理模塊。

數據訪問層是實現對數據的增、刪、改、查功能。將數據操作結果提交給業務邏輯層，同時保存業務邏輯層的數據。系統涉及數據包括：地球大氣數據、火星大氣數據、探測器模型數據、減速傘模型數據、探測器及降落傘氣動數據、風場數據等等。

孩子的感冒與成人有所不同，這是他們的身體特點所決定的。小兒肺脾常不足，而心肝相對有余，所以感受外邪之后，也容易影響心肝脾的正常功能，造成感冒夾滯、夾驚、夾痰3種表現：

2.2 系統模塊組成

依據EDL過程不同階段將系統分為4個動力學仿真子模塊及1個后處理模塊，共5個子模塊，如圖 2所示。其中，4個動力學仿真模塊為：再入六自由度動力學仿真模塊、降落傘彈射拉直過程動力學仿真模塊、降落傘充氣過程動力學仿真模塊、物傘組合體著陸過程動力學仿真模塊。

圖2 系統模塊組成

2.3 系統流程

系統流程如圖3所示。首先給定探測器進入段初始狀態參數：經度、緯度、高度、速度、飛行路徑角、速度方位角、角速度、姿態角等。啟動進入段六自由度動力學模塊進行彈道積分計算，得到每一瞬時探測器的狀態參數。當達到開傘條件時，停止當前模塊的仿真任務。并啟動下一個仿真模塊，即調用降落傘彈射拉直過程動力學仿真模塊進行仿真計算，仿真的初始狀態參數從上一個模塊的末端狀態參數獲取。其余模塊的執行與此類似。

圖3 系統流程及模塊調用示意

動力學仿真模塊執行中均需要調用星球模型：火星重力模型、大氣模型、風場數據模型，以及飛行器模型。飛行器模型包括：氣動模型、質量及附加質量特性模型。飛行器類型包括探測器、降落傘和物傘組合體三種，其中物傘組合體著陸動力學仿真模塊還需調用吊掛模型。

3 “火星探路者”EDL過程仿真計算

3.1 進入段六自由度模塊驗證

以“火星探路者”探測器為研究對象，采用文獻[17-19]數據對進入段六自由度模型進行驗證。“火星探路者”氣動性能數據來自文獻[20]，仿真初始狀態參數如表1所示。

表1 彈道初始條件

Tab.1 Initial conditions of trajectory

圖4給出了探測器高度—速度變化關系與文獻[18]仿真結果的對比，可以看到本文仿真結果與文獻仿真結果一致。圖4中，在60km高度以上，由于大氣密度較低，氣動減速不明顯。主要減速段在60km高度以下。

3.2 物傘組合體動力學模塊驗證

物傘組合體模型由“火星探路者”號探測器、盤縫帶傘和彈性吊掛組成，如圖5所示，其中盤縫帶傘參數見文獻[21]，吊掛系統考慮有中間自由點和無中間自由點兩種模型，連接帶參數見文獻[14]。

圖4 進入段高度速度曲線

圖5 物傘組合體及吊掛模型示意

本文以表2作為初始狀態，進行物傘組合體動力學仿真計算，并與文獻[21]結果進行對比驗證。

表2 仿真初始條件

Tab.2 Initial conditions of simulation

文獻[21]同時給出了計算結果和加速度計測量數據的處理結果。圖 6為本文探測器高度隨時間的變化計算結果與文獻的對比，可以看到，本文仿真結果與文獻仿真結果一致，但與文獻中的測量數據存在一定的偏差。文獻指出誤差可能是由于仿真用的大氣模型和火星實際大氣參數存在一定的差異。另外，吊掛系統為物、傘之間的內部約束，不影響物傘系統整體運動。圖 6中可以看到，兩種吊掛模型下，探測器高度隨時間變化曲線幾乎是重合的。

圖7為探測器下降速率隨時間的變化。同樣，由于吊掛系統為內部約束，圖中兩種吊掛模型對應的探測器垂直下降速率曲線幾乎重合。在前10s，本文結果、文獻結果、測量數據之間相差不大。10s之后出現較大的偏差，總體上垂直速度的計算結果要大于測量結果，這也與圖6中高度仿真結果比測量結果下降更快相一致。與文獻仿真結果的偏差可能是由于本文吊掛系統參數、降落傘結構和性能參數等與文獻存在差別所致。

圖6 探測器高度隨時間的變化

圖7 探測器垂直下降速度隨時間的變化

3.3 火星EDL全過程動力學仿真

進入段初始狀態參數見表1。降落傘充氣過程阻力面積隨時間的變化參見文獻[22]。物傘組合體模型與3.2節相同。

利用本文所建立的動力學模型，通過給定的進入段初始狀態，可獲得EDL全過程的動力學和運動學特性。其中圖8為探測器高度隨時間的變化關系。圖9為探測器速度隨時間的變化關系。兩者變化趨勢均與文獻[18]一致。

圖8 高度隨時間的變化

圖9 速度隨時間的變化

圖10為探測器加速度隨時間的變化關系。圖中，最大過載在進入段69.4s，最大過載值為17.18n，對應的飛行速度為176.36m/s。這與文獻[23]中最大過載17.3n，及對應飛行速度169.5m/s非常接近。另外，本文開傘過程最大過載8n，文獻中開傘過程最大過載不到6n。這一差異可能是由于充氣過程中阻力面積的變化與文獻不同導致。

圖10 加速度隨時間的變化

以上分析表明，本文所開發的動力學集成軟件可以對火星EDL全過程動力學特性進行有效地預測與分析。

4 結束語

本文對火星EDL過程不同階段，分別建立了進入段六自由度動力學模型、降落傘拉直和充氣過程動力學模型、物傘組合體著陸過程動力學模型。并開發了相應的動力學仿真軟件模塊，將各模塊集成形成火星EDL全過程動力學仿真軟件系統。通過對“火星探路者”和盤縫帶傘組成的系統在火星EDL過程中的動力學特性進行仿真分析，驗證了軟件的有效性和可行性。

本文所研制動力學軟件可用于火星EDL過程動力學特性的預測與分析。

[1] 饒偉, 孫澤洲, 孟林智, 等. 火星著陸探測任務關鍵環節技術途徑分析[J]. 深空探測學報, 2016, 3(2): 121-128.

RAO Wei, SUN Zezhou, MENG Linzhi, et al. Analysis and Design for the Mars Entry, Descent and Landing Mission[J]. Journal of Deep Space Exploration, 2016, 3(2): 121-128. (in Chinese)

[2] BRAUN R D, MANNING R M. Mars Exploration Entry, Descent, and Landing Challenges[J]. Journal of Spacecraft and Rockets, 2007, 44(2): 310-323.

[3] 王海濤, 郭鵬, 程文科, 等. 傘艙系統在火星環境中的運動特性分析[J]. 航天返回與遙感, 2010, 31(1): 28-33.

WANG Haitao, GUO Peng, CHENG Wenke, et al. Analysis of Motion Characteristic of Parachute-vehicle System in Mars Environment[J]. Spacecraft Recovery & Remote Sensing, 2010, 31(1): 28-33. (in Chinese)

[4] 魯媛媛, 榮偉, 吳世通. “火星探路者”艙傘系統動力學特性仿真研究[J]. 中國空間科學技術, 2018, 38(2): 67-74.

LU Yuanyuan, RONG Wei, WU Shitong. Dynamic Characteristics Simulation Analysis for Parachute-vehicle System of Mars Pathfinder[J]. Chinese Space Science and Technology, 2018, 38(2): 67-74. (in Chinese)

[5] BALARAM J, AUSTIN R, BANERJEE P, et al. DSENDS a High-fidelity Dynamics and Spacecraft Simulator for Entry, Descent and Surface Landing[C]//IEEE Aerospace Conference, March 9-16, 2002, Big Sky, MT, USA.

[6] CAMERON J M, JAIN A, BURKHART P D, et al. DSENDS: Multi-mission Flight Dynamics Simulator for NASA Missions[C]//American Institute of Aeronautics and Astronautics AIAA SPACE, July 25-27, 2016, Long Beach, California, USA.

[7] 沈剛輝. 火星探測傘降段精確建模與EDL過程Unity3D仿真實現[D]. 北京: 北京理工大學, 2016.

SHEN Ganghui. The Precise Model of Parachute and Unity3D Simulation in Mars EDL[D]. Beijing: Beijing Institute of Thchnology, 2016. (in Chinese)

[8] 張青斌, 豐志偉, 馬洋, 等. 火星EDL過程動力學建模與仿真[J]. 宇航學報, 2017, 38(5): 443-450.

ZHANG Qingbin, FENG Zhiwei, MA Yang, et al. Modeling and Simulation of Mars EDL Process[J]. Journal of Astronautics, 2017, 38(5): 443-450. (in Chinese)

[9] VINH N X. Optimal Trajectories in Atmospheric Flight[M]. New York: Elsevier Scientific Publishing Company, 1981.

[10] 趙漢元. 飛行器再入動力學和制導[M]. 長沙: 國防科技大學出版社, 1997.

ZHAO Hanyuan. Reentry Dynamic and Guidance of Aircraft[M]. Changsha: National University of Defense Technology Press, 1997. (in Chinese)

[11] WOLF D . A Simplified Dynamics Model of Parachute Inflation[C]//American Institute of Aeronautics and Astronautics 4th Aerodynamic Deceleration Systems Conference , May 21-23, 1973, Palm Springs, CA, USA.

[12] 張青斌, 唐乾剛, 彭勇, 等. 飛船返回艙降落傘系統動力學[M]. 北京: 國防工業出版社, 2013.

ZHANG Qingbin, TANG Qiangang, PENG Yong, et al. Dynamics of Parachute-capsule Recovery System[M]. Beijing: National Defense Industry Press, 2013. (in Chinese)

[13] 熊菁, 秦小波, 程文科. 降落傘系統附加質量的研究[J]. 中國空間科學技術, 2002, 7(4): 32-38.

XIONG Jing, QIN Xiaobo, CHENG Wenke. The Added Mass Research in Parachute System[J]. Chinese Space Science and Technology, 2002, 7(4): 32-38. (in Chinese)

[14] 劉濤, 許望晶, 滕海山, 等. 射傘火箭開傘過程動力學仿真與分析[J]. 航天返回與遙感, 2020, 41(4): 45-54.

LIU Tao, XU Wangjing, TENG Haishan, et al. Dynamic Simulation and Analysis of Parachute Deployment Using a Propelling Rocket[J]. Spacecraft Recovery & Remote Sensing, 2020, 41(4): 45-54. (in Chinese)

[15] 宋旭民, 彭勇, 程文科, 等. 一種新的吊掛系統建模方法[J]. 中國空間科學技術, 2005, 10(5): 57-61.

SONG Xumin, PENG Yong, CHENG Wenke, et al. A Novel Modeling Method of Suspension System[J]. Chinese Space Science and Technology, 2005, 10(5): 57-61. (in Chinese)

[16] 宋旭民, 秦子增, 程文科, 等. 具有倒“Y”型吊掛的降落傘系統動力學建模[J]. 國防科技大學學報, 2005, 27(5): 103-106.

SONG Xumin, QIN Zizeng, CHENG Wenke, et al. The Dynamic Model of a Parachute System with the Inverted “Y” Suspension[J]. Journal of National University of Defense Technology, 2005, 27(5): 103-106. (in Chinese)

[17] SPENCER D A, BLANCHARD R C, BRAUN R D. Mars Pathfinder Entry, Descent, and Landing Reconstruction[J]. Journal of Spacecraft and Rockets, 1999, 36(3): 357-367.

[18] BRAUN R D, POWELL R W, ENGELUND W C, et al. Mars Pathfinder Six Degree-of-freedom Entry Analysis[J]. Journal of Spacecraft and Rockets, 1995, 32(6): 670-676.

[19] DESAI P N, SCHOENENBERGER M, CHEATWOOD F M. Mars Exploration Rover Six-degree-of-freedom Entry Trajectory Analysis[J]. Spacecraft and Rockets, 2003, 43(5): 1-5.

[20] MOSS J N, BLANCHARD R C, WILMOTH R G, et al. Mars Pathfinder Rarefied Aerodynamics: Computations and Measurements[J]. Journal of Spacecraft and Rockets, 1999, 36(3): 330-339.

[21] WITKOWSKI A. Mars Pathfinder Parachute System Performance[C]//15th Aerodynamic Decelerator Systems Technology Conference, June 8-11, 1999, Toulouse, France.

[22] EDWARD J F. System Design Overview of the Mars Pathfinder Parachute Decelerator Sub System[C]//14th AIAA Aerodynamic Decelerator System Technology Conference, June3-5, 1997, San Francisco, CA, USA.

[23] DAVID A S, ROBERT D B. Mars Pathfinder Atmospheric Entry: Trajectory Design and Dispersion Analysis[J]. Journal of Spacecraft and Rockets, 1996, 33(5): 670-676.

Mars EDL Process Dynamic Simulation Platform

TANG Mingzhang1WANG Qi1HUANG Wei1QIAN Kai2

（1 Beijing Institute of Space Mechanics & Electricity, Beijing 100094, China）（2 Beijing Dongfang Chuangda Software Co., Ltd, Beijing 100089, China）

Mars EDL process is one of the most important stages in Mars exploration mission, whether the landing on Mars’s surface is successful, determines the success or failure of the mission. This paper concerns about the major dynamic phases in EDL process: entry, parachute deployment, parachute inflation and parachute-payload steady descent. The dynamic models are developed to describe the dynamic phases, which include entry 6-DOF dynamic model, parachute deployment and inflation dynamic model, multi-body dynamic model with elastic coupling. Correspondingly the computational modules based on dynamics are developed. The modules are integrated in sequence of EDL process, the data transfer between these modules is implemented. Finally, the integrated software system is used to simulate the dynamic of Mars Pathfinder’s EDL whole process. The results indicate that the software can be applied to effectively predict and analyze the dynamics of Mars EDL process.

Mars exploration; Entry, Descent, Landing (EDL) dynamics; simulation; software design

V11；V445

A

1009-8518(2021)03-0032-09

10.3969/j.issn.1009-8518.2021.03.004

2021-03-05

國家重大科技專項工程

唐明章, 王奇, 黃偉, 等. 火星EDL過程動力學仿真平臺[J]. 航天返回與遙感, 2021, 42(3): 32-40. TANG Mingzhang, WANG Qi, HUANG Wei, et al. Mars EDL Process Dynamic Simulation Platform[J]. Spacecraft Recovery & Remote Sensing, 2021, 42(3): 32-40. (in Chinese)

唐明章，男，1981年生，2006年獲中國空間技術研究院飛行器設計專業工學碩士學位，高級工程師。研究方向為航天器進入，減速及著陸技術研究。E-mail：tangmz.cast@foxmail.com。

（編輯：陳艷霞）