張延泰，孫建紅，2，侯 斌

（1.南京航空航天大學(xué) 飛行器環(huán)境控制與生命保障工業(yè)和信息化部重點(diǎn)實(shí)驗(yàn)室，南京 210016；2.南京航空航天大學(xué)民航學(xué)院，南京 210016）

馬航370航班失聯(lián)后，南京航空航天大學(xué)民機(jī)救生團(tuán)隊(duì)與中國商用飛機(jī)有限責(zé)任公司（中國商飛）針對民航客機(jī)聯(lián)合研發(fā)了分離式應(yīng)急記錄跟蹤系統(tǒng)，該系統(tǒng)又被稱為報(bào)信者（harbinger，HBG）系統(tǒng)[1-4].在HBG系統(tǒng)中，應(yīng)急飛行數(shù)據(jù)記錄器（emergency flight data recorder，EFDR）在飛機(jī)墜毀前彈射離機(jī).利用降落傘可以有效降低EFDR墜落速度，減小墜落沖擊.降落傘離機(jī)后迅速開傘，其傘衣充氣過程較為復(fù)雜[5]，而針對小型裝置設(shè)計(jì)的降落傘不多，對其展開過程的研究相對較少，其中小型傘是否存在傘頂塌陷及其塌陷機(jī)理均缺乏細(xì)致研究.如果小型傘傘衣在充氣階段發(fā)生塌陷或扭曲變形，盡管是局部也可能會(huì)導(dǎo)致開傘失敗.研究小型降落傘充氣過程，分析傘衣附近的流場特性，判斷其工作充氣過程中的安全性對小型裝置投放等工程應(yīng)用具有重要意義.

降落傘在充氣過程中，流體與傘衣的耦合作用具有非定常、柔性大變形等非線性特性，是流固耦合研究的難點(diǎn)之一.降落傘充氣過程的理論模型最早出現(xiàn)于1927年，是后續(xù)降落傘充氣過程計(jì)算方法的基礎(chǔ).在此基礎(chǔ)上，Mohaghegh等[6]研究了降落傘的充氣過程，發(fā)現(xiàn)充氣時(shí)間與傘型、傘口面積和來流速度等因素相關(guān).這些充氣模型多用于計(jì)算充氣時(shí)間和開傘動(dòng)載，沒有考慮傘衣局部變形等問題.Wolf等[7]研究了群傘的充氣過程，發(fā)現(xiàn)減速傘充滿后發(fā)生塌陷會(huì)嚴(yán)重影響傘衣的減速效果.隨后，NASA在獵戶座計(jì)劃中[8]發(fā)現(xiàn)傘衣在充氣過程中存在過度充氣和部分塌陷的問題.在實(shí)驗(yàn)的基礎(chǔ)上，Yavuz[9]針對物傘系統(tǒng)和傘衣尾流區(qū)速度分布建立了理論模型.彭勇等[10]基于流體連續(xù)方程和動(dòng)量方程，提出了降落傘充氣過程的尾流再附動(dòng)力學(xué)模型.在試驗(yàn)與仿真結(jié)合的研究中，王海濤[11]研究了大型降落傘充氣過程中的傘衣變形和抽打現(xiàn)象，分析了抽打現(xiàn)象的特點(diǎn)和形成原因.

此外，利用實(shí)驗(yàn)和數(shù)值模擬技術(shù)可以進(jìn)一步分析傘衣附近流場特性對傘衣變形的影響.Johari等[12]采用PIV和水洞試驗(yàn)技術(shù)研究了降落傘的尾流場，認(rèn)為渦脫落對傘衣背部有周期性的吸力作用.Fang等[13]采用數(shù)值模擬方法研究了十字傘的充氣過程，對比分析了單傘和雙傘系統(tǒng)，討論了傘衣附近流場特性對傘衣充氣過程的影響.孫建紅等[14]針對阻力傘實(shí)驗(yàn)中的非定常運(yùn)動(dòng)和尾流干擾，結(jié)合理論分析和數(shù)值模擬技術(shù)，分析了附加質(zhì)量模型，建立了傘衣阻力系數(shù)的修正方法，并與實(shí)驗(yàn)結(jié)果進(jìn)行了對比分析.

綜上所述，研究人員在較早的工程應(yīng)用中發(fā)現(xiàn)并討論了尾流再附與傘衣塌陷等問題，但研究主要關(guān)注傘衣投影直徑、充氣時(shí)間和開傘載荷的變化.針對小型降落傘的流場特性與傘衣變形相互作用的研究相對較少.本研究對小型降落傘的充氣和傘衣塌陷過程進(jìn)行數(shù)值模擬，分析了傘衣充氣過程中流場特性和傘衣變形的演化規(guī)律，研究了來流速度和前置體質(zhì)量對傘衣充氣過程的影響，以期給出傘衣安全充氣的工況范圍.

1 數(shù)值計(jì)算方法

降落傘的充氣過程是柔性體大變形流固耦合問題，可采用ALE方法進(jìn)行數(shù)值模擬，通過罰函數(shù)方法進(jìn)行流固耦合計(jì)算.根據(jù)美國國家運(yùn)輸安全委員會(huì)（national transportation safety board，NTSB）的事故報(bào)告，本研究工況中，設(shè)置HBG系統(tǒng)降落傘的開傘速度區(qū)間為40～160 m/s，可壓縮流體控制方程為

式（1）～式（3）中：t為時(shí)間；ρ為流體密度；ui為流體速度張量；p為流體壓力；fi為體積力張量；ε為流體比內(nèi)能；Φ為耗散函數(shù)；T為流體溫度.

由于ALE方法在歐拉和拉格朗日坐標(biāo)之外引入了參考坐標(biāo)，所以在ALE方法中，材料域的節(jié)點(diǎn)信息需要被映射到空間域內(nèi)，其流場變量的輸運(yùn)方程為

式（4）和式（5）中：φ為流場變量；φ0為當(dāng)前時(shí)間步流場變量的初始值；uig為材料速度張量，當(dāng)uig=0時(shí)對應(yīng)歐拉方法，當(dāng)uig=ui時(shí)對應(yīng)拉格朗日方法.該數(shù)值模擬方法在本課題組前期研究中用于分析降落傘和阻力傘的充氣過程[5，15].因此，有理由認(rèn)為該數(shù)值模擬方法適用于本研究.

2 計(jì)算模型

HBG系統(tǒng)的降落傘為平面圓形傘，傘衣名義直徑為1.9 m.根據(jù)平面圓形傘的折疊方法，傘衣首先按照傘衣幅做“平折疊”，隨后做“Z字形”折疊，最后放入傘包.在數(shù)值模擬中采用類似真實(shí)平面圓形傘折疊的方法進(jìn)行建模，圖1為傘衣折疊前后的面積變化和折疊后的傘衣網(wǎng)格.由于降落傘掛載的電子設(shè)備體積較小，所以計(jì)算模型中將前置體簡化為1個(gè)質(zhì)點(diǎn).計(jì)算域尺寸為6D0×6D0×15D0，在傘衣充氣區(qū)域進(jìn)行局部加密，總網(wǎng)格數(shù)約為2.6×106個(gè)，為保證質(zhì)量守恒，計(jì)算域邊界采用速度入口.

圖1 數(shù)值模擬中的降落傘折疊模型Fig.1 Parachute folding model in numerical simulation

3 結(jié)果與分析

3.1 流場分析

由于大型降落傘充氣過程中的劇烈變形可能導(dǎo)致傘衣破損或開傘失敗[16]，而HBG系統(tǒng)的電子設(shè)備質(zhì)量較小，因此可采用小型降落傘完成緩降功能.小型降落傘的剛度較大，其開傘過程與大型降落傘存在明顯區(qū)別.為了闡明HBG系統(tǒng)小型降落傘的充氣過程，首先以來流速度為100 m/s、電子設(shè)備質(zhì)量為4 kg的工況為代表，研究傘衣充氣過程的流場特性.

研究表明，傘衣外形變化與剪切層和近尾跡區(qū)渦結(jié)構(gòu)等的演化密切相關(guān)[14].圖2為來流速度為100 m/s時(shí)，傘衣充氣過程的外形變化與壓力云圖，采用速度等值線的方法描繪出剪切層的位置和形態(tài).

圖2 傘衣對稱面的壓力云圖Fig.2 Pressure contour of symmetrical plane of parachute

由圖2（a）可以看出，當(dāng)傘衣呈“燈泡”狀時(shí)，剪切層在傘衣底邊形成并緊貼于傘衣側(cè)面.當(dāng)剪切層發(fā)展到傘衣肩部時(shí)，傘衣發(fā)生分離并迅速失穩(wěn)（圖2（a））.隨著傘衣擴(kuò)張，剪切層的分離點(diǎn)逐漸向傘衣底邊移動(dòng)（圖2（a）～圖2（d））.同時(shí)，剪切層失穩(wěn)后卷起的尾渦也在向傘衣外側(cè)移動(dòng).在尾渦由中心向外側(cè)移動(dòng)時(shí)，傘衣頂部出現(xiàn)了高壓區(qū)（圖2（d）），說明傘衣在充氣過程中出現(xiàn)了尾流再附現(xiàn)象.尾流再附造成傘衣在傘頂孔附近出現(xiàn)非對稱的局部塌陷，加劇了傘衣近尾跡區(qū)流場的非對稱特性，導(dǎo)致傘衣沿徑向發(fā)生偏轉(zhuǎn)（圖2（f））.

為了進(jìn)一步研究傘衣充氣過程中流場演化的動(dòng)力學(xué)特征，采用Lamb矢量散度描述傘衣附近流場特性[15].Lamb矢量散度定義為Δ·L，其中L=ω×u為Lamb矢量，ω為偽渦矢量，u為流體速度矢量.圖3為流向速度為100 m/s時(shí)利用Δ·L描述的流場拓?fù)浣Y(jié)構(gòu).由圖3（a）～圖3（c）可知，在傘衣側(cè)面的剪切層中，Δ·L出現(xiàn)了明顯的正負(fù)雙層結(jié)構(gòu).該結(jié)構(gòu)的動(dòng)力學(xué)過程包含了強(qiáng)應(yīng)變率區(qū)和渦量區(qū)間的相互作用[15].尾流再附發(fā)生后，傘衣側(cè)面的Δ·L雙層結(jié)構(gòu)演化為正負(fù)交替結(jié)構(gòu)（圖3（d）～圖3（f））.該正負(fù)交替結(jié)構(gòu)一般出現(xiàn)在鈍頭體繞流的近尾跡區(qū)中，說明此時(shí)剪切層在傘衣底邊形成后不能穩(wěn)定地向下游發(fā)展，而是會(huì)迅速失穩(wěn).這種流場演化現(xiàn)象與無限質(zhì)量充氣過程明顯不同.在無限質(zhì)量充氣過程中，當(dāng)傘衣投影面積達(dá)到峰值后，剪切層和尾渦會(huì)進(jìn)入周期性演化[17].而在HBG系統(tǒng)小型降落傘的有限質(zhì)量充氣過程中，傘衣側(cè)面剪切層不會(huì)進(jìn)入周期性演化，而是出現(xiàn)劇烈的渦運(yùn)動(dòng)和流體拉伸變形運(yùn)動(dòng).

圖3 傘衣對稱面Lamb矢量的散度云圖Fig.3 Lamb vector divergence contour of symmetrical plane of parachute

HBG系統(tǒng)的降落傘具有前置體質(zhì)量小、傘衣面積小的特點(diǎn)，其開傘過程屬于有限質(zhì)量充氣.Desabrais[17]從流場特性演化的角度將降落傘充氣過程劃分為3個(gè)階段.與Desabrais[17]研究的無限質(zhì)量充氣過程相比，HBG系統(tǒng)降落傘充氣過程前2個(gè)階段的流場特性和傘衣外形變化特征是非常相似的.而在無限質(zhì)量充氣的第3個(gè)階段，傘衣側(cè)面的剪切層開始周期性演化，傘衣產(chǎn)生周期性的“呼吸”現(xiàn)象.與此不同，HBG系統(tǒng)降落傘在第3階段出現(xiàn)了尾流再附，導(dǎo)致降落傘在充氣過程的后2個(gè)階段中，剪切層不會(huì)進(jìn)入周期性演化，流場的非對稱特性也比無限質(zhì)量充氣更明顯.

圖4為流向速度為100 m/s時(shí)，HBG系統(tǒng)降落傘充氣過程中動(dòng)載和傘衣投影面積的變化曲線.

圖4 傘衣充氣過程中的動(dòng)載與投影面積Fig.4 Load and projected area of parachute during inflation

由圖4可以看出，根據(jù)流場演化的特點(diǎn)，充氣過程可以被分為4個(gè)階段.第1個(gè)階段為0～0.060 s，傘繩和傘衣從折疊狀態(tài)拉直，傘繩的拉直力使動(dòng)載出現(xiàn)第1個(gè)峰值；第2個(gè)階段為0.060～0.080 s，傘衣底邊首先開始充氣擴(kuò)張，隨著氣流進(jìn)入傘衣并在傘頂聚集形成高壓區(qū)，傘衣頂部開始擴(kuò)張出現(xiàn)“燈泡”狀；第3個(gè)階段為0.080～0.130 s，此階段內(nèi)傘衣投影面積迅速增大，開傘動(dòng)載出現(xiàn)峰值，隨后由于尾流再附，傘衣頂部出現(xiàn)局部塌陷；第4個(gè)階段為0.130～0.165 s，傘衣從傘頂局部塌陷的狀態(tài)逐漸回彈，在非對稱流場的影響下，傘衣在沿軸向運(yùn)動(dòng)的同時(shí)沿徑向發(fā)生偏轉(zhuǎn).在大型降落傘的空投實(shí)驗(yàn)中，傘頂在塌陷后的回彈過程中，可能發(fā)生傘衣纏繞造成嚴(yán)重破損[11].由于HBG系統(tǒng)降落傘尺寸較小，雖然傘頂在充氣過程中發(fā)生了塌陷和回彈，但并沒有纏繞其他部分傘衣.

3.2 前置體質(zhì)量和來流速度對降落傘充氣過程的影響

由前文分析可知，在典型工況下，HBG系統(tǒng)的小型降落傘在回彈過程中不會(huì)發(fā)生傘衣纏繞.為了進(jìn)一步分析該小型降落傘在充氣過程中的安全性，需要給出前置體質(zhì)量和來流速度等因素對傘衣充氣過程的影響.根據(jù)環(huán)縫傘相關(guān)研究[18]可知，尾流再附現(xiàn)象與物傘系統(tǒng)的質(zhì)量比和速度比相關(guān)，質(zhì)量比和速度比分別為

式（6）和式（7）中：m為前置體質(zhì)量；ρ為流體密度；D0為傘衣名義直徑；U∞為來流速度；UF為緩降速度.

現(xiàn)役飛行數(shù)據(jù)記錄器質(zhì)量約為6～12 kg，在HBG系統(tǒng)中，離機(jī)緩降的電子設(shè)備質(zhì)量小于傳統(tǒng)飛行數(shù)據(jù)記錄器.電子設(shè)備在飛機(jī)墜毀前彈射離機(jī)，假設(shè)來流密度為1.225 kg/m3不變，當(dāng)來流速度為100 m/s（UR=0.5）、傘衣名義直徑為1.9 m時(shí)，改變電子設(shè)備質(zhì)量，分別取質(zhì)量比為0.2、0.5、0.8、1.1和1.4的工況進(jìn)行分析.圖5為流向速度U∞=100 m/s、不同質(zhì)量比RN時(shí)，傘衣在充氣過程第3和第4階段內(nèi)的外形變化.

圖5 不同質(zhì)量比時(shí)，傘衣充氣過程中的外形變化Fig.5 Shape change of parachute during inflation with different mass ratio

由圖5可以看出，在0～0.120 s的時(shí)間內(nèi)，不同質(zhì)量比工況下，傘衣外形的變化規(guī)律相近.在傘衣充氣過程的第3階段，均出現(xiàn)了傘頂?shù)木植克?當(dāng)RM=0.2時(shí)，頂部局部塌陷后，傘衣塌陷的范圍逐漸變大，最終導(dǎo)致全部傘衣劇烈變形.當(dāng)RM=1.4時(shí)，頂部局部塌陷后，傘衣回彈，沒有出現(xiàn)傘衣完全塌陷的現(xiàn)象.

流向速度U∞=100 m/s時(shí)，開傘過程中的傘衣投影面積如圖6所示.

圖6 質(zhì)量比對傘衣投影面積的影響Fig.6 Effect of mass ratio on projected area of parachute

開傘過程的前0.070 s是傘衣傘繩的拉直過程，由于傘衣底面充氣口面積較小，傘衣充氣速度較慢，導(dǎo)致傘衣投影面積變化較小.隨后傘衣充氣過程進(jìn)入第2個(gè)階段，大量空氣沖入傘衣，造成傘衣投影面積迅速增加.在充氣過程的最后階段，當(dāng)RM=0.2時(shí)，傘衣大面積塌陷使其投影面積明顯大于其他工況.由圖6可以看出，當(dāng)RM=0.2時(shí)，傘衣投影面積峰值明顯大于其他工況.當(dāng)RM=0.5～1.4時(shí)，傘衣投影面積的變化曲線相近，說明在0.5＜RM＜1.4工況中，不會(huì)由于尾流再附導(dǎo)致傘衣大面積塌陷，是較為安全的工況.

為了進(jìn)一步說明前置體質(zhì)量不同時(shí)，傘衣充氣過程的流場特性，圖7給出了流向速度U∞=100 m/s時(shí)，傘衣充氣過程第3和第4階段典型時(shí)刻的壓力云圖.由圖7可以看出，t=0.100 s時(shí)傘衣附近流場特性相近，剪切層在傘衣底邊形成后向下游發(fā)展.在傘衣的近尾跡區(qū)中，剪切層失穩(wěn)并卷起旋渦.隨著傘衣充氣過程向第4階段發(fā)展，傘衣附近流場區(qū)別越來越明顯.在RM=0.2的工況中，由于前置體質(zhì)量較小，在降落傘氣動(dòng)力作用下，前置體向傘衣方向回彈，導(dǎo)致t=0.140 s時(shí)部分傘繩處于松弛狀態(tài)（圖7（a））.由于傘繩對傘衣的拉扯作用減小，傘衣難以保持類似半球殼狀的外形.對比RM=1.4、t=0.140 s時(shí)傘衣對稱面的壓力云圖（圖7（b））可以發(fā)現(xiàn)，RM=0.2時(shí)傘前氣體的高壓區(qū)域幾乎消失，說明此時(shí)傘衣對于來流氣體的滯止效應(yīng)非常微弱.

圖7 不同質(zhì)量比傘衣對稱面的壓力云圖Fig.7 Pressure contour of symmetrical plane of parachute with different mass ratio

本研究中來流速度為40～160 m/s，緩降速度為8 m/s，當(dāng)前置體質(zhì)量為4 kg（RM=0.5）時(shí)，分別取速度比RV為5.00、8.75、12.50、16.25和20.00的工況進(jìn)行分析.圖8中給出了速度比RV為5.00和20.00的工況下，傘衣在充氣過程第3和第4階段的外形變化情況.圖9進(jìn)一步給出了速度比RV為5.00、8.75、12.50、16.25和20.00的工況下，開傘過程中的動(dòng)載和傘衣投影面積變化曲線.

圖8 不同速度比時(shí)，傘衣充氣過程中的外形變化Fig.8 Shape change of parachute during inflation in different speed ratio

由圖8可以看出，在不同速度比的工況下，傘衣在充氣過程中均經(jīng)歷了由充滿到傘頂塌陷再到傘頂回彈的外形變化過程.在傘衣充氣時(shí)，沒有出現(xiàn)傘衣完全塌陷的現(xiàn)象.為了方便對比分析，將動(dòng)載峰值時(shí)刻定義為傘衣充滿時(shí)間，采用無量綱時(shí)間t/tf繪制傘衣投影面積變化曲線（圖9（b））.由圖9（a）可以看到，傘衣充氣時(shí)間隨著來流速度的增加而降低，同時(shí)動(dòng)載峰值隨著來流速度的增加而上升.由圖9（b）可知，速度比不同時(shí)，傘衣投影面積的變化較為相似，說明在RM=0.5、5.00＜RV＜20.00工況下，傘衣在充氣過程中均未出現(xiàn)大范圍的傘衣塌陷現(xiàn)象，是較為安全的工況.

圖9 速度比對開傘動(dòng)載和傘衣投影面積的影響Fig.9 Effect of speed ratio on opening load and projected area of parachute

4 結(jié)論

本文采用任意拉格朗日-歐拉（ALE）方法數(shù)值模擬了HBG系統(tǒng)降落傘的充氣過程，分析了傘衣充氣過程的流場特性，研究了前置體質(zhì)量和來流速度對傘衣塌陷的影響，研究結(jié)果表明：

（1）根據(jù)流場演化的特點(diǎn)，傘衣的充氣過程可分為4個(gè)階段.其中，前2個(gè)階段的流場特性和傘衣外形變化特征與無限質(zhì)量充氣過程非常相似.在充氣過程的后2個(gè)階段，剪切層不會(huì)進(jìn)入周期性演化，流場的非對稱特性也比無限質(zhì)量充氣更明顯.

（2）傘衣在充氣過程中出現(xiàn)尾流再附，并導(dǎo)致傘頂局部塌陷.若前置體質(zhì)量過小，傘衣塌陷的范圍會(huì)逐漸變大，并最終導(dǎo)致全部傘衣劇烈變形，是較為危險(xiǎn)的工況.當(dāng)0.5＜RM＜1.4時(shí)，傘衣頂部局部塌陷后傘衣回彈，不會(huì)由于尾流再附導(dǎo)致傘衣大面積塌陷，開傘失敗的風(fēng)險(xiǎn)較小，是較為安全的工況.

（3）在較大的工況范圍內(nèi)，當(dāng)來流速度不同時(shí)，傘衣投影面積的變化規(guī)律較為相似.當(dāng)5.00＜RV＜20.00時(shí)，傘衣在充氣過程中均未出現(xiàn)大范圍的傘衣塌陷現(xiàn)象，開傘失敗的風(fēng)險(xiǎn)較小，是較為安全的工況.