基于CT的正常與骨質疏松股骨模態分析

王汝良 張凱旋 胡霖霖

摘 ?要： 通過對志愿者人體股骨的CT影像數據進行圖像處理，獲得高精度的股骨三維模型，采用正向工程建立有限元模型，針對股骨的各向異性的材料特性，應用股骨材料屬性經驗公式 獲得正常股骨的材料屬性與骨質疏松的股骨材料和屬性，利用ANSYS MODAL仿真軟件對正常股骨、骨質疏松的股骨進行前6階模態計算，獲得前6階的股骨振動頻率。根據所得數據，可對股骨假體置換手術提供借鑒，應避免與股骨的頻率相近而產生諧振。

關鍵詞： 股骨;模態;正向工程;振動頻率

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

本文著錄格式：王汝良，張凱旋，胡霖霖，等. 基于CT的正常與骨質疏松股骨模態分析[J]. 軟件，2019，40（4）：7780

【Abstract】： Through image processing of CT image data of human femur of volunteers， a high-precision three- dimensional model of femur is obtained. A finite element model is established by forward engineering. According to the anisotropic material properties of femur， the material properties of normal femur and the femoral material and properties of osteoporosis are obtained by using empirical formula of femoral material properties. The first 6 modes of normal femur and osteoporotic femur are calculated by using ANSYS MODAL simulation software to obtain the first 6 modes of femur vibration frequency. According to the obtained data， it can provide reference for femoral prosthesis replacement surgery and should avoid resonance caused by the frequency close to that of femur.

【Key words】： Femur; Modal; Forward engineering; Oscillation frequency

0 ?引言

髖關節和膝關節是人體的重要承重部位，也是關節損傷和炎癥的多發部位，假體置換是緩解人體關節嚴重失效的常用方法之一[1-3]。由于置換假體涉及到非常復雜的生物力學，置換的假體會與髖關節生理載荷發生響應，造成一般的實驗手段無法準確的測定。隨著計算機硬件和有限元理論的發展，利用計算機模擬來建立更接近真實形態股骨生物力學模型，對于置換假體的優化設計提供了很重要的參考[4]。目前，在骨骼建模中，多數學者應用CAD逆向建模技術三維重建出模型，但是由于 CAD 軟件曲面擬合效果帶來很多的誤差，不同軟件之間的切換、導入，容易導致模型的信息丟失，而且建模費時，費力，因此建模結果的準確性會影響后續的計算結果。在材料方面，多數學者采用皮質骨、松質骨分開建模、分別賦予均質、各向同性的材料屬性，甚至有的采用一種材料屬性，這種方法與真實的股骨材料特性差距較大。本文采用的是0.5 mm層厚的股骨CT圖像DICOM數據進行三維建模，利用MIMICS的高級分割功能、計算三維工具，提取出股骨的三維模型，并通過解剖學CAD軟件3-matic對模型進行光順等編輯處理，利用SimLab強大的網格劃分功能進行網格劃分，材料方面，應用基于CT灰度值技術賦予股骨10種各向異性的非均質股骨材料屬性[5，6]，以模擬股骨的實際力學特性。

1 ?材料與方法

1.1 ?基于CT的股骨DICOM數據資料

數據來源：經志愿者本人同意，并簽署知情同意書，采集男性、年齡41周歲的志愿者股骨CT影像DICOM格式數據。掃描設備為東軟公司NMS/ NeuViz 128排CT機，矩陣：512×512;電壓：120 kv;層距：0.5 mm;層數：1206。本實驗經醫院倫理委員會批準，符合倫理學規定;圖像后處理工作站：戴爾Precision T7810：Xeon E5-2609 v3 中央處理器、nVIDIA Quadro k2200圖像顯卡;16G內存;圖像處理及分析軟件：Mimics 21.0、3-matic 13.0（比利時Materialise公司醫學圖像處理軟件和解剖學CAD軟件）;SimLab 2017（美國Altair公司有限元分析前處理軟件）;Ansys Workbench 18.0（美國ANSYS公司大型有限元仿真軟件包）。

1.2 ?方法

本實驗采用正向工程方式建立有限元模型，采用各向異性、非均質的方式賦予材料屬性。具體流程為（圖1）。

1.2.1 ?股骨模型三維重建與網格劃分

利用MIMICS 21.0軟件導入股骨CT影像DICOM數據，使用CT Bone Segmentation功能，將所選擇的股骨提取出來，并初步進行光滑等處理，直接導入3-matic13.0軟件中，對其進行補洞、光滑、修復三角面片缺陷以及初步劃分面網格等。經上述操作后，將股骨模型導入有限元分析前處理軟件SimLab 2017中，進行網格劃分。SimLab軟件網格劃分功能強大，是一個面向流程，基于幾何特征的有限元建模軟件。應用SimLab軟件劃分股骨stl格式的模型，省略了其他逆向工程的有限元建模過程，即采用逆向工程軟件將stl格式的文件通過Geomagic等逆向工程軟件轉換成以點、線、面為結構的CAD軟件格式（如.step、parasolid）。相較于逆向工程的建模方法，本建模模式可有效的保證了模型的結構完整性、更接近于真實的模型結構、減少不同軟件切換而導致的幾何信息丟失問題。使用SimLab 2017劃分股骨四面體網格節點數為390146、單元數為266073，單元類型為Solid187。

1.2.2 ?股骨材料賦予

1）正常股骨的材料賦予

2）骨質疏松股骨的材料賦予

以正常的股骨的彈性模量減少66%，泊松比保持不變為條件賦予股骨材質，模擬形成骨質疏松股骨有限元模型[8-10]。

1.2.3 ?股骨的振動模態計算

本實驗將股骨視為無阻尼的振動系統，同時按照自然股骨的生理狀態，約束股骨頭部所有節點的 3個方向的平動自由度，保留其3個方向的轉動自由度（髖關節位置）;約束股骨髁間窩處節點的3個平動自由度和繞y、z軸的轉動自由度，僅保留其繞 x軸的轉動自由度（膝關節位置）。分別對骨質疏松股骨和正常股骨進行6階模態計算。

2 ?結果

將股骨的模型定位在空間坐標OXYZ的坐標系中，XOY面為股骨的橫斷面，XOZ 面為矢狀面，YOZ面為冠狀面。由圖2、圖3可見，股骨的6階陣型主要是彎曲和扭轉，分布在股骨的第1階振型（矢狀面內）與第2階振型（冠狀面內）都是1階

彎曲;第3階振型（矢狀面內）與第4階振型（冠狀面內）都是2階彎曲;第5階振型為沿股骨長軸的2階扭轉;第6階振型為空間的2階扭轉。從基頻開始股骨振型趨于復雜，彎曲、扭轉的復合使得股骨產生微動位移。股骨的前兩階振型中，股骨的位移最大部位為股骨中段、上段，股骨踝也出現較大位移。從表1的正常股骨6階頻率與骨質疏松股骨6階頻率對比發現，正常的股骨振動頻率都要比骨質疏松股骨的頻率高。

3 ?討論

本文通過對股骨CT影像DICOM數據進行三維重建，獲得股骨的三維模型，利用3-matic軟件的模型修復功能，對模型進行三角面片優化、光順，利用SimLab軟件對模型進行四面體網格劃分，獲得質量較高的有限元模型，應用彈性模量差異性模擬正常股骨和骨質疏松股骨的材料屬性，最后應用ANSYS MODAL分別對股骨進行模態振動仿真計算，提取出股骨的固有6階振型頻率。通過基于灰度值賦予材料屬性、基于正向工程構建有限元模型，可以獲得較為準確的材料分布特性[11-14]，減少因不同軟件的文件導入造成模型信息缺失而降低準確準確性，又可以保證高度不規則模型的網格劃分質量，進而有利于提高模態計算的精度[15-17]。股骨的振動特性為股骨的假體優選提供了指導，在使用置換假體的過程中應避開股骨的固有頻率。

參考文獻

[1] Pyburn E， Goswami T. Finite element analysis of femoral components paper III-hip joints[J]. Materials & Design， 2004， 25（8）： 705-713.

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

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

[4] 韓樹洋， 葛世榮. 人工全髖關節置換術對天然股骨生物力學行為的影響[J]. 醫用生物力學， 2010， 12（6）： 471-474.

[5] 張國棟， 廖維靖， 陶圣祥， 等.股骨有限元分析賦材料屬性的方法[J]. 中國組織工程與臨床康復研究. 2009， 13（43）： 8436-8437.

[6] 董鵬飛， 雷建銀， 劉海波， 等. 基于CT圖像的股骨上段有限元建模及單元尺寸分析[J]. 醫用生物力學. 2016， 31（2）： 129-130.

[7] 方國芳， 林荔軍， 于博， 等. 不同狀態下股骨的應力分布及臨床應用[J]. 中國組織工程與臨床康復研究. 2012， 16（17）： 3045-3047.

[8] 郭文文， 劉靜， 曹慧， 等. 正常與骨質疏松肱骨的三維重建及有限元分析[J]. 中國醫療設備. 2018， 33（4）： 34-36.

[9] 汪金平， 楊天府. 股骨生物力學特性的有限元分析[J]. 中華創傷骨科雜志， 2005， 7（10）： 931-934.

[10] 徐振煬， 胡永成， 葉金鐸， 等. 股骨三維有限元模型的建模方法[J]. 天津理工大學學報， 2015， 10（31）： 52-54.

[11] 陳國棟， 羅羽婕， 王銳英. 有限元分析在股骨生物力學研究中的應用[J]. 實用醫學雜志， 2011， 27（2）： 334-336.

[12] 展影， 雷新瑋， 屈瑾， 等. 不同股骨假體植入角度的髖關節表面置換術后股骨近端力學的有限元研究[J]. 國際醫學放射學雜志， 2010， 33（5）： 451-455.

[13] 呂志軍， 李玲. 不同材料嵌體修復的三維有限元分析[J].齊魯醫學雜志， 2010， 25（5）： 445-446.

[14] 成海平， 柳松楊， 王興偉， 等. 人股骨上段三維有限元模型的建立[J]. 生物醫學工程研究， 2010， 29（2）： 106-108.

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

[16] 崔林江， 張峰， 修樹毅. 三種內固定治療方法在股骨粗隆間骨折的應用[J]. 青島大學醫學院學報， 2007， 43（5）： 457-458.

[17] 雷建銀， 劉海波， 王志華， 等. 站立位下骨盆與骨折內固定穩定性分析[J]. 醫用生物力學， 2014， 29（6）： 517-523.