自落式緩沖氣囊仿真研究

周清艷 言金

（中航工業航宇救生裝備有限公司，襄陽 441003）

0 引言

最早的氣囊裝置是在20世紀40年代因飛機安全保護提出[1]，氣囊保護技術被廣泛應用于國防軍事領域、航天器著陸與回收、武器裝備與物資空投、汽車安全氣囊等領域[2-3]。隨著對緩沖效果的要求提高，結構更復雜、材料更先進、效果更可靠的氣囊正在為大家所關注和研究。

緩沖氣囊種類較多，按充氣方式不同可分為主動充氣式和自然充氣式；按結構形式可分為單氣囊和組合式氣囊；按是否放氣可分為封閉式和排氣式，排氣式又可分為恒定排氣孔式和可控排氣孔式。

對氣囊的緩沖研究主要有兩種：一類是解析分析方法，文獻[4]建立的模型假設氣囊壓縮為等溫過程，給出了重裝空投用自落式緩沖氣囊的一般設計思路和計算方法；另一類是數值模擬即動態非線性有限元仿真[5-6]，其優勢在于可以進行緩沖氣囊的全向緩沖特性分析，模擬各種著陸環境，適用于計算各種柔性織物的結構變形，對其工作過程中的速度、加速度、動量、能量及應力應變的分布情況都能夠有效地表示出來。

國內對于氣囊工作過程的模擬仿真研究已經逐漸發展起來[7-10]，目前較為流行的有控制體積法（Control Volume Method, CV）和流固耦合法（Arbitrary Lagrangian Eulerian, ALE）。CV法研究氣囊緩沖過程，可以得到較為理想的氣囊外形變化情況，計算時間短，占用計算機資源少，大大降低了計算成本。ALE法尚在探索和改進中，相比CV法而言，它可以計算并獲取囊內流場數據，在計算氣囊、降落傘等高速充氣過程中非常有必要[11-12]。因此，氣囊的仿真研究有著很大的發展空間，有待于進一步的努力。

本文針對某空投貨臺緩沖系統開展研究，重點分析氣囊緩沖工作過程。采用CV法，運用LS-DYNA軟件進行仿真計算，結合工程估算，為研究氣囊工作性能提供理論依據，為氣囊設計提供參考。

1 模型建立

本文以某型空投緩沖系統所用氣囊為研究對象，該緩沖氣囊系統由8個獨立的氣囊并聯組成。每個獨立氣囊可以分為主囊和輔囊2部分。主囊底部設有進氣口，兩側面設有排氣口。輔囊縫合在主囊斜面上，通過排口與主囊相通?？胀对O備與主囊上表面相連，為研究方便，本文以同等質量鋼板代替空投設備。

主氣囊依靠底部進氣孔自落充氣，當氣囊與地面接觸后，進氣口被封閉（假設完全封閉）。主氣囊受壓縮后，一方面通過主囊側面排氣口向外排氣，另一方面開始向輔囊充氣，當輔囊內壓超過搭扣帶的粘合力，輔囊排氣口開啟并向外排氣。氣囊卸壓緩沖過程中，空投設備只與主氣囊上表面接觸，輔囊壓縮空氣起延緩放氣作用，并不直接對空投設備施加作用力。

為方便建模和計算，對緩沖系統做如下假設：

1）氣囊為柔性織物，但忽略織物的彈性，在壓縮過程中不產生變形；

2）氣囊的緩沖作用是由囊內氣體壓力的損失來提供能量的，不考慮氣動阻力；

3）緩沖過程空氣流動視為不可壓縮流，為絕熱過程；

4）輔囊搭扣帶之間的粘合力極小，可以假設其排氣孔的爆破壓力與環境大氣壓一致。

氣囊在貨臺橫向方向所占據的尺寸為2 360mm，縱向方向為4 500mm，氣囊高度1 440mm。空投設備質量為7 600kg。由于每個獨立氣囊的排氣口之間無影響，在忽略2個氣囊在緩沖過程中的相互擠壓（試驗顯示2個氣囊之間的擠壓不是很明顯）造成的影響時，可針對其中1個獨立氣囊進行建模并仿真計算。

其邊界條件見表1，幾何尺寸如圖1。

表1 初始邊界條件Tab.1 The initial boundary conditions

圖1 幾何模型Fig.1 The geometrical model of the airbag

2 計算結果與分析

將幾何模型導入有限元分析前處理軟件Hypermesh中，添加相應關鍵字，設置載荷分布，同時，定義氣囊邊界、氣體性質，并添加空投設備與氣囊的連接，進行網格劃分，形成K文件并導出，利用LS-DYNA對上述K文件進行數值計算，運用LS-Prepost軟件對計算結果進行后處理，得到氣囊外形及囊體表面應力仿真結果，見圖2。

圖2 不同時刻氣囊外形變化及表面應力分布云圖Fig.2 Shape changes and stress distribution of the airbag in different time

圖2分別為緩沖氣囊著陸時刻、氣囊囊體所受最大應力時刻和緩沖結束時刻的計算狀態，整個緩沖過程時間大約為0.25s。分析整個緩沖過程應力分布情況，在t=0.118s時主囊外排氣口處出現應力最大值。除排氣口出現應力集中外，氣囊表面應力分布比較均勻。相對主囊而言，輔囊表面應力較小且分布比較平穩。

針對以上計算結果，下面就主囊及輔囊內氣體體積變化、空投設備系統速度、加速度變化情況分別進行分析。

2.1 囊內氣體體積變化情況

分別對主囊、輔囊內氣體體積變化進行分析，得到其變化曲線如圖3和圖4所示。

圖3 主囊內氣體體積變化Fig.3 The air volume change in the main airbag

圖4 輔囊內氣體體積變化Fig.4 The air volume change in the subsidiary airbag

從氣囊著陸開始（t=0.02s），主囊被壓縮，同時向輔囊進行排氣，輔囊搭扣被撕裂，進而打開排氣口，整個系統開始緩沖，囊內氣體體積減少。由于存在微弱反彈，氣體未全部排空。從圖4可以看出，輔囊的收縮現象（呼吸現象）比較明顯，這也體現了其在整個系統中的作用，即延遲囊內氣體排出，穩定緩沖氣囊系統。

2.2 空投設備速度變化

圖5為空投設備速度變化曲線。從圖中可以看出，在緩沖初始階段（0.1s之前）由于空投設備重力作用大于氣囊對其反作用力，空投設備有微量增速現象。在t=0.1s之后，氣囊反作用力大于空投設備重力作用，系統開始迅速減速。在t=0.24s左右，速度降為0，之后Z向速度為正，表明此時系統出現反彈現象，但由于此前氣囊被壓縮，沒有出現離地反彈，僅為氣囊向上膨脹。

圖5 空投設備速度變化曲線Fig.5 The velocity change of the cargo system

2.3 空投設備加速度（過載）變化

圖6為空投設備Z向加速度變化曲線，從圖中可以看出，緩沖過程初始階段（0.1s之前），空投設備由于其重力作用大于氣囊對其反作用力，初始加速度為Z軸相反方向，并逐漸減小。隨著氣囊進一步壓縮，對空投設備形成逐漸增大的反作用力，空投設備加速度開始沿著Z軸正方向逐漸增大。在t=0.21s左右，達到最大值。此時，空投設備所受最大過載為10.8gn。

圖6 空投設備加速度變化Fig.6 The acceleration change of the cargo system

2.4 與試驗結果、工程計算的對比

為更好驗證此氣囊緩沖系統性能，對其進行了空投沖擊試驗，測量系統緩沖過程中空投設備的速度、過載及位移變化情況。同時，基于熱力學原理，對緩沖氣囊建立工程計算模型，設計氣囊緩沖過程計算軟件，對該系統緩沖工作過程進行了工程計算。

如圖7所示，就以上3種研究方法所得空投設備速度變化情況進行對比與分析。結果表明，工程計算結果與試驗所得結果相比，趨勢基本保持一致，但誤差較大。這也是由于計算時對輔囊的簡化假設，沒有考慮主囊向輔囊的排氣過程；且受迭代計算的限制，它不能計算氣囊反彈過程。數值仿真計算與試驗測量值相比誤差較小，且其計算趨勢與試驗測量值基本一致，能有效計算氣囊緩沖過程中出現的反彈，可通過對模型進一步優化，減少計算誤差。

圖7 兩種計算與試驗所得速度對比Fig.7 The comparison between the calculations and tests

3 結束語

通過對該型緩沖氣囊工作過程的仿真與分析，結合工程計算與試驗測試，驗證了此氣囊設計的合理性。同時，也提出一種全新的自落式緩沖氣囊設計思路，即在綜合考慮緩沖效果、材料特性、空投裝備性質等因素下，根據需要運用工程計算確定氣囊初步設計參數，再由仿真計算模擬氣囊緩沖過程，結合沖擊試驗，可為氣囊的設計提供合理可行的研制思路與方法，極大提高了研制品質和效率，降低研制成本。

（

）

[1] 鐘志華, 張維剛, 曹立波, 等.汽車碰撞安全技術[M].北京: 機械工業出版社, 2005.

ZHONG Zhihua, ZHANG Weigang, CAO Libo, et al.The Safety Technology of Automobile Crash [M].Beijing: Machanical Industry Press, 2005.（in Chinese）

[2] 賈賀, 榮偉.火星探測器減速著陸技術分析[J].航天返回與遙感, 2010, 31（3）: 6-14.

JIA He, RONG Wei.Mars Exploration Deceleration Landing Technology Analysis[J].Spacecraft Recovery & Remote Sensing,2010, 31（3）: 6-14.（in Chinese）

[3] 榮偉, 陳國良.火星探測器減速著陸技術特點[J].航天返回與遙感, 2010, 31（4）: 1-6.

RONG Wei, CHEN Guoliang.The Characters of Deceleration and Landing Technology on Mars Explorer[J].Spacecraft Recovery & Remote Sensing, 2010, 31（4）: 1-6.（in Chinese）

[4] 孫曉偉.自落式緩沖氣囊設計計算的一般方法[J].中航救生, 2005（4）: 29-31.

SUN Xiaowei.Method of the Design Calculation for Self-dropping Cushion Airbag[J].China Aeronautical & Astronautical Life-support, 2005（4）: 29-31.（in Chinese）

[5] 衛劍征, 譚惠豐, 萬志敏, 等.緩沖氣囊展開與緩沖著陸過程的仿真分析[J].航天返回與遙感, 2010, 31（5）: 1-8.

WEI Jianzheng, TAN Huifeng, WAN Zhimin, et al.Simulation for Airbag Deployment and Landing Process of Inflatable Landing Vechiles[J].Spacecraft Recovery and Remote Sensing, 2010, 31（5）: 1-8.（in Chinese）

[6] 尹漢鋒, 文桂林, 韓旭.空投設備緩沖氣囊的優化設計[J].系統仿真學報, 2008, 24（5）: 36-40.

YING Hanfeng, WEN Guilin, HAN Xu.Optimal Design of Airbag Impact Attenuation System for Airdropping Equipment[J].Journal of System Simulation, 2008, 24（5）: 36-40.（in Chinese）

[7] 樂莉, 丁剛毅.空投武器設備氣囊緩沖過程的數值模擬[J].航空兵器, 2002（4）: 21-23.

LE li, DING Gangyi.Numerical Simulation of Airbag Impact Attenution System for Airdropping Equipment[J]．Aero Weaponry, 2002（4）: 21-23.（in Chinese）

[8] 何成, 陳國平, 何歡, 等.氣囊著陸緩沖系統的沖擊動力學多目標優化[J].航天返回與遙感, 2012, 33（5）: 1-8.

HE Cheng, CHEN Guoping, HE Huan, et al．Multi-objective Optimization of Impact Dynamics for Airbag Cushion Landing System[J].Spacecraft Recovery & Remote Sensing, 2012, 33（5）: 1-8.（in Chinese）

[9] 馬常亮, 何歡, 陳國平.氣囊著陸緩沖系統的排氣孔面積優化[J].航天返回與遙感, 2009, 30（4）: 6-10.

MA Changliang, HE Huan, CHEN Guoping.The Optimization of Venting Area of Airbag Landing System[J].Spacecraft Recovery and Remote Sensing, 2009, 30（4）: 6-10.（in Chinese）

[10] 王亞偉, 楊春信, 柯鵬.貨臺空投系統氣囊緩沖過程仿真[J].系統仿真學報, 2007, 19（14）: 3176-3179.

WANG Yawei, YANG Chunxin, KE Peng.Airbag Cushion Process Simulation for Cargo Airdrop System[J].Journal of System Simulation, 2007, 19（14）: 3176-3179.（in Chinese）

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

XU Jingjing.Simulation of Airbag Dynamic Mechanical Performance Based on Fluid-structure Interaction Approach[D].Shanghai: Donghua University, 2010.（in Chinese）

[12] 余莉, 明曉.降落傘的計算機仿真研究[J].航天返回與遙感, 2005, 26（4）: 6-9,49.

YU Li, MING Xiao.Investigation of Parachute’s Emulator and Computer Simulation[J].Spacecraft Recovery & Remote Sensing, 2005, 26（4）: 6-9,49.（in Chinese）