胸部鈍性沖擊時刻對心臟損傷的影響

蘭鳳崇　童芳　陳吉清　李雄

摘? ?要：為研究與心動周期有關(guān)的胸部沖擊時刻對鈍性心臟損傷的影響，開發(fā)并驗證了一個精細(xì)人體心臟生物力學(xué)模型，并針對40例處于不同心動周期的心臟進(jìn)行胸部沖擊仿真. 心臟模型幾何來源于醫(yī)學(xué)影像，包含四個心腔、瓣膜以及主動脈、肺動脈等結(jié)構(gòu). 通過Abaqus中基于面的流體腔方法考慮沖擊過程中血液與心臟的流固耦合作用. 心臟的不同狀態(tài)采用各心腔內(nèi)壓力與瓣膜的激活與否來表征. 仿真結(jié)果顯示：（1）建立的心臟模型在胸部正面沖擊下心內(nèi)血壓曲線在試驗范圍內(nèi)，能夠正確反應(yīng)心臟的沖擊響應(yīng). （2）乘員胸部正面沖擊下，左心房的血壓峰值為（164.91±17.33） kPa，明顯高于右心腔，使二尖瓣損傷風(fēng)險高于三尖瓣;右心心肌的應(yīng)力值為（1887.07±168.74） kPa，遠(yuǎn)大于左心，導(dǎo)致右心更容易發(fā)生心肌破裂. （3）心臟處于心室充盈期受到?jīng)_擊時，心肌應(yīng)力為（1901.3±150.7） kPa明顯高于其他時刻，更容易發(fā)生心臟損傷. （4）心腔內(nèi)的初始血壓對碰撞中的峰值血壓影響不明顯，Pearson相關(guān)系數(shù)小于0.2;而瓣膜的開閉狀態(tài)對心肌的損傷風(fēng)險影響較大. 研究結(jié)果有助于進(jìn)一步理解鈍性心臟損傷機(jī)理，為車輛安全設(shè)計提供基礎(chǔ).

關(guān)鍵詞：生物力學(xué);鈍性心臟損傷;沖擊時刻;心動周期;流固耦合

中圖分類號：U467.14? ? ? ? ? ? ? ? ? ? ? ? ? ? ?文獻(xiàn)標(biāo)志碼：A

Effect of Blunt Thoracic Impact Timing on Heart Injury

LAN Fengchong，TONG Fang，CHEN Jiqing LI Xiong

（School of Mechanical &Automotive Engineering，South China University of Technology，Guangzhou 510641，China）

Abstract：To explore the effect of thoracic impact timing related with the cardiac cycle on blunt cardiac injury，a biomechanics model of the heart was developed and validated. Simulations of thoracic impact with 40 hearts in different phases of the cardiac cycle were conducted. The geometry of heart model was obtained from medical imaging，including the structures of four heart chambers，heart valves and vessels etc. The fluid-structure interaction between heart and blood was considered by the method of surface-based fluid cavities in Abaqus. The different phrases of hearts were presented by the intracardiac pressure and valve activation. The simulation results show that： （1）the established heart model can present the correct impact response of the heart under thoracic blunt from the curves of intracardiac blood pressure in the testing ranges. （2）When the passenger thoracic was under the blunt，the peak of intracardiac blood pressure in the left atrium was （164.91±17.33） kPa，which was obviously higher than that of the right heart，causing that the mitral valve is more vulnerable than the tricuspid valve; the stress on the right heart was（1887.07±168.74） kPa，which was higher than the left heart，so the heart rupture is more likely to? occur on the right heart. （3）When the heart was impacted during ventricular filling period， the stress on the myocardium was（1901.3±150.7） kPa， which was higher than that during other impact periods. （4）Initial intracardiac pressure provided little effect on peak pressure， and the Pearson correlation coefficient was less than 0.2; while the state of the atrioventricular valves， open or closed， affects the myocardial injury a lot. This study is helpful to further understand the mechanism of blunt heart injury and provide a basis for vehicle safety design.

Key words：biomechanics;blunt heart injury;impact timing;cardiac cycle;fluid-structure interaction

在交通事故中，乘員胸部極易與方向盤、安全氣囊及汽車內(nèi)飾發(fā)生碰撞[1]，由此引起的乘員鈍性心臟損傷（Blunt cardiac injury， BCI）的死亡率極高[2]. 數(shù)據(jù)顯示，事故導(dǎo)致的心臟鈍性損傷患者有6.3%死于交通事故現(xiàn)場，即使送往醫(yī)院總死亡率也高達(dá)88.5%[3]. 因此，心臟鈍性損傷機(jī)理研究對于車輛設(shè)計中人體安全防護(hù)裝置的優(yōu)化，以及機(jī)械暴力事故中人體心臟的診斷和治療具有重要意義.

胸部鈍性沖擊下影響B(tài)CI的因素有很多，沖擊胸部的物體質(zhì)量和沖擊速度越大，會使胸腔發(fā)生越劇烈的變形，增加了BCI的風(fēng)險[4-5]. 由人種、性別、年齡等因素引起的胸壁順應(yīng)性差別在鈍性胸部沖擊下會對心臟造成不同程度的壓迫[2，6]. 老年人群肋骨斷裂的極限應(yīng)力應(yīng)變較低，在碰撞過程中更容易發(fā)生肋骨骨折，增加心臟破裂的可能性[7]. 不僅如此，心臟是一個充滿血液的四腔室結(jié)構(gòu)，心肌在心電刺激下持續(xù)地有節(jié)律的收縮，實現(xiàn)了血液的體循環(huán)和肺循環(huán)[8]. 在這個過程中，心肌的形態(tài)以及心腔內(nèi)的血流動力學(xué)均處于動態(tài)變化，心動周期內(nèi)不同胸部沖擊時刻對心臟損傷產(chǎn)生的影響不可忽略[9-10]. 許多研究者利用活體動物實驗來研究胸部沖擊時刻對鈍性心臟損傷的影響. Link等[11]利用18只幼豬做不同心動周期時刻的胸部沖擊試驗，發(fā)現(xiàn)心室纖顫的發(fā)生與否與胸部的撞擊時刻有關(guān). 劉維永等[12]利用78只麻醉后的家兔做心室前壁的撞擊實驗，結(jié)果顯示在舒張末期心臟破裂的幾率較高. Stein等人以及Life等人[13-14]通過對狗的胸部沖擊試驗發(fā)現(xiàn)，心臟收縮期時心室更容易受傷，而在心臟舒張期時心房更容易受傷. 然而，人與動物的生理和解剖學(xué)上存在差異，動物試驗不能完全真實地反應(yīng)人體的響應(yīng)[15]. 在這種情形下，基于人體解剖結(jié)構(gòu)建立的生物力學(xué)模型已成為分析人體損傷機(jī)理的重要方法之一.

目前的損傷生物力學(xué)模型如福特的整人模型[16]、豐田的安全人體模型（Total Human Model for Safety， THUMS）[17]、華南理工大學(xué)的中國人體模型（Chinese Human Body Model， CHUBM）[18]、全球人體聯(lián)盟模型（Global Human Body Models Consortium， GHBMC）[19]等，在人體碰撞損傷研究方面都取得了重要進(jìn)展. 然而，目前人體生物力學(xué)模型的建立基本不包含心臟的精細(xì)化建模，對心臟的鈍性沖擊損傷的研究較少. Jiang等[20]以及Iwamoto等[17]建立的人體胸部有限元模型中，心臟被看作一個均質(zhì)的實體，主要通過胸部壓縮量及肋骨斷裂情況研究胸腔整體的響應(yīng). 蘭鳳崇等[21]建立的胸部模型中，基于網(wǎng)格劃分情況將心臟分出兩個腔室，但對心房心室輪廓描述較為粗糙. 方紅榮等[22]采用尸體的心臟斷層圖像數(shù)據(jù)建立了心臟有限元模型，研究了血壓作用下心肌的應(yīng)力分布. 但其中血壓采用恒定的均布壓力載荷來模擬，對于血壓的動態(tài)變化模擬不足，無法反映血液與心肌的相互作用.

為了考慮碰撞過程中血液與心肌的耦合作用，并進(jìn)一步分析與心動周期有關(guān)的碰撞時刻對心臟損傷的影響，本文首先基于成年男性胸部的增強(qiáng)CT影像構(gòu)建了心臟的有限元模型，在Abaqus仿真軟件中對心房心室進(jìn)行血液流體的填充. 將血液-心臟流固耦合模型裝配到團(tuán)隊前期開發(fā)的中國50百分位男性模型胸腔內(nèi)，并利用尸體胸部的碰撞試驗進(jìn)行了模型有效性的驗證. 隨后，對心動周期內(nèi)40個不同狀態(tài)的心臟模型進(jìn)行胸部正面碰撞的仿真，得到了心臟舒張、收縮期胸部碰撞下心肌的應(yīng)力應(yīng)變分布以及四個心腔內(nèi)的血壓變化. 分析了四個腔室在胸部正面碰撞下的損傷風(fēng)險，碰撞時刻對心肌應(yīng)力、血壓的影響，以及瓣膜開閉狀態(tài)對BCI的影響.

1? ?模型與方法

1.1? ?乘員心臟有限元模型的建立

根據(jù)《中國成年人人體尺寸（GB/T 10000-1988）》，18-60歲的中國男性平均身高為167.8 cm，體重為59 kg. 據(jù)此，選擇身高170 cm，體重60 kg的30歲健康男性志愿者，對其注射造影劑后進(jìn)行增強(qiáng)CT掃描. 將獲得的影像數(shù)據(jù)導(dǎo)入Mimics軟件中進(jìn)行心肌及血管組織的點云處理. 利用閾值調(diào)整、區(qū)域增長、3D編輯等功能配合手動調(diào)整蒙罩部分區(qū)域分別獲取心室、心房、主動脈、心肌外輪廓的點云數(shù)據(jù). 提取血管的中心線，輔助調(diào)整血管與組織的局部干涉. 將提取的點云導(dǎo)入Geomagic Studio軟件中進(jìn)行幾何重建，生成的幾何模型通過Hypermesh軟件進(jìn)行網(wǎng)格劃分，得到心臟及血管的有限元模型如圖1所示.

建立的心臟模型包括心包膜、心肌、左心房、左心室、右心房、右心室、主動脈、肺動脈、肺靜脈、上下腔靜脈等. 心包膜和心內(nèi)膜采用三角形殼單元模擬，厚度為0.5 mm[23];血管用三角形殼單元模擬，厚度為1.5 mm[24];心肌因形狀復(fù)雜采用四面體單元劃分. 最終建立的模型具有154 202個單元和34 149個節(jié)點. 網(wǎng)格質(zhì)量顯示所有單元雅各比大于0.8，翹曲度小于4.7°且扭曲都小于37°. 心肌采用粘彈性材料模擬，如公式（1）所示：

G（t） = G_∞+ （G_o- G_∞）e^-βt （1）

式中：G_o是短效剪切模量，G_∞是長效剪切模量，用于控制粘彈性響應(yīng)，β為延遲系數(shù). 在心肌材料中剪切模量分別取0.44以及0.15[16].

血管采用線彈性材料模擬，更多部分的材料參數(shù)見表1. 其中ρ表示材料密度，E表示彈性模量，γ表示材料泊松比，H表示殼單元的厚度.

1.2? ?血液流體單元的填充

心臟是一個充滿血液的四腔室結(jié)構(gòu). 軟件Abaqus中基于面的流體腔算法能夠模擬在一個充滿流體的結(jié)構(gòu)中流體與結(jié)構(gòu)的相互作用. 對腔體內(nèi)的流體，用兩種模型來表達(dá)其性質(zhì)：基于理想氣體方程的氣動模型與表示不可壓縮或近似不可壓縮流體的水力模型. 血液是一種近似不可壓縮流體，故本研究選用水力模型來模擬血液的性質(zhì). 在Abaqus的顯示計算中，流體的壓力是根據(jù)腔體的體積來計算的. 通過設(shè)置流體的體積模量來定義其壓縮性，本研究中血液的體積模量為2 200 MPa.

腔體的體積等于組成腔體的每一個流體單元的體積總和. 單個流體單元是由腔體邊界的面單元與腔內(nèi)的參考點構(gòu)成. 圖2以右心室為例展示了一個流體腔以及內(nèi)部的流體單元. P作為右心室內(nèi)部的參考點，用于計算腔體體積以及輸出流體的壓力等性質(zhì);三角形ABC是組成右心室內(nèi)部邊界的面單元;四面體P-ABC構(gòu)成一個流體單元.

建立的血液-心臟流固耦合模型中，血液填充的四個腔室分別為左心房、左心室、右心房和右心室. 后續(xù)可根據(jù)心臟的狀態(tài)分別在四個心腔的參考點定義腔體初始壓力以及瓣膜的激活狀態(tài).

1.3? ?血液-心臟流固耦合模型的裝配及驗證

本團(tuán)隊前期已開發(fā)了中國50百分位成年男性乘員的生物力學(xué)模型，并對比文獻(xiàn)中的沖擊器試驗以及臺車試驗分別進(jìn)行了頭、頸、胸腹、上下肢等的局部和整體仿真[25-28]，驗證了模型的可靠性. 將建立的心臟模型與整人模型進(jìn)行裝配，如圖3所示. 微調(diào)肺、肝臟等與心臟的相對位置符合解剖學(xué)結(jié)構(gòu). 肺動脈、肺靜脈等與肺的連接方式采用共節(jié)點以提升計算的穩(wěn)定性.

在Kroell[29]的研究中，采用23具不同年齡性別的尸體樣本進(jìn)行正面胸部鈍性沖擊試驗，研究胸部及心臟在此過程中的動態(tài)響應(yīng). 為了模擬體內(nèi)的血壓狀態(tài)，在實驗前給心臟和主動脈填充生理鹽水并縫合，記錄實驗過程中各腔室的時間壓力變化曲線. 實驗中，樣本上肢被膠帶固定在一個5 cm寬的平面上（碰撞前撕開），用一個半徑為15.24 cm的圓柱（內(nèi)附加不同質(zhì)量）以不同的速度撞擊胸部第四根肋骨的位置，如圖4所示. 研究中，有9例試驗記錄了血壓變化，同時無背部約束及死后僵直現(xiàn)象. 按照沖擊器質(zhì)量和沖擊速度將這9例試驗分為2組，分別對模型在不同工況下的仿真響應(yīng)進(jìn)行驗證見表2，其中m代表沖擊器質(zhì)量，v代表沖擊速度.

本文建立了與尸體實驗相同的邊界條件來模擬胸部的鈍性沖擊過程. 將裝配有流固耦合心臟模型的整人模型調(diào)整姿勢到與實驗一致. 人體模型坐于剛性板上，上肢舉起約45°角. 約束沖擊器除水平移動方向的其他5個自由度，保證其不發(fā)生翻轉(zhuǎn).

1.4? ?不同胸部沖擊時刻的碰撞仿真

一個心動周期由四個階段組成，分別是等容收縮期、射血期、等容舒張期和心室充盈期，在每個階段選取10個心臟的瞬時狀態(tài)來進(jìn)行碰撞模擬. 心臟的瞬時狀態(tài)用各心房心室的血壓以及瓣膜的開閉情況來表征. 當(dāng)心臟處于等容收縮期、射血期和等容舒張期時，房室瓣膜處于關(guān)閉狀態(tài)，此時激活瓣膜與血液的接觸作用（圖5（a）），阻止血液在房室間流動;當(dāng)心臟處于心室充盈期時，取消瓣膜與血液的接觸作用（圖5（b）），使得心房和心室連通. 不同時刻的心腔內(nèi)血壓來源于文獻(xiàn)記錄，如圖6所示[8]. 例如，在射血期的仿真a根據(jù)圖6分別賦予四個心腔碰撞時的初始壓力15.8 kPa（LV），0.11 kPa（LA），3.16 kPa（RV）以及0.08 kPa（RV），并激活瓣膜與血液的接觸. 將40個不同狀態(tài)的心臟放入整人模型中，建立如圖4中仿真的約束條件，沖擊器質(zhì)量為18.98 kg，沖擊速度為7.2 m/s. 記錄不同狀態(tài)的心臟在胸部碰撞過程中心肌應(yīng)力應(yīng)變分布及心房心室的血壓輸出曲線.

2? ?結(jié)? ?果

2.1? ?心臟模型的驗證

兩組不同工況的胸部沖擊試驗及對應(yīng)仿真中，左心室血壓隨時間的變化曲線如圖7所示，右心房血壓隨時間的變化曲線如圖8所示. 針對左心室的壓力，實驗和仿真曲線在開始的10 ms變化不大，隨后10 - 20 ms迅速增大. 達(dá)到峰值后的約25 ms在高壓范圍內(nèi)波動，隨后壓力逐漸減小至30 kPa以下. 仿真1中的峰值血壓為72.4 kPa，仿真2中的峰值血壓為125.5 kPa，均在實驗數(shù)據(jù)范圍內(nèi)，符合實驗結(jié)果.

圖8所示右心房壓力曲線，經(jīng)過10 ms的胸廓擠壓階段（此時心臟未被壓縮）后，有波動的上升. 壓力達(dá)到最大之后即開始緩慢下降，在50 - 60 ms間幾乎恢復(fù)到初始值. 仿真1的右心房壓力峰值為72.1 kPa，仿真2中右心房壓力峰值為128.6 kPa，均在實驗數(shù)據(jù)范圍內(nèi)，符合實驗結(jié)果.

2.2? ?胸部正面沖擊下心肌及血壓的整體響應(yīng)

首先分析40個不同心動周期的心臟受到?jīng)_擊時的整體響應(yīng)，如圖9所示. 由圖可見不論心臟處于什么狀態(tài)，胸壁的正面沖擊都大大增加了心內(nèi)血壓及心肌的應(yīng)力. 在碰撞過程中，左心的血壓明顯高于右心的血壓，其中左心房的血壓最大達(dá)到（164.91±17.33） kPa. 對于心肌應(yīng)力而言，右心的心肌比左心的心肌承受了更大的應(yīng)力. 右心室的心肌上應(yīng)力最大，達(dá)到（1887.07±168.74） kPa.

選擇其中一例仿真結(jié)果（心臟處于等容收縮期時的碰撞，心腔初始壓力為左心室6.01 kPa，左心房1.00 kPa，右心室1.20 kPa，右心房0.75 kPa）來分析心肌及血管的應(yīng)力分布，如圖10所示. 當(dāng)胸腔被極度壓縮時（20 ms），Von Mises應(yīng)力主要集中在肺靜脈與左心房血管與心房的連接處，以及右心房與靜脈腔的連接處（圖10（a））. 這是由于胸部在受到撞擊時，心臟與肺等其他組織產(chǎn)生相對運動，導(dǎo)致血管在連接處產(chǎn)生拉扯. Turan[2]的研究中也報道了類似的鈍性心臟損傷機(jī)理：車輛碰撞導(dǎo)致人體突然的減速會使血管與心房的連接處斷裂.

胸部鈍性碰撞過程中，心肌的Von Mises應(yīng)力分布如圖10（b）所示. 當(dāng)胸骨開始接觸心臟時（9 ms），首先對右心產(chǎn)生沖擊，瞬間的作用力容易造成右心室的破裂. 載荷沖擊下胸骨持續(xù)擠壓心臟，20 ms時應(yīng)力集中在左右心室間隔處以及右心的內(nèi)壁，使得右心房心室以及心室間隔容易發(fā)生損傷. Kroell[29]利用尸體做的胸部正面沖擊實驗中，試驗93號，94號以及104號等也都產(chǎn)生了左右心腔間隔的破裂損傷.

2.3? ?胸部沖擊時刻對鈍性心臟損傷的影響

將心臟在不同狀態(tài)受到?jīng)_擊時的心肌最大應(yīng)力與峰值血壓繪制成點狀圖如圖11所示. 四種不同的點標(biāo)記分別代表等容收縮期、射血期、等容舒張期以及心室充盈期. 結(jié)果顯示，在心臟的四個腔室中，不論是峰值血壓還是心肌最大應(yīng)力，心臟處于心室充盈期的碰撞結(jié)果與其它三個狀態(tài)的結(jié)果具有明顯的差異性：在左心室中，心室充盈期的碰撞結(jié)果具有更高的峰值血壓與更大的心肌應(yīng)力;在左心房和右心室中，心室充盈期碰撞時產(chǎn)生了較低的血壓峰值與較大的心肌應(yīng)力;在右心房中，心室充盈期碰撞引起的腔內(nèi)血壓較低. 等容收縮期、射血期和等容舒張期碰撞時，應(yīng)力和血壓結(jié)果在一定范圍內(nèi)的分布較為平均，沒有明顯的規(guī)律.

當(dāng)心臟處于等容收縮期、射血期和等容舒張期時，房室瓣膜均處于關(guān)閉狀態(tài)，心腔內(nèi)的血液壓力是區(qū)分心臟狀態(tài)的主要變量. 圖11的結(jié)果表明，心動周期內(nèi)不同階段的心腔內(nèi)血壓對心臟鈍性損傷的影響不明顯. 為了進(jìn)一步探究每個心腔內(nèi)的初始血壓對該心腔的損傷影響，用Pearson相關(guān)系數(shù)來衡量它們之間的關(guān)系如圖12所示. 第一行的圖12（a）-（d）是四個心腔內(nèi)初始血壓與碰撞期間峰值血壓的相關(guān)系數(shù)計算;第二行的圖12（e）-（h）是四個心腔內(nèi)初始血壓與心肌最大應(yīng)力的相關(guān)系數(shù)計算. 從8組數(shù)據(jù)計算得到的Pearson相關(guān)系數(shù)表明，心腔內(nèi)的初始壓力與胸部撞擊時心內(nèi)的峰值之間幾乎沒有線性關(guān)系;在左心房中，初始血壓和心肌最大應(yīng)力的相關(guān)性系數(shù)為0.47，理論上存在中等偏弱的線性相關(guān)性. 這可能是由于左心房壁厚較薄，且在碰撞過程中直接與脊柱接觸，對血壓的敏感性稍強(qiáng).

當(dāng)心臟處于心室充盈期時，房室瓣膜處于開啟狀態(tài). 當(dāng)心臟處于等容收縮期、射血期和等容舒張期時，房室瓣膜均處于關(guān)閉狀態(tài). 為了研究房室瓣膜開閉狀態(tài)對心臟鈍性損傷的影響，將40個仿真結(jié)果分為兩組：第一組是瓣膜關(guān)閉狀態(tài)的心臟損傷結(jié)果;第二組是瓣膜開啟狀態(tài)的心臟損傷結(jié)果，對比如圖13所示. 房室瓣膜的開啟降低了心臟沖擊過程中左心房、右心房和右心室的血壓，增大了左心室的血壓;對于心肌的應(yīng)力分布而言，開啟的房室瓣膜增加了四個腔室的心肌應(yīng)力最大值. 從數(shù)值上看，第二組的心肌應(yīng)力最大值在左心室、左心房、右心室和右心房分別高于第一組15.99%、13.35%、15.00%以及4.00%;在血壓峰值方面，第二組的左心房（15.59%）、右心室（12.77%）、右心房（12.26%）較第一組低，而左心室的壓力高于第一組17.19%.

3? ?討? ?論

心臟在正常生理條件下的動態(tài)特性使得研究者們開始關(guān)注心動周期對BCI的影響，但目前的研究多是利用動物實驗進(jìn)行的，缺乏對人體心臟解剖學(xué)結(jié)構(gòu)的反映[19-22]. 本文利用CT影像建立了人體的血液-心臟流固耦合模型，與前期開發(fā)的人體模型進(jìn)行裝配后，利用尸體的正面胸部沖擊器試驗對其進(jìn)行了有效性驗證. 通過賦予各心腔不同的血壓及瓣膜的開閉狀態(tài)，模擬了心動周期內(nèi)40個不同狀態(tài)的心臟，并分別進(jìn)行胸部鈍性沖擊仿真，記錄了碰撞過程中各心腔的血壓變化及心肌應(yīng)力分布.

對于心臟整體響應(yīng)而言，胸部正面沖擊下右心心肌應(yīng)力明顯大于左心，說明在給定載荷下右心更容易受傷. 這一結(jié)論與Turk等[7，30]的研究結(jié)果一致，他們的基于尸檢的統(tǒng)計顯示BCI患者中右心破裂的數(shù)量是左心破裂數(shù)量的2-3倍. 造成這一現(xiàn)象的原因是胸部正面沖擊下，胸骨直接將碰撞力傳遞到右心，尤其是右心室，導(dǎo)致接觸區(qū)域出現(xiàn)局部大變形. Leavitt等[10]也描述了這一損傷機(jī)理，他認(rèn)為心腔的破裂常常由于胸骨的直接沖擊引起的. 由于左心與脊柱的接觸面積小于右心與胸壁的接觸面積，心臟受到胸壁和脊柱的擠壓時，左心尤其是左心房內(nèi)會產(chǎn)生更高的血壓，增加了二尖瓣及主動脈瓣的損傷風(fēng)險.

心內(nèi)壓力對鈍性心臟損傷的影響不明顯. 心臟處于等容收縮期、射血期和等容舒張期發(fā)生胸部鈍性沖擊時，心腔內(nèi)的血液壓力峰值與心肌最大應(yīng)力在一定范圍內(nèi)的分布沒有明顯規(guī)律. 不僅如此，pearson相關(guān)系數(shù)顯示各心腔的初始壓力也與相應(yīng)的心腔損傷幾乎沒有關(guān)系. 引起這種現(xiàn)象的一個原因可能是，在給定載荷的胸部沖擊下，心內(nèi)血壓能夠達(dá)到的峰值遠(yuǎn)遠(yuǎn)高于心臟在正常工作下的血壓，所以心臟正常工作的初始血壓對鈍性心臟損傷的影響很小. 如受到?jīng)_擊時，左心室的峰值血壓能夠達(dá)到130 kPa，而正常工作下峰值血壓約為15 kPa.

與心內(nèi)壓力不同，房室瓣膜的開閉狀態(tài)顯著影響了鈍性心臟損傷. 心臟處于心室充盈期（瓣膜開啟）受到鈍性沖擊時，心肌應(yīng)力明顯高于其他時期，具有更高的損傷風(fēng)險. Leavitt等以及Huis等基于BCI病人的損傷結(jié)果分析同樣認(rèn)為心臟破裂更容易發(fā)生在舒張末期[10，32].? 推測由于心臟處于心室充盈期時，心腔被血液完全填滿，碰撞過程中血液對心肌的作用力更大. 不僅如此，開啟的房室瓣膜允許血液在腔內(nèi)流動，使心肌整體在碰撞力下更容易變形，對于室間隔及內(nèi)壁結(jié)構(gòu)的局部擠壓力增加，產(chǎn)生更大的應(yīng)力. 另一方面，瓣膜的精細(xì)結(jié)構(gòu)使得它比心肌更容易受到血壓的影響[31]. 本文的結(jié)果表明心臟受到?jīng)_擊時，開啟的瓣膜在一定程度上能夠平衡房室的血壓差，減小房室瓣膜的損傷風(fēng)險. McDonald等人[33]的研究也表明，房室瓣膜的損傷更容易在收縮早期發(fā)生，此時房室瓣膜處于關(guān)閉狀態(tài).

本文的研究內(nèi)容也存在一定的局限性. 心臟在不同心動周期時刻的差別不僅有血壓、瓣膜狀態(tài)，還有心臟的形態(tài)、心肌在電信號刺激下收縮力的變化. 此外，本研究僅考慮了正面胸部沖擊下心臟的鈍性損傷，不同的外界載荷可能會使心臟呈現(xiàn)不同的損傷特點，這也是下一步研究的方向.

4? ?結(jié)? ?論

本文基于CT影像建立了人體心臟有限元模型，與之前開發(fā)的整人模型進(jìn)行裝配后，利用尸體的胸部沖擊試驗驗證了模型的有效性. 結(jié)合流固耦合數(shù)值算法，進(jìn)行了40例心臟處于不同心動周期的胸部鈍性沖擊仿真，分別研究了胸部鈍性沖擊下心臟的整體響應(yīng)以及心動周期對BCI的影響. 結(jié)果表明：胸腔受到?jīng)_擊時，左心腔的血壓高于右心腔，使二尖瓣的損傷風(fēng)險高于三尖瓣;右心的心肌應(yīng)力大于左心，導(dǎo)致心臟的右心室、右心房比左心更容易破裂. 當(dāng)心臟處于心室充盈期受到?jīng)_擊時，心肌應(yīng)力明顯高于其他時期的碰撞，更容易發(fā)生鈍性心臟損傷. 心腔內(nèi)的血壓對心肌損傷結(jié)果影響不明顯，對瓣膜損傷影響較大;而瓣膜的開閉狀態(tài)對心肌損傷的影響較大，開啟的瓣膜使心肌受到了更大的應(yīng)力，承受了較高的損傷風(fēng)險. 本文的研究結(jié)果為進(jìn)一步理解胸部鈍性沖擊下的心臟損傷提供了基礎(chǔ)，對于BCI的診斷和治療具有積極作用.

參考文獻(xiàn)

[1]? ? 韓勇，楊濟(jì)匡，水野幸治. 行人在與汽車碰撞中胸部動力學(xué)響應(yīng)和損傷機(jī)理的研究[J]. 汽車工程，2015，37（5）：516—519.

HAN Y，YANG J K，KOJI M. A study on the dynamic response and injury mechanism of human chest in vehicle-pedestrian collisions[J]. Automotive Engineering，2015，37（5）：516—519. （In Chinese）

[2]? ? TURAN A A，KARAYEL F A，AKYILDIZ E，et al. Cardiac injuries caused by blunt trauma：an autopsy based assessment of the injury pattern[J]. Journal of Forensic Sciences，2010，55（1）：82—84.

[3]? ? TEIXEIRA P G R，INABA K，ONCEL D，et al. Blunt cardiac rupture：a 5-year NTDB analysis[J]. The Journal of Trauma：Injury，Infection，and Critical Care，2009，67（4）：788—791.

[4]? ? 肖森，楊濟(jì)匡，肖志，等. 基于正面碰撞實驗的胸部損傷有限元分析[J]. 力學(xué)學(xué)報，2017，49（1）：191—201.

XIAO S，YANG J K，XIAO Z，et al. Analysis of chest injury in frontal impact via finite element modelling based on biomechanical experiment[J]. Chinese Journal of Theoretical and Applied Mechanics，2017，49（1）：191—201. （In Chinese）

[5]? ? COOPER G J，PEARCE B P，STAINER M C，et al. The biomechanical response of the thorax to nonpenetrating impact with particular reference to cardiac injuries[J]. The Journal of Trauma，1982，22（12）：994—1008.

[6]? ? 李沛雨. 碰撞載荷下考慮人體差異的胸腔參數(shù)化建模及損傷研究[D]. 北京：清華大學(xué)，2017.

LI P Y. Parametric modeling and injury analysis of occupant thorax under impact conditions concerning variations among populations[D]. Beijing：Tsinghua University，2017. （In Chinese）

[7]? ? 石向南. 老齡和肥胖乘員有限元模型建模方法及損傷機(jī)理研究[D]. 長沙：湖南大學(xué)，2014.

SHI X N. The research on the modeling method of finite element model and injury mechanisms of older and obese occupant[D]. Changsha：Hunan University，2014. （In Chinese）

[8]? ? GUYTON A C，HALL J E. Textbook of Medical Physiology[M]. Philadelphia：Elsevier Inc，2006：106—107.

[9]? ? PEVEC W C，UDEKWU A O，PEITZMAN A B. Blunt rupture of the myocardium[J]. The Annals Thoracic Surgery，1989，48（1）：139—142.

[10]? LEAVITT B J，MEYER J A，MORTON J R，et al. Survival following nonpenetrating traumatic rupture of cardiac chambers[J]. The Annals of Thoracic Surgery，1987，44（5）：532—535.

[11]? LINK M S，WANG P J，PANDIAN N G，et al. An experimental model of sudden death due to low-energy chest-wall impact （commotio cordis）[J]. New England Journal of Medicine，1998，338（25）：1805—1811.

[12]? 劉維永，蔡建輝，易定華. 心臟撞擊傷生物力學(xué)致傷機(jī)制及心肌挫傷分級[J]. 第四軍醫(yī)大學(xué)學(xué)報，2000，21（5）：540—542.

LIU W Y，CAI J H，YI D H. Experimental study on biomechanical mechanisms of blunt cardiac injury and graded criterion of myo-cardial contution[J]. Journal of the Fourth Military Medical University，2000，21（5）：540—542. （In Chinese）

[13]? STEIN P D，SABBAH H N，VIANO D C，et al. Response of the heart to nonpenetrating cardiac trauma[J].The Journal of Trauma，1982，22（5）：364—373.

[14]? LIFE J S，PINCE B W. Response of the canine heart to thoracic impact during ventricular diastole and systole[J]. Journal of? Biomechanics，1968，1（3）：169—173.

[15]? HUELKE D F，NUSHOLTZ G S，KAIKER P S. Use of quadruped models in thoraco-abdominal biomechanics research[J]. Journal of? Biomechanics，1986，19（12）：969—977.

[16]? RUAN J，EL-JAWAHRI R，CHAI L，et al. Prediction and analysis of human thoracic impact responses and injuries in cadaver impacts using a full human body finite element mode[C]//SAE Technical Paper Series.400 Commonwealth Drive，Warrendale，PA，United States：SAE International，2003.

[17]? IWAMOTO M，NAKAHIRA Y，KIMPARA H.Development and validation of the total human model for safety （THUMS） toward further understanding of occupant injury mechanisms in precrash and during crash[J]. Traffic Injury Prevention，2015，16（1）：36—48.

[18]? DU T Y，CHEN J Q，LAN F C. Development of the Chinese 50th percentile male human body model[C]//IRCOBI Conference Proceedings. Malaga，Spain：IRCOBI，2016：570—572.

[19]? SCHOELL S L，WEAVER A A，URBAN J E，et al. Development and validation of an older occupant finite element model of a mid-sized male for investigation of age-related injury risk[J]. Stapp Car Crash Journal，2015，59（11）：359—383.

[20]? JIANG B H，CAO L B，MAO H J，et al. Development of a 10-year-old paediatric thorax finite element model validated against cardiopulmonary resuscitation data[J]. Computer Methods in Biomechanics and Biomedical Engineering，2014，17（11）：1185—1197.

[21]? 蘭鳳崇，蔡志華，陳吉清，等. 汽車碰撞中胸-腹部的生物力學(xué)響應(yīng)與損傷評價[J]. 華南理工大學(xué)學(xué)報（自然科學(xué)版），2012，40（12）：70—78.

LAN F C，CAI Z H，CHEN J Q，et al. Biomechanical responses and injury evaluation of human thorax and abdomen during vehicle collision[J]. Journal of South China University of Technology （Natural Science Edition），2012，40（12）：70—78. （In Chinese）

[22]? 方紅榮，莊茁. 用三維重建和有限元方法對人體心臟進(jìn)行力學(xué)分析[J]. 醫(yī)用生物力學(xué)，2008，23（1）：43—46.

FANG H R，ZHUANG Z. Biomechanics analysis on human heart by establishment of three-dimension model and finite element method[J]. Journal of Medical Biomechanics，2008，23（1）：43—46. （In Chinese）

[23]? SOMMER G，HASPINGER D C，ANDR？魧 M，et al. Quantification of shear deformations and corresponding stresses in the biaxially tested human myocardium[J]. Annals of Biomedical Engineering，2015，43（10）：2334—2348.

[24]? 王志超. 基于超彈模型的主動脈瓣流固耦合分析[D]. 濟(jì)南：山東大學(xué)，2018.

WANG Z C. Fluid structure interaction simulation of the bioprosthetic heart valve based on hyperelastic constitutive model[D]. Jinan：Shandong University，2018. （In Chinese）

[25]? 陳吉清，劉朝陽，蘭鳳崇. 中國人體上肢碰撞損傷有限元模型的開發(fā)[J]. 華南理工大學(xué)學(xué)報（自然科學(xué)版），2017，45（8）：21—27.

CHEN J Q，LIU C Y，LAN F C. Development of finite element model of impact injury to upper limb of Chinese human body[J]. Journal of South China University of Technology （Natural Science Edition），2017，45（8）：21—27. （In Chinese）

[26]? 黃偉，蘭鳳崇，陳吉清，等. 乘員下肢生物力學(xué)模型與KTH正面碰撞損傷分析[J]. 汽車工程學(xué)報，2016，6（2）：106—114.

HUANG W，LAN F C，CHEN J Q，et al. Occupant lower extremity biomechanical model and KTH injury research in frontal crashes[J]. Chinese Journal of Automotive Engineering，2016，6（2）：106—114. （In Chinese）

[27]? 陳吉清，杜天亞，蘭鳳崇，等. 汽車碰撞中乘員腹部損傷機(jī)理及生物力學(xué)響應(yīng)研究進(jìn)展[J]. 汽車工程學(xué)報，2015，5（2）：79—89.

CHEN J Q，DU T Y，LAN F C，et al. Research progress in injury mechanism and biomechanical response of human abdomen during a car crash[J]. Chinese Journal of Automotive Engineering，2015，5（2）：79—89. （In Chinese）

[28]? 蔡志華. 汽車碰撞中胸部生物力學(xué)響應(yīng)與損傷評估研究[D]. 廣州：華南理工大學(xué)，2013.

CAI Z H. Study on Biomechanical responses and injury assessment of human thorax in vehicle[D]. Guangzhou：South China University of Technology，2013. （In Chinese）

[29]? KROELL C K，SCHNEIDER D C，NAHUM A M. Impact tolerance and response of the human thorax II[C]//SAE Technical Paper Series.400 Commonwealth Drive，Warrendale，PA，United States：SAE International，1974.

[30]? TURK E E，TSANG Y W，CHAMPANERI A，et al. Cardiac injuries in car occupants in fatal motor vehicle collisions--an autopsy-based study[J]. Journal of Forensic & Legal Medicine，2010，17（6）：339—343.

[31]? PATHI V，JONES B，DAVIDSON K G. Mitral valve disruption following blunt trauma：case report and review of the literature[J]. European Journal of Cardio-Thoracic Surgery，1996，10（9）：806—808.

[32]? HUIS M A，CRAFT C A，HOOD R E. Blunt cardiac trauma review[J]. Cardiology Clinics，2018，36（1）：183—191.

[33]? MCDONALD M L，ORSZULAK T A，BANNON M P，et al. Mitral valve injury after blunt chest trauma[J]. The Annals of Thoracic Surgery，1996，61（3）：1024—1029.