連亮張紅英李方童明波

（1南京航空航天大學航空宇航學院，南京 210016）

（2空降兵研究所，孝感 432100）

1 引言

火星探測器要成功登陸火星，需要一套高精度進入、減速和著陸系統(tǒng) （Entry,Descent,and Landing,EDL），而降落傘減速系統(tǒng)（Parachute Deceleration System,PDS）是EDL系統(tǒng)的重要組成部分，給探測器提供非常有效的大氣阻力，使探測器的下降速度從超聲速減到亞聲速，為軟著陸創(chuàng)造有利條件。由于火星大氣密度非常稀薄，降落傘的開傘條件和工作性能都與地球上的降落傘有所不同。因此，從火星著陸探測計劃開始，降落傘的研究就一直是整個探測計劃的工作重點和難點之一。

從“海盜號”到“火星探測漫游者”的幾次火星探測著陸任務(wù)中，降落傘均采用了盤縫帶傘[1-3]。如“火星探路者”(MarsPathfinder,MPF)使用的盤縫帶傘是在“海盜號”的基礎(chǔ)上，根據(jù)進入質(zhì)量、開傘高度、動壓等技術(shù)要求作了相應改進[4-6]。參數(shù)與“海盜號”有顯著不同，主要是帶寬明顯增加，達到名義直徑的0.234倍；使阻力系數(shù)從0.67減小到0.41，而降落傘在超聲速時的穩(wěn)定性卻得到明顯改善[7-8]。因此降落傘結(jié)構(gòu)參數(shù)對開傘性能有著明顯的影響，本文以“火星探路者”所用傘為基礎(chǔ)，對其結(jié)構(gòu)作相應改變后進行充氣過程仿真計算，研究結(jié)構(gòu)參數(shù)對降落傘充氣性能的影響。希望能夠為我國未來的火星探測研究工作提供參考。

2 數(shù)學模型

進行數(shù)值模擬計算時，一般把降落傘充氣過程分為兩個階段：初始充氣階段和主充氣階段。初始充氣階段是指氣流從傘衣底部沖到傘衣頂部的階段，主充氣階段是從氣流沖到傘衣頂部開始一直到傘衣第一次完全充滿。初始充氣階段時間非常短，通常情況下對開傘動載和降落傘的充氣性能影響不大，本文只對主充氣階段進行計算。

由于盤縫帶傘是軸對稱傘，傘衣由一定數(shù)量的相同傘衣幅彼此縫合連接而成，相鄰傘衣幅之間由加強帶相連。在對稱充氣情況下，選取徑向?qū)ΨQ的2個加強帶來表示傘衣充氣過程的形狀變化，將研究對象（1根加強帶及其兩邊的2個半幅傘衣）離散為一系列用阻尼彈簧連接的質(zhì)點，分布在加強帶上，各質(zhì)點之間用傘繩、加強帶和傘衣連接，這樣質(zhì)點在彈性力、阻尼力和氣動力（忽略質(zhì)點重力）的作用下運動，傘衣各質(zhì)點的結(jié)構(gòu)動力學方程參見參考文獻[9]。

降落傘傘衣的厚度很小，遠小于傘衣的幾何尺度。因此從流體力學的角度來考慮，可以把傘衣看成是薄膜材料，其厚度對流場本質(zhì)不會產(chǎn)生影響，在對降落傘進行流體力學數(shù)值模擬時可以忽略傘衣厚度的影響。采用準定常假設(shè)，把充氣過程中的降落傘視為剛體，通過對不同充氣時刻的不同傘衣外形進行Navier-Stokes(N-S)方程求解來模擬降落傘的流場特性?？刂品匠滩捎每蓧嚎s雷諾平均N-S方程，其積分形式可以寫成：

式中 t是時間；S是面積； 是控制體體積； 是控制體邊界； 是守恒變量； 和 分別是對流通量和粘性通量。

在這基礎(chǔ)上對結(jié)構(gòu)模型和流場模型進行耦合計算：先根據(jù)初始數(shù)據(jù)，在每一時間步長開始時，降落傘到達一個新位置，將傘面附近的流場網(wǎng)格點移動至傘面上；通過修改這些點相應的動量方程，更新源項，利用CFD程序計算得出交界處質(zhì)點的壓差，將計算結(jié)果傳給結(jié)構(gòu)動力學方程，進行表面力計算和進行傘面變形運動計算，得到下一形狀；依據(jù)下一形狀，流場計算程序進入下一個時間步長，開始新一輪的計算。

3 數(shù)學模型驗證

為了校驗仿真結(jié)果的可信度，首先以“火星探路者”降落傘為例，進行火星探測用降落傘的充氣過程仿真計算。根據(jù)查找的參考文獻[8,10]所對應的計算狀態(tài)進行數(shù)值模擬，計算參數(shù)如表1所示，將計算得到的充滿時間和開傘過載曲線與試驗結(jié)果進行對比，以驗證仿真程序在火星大氣環(huán)境下的可靠性。

?

[8]中，試驗得到的充滿時間是0.53s，開傘動載曲線如圖1中的虛線所示。本文通過仿真計算，得到“火星探路者”探測器所用的盤縫帶傘首次達到穩(wěn)降時，47m2的投影面積所用的時間為0.54s，與實測充滿時間0.53s誤差僅0.01s，計算得到的充氣距離為198m。在圖1的開傘過載曲線中，實線為計算結(jié)果，虛線為試驗結(jié)果。仿真數(shù)據(jù)顯示：降落傘充氣過程中最大動載出現(xiàn)在0.58s，峰值3950N，則開傘過載為6.89gn，開傘峰值出現(xiàn)在傘衣完全充滿后。從圖中可知無論是從計算得到的最大過載，還是最大過載出現(xiàn)的時間點，計算與實測曲線均比較一致。說明該數(shù)值模擬方法能較準確地模擬火星探測用降落傘的開傘過程，可利用本文的仿真手段開展進一步的分析。

圖1 過載與充氣時間曲線Fig.1 Load factor with in flat ion time

4 計算結(jié)果分析

表2給出了“火星探路者”與“海盜號”降落傘的結(jié)構(gòu)參數(shù)。本文以“火星探路者”所用傘為基礎(chǔ)，對其結(jié)構(gòu)作5種變化后進行仿真計算，研究結(jié)構(gòu)參數(shù)對降落傘性能的影響。修改后5種降落傘（A，B，C，D，E）的結(jié)構(gòu)參數(shù)如表3所示，計算參數(shù)如表4所示，其中A傘即MPF傘。

表2 “火星探路者”與“海盜號”傘結(jié)構(gòu)參數(shù)Tab. 2 Structure parameters of parachutes used for MPF and VINKING

表3 修改后的降落傘結(jié)構(gòu)參數(shù)Tab. 3 Structure parameters of parachutes after modification

表4 修改后的降落傘計算參數(shù)Tab. 4 Calculation parameters of parachutes after modification

修改后的降落傘計算結(jié)果值如表5所示。表中，V為計算結(jié)束時的速度；Tm為傘衣投影面積第一次達到充滿面積47m2的時間；T1為開傘載荷達到最大值的時間；Cs為計算結(jié)束時的阻力系數(shù)；n為開傘最大過載。

表5 修改后的降落傘計算結(jié)果Tab. 5 Calculation results of parachutes after modification

序號A，B和C只是帶寬不同，帶寬由小到大的順序是B，A，C。從表5可看出，對應的傘衣充滿后的阻力系數(shù)分別為0.587、0.520、0.435，即隨著帶寬增長阻力系數(shù)減小。投影面積隨充氣時間的變化圖如圖2所示，由圖可知，投影面積的波動隨帶寬增長變得緩和，可見開傘過程中傘衣的顫動減小。

圖2 A，B，C三傘投影面積隨充氣時間變化圖Fig.2 Projected area of parachutes A,B and Cw ith inflation time

圖3 B傘開傘過程傘衣外形變化圖 Fig.3 Configuration during opening of parachute B

圖4 A傘開傘過程傘衣外形變化圖Fig.4 Configuration during opening of parachute A

圖5 C傘開傘過程傘衣外形變化圖Fig.5 Configuration during opening of parachute C

圖3 、圖4和圖5分別為B，A，C開傘過程的傘衣外形變化圖，可以看出B傘開傘過程中外形波動最為明顯，傘衣底邊有嚴重的擺動，而傘頂孔則存在忽上忽下的脈動，這種現(xiàn)象就是超聲速開傘時的喘振現(xiàn)象。該現(xiàn)象在A傘上得到緩解，而C傘開傘過程中傘衣外形則比較光滑，喘振現(xiàn)象已經(jīng)變得不明顯。

圖6、圖7和圖8分別是B，A，C三傘接近充滿狀態(tài)時的某一時刻傘衣周圍流場圖譜。從圖6可以看出B傘由于縫寬最小，使得傘衣頂部的氣流分離最劇烈，漩渦尺寸最大，也正是由于傘頂尾流的作用，開傘過程中傘衣的波動最為明顯。從圖7和圖8可看出，隨著帶寬增加，傘頂分離減弱，尾流區(qū)的漩渦尺寸也減小。由此也可得出，超聲速開傘時，降落傘充氣過程中傘衣的喘振現(xiàn)象隨著帶寬增加而減弱，這對盤縫帶傘的開傘安全性和可靠性有利。

圖6 B傘充滿時刻傘衣周圍流場圖譜Fig.6 Stream linesof inflated parachute B

圖7 A傘充滿時刻傘衣周圍流場圖譜Fig.7 Stream linesof inflated parachute A

圖8 C傘充滿時刻傘衣周圍流場圖譜Fig.8 Stream linesof inflated parachute C

圖9 A,D、E三傘投影面積隨充氣時間變化圖Fig.9 Projected area of parachutesA,D and Ew ith inflation time

序號A，D和E只是縫寬不同，A傘的縫寬最小，E的縫寬最大，它們的阻力系數(shù)分別為0.52、0.504、0.496，即隨縫寬增長阻力系數(shù)減?。ㄓ捎谶x取的縫寬相差不大，計算結(jié)果變化規(guī)律不明顯）。若從盤縫帶傘設(shè)計的角度考慮，相對而言D傘的縫寬0.54m比較合適，因為它在三者中過載最小，為11gn。另外，從圖9的A，D，E三傘投影面積隨充氣時間變化圖中可看出，D傘的投影面積變化也比較緩和。

5 結(jié)束語

本文以火星探測用傘（盤縫帶傘）為研究對象，研究其結(jié)構(gòu)參數(shù)對降落傘充氣性能的影響，為盤縫帶傘的優(yōu)化設(shè)計提供參考。通過計算得到如下結(jié)論：

1）隨著帶寬增長，盤縫帶傘阻力系數(shù)減小。

2）超聲速開傘時，降落傘充氣過程中傘衣的喘振現(xiàn)象隨著帶寬的增加而減弱，投影面積的波動隨帶寬增長變得緩和，可見開傘過程中傘衣的顫動減小。這對盤縫帶傘的開傘安全性和可靠性有利。

3）隨著縫寬增長，盤縫帶傘阻力系數(shù)減小。

參考文獻（References）

[1]Cruz JR,Stephen L J.Aerodynamic Decelerators for Planetary Exploration:Past,Present,and Future[R].AIAA 2006-6792,NASA Langley Research Center,2006.

[2]Jinwook L,Reuben RR,StephenM R.Fluid-structure AnalysisofaClamped Ballute in Titan's Atmosphere[C]//AIAA 2007-4308.M iami,USA:AIAA Computational Fluid DynamicsConference,2007.

[3]Edward JF.System Design Overview of theMarsPathfinder Parachute Decelerator Subsystem[R].AIAA-97-1511,1997.

[4]A llenW.MarsExploration Rover Parachute Decelerator System Program Overview[R].AIAA 2003-2100,2003.

[5]Charles H W.Inflation and Performance of Three Parachute Con Figurations from Supersonic Flight Testsn a Low-density Environment[R].NASA TN D-5296,1969.

[6]榮偉,陳旭.火星探測用降落傘研制試驗簡介[J].航天返回與遙感,2007,28(1):12-17.RONGWei,CHEN Xu.Introduction to Testsof Parachute Development forMars Exploration[J].Spacecraft Recovery amp;Remote Sensing,2007,28(1):12-17.(in Chinese)

[7]Anita S.An Overview of theMars Science Laboratory Parachute Decelerator System[C].//Aerospace Conference,2007 IEEE.

[8]A llenW.MarsPathfinder Parachute System Performance[R].AIAA 99-1701,11-19.

[9]張紅英,劉衛(wèi)華,秦福德,等.降落傘充氣過程中傘衣外形及流場變化研究[J].空氣動力學學報,2011,29(3):288-294.ZHANG Hongying,LIUWeihua,QIN fude,etal.Study on Canopy Shape and Flow Field During Parachute Inflation Process[J].Acta Aerodynamica Sinica,2011,29(3):288-294.(in Chinese)

[10]Norio A,Kensaku S.Inflation Process and Free Oscillation of Flexible Parachute-like Body[R].AIAA 2005-1610,2005.