肖志　張云飛　龐通　李順峰　劉國斌　王方

摘? ?要：為了提高摩托車頭盔防護性能降低頭部損傷風險，本文引入功能梯度仿生泡沫結構替代傳統(tǒng)頭盔的均勻密度泡沫襯墊. 通過頭盔耦合生物力學頭部有限元模型，獲取撞擊過程中頭部質心加速度、生物力學響應和頭盔結構響應，綜合分析不同密度梯度方案對頭盔防護能力的影響. 結果表明，功能梯度泡沫結構在中高速沖擊下有更明顯的優(yōu)勢. 相比傳統(tǒng)頭盔吸能緩沖層的均勻泡沫密度方案和新型頭盔的正/負密度梯度方案，最大密度為80 kg/m3的負密度梯度泡沫設計方案可以更有效改善頭盔結構碰撞響應并降低碰撞過程中頭部損傷，且隨著密度差的增大，負密度梯度頭盔的防護性能得到了進一步提高.

關鍵詞：功能梯度泡沫;仿生頭盔;頭部損傷;有限元模型;碰撞

中圖分類號：U467.14? ? ? ? ? ? ? ? ? ? ? ? ? ? ?文獻標志碼：A

Analysis on Head Injury and Protective Performance

of Functionally Graded Bionic Helmet

XIAO Zhi ZHANG Yunfei PANG Tong LI Shunfeng LIU Guobin WANG Fang

（1. State Key Laboratory of Advanced Design and Manufacturing for Vehicle Body，Hunan University，Changsha 410082，China;

2. Dep of Mechanical，Aerospace and Civil Engineering，University of Manchester，Manchester M13 9PL，UK;

3. The Sixth Research Institute of China Aerospace Science and Industry Corporation，Hohhot，010010，China;

4. College of Automotive and Mechanical Engineering，Changsha University of Science & Technology，Changsha 410082，China）

Abstract：Functionally graded foam is introduced in this paper to replace the uniform-density energy-absorbing liner of a conventional motorcycle helmet for improving its protection performance and further reducing head injuries during impact scenarios. The helmet finite element model coupled with biomechanical head model was developed，the acceleration transmitted to the gravity center of head，biomechanical responses of the head and crushing behaviors of the helmet under the impacts are obtained to comprehensively investigate the effects of functionally graded foam on helmet impact protection. According to the results，functionally graded foam has more advantages under medium/high speed impact than uniform liner foam. In comparison with the uniform-density design for the conventional helmet liner and the positive/negative functionally graded foam design for the novel helmet liner，the negative functionally graded foam with a maximum density of 80 kg/m3 is of the best crushing responses and the severities of head injuries can be reduced more effectively. With the increase of density difference，the impact protection of novel helmet with negative functionally graded foam design can be further improved.

Key words：functionally graded foam;bionic helmet;head injury;finite element model;impact

2020年，“一盔一帶”安全守護行動的開展對摩托車頭盔防護性能提出了更高的要求. 在摩托車碰撞事故中，頭部損傷是駕乘人員重傷甚至死亡的主要原因，而佩戴頭盔可以有效降低頭部損傷程度. 在實際碰撞中，頭盔主要作用是耗散沖擊過程中的碰撞能量，吸收碰撞過程中頭部動能，緩沖碰撞過程.

吸能緩沖層作為摩托車頭盔最主要的吸能結構，通過碰撞過程中泡沫的塑性變形吸收大部分碰撞能量. 近年來，為進一步提高吸能緩沖層的能量耗散水平，一批諸如蜂窩結構[1]、波紋形多孔結構[2]和桁架結構[3]等多種吸能結構被引入頭盔結構設計中，并在一定程度上提高了頭盔的防護能力.

功能梯度泡沫作為一種輕質高效的吸能結構，首先在生物材料中被觀察到，隨后引入到耐撞性結構研究和馬術頭盔中. Fischer等人[4]和Fleck等人[5]觀察了柚子皮的微觀結構，發(fā)現(xiàn)了明顯的梯度分布現(xiàn)象，并通過柚子自由落體實驗驗證了柚子皮良好的抗沖擊能力. 隨后，Cui等人[6]探究了EPS（Expanded Polystyrene）泡沫梯度布置方案在低速沖擊下的吸能特性，在此基礎上，將其引入到馬術頭盔，通過分析頭盔碰撞響應和頭部加速度響應，探究了功能梯度結構對馬術頭盔防護能力的影響[7]. 但截至目前，功能梯度泡沫在摩托車頭盔中的應用及頭部生物力學響應并沒有得到很好的研究.

為進一步提高頭盔的防護吸能，并填補功能梯度頭盔碰撞響應與頭部損傷生物力學響應耦合分析的研究空白，本研究首先建立了摩托車頭盔有限元模型，并通過標準試驗驗證頭盔模型;基于泡沫材料相關理論，獲得一系列不同密度EPS泡沫本構模型，并將其應用到功能梯度頭盔模型中;最后，通過頭盔耦合生物力學頭部模型，探究了功能梯度結構在不同沖擊速度下的防護性能，在此基礎上研究了7.5 m/s沖擊速度下不同梯度方案的頭部加速度響應、生物力學響應和頭盔結構碰撞響應，綜合評價功能梯度泡沫結構對頭盔防護能力的影響.

1? ?有限元模型

本節(jié)介紹了所用的兩種頭部有限元模型，建立并驗證了摩托車頭盔有限元模型，利用相關理論獲取了所需的EPS泡沫材料本構模型.

1.1? ?頭部模型

本文采用由LSTC公司（Livermore Software Technology Corporation）開發(fā)的Hybrid Ⅲ剛性頭部模型輸出頭部質心的加速度曲線來和實驗對標[8]，同時采用由日本豐田汽車研究院開發(fā)的THUMS（Total Human Model for Safety）生物力學頭部模型輸出頭部生物力學響應和質心加速度分析頭部損傷[9]. 其中THUMS頭部模型包括人體腦部、顱骨、腦膜、腦脊液及面部等組織，能夠很好地模擬碰撞工況下的頭部生物力學響應.

1.2? ?頭盔模型

摩托車頭盔主要由硬質外殼、吸能緩沖層、舒適襯墊、系帶和護目鏡組成，其中，頭盔的防護吸能特性主要是由硬質外殼和吸能緩沖層決定，絕大多數(shù)摩托車頭盔外殼和吸能緩沖層分別由ABS（Acrylonitrile Butadiene Styrene）塑料和EPS泡沫制成.

如圖1所示，選取AGV-K3摩托車頭盔通過GOM三維掃描儀掃描得到相應的Stl（Stereolithography）模型，再經(jīng)過Geomagic軟件逆向重構得到頭盔幾何模型，隨后在Hyperworks軟件中完成模型的前后處理，最后選用LS-DYNA971 R4.2求解器進行模型計算.

本文中僅對頭盔主要部件（即：ABS外殼、EPS泡沫吸能緩沖層、系帶）進行建模，忽略對緩沖吸能影響很小的其余部件. 頭盔外殼厚度為4 mm，選用LS-DYNA中*MAT_3塑性隨動硬化材料模型和3.5 mm全積分四邊形殼單元網(wǎng)格進行仿真. 吸能緩沖層厚度范圍為6 mm - 45 mm，選用*MAT_63壓縮泡沫材料模型進行模擬，本構曲線在下節(jié)1.3中推導得出. 在撞擊區(qū)域選擇4 mm全積分六面體網(wǎng)格平均劃分泡沫厚度，每層網(wǎng)格定義為一個部件，定義從最外層到最內(nèi)層分別為第1層至第8層，各層之間采用共節(jié)點連接，對每一層賦予不同的材料本構參數(shù)，層內(nèi)密度均勻，但在厚度方向上實現(xiàn)密度梯度變化，襯墊的其他區(qū)域選擇單點積分四面體網(wǎng)格. 系帶總長度為200 mm，選用*MAT_1彈性材料模型和4 mm四邊形網(wǎng)格進行模擬. 外殼與吸能緩沖層之間采用Tied_node_to_surface接觸;吸能緩沖層與頭部模型之間，外殼與剛性平砧之間均采用Automatic_surface_to_surface接觸. 頭盔模型中各部分具體的材料參數(shù)如表1所示.所建立的頭盔有限元模型如圖2所示.

1.3? ?EPS泡沫材料應力-應變曲線

本文利用Schraad和Harlow[11]提出的隨機多孔材料理論推導出一系列泡沫壓縮應力-應變曲線，并將這些曲線應用到頭盔模型吸能緩沖層中.

1.4? ?模型驗證

根據(jù)摩托車頭盔法規(guī)GB 24429-2009[13]和ECE R22.05[14]，選取頭盔B點位置進行跌落碰撞實驗，B點所在位置如圖5所示. 該實驗在頭盔性能測試試驗機中完成，如圖6所示. 在實驗中頭盔和M號剛性金屬頭型固定在一起，將頭盔拉升至距離剛性平砧2.9 m高度處釋放，使其自由下落產(chǎn)生7.5 m/s的撞擊速度，實驗重復三次，實驗數(shù)據(jù)通過內(nèi)置于剛性頭型質心位置的加速度傳感器測得. 如圖7所示，頭盔跌落實驗結果和對應仿真（圖2（a））的加速度-時間曲線顯示了良好的一致性，頭盔模型的準確性得到驗證.

2? ?功能梯度頭盔梯度參數(shù)選擇

如圖8所示，基于柚子皮組織梯度分布現(xiàn)象，本文設計了13種密度梯度方案探究仿生梯度EPS泡沫結構對頭盔防護能力的影響. 各種梯度方案如表2所示，其命名方法如下：Ave80和Max80分別表示泡沫襯墊平均密度和最大密度為80 kg/m³;Δρ20和Δρ-20分別表示襯墊最外層（第1層）與最內(nèi)層（第8層）泡沫密度之差為20 kg/m³和-20 kg/m³. 例如，Ave80_Δρ20代表平均密度為80 kg/m³且最外層與最內(nèi)層泡沫密度差為20 kg/m³的正密度梯度方案.

3? ?功能梯度頭盔防護性能分析

3.1? ?均勻密度頭盔防護性能分析

基于驗證的頭盔耦合生物力學頭部模型（圖2（b））對均勻密度頭盔進行跌落仿真. 頭盔EPS泡沫吸能緩沖層的碰撞變形過程如圖9所示;頭部加速度曲線及對應的腦壓云圖和最大主應變云圖如圖10所示. 在頭盔-頭部模型和剛性砧之間開始接觸碰撞時，傳遞至頭部的加速度迅速升高;當碰撞時間t=2 ms時，加速度增長速率開始降低，各層泡沫撞擊點處應變?yōu)?.12-0.29，此時大部分泡沫處于平臺吸能階段的前中期（圖3和圖4），泡沫吸能效率開始增大;在碰撞時刻t=6 ms，加速度達到峰值，頭部模型運動到極限位置;隨后頭型開始回彈，加速度快速降低. 當加速度到達峰值時，大腦壓力云圖呈現(xiàn)出明顯的梯度分布，最大腦正壓出現(xiàn)在撞擊位置附近（A點），最大主應變在額葉和顳葉附近（C點）取得最大值，這與Zhou等人[15]的研究結果具有較好的一致性.

3.2? ?沖擊速度對頭盔防護性能的影響

本部分在標準沖擊速度的基礎上，進一步探究了5 m/s、7.5 m/s和9 m/s三種沖擊速度下13種功能梯度方案（見表2）頭盔結構響應和頭部損傷響應;其中，Ave80_Δρ0、Max80_Δρ50和 Max80_Δρ-50三種典型梯度方案的大腦壓力云圖和頭盔吸能緩沖層等效應變云圖分別如圖11和圖12所示.

從壓力分布看，當沖擊速度增大時，最大腦正壓的峰值和高壓力區(qū)分布面積均逐漸增大. 和均勻密度方案相比，低速沖擊兩種典型功能梯度方案的最大腦正壓分布面積并沒有明顯變化，中高速沖擊兩種功能梯度方案的高壓力區(qū)分布面積顯著降低，且負密度梯度方案在降低頭部損傷方面的優(yōu)勢更加明顯.

從吸能緩沖層應變云圖看，隨著沖擊速度的增大，三種梯度方案峰值應變逐漸增大，高應變分布區(qū)域也逐漸擴大. 低速碰撞時，均勻密度方案的內(nèi)外層泡沫變形均不充分，正梯度方案可以進一步增大內(nèi)層部分泡沫的變形量，同時基本保持外層泡沫變形量，但由于應變顯著增大區(qū)域體積較小，優(yōu)勢并不明顯;而負梯度方案不僅沒有顯著改善外層泡沫的變形吸能，反而抑制了內(nèi)層泡沫的變形. 中高速沖擊下均勻密度方案的內(nèi)層泡沫變形已比較充分，此時的正梯度方案引起的內(nèi)層泡沫進一步變形增大了峰值應變和高應變區(qū)域體積，但抑制了外層泡沫的變形吸能;負密度梯度方案在保證內(nèi)層應變足夠大（平臺階段后期）的基礎上極大地改善了外層泡沫的變形吸能. 考慮到低速沖擊下頭部損傷較輕，損傷概率較低，功能梯度結構在中高速沖擊下的應用潛力更大.

3.3? ?標準沖擊速度頭部損傷結果分析

頭部損傷標準HIC（Head Injury Criterion）作為當前應用最廣泛的頭部損傷評價指標，可以在一定程度上反映頭部的整體損傷，如公式（4）所示[16]. 然而，在實際碰撞中，常見的損傷有腦挫傷、腦血腫、彌漫性軸索損傷等，HIC并不足以對各種具體損傷做精準預測. 相關研究表明，腦部最大腦正壓和最大腦負壓可以較好地預測腦挫傷和腦血腫，同時結合大腦最大主應變可以對彌漫性軸索損傷做較為精準的預測[17]. 其中，最大腦正壓和最大腦負壓為撞擊過程中腦組織受慣性影響產(chǎn)生顱骨和大腦相互碰撞，并在撞擊側和撞擊對側分別產(chǎn)生的峰值正負腦壓，是常用的腦損傷評價指標[18].

3.4? ?頭盔-頭部碰撞力學響應分析

3.4.1? ?頭部接觸力響應

3.4.2? ?吸能緩沖層變形分析

圖15對比了各層EPS泡沫撞擊點位置的最大變形量. 對于均勻密度泡沫，第1層（最外層）泡沫的變形量最小，而第8層（最內(nèi)層）泡沫的變形量最大，內(nèi)層泡沫吸能潛能發(fā)揮得較為充分. 而六種正梯度方案內(nèi)層泡沫變形進一步減小，不僅不利于能量吸收，反而會使頭部接觸力超出相同密度差的負密度方案，甚至均勻密度方案，不利于頭部損傷降低. Ave80的負梯度方案雖然增大了外層泡沫的變形，但總變形沒有明顯改變，而Max80的負梯度方案在保持各層泡沫變形比例均衡的前提下提升了總變形量，泡沫利用率進一步提升. 尤其是梯度方案Max80_Δρ_50，各層泡沫變形量和總變形量提升程度達到最大，更有利于頭盔碰撞吸能能力的提升.

3.4.3? ?頭盔吸能分析

圖16對比了功能梯度頭盔各部分能量吸收情況. 6種正梯度方案頭盔總吸能略低于均勻密度方案，且密度差越大，總吸能越小，主要原因是外層泡沫變形量的降低使得外殼吸能出現(xiàn)較大程度的降低. 6種負梯度方案總吸能均高于均勻密度方案和6種正梯度方案，且隨著密度差的增大，頭盔吸能逐漸增大，其主要原因是外層泡沫變形量增大使得外殼吸能出現(xiàn)較大幅度的增大，此時外殼吸能能力成為頭盔防護能力的主要影響因素. 此外，雖然Max80和Ave80梯度方案頭盔吸能相差不大，但是由于Max80方案的質量降低引起的初始動能降低，頭部剩余動能遠小于Ave80梯度方案，使得Max80梯度方案頭部損傷遠低于Ave80梯度方案.

3.5? ?參數(shù)相關性分析

為探究頭盔結構響應與頭部損傷參數(shù)（HIC，最大腦正壓，最大腦負壓，最大主應變）之間的對應關系，生物力學損傷參數(shù)分別對頭部接觸力、各層泡沫吸能標準差、各層泡沫變形標準差進行了相關性分析，如圖17和圖18所示. 由圖17可知，四種生物力學損傷參數(shù)相對減小量均隨著接觸力的增大而逐漸降低，對四種參數(shù)進行線性擬合進一步推導得到損傷參數(shù)與峰值接觸力之間的線性表達式如圖所示，其決定系數(shù)R2分別為0.957、0.961、0.925、0.993，線性關系明顯.

在圖18中，相對標準差為一種梯度方案下8層泡沫變形（或吸能）的標準差相對均勻密度方案標準差的歸一化處理，來表征各層泡沫變形（或吸能）的差異性. 由圖18（a）知，六種正梯度方案吸能標準差遠大于負梯度方案，且六種負梯度方案的吸能標準差隨著密度差的增大而逐漸減小，而正梯度方案則呈現(xiàn)相反的趨勢，這說明大密度差負梯度方案在改善頭盔各部分吸能比例方面優(yōu)勢更明顯. 由圖18（b）知，負梯度方案的變形標準差小于正梯度方案，但負梯度方案對密度差的相關性不明顯，這與圖18（a）的結果存在差異，可能原因是在撞擊過程中泡沫存在較大的剪切變形，吸收了另一部分碰撞能量.

4? ?結? ?論

本文通過驗證的頭盔-生物力學頭部耦合模型對比研究了均勻密度泡沫的傳統(tǒng)頭盔與功能梯度泡沫結構的新型頭盔的防護性能. 碰撞過程中，對不同頭盔保護下的頭部損傷進行了研究，同時還分析了頭盔結構碰撞響應，得出的主要結論如下：

1）功能梯度結構在中高速沖擊下對頭盔防護性能的改善效果明顯優(yōu)于低速沖擊;

2）中高速沖擊下，對比傳統(tǒng)頭盔吸能緩沖層的均勻泡沫密度方案和新型頭盔的正/負密度梯度方案，具有負密度梯度泡沫設計方案具有最佳的結構碰撞響應并能有效降低碰撞過程中頭部損傷，且隨著密度差的增大，負密度梯度頭盔的防護性能得到了進一步提高;

3）在13種密度梯度方案中，最大密度為80 kg/m³的方案（Max80）在降低頭部損傷方面優(yōu)于平均密度為80 kg/m³的梯度方案（Ave80）;

4）在頭盔結構響應與頭部生物損傷參數(shù)相關性分析中，頭部接觸力與生物力學損傷參數(shù)之間具有很好的線性關系，且較大密度差的負梯度泡沫可以進一步提高吸能緩沖層的整體變形吸能程度，從而提升頭盔防護性能并降低頭部損傷的發(fā)生.

參考文獻

[1]? ? CASERTA G D，IANNUCCI L，GALVANETTO U.Shock absorption performance of a motorbike helmet with honeycomb reinforced liner[J].Composite Structures，2011，93（11）：2748—2759.

[2]? ? BLIVEN E，ROUHIER A，TSAI S，et al. Evaluation of a novel bicycle helmet concept in oblique impact testing[J]. Accident Analysis & Prevention，2019，124：58—65.

[3]? ? KHOSROSHAHI S F，TSAMPAS S A，GALVANETTO U. Feasibility study on the use of a hierarchical lattice architecture for helmet liners[J].Materials Today Communications，2018，14：312—323.

[4]? ? FISCHER S F，THIELEN M，LOPRANG R R，et al. Pummelos as concept generators for biomimetically inspired low weight structures with excellent damping properties[J]. Advanced Engineering Materials，2010，12（12）：B658—B663.

[5]? ? B？譈HRIG-POLACZEK A，F(xiàn)LECK C，SPECK T，et al. Biomimetic cellular metals-using hierarchical structuring for energy absorption[J].Bioinspiration & Biomimetics，2016，11（4）：045002.

[6]? ? CUI L，KIERNAN S，GILCHRIST M D. Designing the energy absorption capacity of functionally graded foam materials[J]. Materials Science and Engineering：A，2009，507（1/2）：215—225.

[7]? ? CUI L，F(xiàn)ORERO RUEDA M A，GILCHRIST M D. Optimisation of energy absorbing liner for equestrian helmets.Part II：Functionallygraded foam liner[J]. Materials & Design，2009，30（9）：3414—3419.

[8]? ?MOHAN P，PARK C K，MARZOUGUI D，et al. LSTC /NCAC Dummy Model Development[C]// 11th International LS-DYNA Users Conference，Livermore，California，USA，2013：53—64.

[9]? ?MASAMI I，YUKO N，ATSUTAKA T，et al. Development of advanced human models in THUMS[C]// 6th European LS-DYNA Users' Conference. Livermore，California，USA，2015：47—56.

[10]? PINNOJI P K，MAHAJAN P，BOURDET N，et al. Impact dynamics of metal foam shells for motorcycle helmets：Experiments & numerical modeling[J]. International Journal of Impact Engineering，2010，37（3）：274—284.

[11]? SCHRAAD M W，HARLOW F H. A stochastic constitutive model for disordered cellular materials：Finite-strain uni-axial compression[J]. International Journal of Solids and Structures，2006，43（11/12）：3542—3568.

[12]? KIERNAN S，GILCHRISTM D. Towards a virtual functionally graded foam：Defining the large strain constitutive response of an isotropic closed cell polymeric cellular solid[J]. International Journal of Engineering Science，2010，48（11）：1373—1386.

[13]? 運動頭盔 自行車、滑板、輪滑運動頭盔的安全要求和試驗方法： GB 24429—2009[S]. 北京：中國標準出版社，2010.

Sportshelmets-safety requirements and testing methods for sports helmets for cyclists and users of skateboards and roller skates：GB24429—2009[S]. Beijing：Standards Press of China，2010.（In Chinese）

[14]? Uniform provisions concerning the approval of protective helmets ECE R22.05[S]. Geneva：Economic Commission of Europe，2002.

[15]? ZHOU Z，JIANG B H，CAO L B，et al. Numerical simulations of the 10-year-old head response in drop impacts and compression tests[J]. Computer Methods and Programs in Biomedicine，2016，131：13—25.

[16]? 陳勇，楊濟匡. 轎車與兒童行人碰撞事故重建研究[J]. 湖南大學學報（自然科學版），2009，36（1）：35—40.

CHEN Y，YANG J K. Reconstruction of a passenger car-to-child pedestrian accident[J]. Journal of Hunan University （Natural Sciences），2009，36（1）：35—40. （In Chinese）

[17]? 汪黎.摩托車與汽車碰撞中騎車人的頭部響應和顱腦損傷研究[D]. 長沙：湖南大學，2019.

WANG L. A study on motorcyclists head responses and brain injuries during impacts against vehicle[D]. Changsha：Hunan University，2019. （In Chinese）

[18]? 水野幸治. 汽車碰撞安全[M]. 韓勇，陳一唯譯 . 北京：人民交通出版社，2016.

YUKI MIZUNO. Crash safety of passenger vehicles[M]. Beijing：China Communications Press，2016.