超大型組合翼傘開傘性能研究

吳瓊 張紅英 段非 胡雪鵬

超大型組合翼傘開傘性能研究

吳瓊1張紅英1段非1胡雪鵬2

（1 南京航空航天大學(xué)航空學(xué)院，南京 210016）（2 航空工業(yè)航宇救生裝備有限公司，襄陽 441003）

為深入了解超大型組合翼傘與連續(xù)翼傘開傘過程的差異，文章通過數(shù)值模擬計算對超大型組合翼傘和連續(xù)翼傘的開傘性能進(jìn)行了對比研究。首先，基于自由曲面變形（Free-Form-Deformation，F(xiàn)FD）理論建立了超大型翼傘的展向折疊模型，然后應(yīng)用任意拉格朗日-歐拉（Arbitrary Lagrange-Euler，ALE）方法對翼傘折疊模型進(jìn)行了充氣過程的流固耦合仿真計算。分析了超大型組合翼傘和連續(xù)翼傘充氣過程中的傘衣外形、傘衣應(yīng)力和氣室充氣規(guī)律，深入對比了超大型組合翼傘和連續(xù)翼傘的充滿時間和開傘動載。研究表明，連續(xù)翼傘和組合翼傘的開傘性能差異很小，兩者的充滿時間相差不大，組合翼傘展開效果更好，動載之差在5%左右。與連續(xù)翼傘相比，組合翼傘更易于制造和使用。研究將為超大型翼傘的設(shè)計優(yōu)化提供一定的參考。

超大型組合翼傘 充氣展開 自由曲面變形理論 任意拉格朗日-歐拉方法 開傘動載

0 引言

翼傘，相比于傳統(tǒng)降落傘，具有高滑翔比、可操縱性強(qiáng)、著陸精準(zhǔn)及可靠性強(qiáng)等優(yōu)勢，在航空、航天、民用等領(lǐng)域得到了廣泛應(yīng)用，近些年來更成為國內(nèi)外研究重點[1]。翼傘最具有代表性的應(yīng)用是在20世紀(jì)90年代，美國對X-38原型機(jī)開展了無人駕駛的投放試驗[2-4]。他們利用翼傘的減速與滑翔特性，使樣機(jī)最后安全定點著陸。此次試驗成功后，其他各國航天機(jī)構(gòu)紛紛開始加入到翼傘技術(shù)的研究與試驗中。

翼傘按上下翼面的傘衣是否為連續(xù)曲面可以分為連續(xù)翼傘與組合翼傘。連續(xù)翼傘即常規(guī)翼傘，上下翼面的傘衣均為一個連續(xù)曲面。組合翼傘上下翼面的傘衣均為非連續(xù)曲面，由多個模塊組合而成，每個模塊包含多個氣室，相鄰模塊的組合方式通過上下翼面設(shè)置連接點進(jìn)行連接，在連接點位置設(shè)置連接扣，再通過快卸式連接繩進(jìn)行連接，以便能夠快速固定和分離。

近些年，我國對連續(xù)翼傘進(jìn)行了諸多數(shù)值仿真研究。汪龍芳計算了翼傘從自然展開狀態(tài)到充氣鼓包狀態(tài)的充氣過程[5]；聶帥提出了一種基于R-N流體傳輸?shù)姆謪^(qū)耦合迭代方法，用于翼傘展開過程的流固耦合數(shù)值模擬[6]；張思宇等人開展了翼傘充氣過程中的流固耦合動力學(xué)特性研究[7]；陸偉偉等人對某大型翼傘進(jìn)行了三維定常數(shù)值模擬[8]；田似營等人采用流場、結(jié)構(gòu)松耦合方法數(shù)值模擬翼傘的充氣過程[9]；熊宇超等人通過仿真證明了穩(wěn)定傘和減速傘可以使火箭助推器穩(wěn)定減速，為翼傘開傘提供前提條件[10]；王立武通過調(diào)研國際上航天器翼傘精確回收領(lǐng)域發(fā)展情況并差距分析，初步規(guī)劃了技術(shù)路線，擬通過設(shè)計、仿真、試驗互相迭代實現(xiàn)關(guān)鍵技術(shù)的理論突破[11]。上述研究中，并沒有人使用流固耦合方法仿真計算研究組合翼傘的開傘性能。

由于大型連續(xù)翼傘具有面積大、質(zhì)量大的特點，在加工生產(chǎn)、包裝使用、回收維護(hù)等方面都存在一系列問題。國外已開展大型組合翼傘的探索研究，借鑒常規(guī)重裝空投系統(tǒng)中群傘設(shè)計理念，使用模塊組合的方式來進(jìn)行大型翼傘的設(shè)計。美國研制的組合翼傘Megafly成功開傘驗證了組合翼傘設(shè)計的可行性[12]。Jean-Christophe Berland等人將組合翼傘進(jìn)行了40次的空投實驗[13]。國內(nèi)楊鋒魁對組合翼傘也開展了風(fēng)洞試驗[14]。該思路可以有效解決大型翼傘在加工、使用、維護(hù)等方面存在的一系列問題。

從總體上看，國內(nèi)關(guān)于超大型組合翼傘的仿真研究尚屬空白，本文建立超大型組合翼傘的流固耦合仿真模型，通過仿真計算與分析，對比分析超大型組合翼傘和超大型連續(xù)翼傘的開傘性能，以期為超大型翼傘的優(yōu)化設(shè)計與開傘性能研究提供一定的參考依據(jù)。

1 理論基礎(chǔ)

1.1 流固耦合控制方程

翼傘展開過程是一個柔性體大變形的流固耦合過程，使用拉格朗日法或歐拉法進(jìn)行結(jié)構(gòu)運動的描述會導(dǎo)致網(wǎng)格畸變嚴(yán)重，計算終止。任意拉格朗日-歐拉（Arbitrary Lagrang-Euler，ALE）方法結(jié)合了拉格朗日法和歐拉法的優(yōu)點，非常適合對大變形結(jié)構(gòu)的運動進(jìn)行描述。本文在ALE法基礎(chǔ)上采用結(jié)構(gòu)化任意拉格朗日歐拉（Structured Arbitrary Lagrange-Euler，S-ALE）方法進(jìn)行流固耦合仿真，采用結(jié)構(gòu)化的網(wǎng)格結(jié)構(gòu)，使得單元和節(jié)點的連通性更加簡單，這將極大簡化ALE耦合方法中的搜索算法。

S-ALE方法的流體控制方程由質(zhì)量守恒方程、動量守恒方程和能量守恒方程組成，如下所示：

網(wǎng)格控制方程為

翼傘的結(jié)構(gòu)控制方程為

1.2 翼傘展向折疊方法

翼傘為非展平曲面所構(gòu)成的多氣室-雙翼面結(jié)構(gòu)，并且通過開孔的肋片將內(nèi)部分割成多個氣室，正是因為翼傘結(jié)構(gòu)的復(fù)雜性，所以無法用常規(guī)的折疊方法來建立展向折疊模型。針對這一難點，借鑒計算機(jī)進(jìn)行三維視覺特效的曲面變形的技術(shù)，再結(jié)合翼傘實際折疊包裝的工序，本文提出了基于自由曲面變形（Free-Form-Deformation，F(xiàn)FD）理論的柔性翼傘折疊建模技術(shù)。該方法主要是通過控制幾何模型移動方向與位移，對曲面進(jìn)行相應(yīng)的變形處理，然后再通過齊次坐標(biāo)法的張量矩陣進(jìn)行展向的移動與復(fù)制等相關(guān)變換，從而順利解決了多氣室沖壓翼傘柔性傘衣織物的展向折疊建模問題，對進(jìn)一步基于數(shù)值模擬耦合方法開展翼傘充氣機(jī)制的研究具有重要意義。

圖1 FFD六面體控制晶格

翼傘任意一點對應(yīng)的全局坐標(biāo)為

Bernstein多項式

2 計算模型

2.1 翼傘幾何模型

參考Megafly翼傘的尺寸[12-13]，本文研究的836 m2的翼傘幾何模型如圖2所示，具體幾何參數(shù)如表1所示。翼傘肋片材料均采用的是具有大透氣性特性的紡織物，將中間氣室沿展向的中心面定義為中軸面，于是整具翼傘關(guān)于中軸面對稱，且傘繩匯交點在中軸面上。翼傘各氣室的肋片和左右兩側(cè)的傘衣上的傘繩展向均勻分布64排、弦向均勻分布9排，合計576根傘繩；操縱繩左右各有1根，每根操縱繩與9根傘繩相連，這9根傘繩分別連接在除了中間氣室以外的側(cè)邊氣室的所有肋片處。以此幾何模型為例，分別對超大型連續(xù)翼傘和超大型組合翼傘進(jìn)行建模仿真。

圖2 翼傘結(jié)構(gòu)示意

表1 翼傘的幾何參數(shù)

2.2 展向折疊模型

進(jìn)行超大型連續(xù)翼傘與超大型組合翼傘開傘過程研究的第一步就是需要進(jìn)行折疊建模。從幾何模型來說，翼傘的上下兩翼面是柔性非展平曲面，并且通過開孔的肋片將內(nèi)部分割成多個氣室，正是因為翼傘結(jié)構(gòu)的復(fù)雜性，所以無法用常規(guī)的折疊方法來建立折疊模型。針對這一難點，本文提出了基于FFD法的柔性翼傘折疊建模技術(shù)。該方法主要是通過控制幾何模型上移動的方向與位移，將曲面進(jìn)行相應(yīng)的變形處理，然后再通過齊次坐標(biāo)法的張量矩陣進(jìn)行展向的移動與復(fù)制等相關(guān)變換，從而順利解決了多氣室沖壓翼傘柔性傘衣織物的展向折疊建模問題，對進(jìn)一步基于數(shù)值模擬耦合方法開展翼傘充氣機(jī)制的研究具有重要意義。

為保證折疊建模的高效性，傘衣折疊過程僅對參與折疊過程的上下翼面和肋片進(jìn)行操作，再與傘載系統(tǒng)的其他結(jié)構(gòu)合并成完整模型。

翼傘每個氣室相同，且每個氣室均沿中軸面對稱，為了保持各氣室網(wǎng)格的一致性，首先建立半個氣室模型，再采用對稱鏡像復(fù)制的方法補齊氣室。

基于FFD的控制函數(shù)，分別進(jìn)行上翼面和下翼面的折疊。折疊好上下翼面之后，對肋片進(jìn)行復(fù)制后平移到相應(yīng)的位置上共節(jié)點，以保證翼面與肋片之間，翼面與翼面之間是連續(xù)的，這樣就簡單形成了一個氣室，再對單一氣室進(jìn)行復(fù)制鏡像，之后對形成的兩個氣室進(jìn)行復(fù)制后移動，經(jīng)過反復(fù)的這種操作后，最終便形成如圖3所示的單個氣室和多個氣室，再進(jìn)行重合點合并連接。

圖3 單個氣室與多個氣室折疊

最后連接上下翼面，并對尾緣進(jìn)行內(nèi)嵌處理，形成翼傘完整傘衣結(jié)構(gòu)的展向折疊模型如圖4所示。連續(xù)翼傘上下翼面的傘衣均為連續(xù)曲面，共有63個氣室，由62個承載肋片連接上下翼面，如圖4（a）和圖4（c）所示；而組合翼傘上下翼面的傘衣均為非連續(xù)曲面，分為5個模塊，中間是一個13氣室翼傘，緊接著是兩個14氣室翼傘，最外端是兩個11氣室翼傘，共有63個氣室，由58個承載肋片連接上下翼面，在連接處上翼面共有11個連接點，其中下翼面有9個連接點，通過連接點將五個模塊連接縫合成一個整體，如圖4（b）和圖4（d）所示。

圖4 完整傘衣結(jié)構(gòu)的展向折疊模型

2.3 仿真計算模型

翼傘傘衣折疊模型建立完以后，再添加傘載系統(tǒng)的其他結(jié)構(gòu)并劃分網(wǎng)格，合并成完整的仿真計算模型。仿真計算模型如圖5所示。傘衣主要分為上翼面、下翼面、肋片（包括承載肋片與非承載肋片）這三個部分，均采用的是四邊形的殼單元，傘繩采用的是梁/繩索單元，采用共節(jié)點的方式加以連接。整體黃色網(wǎng)格為流場網(wǎng)格，流場域的尺寸為72.8 m×40.4 m×52.4 m，在傘衣面區(qū)域?qū)α鲌鼍W(wǎng)格進(jìn)行局部加密。采用有限質(zhì)量方法充氣，水平方向速度為17 m/s，豎直方向速度為20 m/s，開傘速度是最大開傘動載的主要影響因素，載荷質(zhì)量為13 700 kg，載荷的質(zhì)量會影響翼傘的下降速度。

圖5 計算模型網(wǎng)格

3 開傘性能分析

基于軟件的S-ALE流固耦合方法，進(jìn)行同樣的材料、算法、邊界條件等的關(guān)鍵字設(shè)置，數(shù)值模擬了超大型連續(xù)翼傘和超大型組合翼傘的開傘過程，分別獲得了傘衣充氣外形、傘衣應(yīng)力、氣室充氣規(guī)律以及開傘性能。開傘性能主要通過分析開傘時間與開傘動載來探究兩種傘型在開傘過程的異同。

3.1 充氣外形分析

超大型連續(xù)翼傘與組合翼傘充滿狀態(tài)下的外形如圖6所示。從圖6（a）和圖6（b），可以看出與連續(xù)翼傘相比，組合翼傘的前緣張開得更充分一些。由于尾緣處于封閉空間，膨脹速率更快，外形也更為飽滿；前緣受傘繩的約束作用，導(dǎo)致前緣膨脹程度較弱，且向中間氣室方向靠攏。此時上下翼面呈現(xiàn)出凹凸不平的鼓包狀態(tài)，主要是由于氣室內(nèi)充滿氣流，上下翼面受到內(nèi)腔氣流壓力向氣室外側(cè)膨脹，同時又因為上下翼面受到承載肋片的約束，從而在承載肋片與傘衣連接處呈現(xiàn)出向內(nèi)腔凹陷的外形，故在上下翼面表現(xiàn)為凹凸不平的波浪狀；分別對比圖6（a）與圖6（c）、圖6（b）與6（d）發(fā)現(xiàn)下翼面向氣室外側(cè)的凸出程度相對于上翼面較弱，這是由于下翼面連接有傘繩，受傘繩下拉作用，使得承載肋片向下移動，所以下翼面整體的凸出程度相對較弱。

3.2 傘衣應(yīng)力分析

經(jīng)仿真計算，連續(xù)翼傘和組合翼傘充氣過程的的應(yīng)力云圖如圖7所示。在氣室充氣過程中，翼傘前緣所受的約束較少，前緣切口處存在一定角度，且一直遭受氣流吹襲，擺動頻繁，0.3 s時連續(xù)翼傘和組合翼傘前緣均產(chǎn)生明顯的應(yīng)力集中現(xiàn)象，同時組合翼傘在各縫合接口處也出現(xiàn)應(yīng)力集中現(xiàn)象，需要在此處適當(dāng)?shù)卦黾訌?qiáng)度。各氣室受內(nèi)外氣流的壓力與肋片約束力的共同作用，使得翼面呈現(xiàn)出凹凸不一的外形。由于展向每排的傘繩長度一致，導(dǎo)致傘衣充滿狀態(tài)下下翼面受傘繩拉力成一定弧度的外形。初期傘衣受高壓氣流吹襲，承受較大氣動力，應(yīng)力峰值大約出現(xiàn)在0.30 s～0.76 s之間。同時受到傘繩的下拉作用，傘衣應(yīng)力主要集中在承載肋片周邊，隨著開傘過程阻力面積的增大，0.88 s高壓區(qū)向四周擴(kuò)散，傘衣載荷減小。

圖6 超大型翼傘充滿外形

圖7 充氣過程傘衣應(yīng)力云圖

3.3 氣室充氣規(guī)律

翼傘前緣、中部、尾緣三部分充氣過程的展長變化見圖8。連續(xù)翼傘和組合翼傘的前緣分別在0.75 s、0.76 s達(dá)到展長最大值24.8 m、28.0 m；中部在0.63 s、0.65 s達(dá)到展長最大值35.8 m、38.7 m；尾緣在0.58 s、0.57 s達(dá)到展長最大值43.0 m、43.3 m。由于傘繩和傘衣自身具有彈性，各部位在充氣充滿后產(chǎn)生了一定的回彈并伴有氣室“鼓包”現(xiàn)象?？梢钥闯鲇捎跉饬鞯耐苿?，氣體首先進(jìn)入氣室尾部，由于尾緣處于封閉空間，膨脹速率最快，前緣膨脹速率最慢。充滿后尾緣的充氣效果最佳，氣室最飽滿。通過對比連續(xù)翼傘和組合翼傘的充滿時間，發(fā)現(xiàn)兩者的充滿時間相差不大，充滿時組合翼傘的展長比連續(xù)翼傘的展長更長，組合翼傘展開效果更好。

圖8 充氣過程的展長變化

圖9 開傘動載隨時間的變化曲線

3.4 開傘動載對比

圖9為連續(xù)翼傘和組合翼傘開傘動載的對比情況，從圖中可以看出：兩種傘形得到的降落傘開傘動載曲線變化較為接近，連續(xù)翼傘開傘動載的峰值時間為0.266 s，載荷峰值為570 kN；組合翼傘開傘動載的峰值時間為0.222 s，載荷峰值為542 kN。由于組合翼傘是多塊傘衣拼接而成的，傘衣之間存在縫隙，部分氣流從縫隙流出未經(jīng)過傘衣，故組合翼傘相對連續(xù)翼傘動載較小一點，動載之差在5%左右。

4 結(jié)論

本文對超大型連續(xù)翼傘和超大型組合翼傘的折疊充氣展開過程進(jìn)行了數(shù)值模擬，研究了超大型組合翼傘的開傘性能，主要得出以下結(jié)論：

1）充氣外形。連續(xù)翼傘和組合翼傘在充滿時刻，上下翼面均呈現(xiàn)凹凸不平的鼓包狀態(tài)，但下翼面整體的凸出程度較弱，這是由于下翼面連接有傘繩，會受到傘繩的下拉作用，承載肋片向下移動。與連續(xù)翼傘相比，組合翼傘的前緣張開的更充分一些。

2）傘衣應(yīng)力。連續(xù)翼傘和組合翼傘前緣均產(chǎn)生明顯的應(yīng)力集中現(xiàn)象，同時組合翼傘在各縫合接口處也出現(xiàn)應(yīng)力集中現(xiàn)象，翼面受傘繩拉力成一定弧度的外形。初期傘衣應(yīng)力主要集中在承載肋片周邊，隨著開傘過程阻力面積的增大，高壓區(qū)向四周擴(kuò)散，傘衣載荷減小。

3）氣室充氣規(guī)律。由于氣流的推動，氣體首先進(jìn)入氣室尾部，由于尾緣處于封閉空間，膨脹速率最快，前緣膨脹速率最慢。充滿后尾緣的充氣效果最佳，氣室最飽滿。通過對比連續(xù)翼傘和組合翼傘的充滿時間，發(fā)現(xiàn)兩者的充滿時間相差不大，充滿時組合翼傘的展長比連續(xù)翼傘的展長更長，因此組合翼傘展開效果更好。

4）開傘性能。連續(xù)翼傘和組合翼傘的開傘性能差別很小，充滿時間相差不大，組合翼傘開傘動載相對較小一點，動載之差在5%左右。這是由于組合翼傘是多塊傘衣拼接而成的，傘衣之間存在縫隙，部分氣流從縫隙流出未經(jīng)過傘衣。

[1] 王利榮. 降落傘理論與應(yīng)用[M]. 北京: 宇航出版社, 1997: 528-600. WANG Lirong. Theory and Application of Parachutes[M]. Beijing: China Astronautic Publishing House, 1997: 528-600. (in Chinese)

[2] JOHN S. Development of the NASA X-38 Parafoil Landing System[C]//15th CEAS/AIAA Aerodynamic Decelerator Systems Technology Conference, June 8-11, 1999, Toulouse, France.

[3] SYGULLA D, KL?DTKE R, Püttmann N, et al. X-38, CRV and CRV Evolution Program Overview and European Role[J]. Acta Astronautica, 2001, 48(5): 859-868.

[4] BENNETT T, FOX R. Design, Development & Flight Testing of the NASA X- 38 7500 ft2 Parafoil Recovery System[C] //17th AIAA Aerodynamic Decelerator Systems Technology Conference and Seminar, May 19-22, 2003, Monterey, CA, USA.

[5] 汪龍芳, 賀衛(wèi)亮. 基于索膜有限元模型的翼傘氣動變形仿真[J]. 北京航空航天大學(xué)學(xué)報, 2017, 43(1): 47-52. WANG Longfang, HE Weiliang. Parafoil Aerodynamic Deformation Simulation Based on Cable-membrane Finite Element Model[J]. Journal of Beijing University of Aeronautics and Astronautics, 2017, 43(1): 47-52. (in Chinese)

[6] NIE S, CAO Y, WU Z. Numerical Simulation of Parafoil Inflation via a Robin-neumann Transmission-based Approach[J]. Proceedings of the Institution of Mechanical Engineers, 2018, 232(4): 797-810.

[7] 張思宇, 余莉, 劉鑫. 翼傘充氣過程的流固耦合方法數(shù)值仿真[J]. 北京航空航天大學(xué)學(xué)報, 2020, 46(6): 1108-1115. ZHANG Siyu, YU Li, LIU Xin. Numerical Simulation of Parafoil Inflation Process Based on Fluid-structure Interaction Method[J]. Journal of Beijing University of Aeronautics and Astronautics, 2020, 46(6): 1108-1115. (in Chinese)

[8] 陸偉偉, 張紅英, 連亮. 大型翼傘的三維氣動性能分析[J]. 航天返回與遙感, 2015, 36(3): 1-10. LU Weiwei, ZHANG Hongying, LIAN Liang. A Three-dimensional Analysis on Aerodynamic Performance of a Large Parafoil[J]. Spacecraft Recovery & Remote Sensing, 2015, 36(3): 1-10. (in Chinese)

[9] 田似營, 常浩, 凌蕾, 等. 沖壓翼傘充氣過程的數(shù)值模擬[J]. 航空科學(xué)技術(shù), 2019, 30(8): 28-34. TIAN Siying, CHANG Hao, LING Lei, et al. Numerical Simulation of Ram-air Parachute in Inflation Stage[J]. Aeronautical Science & Technology, 2019, 30(8): 28-34. (in Chinese)

[10] 熊宇超, 張紅英, 陳建平, 等. 火箭助推器回收穩(wěn)定減速階段動力學(xué)仿真分析[J]. 飛行力學(xué), 2021, 39(5): 63-70. XIONG Yuchao, ZHANG Hongying, CHEN Jianping, et al. Dynamic Simulation Analysis on Stabilization and Deceleration Stage of Rocket Booster Recovery[J]. Flight Dynamics, 2021, 39(5): 63-70. (in Chinese)

[11] 王立武, 許望晶, 劉濤, 等. 航天器翼傘精確回收技術(shù)發(fā)展及展望[J]. 航天返回與遙感, 2020, 41(4): 21-30.WANG Liwu, XU Wangjing, LIU Tao, et al. Development and Prospect of Parafoil Precise Recovery Technology for Spacecraft[J]. Spacecraft Recovery & Remote Sensing, 2020, 41(4): 21-30. (in Chinese)

[12] DUNKER S, BERLAND J C. Modularity Concepts for a 30,000 lb Capacity Ram-air Parachute[C]// AIAA Aerodynamic Decelerator Systems Technology Conference and Seminar, May 21-24, 2007, Williamsburg, VA, USA.

[13] BERLAND J C, BARBER J, GARGANO B. Autonomous Precision Delivery of 42000 Pounds (19 000 kg) under One Parachute[C]//20th AIAA Aerodynamic Decelerator Systems Technology Conference and Seminar, May 4-7, 2009, Seattle, Washington, USA.

[14] 楊鋒魁, 劉琦. 基于風(fēng)洞試驗的組合翼傘氣動特性研究[J]. 測控技術(shù), 2022, 41(2): 70-74. YANG Fengkui, LIU Qi. Aerodynamic Characteristics of Combined Parafoil Based on Wind Tunnel Experiment[J]. Measurement & Control Technology, 2022, 41(2): 70-74. (in Chinese)

[15] 張湘玉. 基于細(xì)分的曲線曲面變形技術(shù)研究[D]. 南京: 南京航空航天大學(xué), 2010. ZHANG Xiangyu. Research on Curve and Surface Deformation Technology Based on Subdivision[D]. Nanjing: Nanjing University of Aeronautics and Astronautics, 2010. (in Chinese)

[16] 呂狀. 三維網(wǎng)格模型編輯變形方法的研究與實現(xiàn)[D]. 西安: 西北大學(xué), 2011. LYU Zhuang. Research and Implementation of 3D Mesh Model Editing and Deformation Method[D]. Xi’an: Northwest University, 2011. (in Chinese)

[17] SEDERBERG T W, PARRY S R. Free-form Deformation of Solid Geometric Models[J]. ACM Siggraph Computer Graphics, 1986, 20(4): 151-160.

[18] 劉鑫. 基于流固耦合方法的翼傘開傘過程仿真研究[D]. 南京: 南京航空航天大學(xué), 2019. LIU Xin. Simulation of Parachute Opening Process Based on Fluid-solid Coupling Method[D]. Nanjing: Nanjing University of Aeronautics and Astronautics, 2019. (in Chinese)

[19] 張雄, 陸明萬, 王建軍. 任意拉格朗日歐拉描述法研究進(jìn)展[J]. 計算力學(xué)學(xué)報, 1997, 14(1): 91-102. ZHANG Xiong, LU Mingwan, WANG Jianjun. Research Progress in Arbitrary Lagrangian-eulerian Method[J]. Chinese Journal of Computational Mechanics, 1997, 14(1): 91-102. (in Chinese)

[20] 簡相輝, 金哲巖. 降落傘工作過程數(shù)值模擬研究綜述[J]. 航空科學(xué)技術(shù), 2016, 27(10): 1-7. JIAN Xianghui, JIN Zheyan. Review on the Development of Numerical Simulations on Parachutes[J]. Aeronautical Science & Technology, 2016, 27(10): 1-7. (in Chinese)

[21] 張宇晨, 胡雪鵬, 張紅英, 等. 沖壓式翼傘滑布控制開傘技術(shù)研究[J]. 航天返回與遙感, 2023, 44(2): 14-23. ZHANG Yuchen, HU Xuepeng, ZHANG Hongying, et al. Research on Slider-reefed Control Technology of Ram-air Parafoil[J]. Spacecraft Recovery & Remote Sensing, 2023, 44(2): 14-23. (in Chinese)

[22] FOGELL N, SHERWIN S J, COTTER C J, et al. Fluid-structure Interaction Simulation of the Inflated Shape of Ram-air Parachutes[C]//AIAA Aerodynamic Decelerator Systems Technology Conference, March 25-28, 2013, Daytona Beach, Florida, USA.

[23] STEIN K, BENNEY R, KALRO V. Parachute Fluid-structure Interactions: 3-D Computation[J]. Computer Methods in Applied Mechanics and Engineering, 2000, 190(3/4): 373-386.

[24] TEIXEIRA P R, AWRUCH A M. Numerical Simulation of Fluid-structure Interaction Using the Finite Element Method[J]. Computers & Fluids, 2005, 34(2): 249-273.

Research on Opening Performance of Super Large Combined Parafoil

WU Qiong1ZHANG Hongying1DUAN Fei1HU Xuepeng2

（1 Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China）（2 AVIC Aerospace Life-Support Industries Ltd., Xiangyang 441003, China）

In order to understand the difference of opening process between super large combined parafoil and continuous parafoil, the opening performance of super large combined parafoil and continuous parafoil in the process of inflation is studied and analyzed through numerical simulation. Firstly, the span-folding model of super large parafoil is established based on Free-Form-Deformation (FFD) theory, and then the inflatable expansion of the airfoil model is simulated with the Arbitrary Lagrangian Euler (ALE) method. The full shape, load distribution and flow field variation of super large combined and continuous parachutes in the process of inflation are analyzed. The filling time and opening dynamic load of super large combined and continuous parafoil are compared and analyzed deeply. The difference in opening performance between continuous parafoil and combined parafoil is very small, there is little difference in the filling time between the two types of parafoils, expansion effect is equivalent, and the difference of opening load is about 3%. Compared with the continuous parafoil, the combined parafoil is easier to manufacture and use. The above research will provide a certain reference for the design and optimization of super large parafoil.

super large combined parafoil; inflation; Free-Form-Deform (FFD) theory; Arbitrary Lagrange-Euler (ALE); opening dynamic load

V244.21+6

A

1009-8518(2023)05-0001-10

10.3969/j.issn.1009-8518.2023.05.001

吳瓊，女，1998年生，2021年獲河南理工大學(xué)工程力學(xué)專業(yè)學(xué)士學(xué)位，現(xiàn)在南京航空航天大學(xué)工程力學(xué)專業(yè)攻讀碩士學(xué)位。主要研究方向為流固耦合力學(xué)。E-mail：1335294618@qq.com。

2022-08-25

航空科學(xué)基金（20200029052004）

吳瓊, 張紅英, 段非, 等. 超大型組合翼傘開傘性能研究[J]. 航天返回與遙感, 2023, 44(5): 1-10.

WU Qiong, ZHANG Hongying, DUAN Fei, et al. Research on Opening Performance of Super Large Combined Parafoil[J]. Spacecraft Recovery & Remote Sensing, 2023, 44(5): 1-10. (in Chinese)

（編輯：陳艷霞）