曹立波，張 愷，顏凌波，湯 駿

(湖南大學，汽車車身先進設計制造國家重點實驗室，長沙 410082)

2016132

基于中國人體測量學尺寸的假人頭部跌落試驗的仿真研究*

曹立波，張 愷，顏凌波，湯 駿

(湖南大學，汽車車身先進設計制造國家重點實驗室，長沙 410082)

為探究基于中國人體測量尺寸的假人頭部在跌落試驗中的響應與混合III型假人頭部響應的差異，采用兩種常用的假人縮放方法建立了中國5百分位、50百分位和95百分位假人頭部有限元模型，并在3種不同沖擊速度下進行跌落試驗仿真。結果表明，非均一縮放方法獲得的中國假人頭部模型幾何尺寸與中國人體頭部測量學尺寸更為吻合，而采用該方法獲得的中國人3種體型的假人頭部合成加速度峰值均比對應的混合III型假人頭部模型高3%～6%。

假人頭部；中國人體測量學尺寸；跌落試驗仿真；縮放方法

前言

據統計，在美國每年約有30%的致命傷是由創(chuàng)傷性顱腦損傷(TBI)所引起的[1]，盡管研究人員對頭部損傷機理的研究不斷深入，在此基礎上采取了許多措施來降低交通事故中的頭部損傷程度，例如戴頭盔和乘員約束系統等，頭部損傷的數量和嚴重程度也得到了明顯降低，然而，交通事故中的頭部損傷仍然是導致人員死亡和傷殘的一個最主要原因[2]。汽車碰撞假人是根據人體生物力學研究成果而設計的模擬人體結構與力學特性的機械裝置，是評估乘員在碰撞載荷下遭受傷害風險的重要工具[3]。目前，在世界范圍內，混合III型系列假人是應用最為廣泛的汽車正面碰撞測試假人，其中，50百分位男性假人表征美國中等體型成年男性，5百分位女性假人和95百分位男性假人通過縮放50百分位男性假人獲得，分別表征小體型成年女性和較大體型成年男性。研究表明，混合III型50百分位假人頭部輪廓與美國中等體型人體頭部輪廓十分吻合[4]，同時，在剛性碰撞中可以準確地模擬人體頭部在相同工況下的力學響應[5-6]。

由于混合III型系列假人的人體統計學參數參考的是美國20世紀80年代的研究成果，該參數與中國人體參數差異較大，2011年的統計數據顯示，中國50百分位男性身高與混合III型假人身高的差異達到了3.48%，體質量的差異更是達到了11.89%[7]。文獻[8]和文獻[9]中對中國假人模型與混合III型50百分位假人模型(標準假人模型)正面碰撞響應的差異進行了研究。研究表明，中國假人模型與標準假人模型在車輛正面碰撞中的損傷風險差異很大，開發(fā)具有中國人體特征的假人具有重要意義。但是，上述研究是利用MADYMO軟件縮放得到的符合中國50百分位人體身材的假人模型，并未對縮放方法進行探究。同時，上述研究討論的是兩種假人在約束系統中的響應差異，由于在正面碰撞中假人的動力學響應與乘員約束系統的設計參數直接相關，這樣不利于從根本上研究兩種假人的損傷差異。此外，之前的研究大多集中于假人的整體響應，并未對某一特定部位的損傷差異進行研究。

由此可見，深入研究基于中國人體尺寸的假人頭部在跌落試驗中的響應對研究中國人體頭部損傷及開發(fā)中國汽車碰撞測試假人具有重要意義。本文中首先探究了兩種縮放方法獲得中國人體假人頭部有限元模型的準確性，建立了5百分位、50百分位和95百分位3種假人頭部的跌落試驗仿真模型，在此基礎上研究了中國假人頭部模型與混合III型假人頭部模型在跌落試驗中的加速度差異。

1 假人頭部模型縮放方法

研究人員常常采用縮放的方法通過現有成熟的假人模型獲得其他百分位假人模型，這樣既降低了假人模型的開發(fā)周期與開發(fā)成本，也在一定程度上保證了假人的生物逼真度。在處理成年人的縮放系數時，常常假設成年人的生物組織材料密度、剛度是一致的[10]，所以使用縮放方法獲得中國成年假人模型的過程中，幾何縮放因數是一個重要的參數。

混合III型50百分位假人頭部輪廓是通過大量樣本均一化處理后獲得的，為獲取其他體型人體的假人頭部輪廓，可以根據目標人體頭部的幾何參數通過線性縮放獲得[11]。通過文獻研究，常用的假人頭部模型的縮放方法有兩種，一種是均一縮放方法，在開發(fā)混合III型5百分位女性假人頭部和混合III型95百分位男性假人頭部時便采用了這種方法[12]，即在假人頭部局部坐標系下X向、Y向和Z向的縮放系數λx，λy和λz相等，其中頭部局部坐標系如圖1所示，頭部縮放系數的具體確認方法為

(1)

λm=λx·λy·λz

(2)

式中：C為頭部周長；W為頭部寬度(Y方向)；L為頭部長度(X方向)；m為頭部質量；λm為頭部質量縮放系數；下標1表示中國人體參數，下標0表示美國人體參數。

另一種縮放方法為非均一縮放方法，在開發(fā)50百分位女性后碰撞假人—EvaRID時，瑞典查爾莫斯科技大學的研究人員便采用了該方法[13]，頭部縮放系數的具體確認方法為

λx=L1/L0

(3)

λy=W1/W0

(4)

λz=λm/(λy·λx)

(5)

2 假人頭部模型的建立

為更加深入研究基于中國人體尺寸的假人頭部模型與混合III型假人頭部模型的差異，本文中選取了5百分位、50百分位和95百分位3種體段人體的參數進行研究，如表1所示，其中中國人體參數來自GB/T 10000—1988，開發(fā)假人的參考人體參數來自AMVO(the anthropometry for motor vehicle occupants)中的統計數據[11]。

表1 頭部尺寸參數

根據上述的兩種縮放方法，得到了3種體型中國假人頭部的縮放系數，如表2所示，其中系數大于1表示在原模型的基礎上進行放大，系數小于1表示在原模型的基礎上縮小。

表2 頭部縮放系數

用于開發(fā)中國5百分位、50百分位和95百分位假人頭部有限元模型的原模型為LSTC公司與美國國家碰撞分析中心(NCAC)共同開發(fā)的混合III型5百分位、50百分位和95百分位假人有限元模型。

以50百分位假人頭部模型為例，圖2中從左至右依次為混合III型50百分位假人頭部有限元模型、均一縮放中國50百分位假人頭部有限元模型和非均一縮放中國50百分位假人頭部有限元模型。為深入比對通過兩種縮放方法得到的中國假人頭部的幾何差異，分別抽取了3種50百分位假人頭部模型矢狀面和冠狀面的輪廓線，如圖3所示。結合表1和圖3可以發(fā)現，中國50百分位人體頭部寬度尺寸基本一致，長度和周長尺寸差異較大，若采用均一系數的縮放方法，得到的中國50百分位假人頭部模型在矢狀面和冠狀面的輪廓線與中國人體的頭部幾何參數均有明顯偏差。由于絕大部分的創(chuàng)傷性顱腦損傷都發(fā)生于頭部與其他物體的碰撞過程中，所以假人頭部輪廓與人體頭部輪廓高度一致是保證假人生物逼真度的重要前提，相對而言，非均一縮放方法更適合用于開發(fā)中國尺寸假人模型。

3 假人頭部跌落試驗仿真

雖然LSTC與NCAC共同開發(fā)的混合III型假人有限元模型得到了充分的驗證[14]，在進行中國假人頭部與混合III型假人頭部跌落試驗差異研究之前，本文中對混合III型50百分位假人頭部有限元模型在高速、中速和低速3種試驗工況下的質心加速度與試驗數據進行對標。第一種試驗設置為：調整假人的頭部，使假人額頭最低點處于假人鼻子最低點的下方12.7mm處，將假人頭部從376mm的高處自由跌落到剛性平面上，即假人頭部以2.72m/s的初速度撞擊剛性平面，該試驗也就是假人頭部標定試驗，是檢驗假人生物逼真度與驗證假人有限元模型必不可少的試驗。后兩種試驗設置與第一種試驗類似，只是將撞擊初速度調整為1.9和4m/s，試驗過程中采集頭部加速度曲線，并使用SAE 1000等級進行濾波。參考試驗設置建立假人頭部跌落試驗仿真模型，如圖4所示，調用LS-DYNA有限元求解器，仿真后輸出頭部質心合成加速度曲線，同樣使用SAE 1000等級進行濾波。

仿真與試驗得到的混合III型50百分位假人頭部在3種跌落試驗中的頭部質心合成加速度曲線如圖5所示，其中撞擊初速度為4m/s的試驗曲線來自文獻[15]中的研究,加速度峰值如表3所示。從圖中可以看出，3種試驗工況下有限元仿真模型與試驗曲線吻合較好，但是仿真曲線的波形較寬，合成加速度峰值的差異均小于3%。

撞擊初速度v0/(m·s-1)合成加速度峰值/g試驗仿真差異/%1.9159.43157.231.382.72260.12258.70.554437448.432.62

對利用非均一縮放方法獲得的中國人體5百分位、50百分位和95百分位假人頭部建立3種碰撞速度的跌落試驗仿真模型，仿真結果如圖6所示。在3種試驗工況下，基于中國人體尺寸建立的3種體型的假人頭部的合成加速度峰值均高于對應的混合III型假人頭部模型的合成加速度峰值，差異在3%～6%，也就是說在相同碰撞強度下，中國假人頭部會遭受更加嚴重的損傷。與此同時，隨著體型的增大，假人的頭部合成加速度峰值呈現降低趨勢。

造成中國假人頭部合成加速度峰值高于混合III型假人頭部模型的原因是多方面的。首先，中國人體頭部與用于混合III型假人開發(fā)的美國人體頭部在幾何尺寸上存在差異。其次就是頭部質量方面的差異，在分析假人頭部跌落試驗中的頭部質心加速度時，可將假人頭部簡化為一個單一的彈簧-質量模型[12]，根據能量守恒定律得

(6)

式中：v為碰撞速度；d為質量塊最大移動量；k為彈簧剛度。根據牛頓第二定律和胡克定律可得

F=ma

(7)

F=kd

(8)

式中：F為碰撞力；a為碰撞中質量塊加速度。

綜合式(6)～式(8)，可得

(9)

由于試驗條件相同，即兩組試驗的碰撞速度v相同，同時可認為彈簧-質量模型的彈簧剛度近似，則碰撞過程中加速度峰值的平方反比于質量，即頭部質量越小，加速度峰值越大，通過改變混合III型假人頭部質量，在標定試驗中頭部質心加速度峰值的變化規(guī)律也印證了此點[16]。再次，縮放導致假人頭部皮膚厚度的變化也是引起頭部質心加速度峰值差異的重要原因。通過仿真分析，將中國假人頭部皮膚的厚度調整為與混合III型假人頭部皮膚的厚度一致時，二者在跌落試驗中加速度峰值的差異僅為1%～3%。當然，在縮放過程中頭部形狀、頭部質量和頭皮厚度這些因素都會隨之變化，三者共同作用導致了中國假人頭部質心加速度峰值高于混合III型假人頭部加速度峰值的結果。

4 結論

歸納總結兩種開發(fā)假人常用的幾何縮放方法，結合中國人體與混合III型假人開發(fā)時期美國人體在頭部幾何參數的差異，通過縮放方法獲得了一系列符合中國人體尺寸的假人頭部有限元模型，通過3種工況下頭部跌落試驗的仿真研究，得到如下結論。

(1) 通過非均一縮放方法獲得中國假人頭部模型是一種更為理想的開發(fā)策略。

(2) 在3種試驗工況下，基于中國人體尺寸建立的3種體型的假人頭部的合成加速度峰值均比對應的混合III型假人頭部模型的仿真結果高3%～6%。

(3) 隨著體型的增大，假人的頭部合成加速度峰值呈下降趨勢。

目前，我國現行的汽車安全法規(guī)與新車評價規(guī)程中汽車正面碰撞試驗均使用混合III型假人，由于我國人體在體質量、身高與假人具有明顯差異，在相同沖擊環(huán)境下，中國假人頭部會承受更大的加速度，使用混合III型假人替代中國人體是否合理仍存在不小爭議。因此，建議我國在制定法規(guī)以及開發(fā)新車型時考慮中國人體與現行假人在幾何尺寸與損傷風險的差異，為我國乘員提供更好的保護，同時，對于符合中國人體尺寸假人的研究工作必須穩(wěn)步推進，在沒有足夠充分的人體生物力學試驗數據的前提下，采用非均一縮放方法獲得中國假人頭部模型具有參考價值。

[1] FAUL M, XU L, WALD M M, et al. Traumatic brain injury in the United States[C]. Atlanta, GA: Centers for Disease Control and Prevention, National Center for Injury Prevention and Control,2010.

[2] 汽車與運動損傷生物力學[M].北京:機械工業(yè)出版社,1992.

[3] FOSTER J, KORTGE J, WOLANIN M. Hybrid III-a biomechanically-based crash test dummy[C]. SAE Paper 770938.

[4] LEE C, LOYD A M, NIGHTINGALE R, et al. Three-dimensional adult male head and skull contours[J]. Traffic Injury Prevention,2014,15(4):402-406.

[5] HUBBARD R P, MCLEOD D G. Definition and development of a crash dummy head[C]. Proceedings: Stapp Car Crash Conference. Society of Automotive Engineers SAE,1974,18.

[6] LOYD A M, NIGHTINGALE R W, SONG Y, et al. The response of the adult and ATD heads to impacts onto a rigid surface[J]. Accident Analysis & Prevention,2014,72:219-229.

[7] CAO L, ZHANG K, Lü X, et al. Comparison of current ATDs with Chinese adults in anthropometry[J]. Traffic Injury Prevention,2016,17(4):430-433.

[8] 曹立波,余志剛,白中浩,等.中國假人模型與 HybridⅢ 假人模型碰撞響應的差異[J].中國機械工程,2009,10(3):370-374.

[9] 白中浩,曹立波,余志剛.中國 50 百分位人體與標準試驗假人正面碰撞響應差異的研究[J].汽車工程,2008,30(11):993-997.

[10] SUN E, EPPINGER R, KUPPA S, et al. Development of improved injury criteria for the assessment of advanced automotive restraint systems[M]. Washington, DC: National Highway Traffic Safety Administration,1998.

[11] LOYD A M, DANIEL C, MUKUNDAN S, et al. 3-D average skull contour data for 3-year-olds and 5-year-olds: a pilot study[C]. Injury Biomechanics Research Thirty-Fourth International Workshop,2006.

[12] MERTZ H J, IRWIN A L, MELVIN J W, et al. Size, weight and biomechanical impact response requirements for adult size small female and large male dummies[C]. SAE Paper 890756.

[13] CARLSSON A, CHANG F, LEMMEN P, et al. EvaRID-A 50th percentile female rear impact finite element dummy model[C]. Proceedings: IRCOBI Conference,2012:12-14.

[14] Livermore Software Technology Corporation. LSTC_NCAC hybrid III 50th dummy positioning & post-processing[R]. Michigan USA,2013.

[15] CHOU C C, ZHAO Y, HUANG Y, et al. Development and validation of a deformable featureless headform model using LS-DYNA3D[C]. Proceedings: International Technical Conference on the Enhanced Safety of Vehicles. National Highway Traffic Safety Administration,1996:1313-1328.

[16] 曹立波,張愷,呂鑫,等.混合Ⅲ型第50百分位假人頭部跌落試驗影響因素探究[J].汽車工程,2016,38(3):302-307.

A Simulation Study on the Drop Test of Dummy HeadBased on Chinese Anthropometric Dimensions

Cao Libo, Zhang Kai, Yan Lingbo & Tang Jun

HunanUniversity,StateKeyLaboratoryofAdvancedDesignandManufacturingforVehicleBody,Changsha410082

In order to investigate the difference of head responses between the dummy based on Chinese anthropometric dimensions and Hybrid-III dummy in drop simulation, two common dummy scaling methods are adopted to create the finite element models for the 5th, 50th and 95th percentile Chinese dummy head, on which a drop-test simulation is conducted at three different impact speeds. The results show that the geometric dimensions of Chinese dummy head model obtained by inhomogeneous scaling method are more consistent with Chinese anthropometric dimensions, and the resultant acceleration peak values of three percentile Chinese dummy heads obtained by that method are 3%～6% higher than corresponding Hybrid III dummy head model.

dummy head; Chinese anthropometric dimensions; drop-test simulation; scaling method

*湖南大學汽車車身先進設計制造國家重點實驗室開放基金(31575005)和國家自然科學基金(11172099)資助。

原稿收到日期為2015年5月28日。