不同燃氣流量對空間氣囊充氣展開過程的影響

信志濤，王 浩，陶如意

(南京理工大學 能源與動力工程學院， 南京 210094)

空間氣囊作為一種充氣膜結構，在空間天線、太陽帆板等大型空間結構上均得到了廣泛應用，并且成為未來大型空間結構應用的一種趨勢[1]。目前，國內外對汽車用安全氣囊、船舶用氣囊和航空航天中使用的緩沖氣囊等類型的氣囊實驗和理論研究比較成熟[2-6]。

某空間氣囊采用火藥在燃氣發生器內燃燒產生燃氣作為動力源，燃氣流場的流動特性和溫度變化特性是影響氣囊應力和應變的關鍵因素，囊內流場的形成和變化又和燃氣發生器內的燃氣狀態密不可分，因此將內彈道計算模型和有限元分析方法相結合，為空間氣囊充氣展開的研究提供有效手段。

為驗證不同的燃氣流量對空間氣囊充氣展開過程影響，本文設計了空間氣囊充氣展開系統，通過不同的裝藥量在燃氣發生器內燃燒產生不同的燃氣流量作為氣囊展開的動力源，分析對氣囊展開過程的影響。

1 燃氣發生器內彈道過程數值模擬

某空間氣囊展開系統由點火系統、燃氣發生系統及氣囊等組成，整個裝置如圖1所示。點火藥采用2#小粒黑，主裝藥為3/1樟火藥，所用火藥為顆粒狀制式火藥，采用經典內彈道計算時不需要考慮裝藥結構設計，不考慮填充過程，計算只需裝藥量和燃氣發生器參數。氣囊展開過程以點火藥被點燃為起點，通常在充氣30 ms內氣囊充氣展開完全并趨于穩定，因此以30 ms為計算終點，期間包含了點火藥的點燃、主裝藥的燃燒、拋放彈破膜、火藥燃氣流動、氣囊膨脹等復雜變化過程。

1.1 基本假設

根據經典內彈道理論及設計方案，針對燃氣發生系統內的火藥燃燒、燃氣流動及氣囊充氣膨脹等過程，在建立內彈道數理模型的過程中有以下基本假設[7-8]：

1) 火藥燃燒模型：點火具與點火藥瞬間燃完，同時主裝藥瞬間全面點燃，火藥的燃燒過程遵循幾何燃燒定律，燃燒速度與燃氣發生室內的平均壓力呈指數函數關系；

2) 火藥燃氣模型：火藥燃燒產生組份不變，火藥燃氣服從諾貝爾-阿貝爾狀態方程，有關參數如火藥力f，氣體絕熱指數k=Cp/Cv，余容α等為常數；火藥燃氣流入氣囊內部的過程視為等熵流動，噴嘴處火藥燃氣流速為當地音速；

3) 燃氣發生系統模型：當拋放彈內壓力達到破膜壓力后，限壓膜片瞬間打開，火藥燃氣及未燃完的火藥顆粒瞬時充滿并均勻分布于整個燃氣發生器。

1.2 內彈道過程控制方程

基于以上基本假設和內彈道學理論，燃氣發生器內彈道過程劃分為兩個階段：① 從點火藥點燃到拋放彈內限壓膜片破裂瞬間，認為點火藥瞬間燃完并點燃主裝藥，這一階段為火藥在拋放彈內定容燃燒；② 從拋放彈限壓膜片破裂到火藥燃燒結束，這一階段包含了火藥燃燒、燃氣流動、氣囊充氣膨脹等過程，是最為復雜的階段。

依據以上各階段的劃分，可建立內彈道方程組[9]：

1) 燃速方程

(1)

式中：Z為已燃厚度百分比；u1為燃速系數；p為燃氣發生器內壓力；n為燃速指數；e1為火藥弧厚的一半。

2) 形狀函數

(2)

式中：ψ為火藥已燃質量百分比；χ、λ、μ為火藥形狀特征量；Zk為火藥已燃相對厚度。

3) 燃氣發生器內火藥氣體狀態方程

(3)

式中：f為主裝藥火藥力；τ=T/T1，其中T為燃燒室溫度，T1為火藥爆燃溫度；ω為主裝藥質量；f1為點火藥力；ω1為點火藥質量；V0為拋放彈藥室容積；V1為燃氣發生器容積；ρp為火藥密度；α為主裝藥氣體余容；η為燃氣發生器內火藥氣體流入氣囊的相對流出量；α1為點火藥氣體余容。

4) 氣體流量方程

(4)

式中：st為燃氣發生器內噴孔面積；φ為流量損耗系數，k為火藥燃氣比熱比。

5) 燃氣發生器能量守恒方程

(5)

式中，θ=k-1。

1.3 計算結果

根據前面建立的模型，通過數值計算得出裝藥量為4 g、4.5 g和5 g時燃氣發生器噴嘴質量流量和燃燒室溫度曲線，如圖2和圖3所示，從圖中可以看出，裝藥量越大，噴嘴流量上升越快，因為裝藥越大燃燒室壓力上升越快，火藥燃燒也越快；破膜后燃燒室溫度隨著火藥流出迅速下降，之后隨著火藥燃燒加快，燃燒室溫度下降較慢。在裝藥量為5 g時，噴嘴流量和燃燒室溫度在27.05 ms時刻都經歷了迅速下降的階段，因為此時火藥已燃盡，由燃氣發生器能量守恒方程和氣體狀態方程可知，此時燃氣的流出會導致燃氣發生器溫度和壓力迅速下降。30 ms時4 g、4.5 g和5 g裝藥量下，流入氣囊的燃氣總質量分別為2.31 g、3.48 g和4.73 g。計算結果將代入LS-DYNA軟件作為氣囊展開計算的入口邊界條件。

2 氣囊展開數值理論

控制容積法(CV法)也稱為均壓法(Uniform Pressure Method)，在充氣結構展開數值計算中應用比較多，如圖4所示的氣囊模型，首先把氣囊看作是一個由表面膜單元圍成的不斷擴大的控制容積，燃氣發生器及其囊內火藥燃氣壓力遠遠低于600 MPa，可以視為理想氣體，且認為燃氣的比熱容不變，囊內各處壓力處處相等，在不考慮氣囊與外界傳熱的前提下，囊內火藥燃氣滿足下列關系式

(6)

式中：P為囊內壓力；V為囊內體積；m為流入氣囊的火藥燃氣質量；M為火藥燃氣的摩爾分子質量；e為比內能；ρ為火藥燃氣的密度；CV為火藥燃氣的定壓比熱容。

3 氣囊有限元模型

采用自底向上的建模方式，在LS-PREPOST中建立如圖5所示的氣囊未折疊模型，該氣囊模型由上下兩片囊壁和底部的噴嘴組成，這里不需要對燃氣發生器進行建模，只需要建立噴嘴的一部分模型即可。氣囊直徑為800 mm，由上下兩片間距為1 mm的圓形織物囊壁通過共節點方式連接而成，氣囊與噴嘴之間也是通過共節點方式連接，模型共有 10 289個節點，10 152個膜單元。采用系統自帶的折疊工具進行折疊，對其進行折疊后的模型如圖6所示。

CV法氣囊模型的主要關鍵字如下：

CONTACT_AIRBAG_SINGLE_SURFACE：用于設置氣囊的自接觸；

*BOUNDARY_SPC_SET：用于噴嘴的固定；

*DEFINE_CRUVE：用于實現噴嘴流量和燃氣溫度曲線的輸入；

*DATABASE_ABSTATE：用于輸出氣囊內部壓力、體積等數據；

*AIRBAG_WANG_NEFSKE：用于實現氣囊模型的定義，確定采用的計算方法為CV法，可以輸入環境參數及燃氣數據曲線編號，實現對氣囊的計算；

*CONTROL_CONTACT：用于對接觸進行控制，IGNORE設為1，可以消除初始穿透對氣囊計算的影響。

氣囊采用織物材料，不考慮囊內氣體的泄漏，不考慮燃氣與氣囊壁之間的導熱，氣囊材料為各向同性材料，單元算法采用Belytschko-Tsay全積分膜單元，氣囊的材料參數如表1所示。

表1 氣囊的材料參數

4 模擬結果及分析

4.1 氣囊動力學分析

不同裝藥量條件下氣囊展開過程如圖7、圖8和圖9所示，應力單位為GPa，分別取氣囊展開過程中3 ms、10 ms、20 ms和30 ms時刻的展開圖形對比分析，為了便于觀察氣囊應力分布，尤其是噴嘴附近的應力分布，考慮到折疊氣囊結構的對稱性，選取氣囊一半的模型分析應力和氣囊展開過程。從中可以看出，3 ms時氣囊均處于展開初期，應力集中在噴嘴附近，應力相對較小。10 ms時氣囊均處于展開中期，氣囊的展開程度隨著裝藥量的增加而增加，此時氣囊受到的應力集中在噴嘴附近。裝藥量為4 g、4.5 g和5 g時的最大應力分別為60.5 MPa、113.6 MPa和89.7 MPa。20 ms時不同裝藥量的氣囊均已基本展開，4 g裝藥量條件下氣囊仍有部分褶皺未展開，最大應力位于氣囊褶皺處；4.5 g裝藥量條件下氣囊已展開完全，隨著褶皺的展開，應力主要集中在噴嘴處；5 g裝藥量下氣囊完全展開，應力集中分布于氣囊噴嘴處，隨著距離噴嘴距離增加而減小。30 ms時氣囊均已完全展開，應力集中于噴嘴和展開軸線垂直的方向上，但4 g裝藥量下有一處褶皺未展開，導致褶皺處應力較大。完全展開時氣囊的最大應力呈現隨裝藥量增大而增大的趨勢，從圖形中可以看出此時最大應力又隨著氣囊膨脹程度增大而減小。從氣囊全部展開過程看，4 g裝藥量條件下氣囊展開過程中最大應力為200.2 MPa，位于氣囊褶皺處；4.5 g裝藥量時最大應力為215.4 MPa，位于氣囊噴嘴連接處；5 g裝藥量時最大應力為228.0 MPa，位于氣囊噴嘴連接處。

噴嘴正上方節點ID為7 685點，該節點的運動可以反應氣囊的外形變化，其沿噴嘴軸線位移曲線如圖10所示，從圖10可以看出，7 685點先是隨著氣囊膨脹沿噴嘴軸線向上運動，到達頂端后隨著氣囊的回彈又向下運動。裝藥量越大該點沿噴嘴軸線位移越大，并且達到頂端用時越少，說明裝藥量越大氣囊變形程度越大，且氣囊膨脹速度越快。由于噴嘴固定和氣囊織物材料的彈性影響，氣囊向頂端位移越大，受到的向下拉力越大，因此氣囊回彈越快。當氣囊回彈結束時，氣囊又會出現反方向回彈，如5 g裝藥量下氣囊在27.9 ms時出現了反方向的回彈，由于氣囊已完全展開，回彈速度變慢。

4.2 氣囊容積及壓力分析

不同裝藥量下氣囊展開容積變化曲線如圖11所示，可以看出氣囊展開過程中容積隨著燃氣流量的增加而增加，達到最大容積后小幅波動，且裝藥量越多，氣囊容積越大，這是由于氣囊的織物材料有一定的彈性，裝藥量越大囊內壓力越大，織物材料拉伸越大。4 g裝藥量時氣囊容積在23.6 ms時接近最大值8.35×107mm3，4.5 g裝藥量時氣囊容積在20.2 ms時接近了最大值8.56×107mm3，5 g裝藥量時氣囊容積在18.1 ms時接近了最大值8.98×107mm3。當氣囊完全展開后，其容積一直處于波動狀態，這是由于氣囊由向上運動時受到噴嘴限制變為向下運動，向下運動到一定程度又會受到噴嘴的限制變為向上運動。當氣囊改變運動方向時，會對燃氣有一定的壓縮作用，導致氣囊容積處于波動狀態。

不同裝藥量下氣囊壓力變化曲線如圖12所示，從圖12可以看出，在氣囊展開初期，囊內壓力有一個躍升階段，因為此時氣囊從靜止到充氣展開，體積增加較慢。之后囊內壓力又迅速增加，在約7 ms左右達到各自的峰值，這個階段壓力上升主要是因為氣囊的折疊結構使氣囊容積上升較慢，而燃氣流入不受影響。隨著氣囊的展開，囊內壓力有所下降，等氣囊完全展開后，氣囊的容積變化很小，而燃氣不斷流入，壓力不斷升高。4 g、4.5 g和5 g裝藥量下囊內壓力的峰值出現在30 ms計算結束時刻，并且隨著裝藥量的增加，囊內最大壓力也隨之增加。

5 結論

1) 可以根據燃氣發生器經典內彈道數理模型，理論計算得到不同裝藥量下的噴嘴流量曲線和燃燒室溫度曲線。

2) 采用CV法可以模擬真空下氣囊展開過程，但由于CV法原理的限制，與實際展開過程不完全相同。

3) 燃氣流量越大，氣囊受到的應力越大，展開越快，展開后容積越大，囊內壓強越大。調整裝藥量可以改變燃氣流量，從而改變氣囊展開進程，可以根據裝藥量合理選擇氣囊材料。

4) 氣囊完全展開后，由于受到噴嘴的限制以及氣囊材料的影響，會沿噴嘴軸線出現反復回彈運動。

[1] 韓克良,關富玲,曹莉,等.充氣膜結構展開過程仿真分析[J].建筑結構學報,2010(S1):288-292.

[2] 徐靜靜.基于流—固耦合模型的氣囊織物動態力學仿真研究[D].上海：東華大學,2010.

[3] 衛劍征.空間折疊薄膜管充氣展開過程氣固耦合問題研究[D].哈爾濱：哈爾濱工業大學,2008.

[4] GAO X Y,ZHANG P.Testing of the Dynamic Permeability and Deformation of Airbag Fabrics[J].Advanced Materials Research,2014,941-944:1345-1349.

[5] NIE M,XU C,LI T,et al.Design and Implementation of Development and Testing Platform for Airbag Control System[J].Applied Mechanics and Materials,2013,321-324:1554-1557.

[6] SHAO,ZHIJIAN,CHENG,et al.Multi-objective Optimization Design of Vented Cylindrical Airbag Cushioning System for Unmanned Aerial Vehicle[J].Transactions of Nanjing University of Aeronautics and Astronautics,2016,33(2):208-214.

[7] 李嗄.氣囊式子母彈拋撒過程數值模擬[D].南京：南京理工大學,2012.

[8] 王浩.子母彈內燃式氣囊拋撒模型及計算機仿真[J].兵工學報,2001,22(2):178-181.

[9] 張小兵.槍炮內彈道學[M].北京：北京理工大學出版社,2014.

[10] 陳東雷.火藥氣體成分與簡化狀態方程[J].南京理工大學學報,1987(1):37-46.

[11] 賀增弟，蔡鎖章.火藥燃燒平衡組成計算中的最優化數學模型[J].華北工學院學報,2003(5):332-335.