王爽 楊寶軍 陳廣新

摘 ?要： 按照股骨在站立位、膝關節(jié)先著地條件下進行分類，分析股骨在不同狀態(tài)下的應力分布情況，并對股骨骨折的臨床應用提供借鑒。利用MIMICS軟件對志愿者股骨CT的DICOM數(shù)據(jù)進行三維重建模型，利用3-matic進行模型修復、光順，利用有限元前處理軟件SimLab進行有限元網(wǎng)格劃分、并應用MIMICS基于股骨CT灰度值賦予材料屬性，最后用ANSYS Workbench結(jié)構靜力學軟件進行股骨應力模擬分析。計算結(jié)果表明股骨在站立位，應力主要集中在股骨的上中段，在膝關節(jié)著地時，股骨下段受力增大，且股踝間受力出現(xiàn)集中現(xiàn)象。計算結(jié)果說明不同受力狀態(tài)下的股骨骨折的可能性。

關鍵詞： DICOM;CT灰度值;股骨;應力分布

中圖分類號： TP319 ? ?文獻標識碼： A ? ?DOI：10.3969/j.issn.1003-6970.2019.04.010

本文著錄格式：王爽，楊寶軍，陳廣新，等. 基于正向工程建模和CT灰度值賦予材料屬性的人體股骨受力分析[J]. 軟件，2019，40（4）：5155

【Abstract】： According to the classification of femur in standing position and knee joint landing first， the stress distribution of femur in different states is analyzed， and the clinical application of femoral fracture is provided for reference. MIMICS software was used to reconstruct the DICOM data of volunteer femoral CT accurately in three dimensions. 3-matic software was used to repair and smooth the model. SimLab software was used to divide the finite element mesh， load the force， and MIMICS was used to assign material attributes based on the gray value of femoral CT. Finally， ANSYS WORKBENCH structural statics software was used to simulate and analyze the femoral stress. When the femur is standing， the stress mainly concentrates on the upper and middle part of the femur. When the knee joint touches the ground， the force on the lower part of the femur increases， and the force between the femur and ankle appears to be concentrated. It shows that different stress states determine the possibility of femoral fracture.

【Key words】： DICOM; CT Gray Value; Femur; Stress distribution

0 ?引言

股骨是人體中最大的長管狀骨，股骨骨折是骨科常見疾病，股骨通過股骨頭與上方的髖臼相連，組成了人體全身最大且最穩(wěn)定的關節(jié)，在股骨頭的外下方是股骨頸，股骨頸的內(nèi)部全部由松質(zhì)骨組 ?成[1]。近年來人股骨骨折的發(fā)病率日益增多，已經(jīng)成為現(xiàn)代骨科常見疾病，股骨骨折方式有兩種，一種是股骨近端骨折，股骨近端骨折又可以分為股骨頭骨折、股骨頸骨折、粗隆間骨折等等幾種;另一種是股骨干骨折[2]。建模與仿真是研究骨骼生物力學的重要手段，在探索骨骼發(fā)病機理和完善治療方案中發(fā)揮重要的作用。一些學者應用逆向工程方法，采用皮質(zhì)骨、松質(zhì)骨分開建模，以各向同性方式分別給皮質(zhì)骨、松質(zhì)骨賦予材料屬性。這種建模和賦予材料屬性方法存在很多缺點，如模型信息丟失較多，效率低下，計算結(jié)果可能會與真實情況相差較多，還有一部分學者采用局部分析法，只分析股骨近端，這種研究不夠全面。本實驗研究基于正向工程方式，利用比利時Materialise公司的MIMICS 21.0軟件對人體股骨DICOM數(shù)據(jù)進行三維重建和基于灰度值方法賦予模型材料屬性，應用解剖學CAD軟件3-matic13.0進行模型表面光順等處理，應用美國Altair公司的SimLab 2017軟件對模型進行有限元分析前處理，應用美國ANSYS公司的ANSYS 18.0進行結(jié)構靜力學仿真計算，通過兩種不同受力狀態(tài)分析，探討股骨骨折與受力的因果關系。

1 ?對象和方法

對象：選取男性健康志愿者1名，42歲，身高170cm，應用CT檢查排除病變、損傷等情況，受試者知情同意并簽署知情同意書，得到醫(yī)院醫(yī)學倫理學委員會批準。

設備與軟件：掃描設備：中國東軟公司的NeuViz 128排螺旋CT;計算機：Dell 7810圖形工作站，16G內(nèi)存，E5處理器，NVIDIA 2G顯卡，21寸液晶顯示器;軟件：比利時Materialise公司醫(yī)學圖像處理軟件MIMICS 21.0、正向工程軟件3-matic 13.0、美國有限元前處理軟件SimLab 2017、美國ANSYS公司有限元仿真分析軟件ANSYS 18.0。

方法：（1）三維重建模型：使用東軟公司NeuViz 128排螺旋CT機薄層掃描志愿者下肢，掃描電壓為120 kv，掃描電流為300 mA，層厚為1 mm。共獲得1206層斷層圖像，以DICOM 3.0格式保存。將獲得1206層斷層圖像以Dicom 3.0格式導入MIMICS21.0軟件中[3]使用高級分割功能CT Bone對股骨斷層圖像進行Bone selection、Mask Threhold、Calculate Part處理，獲得右側(cè)股骨的三維重建模型，從左往右依次為股骨正面觀、側(cè)面觀、后面觀（圖1）。

（2）股骨有限元模型的網(wǎng)格劃分：由于MIMICS獲得股骨初步模型存在表面粗糙、細節(jié)丟失、三角面片不精確等問題，無法直接用于有限元計算，故需要在將MIMICS中獲得的股骨三維初步模型導入3-matic中進行光順（Smooth）、包裹（Wrap）、補洞（Fill hole）等處理，確保模型足夠光滑。應用3-matic Remsh模塊對模型進行面網(wǎng)格、將3-matic畫完的面網(wǎng)格導入SimLab軟件進行體網(wǎng)格劃分。體網(wǎng)格劃分結(jié)果為：四面體10節(jié)點，單元數(shù)201287，節(jié)點數(shù)352691。

（3）股骨有限元模型的材料賦值：將模型網(wǎng)格文件導入MIMICS軟件中，基于CT斷層圖像的灰度值，使用經(jīng)驗公式按十種材料屬性對股骨有限元模型進行材料賦值（圖2）[4]，賦值之后模型以solid187單元類型導出：為避免密度出現(xiàn)負值，本文定義兩種材料類型：

（4）股骨仿真計算：對股骨在站立位、膝關節(jié)著地情況下的結(jié)構靜力學仿真計算，對生理狀態(tài)下的股骨頭垂直加載1000 N的力，約束條件分別為：股骨在站立位情況下：對脛骨關節(jié)面進行約束;膝關節(jié)先著地，在股骨下段的前下方進行約束。

主要觀察指標：不同受力條件下的股骨應力與變形分布情況。

2 ?結(jié)果

2.1 ?股骨在站立情況下的應力與變形分析

由股骨的應力云圖（圖3）可知，股骨在站立情況下，股骨前、中、后的應力分布不均勻，其中最大應力集中在股骨中段附近。股骨頸存在應力集中現(xiàn)象，最大值為25.379，遠小于股骨中段附近的最大值，但股骨頸應力比股骨頭大，轉(zhuǎn)子間也存在應力集中且比大小轉(zhuǎn)子大（圖4）。股骨頭變形最大（圖5），變形程度沿股骨軸向由近端到股骨的中段

部位逐漸減小，股骨中段到股骨遠端幾乎無變形。計算結(jié)果驗證了在臨床上垂直墜落傷時，股骨骨折以股骨中段骨折為主，下段骨折則以螺旋形骨折多見，而股骨頸骨折在這個受傷機制相對較少一些，除了一些骨質(zhì)疏松患者，特別是股骨頸骨質(zhì)疏松為主的，可以出現(xiàn)股骨頸及轉(zhuǎn)子間骨折，頭及大小轉(zhuǎn)子間骨折[5-7]。

2.2 ?股骨在膝關節(jié)著地情況下的應力與變形分析

股骨在膝關節(jié)著地時的應力分布明顯與站力情況不同（圖6）。股骨的中段仍然是受力最集中的區(qū)域，股骨下段應力增大，股骨髁間應力也增大[8-12]，股骨頸的應力也有部分區(qū)域達到了25MPa，轉(zhuǎn)子間的應力也增大不少。由股骨的變形云圖（圖7）可見，股骨變形程度最大的區(qū)域為外側(cè)踝、內(nèi)側(cè)踝、踝間窩、內(nèi)側(cè)踝等區(qū)域。綜上所述，仿真計算結(jié)果與臨床上股骨中段骨折、股骨頸、轉(zhuǎn)子間骨折發(fā)生率相吻合。

3 ?討論

在骨骼生物力學建模與仿真中，采用基于CT的DICOM數(shù)據(jù)通過MIMICS等醫(yī)學圖像處理軟件進行建模方法日趨成為主流[13-17]，但在建模中的有限元模型構建存在很大的差異，目前主流的方法是才用逆向工程軟件比如Geomagic等軟件將人體骨骼模型由STL格式轉(zhuǎn)換成CAD格式，再構建出骨骼模型[18-25]，其中，對皮質(zhì)骨、松質(zhì)骨分別建模，在賦予材質(zhì)上采用各向同性的材料賦值法。此種方法的優(yōu)點是非常有利于有限元模型的裝配、網(wǎng)格劃分、有限元軟件的導入，計算收斂方面也有很多的優(yōu)勢，成功率高。但是，由于人體骨骼的特殊性，使用多個軟件進行往復導入，必然導致模型的信息丟失，而且容易造成模型的結(jié)構與真實情況相差很大。3-matic是比利時Materialise公司開發(fā)的應用于“解剖學CAD”的軟件，可以直接以stl格式進行模型處理、網(wǎng)格劃分，減少逆向工程的往復導入，減少因往復導入導致的信息缺失的錯誤，此外采用MIMICS自有的各向異性的基于灰度值的材料屬性賦予方法，相較于均質(zhì)法，更真實、更科學，更貼近人體真實情況。

參考文獻

[1] 樊瑜波， 王麗珍. 骨肌系統(tǒng)生物力學建模與仿真. 人民衛(wèi)生出版社， 2018. 6： 327-329.

[2] 陳龍， 張軍洋. 基于CT股骨有限元模型精確重建及模態(tài)分析[J]. 計算機仿真， 2014. 31（2）： 280-281.

[3] 陳廣新， 董默， 宋海南， 等. 股骨3D模型構建在股骨頭置換手術中的應用研究[J]. 軟件， 2017. 38（12）： 74-77.

[4] 李峰， 吳華. 基于DICOM數(shù)據(jù)構建股骨三維有限元模型的精確力學分析. 中國組織工程研究[J]， 2013. 17（30）： 5483- 5484.

[5] 方國芳， 林荔軍， 于博， 等. 不同狀態(tài)下的股骨應力分布及臨床應用. 中國組織工程研究[J].

[6] Lubin I M， Aziz N， Babb L J， et al. Principles and Recommendations for Standardizing the Use of the Next- Generation Sequencing Variant File in Clinical Settings[J]. The Journal of Molecular Diagnostics， 2017， 19（3）： 417-426.

[7] Rajapakse C S， Chang G. Micro-Finite Element Analysis of the Proximal Femur on the Basis of High-Resolution Magnetic Resonance Images[J]. Current Osteoporosis Reports， 2018， 16（6）： 657-664.

[8] Mirzaei M， Alavi F， Allaveisi F， et al. Linear and nonlinear analyses of femoral fractures： Computational/experimental study[J]. Journal of Biomechanics， 2018： 155-163.

[9] Khor F， Cronin D S， Watson B， et al. Importance of asymmetry and anisotropy in predicting cortical bone response and fracture using human body model femur in three-point bending and axial rotation[J]. Journal of The Mechanical Behavior of Biomedical Materials， 2018： 213-229.

[10] Dagneaux L， Allal R， Pithioux M， et al. Femoral malrotation from diaphyseal fractures results in changes in patellofemoral alignment and higher patellofemoral stress from a finite element model study[J]. Knee， 2018， 25（5）： 807-813.

[11] Completo A， Fonseca F， Simoes J A， et al. Experimental validation of intact and implanted distal femur finite element models[J]. Journal of Biomechanics， 2007， 40（11）： 2467- 2476.

[12] Hambli R， Lespessailles E， Benhamou C， et al. Integrated remodeling-to-fracture finite element model of human proximal femur behavior[J]. Journal of The Mechanical Behavior of Biomedical Materials， 2013： 89-106.

[13] Synek A， Pahr D H. Plausibility and parameter sensitivity of micro-finite element-based joint load prediction at the proximal femur[J]. Biomechanics and Modeling in Mechanobiology， 2018， 17（3）： 843-852

[14] Chang G， Hotcacho A， Rusinek H， et al. Measurement reproducibility of magnetic resonance imaging-based finite element analysis of proximal femur microarchitecture for in vivo assessment of bone strength[J]. Magnetic Resonance Materials in Physics Biology and Medicine， 2015， 28（4）： 407-412

[15] Klein K F， Hu J， Reed M P， et al. Development and Validation of Statistical Models of Femur Geometry for Use with Parametric Finite Element Models. [J]. Annals of Biomedical Engineering， 2015， 43（10）： 2503-2514.

[16] 王俊然， 杜瑋瑾， 王長江， 郭媛， 陳維毅. 有限元分析不同屈曲狀態(tài)下脛-股關節(jié)的生物力學變化[J]. 中國組織工程研究， 2018， 22（31）： 4975-4981.

[17] 產(chǎn)啟平， 陳龍. 基于等幾何分析的非均質(zhì)股骨近端模型的受力分析[J]. 電子科技， 2018， 31（02）： 15-19.

[18] Grassi L， Vaananen S P， Yavari S A， et al. Experimental Validation Of Finite Element Model For Proximal Composite Femur Using Optical Measurements[J]. Journal of The Mechanical Behavior of Biomedical Materials， 2013： 86-94.

[19] 汪浩， 陳彩鳳， 仲梁維. 股骨有限元模型構建及鈍力性損傷研究[J]. 軟件導刊， 2017， 16（07）： 137-139.

[20] Duarte R J， Ramos A， Completo A， et al. The importance of femur/acetabulum cartilage in the biomechanics of the intact hip： experimental and numerical assessment[J]. Computer Methods in Biomechanics and Biomedical Engineering， 2015， 18（8）： 880-889.

[21] Bryan R， Mohan P S， Hopkins A R， et al. Statistical modelling of the whole human femur incorporating geometric and material properties[J]. Medical Engineering & Physics， 2010， 32（1）： 57-65.

[22] Ennsbray W S， Ariza O， Gilchrist S， et al. Morphology based anisotropic finite element models of the proximal femur validated with experimental data[J]. Medical Engineering & Physics， 2016， 38（11）： 1339-1347.

[23] 常文舉， 丁海. 股骨近端解剖與生物力學研究進展[J]. 醫(yī)用生物力學， 2016， 31（02）： 188-192.

[24] Machado M M， Fernandes P R， Zymbal V， et al. Human Proximal Femur Bone Adaptation to Variations in Hip Geometry[J]. Bone， 2014： 193-199.

[25] Klein K F， Hu J， Reed M P， et al. Validation of a Parametric Finite Element Human Femur Model[J]. Traffic Injury Prevention， 2017， 18（4）： 420-426.