蘭鳳崇+黃偉+陳吉清+吳凱

摘 要：為了深入研究汽車行人碰撞過程中下肢的生物力學(xué)響應(yīng)和損傷機(jī)理，基于人體解剖學(xué)結(jié)構(gòu)建立了具有完整下肢組織結(jié)構(gòu)和高仿生精度的成年男性行人下肢有限元模型.包括股骨、脛骨、腓骨、髕骨等下肢骨骼以及皮膚、肌肉、韌帶、關(guān)節(jié)囊、半月板等重要軟組織.針對(duì)長骨骨干斷面幾何不均勻的特征，提出以CT斷面影像數(shù)據(jù)為依據(jù)，建立以真實(shí)皮質(zhì)骨內(nèi)外表面為邊界實(shí)現(xiàn)皮質(zhì)骨斷面厚度和形狀連續(xù)變化的長骨數(shù)值模型，對(duì)不同建模方式進(jìn)行了對(duì)比，采用兩層實(shí)體單元模擬皮質(zhì)骨以獲得相對(duì)較高的計(jì)算精度和計(jì)算效率.通過模擬相關(guān)生物力學(xué)實(shí)驗(yàn)，獲得了行人下肢各部位的損傷參數(shù)，分析了皮質(zhì)骨厚度變化以及不同撞擊方向?qū)ο轮珦p傷參數(shù)和損傷機(jī)理的影響.上述損傷參數(shù)的獲得可為我國汽安全性設(shè)計(jì)提供參考.

關(guān)鍵詞：行人安全；下肢；有限元模型；生物力學(xué)；損傷機(jī)理

中圖分類號(hào)：U461 文獻(xiàn)標(biāo)識(shí)碼：A

文章編號(hào)：1674-2974（2016）10-0042-10

Abstract： In order to study the biomechanical response and injury mechanisms of pedestrian`s lower limb during impact， a three-dimensional FE model of lower limb for adult pedestrians with high precision was developed on the basis of human anatomical structure. The lower limb model included complete anatomical structure of femur， tibia， fibula， patella as well as soft tissues such as skin， fresh， ligaments， capsule meniscus and cartilage. Considering the nonuniformity of the section of cortical bone， a long bone model was developed based on CT data with the thickness and the shape of the cortical bone section varying continuously. In comparison with other modeling methods， cortical bone was modeled with two-layer solid elements to obtain higher precision and efficiency. The injury criterions of pedestrians′ lower limb such as the ultimate bending moment of femur and tibia were obtained by modeling related biomechanical experiments. In addition， the influence of the thickness of cortical bone and the impact direction to injury mechanism and injury parameters of lower limb were analyzed. The injury parameters can provide important reference for the design of cars.

Key words：pedestrian safety； lower limb； finite element model； biomechanics； injury mechanism

中國是交通事故死亡人數(shù)最多的國家.相對(duì)于汽車乘員，行人是道路使用者中的弱勢群體.據(jù)世界衛(wèi)生組織統(tǒng)計(jì)，中國行人死亡人數(shù)占交通事故死亡人數(shù)的26%，而美國只占11%，德國為14%[1].由此可見，中國在行人安全方面與發(fā)達(dá)國家還存在較大的差距.如何避免和降低交通事故中行人的死傷，是中國道路安全領(lǐng)域面臨的重要課題，對(duì)其進(jìn)行研究，具有重大的社會(huì)意義.

在車人碰撞交通事故中，行人最易受傷的部位為頭部和下肢，所占比例分別為31.3%和32.4%[2].雖然頭部損傷是造成行人死亡的主要原因，但一般情況下，下肢往往首先與汽車發(fā)生碰撞，先于頭部發(fā)生損傷，下肢損傷是交通事故中使人致殘的主要因素.因此，下肢防護(hù)是行人保護(hù)研究的重點(diǎn).

鑒于交通事故中行人下肢損傷的頻發(fā)性和嚴(yán)重性，國內(nèi)外學(xué)者進(jìn)行了下肢損傷的相關(guān)研究.隨著計(jì)算機(jī)和醫(yī)學(xué)成像技術(shù)的不斷發(fā)展，有限元法已逐漸成為研究下肢損傷生物力學(xué)的重要工具，使下肢損傷和生物力學(xué)響應(yīng)的數(shù)值模擬成為可能.Yang等[3]建立了由股骨、脛骨和膝關(guān)節(jié)構(gòu)成的人體下肢有限元模型，用來預(yù)測下肢長骨和膝關(guān)節(jié)的損傷風(fēng)險(xiǎn).Takahashi等[4-5]分別基于Viewpoint Datalab和志愿者M(jìn)RI數(shù)據(jù)，建立了具有較高生物仿真度的行人下肢模型，并采用非線性、應(yīng)變率材料模型定義下肢長骨和膝關(guān)節(jié)韌帶.Untaroiu等[6]根據(jù)Visible Human Male Project（VHMP）的幾何數(shù)據(jù)，在縱、橫向上分別以0.95和0.85的縮放系數(shù)進(jìn)行縮放，建立了50百分位（美國）男性行人的下肢有限元模型，并模擬相關(guān)生物力學(xué)實(shí)驗(yàn)對(duì)下肢模型的可靠性進(jìn)行了全面的驗(yàn)證.張冠軍等[7]基于THUMS模型改進(jìn)了行人下肢模型，并對(duì)其進(jìn)行了整體彎曲驗(yàn)證，具有較好的生物逼真度.

中國在人體損傷生物力學(xué)方面的研究起步較晚，目前汽車企業(yè)和研究機(jī)構(gòu)自主開發(fā)可用于汽車安全性研究的相關(guān)生物力學(xué)模型還較少；且對(duì)于能夠準(zhǔn)確、深入分析行人下肢的生物力學(xué)響應(yīng)和損傷機(jī)理而言，相關(guān)模型在結(jié)構(gòu)和材料等方面做了較多簡化，如將下肢長骨模型皮質(zhì)骨各斷面賦予相同厚度，不符合骨骼的解剖學(xué)特征，會(huì)影響損傷的分析結(jié)果，因此建立精細(xì)化的行人下肢模型成為迫切的需要.本文根據(jù)人體解剖學(xué)結(jié)構(gòu)，建立了具有精細(xì)解剖學(xué)結(jié)構(gòu)的成年行人下肢有限元模型，通過模擬相關(guān)生物力學(xué)實(shí)驗(yàn)進(jìn)行了驗(yàn)證，獲得了下肢各部位的損傷參數(shù)，為中國汽車安全性設(shè)計(jì)提供參考.

1 行人下肢有限元模型的構(gòu)建

1.1 行人下肢有限元模型的建立

根據(jù)國標(biāo)GB 10000中50百分位中國成年男性的身體尺寸標(biāo)準(zhǔn)（身高1 678 mm，體重59 kg），本文選定一位30歲、身高1 680 mm、體重約60 kg的中國成年男性志愿者進(jìn)行下肢螺旋CT掃描，為獲得較高的掃描精度，掃描厚度定為0.6 mm.根據(jù)掃描得到的CT斷層文件，采用醫(yī)學(xué)影像技術(shù)，通過mimics軟件進(jìn)行圖像分割，提取得到下肢骨骼的三維點(diǎn)云數(shù)據(jù)，并利用逆向開發(fā)軟件重建其三維曲面模型.在幾何模型的基礎(chǔ)上，利用ANSYS ICEM CFD軟件及其特有的Block-Controlled網(wǎng)格劃分法構(gòu)建下肢骨骼三維有限元模型；由于韌帶、關(guān)節(jié)囊等下肢軟組織難以從CT影像中提取，本文通過研究其解剖學(xué)結(jié)構(gòu)確定下肢各軟組織的形態(tài)特征，在下肢骨骼模型的基礎(chǔ)上利用HYPERMESH有限元前處理軟件構(gòu)建下肢軟組織的有限元模型.

行人下肢有限元模型如圖1所示，該模型具備完整的下肢解剖學(xué)結(jié)構(gòu)，包括股骨、脛骨、腓骨、髕骨等下肢骨骼以及皮膚、肌肉和膝關(guān)節(jié)軟組織等.行人汽車碰撞事故中，下肢損傷的主要類型包括：大腿損傷（股骨骨折）、小腿損傷（脛骨、腓骨骨折）及膝關(guān)節(jié)損傷[8].因此，下肢長骨和膝關(guān)節(jié)是本文的建模重點(diǎn).

本文對(duì)長骨不同建模方式進(jìn)行了比較和分析，采用六面體實(shí)體單元模擬皮質(zhì)骨和松質(zhì)骨，其中骨干區(qū)域皮質(zhì)骨采用2層實(shí)體單元模擬，后續(xù)分析表明采用2層實(shí)體單元模擬長骨骨干能夠兼顧計(jì)算精度和計(jì)算效率；骨骺區(qū)域皮質(zhì)骨較薄，為避免單元尺寸過小，僅用1層單元模擬.骨干和骨骺皮質(zhì)骨之間采用五面體單元實(shí)現(xiàn)平滑過渡.以股骨模型為例，其有限元模型如圖2所示.

下肢長骨皮質(zhì)骨各斷面具有形狀不規(guī)則性和厚度不均勻性，同時(shí)各斷面形狀、厚度變化具有連續(xù)性.前人建立的長骨皮質(zhì)骨模型多是通過外表面向內(nèi)偏移而成，因而皮質(zhì)骨各斷面厚度相同；或采用殼單元模擬分段賦予皮質(zhì)骨不同的厚度值，但段間皮質(zhì)骨厚度的變化不連續(xù)且段內(nèi)厚度相同，因此上述方法均不能準(zhǔn)確地反映長骨皮質(zhì)骨的解剖學(xué)特征.本文通過CT影像技術(shù)，采用較小的掃描層厚以獲得精確反映下肢長骨復(fù)雜幾何特征的長骨點(diǎn)云數(shù)據(jù)，并擬合生成皮質(zhì)骨內(nèi)、外表面的自由曲面數(shù)學(xué)模型，以此為邊界建立皮質(zhì)骨厚度、形狀連續(xù)變化的長骨有限元模型，真實(shí)表征下肢長骨皮質(zhì)骨形狀的不規(guī)則性、厚度的不均勻性以及形狀、厚度變化的連續(xù)性.以脛骨為例，其CT影像和有限元模型分別如圖3，圖4所示.

膝關(guān)節(jié)是人體內(nèi)最大最復(fù)雜的關(guān)節(jié)，其有限元模型如圖5所示.膝關(guān)節(jié)軟組織模型包括韌帶、關(guān)節(jié)囊、半月板和關(guān)節(jié)軟骨等，對(duì)膝關(guān)節(jié)的穩(wěn)定性具有重要作用.影響膝關(guān)節(jié)運(yùn)動(dòng)學(xué)響應(yīng)的韌帶主要有4條，

分別為：脛側(cè)副韌帶（MCL），腓側(cè)副韌帶（LCL），前交叉韌帶（ACL）和后交叉韌帶（PCL）.上述韌帶均采用2層實(shí)體單元模擬，厚度約為4 mm[9].關(guān)節(jié)囊采用殼單元模擬，包圍膝關(guān)節(jié)四周（為顯示膝關(guān)節(jié)內(nèi)部軟組織結(jié)構(gòu)，圖6（b）隱藏部分關(guān)節(jié)囊單元），厚度定義為0.5 mm[6].在股骨、脛骨及髕骨的關(guān)節(jié)面上，分別附上一層實(shí)體單元模擬關(guān)節(jié)軟骨，厚度為2～3 mm[10].內(nèi)、外側(cè)半月板采用一層實(shí)體單元模擬，位于股骨和脛骨關(guān)節(jié)面之間.膝關(guān)節(jié)各關(guān)節(jié)面之間定義為自動(dòng)面面接觸方式（*CONTACT_AUTOMATIC_SURFACE_TO_SURF -ACE），關(guān)節(jié)各軟組織之間定義為單面接觸方式（*CONTACT_AUTOMATIC_SINGLE_SURFACE）.

肌肉與皮膚能夠吸收碰撞能量，改變骨骼的受力分布情況，具有一定的緩沖作用.肌肉采用六面體單元模擬，與長骨模型共節(jié)點(diǎn)連接，并在表面附上一層殼單元模擬皮膚，厚度定義為1 mm.

建立的行人下肢有限元模型共包括68 830個(gè)單元，73 440個(gè)節(jié)點(diǎn)，最小單元尺寸0.7 mm，最小雅各比0.5，滿足計(jì)算要求.

1.2 行人下肢有限元模型的材料參數(shù)

盡管不同人種身體在尺寸上存在差異，但其生物組織材料參數(shù)仍具有一定的借鑒意義，本文下肢模型的材料參數(shù)主要參考國內(nèi)外相關(guān)文獻(xiàn)獲得，如表1～表3所示.其中皮質(zhì)骨和松質(zhì)骨采用各向同性的彈塑性材料模擬；膝關(guān)節(jié)主要韌帶定義為線黏彈性材料，以模擬韌帶的非線性和黏彈性特征，相關(guān)參數(shù)參考Van Dommelen等不同應(yīng)變率下的韌帶拉伸實(shí)驗(yàn)反求獲得[15]；下肢肌肉參考文獻(xiàn)定義為粘彈性材料，皮膚、關(guān)節(jié)囊、關(guān)節(jié)軟骨和半月板等軟組織均采用彈性材料模擬，不考慮其損傷.

2 行人下肢損傷參數(shù)研究

損傷閾值的獲取是研究人體損傷風(fēng)險(xiǎn)和機(jī)制的基礎(chǔ).由于國內(nèi)在人體損傷生物力學(xué)方面的研究起步較晚，目前已發(fā)表的文獻(xiàn)資料中，缺少與中國人體下肢損傷閾值相關(guān)的基礎(chǔ)性實(shí)驗(yàn)研究.在此背景下，本文參考?xì)W美人體生物力學(xué)實(shí)驗(yàn)方法，利用本文建立的高精度下肢生物力學(xué)模型，對(duì)下肢各部位的損傷參數(shù)進(jìn)行了仿真研究.

2.1 下肢長骨模型準(zhǔn)靜態(tài)條件下的損傷參數(shù)研究

本文參考Yamada，Asang等[17-18]的尸體實(shí)驗(yàn)，在準(zhǔn)靜態(tài)加載條件下，對(duì)下肢長骨（股骨、脛骨和腓骨）的損傷參數(shù)進(jìn)行了仿真研究，如圖6所示.將下肢長骨模型水平放置于剛性支撐面上，定義接觸面摩擦系數(shù)為0.3，直徑為25 mm的剛性圓柱體以0.01 m/s的恒定速度分別從A-P方向?qū)﹂L骨模型中部進(jìn)行加載，直至骨折斷裂.

準(zhǔn)靜態(tài)三點(diǎn)彎曲條件下股骨、脛骨和腓骨模型的生物力學(xué)響應(yīng)分別如圖7～圖9所示.準(zhǔn)靜態(tài)加載作用下，接觸力隨加載處位移的增加而增大，并近似成正比例關(guān)系.股骨、脛骨和腓骨模型骨折的耐受極限分別為4.50 kN，4.03 kN和0.41 kN，單位位移下增加的載荷分別為0.39 kN/mm，0.39 kN/mm和0.03 kN/mm，其中股骨的耐受極限最高，約是腓骨的10倍，而脛骨的損傷閾值略低于股骨.仿真曲線與文獻(xiàn)中實(shí)驗(yàn)曲線對(duì)比可知，兩者具有相同趨勢，盡管存在一定的誤差，但考慮到實(shí)驗(yàn)樣本（歐美）與仿真模型（中國）在尺寸上的差異以及計(jì)算誤差，上述誤差在合理范圍內(nèi).由于模型仿真與尸體實(shí)驗(yàn)具有類似的生物力學(xué)響應(yīng)，因此認(rèn)為下肢長骨模型能夠反映長骨的準(zhǔn)靜態(tài)生物力學(xué)特性.

2.2 下肢長骨動(dòng)態(tài)加載條件下的損傷參數(shù)研究

下肢長骨動(dòng)態(tài)加載條件下的損傷參數(shù)研究參考文[19]的尸體實(shí)驗(yàn)方法進(jìn)行.仿真設(shè)置如圖10所示，將長骨模型兩端塞入下方帶有弧形金屬板的金屬方盒內(nèi)，以提供一個(gè)簡單的支撐條件，并保持模型姿態(tài)與其在人體中的姿態(tài)一致.實(shí)驗(yàn)中金屬盒采用聚氨酯泡沫填充塞實(shí)，仿真時(shí)通過定義長骨模型兩端與金屬盒的剛性連接（*CONSTRAINED_EXTRA_NODES_SET）來模擬.前端為弧形的剛性沖擊器以恒定速度從L-M方向加載于股骨（1.2 m/s）和脛骨（1.45 m/s）中部直至骨折.脛骨彎曲實(shí)驗(yàn)中，沖擊器前端包裹有25 mm的conforTM泡沫，其仿真材料特性來自文獻(xiàn)[20].

仿真輸出長骨模型中部彎矩與加載處位移的關(guān)系曲線，并與實(shí)驗(yàn)結(jié)果相對(duì)比，如圖11，圖12所示.在動(dòng)態(tài)加載條件下，股骨和脛骨模型中部承受的彎矩隨加載處位移的增加而增加，直至骨折，極限彎矩分別為417 N·m和293 N·m，分別與Kerrigan等股骨極限彎矩為412±102 N·m和Nyquist等脛骨極限彎矩為317±88 N·m[21]的實(shí)驗(yàn)結(jié)果相近.

隨著年齡的增長，人體皮質(zhì)骨內(nèi)側(cè)將出現(xiàn)皮質(zhì)骨小梁化，導(dǎo)致皮質(zhì)骨厚度不斷減薄，Ward等研究表明，50歲以上成年人皮質(zhì)骨減薄速度約為14%/10年[22].皮質(zhì)骨厚度減薄導(dǎo)致下肢長骨的耐受極限明顯下降.為了研究皮質(zhì)骨厚度對(duì)下肢長骨損傷參數(shù)的影響，本文借助hypermesh有限元前處理軟件，將皮質(zhì)骨模型內(nèi)側(cè)面節(jié)點(diǎn)向外等比例偏移，分別建立了不同皮質(zhì)骨厚度的股骨和脛骨模型，并進(jìn)行動(dòng)態(tài)三點(diǎn)彎曲仿真，得到不同皮質(zhì)骨厚度下肢長骨骨折的耐受極限，如表4所示.

為得到股骨和脛骨皮質(zhì)骨厚度減薄與其耐受極限下降程度的對(duì)應(yīng)關(guān)系，采用最小二乘法對(duì)上述數(shù)據(jù)進(jìn)行擬合，結(jié)果如圖13所示.股骨/脛骨擬合曲線的斜率分別為0.74/0.67，標(biāo)準(zhǔn)差分別為0.016/0.015.由此可以得出，下肢長骨皮質(zhì)骨厚度每減薄1%，其抗彎極限將下降約0.7%.

本文對(duì)比了動(dòng)態(tài)載荷下長骨皮質(zhì)骨不同建模方式對(duì)其力學(xué)響應(yīng)的影響，不同建模方式如圖14所示.圖15所示為股骨模型在動(dòng)態(tài)3點(diǎn)彎曲載荷條件下3種不同建模方式的力學(xué)特性對(duì)比.當(dāng)采用兩層和三層實(shí)體單元模擬股骨骨干皮質(zhì)骨時(shí)，兩者的彎矩與位移曲線具有相近的斜率，而當(dāng)僅采用一層單元模擬時(shí)，在相同的加載位移下，骨干中部的彎矩要低于前兩者.由此可見，僅采用一層實(shí)體單元不能較好的模擬長骨骨干皮質(zhì)骨的力學(xué)性能，而采用兩層和三層實(shí)體單元模擬具有相近的計(jì)算精度.因此，本文采用兩層實(shí)體單元模擬骨干皮質(zhì)骨，兼顧了計(jì)算精度和計(jì)算效率，具有一定的合理性.

2.3 大腿和小腿動(dòng)態(tài)加載條件下的損傷參數(shù)研究

本文參考Kerrigan等的尸體實(shí)驗(yàn)方法由側(cè)向（L-M）對(duì)帶有肌肉的大腿和小腿模型進(jìn)行了動(dòng)態(tài)三點(diǎn)彎曲仿真，獲得了大腿和小腿動(dòng)態(tài)加載條件下的損傷參數(shù)，仿真設(shè)置與下肢長骨模型的動(dòng)態(tài)三點(diǎn)彎曲仿真相同，如圖16所示.

仿真得到加載處接觸力與位移的關(guān)系曲線，并與文獻(xiàn)中實(shí)驗(yàn)結(jié)果進(jìn)行對(duì)比，如圖17，圖18所示.大腿動(dòng)態(tài)三點(diǎn)彎曲加載時(shí)，由于大腿皮膚與肌肉的緩沖作用，起始階段曲線上升較為平緩，隨著股骨受力，接觸力上升加快，當(dāng)接觸力達(dá)到耐受極限4.82 kN時(shí)，股骨斷裂；小腿動(dòng)態(tài)三點(diǎn)彎曲加載時(shí)，由于腓骨骨折，仿真曲線上升到2.22 kN后接觸力略有下降，然后繼續(xù)上升，直至脛骨骨折，此時(shí)接觸力達(dá)到最大值4.02 kN.大腿仿真曲線位于文獻(xiàn)中實(shí)驗(yàn)上下邊界之間，小腿仿真曲線初始階段存在一定偏差，這可能由腓骨的動(dòng)態(tài)特性差異所引起，但基本與實(shí)驗(yàn)上邊界相近.盡管上述實(shí)驗(yàn)結(jié)果是基于歐美人體實(shí)驗(yàn)所得，但從載荷上升趨勢和損傷閾值上分析，本文建立的大腿和小腿模型能真實(shí)地反映相關(guān)生物體的生物力學(xué)特性.

由于位于小腿外側(cè)的腓骨其損傷閾值明顯較低，同時(shí)考慮到小腿不同位置肌肉對(duì)骨骼的包裹厚度不均勻（如小腿前側(cè)基本沒有肌肉包覆），不同方向加載小腿的生物力學(xué)響應(yīng)有所不同.為研究不同加載方向?qū)π⊥葥p傷機(jī)理和損傷參數(shù)的影響，本文對(duì)比了外側(cè)、前側(cè)和后側(cè)3種不同方向加載條件下小腿的生物力學(xué)響應(yīng)，如圖19所示，仿真設(shè)置與上文相同.

撞擊塊從前側(cè)對(duì)小腿模型進(jìn)行加載時(shí)，由于小腿前側(cè)基本沒有肌肉包覆，撞擊塊直接作用于位于前側(cè)的脛骨，相比于側(cè)向加載，接觸力上升明顯較快，當(dāng)接觸力達(dá)到5.24 kN時(shí)，脛骨發(fā)生骨折，由于腓骨損傷閾值相對(duì)較低，腓骨隨之也發(fā)生骨折.而撞擊塊從后側(cè)對(duì)小腿模型進(jìn)行加載時(shí)，由于小腿后側(cè)肌肉包覆較厚，接觸力上升較為緩慢，當(dāng)接觸力達(dá)到3.52 kN時(shí)，位于小腿側(cè)后方腓骨首先發(fā)生骨折，接觸力經(jīng)歷短暫下降后隨著脛骨的承載而繼續(xù)上升，當(dāng)接觸力達(dá)到3.98 kN時(shí)，脛骨發(fā)生骨折.對(duì)比3種加載方式，由前側(cè)進(jìn)行加載時(shí)，盡管前側(cè)肌肉包覆較少，但脛骨首先承受載荷，因此損傷閾值要明顯高于其他兩種加載方式；而側(cè)向和后向加載時(shí)，腓骨首先承載并先于脛骨發(fā)生骨折，損傷閾值較低，但相比之下，后向加載時(shí)腓骨對(duì)應(yīng)的受力極限要高于前者，這可能是由于后側(cè)肌肉包覆較厚所致.由此可以看出，不同撞擊方向?qū)π⊥鹊膿p傷機(jī)理和損傷閾值均具有較大影響.

2.4 膝關(guān)節(jié)損傷參數(shù)研究

彎曲和剪切載荷被認(rèn)為是行人膝關(guān)節(jié)損傷的主要機(jī)理[23].本文參考文獻(xiàn)[24]的膝關(guān)節(jié)四點(diǎn)彎曲實(shí)驗(yàn)，對(duì)行人膝關(guān)節(jié)損傷參數(shù)進(jìn)行了研究.試驗(yàn)裝置如圖20所示，膝關(guān)節(jié)近心端和遠(yuǎn)心端分別與伸長桿固定，伸長桿通過轉(zhuǎn)動(dòng)軸承與支座連接，其中近心端支座為滑動(dòng)支座，遠(yuǎn)心端支座為固定支座.帶轉(zhuǎn)動(dòng)鉸的叉形沖擊桿加載于左右伸長桿上，使膝關(guān)節(jié)以1°/ms的速率彎曲，模擬行人下肢在一定沖擊載荷下的膝關(guān)節(jié)彎曲響應(yīng).

仿真參考實(shí)驗(yàn)設(shè)置，如圖21所示.膝關(guān)節(jié)模型近心端和遠(yuǎn)心端分別與左右伸長桿剛性連接，叉形桿、伸長桿和支座均定義為剛體，約束遠(yuǎn)心端支座所有自由度，僅釋放近心端支座左右移動(dòng)自由度.通過加載裝置對(duì)上述機(jī)構(gòu)進(jìn)行加載，使膝關(guān)節(jié)模型彎曲速度與實(shí)驗(yàn)一致，如圖22所示.

仿真輸出膝關(guān)節(jié)彎矩與膝關(guān)節(jié)彎曲角度的關(guān)系曲線，如圖23所示.

由于韌帶的約束作用，隨著膝關(guān)節(jié)彎曲角度的增大，其所受的彎矩也逐漸上升，當(dāng)膝關(guān)節(jié)彎曲角度達(dá)到14.5°時(shí)，所受彎矩達(dá)到最大值122 N·m，MCL韌帶隨之?dāng)嗔眩鐖D24所示，與文獻(xiàn)中實(shí)驗(yàn)結(jié)果相吻合.由于仿真時(shí)沒有考慮皮膚、肌肉以及關(guān)節(jié)囊對(duì)膝關(guān)節(jié)的約束作用，MCL韌帶斷裂后膝關(guān)節(jié)所受彎矩與實(shí)驗(yàn)相比下降較為迅速.仿真曲線與實(shí)驗(yàn)曲線無論是趨勢還是峰值都能較好地吻合，因此本文建立的膝關(guān)節(jié)模型可用來模擬人體膝關(guān)節(jié)的生物力學(xué)響應(yīng)和損傷細(xì)節(jié).由膝關(guān)節(jié)模型仿真結(jié)果可知，在1°/ms的加載速率下，行人膝關(guān)節(jié)的損傷閾值如極限彎矩和極限彎曲角度分別約為122 N·m和15°.

本文利用下肢生物力學(xué)模型，獲得的中國人體下肢（50百分位）各部位的損傷參數(shù)如表5所示.損傷參數(shù)的獲取對(duì)于中國汽車安全性設(shè)計(jì)與虛擬安全評(píng)價(jià)具有一定的參考意義.

3 結(jié) 語

1）基于人體解剖學(xué)特征建立了具有完整和精細(xì)解剖學(xué)結(jié)構(gòu)的成年行人下肢生物力學(xué)有限元模型，包括股骨、脛骨、腓骨、髕骨等下肢骨骼以及皮膚、肌肉和膝關(guān)節(jié)軟組織等，其中膝關(guān)節(jié)軟組織包括內(nèi)外側(cè)副韌帶、前后交叉韌帶、髕韌帶、半月板、關(guān)節(jié)囊以及關(guān)節(jié)軟骨等.

2）建立的下肢有限元模型具有較高精度.針對(duì)長骨骨干斷面幾何不均勻的特征，提出以骨骼斷面影像數(shù)據(jù)為依據(jù)，建立以皮質(zhì)骨內(nèi)外表面為邊界實(shí)現(xiàn)皮質(zhì)骨厚度、形狀連續(xù)變化的長骨數(shù)值模型；通過對(duì)不同建模方式進(jìn)行對(duì)比，采用兩層實(shí)體單元模擬長骨骨干皮質(zhì)骨，在保證計(jì)算效率的前提下獲得較高的建模精度.

3）通過模擬相關(guān)生物力學(xué)實(shí)驗(yàn)，獲得了50百分位中國成年人體下肢各部位在動(dòng)、靜態(tài)加載條件下的損傷參數(shù)，并分析了皮質(zhì)骨厚度對(duì)長骨損傷參數(shù)的影響：皮質(zhì)骨厚度每減薄1%，長骨抗彎極限下降約0.7%；此外，對(duì)比了行人小腿在不同方向撞擊下的損傷機(jī)理和損傷參數(shù)：前側(cè)加載條件下，小腿承載能力較好，而在側(cè)向及后向加載條件下，小腿損傷閾值較低且腓骨具有較高的損傷風(fēng)險(xiǎn).下肢各部位損傷參數(shù)的獲得可為汽車安全性設(shè)計(jì)提供參考.

參考文獻(xiàn)

[1] 劉庭志， 陳吉清.汽車行人保護(hù)開發(fā)與研究進(jìn)展[J].汽車實(shí)用技術(shù)，2012（1）：17-22.

LIU Ting-zhi， CHEN Ji-qing. The latest progress of the research for pedestrian protection of car[J]. Automobile Technology， 2012（1）： 17-22. （In Chinese）

[2] 鄭巍.行人碰撞腿部保護(hù)研究[J].城市車輛，2004（3）：26-29.

ZHENG Wei. A study on leg protection during collision with pedestrians[J]. Urban Vehicles， 2004（3）： 26-29. （In Chinese）

[3] YANG J K， WITTEK A， KAJZER J. Finite element model of the human lower extremity skeleton in a lateral impact [C]//Proceedings of the 1996 International IRCOBI Conference on the Biomechanics of Impact. Dublin Ireland， 1996： 377-388.

[4] TAKAHASHI Y， KIKUCHI Y， KONOSU A， et al. Development and validation of the finite element model for the human lower limb of pedestrians [J]. Stapp Car Crash Journal， 2000，（44）44：335-355.

[5] TAKAHASHI Y， KIKUCHI Y， MORI F， et al. Advanced FE lower limb model for pedestrians [C]//Proceeding of the 18th International Technical Conference on the Enhanced Safety of Vehicles（ESV）， Nagoya， Japan： US Department of Transportation， National Highway Traffic Safety Administration， 2003：1-14.

[6] UNTAROIU C， DARVISH K， CRANDALL J， et al. A finite element model of the lower limb for simulating pedestrian impacts [J]. Stapp Car Crash Journal， 2005（49）：157-181.

[7] 張冠軍， 魏嵬， 曹立波，等. 行人下肢有限元模型的建立與驗(yàn)證研究[J]. 湖南大學(xué)學(xué)報(bào)：自然科學(xué)版， 2015， 42（2）：48-54.

ZHANG Guan-jun， WEI Wei， CAO Li-bo， et al. Development and validation of finite element models for the lower limb of pedestrians [J]. Journal of Hunan University： Natural Sciences， 2015， 42（2）：48-54. （In Chinese）

[8] 韓勇， 楊濟(jì)匡， 李凡，等.汽車-行人碰撞中人體下肢骨折的有限元分析[J]. 吉林大學(xué)學(xué)報(bào)：工學(xué)版， 2011，41（1）：6-11.

HAN Yong， YANG Ji-kuang， LI Fan， et al. Finite element analysis of lower extremity fractures in vehicle pedestrian collision[J]. Journal of Jilin University： Engineering and Technology Edition， 2011，41（1）：6-11. （In Chinese）

[9] 張冠軍， 曹立波， 官鳳嬌，等. 行人膝關(guān)節(jié)韌帶的建模及驗(yàn)證研究[J].汽車工程， 2012，34（1）：57-61.

ZHANG Guan-jun， CAO Li-bo， GUAN Feng-jiao， et al. A study on the modeling and its validation for pedestrian`s knee ligaments[J]. Automotive Engineering， 2012， 34（1）：57-61. （In Chinese）

[10]董躍福，董英海，胡廣洪，等.準(zhǔn)確構(gòu)建個(gè)體化膝關(guān)節(jié)有限元解剖模型的方法[J].臨床骨科雜志， 2011，14（2）：198-203.

DONG Yue-fu， DONG Ying-hai， HU Guang-hong， et al. A method for accurate reconstruction of subject-specific anatomical finite element model of knee joint[J]. Journal of Clinical Orthopaedics， 2011， 14（2）：198-203. （In Chinese）

[11]張冠軍.行人下肢的碰撞損傷特性及相關(guān)參數(shù)研究[D].長沙：湖南大學(xué)機(jī)械與運(yùn)載工程學(xué)院，2009.

ZHANG Guan-jun. A study on characteristics of lower extremity and related parameters in vehicle-pedestrian crashes[D]. Changsha： College of Mechanical and Vehicle Engineering， Hunan University， 2009. （In Chinese）

[12]李海巖，顧玉龍，阮世捷，等.人體膝關(guān)節(jié)有限元?jiǎng)恿W(xué)分析模型的建立與驗(yàn)證[J].生物醫(yī)學(xué)工程學(xué)雜志，2012，29（1）：97-101.

LI Hai-yan， GU Yu-long， RUAN Shi-jie， et al. Development and validation of a finite element model of human knee joint for dynamic analysis[J]. Journal of Biomedical Engineering， 2012，29（1）：97-101. （In Chinese）

[13]DOMMELEN J A W， IVARSSON B J， JOLANDAN M M， et al. Characterization of the rate-dependent mechanical properties and failure of human knee ligaments[R].Warrenda PA：SAE International， 2005-01-0293.

[14]MO F， ARNOUX P J， MASSON D C C. The failure modelling of knee ligaments in the finite element model [J]. International Journal of Crashworthiness， 2012， 17（6）：630-636.

[15]LAN Feng-chong， HUANG Wei， CHEN Ji-qing， et al. Pedestrian knee joint modeling and viscoelastic material parameters identification for ligaments[C]// Eighth Conference on Measuring Technology and Mechatronics Automation.IEEE，2016：461-466.

[16]SNEDEKER J G， MUSER M H， WALZ F H. Assessment of pelvis and upper leg injury risk in car-pedestrian collisions： comparison of accident statistics， impactor tests and a human body finite element model [J]. Stapp Car Crash Journal， 2003，47（47）：437-57.

[17]YAMADA H. Strength of biological material[M]. Edited by EVANSFG，Baltimore： The Williams & Wilkins Company， 1970：19-280.

[18]ASANG E. Experimental biomechanics of the human leg： A basis for interpreting typical skiing injury mechanisms[J]. Orthopedic Clinics of North America， 1976， 7（1）：63-73.

[19]KERRIGAN J R， BHALLA A K， MADELEY N J， et al. Experiments for establishing pedestrian-impact lower limb injury criteria[C]//SAE 2004 World Erhibition.Detroit， USA：2003：1-19.

[20]DARVISH K K， TAKHOUNTS E G， MATHEWS B T， et al. A nonlinear viscoelastic model for polyurethane foams[R]. Washington DC：SAE International，1999-01-0299.

[21]NYQUIST G W， CHENG R， EL-BOHY A， et al. Tibia bending： strength and response [C]//Proc Warrendale.1985：99-112.

[22]WARD K A， ADAMS J E， HANGARTNER T N. Recommendations for thresholds for cortical bone geometry and density measurement by peripheral quantitative computed tomography[J]. Calcified Tissue International， 2005， 77（5）： 275-280.

[23]楊濟(jì)匡.汽車與行人碰撞中的損傷生物力學(xué)研究概覽[J].汽車工程學(xué)報(bào)，2011， 1（2）：81-93.

YANG Ji-kuang. Overview of research on injury biomechanics in car-pedestrian collisions[J]. Chinese Journal of Automotive Engineering， 2011， 1（2）：81-93. （In Chinese）

[24]BOSE D， BHALLA K， ROOIJ L， et al. Response of the knee joint to the pedestrian impact loading environment [R].Washington DC：SAE International，2004-01-1608.