方世興 黃偉 榮偉

?

盤(pán)縫帶傘細(xì)化結(jié)構(gòu)的仿真影響研究

方世興 黃偉 榮偉

（北京空間機(jī)電研究所，北京100094）

在火星探測(cè)任務(wù)中，盤(pán)縫帶傘對(duì)航天器在低密度大氣條件下的減速起到了重要作用。文章采用ProE軟件進(jìn)行幾何建模；利用Hyperwork軟件進(jìn)行傘衣有限元建模及前處理；利用LS-DYNA軟件的任意拉格朗日–歐拉方法著重研究了徑向加強(qiáng)帶細(xì)化結(jié)構(gòu)對(duì)于盤(pán)縫帶傘無(wú)限質(zhì)量下，流固耦合仿真結(jié)果的影響。通過(guò)仿真研究發(fā)現(xiàn)，在考慮徑向加強(qiáng)帶細(xì)化結(jié)構(gòu)后，縮比的盤(pán)縫帶傘的阻力系數(shù)值與實(shí)際試驗(yàn)結(jié)果相近，傘衣外形、傘衣擺角與實(shí)際試驗(yàn)結(jié)果更為接近。此外，考慮徑向加強(qiáng)帶后，傘衣應(yīng)力的仿真結(jié)果也與真實(shí)物理現(xiàn)象更為貼近。文章的建模方法及仿真結(jié)果對(duì)盤(pán)縫帶傘的設(shè)計(jì)與仿真驗(yàn)證具有一定的參考價(jià)值。

盤(pán)縫帶傘 細(xì)化結(jié)構(gòu) 流固耦合 有限元方法 回收著陸

0 引言

降落傘在飛船返回著陸、火箭助推器回收、返回式衛(wèi)星回收以及火星探測(cè)器進(jìn)入著陸等方面，都發(fā)揮著重要的作用[1]。深空探測(cè)任務(wù)已經(jīng)作為新世紀(jì)的探究熱點(diǎn)，而火星探測(cè)則是目前世界各國(guó)深空探測(cè)領(lǐng)域的研究重點(diǎn)?；鹦潜砻婢哂邢”〉拇髿?，探測(cè)器登陸火星時(shí)常采用氣動(dòng)外形與降落傘相結(jié)合的減速方案[2]。目前，成功的火星探測(cè)任務(wù)，大部分采用的降落傘是盤(pán)縫帶（Disk-Gap-Band，DGB）傘。特殊的火星環(huán)境與超聲速的開(kāi)傘條件，使得DGB傘的工作狀態(tài)更為復(fù)雜，難以在地球上直接進(jìn)行模擬試驗(yàn)[2]。目前國(guó)內(nèi)外主要通過(guò)風(fēng)洞試驗(yàn)、空投試驗(yàn)、高空開(kāi)傘試驗(yàn)來(lái)驗(yàn)證DGB傘的氣動(dòng)性能[3]。但這些試驗(yàn)往往費(fèi)用較高，又無(wú)法涵蓋所有的試驗(yàn)工況。因此，基于計(jì)算機(jī)技術(shù)的數(shù)值仿真方法就成為關(guān)注的熱點(diǎn)。

目前有許多適用于降落傘氣動(dòng)性能仿真的數(shù)值方法，LS-DYNA軟件的任意拉格朗日–歐拉方法（Arbitrary Lagrange Euler，ALE）能有效模擬出低速狀態(tài)下，降落傘在風(fēng)洞試驗(yàn)、空投試驗(yàn)中的流固耦合現(xiàn)象，其對(duì)應(yīng)的無(wú)限質(zhì)量[4-5]和有限質(zhì)量[6]的仿真計(jì)算結(jié)果也與試驗(yàn)結(jié)果相近，驗(yàn)證了ALE方法的有效性。在真實(shí)情況下，降落傘的加強(qiáng)帶、傘繩及吊帶等結(jié)構(gòu)，都是降落傘設(shè)計(jì)時(shí)需要重點(diǎn)關(guān)注的對(duì)象，也是保證降落傘強(qiáng)度的重要部分[7]；而以往國(guó)內(nèi)外的降落傘流固耦合仿真研究中，都對(duì)降落傘傘繩、加強(qiáng)帶及吊帶等結(jié)構(gòu)進(jìn)行了一定程度的簡(jiǎn)化[8-10]。在實(shí)際的風(fēng)洞試驗(yàn)中，由于縮比降落傘模型的尺寸較小，這些結(jié)構(gòu)的占比不可忽略，但對(duì)試驗(yàn)和仿真結(jié)果造成的影響目前尚不明確。國(guó)外的Tutt和Lingard團(tuán)隊(duì)利用LS-DYNA軟件，針對(duì)DGB傘模型進(jìn)行超聲速條件下的流固耦合仿真研究，并得到其不同馬赫數(shù)值條件下的氣動(dòng)特性，但都簡(jiǎn)化了DGB傘的傘繩及加強(qiáng)帶結(jié)構(gòu)[11-13]。文獻(xiàn)[14]曾利用浸入式邊界法，研究了超聲速條件下，柔性及剛性降落傘的氣動(dòng)特性。文獻(xiàn)[15]研究了傘繩在超聲速條件下剛體化降落傘氣動(dòng)仿真結(jié)果的影響，發(fā)現(xiàn)傘繩引起的激波與前體及降落傘本身的激波發(fā)生耦合。

針對(duì)DGB傘的傘繩及加強(qiáng)帶結(jié)構(gòu)對(duì)DGB傘的流固耦合仿真結(jié)果影響，目前國(guó)內(nèi)外的學(xué)者還未有過(guò)多的研究。因此，本文著重探究DGB傘細(xì)化結(jié)構(gòu)對(duì)于流固耦合仿真結(jié)果的影響。

本文采用ProE軟件進(jìn)行DGB傘縮比模型幾何建模，而后利用Hyperwork軟件對(duì)其進(jìn)行有限元建模及前處理，并給出了徑向加強(qiáng)帶細(xì)化結(jié)構(gòu)的建模方法；采用LS-DYNA軟件的ALE方法，對(duì)縮比的普通傘模型（后文稱(chēng)“普通傘”）和含細(xì)化結(jié)構(gòu)的細(xì)化傘模型（后文稱(chēng)“細(xì)化傘”）進(jìn)行無(wú)限質(zhì)量下的流固耦合仿真研究。文章對(duì)比了仿真計(jì)算與試驗(yàn)的阻力系數(shù)，并在此基礎(chǔ)上，對(duì)DGB傘徑向加強(qiáng)帶細(xì)化結(jié)構(gòu)帶來(lái)的影響進(jìn)行了研究，考察了充氣過(guò)程傘衣外形、傘衣擺角及傘衣應(yīng)力的變化情況，結(jié)果可為DGB傘的設(shè)計(jì)與仿真提供參考。

1 DGB傘的流固耦合計(jì)算機(jī)理

本文所采用的數(shù)值仿真軟件為L(zhǎng)S-DYNA，利用ALE方法[16-18]對(duì)DGB傘充氣過(guò)程進(jìn)行流固耦合仿真計(jì)算。

1.1 結(jié)構(gòu)動(dòng)力學(xué)模型

降落傘傘衣是柔性大變形體，具有典型的非線性動(dòng)力學(xué)特性，而且還是具有透氣量的多孔滲透性薄膜結(jié)構(gòu)，則結(jié)構(gòu)邊界的動(dòng)力學(xué)模型為[16]：

式中為位移矢量；s為材料密度；為作用在結(jié)構(gòu)上的體積力；s為結(jié)構(gòu)的Cauchy應(yīng)力張量；為積分時(shí)間；下標(biāo)s表示結(jié)構(gòu)。

1.2 流體動(dòng)力學(xué)模型

在ALE方法中，涉及流體的Navier-Stokes控制方程可用質(zhì)量守恒方程、動(dòng)量守恒方程和能量守恒方程來(lái)表示[16]。

質(zhì)量守恒方程為：

式中代表流體的密度；v代表相對(duì)空間域的流體速度；w表示相對(duì)參考域的流體速度；q表示流體的歐拉坐標(biāo)；=,,，表示節(jié)點(diǎn)的方向。

動(dòng)量守恒方程為：

式中為流體的應(yīng)力張量；b代表流體的體積力；,,表示的,,分量。

的具體形式可表述為：

式中為Kronecker函數(shù)；為流體的壓力；表示流體動(dòng)力粘性系數(shù)。

能量守恒方程為：

式中表示流體的內(nèi)能。

1.3 流固耦合算法

本文的數(shù)值算法是ALE方法，其主要特點(diǎn)是在建立幾何模型和進(jìn)行有限元網(wǎng)格劃分時(shí)，固體與流體的幾何外形以及對(duì)應(yīng)的網(wǎng)格可以重疊在一起，如圖1所示。計(jì)算中則通過(guò)罰函數(shù)約束方法將固體域與流體域耦合在一起，以實(shí)現(xiàn)力學(xué)參量的傳遞。

本文的研究對(duì)象是DGB傘，其常用的約束方法是罰函數(shù)法。罰函數(shù)耦合方法是通過(guò)數(shù)值計(jì)算追蹤結(jié)構(gòu)的從節(jié)點(diǎn)和流體的主物質(zhì)之間的相對(duì)距離，如圖2所示。檢查每一個(gè)從節(jié)點(diǎn)對(duì)主物質(zhì)表面的貫穿，如果從節(jié)點(diǎn)不出現(xiàn)貫穿，就不進(jìn)行任何操作；如果發(fā)生從節(jié)點(diǎn)對(duì)主物質(zhì)表面的貫穿，則界面接觸耦合力F就會(huì)分布到流體的節(jié)點(diǎn)上，界面接觸耦合力的大小與發(fā)生貫穿的數(shù)量呈正比，即：

式中k表示基于主從節(jié)點(diǎn)質(zhì)量模型特性的剛度系數(shù)，該值與實(shí)際的相對(duì)距離有關(guān)。

圖1 流體與固體的網(wǎng)格示意

Fig.1 Mesh of fluid and structure

圖2 罰函數(shù)耦合算法示意

2 仿真模型

本文仿真所采用的DGB傘參考“海盜號(hào)”（Viking）降落傘，其傘型參數(shù)比例見(jiàn)表1所示[19]。但考慮到計(jì)算成本的需要，本文根據(jù)某地面風(fēng)洞試驗(yàn)參數(shù)，對(duì)DGB傘的模型進(jìn)行了縮比，見(jiàn)后文所述。

表1 “海盜號(hào)”型DGB傘參數(shù)比例關(guān)系

Tab.1 The parameter relationship of the Viking type DGB parachute

2.1 DGB傘縮比模型參數(shù)

DGB傘縮比模型的結(jié)構(gòu)參數(shù)如表2所示，其中細(xì)化傘在幾何建模的時(shí)候僅考慮徑向加強(qiáng)帶，其余細(xì)化結(jié)構(gòu)均在有限元建模中予以展示。本文普通傘模型和細(xì)化傘模型的傘衣幾何建模結(jié)果如圖3所示。

表2 仿真DGB傘結(jié)構(gòu)參數(shù)

Tab.2 The structural parameter of the simulation parachute

（a）普通傘

（a）Normal parachute

（b）細(xì)化傘

通過(guò)圖3可以得知，細(xì)化傘模型與普通傘模型最大的區(qū)別，在于細(xì)化傘模型具有從傘頂孔到帶底部的徑向加強(qiáng)帶。徑向加強(qiáng)帶的寬度為3mm，且徑向加強(qiáng)帶與傘衣完全縫合。

為了提供后續(xù)有限元計(jì)算的輸入條件，本文利用Hyperwork軟件對(duì)兩種DGB傘模型進(jìn)行網(wǎng)格劃分及前處理。單幅傘衣的頂孔邊劃分為1個(gè)單元；盤(pán)底邊劃分為7個(gè)單元，傘衣盤(pán)部分沿軸向劃分為16個(gè)單元；傘衣帶部分底邊劃分為7個(gè)單元，沿軸向劃分為6個(gè)單元；為盡快獲得穩(wěn)態(tài)數(shù)據(jù)，仿真模型的傘繩劃分為1個(gè)單元。

為了探究徑向加強(qiáng)帶對(duì)于DGB傘流固耦合仿真結(jié)果的影響，針對(duì)LS-DYNA軟件流固耦合仿真的要求，本文對(duì)普通傘模型的加強(qiáng)帶采用1D繩索單元?jiǎng)澐?，而?duì)細(xì)化傘模型的徑向加強(qiáng)帶則采用2D殼單元?jiǎng)澐?。?jīng)過(guò)網(wǎng)格劃分后，DGB傘的網(wǎng)格信息如表3所示。

表3 DGB傘網(wǎng)格參數(shù)

Tab.3 Mesh parameters of the DGB parachute

2.2 流場(chǎng)模型參數(shù)

LS-DYNA軟件在進(jìn)行流固耦合仿真時(shí)，其計(jì)算總時(shí)長(zhǎng)與流場(chǎng)網(wǎng)格大小、網(wǎng)格總量有關(guān)。為了縮減計(jì)算總時(shí)長(zhǎng)，流場(chǎng)網(wǎng)格設(shè)置時(shí)僅在傘衣附近區(qū)域進(jìn)行加密，其余區(qū)域按線性變化進(jìn)行設(shè)置。為了盡可能模擬出風(fēng)洞條件，流場(chǎng)域選定為長(zhǎng)方體，其中分為流場(chǎng)入口、流場(chǎng)出口、風(fēng)洞壁面和中心試驗(yàn)段四部分，流場(chǎng)的幾何和網(wǎng)格參數(shù)見(jiàn)表4所示。為了進(jìn)一步模擬出無(wú)限質(zhì)量的試驗(yàn)條件，流場(chǎng)入口條件為恒定速度入口條件，流場(chǎng)出口設(shè)定為無(wú)反射出流條件，風(fēng)洞壁面設(shè)定為無(wú)滑移邊界條件，中心試驗(yàn)部分設(shè)定為滿(mǎn)足理想氣體狀態(tài)方程條件。

表4 流場(chǎng)域參數(shù)

Tab.4 The parameter of the fluid domain

2.3 仿真工況參數(shù)

在火星探測(cè)任務(wù)中，DGB傘是在超聲速條件下進(jìn)行充氣展開(kāi)的。而在地球環(huán)境中，則可以通過(guò)風(fēng)洞試驗(yàn)來(lái)測(cè)試DGB傘的阻力特性，其測(cè)試涵蓋低速、亞聲速、跨聲速和超聲速范圍[3]。考慮到LS-DYNA軟件在低速情況下的仿真結(jié)果是具有較高的可信度。因此，為了研究細(xì)化結(jié)構(gòu)對(duì)于降落傘流固耦合仿真結(jié)果的影響，本文選取地球大氣條件下的某低速風(fēng)洞試驗(yàn)參數(shù)作為仿真輸入條件，并將仿真結(jié)果與實(shí)際試驗(yàn)值進(jìn)行對(duì)比。仿真對(duì)應(yīng)的風(fēng)洞試驗(yàn)參數(shù)如表5所示，仿真過(guò)程中僅考慮傘衣織物透氣性，不考慮細(xì)化傘徑向加強(qiáng)帶的透氣性。

表5 DGB傘低速風(fēng)洞試驗(yàn)參數(shù)

Tab.5 Low-speed wind tunnel test parameters of DGB parachute

3 仿真結(jié)果與分析

3.1 阻力系數(shù)的變化分析

仿真所用的DGB傘模型的名義直徑為500mm，其名義面積為0.182m2，通過(guò)仿真得到的傘衣阻力可以求得阻力系數(shù)隨著時(shí)間的變化，如圖4所示。由圖4可知，兩種DGB傘從仿真初始狀態(tài)到充滿(mǎn)狀態(tài)的充氣時(shí)間都相對(duì)較短，細(xì)化傘在0.2s左右就達(dá)到了穩(wěn)定狀態(tài)，而普通傘的穩(wěn)定時(shí)間比細(xì)化傘更長(zhǎng)一些，接近0.35s。細(xì)化傘的呼吸頻率比普通傘的要大一些，細(xì)化傘的阻力系數(shù)值比普通傘的要偏大，這與實(shí)際的物理建模是符合的。由于徑向加強(qiáng)帶為2D單元，實(shí)際上是增加了傘衣的阻力面積，導(dǎo)致傘衣的阻力增大。經(jīng)過(guò)計(jì)算求得，細(xì)化傘阻力系數(shù)是0.687，與試驗(yàn)結(jié)果的誤差值為+8.36%，這說(shuō)明本文細(xì)化結(jié)構(gòu)的建模方法和計(jì)算結(jié)果是可信的。

圖4 DGB傘阻力系數(shù)與時(shí)間的關(guān)系曲線

3.2 傘衣外形的變化分析

為了更好地得到細(xì)化結(jié)構(gòu)對(duì)于DGB傘流固耦合仿真結(jié)果的影響，給出了3個(gè)任意時(shí)刻的傘衣變形圖，如圖5所示。圖中為仿真充氣時(shí)間。通過(guò)對(duì)比分析兩種不同DGB傘的變形情況，可以發(fā)現(xiàn)，在降落傘從初始外形到充氣張滿(mǎn)的過(guò)程中，由于徑向加強(qiáng)帶剛度較大，細(xì)化傘的變形量沒(méi)有普通傘的大，故外形變化并不劇烈。此外，隨著充氣時(shí)間的推移，細(xì)化傘的外形更為飽滿(mǎn)，但呼吸頻率略比普通傘要大一些。

圖5 DGB傘變形示意

圖6中顯示的是考慮徑向加強(qiáng)帶細(xì)化結(jié)構(gòu)后，DGB傘穩(wěn)定狀態(tài)下的變形對(duì)比圖。通過(guò)仿真結(jié)果圖片可以看出，在未考慮細(xì)化結(jié)構(gòu)時(shí)，DGB傘的盤(pán)與帶之間間隙較大，這也說(shuō)明了加強(qiáng)帶對(duì)降落傘充氣展開(kāi)的外形會(huì)有一定的影響。根據(jù)測(cè)量的結(jié)果顯示，考慮細(xì)化結(jié)構(gòu)的情況下，細(xì)化傘的盤(pán)直徑為317mm，普通傘的盤(pán)直徑為315mm；細(xì)化傘的實(shí)際盤(pán)直徑較大一些，導(dǎo)致傘衣的阻力面積增大，因此也對(duì)仿真結(jié)果也造成一定的影響。

（a）普通傘

（a）Normal parachute

（b）細(xì)化傘

本文所建立的DGB傘模型尺寸較小，而徑向加強(qiáng)帶細(xì)化結(jié)構(gòu)，除了在尺寸上會(huì)對(duì)建模結(jié)果產(chǎn)生一定的影響，也從材料剛度上對(duì)DGB傘產(chǎn)生較大的影響。從理論上來(lái)說(shuō)，當(dāng)考慮了降落傘加強(qiáng)帶后，一定程度上是減小了結(jié)構(gòu)的透氣量，而徑向加強(qiáng)帶剛度較大，變形時(shí)會(huì)進(jìn)一步約束細(xì)化傘的外形，導(dǎo)致進(jìn)一步增加了DGB傘的阻力面積。而如上文所述，展開(kāi)后考慮加強(qiáng)帶的降落傘阻力面積更大。但是通過(guò)對(duì)比實(shí)際試驗(yàn)外形結(jié)果和仿真結(jié)果的外形，如圖7所示，細(xì)化傘更接近于實(shí)際結(jié)果；這也說(shuō)明了在仿真時(shí)，需要考慮細(xì)化結(jié)構(gòu)帶來(lái)的影響。

（a）仿真結(jié)果

（a）Simulation result

（b）試驗(yàn)結(jié)果

3.3 傘衣擺角的變化分析

通過(guò)進(jìn)一步分析仿真計(jì)算得到的結(jié)果，可以獲得DGB傘兩種傘模型對(duì)應(yīng)的側(cè)向力系數(shù)C（=,），如圖8所示。通過(guò)圖8可以看出，普通傘的側(cè)向力系數(shù)明顯要比細(xì)化傘的大很多，且其值出現(xiàn)類(lèi)似周期性的變化，這也說(shuō)明普通傘的擺動(dòng)較為劇烈。通過(guò)圖8可以得知，普通傘的向側(cè)向力系數(shù)較大，對(duì)降落傘的穩(wěn)定性造成較大的影響。

本文通過(guò)分析計(jì)算傘頂孔中心位置到傘繩匯交點(diǎn)的向量與DGB傘初始狀態(tài)軸向向量的夾角，進(jìn)一步探究細(xì)化結(jié)構(gòu)對(duì)于仿真結(jié)果的影響。傘衣的擺角隨著時(shí)間變化的曲線如圖9所示。風(fēng)洞試驗(yàn)結(jié)果顯示，DGB傘在低速情況下的擺角為2°~5°。通過(guò)圖9可以看出，普通傘的擺角是屬于正弦周期的擺動(dòng)，且擺角在12°~25°之間，擺角較大，與實(shí)際的試驗(yàn)值相差較大。而細(xì)化傘的擺角則在穩(wěn)定后基本不變，其擺角在2°~3°之間，與實(shí)際的試驗(yàn)值相符。通過(guò)對(duì)比可以發(fā)現(xiàn)，考慮了徑向加強(qiáng)帶的仿真結(jié)果，在降落傘的傘衣擺角上與實(shí)際結(jié)果更為相近，細(xì)化結(jié)構(gòu)在流固耦合仿真時(shí)應(yīng)該予以考慮。

圖9 傘衣擺角與時(shí)間的關(guān)系

3.4 傘衣應(yīng)力的變化分析

圖10是普通傘和細(xì)化傘在傘衣應(yīng)力分布的變化情況。通過(guò)對(duì)比結(jié)果可以發(fā)現(xiàn)，細(xì)化傘的應(yīng)力集中分布區(qū)域與普通傘的應(yīng)力集中分布區(qū)域有較大的區(qū)別。

圖10 傘衣應(yīng)力分布

普通傘的應(yīng)力集中區(qū)域主要是在傘衣盤(pán)部分的中間部分和底邊部分；而細(xì)化傘的應(yīng)力集中區(qū)域則主要在傘頂孔附近的加強(qiáng)帶及傘衣徑向加強(qiáng)帶的鄰近部分。通過(guò)對(duì)比兩者可以看出，細(xì)化傘的應(yīng)力集中區(qū)域與實(shí)際實(shí)驗(yàn)中應(yīng)力集中區(qū)域更接近；而實(shí)際風(fēng)洞試驗(yàn)中降落傘在遭受破壞時(shí)，也主要是在加強(qiáng)帶附近出現(xiàn)撕裂或破壞的情況[7]。此外，細(xì)化傘的應(yīng)力計(jì)算值比普通傘更大，能夠?yàn)榻德鋫銖?qiáng)度設(shè)計(jì)提供更高的閾值，從而保證降落傘在實(shí)際試驗(yàn)中不易遭受破壞。因此，在降落傘的流固耦合仿真中，徑向加強(qiáng)帶應(yīng)著重予以考慮，不應(yīng)忽略。

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

本文采用軟件（ProE和Hyperwork）進(jìn)行傘衣結(jié)構(gòu)建模，利用ALE方法，著重研究了徑向加強(qiáng)帶細(xì)化結(jié)構(gòu)對(duì)于盤(pán)縫帶傘無(wú)限質(zhì)量下流固耦合仿真結(jié)果的影響，并得出以下結(jié)論：

1）徑向加強(qiáng)帶的簡(jiǎn)化結(jié)構(gòu)和細(xì)化結(jié)構(gòu)對(duì)降落傘的阻力系數(shù)影響不大，細(xì)化傘模型的仿真結(jié)果在誤差允許范圍內(nèi)，說(shuō)明本文的建模方法與仿真計(jì)算方法可信。

2）考慮徑向加強(qiáng)帶細(xì)化結(jié)構(gòu)后，傘衣的變形量較?。患?xì)化傘衣外形的仿真結(jié)果與實(shí)際風(fēng)洞試驗(yàn)結(jié)果更相近，仿真時(shí)應(yīng)考慮細(xì)化結(jié)構(gòu)。

3）徑向加強(qiáng)帶的簡(jiǎn)化結(jié)構(gòu)對(duì)降落傘的擺角產(chǎn)生較大的影響。在考慮細(xì)化結(jié)構(gòu)后，盤(pán)縫帶傘的側(cè)向力系數(shù)更小，傘衣擺角更小。因此，在研究降落傘擺動(dòng)性能時(shí)，應(yīng)考慮徑向加強(qiáng)帶。

4）考慮細(xì)化結(jié)構(gòu)后，盤(pán)縫帶傘的應(yīng)力集中主要在徑向加強(qiáng)帶的細(xì)化結(jié)構(gòu)上；對(duì)降落傘結(jié)構(gòu)應(yīng)力的模擬結(jié)果更加真實(shí)，能夠更有效地指導(dǎo)降落傘強(qiáng)度設(shè)計(jì)工作。

本文僅針對(duì)盤(pán)縫帶傘的風(fēng)洞試驗(yàn)縮比模型進(jìn)行了研究，全尺寸條件下細(xì)化結(jié)構(gòu)的影響是未來(lái)進(jìn)一步研究的重點(diǎn)。此外，本文研究?jī)?nèi)容僅考慮了徑向加強(qiáng)帶細(xì)化結(jié)構(gòu)對(duì)降落傘仿真結(jié)果的影響，未來(lái)也可進(jìn)一步研究傘繩等對(duì)仿真結(jié)果的影響。在未來(lái)的工作中，還可在不同速度條件下，開(kāi)展更多的降落傘細(xì)化結(jié)構(gòu)的流固耦合仿真影響研究，以明確定量的規(guī)律。

(References)

[1] 王海濤. 大型降落傘抽打現(xiàn)象及運(yùn)動(dòng)穩(wěn)定性研究[D]. 長(zhǎng)沙: 國(guó)防科學(xué)技術(shù)大學(xué), 2011. WANG Haitao. Research on Bull Whipping and Dynamic Stability of Large Parachute System[D]. Changsha: National University of Defense Technology, 2011. (in Chinese)

[2] 榮偉. 火星探測(cè)器減速著陸技術(shù)研究[D]. 北京: 中國(guó)空間技術(shù)研究院, 2008. RONG Wei. Mars Exploration Deceleration Landing Technology Research[D]. Beijing: China Academy of Science and Technology, 2008. (in Chinese)

[3] GILLIS G. The Viking Decelerator System an Overview[C]. 4th Aerodynamic Deceleration Systems Conference, Palm Springs, CA, USA, AIAA 1973-442, 1973.

[4] TUTT B, TAYLOR A. The Use of LS-DYNA to Simulate the Inflation of a Parachute Canopy[C]//AIAA Aerodynamic Decelerator Systems Technology Conference and Seminar, Munich, Germany, AIAA 2005-1608, 2005.

[5] 高興龍, 高慶玉, 張青斌, 等. 基于ALE方法的開(kāi)縫降落傘充氣過(guò)程研究[J]. 航天返回與遙感, 2013, 34(1): 13-20. GAO Xinglong, GAO Qingyu, ZHANG Qingbin, et al. A Study on the Slots-parachute Inflation Based on ALE Method[J]. Spacecraft Recovery & Remote Sensing, 2013, 34(1): 13-20. (in Chinese)

[6] 程涵, 余莉, 楊雪松, 等. 有限質(zhì)量情況下降落傘開(kāi)傘過(guò)程數(shù)值仿真研究[J]. 空氣動(dòng)力學(xué)學(xué)報(bào), 2014, 32(2): 259-263. CHEN Han, YU Li, YANG Xuesong, et al. Numerical Simulation of Parachute Opening Process in Finite Mass Situation[J]. Acta Aerodynamica Sinca, 2014, 32(2): 259-263. (in Chinese)

[7] 王利榮. 降落傘理論與應(yīng)用[M]. 北京: 宇航出版社, 1997. WANG Lirong. Parachute Theory and Application[M]. Beijing: Aerospace Press, 1997. (in Chinese)

[8] 薛曉鵬, 溫志涌. 超音速降落傘系統(tǒng)的氣動(dòng)干擾數(shù)值模擬研究[J]. 航天返回與遙感, 2016, 37(3): 9-18. XUE Xiaopeng, WEN Chihyung. Numerical Simulation of Aerodynamic Interaction of a Supersonic Parachute System[J]. Spacecraft Recovery & Remote Sensing, 2016, 37(3): 9-18. (in Chinese)

[9] 王中陽(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)

[10] 賈賀, 榮偉, 陳國(guó)良. 基于LS-DYNA軟件的降落傘充氣過(guò)程仿真研究[J]. 航天器環(huán)境工程, 2010, 27(3): 367-373. JIA He, RONG Wei, CHEN Guoliang. The Simulation Research of Parachute Inflation Process Based on LS-DYNA Software[J]. Spacecraft Environment Engineer, 2010, 27(3): 367-373. (in Chinese)

[11] LINGARD J, DARLEY M. Simulation of Parachute Fluid Structure Interaction in Supersonic Flow[R]. AIAA 2005-1607, 2005.

[12] LINGARD J, DARLEY M, UNDERWOOD J. Simulation of Mars Supersonic Parachute Performance and Dynamics[R]. AIAA 2007-2507, 2007.

[13] STEPHEN J L, UNDERWOOD J. Supersonic Tests of the Double Gap Disk-gap Band Parachute and Fluid Structure Interaction Simulation[R]. AIAA 2015-2109, 2015.

[14] XUE Xiaopeng, HIROTO K, YOSHIAKI N. Numerical Simulation on Supersonic Aerodynamic Interference for Rigid and Flexible Parachutes [C]. 42nd AIAA Fluid Dynamics Conference and Exhibit, Fluid Dynamics and Co-located Conferences, New Orleans, Louisiana, USA, AIAA 2012-3269, 2012.

[15] XUE Xiaopeng, HIROTO K, YOSHIAKI N, et al. Effects of Suspension Line on Flow Field Around a Supersonic Parachute [J]. Aerospace Science and Technology, 2015, 43: 63-77.

[16] 李裕春, 時(shí)黨勇, 趙遠(yuǎn). ANSYS 11.0/LS-DYNA基礎(chǔ)理論與工程實(shí)踐[M]. 北京: 中國(guó)水利水電出版社, 2008. LI Yuchun, SHI Dangyong, ZHAO Yuan. ANSYS 11.0/LS-DYNA Basic Theory and Engineering Application[M]. Beijing: China Waterpower press, 2008. (in Chinese)

[17] 白金澤. LS-DYNA3D理論基礎(chǔ)與實(shí)例分析[M]. 北京: 科學(xué)出版社, 2005. BAI Jinze. LS-DYNA 3D Basic Theory and Application analysis[M]. Beijing: Science Press, 2005. (in Chinese)

[18] 趙海鷗. LS-DYNA動(dòng)力分析指南[M]. 北京: 兵器工業(yè)出版社, 2003. ZHAO Haiou. LS-DYNA Dynamic Analysis Manual[M]. Beijing : Ordnance Industry Press, 2003.(in Chinese)

[19] 榮偉, 陳旭. 火星探測(cè)用降落傘研制試驗(yàn)簡(jiǎn)介[J]. 航天返回與遙感, 2007, 28(1): 12-17. RONG Wei, CHEN Xu. Introduction of MARS Exploration Parachute Research and Test[J]. Spacecraft Recovery & Romote Sensing, 2007, 28(1): 12-17. (in Chinese)

[20] STELZNERA, CRUZJ, BRUNOR, et al. Opportunities and Limitations in Low-earth Subsonic Testing and Qualification of Extraterrestrial Supersonic Parachute Designs [R].Lanley: NASA Jet Propulsion Laboratory, AIAA2003-2135.

（編輯：陳艷霞）

Study on the Detailed Structure of Disk-Gap-Band Parachute Simulation

FANG Shixing HUANG Wei RONG Wei

（Beijing Institute of Space Mechanics & Electricity, Beijing 100094, China）

In the Mars exploration mission, Disk-Gap-Band (DGB) parachute plays an important role in decreasing spacecraft's speed under low density atmosphere conditions. In this paper, ProE software is used to model the geometry of the canopy while Hyperwork software is used for its finite element modeling and pre-processing, and the effect of the radial reinforcement detailed structure of the DGB parachute on the wind tunnel simulation is analyzed by Arbitrary Lagrange Euler (ALE) method of the LS-DYNA software. The simulation results show that when considering the radial reinforcement detailed structure, drag coefficient of the scale canopy is closed to the experiment result, and the shape and swing angle also meet the experiment data. Furthermore, the stress of canopy are also analyzed according to the simulation results after considering the reinforcement detailed structure, which is more closed to actual physical phenomenons. The modeling method and the simulation results in this paper have certain reference value for the design and simulation verification of the DGB parachute.

Disk-Gap-Band parachute; detailed structure; fluid-structure interaction; finite element method; recovery and landing

V445

A

1009-8518(2017)02-0017-10

10.3969/j.issn.1009-8518.2017.02.003

2017-01-10

方世興，男，1992年生，2014年獲北京航空航天大學(xué)飛行器設(shè)計(jì)與工程（航天）專(zhuān)業(yè)學(xué)士學(xué)位，現(xiàn)在中國(guó)空間技術(shù)研究院飛行器設(shè)計(jì)專(zhuān)業(yè)攻讀碩士學(xué)位。研究方向?yàn)楹教炱鞣祷嘏c著陸技術(shù)。E-mail：fs_xing@163.com。