基于定量CT的脊柱有限元分析現狀

湯凱文，綜述，熊波，劉社庭審校

（南華大學附屬郴州市第一人民醫院脊柱外科，湖南郴州423000）

基于定量CT的脊柱有限元分析現狀

湯凱文，綜述，熊波，劉社庭審校

（南華大學附屬郴州市第一人民醫院脊柱外科，湖南郴州423000）

脊柱；體層攝影術，X線計算機；有限元分析；三維定量計算機斷層掃描；椎體

骨質疏松容易誘發脊柱胸、腰段骨折，相關研究人員花費了大量時間及精力研發能監測人體脊柱中骨礦物質含量的影像設備。二維雙能X射線骨密度儀（dualenergy X-ray absorptiometry，DXA）是應用于骨質疏松癥的診斷及監測骨礦物質含量最常用的工具[1-2]。但因為標準DXA應用的是二維影像技術，其不能分析骨質的三維空間結構和區分在骨皮質與骨小梁之間礦物質密度變化，三維定量計算機斷層掃描（quantitative computed tomography，QCT）提供了骨質結構的三維影像及骨礦物質的空間分布情況，因此，也被應用于分析骨皮質與骨小梁各自的礦物質密度[3]。一種源于在機械及土木工程領域的計算機技術——有限元分析法能采用從QCT中獲得的骨質圖像創建使用者所感興趣的脊柱有限元模型，然后，進行有限元分析，評估在不同影響因素情況下骨骼的強度變化[4-6]。

有限元法需要將一個整體結構分割成許多小的部分，而將小的部分稱之為單元。當外界的應力加載于物體結構上時通過有限元法能分析施加在該結構上每個單元的機械應力及物體的形變。數字QCT影像被劃分為許多小體，這些小體均具有所知的骨密度。因為在之前的實驗研究中已有學者發現了在基于QCT的骨礦物質密度與骨質的機械性質之間的關系，在這樣的情形下，骨骼的三維醫學影像轉化為三維有限元模型是非常便利的?；赒CT影像的有限元分析作為一種適用于生物力學分析的工具，可用來評估骨骼的機械載荷。本文主要概述了采用脊柱的QCT圖像創建脊柱有限元模型并對其進行有限元分析，通過該法能確定一些因素對脊柱生物力學的影響，如年齡、性別、骨質疏松、抗骨質疏松藥物治療等，攻克了許多臨床上的難題。

1 有限元模型創建及分析的基本步驟

對脊柱進行有限元分析需要6個步驟：（1）圖像分割；（2）網格劃分；（3）骨礦物質的賦值；（4）邊界條件的約束和加載；（5）求解；（6）后處理。每一個步驟對最后的結果均會產生重要的影響，而且對于每一個需要用有限元分析解決的臨床問題，每一個步驟均是特定的。

1.1 圖像分割通過獲得所需要的解剖部位的CT影像，從中提取骨骼的幾何信息及骨密度，創建基于QCT的骨骼幾何模型，該過程稱為圖像分割。目前的研究正致力于使圖像分割技術盡可能自動化，其能幫助減少手動操作產生的誤差及縮短建立有限元模型的時間。通過利用解剖結構的信息進行圖像分割是當前最常用的方法之一[7-9]。

1.2 網格劃分在創建了需要分析的骨骼幾何模型后需要劃分為更小的有限元單元，再將其單元組合，該過程稱為網格劃分。模型經過網格劃分后包含了許多不同形狀的單元和節點。四面體單元由4個角和4個面組成，是最常用的單元之一。四面體單元的優點在于能提供更光滑的骨皮質表面，有利于網格劃分，而且使有限元分析結果更加精確[10]。

1.3 骨礦物質的賦值每個單元在劃分網格時均必須被賦予材料屬性，當載荷施加在這些單元時就不會被壓縮變形。通過賦予模型基于QCT的骨礦物質濃度的材質并且進行力學測試，結果顯示，與實體生物力學測試比較，在誤差范圍內。所以，基于QCT的骨礦物質濃度數據可直接定義有限元模型的材料屬性。

1.4 邊界條件的約束和加載當模型網格劃分及材料屬性被定義后，接下來需要對模型賦予正確的邊界條件及機械載荷。邊界條件的約束需要施加在有限元模型的節點上。機械載荷包括面載荷及集中載荷，均能直接施加在節點及單元上。當邊界條件及載荷一起施加于有限元模型時能模擬生理情況下力的分布形式。對脊柱進行有限元分析時最大難點是定義最適合的邊界條件和機械載荷，以便得到最接近生理的生物力學分析。近年來，有一個權威的OrthoLoad數據庫網站能提供施加在脊柱、骨盆等部位力量的大小及方向的數據，且經常被國內外學者引用[11-12]。

1.5 求解完成有限元模型的建立、材料的賦值、邊界條件的約束和載荷的施加步驟后，再對有限元模型進行求解，其能計算出每個節點的位移和應力。對于線性模型，求解只需一步就能完成，而且相對于非線性模型，線性模型的計算速度快且處理技術簡便。

1.6 后處理在建立有限元模型并求取結果后便進入后處理階段，其是有限元分析的重要步驟，能將節點自由度的變化參數轉換為分析需要的數據，這些數據也可通過多種形式的圖形直觀地顯示出來。

基于QCT建立的有限元模型進行有限元分析，6個步驟中的每一步對最后的結果均有重要影響。對復雜的臨床問題采取有限元分析法時依靠的是每一個步驟的正確實施。通過機械測試對有限元模型進行驗證能證明模型建立的有效性和有限元分析結果的可靠性[13-16]。脊柱有限元模型建立、模型材料的賦值、邊界條件的約束、接近生理狀態的加載、模型后處理及最后對模型的有效性驗證等能使脊柱外科醫生在脊柱的生物力學領域中有更加深刻的見解。

2 脊柱的有限元分析

因骨質疏松引發的骨折最常見的是胸、腰椎骨折，所以，研究骨質的流失及藥物治療對椎體強度的生物力學影響很有必要。1990年初就有學者利用QCT對脊柱進行了有限元分析，對椎體強度進行了有效評估。隨著CT影像技術和計算機技術的發展，促使基于QCT的有限元分析法應用于更加廣泛的脊柱研究中。這些研究有助于闡明椎體結構對椎體強度和骨折風險的影響，老年人中男、女性骨質流失的差異，不同抗骨質疏松藥物的療效等。

2.1 有限元分析在評估骨折風險中的應用骨的微結構及骨的力學特性是決定骨折危險性的重要因素[17]，在對男性胸、腰椎骨折的大量研究中發現，基于QCT的胸、腰椎有限元分析對評估骨折風險具有重要價值[18]。L1椎體的強度通過有限元分析顯示，骨折風險比率為8.5，甚至在校正年齡、種族、體質量指數、骨面密度后。近年來，有研究在服用糖皮質激素誘發骨質疏松的男性群體中經基于脊柱QCT的有限元分析發現，因骨質疏松導致的T12椎體骨折中，T12椎體的抗壓強度較沒有骨折椎體的抗壓強度小29.0%[19]。分析骨微結構、骨密度、骨力學特性可進一步評估骨強度[20-21]，基于QCT的有限元分析是一個預估骨折風險的工具。

2.2 老年人中男、女性骨質流失的差異隨著年齡增大所致的骨骼量流失，對女性脊柱方面引起的生物力學影響可能大于男性。Christiansen等[22]應用基于T10及L3椎體QCT的有限元分析發現，在兩性群體中椎體骨小梁丟失量是類似的，但女性骨皮質層骨量丟失更多。在平均年齡78歲的女性老年人中L3椎體骨皮質強度較平均年齡40歲的中年人低37.0%，T10椎體骨皮質強度低33.0%。但38～78歲男性群體椎體強度比較，無明顯差別。表明在女性群體中椎體皮質層骨量丟失可能是導致骨折發生率增加的一個因素。

2.3 不同抗骨質疏松藥物療效的對比基于脊柱的QCT有限元分析法能評估骨吸收抑制劑（抗骨質疏松藥物）的療效。女性骨質疏松癥患者每天口服阿倫磷酸鈉10 mg，6個月后對患者L3椎體進行有限元分析，通過6個月前后對比發現，服藥后患者椎體強度較用藥前增加了4.9%[23]。有一項類似的研究對骨質疏松癥患者服用抗骨質疏松藥物——伊班膦酸，每次口服150 mg，每個月1次，服用1年。1年后對患者采用有限元分析法發現，患者椎體強度也較治療前增加[24]。在該項研究中用伊班膦酸鹽治療的女性骨質疏松癥患者，腰椎強度增加了4.0%，但服用安慰劑者椎體強度反而降低了3.0%。

特立帕肽是一種重組人甲狀旁腺激素1-34，具有增加骨量、骨形成、再吸收等促合成效應。評估該抗骨質疏松藥物對骨質的生物力學影響，采用有限元方法對患者腰椎進行分析，結果顯示，特立帕肽增強骨質強度的效果強于骨吸收抑制劑。如女性骨質疏松癥患者每天服用特立帕肽20 μg，在服用藥物后相應階段對患者腰椎進行有限元分析顯示，L3椎體在服用藥物6個月后強度增加13.0%，而在18個月后增強了21.0%。

應用脊柱有限元分析法可幫助發現隨著年齡的增大脊柱強度在男、女性之間的差別，以及骨吸收抑制劑和重組人甲狀旁腺激素1-34治療骨質疏松患者的療效對比。然而，對于還有很多影響骨質代謝和骨質量的因素，如慢性腎臟疾病、糖尿病、放療等關于其對脊柱生物力學的影響所了解的知識還很局限。

3 小結與展望

在過去的20年中基于脊柱的QCT有限元分析法是一種很有價值的生物力學研究工具，其首先通過QCT獲取感興趣的影像，然后，再通過相關軟件建立有限元模型，賦予相對應的材質，模擬所需要的力學狀態，進行有限元分析，其能評估患者因年齡增大、抗骨質疏松藥物的治療所導致骨質強度的變化情況，以此進行對骨質疏松癥的治療及預防。隨著利用基于脊柱QCT的有限元分析法在研究脊柱相關問題方面日漸增多，越來越多關于脊柱病理和臨床的問題均會迎刃而解，也會推動脊柱外科在生物力學領域的發展。

[1]Adams JE.Advances in bone imaging for osteoporosis[J].Nat Rev Endocrinol，2012，9（1）：28-42.

[2]Blake GM，Fogelman I.The clinical role of dual energy X-ray absorptiometry[J].Eur J Radiol，2009，71（3）：406-414.

[3]Lang TF.Quantitative computed tomography[J].Radiol Clin North Am，2010，48（3）：589-600.

[4]Amin S，Kopperdhal DL，Melton LJ，et al.Association of hip strength estimates by finite-element analysis with fractures in women and men[J].J Bone Miner Res，2011，26（7）：1593-1600.

[5]Shim VB，Pitto RP，Anderson IA.Quantitative CT with finite element analysis：towards a predictive tool for bone remodelling around an uncemented tapered stem[J].Int Orthop，2012，36（7）：1363-1369.

[6]Zysset PK，Dall′ara E，Varga P，et al.Finite element analysis for prediction of bone strength[J].Bonekey Rep，2013，2：386.

[7]Cerveri P，Manzotti A，Marchente M，et al.Mean-shifted surface curvature algorithm for automatic bone shape segmentation in orthopedic surgery planning：a sensitivity analysis[J].Comput Aided Surg，2012，17（3）：128-141.

[8]Wu C，Murtha PE，Jaramaz B.Femur statistical atlas construction based on two-level 3D non-rigid registration[J].Comput Aided Surg，2009，14（4/5/6）：83-99.

[9]Dai Y，Niebur GL.A semi-automated method for hexahedral mesh construction of human vertebrae from CT scans[J].Comput Methods Biomech Biomed Engin，2009，12（5）：599-606.

[10]Ramos A，Simoes JA.Tetrahedral vs hexahedral finite elements in numerical modelling of the proximal femur[J].Med Eng Phys，2006，28（9）：916-924.

[11]Rohlmann A，Gabel U，Graichen F，et al.An instrumented implant for vertebral body replacement that measures loads in the anterior spinal column[J].Med Eng Phys，2007，29（5）：580-585.

[12]Rohlmann A，Petersen R，Schwachmeyer V，et al.Spinal loads during position changes[J].Clin Biomech（Bristol，Avon），2012，27（8）：754-758.

[13]Dickinson AS，Taylor AC，Ozturk H，et al.Experimental validation of a finite element model of the proximal femur using digital image correlation and a composite bone model[J].J Biomech Eng，2011，133（1）：014504.

[14]Karunratanakul K，Kerckhofs G，Lammens J，et al.Validation of a finite element model of a unilateral external fixator in a rabbit tibia defect model[J].Med Eng Phys，2013，35（7）：1037-1043.

[15]Ghosh R，Gupta S，Dickinson A，et al.Experimental validation of finite element models of intact and implanted composite hemipelvises using digital image correlation[J].J Biomech Eng，2012，134（8）：081003.

[16]Tseng ZJ，McNitt-Gray JL，Flashner H，et al.Model sensitivity and use of the comparative finite element method in mammalian jaw mechanics：mandible performance in the gray wolf[J].PLoS One，2011，6（4）：e19171.

[17]陳惠恩，鐘永青.QCT應當在骨質疏松癥的診斷、治療中發揮更大的作用[J].中國骨質疏松雜志，2008，14（11）：803-807.

[18]Wang X，Sanyal A，Cawthon PM，et al.Prediction of new clinical vertebral fractures in elderly men using finite element analysis of CT scans[J].J Bone Miner Res，2012，27（4）：808-816.

[19]Graeff C，Marin F，Petto H，et al.High resolution quantitative computed tomography-based assessment of trabecular microstructure and strength estimates by finite-element analysis of the spine，but not DXA，reflects vertebral fracture status in men with glucocorticoid-induced osteoporosis [J].Bone，2013，52（2）：568-577.

[20]唐海，李永剛，陳浩.骨密度測量技術的臨床應用及評價[J].實用糖尿病雜志，2006，2（4）：9-11.

[21]Baqi CM，Hanson N，Andresen C，et al.The use of micro-CT to evaluate cortical bone geometry and strength in nude rats：correlation with mechanical testing，pQCTand DXA[J].Bone，2006，38（1）：136-144.

[22]Christiansen BA，Kopperdahl DL，Kiel DP，et al.Mechanical contributions of the cortical and trabecular compartments contribute to differences in age-related changes in vertebral body strength in men and women assessed by QCT-based finite element analysis[J].J Bone Miner Res，2011，26（5）：974-983.

[23]Keaveny TM，Donley DW，Hoffmann PF，et al.Effects of teriparatide and alendronate on vertebral strength as assessed by finite element modeling of QCT scans in women with osteoporosis[J].J Bone Miner Res，2007，22（1）：149-157.

[24]Lewiecki EM，Keaveny TM，Kopperdahl DL，et al.Once-monthly oral ibandronate improves biomechanical determinants of bone strength in women with postmenopausal osteoporosis[J].J Clin Endocrinol Metab，2009，94（1）：171-180.

10.3969/j.issn.1009-5519.2015.04.023

：A

：1009-5519（2015）04-0546-03

2014-11-13）

湯凱文（1987-），男，湖南郴州人，碩士研究生，主要從事脊柱外科臨床研究；E-mail：48844879@qq.com。

熊波（E-mail：boxiongcz@sina.com）。