基于Sobol法的載人空降人體頭－頸部損傷響應特性研究

劉 鑫，張志勇，劉桂萍

(1.長沙理工大學 汽車與機械工程學院，長沙 410114;2.湖南大學 汽車車身先進設計制造國家重點實驗室，長沙 410082)

基于Sobol法的載人空降人體頭－頸部損傷響應特性研究

劉 鑫1，張志勇1，劉桂萍2

(1.長沙理工大學 汽車與機械工程學院，長沙 410114;2.湖南大學 汽車車身先進設計制造國家重點實驗室，長沙 410082)

建立了載人著陸沖擊的“假人－座椅”數值模型，并進行數值模擬和實驗驗證。模型經驗證后，其數值結果與實驗數據非常接近，并基于Sobol法對人體不同部位角度對人體頭－頸部損傷響應的影響進行了靈敏度分析。結果表明：人體頭－頸部角度和體位角對人體損傷響應有著重要影響。根據分析結果，可為最佳人體著陸姿勢的研究提供變量選取依據，在載人空降防護技術領域具有一定的實際工程意義。

著陸沖擊;人體損傷響應;Sobol法;靈敏度分析

載人空降過程中，乘員的著陸姿勢將直接影響到人體損傷的動態響應。目前，許多學者［1－4］已對人體著陸姿勢的問題進行了相關探討。Brown［5］針對Apollo飛船著陸的工況，對人體的24種不同著陸姿勢進行了相關實驗研究，為Apollo飛船的設計提供了相關科學依據;Weis［6］將20位志愿者以7種不同的體位角姿勢在垂直降落塔上做了人體沖擊載荷實驗，并歸納了人體對沖擊載荷的耐受區間;劉炳坤等［7］也對不同體位角下的人體沖擊響應做了相關研究。他們發現人體的損傷響應與沖擊加速度的峰值、持續時間以及人體的體位角度有關。在他們的實驗中，都只考慮了人體的體位角和加速度對人體損傷的影響。而在實際著陸沖擊過程中，人體的損傷響應并不僅僅只與體位角有關，人體每個部位的角度都對人體的損傷響應存在貢獻。只單一考慮某個部位雖能在模型上得到簡化，但無法全面評估人體著陸姿勢對人體損傷響應的影響。

本文通過建立人體損傷響應的數值模型，提出一種基于Sobol法的人體著陸姿勢變量靈敏度分析方法，研究乘員在空降過程中人體不同部位角度對人體損傷響應的影響。根據變量靈敏度分析結果，可為最佳人體著陸姿勢的研究提供變量選取依據，從而有利于提高載人空降乘員的安全性。

1 載人空降過程中的人體損傷響應評價

在載人空降過程中，乘員仰臥在防護座椅裝置里面，除了承受作用于人體且方向垂直向下的重力加速度場之外，還將承受到在地面沖擊載荷作用下造成的沿人體背—胸方向的加速度場。根據汽車碰撞安全損傷原則以及載人空降的實際情況，選擇損傷評價準則。

1.1 頭部損傷評價

著陸沖擊下的頭部傷害屬于閉合性腦損傷，主要有頭骨骨折和腦損傷(腦震蕩和腦挫裂傷)等。本文采用如下公式進行頭部損傷(head injury criterion，HIC)評價［8］：

式中：a(t)表示碰撞過程中頭部質心合成加速度;t2－t1表示HIC達到最大值時的時間間隔，在實際應用中最大時間間隔取36 ms;HIC值的耐受限度為1 000。

1.2 頸部損傷評價

頸部損傷準則NIC的基礎是脊髓神經根損傷機理和生物力學試驗。該準則基于頸部上下兩端的相對加速度和相對速度，其計算公式為：

其中，αrel(t)為頭部和胸部T1的相對加速度，Vrel(t)為頭部和胸部T1的相對速度。Bostom等［9］推薦NIC耐受限度值15來判斷頸部是否發生損傷。

2 人體損傷響應模型的數學描述

乘員由于受到地面沖擊載荷的作用，頭部和頸部極易發生損傷，這些部位的損傷值又和人體的著陸姿勢緊密相關。所以研究不同人體著陸姿勢對人體損傷響應的影響是非常必要的。需要說明的是，本文所研究的人體著陸姿勢是針對在同一著陸沖擊載荷作用下，暫不考慮安全帶對乘員胸部損傷的影響，選取頭部和頸部兩處損傷響應值作為評價人體著陸姿勢優劣的標準［10］。

考慮到不同的人體著陸姿勢對人體損傷響應將產生不同的影響，而影響人體著陸姿勢的參數主要有四個，如圖1所示，分別為：頭－頸角度、體位角度、大腿角度和小腿角度。所以選取這四個參數作為設計變量，其范圍如表1所示。

圖1 人體著陸姿勢Fig.1 Human landing position

表1 人體各部位角度范圍Tab.1 The bounds of the design variables

因此，載人著陸沖擊的人體損傷響應模型可用下式表示為：

式中：fH(α)表示頭部損傷響應的HIC值;fN(α)表示頸部損傷響應的 NIC 值;α =［α1，α2，α3，α4］T則是設計變量，分別為：頭－頸角度、體位角度、大腿角度和小腿角度。式(3)中的函數關系較為復雜，其精確表達式難以寫出，本文利用MADYMO軟件建立載人著陸沖擊的數值模型，并根據數值模擬試驗結果結合Sobol法對設計變量參數進行靈敏度分析，形成從載人空降人體損傷響應評價指標的提出、動力學模型的建立到其參數靈敏度分析的一套分析方法。其分析流程如圖2所示。

圖2 人體著陸姿勢分析流程圖Fig.2 The analysis procedure of human landing position

3 Sobol靈敏度分析方法

Sobol靈敏度分析方法是一種基于方差的蒙特卡羅法［11］。針對以上載人空降問題，首先定義一個 k(k=4)維的單元體Ωk作為變量的空間域，表示為：

基于Sobol法可以用總的方差表示所有參數對人體損傷的影響程度：

用偏方差表示人體各部位角度變化對人體損傷的影響程度：

用偏方差表示人體各部位角度相互之間的交叉作用對人體損傷的影響程度：

這里，Si稱為因素xi的一階靈敏度系數，表示xi對輸出的主要影響，即人體單個部位角度xi對人體損傷值f(α)的影響。Sij(i≠j)為二階靈敏度系數，估計兩參數的交叉影響(αi，αj所引起的變化不能由 αi，αj單獨影響之和直接表示)等等。并且：

人體各部位角度因素的總靈敏度系數由該因素各階靈敏度系數之和表示：

式中：S(i)指所有包含人體部位角度αi的靈敏度。

4 人體著陸姿勢變量靈敏度分析

4.1 “假人－座椅”數值模型的建立

某空降設備質量為8 000 kg，著地垂直速度為6 m/s，著陸角為 0，氣囊座椅內壓為1.4個大氣壓。整個“假人－座椅”數值模型包括氣囊緩沖座椅、假人、安全帶和戰車底板。其中，假人和戰車底板由MADYMO軟件中的多剛體構成，氣囊座椅和安全帶為有限元模型，如圖3所示。

根據氣囊座椅的三維結構模型，使用 Altair HyperMesh軟件作為前處理工具對座椅幾何模型進行最優網格劃分，并把節點和單元信息導入到MADYMO文件中，為MADYMO仿真運算準備其有限元模型。氣囊座椅采用尼龍66材料，其物理參數如表2所示。單元類型選擇三節點線性膜單元—MEM3單元。為了控制氣囊充氣后的形狀，氣囊內部通過拉帶確保氣囊座椅充氣后不會造成體積過度增大，達到控制氣囊座椅形狀的效果。整個座椅模型包括16506個單元。

圖3 “假人－座椅”數值模型Fig.3 The“dummy－seat”numerical model

表2 座椅材料尼龍66的物理參數Tab.2 The physical parameters of nylon 66

安全帶采用混合三點式安全帶模型，結合使用了一維多剛體安全帶模型和二維膜單元有限元安全帶模型。多剛體安全帶模型便于模擬安全帶在卷收器的拉入拉出和在滑環兩側的滑動，有利于準確的模擬安全帶與假人身體的接觸。

模型包括兩個輸入的載荷條件：重力加速度場和著陸沖擊加速度場。著陸沖擊過程中，戰車和座椅約束系統是連接在一起的。由于乘員的質量參數與戰車相比非常小，因此乘員運動對戰車響應的影響可以忽略。在計算機建模時可以將乘員的運動從戰車環境中隔離出來，將戰車環境視為慣性空間，而對乘員的外部作用定義為兩個加速度場：一個是作用于人體且方向垂直向下的重力加速度場，另一個是與實際碰撞加速度方向相反、大小相等的加速度場，該加速度即為戰車空降實驗中測得的車體減速度的反向曲線，如圖4所示。將兩加速度場進行疊加，即得到著陸沖擊過程中乘員相對戰車和座椅約束系統的加速度場。

圖4 沖擊加速度曲線Fig.4 The impact acceleration curves

4.2 數值模型的驗證

將同種工況下的某空降設備進行空降實驗，并通過人體頭部實驗數據來驗證“假人—座椅”數值模型的有效性。

從圖5可以看出，人體頭部加速度曲線的數值計算結果與實驗結果的變化趨勢吻合較好，兩者曲線峰值的大小和峰值出現的時間基本一致，能夠基本反應著陸沖擊時人體頭部的加速度響應，表明該數值模型是有效的，滿足要求可以用于人體著陸姿勢變量靈敏度分析。

圖5 假人頭部加速度實驗結果和數值結果對比Fig.5 The comparison of dummy head acceleration experimental results and numerical results

4.3 靈敏度分析

Sobol靈敏度方法一個顯著的特點是(4)式和(5)式的積分可由蒙特卡羅積分法［12］求得。因此根據設計變量的范圍，由拉丁超立方采樣技術采集50個不同的人體著陸姿勢，并通過MADYMO軟件分別建立與各人體姿勢相匹配的氣囊座椅模型，如圖6所示。這些模型中除了人體著陸姿勢不同之外，氣囊座椅的體積、初始內壓、充氣質量和溫度都相同，并施加相同的沖擊載荷。目的是盡可能的在同等條件下，考慮不同著陸姿勢對人體損傷響應值的影響。

圖6 不同人體姿勢下的“假人－座椅”數值模型Fig.6 Different“dummy－seat”numerical models

以上“假人－座椅”數值模型建立好之后，可對模型中假人的頭部損傷HIC值和頸部損傷NIC值進行計算，并采用Sobol法計算單個變量以及變量之間交互作用下對人體損傷響應值的影響。

如圖7所示，當考慮人體各部位角度單一變化對HIC值的影響程度時，頭－頸角度對人體頭部損傷響應變化的影響最大，靈敏度值為0.190 25。其次是體位角(0.154 56)、大腿角度(0.084 97)、和小腿角度(0.039 967)。如果考慮人體各部位變量之間的交互作用對HIC值的影響，則靈敏度值如圖8所示，影響最大的是體位角(0.749 65)，其次是頭－頸角度(0.614 87)、大腿角度(0.385 08)、小腿角度(0.135 32)。

由圖9可知，如果只考慮單個變量對人體損傷響應NIC值的影響，頭－頸角度對人體頸部損傷值變化的影響最大，靈敏度值為0.190 11。其次是體位角(0.164 55)、大腿角度(0.097 978)、小腿角度(0.071 361)。但是，考慮變量之間交互作用的影響的話，順序則發生了改變，如圖10所示。對人體頸部損傷值變化影響最大的是體位角(0.637 6)。隨后是頭－頸部角度(0.578 07)、大腿角度(0.381 77)、小腿角度(0.172 36)。

圖7 單個變量對HIC的影響Fig.7 The sensitivity of the single variable to HIC

圖8 變量交互作用對HIC的影響Fig.8 The sensitivity of variables interaction with others to HIC

圖9 單個變量對NIC的影響Fig.9 The sensitivity of the single variable to NIC

圖10 變量交互作用對NIC的影響Fig.10 The sensitivity of variables interaction with others to NIC

下面分別討論人體各個部位角度的變化對人體損傷響應的影響。

(1)頭－頸部角度

當僅考慮單個變量對人體損傷響應的影響時，頭－頸部角度不管是對HIC值，還是對NIC值的靈敏度都是最大的。這是因為當人體其它部位角度固定不變，頭－頸部角度發生變化時，頭部所受到的沖擊力方向也在發生改變。在沖擊載荷作用下，頭部的速度和加速度將直接影響到HIC和NIC值。

(2)體位角

當僅考慮單個變量對人體損傷響應的影響時，體位角排在第二。而考慮體位角和其它部位交互作用下對HIC和NIC值的影響時，靈敏度值則排在了第一。人體的體位角發生變化時，勢必會影響到人體頭部的受力方向。從頸部損傷式(2)可知，胸部T1的x方向加速度和速度對NIC值有著顯著作用，而體位角的改變勢必會對胸部的運動狀態產生重要影響。

(3)大腿角度和小腿角度

大腿角度和小腿角度無論是對HIC值，還是對NIC值的影響，相比頭－頸部角度和體位角，都要小得多。但當人體遭受到著陸沖擊載荷時，大腿和小腿會起到支撐的作用。所以它們對人體損傷響應的影響也有一部分貢獻。

5 結論

本文建立了載人著陸沖擊的“假人－座椅”數值模型，并進行了實驗驗證。利用Sobol法對人體不同部位角度對人體損傷響應的影響進行了分析。結果表明：人體頭—頸部角度和體位角對人體損傷響應的HIC值和NIC值起著舉足輕重的作用，而小腿角度不管是考慮它單獨對人體損傷的影響還是考慮它與其它參數的交互作用，都是所有靈敏度中最小的。因此，可根據變量靈敏度分析結果，為最佳人體著陸姿勢的研究提供變量選取依據，在載人空降防護技術領域具有一定的實際工程意義。

［1］譚 軍，韓 旭，劉 鑫.載人防護緩沖氣囊座椅的優化設計［J］.振動與沖擊，2011，30(2)：222－225.

［2］謝 燕，蔣建平，雷勇軍，等.飛船返回艙座椅著陸緩沖系統的局部改進模型［J］.振動與沖擊，2006，25(8)：785－787.

［3］劉炳坤，王憲明，王玉蘭，等.不同體位著陸沖擊時人體的動態響應［J］.航天醫學與醫學工程，2001，14(2)：120－122.

［4］ Cheng Z L，Han Y F，Zeng W Y.Human tolerance to landing impact(+Gz)in sitting position ［J］.Space Medicine＆ Medical Engineering，1997，10(5)：340－343.

［5］Brown W K，Rothstein J D，Foster L P.Human response to predicted apollo landing impacts in selected body orientations［J］.Aerospace Med，1966，37(4)：394－398.

［6］ Weis E B，Clarke N P，Brinkley J W.Human response to several impact orientations and patterns ［J］.Aerospace Med，1963，34(3)：1122－1129.

［7］Liu B K，Ma H L，Jiang S Z.Dynamic responses to landing impact at different key segments in selected body positions［J］.Aerospace Science and Technology，2008，12：331－336.

［8］杜匯良，黃世霖，張金換.高G值著陸沖擊下頭、頸、胸、脊椎損傷因子及評價指標的研究［J］.中國臨床康復，2004，8(8)：1474－1475.

［9］Bostrom O，Fredriksson R，Haland Y，et al.Comparison of car seats in low speed rear－end impacts using the BioRID dummy and the new neck injury criterion(NIC)［J］.Accident Analysis and Prevention，2000，32(2)：321－328.

［10］肖 志.汽車后碰撞中乘員頸部損傷防護的研究［D］.長沙：湖南大學，2007.

［11］ Sobol I.Sensitivity estimates for non－linear mathematical models［J］.Math Modeling Comput Exp，1993，4：407－414.

［12］ Nestor V Q，Raphael T H，Wei S，et al.Surrogate－based analysisand optimization ［J］. Progressin Aerospace Sciences，2005，41：1－28.

Characteristics of human head-neck injury response to manned airdrop based on sobol method

LIU Xin1，ZHANG Zhi－yong1，LIU Gui－ping2

(1.College of Automobile and Mechanical Engineering，Changsha University of Science and Technology，Changsha 410114，China;2.State Key Laboratory of Advanced Design and Manufacturing for Vehicle Body，Hunan University，Changsha 410082，China)

The “dummy－seat”numerical model was constructed and the results of numerical simulation are approaching to the experimental data.The sensitivity degrees of the angles of various body parts to the human injury were investigated based on the Sobol method.The results demonstrate that the head－neck angle and body angle have great influence on human injury responses.According to the analysis results，the engineers are able to datermine the optimal astronaut landing attitude.The method can also be used in the field of other airdrop protections.

landing impact;human injury response;sobol method;sensitivity analysis

V224

A

2011－06－09 修改稿收到日期：2011－07－22

劉 鑫 男，博士，講師，1981年生