盤(pán)–縫–帶傘超聲速充氣過(guò)程仿真研究

王祁 曹義華

?

盤(pán)–縫–帶傘超聲速充氣過(guò)程仿真研究

王祁 曹義華

（北京航空航天大學(xué)，北京 100083）

隨著火星探測(cè)著陸研究的發(fā)展，對(duì)其主要的減速方式——盤(pán)–縫–帶傘應(yīng)用的研究也日趨火熱。而在火星探測(cè)中，對(duì)傘衣充氣過(guò)程的分析又是最為復(fù)雜的難題之一。文章介紹了常用于降落傘流場(chǎng)求解的數(shù)值模擬方法，并對(duì)它們的利弊做了簡(jiǎn)要的描述。文章使用守恒元/解元方法對(duì)盤(pán)–縫–帶傘在超聲速條件下，從收束狀態(tài)開(kāi)始的充氣過(guò)程進(jìn)行了數(shù)值模擬。首先對(duì)前置體進(jìn)行單獨(dú)仿真，之后在不同的馬赫數(shù)下通過(guò)對(duì)降落傘系統(tǒng)的整體仿真完成對(duì)降落傘開(kāi)傘過(guò)程的模擬。仿真數(shù)據(jù)與實(shí)驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比，獲得了較為一致的結(jié)果，驗(yàn)證了方法的可行性。另外，文章對(duì)兩種工作狀態(tài)下的充氣過(guò)程進(jìn)行了對(duì)比，分析了盤(pán)–縫–帶傘在兩者中充氣過(guò)程的異同，以及前置體阻力貢獻(xiàn)隨著馬赫數(shù)的變化。最后通過(guò)開(kāi)傘云圖對(duì)充氣過(guò)程中的開(kāi)傘失敗現(xiàn)象進(jìn)行了分析。研究結(jié)果可為火星探測(cè)提供參考。

盤(pán)–縫–帶傘 流固耦合 守恒元/解元方法 超聲速 探測(cè)著陸

0 引言

在常見(jiàn)的火星探測(cè)器著陸減速系統(tǒng)中，盤(pán)–縫–帶傘是典型的超聲速、低密度降落傘傘型，在超聲速、低密度的工作環(huán)境中穩(wěn)定性好，充氣性能優(yōu)良。國(guó)外已成功著陸的Viking1、Viking2、MPF等火星探測(cè)器均采用了盤(pán)–縫–帶傘[1-2]。降落傘的充氣過(guò)程是整個(gè)工作過(guò)程中最復(fù)雜的一個(gè)階段，充氣過(guò)程處于一個(gè)快速的大變形狀態(tài)，涉及到結(jié)構(gòu)動(dòng)力學(xué)與流體力學(xué)耦合的相關(guān)問(wèn)題[3]，這些在文獻(xiàn)[4-7]的試驗(yàn)結(jié)果中均有體現(xiàn)。而對(duì)于火星降落傘來(lái)說(shuō)，低密度高馬赫數(shù)的工作環(huán)境又很難通過(guò)試驗(yàn)完成，因此對(duì)這一領(lǐng)域進(jìn)行仿真研究顯得尤為重要[8-10]。

降落傘流場(chǎng)的數(shù)值模擬主要有兩種方法。一種是渦元法。它在建立流場(chǎng)方程的基礎(chǔ)上，通過(guò)傘衣型面及尾渦面上的邊界條件，求解出繞傘衣流場(chǎng)的數(shù)值解。渦元法最大的好處是容易收斂，在軸對(duì)稱傘衣上的計(jì)算結(jié)果非常好，但這種方法通常用在不可壓流場(chǎng)的計(jì)算中。另一種是基于網(wǎng)格的流固耦合方法。常用的方法中，又可以分為傳統(tǒng)的松散耦合方法（詳見(jiàn)文獻(xiàn)[11]，由于工作量大且精度無(wú)法保證，不建議采用）、任意拉格朗日歐拉（Arbitrary Lagrange Euler，ALE）方法和守恒元/解元（Conservation Element and Solution Element Method，CE/SE）方法。

ALE方法最初用于數(shù)值模擬流體力學(xué)的有限差分方法中，兼有Lagrage方法和Euler方法的特長(zhǎng)，使得網(wǎng)格不至于出現(xiàn)嚴(yán)重的畸變[12]。國(guó)外的文獻(xiàn)中無(wú)論是理論研究還是工程研究都已經(jīng)頻繁的使用過(guò)ALE方法[13-15]。這種方法在處理大變形問(wèn)題時(shí)非常的有利，但是由于在超聲速條件下，ALE方法計(jì)算耗時(shí)長(zhǎng)，難收斂，往往得不到正確的結(jié)果，還可能出現(xiàn)計(jì)算終止，所以其更適用于計(jì)算亞聲速[11,16]。

CE/SE方法，又稱為守恒元/解元方法，是一種全新的守恒方程計(jì)算方法，國(guó)外文獻(xiàn)對(duì)于這一方法的理論研究較為成熟[17-19]。它可以與浸潤(rùn)邊界法[20]一起使用，將時(shí)間和空間統(tǒng)一起來(lái)同等對(duì)待，構(gòu)造比較簡(jiǎn)單。除了簡(jiǎn)單的泰勒展開(kāi)之外，CE/SE方法不需要采用其他的數(shù)值方法，尤其是不用通過(guò)其他的特征分析數(shù)值方法（如Riemann求解器）來(lái)捕捉激波、抑制振蕩等[21]。

在研究降落傘領(lǐng)域，基于ALE方法的流固耦合數(shù)值仿真已經(jīng)比較成熟，而CE/SE方法的工程化應(yīng)用較晚，針對(duì)這一方法在降落傘領(lǐng)域的研究并不多。國(guó)內(nèi)對(duì)于新方法的研究則還處于起步階段，包括對(duì)降落傘充氣過(guò)程的簡(jiǎn)化以及充氣馬赫數(shù)的選取都不能較為準(zhǔn)確的反應(yīng)盤(pán)–縫–帶傘的真實(shí)充氣過(guò)程。因此，對(duì)盤(pán)–縫–帶傘超聲速充氣過(guò)程的仿真研究具有一定的理論意義和實(shí)踐價(jià)值。

1 仿真過(guò)程

1.1 仿真模型

盤(pán)–縫–帶傘是開(kāi)縫傘的一種，如圖1（a）所示，主要由傘衣、傘繩和吊帶等組成。其中傘衣由平面圓形“盤(pán)”和圓筒型“帶”組成，中間由縫隙將兩者垂直分開(kāi)。傘衣幅呈扇形，盤(pán)的中間有通氣孔。在盤(pán)和帶的周向徑向均有加強(qiáng)條帶，以防止傘衣過(guò)度張開(kāi)[1]。圖1（b）為傘衣的網(wǎng)格示意圖。為了能夠和文獻(xiàn)[12]中的實(shí)驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比，本文選用了文獻(xiàn)中相同的傘型。這是一種縮比的維京傘，相關(guān)的數(shù)據(jù)在表1中列出。

仿真模型中使用的前置體為具有70°鈍角的減速體，整體采用結(jié)構(gòu)網(wǎng)格進(jìn)行劃分。相關(guān)的數(shù)據(jù)體現(xiàn)在圖2中，圖中半徑1=0.085m，2=0.852m。

如圖3所示，選用圓柱形流場(chǎng)，這樣既可以很好地捕捉傘衣的尾流，又可以減少流域邊界對(duì)仿真結(jié)果的影響。由于降落傘帶有前置體，所以流場(chǎng)的周向和徑向都進(jìn)行了加密。為了簡(jiǎn)化模型，外流域選用了相對(duì)稀疏的網(wǎng)格，但是相鄰網(wǎng)格的體積增長(zhǎng)控制在30%以內(nèi)。CE/SE方法所使用的網(wǎng)格和ALE方法比較接近，均使用不兼容的網(wǎng)格，即流場(chǎng)和降落傘為兩套互無(wú)搭接的網(wǎng)格。所不同的是CE/SE方法可以在初始運(yùn)算中將氣流靜壓帶入到流場(chǎng)中，而ALE方法則需要外包一層提供流場(chǎng)靜壓網(wǎng)格，同時(shí)來(lái)流也需要單獨(dú)設(shè)置，所以兩者的流場(chǎng)不甚相同。本文中接觸邊界的力傳遞采用了浸潤(rùn)邊界法，所以也不需要生成特殊的邊界層網(wǎng)格。實(shí)際上在大變形情況下采用傳統(tǒng)方法處理邊界層更容易產(chǎn)生網(wǎng)格畸變。

表1 盤(pán)–縫–帶傘模型參數(shù)

Tab.1 Disk–gap–band parachute model parameters

圖2 前置體模型參數(shù)

圖3 流場(chǎng)模型

1.2 流場(chǎng)仿真條件

流場(chǎng)的仿真條件來(lái)自于文獻(xiàn)[12]，已整理在表2中。因?yàn)槭腔鹦堑拇髿鈼l件，所以密度要比地球低很多，這種低密度高馬赫數(shù)的開(kāi)傘條件也使仿真結(jié)果更加接近真實(shí)充氣環(huán)境。

表2 流場(chǎng)仿真條件

Tab.2 Simulation Condition of fluid field

2 仿真結(jié)果的試驗(yàn)驗(yàn)證

由于降落傘的阻力是在不停的變化中，而且在高馬赫數(shù)條件下降落傘充氣過(guò)程中的阻力劇烈振蕩，仿真的準(zhǔn)確性很難通過(guò)阻力衡量。所以仿真分為兩部分，首先針對(duì)沒(méi)有傘衣作用的前置體進(jìn)行仿真，通過(guò)仿真結(jié)果與實(shí)驗(yàn)數(shù)據(jù)的對(duì)比來(lái)驗(yàn)證仿真對(duì)阻力捕捉的準(zhǔn)確性。在此基礎(chǔ)上再對(duì)傘衣的充氣過(guò)程進(jìn)行仿真。由于馬赫數(shù)=2.8時(shí)，缺少傘衣的試驗(yàn)數(shù)據(jù)，將選用=1.2，1.6，2.0，2.4作為數(shù)據(jù)驗(yàn)證。試驗(yàn)數(shù)據(jù)來(lái)源于文獻(xiàn)[12]。如圖4所示，在前置體的仿真結(jié)果中，實(shí)驗(yàn)數(shù)據(jù)和仿真結(jié)果擬合的非常接近，這說(shuō)明CE/SE方法和浸潤(rùn)邊界法在獲取靜態(tài)物體的阻力系數(shù)時(shí)是有效的，對(duì)邊界層的處理也是符合物理實(shí)際的。

仿真結(jié)果顯示，傘衣的阻力系數(shù)也和實(shí)驗(yàn)結(jié)果很接近，如圖5所示。由于降落傘在仿真過(guò)程中未考慮傘衣透氣性，所以阻力系數(shù)會(huì)在部分時(shí)刻大于實(shí)驗(yàn)值。

圖4 前置體仿真結(jié)果與實(shí)驗(yàn)數(shù)據(jù)對(duì)比

圖5 傘衣仿真結(jié)果與實(shí)驗(yàn)數(shù)據(jù)對(duì)比

需要說(shuō)明的是，由于“呼吸”現(xiàn)象（傘衣第一次充滿后，傘衣繼續(xù)向外擴(kuò)張充氣，導(dǎo)致降落傘阻力面積峰值大于設(shè)計(jì)的全充滿阻力面積，使開(kāi)傘充氣過(guò)程中最大開(kāi)傘力預(yù)測(cè)精度降低[22]），傘衣的阻力是不穩(wěn)定的，而耗費(fèi)大量的時(shí)間求取傘衣穩(wěn)定狀態(tài)下的阻力系數(shù)也不是本文的研究重點(diǎn)，所以這里的阻力系數(shù)近似選取是傘衣第一次完全張滿的狀態(tài)，這里的結(jié)果僅作為仿真有效性的參考。

3 仿真結(jié)果分析

3.1 不同馬赫數(shù)下的充氣過(guò)程

由于模型的傘衣收束狀態(tài)是完美的，所以盤(pán)–縫–帶傘的充氣會(huì)首先在頂盤(pán)處開(kāi)始。呈錐形收束的傘盤(pán)通過(guò)充氣膨脹，繼而帶動(dòng)傘帶擴(kuò)張，因此在充氣過(guò)程中傘帶的投影面積始終小于傘盤(pán)。同時(shí)，通過(guò)仿真結(jié)果可以看出，降落傘在第一次充氣之后，其投影面積會(huì)在一定范圍內(nèi)波動(dòng)，這種“呼吸”現(xiàn)象在真實(shí)的試驗(yàn)過(guò)程中多次出現(xiàn)[22]，可見(jiàn)仿真結(jié)果比較真實(shí)的再現(xiàn)了降落傘的實(shí)際充氣過(guò)程。盤(pán)–縫–帶傘主要應(yīng)用于超聲速環(huán)境下，而且在=2.7以下性能是比較好的，在達(dá)到=2.7時(shí)會(huì)呈現(xiàn)較大的振動(dòng)現(xiàn)象[23]?；谝陨显?，在本小結(jié)中主要分析=1.2和=2.8兩種速度下的充氣過(guò)程，其中=1.2表征了降落傘的正常工作狀態(tài)，=2.8表征了降落傘的非正常工作狀態(tài)。圖6所示即為盤(pán)–縫–帶傘在=1.2條件下的充氣過(guò)程?？梢钥闯銎湓谡９ぷ鳡顟B(tài)下可以很好的完成充氣，其傘衣振蕩幅度也比較小。

圖7所示為=1.2條件下，盤(pán)–縫–帶傘的開(kāi)傘充氣過(guò)程速度云圖。在時(shí)間=0.004s，即充氣的初始時(shí)刻，氣流從傘衣下部入口進(jìn)入，在傘衣內(nèi)壁阻滯，使得傘衣開(kāi)始變形。之后隨著進(jìn)入氣體的增多，傘衣變形逐漸加大，在=0.3s高速氣流到達(dá)頂孔并溢出。由于頂孔的阻滯作用，進(jìn)入的氣流大于頂孔流出氣流，傘衣開(kāi)始快速膨脹，多余的氣體迅速充滿傘衣內(nèi)腔，開(kāi)始從傘縫區(qū)域外流。最終進(jìn)入的氣體和流出的氣體達(dá)到平衡狀態(tài)，傘衣進(jìn)入完全張滿的穩(wěn)定狀態(tài)。

盤(pán)–縫–帶傘在不同工作環(huán)境下的性能表現(xiàn)是不同的。由于其有能夠正常工作的馬赫數(shù)區(qū)間，故而在高馬赫數(shù)條件下降落傘會(huì)發(fā)生劇烈的抖動(dòng)。這些抖動(dòng)對(duì)降落傘的影響是不確定的，降落傘有可能會(huì)最終給張開(kāi)并穩(wěn)定，但是也有可能因?yàn)閯×叶秳?dòng)而發(fā)生開(kāi)傘失敗，如圖8所示。根據(jù)仿真結(jié)果來(lái)看，=1.2時(shí)，降落傘開(kāi)傘成功，形狀和阻力最終都趨于穩(wěn)定，而=2.8時(shí)則會(huì)發(fā)生劇烈振蕩，最終導(dǎo)致降落傘的開(kāi)傘失敗。

3.2 不同馬赫數(shù)下的阻力變化過(guò)程

阻力變化和降落傘的充氣過(guò)程是緊密聯(lián)系的，降落傘過(guò)度擴(kuò)張的時(shí)間點(diǎn)也恰好是阻力值達(dá)到最大的時(shí)刻。同時(shí)，=2.8條件下的阻力變化要遠(yuǎn)遠(yuǎn)劇烈于=1.2，可見(jiàn)馬赫數(shù)對(duì)降落傘的充氣表現(xiàn)至關(guān)重要。同樣的傘型在不同的馬赫數(shù)下也會(huì)有不同的阻力表現(xiàn)。

圖9所示為=1.2條件下的阻力變化過(guò)程。可以明顯的看出，雖然阻力在很小的范圍內(nèi)振蕩，但是整體是較為平滑的。雖然在圖中看不出比較明顯的穩(wěn)定跡象，但是根據(jù)開(kāi)傘過(guò)程動(dòng)畫(huà)可以看出，降落傘在=1s時(shí)已經(jīng)趨于穩(wěn)定。

圖10為=2.8條件下的降落傘阻力變化曲線。通過(guò)和=1.2的仿真結(jié)果對(duì)比可以看出，在相同時(shí)間內(nèi)阻力發(fā)生了更為劇烈的振蕩，而且在阻力的變化過(guò)程也是極其的不穩(wěn)定。這樣的結(jié)果也可以在開(kāi)傘動(dòng)畫(huà)中得到體現(xiàn)，傘衣的張開(kāi)過(guò)程很不穩(wěn)定，會(huì)出現(xiàn)劇烈的抖動(dòng)。這些抖動(dòng)使傘衣的形狀變化不再具有規(guī)律，過(guò)度的擴(kuò)張導(dǎo)致過(guò)度的收縮，進(jìn)而影響結(jié)構(gòu)透氣性，導(dǎo)致最終的開(kāi)傘失敗。

圖9 Ma=1.2時(shí)，傘衣阻力變化過(guò)程

圖10 Ma=2.8時(shí)，傘衣阻力變化過(guò)程

3.3 充氣過(guò)程中前置體阻力貢獻(xiàn)

圖11所示為=2.8時(shí)，降落傘充氣過(guò)程中的流場(chǎng)云圖?？梢钥闯銮爸皿w處對(duì)氣流的影響范圍要比傘衣大一些，但是對(duì)氣流的減速效果遠(yuǎn)小于傘衣，激波范圍也小于傘衣。同時(shí)，圖中可以看出前置體的存在會(huì)對(duì)流經(jīng)傘衣的氣流有減速效果，但是減速效果有限。

圖12所示為=2.8條件下前置體阻力的求解過(guò)程，從圖中可以發(fā)現(xiàn)，對(duì)于固定的物體來(lái)說(shuō)求解是十分穩(wěn)定和迅速的。通過(guò)將前置體的阻力占傘衣張滿狀態(tài)下的整體阻力的百分比求出即可得到前置體在降落傘整體減速效果中的貢獻(xiàn)量。可以看出，隨著馬赫數(shù)的增長(zhǎng)，前置體所占百分比逐漸增多，如表3所示。由于這里的阻力為降落傘完全張滿的阻力，所以在充氣過(guò)程中前置體的減速效果更加巨大，這一部分的初始減速直接決定了傘衣的開(kāi)傘和張滿狀態(tài)的工作馬赫數(shù)。

圖11 Ma=2.8時(shí)，降落傘整體速度云圖

圖12 Ma=2.8時(shí)，前置體阻力變化曲線

表3 前置體阻力貢獻(xiàn)比率

Tab.3 Contribution ratio of the probe drag

3.4 開(kāi)傘失敗的過(guò)程分析

通過(guò)對(duì)比=1.2（在3.1節(jié)中已通過(guò)圖7進(jìn)行過(guò)分析）和=2.8下開(kāi)傘過(guò)程云圖的異同，可以得出為什么高馬赫數(shù)下降落傘更容易開(kāi)傘失敗。在=2.8時(shí)，從圖10中可以看出，降落傘前后共經(jīng)歷兩次阻力峰值，分別在=0.243s左右與=0.584s左右。在=0.243s后降落傘經(jīng)歷多次“呼吸”現(xiàn)象，但整體是收縮狀態(tài)，在=0.5s時(shí)刻恢復(fù)擴(kuò)張狀態(tài)。而在=584s之后則經(jīng)歷急速收縮，直到充氣失敗。

圖13所示即為=2.8時(shí)，=0.243s前后降落傘一次“呼吸”現(xiàn)象的速度云圖，其中色溫越高表明氣流的速度越大?？梢悦黠@看出，在=0.199s時(shí)，降落傘內(nèi)部氣流是正常的充氣狀態(tài)，而=0.224s時(shí)開(kāi)始傘衣內(nèi)部回流氣體增多，此時(shí)降落傘開(kāi)始擴(kuò)張。但是回流的氣體并不能穩(wěn)定的維持在傘衣內(nèi)部，而是在傘衣內(nèi)部呈現(xiàn)不規(guī)則的變化。

圖14為與圖13同一時(shí)刻的密度云圖。綜合圖13和圖14，可以看出，在=0.224s時(shí)降落傘內(nèi)部的高壓區(qū)擴(kuò)大，由于內(nèi)部氣體受壓，不能及時(shí)從頂孔排出而產(chǎn)生阻滯，阻滯的氣流逐漸增多影響了激波的形狀，激波影響域加大，使得進(jìn)入降落傘的氣流減少，內(nèi)部氣壓降低，在彈性傘衣的回復(fù)力作用下，=0.224s之后傘衣開(kāi)始塌縮。而在=0.264s和=0.299s之間，阻滯在傘衣邊緣的回流氣體將整體的正激波擠壓分流，給了外部氣流流入的機(jī)會(huì)，在=0.329時(shí)刻傘衣已經(jīng)逐漸恢復(fù)擴(kuò)張狀態(tài)。

而在圖15中，=0.579s附近，降落傘首先經(jīng)歷了和圖13中一樣的過(guò)度擴(kuò)張，然而之后傘衣收縮嚴(yán)重，最終沒(méi)有能夠使傘衣重新張開(kāi)。從圖14中可以看出，傘衣的擴(kuò)張依靠的是內(nèi)部的阻滯氣體。我們知道，氣體的阻滯是因?yàn)閭阋挛樟藖?lái)流的動(dòng)量，轉(zhuǎn)化為了阻力；阻滯氣體量表征了降落傘傘衣吸收動(dòng)量的能力。所有降落傘擴(kuò)張的過(guò)程中，阻滯氣體都布滿了傘衣內(nèi)壁；在圖13中=0.329s時(shí)刻就存在這種現(xiàn)象，這也就是傘衣在=0.329s之后仍然可以恢復(fù)擴(kuò)張的原因。而在圖15中=0.689s，可以明顯的看出，阻滯氣體已經(jīng)不占有優(yōu)勢(shì)。這說(shuō)明此時(shí)傘衣并沒(méi)有將來(lái)流的動(dòng)量有效轉(zhuǎn)化為阻力。在=0.744s，被分離的激波重新合為一個(gè)，但激波形狀已經(jīng)改變，無(wú)法正常工作了。頂孔此時(shí)已經(jīng)偏離中軸線，所以在云圖中無(wú)法看到，實(shí)際上頂孔仍然在劇烈的噴涌氣體。綜合來(lái)看，應(yīng)該是過(guò)度擴(kuò)張導(dǎo)致的恢復(fù)力大于了阻力，致使收縮過(guò)度，氣體進(jìn)入量減少，最終導(dǎo)致了開(kāi)傘的失敗。

通過(guò)=1.2和=2.8時(shí)充氣過(guò)程的云圖對(duì)比，可以看出，成功充氣的盤(pán)–縫–帶傘內(nèi)部的低速氣體是穩(wěn)定的，同時(shí)速度也是一致的。而在高馬赫數(shù)時(shí)由于進(jìn)氣量大于出氣量，內(nèi)部氣流不穩(wěn)定，容易導(dǎo)致激波的破裂，對(duì)氣流的減速效果下降，傘衣整體的阻力也隨之下降。氣流的不穩(wěn)定還會(huì)使傘衣形狀發(fā)生劇烈變化，外型的改變是不可控的，傘衣在之后是否能回歸正常充氣則成了隨機(jī)事件。尤其對(duì)于盤(pán)–縫–帶傘來(lái)說(shuō)，傘盤(pán)和傘帶是由傘繩連接的，這就更加加劇了傘衣的變形幅度，傘帶的坍縮會(huì)導(dǎo)致嚴(yán)重的后果?？梢?jiàn)，對(duì)于工作于不同馬赫數(shù)下的盤(pán)–縫–帶傘，透氣性的設(shè)計(jì)至關(guān)重要。

4 結(jié)束語(yǔ)

本文通過(guò)CE/SE方法對(duì)盤(pán)–縫–帶傘在低密度超聲速條件下的充氣過(guò)程進(jìn)行了仿真。通過(guò)對(duì)仿真結(jié)果的分析得出了以下結(jié)論：

1）CE/SE方法在對(duì)超聲速流場(chǎng)的仿真求解過(guò)程中對(duì)計(jì)算資源的占用較小，耗時(shí)短，且與試驗(yàn)結(jié)果比較接近，是比較適合工程應(yīng)用的新方法。

2）盤(pán)–縫–帶傘可以在超聲速條件下完成充氣的過(guò)程，但是對(duì)馬赫數(shù)有要求。在較高的馬赫數(shù)條件下有開(kāi)傘失敗的風(fēng)險(xiǎn)。

3）前置體在盤(pán)–縫–帶傘的充氣過(guò)程中對(duì)傘衣的流場(chǎng)有一定影響，且會(huì)提供可觀的阻力貢獻(xiàn)。

4）傘衣內(nèi)部的氣體穩(wěn)定性對(duì)盤(pán)–縫–帶傘的整體穩(wěn)定性有重要影響。傘衣內(nèi)部氣流的不穩(wěn)定會(huì)使傘衣振蕩劇烈，可導(dǎo)致開(kāi)傘失敗的嚴(yán)重后果。為了消解這一不利影響，盤(pán)–縫–帶傘的結(jié)構(gòu)透氣性要謹(jǐn)慎選取。

[1] 高樹(shù)義, 戈嗣誠(chéng), 梁艷. 火星盤(pán)–縫–帶傘跨聲速風(fēng)洞試驗(yàn)研究[J]. 中國(guó)空間科學(xué)技術(shù), 2015(4): 69-75. GAO Shuyi, GE Sicheng, LIANG Yan. Research on Transonic Wing Tunnel Tests of Mars Disk-gap-band Parachutes[J]. Chinese Space Science and Technology, 2015(4): 69-75. (in Chinese)

[2] 于瑩瀟, 田佳林. 火星探測(cè)器降落傘系統(tǒng)綜述[J]. 航天返回與遙感, 2007, 28(4): 12-16. YU Yingxiao, TIAN Jialin. Mars Explorer’s Parachute System Overview[J]. Spacecraft Recovery & Remote Sensing, 2007, 28(4): 12-16. (in Chinese)

[3] 榮偉, 魯媛媛, 包進(jìn)進(jìn), 等. 火星探測(cè)器減速著陸過(guò)程中若干問(wèn)題的研究[J]. 南京航空航天大學(xué)學(xué)報(bào), 2016, 48(4): 445-453. RONG Wei, LU Yuanyuan, BAO Jinjin, et al. Several Problems Relative with Descent and Landing Process of Mars Explorer[J]. Nanjing University of Aeronautics and Astronautics, 2016, 48(4): 445-453. (in Chinese)

[4] REICHENAU D E A. Aerodynamic Characteristics of Disk-gap-band Parachutes in the Wake of Viking Entry Forebodies at Mach Numbers from 0.2 to 2.6[R]. Arnold Engineering Development Center Arnold AFB TN, AEDC-TR-72-78, 1972.

[5] CRUZ J R, MINECK R E, KELLER D F, et al. Wind Tunnel Testing of Various Disk-gap-band Parachutes[C]//17th AIAA Aerodynamic Decelerator Systems Technology Conference and Seminar, Monterey, California, AIAA 2003-2129, 2003: 1-16.

[6] LINGARD J, BARNARD S, KEARNEY P. Comparative Study of the Performance of Parachutes at Mach 0.5 to Mach 4.35 with Reference to Suitability for Use with the Hermes Crew Escape capsule[C]//10th Aerodynamic Decelerator Systems Technology Conference, Cocoa Beach, AIAA 1989-889, 1989: 1-11.

[7] WITKOWSKI A, KANDIS M, ADAMS D. Inflation Characteristics of the MSL Disk Gap Band Parachute[C]//20th AIAA Aerodynamic Decelerator Systems Technology Conference and Seminar
, Seattle, Washington, AIAA 2009-2915, 2009: 1-9.

[8] BENDURA R J, COLTRANE L C, HUCKINS, et al. Performance of a 19.7 Meter Diameter Disk Gap Band Parachute in A Simulated Martian Environment[C]//National Aeronautics and Space Administration, Washington, D. C., 1968: 1-36.

[9] GAO X, ZHANG Q, TANG Q. Numerical Modelling of Mars Supersonic Disk-gap-band Parachute Inflation[J]. Advances in Space Research, 2016, 57(11): 2259-2272.

[10] TAYLOR A. An Investigation of the Apparent Mass of Parachutes under Post Inflation Dynamic Loading through the Use of Fluid Structure Interaction Simulations[C]//17th AIAA Aerodynamic Decelerator Systems Technology Conference and Seminar, Monterey, California, AIAA 2003-2014, 2003: 1-9.

[11] 楊璐瑜, 張紅英, 陸偉偉, 等. 盤(pán)–縫–帶傘超聲速開(kāi)傘過(guò)程研究[J]. 航天返回與遙感, 2016, 37(3): 29-38. YANG Luyu, ZHANG Hongying, LU Weiwei, et al. Study on the Deployment of Disk-gap-band Parachute in Supersonic Flow[J]. Spacecraft Recovery and Remote Sensing, 2016, 37(3): 29-38. (in Chinese)

[12] 王中陽(yáng). 降落傘充氣過(guò)程動(dòng)力學(xué)數(shù)值模擬[D]. 南京: 南京航空航天大學(xué), 2013. WANG Zhongyang. The Dynamic Numerical Simulation of the Parachute Inflation[D]. Nanjing: Nanjing University of Aeronautics and Astronautics, 2013. (in Chinese)

[13] LINGARD J, DARLEY M, UNDERWOOD J C. Simulation of Mars Supersonic Parachute Performance and Dynamics[C]//19th AIAA Aerodynamic Decelerator Systems Technology Conference and Seminar, Williamsburg, VA, AIAA 2007-2507, 2007: 1-11.

[14] LINGARD J, DARLEY M. Simulation of Parachute Fluid Structure Interaction in Supersonic Flow[C]//18th AIAA Aerodynamic Decelerator Systems Technology Conference and Seminar, Munich Germany, AIAA 2005-1607, 2005: 1-9.

[15] SOULI M, OUAHSINE A, LEWIN L. ALE Formulation for Fluid–structure Interaction Problems[J]. Computer Methods in Applied Mechanics & Engineering, 2000, 190(5): 659-675.

[16] 楊璐瑜, 陸偉偉, 張紅英, 等. 速度對(duì)火星用盤(pán)–縫–帶傘超聲速開(kāi)傘性能影響[J]. 航空計(jì)算技術(shù), 2016, 46(5): 34-37. YANG Luyu, LU Weiwei, ZHANG Hongying, et al. Effect of Velocity on Performance of Mars Disk-gap-gand Parachutes in Supersonic Flow[J]. Aeronautical Computing Technique, 2016, 46(5): 34-37. (in Chinese)

[17] ZHANG Zengchan, JOHN Y S T. A Generalized Space-time CE/SE Method for the Euler Equations on Quadrilateral and Hexagonal Meshes[C]//15th AIAA Computational Fluid Dynamics Conference, Anaheim, CA, AIAA 2001-2592, 2001: 1-11.

[18] COOK G, ZHANG Z C, IM K. Applications of the CE/SE Method in LS-DYNA[C]//Applications of the CE/SE Method in LS-DYNA, San Diego, CA, AIAA 2013-3070, 2013: 1-10.

[19] TEZDUYAR T E, SATHE S, STEIN K. Solution Techniques for the Fully Discretized Equations in Computation of Fluid–structure Interactions with the Space–time Formulations[J]. Computer Methods in Applied Mechanics and Engineering, 2006, 195(41): 5743-5753.

[20] MITTAL R, IACCARINO G. Immersed Boundary Methods[J]. Annual Review of Fluid Mechanics, 2005, 37(37): 239-261.

[21] 劉凱欣, 王景燾, 王剛, 等. 時(shí)–空守恒元解元(CE/SE)方法綜述[J]. 力學(xué)進(jìn)展, 2011, 41(4): 447-461. LIU Kaixin, WANG Jingtao, WANG Gang, et al. A Review on the CE/SE Method[J]. Advances in Mechanics, 2011, 41(4): 447-461. (in Chinese)

[22] 郭叔偉, 王海濤, 董楊彪, 等. 降落傘“呼吸”現(xiàn)象研究[J]. 航天返回與遙感, 2010, 31(1): 18-23. GUO Shuwei, WANG Haitao, DONG Yangbiao, et al. Research on Parachute BREATH Behavior[J]. Spacecraft Recovery & Remote Sensing, 2010, 31(1): 18-23. (in Chinese)

[23] 韓晉陽(yáng), 徐宏, 高峰. 超聲速半流傘設(shè)計(jì)與分析[J]. 航天返回與遙感, 2013, 34(5): 20-28. HAN Jinyang, XU Hong, GAO Feng. Design and Analysis of Supersonic Half-flow Parachute[J]. Spacecraft Recovery & Remote Sensing, 2013, 34(5): 20-28. (in Chinese)

（編輯：陳艷霞）

Study on the Simulation of the Inflating Process of Disk-gap-band Parachute in Supersonic Flow

WANG Qi CAO Yihua

（Beihang University, Beijing 100083, China）

With the development of Mars exploration and landing, the research on the decelerator, the disk-gap-band parachute, is becoming more and more popular, and the analysis of the inflation process is one of the most complicated problems. In this paper, the numerical simulation methods used to solve the flow field of parachute are introduced, and their advantages and disadvantages are briefly described. This paper use conservation element and solution element (CE/SE) method to simulate the inflation process of the parachute from folding state in supersonic speed. Firstly, the probe is simulated separately, then the inflating process of the whole system is simulated under different Mach numbers. The simulation result are compared with the experimental data, and the comparison results are consistent, which verified the feasibility of the method. In this paper, the inflating process under two working conditions is compared, and the similarities and differences are analyzed. The resistance contribution of the probe changing along with the Mach number is also analyzed. Finally, the phenomenon of failure in the inflation process was analyzed by the contours of the parachute. The research result can provide reference for Mars exploration.

disk-gap-band parachute; fluid-structure interaction (FSI); CE/SE method; supersonic; probe landing

V445.2+3

A

1009-8518(2018)01-0035-10

10.3969/j.issn.1009-8518.2018.01.005

王祁，男，1992年生。2015年獲北京航空航天大學(xué)飛行器設(shè)計(jì)專業(yè)學(xué)士學(xué)位。現(xiàn)為北京航空航天大學(xué)航空工程專業(yè)碩士研究生。研究方向?yàn)榻德鋫慊厥障到y(tǒng)流固耦合數(shù)值模擬。E-mail: shuimohaichen@163.com。

曹義華，男，1962年生。1990年獲南京航空航天大學(xué)飛行器設(shè)計(jì)專業(yè)空氣動(dòng)力學(xué)博士學(xué)位，1992年南京理工大學(xué)氣體動(dòng)力學(xué)博士后出站。教授，研究方向?yàn)橛?jì)算流體力學(xué)研究和飛行器設(shè)計(jì)。E-mail: yihuacaobu@126.com。

2017-08-22