股骨成形術預防老年髖部骨折的斷裂分析

孫文濤++++++李鵬飛++++++林梓凌++++++何祥鑫++++++李紅庚++++++魯榮貴++++++馮蓮影

[摘要] 目的 基于斷裂力學模擬股骨成形術前后股骨頸的斷裂過程，分析骨水泥干預對老年骨質疏松性髖部骨折的療效。 方法 選取1例骨質疏松股骨頸骨折的健側髖部CT資料以Dicom格式保存，導入Mimics三維建模軟件重建髖部骨骼有限元模型，將模型導入工程力學軟件HYPERMESH劃分網格、定義材料屬性、設置邊界條件，生成的為對照模型；在邊界條件不變的情況下，對股骨頸區域內部的有限元網格重新定義材料屬性，生成的為實驗模型；分別將兩個模型導入動力求解軟件LS-DYNA中進行運算，比較動態斷裂過程以及Von Mises云圖。 結果 對照模型股骨頸斷裂，骨折線由內下向外上方走行，Von Mises云圖上顯示股骨頸骨折部位應力較高；實驗模型股骨頸部位未見明顯骨折線，Von Mises云圖上顯示髖部骨骼周圍應力較對照組明顯降低。 結論 從生物力學角度上，股骨成形術對于預防骨質疏松性股骨頸骨折具有可行性、實用性，為臨床運用提供參考依據。

[關鍵詞] 股骨成形術；骨水泥；預防；骨質疏松；髖部骨折；斷裂分析

[中圖分類號] R683 [文獻標識碼] A [文章編號] 1673-7210（2017）07（c）-0094-03

[Abstract] Objective To analyze effect of bone cement intervention on osteoporotic hip fracture in elderly patients based on biomechanical simulation of femoral neck fracture before and after femoroplasty. Methods CT image data of one case of femoral neck fracture was collected. Mimics software was used to rebuild the three-dimensional finite element model of femoral. Model was imported in HYPERMESH software for meshing， defining material properties， setting the boundary conditions， loading the force curve and the role of points such as pre-processing， generated model for the control model. In the case of the same boundary conditions， the finite element grid inside the femoral neck region was redefined the material properties， and the generated model was selected as the experimental model. Then the two models were imported into the solver LS-DYNA software for operating respectively， then dynamic fracture processes and Von Mises cloud maps were observed. Results The Von Mises cloud showed a high stress on the femoral neck fracture. There was no obvious fracture line in the femoral neck of the experimental model， and the Von Mises cloud showed a hip bone around the hip stress was significantly lower than the control model. Conclusion From the perspective of biomechanics， femoral angioplasty is feasible and practical for the prevention of osteoporotic femoral neck fractures， and provides reference for clinical application

[Key words] Femoroplasty； Cemented； Prevention； Osteoporosis； Hip fracture； Fracture analysis

老年髖部骨折是最常見的損傷之一，隨著全球老齡化人口的增加，這一問題更加突出。在美國，與骨質疏松性骨折相關的費用估計高達10億～20億美元/年，其中約30萬人由于髖關節骨折入院[1-3]，估計有90億美元的直接醫療費將用于髖部骨折的治療[4-5]。據推測，到2025年，全球髖骨骨折的年發病率估計將增至300萬增至1000萬[6]。顯然，預防髖部骨折對于降低發病率、死亡率和減輕社會經濟負擔具有戰略性意義[7]。

骨質疏松主要表現為骨小梁減少，且骨小梁的抗斷裂減弱，輕微跌倒暴力即可導致骨質疏松性骨折。目前預防髖部骨折主要的措施是一方面口服刺激骨生成、抑制骨吸收的藥物[8]，另一方面適當運動鍛煉肌肉協調性防止跌倒以及佩戴髖保護墊[9]，臨床上有取得一定療效，但仍沒達到理想的標準。因此，人們開始著手進行外科微創手術進行預防髖部骨折的探索，國外研究的“股骨成形術”（股骨近端骨水泥強化增強骨小梁抗壓抗剪切力）[10]，它具有與骨組織良好的結合性，并表現出強大的抗斷裂作用[11-13]，可有效預防髖部骨折，但其動態的斷裂力學研究國內外尚未見報道。

團隊成員李鵬飛等[14]在髖部骨折的有限元仿真模擬中，開創性的運用LS-DYNA動態分析軟件模擬跌倒狀態下的股骨頸骨折過程，并且驗證仿真度高達83%，但股骨轉子間骨折由于邊界條件的設置以及技術原因等限制尚未模擬。基于此，本文僅以股骨頸骨折為模型探討股骨頸局部骨水泥強化前后的斷裂過程并作相關分析。

1 材料與方法

1.1 材料

選取2016年12月11日于廣州中醫藥大學第一附屬醫院（以下簡稱“我院”）一骨科經診斷為老年骨質疏松性股骨頸骨折1例（女性，68歲），癥狀：左側髖部腫痛，雙髖關節正位+左側髖部側位X線片報告：左側股骨頸骨折（Garden Ⅲ型）。診斷標準：①明顯外傷史；②有髖部疼痛、腫脹、活動受限等癥狀；③腹股溝中點處壓痛、下肢縱向叩擊痛，患肢呈外旋、短縮畸形；④X線片證實股骨頸骨折。經我院倫理委員會批準并與患者簽署知情同意書后，采用GE 64排螺旋CT對右側股骨全段螺旋掃描并三維重建（患者仰臥位），然后數據以Dicom格式保存。掃描條件：掃描電壓120 kV，掃描電流250 mA，自髖臼至股骨中上段層厚2 mm，層距5 mm，每個掃描層的象素矩陣密度大小為512×512。

1.2 方法

1.2.1 實驗軟件與設備 Mimics 19.0軟件（Materialise公司，比利時），由廣州中醫藥大學國家重點學科中醫骨傷科學數字骨科與生物力學實驗室提供；GE 64排螺旋CT，由廣州中醫藥大學第一附屬醫院提供；Hypermesh 14.0軟件（Altair公司，美國），由廣州有道計算機有限公司提供；大型有限元分析軟件LS-DYNA軟件（Lstc公司，美國），由湖南大學提供。

1.2.2 三維重建 將保存的Dicom格式的影像資料導入Mimics19.0，通過區域增長、編輯蒙罩、Calculate 3D重建，并通過包裹、光滑等工具對模型進行光順，得到理想的髖部三維有限元模型，導出為igs格式文件。

1.2.3 有限元網格劃分與定義材料屬性、邊界條件 將igs文件在Hypermesh中進行體網格劃分，并參考李鵬飛等[14]的模型定義材料屬性和邊界條件模型設定好以后導出為對照模型。

1.2.4 骨水泥干預 新建一個單元體，參考邱興等[15]定義骨水泥的材料屬性；隱藏股骨頸表面以及股骨頭和股骨干，將股骨頸內部的有限元網格通過Organize-move，使其轉移到新建的單元體，導出為實驗模型。見圖1，封三。

1.2.5 運算處理 將生成的兩個模型依次在LS-DYNA中進行運算，運算結果可在Hyperview中動態展示。

1.3 觀察指標

股骨頸骨折的裂紋走向，髖部骨骼模型的Von Mises云圖。

2 結果

對照模型在垂直股骨頭向下的壓力加載時，股骨頸斷裂并發生骨折，且骨折裂紋走向為內下向外上方向（圖2，封三），實驗模型股骨頸部位未見明顯的骨折裂紋，且股骨轉子間亦無明顯骨折裂紋（圖3，封三）。比較兩種模型股骨近端的Von Mises云圖（應力越高，云圖中顏色越鮮艷，紅色區域應力最高，為骨折危險區），可以發現：在相同的載荷等邊界條件下，對照模型股骨頸部位應力集中，為骨折斷裂區域；實驗模型骨近端（包括股骨轉子間）應力均較低；由此，股骨頸骨水泥填充可以有效預防骨質疏松性股骨頸骨折，而且可以減輕股骨轉子間、股骨干等部位的應力。見圖3～4，封三。

3 討論

張國棟等[16]根據Mimics提供的經驗公式[17]：Density=-13.4+1017×Grayvalue，E-Modulus=-388.8+5925×Density對骨骼定義多種材料屬性（密度、彈性模量、泊松比），并驗證了10種材料屬性即可滿足有限元分析要求，但是并未劃分股骨10種材料屬性的區域。林梓凌等[18]基于靜態力學探討了區域骨密度的老年性骨質疏松性的相關性，并建立了基于Abaqus靜態力學分析的有限元骨質疏松，為分析老年髖部骨折奠定基礎。運用斷裂力學研究分析問題，要求根據失效準則明確材料失效參數。然而Mimics經驗公式尚不能定義10種材料失效參數，無法精確導入LS-DYNA動態力學軟件分析，因此無法精確進行骨質疏松建模。本研究中進行簡化模型，視股骨為均質性材料，將骨密度定義為正常骨密度的2/3進行仿真骨質疏松模型，陳振沅[19]在有限元分析中使用該種模型模擬股骨轉子間骨折，并驗證了其有效性。

股骨成形術預防髖部骨折尚處于研究階段，因為技術的不成熟以及可能導致并發癥的發生，致使其尚未在臨床全面應用。由于骨水泥填充股骨近端造成的局部壓力增大可能會出現微血管栓塞，以及骨水泥放熱出現局部骨組織壞死等并發癥的問題[13]，Sutter等[20]研究發現40～50 mL骨水泥填充股骨干、股骨頸、股骨頭、股骨轉子間都具有抗骨脆性骨折的作用，但是填充股骨頸部位最具有顯著性意義；通過計算機軟件編程智能化控制經皮股骨近端骨水泥的用量與部位，可以達到增強骨密度的同時也減少了因骨水泥的熱效應產生的并發癥，也只是處于初步研究階段[21]。

本文的不足之處以及改進措施：首先，僅僅只模擬了股骨頸局部骨水泥強化對于抗骨質疏松骨折的斷裂分析，實驗尚不完整，后期再深入的基礎上可以進一步對股骨轉子間的骨水泥干預進行斷裂模擬，觀察其抗骨質疏松效果；第二，填充股骨近端骨水泥的容量是股骨成形術的一個研究重點，而生物有限元無法進行精確的定量模擬，或者說以筆者現有的有限元技術還不能進行此類操作，股骨成形的有限元斷裂研究尚處于起步階段，隨著計算機有限元的發展以及工程力學軟件性能的提升，對骨水泥進行精準定量模擬肯定能成為現實；第三，將股骨近端有限元模型簡化為均質的材料屬性，運算方便快捷，然而卻忽略了Wardˊ三角區對于髖部骨折的重大意義[22]，因此，將股骨材料屬性定義精準化是提高仿真度的重要一步。

[參考文獻]

[1] Elffors L. Are osteoporotic fractures due to osteoporosis？Impacts of a frailty pandemic in an aging world [J]. Aging （Milano），1998，10（3）：191-204.

[2] Ray WA，Taylor JA，Meador KG，et al. A randomized trial of a consultation service to reduce falls in nursing homes [J]. JAMA，1997，278（7）：557-562.

[3] Johnell O. The socioeconomic burden of fractures：today and in the 21st century [J]. Am J Med，1997，103（2A）：20S-25S.

[4] Stevens JA，Olson S. Reducing falls and resulting hip fractures among older women [J].MMWR Recomm Rep，2000， 5（14）：134-139.

[5] Wehren LE，Magaziner J. Hip fracture：risk factors and outcomes [J]. Curr Osteoporos Rep，2003，1（2）：78-85.

[6] Resch H，Gollob E，KudlacekS，et al. Osteoporosis in the man [J]. Wien Med Wochenschr，2001，151（18-20）：457-463.

[7] Melton LR. How many women have osteoporosis now？[J] J Bone Miner Res，1995，10（2）：175-177.

[8] 劉飛，史曉林，吳連國，等.唑來膦酸治療骨質疏松性股骨粗隆間骨折的臨床研究[J].中國中醫骨傷科雜志，2015（9）：19-21.

[9] Kannus P，Parkkari J，Niemi S，et al. Prevention of hip fracture in elderly people with use of a hip protector [J]. N Engl J Med，2000，343（21）：1506-1513.

[10] Oden ZM，Selvitelli DM， Bouxsein ML. Effect of local density changes on the failure load of the proximal femur [J]. J Orthop Res，1999，17（5）：661-667.

[11] Courtney AC，Wachtel EF，WyersER，et al. Effects of loading rate on strength of the proximal femur [J]. Calcif Tissue Int，1994，55（1）：53-58.

[12] Beckmann J，Springorum R，VettorazziE，et al. Fracture prevention by femoroplasty--cement augmentation of the proximal femur [J]. J Orthop Res，2011，29（11）：1753-1758.

[13] Beckmann J，Ferguson SJ，GebauerM，et al. Femoroplasty—augmentation of the proximal femur with a composite bone cement——feasibility，biomechanical properties and osteosynthesis potential [J]. Med Eng Phys，2007，29（7）：755-764.

[14] 李鵬飛，杜根發，林梓凌，等.基于LS-DYNA模擬老年股骨頸骨折的有限元分析[J].中國組織工程研究，2016， 20（44）：6606-6611.

[15] 邱興，楊圣，蘆健民，等.骨質疏松骨折椎體及相鄰椎體骨水泥注入前后的有限元分析[J].中國組織工程研究，2013，17（17）：3041-3048.

[16] 張國棟，廖維靖，陶圣祥，等.股骨頸有限元分析的賦材料屬性方法探討及有效性驗證[J].中國組織工程研究與臨床康復，2009，13（52）：10263-10268.

[17] Hobatho MC，Rho JY，Ashman RB. Anatomical variation of human cancellous bone mechanical properties in vitro [J]. Stud Health Technol Inform，1997，40（40）：157-173.

[18] 林梓凌，李鵬飛，龐智暉，等.骨密度與老年髖部骨折股骨近端三維有限元模型密度的關系[J].中國老年學雜志，2015，35（11）：3069-3070.

[19] 陳振沅.股骨轉子部骨折六部分骨折分型產生機制的有限元分析[D].遵義：遵義醫學院，2015.

[20] Sutter EG，Mears SC，Belkoff SM. Abiomechanical evaluation of femoroplasty under simulated fall conditions [J]. Orthop Trauma，2010，24（2）：95-99.

[21] Santana AM，Venetsanos DT. A new evolutionary optimization method for osteoporotic bone augmentation [J]. Comput Methods Biomech Biomed Engin，2017，20（7）：691-700.

[22] 唐瀟，錢競光，宋雅偉，等.股骨頸頸干角、骨密度與股骨頸骨折的關聯分析[J].南京體育學院學報：自然科學版，2006，5（1）：29-32.