有限元方法在法醫(yī)損傷鑒定中的應(yīng)用

陳憶九，鄒冬華，邵 煜

（司法部司法鑒定科學(xué)技術(shù)研究所 上海市法醫(yī)學(xué)重點(diǎn)實(shí)驗(yàn)室，上海200063）

有限元方法在法醫(yī)損傷鑒定中的應(yīng)用

陳憶九，鄒冬華，邵 煜

（司法部司法鑒定科學(xué)技術(shù)研究所 上海市法醫(yī)學(xué)重點(diǎn)實(shí)驗(yàn)室，上海200063）

有限元方法（finite element method，F(xiàn)EM）是一種與現(xiàn)代計(jì)算機(jī)技術(shù)相聯(lián)系的計(jì)算應(yīng)力分析方法，已逐漸成為生物力學(xué)領(lǐng)域仿真人體結(jié)構(gòu)力學(xué)功能研究的一個(gè)重要實(shí)驗(yàn)手段，可以有效地分析人體結(jié)構(gòu)的物理性質(zhì)，如受外力作用時(shí)人體組織的外部沖擊響應(yīng)及內(nèi)部應(yīng)力分布等，在法醫(yī)損傷學(xué)研究及實(shí)踐鑒定過(guò)程中具有較高的應(yīng)用價(jià)值。文章綜述了利用有限元方法分析高墜、交通及火器事故中顱腦損傷、脊柱損傷、骨盆損傷等的成傷機(jī)制、損傷生物力學(xué)等研究領(lǐng)域的成果，探討其在法醫(yī)損傷學(xué)研究以及法醫(yī)檢案中的應(yīng)用價(jià)值。

法醫(yī)損傷學(xué)；有限元方法；生物力學(xué)；綜述

Abstract：The finite element method is a mathematical technique using modern computer technology for stress analysis，and has gradually become an important experimental tool to simulate human body structures in the biomechanical field.It provides effective analysis of physical properties of human structures，such as the response to external impacts and inner dispersion of stress；and shows great value in the research and practice of forensic trauma examination.This paper reviews the studies of the finite element method on mechanism of craniocerebral，vetebral，and pelvic injuries in fall，traffic accidents and firearm cases，and discusses its values in forensic research and examination.

Key words：forensic traumatology；finite element method；biomechanics；review

有限元方法（finite element method，F(xiàn)EM）是一種求解數(shù)學(xué)物理問(wèn)題的方法，最初應(yīng)用于工程學(xué)，距今已有60多年的發(fā)展歷史。近年來(lái)，隨著計(jì)算機(jī)技術(shù)的發(fā)展和醫(yī)學(xué)影像水平的提高，有限元方法與CT三維重建技術(shù)及其它虛擬現(xiàn)實(shí)技術(shù)相結(jié)合已逐漸應(yīng)用于生物工程學(xué)、醫(yī)學(xué)等多個(gè)學(xué)科領(lǐng)域，并日漸成為生物力學(xué)領(lǐng)域仿真人體結(jié)構(gòu)力學(xué)功能研究的一個(gè)重要實(shí)驗(yàn)手段。目前，隨著專業(yè)的三維有限元分析軟件和模型模擬軟件的研發(fā)和應(yīng)用，三維有限元分析法的過(guò)程得到簡(jiǎn)化，同時(shí)計(jì)算結(jié)果更接近真實(shí)。通過(guò)基于醫(yī)學(xué)圖像通訊標(biāo)準(zhǔn) （digital imaging and communcations in medicine，DICOM）存儲(chǔ)格式的醫(yī)學(xué)圖像建立的人體有限元分析模型，賦予模型物理材料性質(zhì)并合理模擬人體條件，可以有效地分析人體結(jié)構(gòu)的物理性質(zhì)，例如結(jié)構(gòu)的應(yīng)力/應(yīng)變、模態(tài)分析、外部沖擊響應(yīng)等，這為有限元方法在法醫(yī)學(xué)領(lǐng)域運(yùn)用提供了良好的研究前景，對(duì)深入法醫(yī)學(xué)損傷機(jī)制研究具有積極的作用。

1 有限元方法的原理

有限元方法的基本思想是將連續(xù)的求解區(qū)域離散為一組有限個(gè)、按一定方式相互連結(jié)在一起的單元的組合體。由于單元能按不同的連結(jié)方式進(jìn)行組合，且單元本身又可以有不同形狀，因此可以模型化幾何形狀復(fù)雜的求解域。即將需計(jì)算的幾何結(jié)構(gòu)行有限單元?jiǎng)澐?，用網(wǎng)格實(shí)現(xiàn)趨近其真實(shí)結(jié)構(gòu)，用每個(gè)單元內(nèi)假設(shè)的近似函數(shù)來(lái)分片表示全求解域上待求的未知場(chǎng)函數(shù)，通過(guò)對(duì)各個(gè)單元內(nèi)場(chǎng)函數(shù)進(jìn)行計(jì)算，計(jì)算結(jié)果通過(guò)方程組整合后用來(lái)描述整個(gè)結(jié)構(gòu)的情況[1]。其單元尺寸越小，網(wǎng)格密度越高，結(jié)構(gòu)模擬就越真實(shí)，求解計(jì)算的誤差也越小。最小位能原理和基于它的與結(jié)點(diǎn)位移為基本未知量的位移元是固體力學(xué)有限元方法中應(yīng)用最為廣泛、也是最成熟的一種選擇。現(xiàn)行的有限元方法的通用程序幾乎無(wú)例外的都以位移元作為它最主要的甚至唯一的單元形式。通常有限元方法求解都遵循以下基本步驟[2]：（1）物體的離散化。即選擇合適的單元形狀，將擬分析的物體假想分成有限個(gè)分區(qū)或分塊的集合體。單元結(jié)點(diǎn)的位移、節(jié)點(diǎn)力即為所需求解問(wèn)題的基本未知參數(shù)。（2）挑選形函數(shù)或插值函數(shù)。選擇一組函數(shù)，通過(guò)函數(shù)確保每個(gè)“有限單元”節(jié)點(diǎn)位移、節(jié)點(diǎn)力的唯一性，且能通過(guò)節(jié)點(diǎn)位移、節(jié)點(diǎn)力來(lái)確定該單元的位移、受力狀態(tài)。（3）確定單元的性質(zhì)。對(duì)單元的力學(xué)性質(zhì)進(jìn)行描述。一般用單元的剛度矩陣來(lái)描述單元的性質(zhì)，確定單元節(jié)點(diǎn)力與位移的關(guān)系。（4）組成物體的總性質(zhì)方程組。利用最小能原理建立有限元方程。組成物體的總性質(zhì)方程組集成了單元?jiǎng)偠染仃嚭蛦卧刃ЫY(jié)點(diǎn)載荷列陳來(lái)表示整個(gè)物體性質(zhì)的結(jié)構(gòu)剛度矩陣和結(jié)構(gòu)載荷列陣，建立起整個(gè)總節(jié)點(diǎn)載荷與總結(jié)點(diǎn)位移的關(guān)系。（5）解有限元方程和輔助計(jì)算。引入強(qiáng)制邊界條件，解方程得到結(jié)點(diǎn)位移、節(jié)點(diǎn)力。用一定的數(shù)學(xué)方法計(jì)算上述所形成的總性質(zhì)方程組，解出其未知量。然后，根據(jù)實(shí)際問(wèn)題進(jìn)行必要的輔助計(jì)算。

2 有限元方法在法醫(yī)損傷鑒定中的應(yīng)用

決定有限元生物力學(xué)研究精細(xì)度和準(zhǔn)確度的關(guān)鍵在于如何建模，傳統(tǒng)的磨片切片法、三維測(cè)量法、CT圖像處理法建模等，經(jīng)證明存在較多缺點(diǎn)，特別是整個(gè)過(guò)程是人為識(shí)別組織結(jié)構(gòu)，致使誤差在所難免，從而影響模型的精確度。近年來(lái)，MSCT和MRI等高質(zhì)量的掃描圖像使有限元建模的技術(shù)瓶頸得以突破，硬件和軟件方面的革新使MSCT和MRI圖像可方便地進(jìn)行格式轉(zhuǎn)化，產(chǎn)生符合有限元建模所需的基于醫(yī)學(xué)圖像通訊標(biāo)準(zhǔn)存儲(chǔ)格式的醫(yī)學(xué)圖像，目前該方法已被認(rèn)為是三維有限元模型建模的主流方法。隨著研究的深入，該建模方法逐步得到完善、逐漸成熟，再通過(guò)借助一些第三方的軟件，如MIMICS、SIMPLEWARE等，不同學(xué)者建立的不同部位有限元模型[4-8]已日趨完美。利用上述模型，依據(jù)撞擊壓——對(duì)撞壓產(chǎn)生理論（coup-contre coup pressure theory）將有限元模型在PAM-CRASH、MARC、RADIOSS、ANSYS 和 NISA 等軟件中模擬力作用，不僅可真實(shí)地動(dòng)態(tài)再現(xiàn)鈍力性損傷全過(guò)程，還可對(duì)力作用動(dòng)力學(xué)響應(yīng)進(jìn)行分時(shí)分析，并通過(guò)計(jì)算諸如頭部損傷標(biāo)準(zhǔn)（Head Injury Criterion，HIC）[9]等對(duì)力作用部位形變、位移、應(yīng)力等生物力學(xué)響應(yīng)及其程度進(jìn)行研究。因此，有效建立符合真實(shí)損傷情況的有限元模型及開展虛擬狀態(tài)下的應(yīng)力響應(yīng)研究，已成為國(guó)內(nèi)外學(xué)者研究的熱點(diǎn)和難點(diǎn)[4-8，10-32]，也是目前法醫(yī)損傷學(xué)研究的新趨勢(shì)。

2.1 在顱腦損傷中的應(yīng)用

有限元方法在顱腦損傷中的應(yīng)用開展最早亦最為廣泛。早在上世紀(jì)70年代初，國(guó)際上就開始構(gòu)建用于顱腦損傷研究的有限元模型，早期構(gòu)建的人顱腦有限元模型僅涉及顱骨、腦和頸椎等結(jié)構(gòu)，之后逐漸增加了大腦鐮、小腦幕、蛛網(wǎng)膜下腔，腦實(shí)質(zhì)被劃分為大腦半球、小腦及腦干，并區(qū)分出灰、白質(zhì)、腦室，頂顱骨分為外板、板障和內(nèi)板等。目前，有限元方法在法醫(yī)學(xué)上廣泛應(yīng)用于高墜傷、交通傷、火器傷等顱腦損傷研究中，在損傷的發(fā)生預(yù)測(cè)、損傷程度評(píng)估、不同致傷方式造成的損傷辨別方面展現(xiàn)了其高精確度及實(shí)用性[5-8，10-13]。 Raul等[5]將頭部有限元方法應(yīng)用于經(jīng)連續(xù)兩次高墜造成顱腦損傷的案件中，并成功區(qū)分了每一次高墜所致的損傷，該研究對(duì)于沖擊傷、對(duì)沖傷及顱骨骨折的損傷部位及性質(zhì)、程度預(yù)測(cè)均與檢驗(yàn)結(jié)果相吻合；結(jié)果提示該腦損傷有限元模型可以有效地分析法醫(yī)案例中觀測(cè)到的不同腦損傷的成傷機(jī)制。次年，Raul等[6]還運(yùn)用有限元建模方法重現(xiàn)了一起自殺案件的案發(fā)過(guò)程，確認(rèn)了死者自殺導(dǎo)致尸體頭部?jī)商帢審梽?chuàng)發(fā)生的可能性。同時(shí)他們還通過(guò)研究指出有限元分析已被廣泛地運(yùn)用于汽車碰撞評(píng)估試驗(yàn)中，并闡述了現(xiàn)代影像學(xué)技術(shù)和軟件技術(shù)的發(fā)展為有限元方法在法醫(yī)學(xué)領(lǐng)域應(yīng)用提供了良好的研究前景，提出有限元方法可以在法醫(yī)損傷領(lǐng)域逐步開展應(yīng)用并成為常規(guī)檢驗(yàn)方法之一[7]。J.Motherway等[10]運(yùn)用頭部有限元模型模擬外力擊打頭部的情況，其預(yù)測(cè)的骨折部位與實(shí)際解剖相符，體現(xiàn)了有限元方法在鈍力性顱腦損傷預(yù)測(cè)中的高度準(zhǔn)確性與實(shí)用價(jià)值。Magnus Aare等[11]通過(guò)建立頭戴鋼盔的有限元模型模擬了子彈沖擊鋼盔后鋼盔與頭部相互作用造成的損傷，展現(xiàn)了不同的入射角度、鋼盔材料的剛度對(duì)于頭部造成的損傷情況。國(guó)內(nèi)湖南大學(xué)許偉等[12]基于人體解剖學(xué)結(jié)構(gòu)的HBM（Human body model）建立了人體頭部三維有限元模型，采用人體頭部碰撞實(shí)驗(yàn)數(shù)據(jù)，比較了實(shí)驗(yàn)與仿真中頭部的動(dòng)力學(xué)響應(yīng)和顱內(nèi)壓力分布參數(shù)，對(duì)頭部有限元模型進(jìn)行了驗(yàn)證，可以用來(lái)分析研究車輛交通事故中顱腦創(chuàng)傷和損傷機(jī)制。天津科技大學(xué)何培等[13]對(duì)建立的人體顱腦三維有限元模型進(jìn)行碰撞模擬，用剛性圓柱體撞擊眉心處，撞擊速度為6.33m/s，撞擊方式為水平自由撞擊，均取得了相關(guān)的實(shí)驗(yàn)數(shù)據(jù)。

2.2 在脊柱損傷中的應(yīng)用

建立脊柱的有限元模型（包括椎體、椎間盤、韌帶等），施以壓縮、牽張、扭轉(zhuǎn)、剪切等載荷，通過(guò)觀察并計(jì)算脊柱的位移、形變、破壞等情況，可以研究高墜及交通事故所致的諸多脊柱損傷（如揮鞭樣損傷、安全帶損傷、屈曲、爆裂性骨折等）的發(fā)生機(jī)制。Hakim N S等[14]及Yang KH等人[15]通過(guò)對(duì)脊柱的三維有限元分析表明：脊柱承受壓縮載荷時(shí)，以椎體前方或者前下方、終板的中央部位、椎弓根處產(chǎn)生最大主應(yīng)力，損傷易產(chǎn)生于上述部位。劉雷等[16]建立胸腰椎脊柱有限元模型，分別施加垂直壓縮、壓縮屈曲、分離屈曲載荷，成功解釋了壓縮性骨折、爆裂性骨折及安全帶損傷的的力學(xué)機(jī)制。戴立揚(yáng)等[17]對(duì)于揮鞭樣損傷的研究表明，揮鞭樣損傷機(jī)制并不是頸椎屈曲性損傷，而是碰撞時(shí)頸椎的過(guò)伸性損傷，這種過(guò)伸作用下產(chǎn)生的壓縮載荷并不會(huì)導(dǎo)致椎體的壓縮性骨折，而是使椎體產(chǎn)生異常曲線。張建國(guó)等[18]建立頸部有限元模型研究后發(fā)現(xiàn)，在揮鞭樣損傷過(guò)程中，整個(gè)頸椎的最大后伸角度出現(xiàn)在C7-T1關(guān)節(jié)段，其相應(yīng)的椎間盤達(dá)到了最大的馮米塞斯應(yīng)力值，故此部位極易損傷。

2.3 在骨盆損傷中的應(yīng)用

目前，鈍力性作用（如高墜、交通事故等）致骨盆損傷的有限元建模研究報(bào)道較少，骨盆三維有限元模型主要應(yīng)用于臨床分析骨折內(nèi)固定術(shù)后應(yīng)力及應(yīng)變分布。Dalstra等[19]將骨小梁的材料屬性及皮質(zhì)骨的厚度引入有限元，建立了骨盆的三維有限元模型，該模型能有效的描述骨盆皮質(zhì)骨的應(yīng)力分布情況。Anderson等[20]對(duì)個(gè)體化骨盆三維有限元模型進(jìn)行了研究，認(rèn)為不同的骨皮質(zhì)厚度及彈性模量的設(shè)置對(duì)應(yīng)力集中會(huì)產(chǎn)生很大影響。Li等[21]將恥骨間盤與韌帶引入骨盆有限元模型，通過(guò)分析人坐姿時(shí)外力作用于股骨或者大粗隆時(shí)恥骨聯(lián)合的壓力及剪切力的變化而模擬交通事故所致骨盆恥骨聯(lián)合損傷的受傷機(jī)制。通過(guò)模擬落塔式撞擊器對(duì)于骨盆側(cè)方的撞擊，顯示了在骨盆受到側(cè)方?jīng)_擊過(guò)程中的受力分布及損傷情況，結(jié)論為：在對(duì)側(cè)骨盆固定的情況下，側(cè)方?jīng)_擊使恥骨聯(lián)合同時(shí)受到側(cè)方壓縮、向后彎曲、前后方向的以及上下方向的剪應(yīng)力，當(dāng)各力超過(guò)骨盆耐受極限時(shí)即會(huì)在對(duì)側(cè)的恥骨上支發(fā)生骨折；而當(dāng)解除對(duì)側(cè)骨盆的固定后，骨盆受到上述各力程度明顯減小，骨折最易發(fā)生于沖擊側(cè)的恥骨上支，而沖擊帶來(lái)的韌帶的應(yīng)變、恥骨聯(lián)合間的軟骨盤的壓縮及剪切應(yīng)力、恥骨聯(lián)合關(guān)節(jié)的向后彎曲亦對(duì)骨盆損傷有潛在的影響。國(guó)內(nèi)蘇佳燦等人[22-23]利用影像學(xué)數(shù)據(jù)建立了正常的骨盆有限元模型，分別進(jìn)行了沖擊載荷、靜載荷作用下骨盆三維有限元分析及其生物力學(xué)研究，并通過(guò)分析靜載荷及沖擊載荷作用下骨盆各部位應(yīng)力分布以及骨盆各個(gè)單元在應(yīng)力作用下的位移變化等，為臨床分析及判斷骨盆骨折類型、力學(xué)分布、沖擊載荷影響等提供力學(xué)基礎(chǔ)。

2.4 在四肢長(zhǎng)骨損傷中的應(yīng)用

四肢長(zhǎng)骨的有限元構(gòu)建相對(duì)顱腦、脊柱及骨盆建模較為簡(jiǎn)單。國(guó)內(nèi)薛召軍等[24]以肱骨橫截面CT圖像為解剖學(xué)基礎(chǔ)，提取骨端面的幾何形狀，建立肱骨的三維有限元模型；對(duì)于骨折斷面與肱骨軸線成30°、60°、90°的情況，運(yùn)用國(guó)際上通用的大型有限元分析程序ANSYS進(jìn)行計(jì)算，分析肱骨在軸向受壓時(shí)的應(yīng)力分布及大小。蘇佳燦等[25]應(yīng)用三維有限元軟件分析肱骨骨折后應(yīng)力加載方向與加載力值，明確了肱骨斷面與非斷面處應(yīng)力差異較大，斷面處應(yīng)力較集中，且以斷面為中心呈非嚴(yán)格對(duì)稱分布，應(yīng)力較大區(qū)域分布在距骨折部位10mm左右，是其他應(yīng)力區(qū)的兩三倍。樊黎霞等[26]進(jìn)行了人在摔倒時(shí)股骨上端承載能力的有限元分析，對(duì)建立的有限元模型進(jìn)行了驗(yàn)證并對(duì)摔跤姿態(tài)和載荷位置對(duì)股骨承載能力的影響進(jìn)行了深入研究。劉安慶等人[27]應(yīng)用有限元分析軟件ABAQUS建立人體股骨的三維有限元仿真模型，模擬暴力以4000N情況下的應(yīng)力分布，以了解人體股骨在不同情況下的生物力學(xué)特性和應(yīng)力應(yīng)變分布特點(diǎn)。

2.5 在其它組織損傷中的應(yīng)用

除上述各種損傷之外，有限元方法還被運(yùn)用于人體其他體表結(jié)構(gòu)及體內(nèi)器官的損傷評(píng)估中。如Rangarajan等[28]通過(guò)建立有限元模型闡述了嬰兒搖晃綜合癥中視網(wǎng)膜鋸齒緣處出血及視網(wǎng)膜點(diǎn)狀脫落的機(jī)制，并提出了區(qū)分視網(wǎng)膜交通致傷與虐待致傷的方法。Li等[29]建立有限元肋骨模型，模擬對(duì)其施加前后方向的動(dòng)態(tài)彎曲載荷，預(yù)測(cè)的肋骨骨折部位與實(shí)際尸體實(shí)驗(yàn)結(jié)果相符。Richens等[30]通過(guò)建立胸腔有限元模型（包括心臟、主動(dòng)脈及其他相關(guān)胸腔內(nèi)結(jié)構(gòu)）并對(duì)其施加胸部的撞擊，通過(guò)計(jì)算得到了創(chuàng)傷性主動(dòng)脈斷裂最易發(fā)生于主動(dòng)脈峽部及肺動(dòng)脈分支處的結(jié)論。Huang等[31]通過(guò)建立有限元模型研究了上矢狀竇旁橋靜脈進(jìn)入矢狀竇的角度與沖擊方向?qū)蜢o脈所受應(yīng)力的影響，發(fā)現(xiàn)角加速度對(duì)此處橋靜脈產(chǎn)生的應(yīng)力比線加速度產(chǎn)生的更顯著，更易使橋靜脈破裂致急性硬膜下血腫。

3 有限元方法的優(yōu)缺點(diǎn)

有限元方法是根據(jù)變分法原理求解數(shù)學(xué)上可描述的物理問(wèn)題的一種數(shù)值計(jì)算方法，它可以定量地模擬計(jì)算出模型的應(yīng)力變化，屬于一種理論性分析，其結(jié)果只有與其他研究結(jié)果或臨床檢測(cè)結(jié)果相結(jié)合，才能更為真實(shí)的反映標(biāo)本的受力狀況。它的主要優(yōu)點(diǎn)是能夠解決結(jié)構(gòu)、材料性質(zhì)和載荷情況都比較復(fù)雜的問(wèn)題，且通過(guò)有限元方法建立的人體損傷部位有限元模型幾何學(xué)精確度高、生物力學(xué)仿真性強(qiáng)，具有虛擬性、非損傷性、成本低、可重復(fù)操作、實(shí)驗(yàn)條件（單元、節(jié)點(diǎn)、材料特性、載荷、邊界條件）可控制等特點(diǎn)，其在精確模擬人體組織的生物材料力學(xué)特性等方面是臨床檢查、尸體解剖、假人模型等研究方法無(wú)法比擬的。但受有限元建模的個(gè)體化差異及邊界條件（有限元模型的邊界、人體組織與外界的接觸面、人體不同組織間的接觸面）等影響，直接關(guān)系到不同組織結(jié)構(gòu)間相對(duì)運(yùn)動(dòng)，從而影響到生物力學(xué)模擬過(guò)程最終的評(píng)估結(jié)果。另外由于人體某些組織結(jié)構(gòu)的損傷閾值數(shù)據(jù)尚未測(cè)得，因此對(duì)于這些結(jié)構(gòu)的有限元分析結(jié)果只能顯示其對(duì)于外力響應(yīng)的生物力學(xué)指標(biāo)的大小、高低，并預(yù)測(cè)最易損傷的部位，而無(wú)法評(píng)估外力屬性與損傷的關(guān)聯(lián)性。

4 展望

隨著MSCT和MRI等高質(zhì)量的掃描圖像數(shù)據(jù)轉(zhuǎn)化為符合有限元建模所需的基于醫(yī)學(xué)圖像通訊標(biāo)準(zhǔn)存儲(chǔ)格式的醫(yī)學(xué)圖像這一技術(shù)瓶頸的突破，以此建立起來(lái)的符合真實(shí)情況的損傷三維有限元模型已成為可能，在損傷應(yīng)力分析方面已顯示了其優(yōu)勢(shì)和價(jià)值。通過(guò)對(duì)人體組織不同幾何特征、材料屬性和邊界條件進(jìn)行賦值和約束，建立用于損傷生物力學(xué)分析的有限元人體模型，可以對(duì)各種作用力加載人體后產(chǎn)生的力學(xué)變化和形變特征進(jìn)行仿真研究，實(shí)時(shí)分析計(jì)算其應(yīng)力、應(yīng)變、位移、速度、加速度等應(yīng)力響應(yīng)參數(shù)變化，在法醫(yī)實(shí)踐過(guò)程中分析不同致傷方式、不同應(yīng)力條件下的損傷形態(tài)學(xué)、損傷生物力學(xué)變化等具有很好的應(yīng)用前景。

但將有限元方法應(yīng)用于法醫(yī)損傷學(xué)研究及鑒定實(shí)踐，涉及法醫(yī)學(xué)、診斷影像學(xué)、磁共振波譜學(xué)、計(jì)算機(jī)仿真學(xué)、三維圖像學(xué)以及生物力學(xué)等多學(xué)科的交叉融合，現(xiàn)階段國(guó)內(nèi)法醫(yī)學(xué)領(lǐng)域囿于學(xué)科專業(yè)的限制，對(duì)于有限元方法在科研及檢案實(shí)踐中的應(yīng)用較少。相信，隨著學(xué)科的融合和交叉研究的深入，充分利用有限元方法動(dòng)態(tài)地分析損傷應(yīng)力分布及定量分析應(yīng)力集中等損傷生物力學(xué)研究，其研究結(jié)果應(yīng)用于鑒定實(shí)踐必將推動(dòng)司法鑒定行業(yè)的技術(shù)進(jìn)步。

[1]謝江波，劉亞青，張鵬飛，等.有限元方法概述[J].現(xiàn)代制造技術(shù)與裝備，2007，（5）：29-30.

[2]成思源.有限元方法的方法論[J].重慶大學(xué)學(xué)報(bào)（社會(huì)科學(xué)版），2001，7（4）：61-63.

[3]王元豐.有限元方法的方法論初探[J].科學(xué)技術(shù)與辯證法，1990（6）：16-20.

[4]A.Mota，W.S.Klug，M.Ortiz，et al.Finite-element simulation of firearm injury to the human cranium[J].Computational Mechanics，31（2003）：115-121.

[5]Raul J S，Baumgartner D，Willinger R，et al.Finite element modelling of human head injuries caused by a fall[J].Int J Legal Med，2006，120（4）：212-218.

[6]Raul JS，Deck C，Meyer F，et al.A finite element model investigation of gunshot injury[J].Int J Legal Med，2007，121（2）：143-146.

[7]Raul J S，Baumgartner D，Willinger R，et al.Finite-element models of the human head and their applications in forensic practice[J].Int J Legal Med，2008，122（4）：359-366.

[8]朱偉民，付雙林，賈昭，等.三維顱骨模型的建立和有限元及實(shí)驗(yàn)應(yīng)力分析的研究[J].生物醫(yī)學(xué)工程研究，2007，（3）：244-247.

[9]FMVSS208.Occupant Crash Protections[S].

[10]Motherway J，Doorly M C，Curtis M，et al.Head impact biomechanics simulations：A forensic tool for reconstructing head injury?[J].Legal Medicine，2009，11（S1）：S220-S222.

[11]Magnus Aare，Svein Kleiven.Evaluation of head response to ballistic helmet impacts using the finite element method[J].Int J Impact Eng，2007，34（3）：596-608.

[12]許偉，楊濟(jì)匡.研究顱腦交通傷的有限元模型的建立及驗(yàn)證[J].生物醫(yī)學(xué)工程學(xué)雜志，2008，25（3）：556-561.

[13]何培，張建國(guó)，阮世捷，等.基于螺旋CT的人體顱腦三維有限元模型構(gòu)建[J].天津科技大學(xué)學(xué)報(bào)，2006，21（2）：36-38.

[14]Hakim N S，King A I.A three dimensional finite element dynamic response analysis of a vertebra with experimental verification[J].Journal of Biomechanics.1979，12（4）：277-285.

[15]Yang KH，et al.Finite element model of a functional spinal unit[M].American society of mechanical engineers biomechanics symposium 1983：137-140.

[16]劉雷，沈根標(biāo).胸腰椎損傷的三維有限元模型的應(yīng)力分析[J].中華創(chuàng)傷雜志，1995，11（6）：343-344.

[17]戴力揚(yáng).揮鞭樣損傷的生物力學(xué)[J].醫(yī)用生物力學(xué)，2002，17（4）：248-252.

[18]張建國(guó)，周蕊，薛強(qiáng).基于揮鞭樣損傷研究的頸部有限元模型的建立及驗(yàn)證[J].中國(guó)生物醫(yī)學(xué)工程學(xué)報(bào).2008，27（3）：389-392.

[19]Dalstra M，Huiskes R，van Erning L.Development and validation of a three-dimensional finite element model of the pelvic bone[J].J Biomech Eng 1995；117（3）：272-278.

[20]Anderson AE，Peters CL，Tuttle BD，et al.Subject-specific finite element model of the pelvis：development，validation and sensitivity studies.J Biomech Eng，2005，127（3）：364-373.

[21]Li Z，Kim JE，Davidson JS，et al.Biomechanical response of the pubic symphysis in lateral pelvic impacts：a finite element study[J].J Biomech，2007，40（12）：2758-2766.

[22]Li Z，Kim J E，Davidson J S，et al.Biomechanical response of the pubic symphysis in lateral pelvic impacts：A finite element study[J].Journal of biomechanics.2007，40（12）：2758-2766.

[23]蘇佳燦，管華鵬，張春才，等.沖擊載荷作用下骨盆三維有限元分析及其生物力學(xué)意義.中國(guó)骨傷，2007，20（7）：455-457.

[24]薛召軍，吳建國(guó)，丁祖泉，等.肱骨骨折重建及三維有限元分析[J].醫(yī)用生物力學(xué)，2002，17（1）：24-27.

[25]蘇佳燦，張春才，薛召軍，等.應(yīng)用三維有限元分析肱骨骨折后應(yīng)力加載方向與加載力值[J].中國(guó)臨床康復(fù)，2005，9（14）：234-235.

[26]樊黎霞，Eric Wang.人在摔倒時(shí)股骨上端承載能力的有限元分析——（II）摔跤姿態(tài)和載荷位置對(duì)股骨承載能力的影響研究[J].生物醫(yī)學(xué)工程學(xué)雜志，2006，23（6）：1245-1249.[27]劉安慶，張銀光，王春生，等.人股骨生物力學(xué)特性的三維有限元分析[J].西安醫(yī)科大學(xué)學(xué)報(bào)，2001，22（3）：242-244.

[28]Rangarajan N，Kamalakkannan S B，Hasija V，et al.Finite element model of ocular injury in abusive head trauma[J].Journal of American Association for Pediatric Ophthalmology and Strabismus.2009：1-6.

[29]Li Z，Kindig M W，Kerrigan J R，et al.Rib fractures under anterior-posterior dynamic loads：Experimental and finite-element study[J].Journal of Biomechanics.2010，43：228-234.

[30]Richens D，F(xiàn)ield M，Hashim S，et al.A finite element model of blunt traumatic aortic rupture[J].European journal of cardio-thoracic surgery.2004，25（6）：1039-1047.

[31]Huang H M，Lee M C，Chiu W T，et al.Three-dimensional finite element analysis of subdural hematoma[J].The Journal of trauma.1999，47（3）：538-544.

（本文編輯：張建華）

Application of the Finite Element Method in Forensic Traumatology

CHEN Yi-jiu，ZOU Dong-hua，SHAO Yu
（Shanghai Key Laboratory of Forensic Medicine，Institute of Forensic Sciences，Ministry of Justice，Shanghai 200063，China）

DF795.4

A

10.3969/j.issn.1671-2072.2010.04.005

1671-2072-（2010）04-0028-05

2010-04-20

國(guó)家自然科學(xué)基金資助項(xiàng)目（30872920、81072507）；上海市自然科學(xué)基金資助項(xiàng)目（10ZR1431200）；科技部科研院所社會(huì)公益研究專項(xiàng)（GY0906）。

陳憶九（1961-），男，研究員，碩士研究生導(dǎo)師，主要從事法醫(yī)病理學(xué)研究。E-mail：yijiuchen@yahoo.com.cn。