胸腰段血管瘤有限元模型建立及生物力学分析

doi:10.3969/j.issn.2095-4344.2017.11.020

中国组织工程研究 ›› 2017, Vol. 21 ›› Issue (11): 1753-1757.doi: 10.3969/j.issn.2095-4344.2017.11.020

• 骨与关节生物力学 bone and joint biomechanics • 上一篇下一篇

胸腰段血管瘤有限元模型建立及生物力学分析

江仁兵，董磊，刘启震，武峻申，白靖平

新疆医科大学附属肿瘤医院骨与软组织病区，新疆维吾尔自治区乌鲁木齐市新市 830000

出版日期:2017-04-18 发布日期:2017-05-06
通讯作者: 白靖平，博士，教授，新疆医科大学附属肿瘤医院骨与软组织病区，新疆维吾尔自治区乌鲁木齐市新市 830000
作者简介:江仁兵，男，1968年生，新疆维吾尔自治区乌鲁木齐市新市人，汉族，2014年新疆医科大学毕业，博士，主任医师，副教授，主要从事骨与软组织肿瘤的研究。
基金资助:
国家自然科学基金项目(81360282)

Finite element model establishment of thoracolumbar hemangioma and biomechanical analysis

Jiang Ren-bing, Dong Lei, Liu Qi-zhen, Wu Jun-shen, Bai Jing-ping

Ward of Bone and Soft Tissue, Tumor Hospital Affiliated to Xinjiang Medical University, Urumqi 830000, Xinjiang Uygur Autonomous Region, China

Online:2017-04-18 Published:2017-05-06
Contact: Bai Jing-ping, M.D., Professor, Ward of Bone and Soft Tissue, Tumor Hospital Affiliated to Xinjiang Medical University, Urumqi 830000, Xinjiang Uygur Autonomous Region, China
About author:Jiang Ren-bing, M.D., Chief physician, Associate professor, Ward of Bone and Soft Tissue, Tumor Hospital Affiliated to Xinjiang Medical University, Urumqi 830000, Xinjiang Uygur Autonomous Region, China
Supported by:
the National Natural Science Foundation of China, No. 81360282

摘要/Abstract

摘要：

文章快速阅读：

文题释义：

有限元分析：有限元法是根据变分法原理求解数学上可描述的物理问题的一种数值计算方法。它的主要优点是能够解决结构、材料性质和载荷情况比较复杂的问题。但有限元方法结果受诸多因素的影响，研究的基础和前提是建立精确的模型。

等效应力：等效应力是弹塑性力学里的VonMises应力，主要考察的是材料在各个方向上的应力差值，因为在实验室里获得的材料强度都是单向载荷作用下的强度，所以有时候材料所受的单向载荷可能很大，但并没有造成破坏，这就是看其等效应力。

摘要

背景：关于胸椎段肿瘤生物力学方面的研究比较多，但是有关有限元分析胸腰段椎体血管瘤骨折风险的研究不多。

目的：建立胸腰段椎体血管瘤有限元模型，分析其生物力学分布特点，评估脊椎血管瘤的骨折风险。

方法：建立T₁₂-L₂椎体的正常椎体、血管瘤椎体(血管瘤占椎体松质骨比例分别为20%，40%，60%，80%)、骨水泥填充椎体三维有限元模型，并进行力学分析。分别予以600 N垂直静态压力下正常椎体、血管瘤椎体、骨水泥填充的椎体等效应力分布及特点。

结果与结论：①建立了有效的T₁₂-L₂椎体有限元模型。在静态压力作用下3组模型L₁椎体皮质骨应力传导无明显差异，应力最大部均出现在椎弓根基底部、关节突关节及峡部；②在骨皮质完整情况下，椎体血管瘤占椎体松质骨20%-40%时病变椎体与正常椎体应力分布无统计学差异；当比例为60%-80%时，病变椎体与正常椎体应力分布有统计学差异；③建立的胸腰段三维有限元模型具有可用性；生物力学实验显示，椎体皮质骨完整性、松质骨破坏比例是椎体血管瘤骨折风险关键因素。

中国组织工程研究杂志出版内容重点：人工关节；骨植入物；脊柱；骨折；内固定；数字化骨科；组织工程
ORCID: 0000-0003-1069-693X(江仁兵)

关键词: 骨科植入物, 脊柱植入物, 有限元, 椎体血管瘤, 等效应力, 国家自然科学基金

Abstract:

BACKGROUND: Most of studies focus on the biomechanical characteristics of thoracic spine neoplasm, but there is little report on the fracture risk in the patients with vertebral hemangioma through finite element analysis.

OBJECTIVE: To establish a finite element model of vertebral hemangioma, and to analyze its biomechanical characteristics, and assess the risk of vertebral fracture.

METHODS: Three-dimensional finite element models of T₁₂-L₂ vertebrae from normal individuals, the patients with vertebral hemangioma (hemangioma accounting for 20%, 40%, 60%, 80% of the vertebral cancellous bone) and bone cement filling were established, respectively, and then the mechanical characteristics were analyzed. The stress distribution and characters of each model were determined under a vertical static pressure of 600 N.

RESULTS AND CONCLUSION: (1) Three-dimensional finite element models of T₁₂-L₂ vertebrae were established successfully. Under static pressure, the stress distribution of L₁ cortical bone showed no significant difference among models, and the maximum stress all occurred at the base of pedicle, zygapophysial joint and isthmus. (2) The stress distribution did not differ significantly between vertebral hemangioma accounting for 20%-40% of vertebral cancellous bone with complete cortical bone and normal ones, but which differed significantly in hemangioma accounting for 60%-80%. (3) To conclude, the established thoracolumbar three-dimensional model is available. Additionally, biomechanical tests manifest that the completeness of cortical bone and destruction ratio of cancellous bone destruction are key factors for the fracture risk of vertebral hemangioma.

中国组织工程研究杂志出版内容重点：人工关节；骨植入物；脊柱；骨折；内固定；数字化骨科；组织工程

Key words: Finite Element Analysis, Spine, Hemangioma, Stress, Mechanical, Tissue Engineering

中图分类号:

R318

江仁兵，董磊，刘启震，武峻申，白靖平. 胸腰段血管瘤有限元模型建立及生物力学分析[J]. 中国组织工程研究, 2017, 21(11): 1753-1757.

Jiang Ren-bing, Dong Lei, Liu Qi-zhen, Wu Jun-shen, Bai Jing-ping. Finite element model establishment of thoracolumbar hemangioma and biomechanical analysis [J]. Chinese Journal of Tissue Engineering Research, 2017, 21(11): 1753-1757.

图/表（结果） 3

参考文献

[1] 李健,黄海,张亮,等.经皮椎体成形术治疗症状性颈椎血管瘤1例[J].中国微创外科杂志,2011,11(12):1143-1144.

[2] 万宁军,胡永成,王晗.椎板减压联合椎体成形术治疗合并脊髓功能损害的脊椎血管瘤[J].中华骨科杂志,2012,32(11):1001-1004.

[3] 姜万荣,李兵,朱锡旭.射波刀治疗复发性椎体血管瘤1例报告及文献复习[J].现代肿瘤医学,2012,20(5):1034-1036.

[4] 林松青,王彬,张磊,等.有限元分析在骨科中的应用及研究进展[J].中国中医骨伤科杂志,2013,21(4): 69-73.

[5] 梁得华,鲁世保.有限元分析法在胸腰段脊柱骨折的应用进展[J].北京医学,2016,38(7):706-708.

[6] 张振军,李阳,廖振华,等.有限元法在腰椎生物力学应用中的研究进展和展望[J].生物医学工程杂志,2016,33(6):1097-1202

[7] 王海羽,王挺.有限元法分析全膝关节置换对胫骨近端骨重建的影响[J].中国组织工程研究,2016,20(35):5180-5186.

[8] 张海峰,尹爱华,董毅,等.有限元法分析不同负荷下髋臼区的应力分布[J].中国组织工程研究,2016,20(39):5867-5872.

[9] Zauel R, Yeni YN, Bay BK, et al. Comparison of the linear finite element prediction of deformation and strain of human cancellous bone to 3D digital volume correlation measurements. J Biomech Eng. 2006;128(1):1-6..

[10] 李大森,郭卫,杨荣利,等．症状性脊椎血管瘤的术式选择和疗效分析[J]中国脊柱脊髓杂志,2015,25(2):97-102．

[11] 杨保家,杨开舜,姚汝斌.影像及生物力学测试尼古丁剂量与腰椎后外侧植骨融合率的相关性[J].中国组织工程研究, 2016,20(48): 7251-7260.

[12] 沈彬,孟阳,赵卫东,等．症状性椎体血管瘤影像学表现及手术治[J]中国脊柱脊髓杂志,2013,23(3):251-256

[13] Cohen JE,Lylyk P,Ceratto R,et al.Vertebroplasty: technique and results in 192 procedures.Neurol Res.2004;26(1): 41-49.

[14] Peh WC G,Munk PL,Rashid F,et al.Percutaneous Vertebral Augmentation: Vertebroplasty,Kyphoplasty and Skyphoplasty. Radiologic Clinics of North America,2008,46(3): 611-35,vii.

[15] 姜亮,李杰,刘忠军,等.脊柱血管瘤的诊断与治疗[J]中国脊柱脊髓杂志,2011,21(1):38-54．

[16] Fei Q,Zhao F,Yang Y,et al.[Effects of posterior lumbar spinal fusion on the stability of unstable lumbar segment and biomechanical properties of adjacent segments: a finite element study. Zhonghua Yi Xue Za Zhi.2015;95(45):3681.

[17] Su Y,Ren D,Jiang M,et al.T1 finite element model of Kümmell's disease shows changes in the vertebral stress distribution. Int J Clin Exp Med. 2015;8(11):20046-20055.

[18] Yan L,Chang Z,Xu Z,et al.Biomechanical effects of bone cement volume on the endplates of augmented vertebral body: a three-dimensional finite element analysis[J].中华医学杂志:英文版,2014,127(1):79-84.

[19] Erbulut DU,Zafarparandeh I,Hassan C R,et al.Determination of the biomechanical effect of an interspinous process device on implanted and adjacent lumbar spinal segments using a hybrid testing protocol: a finite-element study. J Neurosurg Spine. 2015;23(2):200-208.

[20] Cao KD,Grimm MJ,Yang KH.Load sharing within a human lumbar vertebral body using the finite element method.Spine. 2001;26(12):253-260.

[21] 王宇, 靳安民, 方国芳,等. 腰椎椎弓根螺钉系统断裂的三维有限元分析[J]. 中国组织工程研究, 2008, 12(48):9439-9442.

[22] Crawford RP,Keaveny TM.Relationship between axial and bending behaviors of the human thoracolumbar vertebra. Spine. 2004;29(20):2248-2255.

[23] 鲍春雨,刘晋浩.人体脊柱腰椎节段三维有限元模型的研究与分析[J].机械设计, 2009, 26(9):61-63.

[24] Dall'Ara E,Schmidt R,Pahr D,et al.A nonlinear finite element model validation study based on a novel experimental technique for inducing anterior wedge-shape fractures in human vertebral bodies in vitro. J Biomech. 2010;43(12): 2374-2380.

[25] Oktenoglu T, Erbulut DU, Kiapour A, et al. Pedicle screw-based posterior dynamic stabilisation of the lumbar spine: in vitro cadaver investigation and a finite element study. Comput Methods Biomech Biomed Engin. 2015;18(11): 1252-1261.

[26] 刘治华,徐新伟,管文浩,等.腰椎有限元模型的建立与不同角度牵引条件下的仿真研究[J].郑州大学学报医学版, 2014(1): 119-122.

[27] Iselin L D,Wahl P,Studer P,et al.Associated lesions in posterior wall acetabular fractures: not a valid predictor of failure.J Orthop Traumatol.2013;14(3):179-184.

引言

材料方法

1.1 设计 有限元模型建立。

1.2 时间及地点 2013年10月至2014年4月在新疆医科大学实验室完成。

1.3 对象 身体健康的成年女性志愿者，40岁，身高 165 cm，体质量60 kg，无脊椎疾病及脊椎外伤，无骨质疏松及其他慢性病史。对实验方案了解同意并签署知情同意书书。详细查体及影像检查未见其他病变，其L₁椎体血管瘤。

1.4 实验方法

1.4.1 有限元模型建立采用的是64排螺旋CT平行扫描，床进速度 0.707 mm/s，扫描层厚0.625 mm，螺旋扫描参数为：120 kV，36 mA。

将CT扫描的DICOM文件导入Mimics15.01软件中，形成蒙皮，通过光滑处理模块进行光滑处理，通过调整阈值，提取椎间盘轮廓，研究者手绘椎间盘大体模型。导入光滑处理软件中进行光滑处理，最终得到光滑后的椎间盘三维模型。前纵韧带、后纵韧带、棘上韧带的3D模型的建立依照上面方法生成。使用布尔操作命令生成1 mm厚的皮质骨和内部松质骨的椎体结构，同样得到包含纤维环和髓核的椎间盘结构，最后完成实体模型的建立。

分别建立T₁₂-L₂节段正常、病变、骨水泥填充3种情况下的有限元模型，利用 Meshing Tools中的volume划分工具对骨骼模型进行自动网格划分，局部实施网格细化，最后得到可以进行试验的有限元模型，其中每节段脊椎中的上下2个椎体设定为正常椎体，中间为病变椎体。有限元法可以有效地模拟临床手术行为和术后影响^[6]。

1.4.2 各种结构不同的属性材料参数分别赋予各种结构不同的属性材料参数(见表1)。

1.4.3 静态生物力学实验病变椎体下方、双侧下关节面以及棘突下方予以固定。在已建立的T₁₂-L₂节段有限元模型椎体上表面放有密度均匀的3 mm厚的刚体，覆盖范围包括椎体及后部结构，在椎体及附件定义16个应力测量点(详见图1)。

1.5 主要观察指标 ①检测给予600 N的垂直静态压力下病变椎体、正常椎体及经过骨水泥填充的椎体的等效应力分布及特点；②测量600 N垂直静态压力下病变椎体血管瘤占据椎体松质骨20%，40%，60%，80%时等效应力分布及特点。

1.6 统计学分析 将数据入SPSS 20.0统计软件，进行配对样本的t 检验，检验标准P < 0.05。

中国组织工程研究杂志出版内容重点：人工关节；骨植入物；脊柱；骨折；内固定；数字化骨科；组织工程

讨论

三维有限元分析方法因为其克服了传统生物力学的缺点，所以在脊柱生物力学的研究中广泛应用^[7]。其通过离散化的有限单元代替原有物体，结点连接单元，而力则通过结点传递，从而成为了椎体生物力学研究的常用工具之一。应用骨是一种人体中的结缔组织，有着支撑身体、保护内脏器官及维持姿势等作用，并且可以在肌肉的牵拉下完成各种活动^[8]。Zauel 等^[9]首次利用有限元分析法，分析了单轴向压力作用对人体松质骨骨小梁所产生的形变情况。椎体血管瘤是临床上常见的椎体病变，有研究发现尸检率约为10%，较多发生于胸椎，但大多数胸椎血管瘤患者无症状，无需处理^[10]。而其中0.9%-1.2%的血管瘤会出现膨胀生长，导致压迫神经组织，最终出现神经症状，需要手术治疗。随着近年来内固定技术的发展，后路植骨融合在临床应用中更为普遍^[11]。病椎在MRI中T₁上表现为中等或低信号而T₂上则为高信号，在CT中则表现为骨小梁紊乱增粗，骨质稀疏^[12]。MRI表现为椎体内T₁、T₂高信号，类似脂肪的改变，一般提示血管瘤不活跃^[13]。无症状的患者随访观察即可，而对于严重的椎体血管瘤如果不及时的治疗则可能出现严重的后果^[14]。大量文献报道，对于椎体血管瘤应行经皮椎体成形术用以稳定椎体，因为很多专家认为椎体血管瘤骨折风险远高于此^[15]，是考虑到椎体血管瘤可能造成的椎体压缩性骨折所带来的风险。而椎体血管瘤是否真的会导致骨折风险增加及哪种类型的血管瘤会增加骨折风险至今尚无明确定论。

1975年，Liu首次报告了腰椎的三维有限元模型。Hakim和King模拟了腰椎后部结构特别是小关节的作用，进行了静力学和动力学分析。有专家通过建立胸腰段的有限元模型发现，髓核中心及后偏移是上下终板高应力的部位。Fei等^[16]经过有限元建模并行生物力学实验证实：椎体皮质骨、终板、椎间盘分别承受43%、10%、13%的垂直应力。Su等^[17]提出当强度在3 kN与5 kN之间时会存在不同的骨折风险，而当强度高于5 kN时是没有骨折风险的，在强度低于3 kN的时候骨折风险为100%。而本研究中在给予600 N的正常压力下，正常的椎体及病椎的应力值均低于此数值，而实验结果显示，此时椎体血管瘤仅占椎体松质骨的15%。椎间盘所受应力分布是不均的，各个椎体的前部所受的应力最大，而后部结构所受应力较小，椎体局部出现受力较大情况，尤以椎体相邻部位更为明显。

在椎体应力分布上，Yan等^[18]在椎体应力分布上采用挖空骨松质仅以外壳测量椎体轴向符合的办法已证实椎体骨密质可提供75%的椎体强度。相反，Erbulut等^[19]则在胸腰椎标本中去掉骨皮质后检测松质骨的轴向负荷抗力，实验结果表明男性中皮质的负荷抵抗力为43.8%(29%- 54%)，而女性中则为35.2%(22%-44%)。Cao等^[20]建立了有限元分析模型，分析结果标明34%的椎体压缩骨折发生在上终板或下终板附近，而63%的椎体压缩骨折发生在上下终板之间。王宇等^[21]研究发现椎弓根螺钉各部分所承受的垂直应力远小于其他屈伸时所承受的应力。Crawford等^[22]研究认为椎体屈曲强度对椎体骨折有之间影响，而椎体的前后径有对椎体的屈曲强度有影响。鲍春雨等^[23]研究发现椎间盘内的应力不仅在前屈位大于后伸位并且有着自上而下增大的规律。Dall Ara等^[24]通过非线性解剖特异的定量CT数据位依据建立的有限元分析模型，预测了体内椎体的力量。Oktenoglu等^[25]的研究发现完整的的脊柱与加装稳定系统的脊柱在位移上没有差别。作者研究发现脊柱应力主要集中在终板中央，终板中央又是其结构的薄弱部位，如果持续的增加应力可能导致骨折的可能性增加，而临床上出现终板骨折的情况较多，并与上述文献相符。

为了进一步研究椎体血管瘤造成骨松质破坏对椎体应力所造成的影响，作者实验中把有限元模型设置松质骨缺损，采用同心圆的形式设计了缺损分别为20%，40%，60%，80%，同时将缺损赋值后采用弹性模量及泊松比模拟出了病变占椎体松质骨20%，40%，60%，80%的情况，在不同比例的情况下完成静态加压600 N，检测其中16个点的应力变化，结果与正常椎体在相同压力下做比较分析后证实，病椎与正常椎体终板应力均未超越3.6 MPa，椎体血管瘤在骨皮质完情况下，椎体血管瘤占椎体松质骨20%-40%的时候，正常椎体与病椎承受应力大小是一致的。这与既往研究结果致^[26]。椎体血管瘤在骨皮质完整情况下，椎体血管瘤占椎体松质骨60%-80%的时候，正常椎体与病椎所承受应力大小是不一致的，结果有统计学意义，这也进一步证明了在骨皮质完整时，椎体病变与椎体松质骨的比例与椎体骨折风险相关，同时也证实明了当松质骨缺损大于60%以上时，病椎的骨折风险较正常椎体大。而当上下终板有应力超越3.6 MPa，骨折风险也会增加。

精确的诊断和正规的治疗计划对获得良好功能及较好远期疗效是很重要的^[27]。胸腰段椎体血管瘤三维有限元模型能够较好的模拟胸腰椎血管瘤实际情况，可用于临床对胸腰椎血管瘤患者的生物力学分析，其结果可以指导制定胸腰椎血管瘤的治疗策略。有限元分析能够将复杂问题简单化，减少血管瘤骨折及术后并发症的发生，为胸椎血管瘤骨折的治疗提供新的科学依据。本研究结果显示：血管瘤对椎体的松质骨破坏比例及皮质骨是否完整是导致骨折的危险因素，由于椎体血管瘤很少发生在关节突及椎弓根，这些部位尚未考虑。有限元分析法仍需要较多的临床观察及术后的长期随访才能对胸椎血管瘤的发生发展过程提供参考。

中国组织工程研究杂志出版内容重点：人工关节；骨植入物；脊柱；骨折；内固定；数字化骨科；组织工程

胸腰段血管瘤有限元模型建立及生物力学分析

Finite element model establishment of thoracolumbar hemangioma and biomechanical analysis

PDF

可视化

摘要/Abstract

引用本文

使用本文

图/表（结果） 3

参考文献

相关文章 15

引言

材料方法

讨论

文章快阅

延伸阅读

编辑推荐

Metrics

本文评价

[1]	魏国强, 李云峰, 王一, 牛晓芬, 车丽芳, 王海燕, 李志军, 史国鹏, 白灵, 莫凯, 张晨晨, 许阳阳, 李筱贺. 非均匀材料股骨在不同负荷情况下的生物力学分析[J]. 中国组织工程研究, 2022, 26(9): 1318-1322.
[2]	吕倩忆, 陈芯仪, 郑慧娥, 何灏龙, 李棋龙, 陈楚淘, 田浩梅. 不同角度肘按法对正常人腰椎椎体及后部结构的应力和位移分析[J]. 中国组织工程研究, 2022, 26(9): 1346-1350.
[3]	张玉芳, 吕蒙, 梅钊. 青少年脊柱侧弯全脊柱生物力学模型的构建及验证[J]. 中国组织工程研究, 2022, 26(9): 1351-1356.
[4]	张吉超, 董跃福, 牟志芳, 张震, 李冰言, 徐祥钧, 李佳意, 任梦, 董万鹏. 骨关节炎患者在不同步态角度下膝关节内部生物力学变化的有限元分析[J]. 中国组织工程研究, 2022, 26(9): 1357-1361.
[5]	白子兴, 曹旭含, 孙承颐, 杨艳军, 陈思, 温建民, 林新晓, 孙卫东. 步态周期中踝关节有限元模型的构建及生物力学分析[J]. 中国组织工程研究, 2022, 26(9): 1362-1366.
[6]	刘峰, 冯毅. 步态周期下不同克氏针张力带治疗髌骨横行骨折的有限元分析[J]. 中国组织工程研究, 2022, 26(9): 1367-1371.
[7]	李桂军, 方晓辉, 孔维峰, 袁晓庆, 金荣忠, 杨俊. 微型钢板联合超强缝线弹性固定治疗拇外翻畸形的有限元分析[J]. 中国组织工程研究, 2022, 26(6): 938-942.
[8]	温明韬, 梁学振, 李嘉程, 许波, 李刚 . 两种方式固定Sanders Ⅱ型跟骨骨折后的力学稳定性[J]. 中国组织工程研究, 2022, 26(6): 838-842.
[9]	王海龙, 李龙, 买合木提·亚库甫, 陈洪涛, 刘旭, 伊力哈木·托合提. 髋臼假体在Lewinnek 安全区内应力分布的有限元分析[J]. 中国组织工程研究, 2022, 26(6): 843-847.
[10]	黄浩, 洪嵩, 瓦庆德. 全膝关节置换过程中股骨假体旋转对髌股关节接触压影响的有限元分析[J]. 中国组织工程研究, 2022, 26(6): 848-852.
[11]	郑永泽, 郑利钦, 何兴鹏, 陈心敏, 李木生, 李鹏飞, 林梓凌. 基于ABAQUS软件股骨颈骨折的扩展有限元建模分析[J]. 中国组织工程研究, 2022, 26(6): 853-857.
[12]	李硕, 苏鹏, 张力, 吴秋龙, 胡湘宇, 来钰梁. 三维重建和有限元分析股骨髁上截骨对矫治膝内翻有积极作用[J]. 中国组织工程研究, 2022, 26(6): 858-863.
[13]	魏兵, 常山. 脊柱骨折矢状面不同角度置钉方式的有限元分析[J]. 中国组织工程研究, 2022, 26(6): 864-869.
[14]	刘宇航, 周建强, 徐雪彬, 曲星月, 李梓瑜, 李琨, 王星, 李志军, 李筱贺, 张少杰. 6岁儿童下颈椎有限元模型建立及有效性验证[J]. 中国组织工程研究, 2022, 26(6): 870-874.
[15]	高文波, 马宗民, 李淑娴, 聂秀吉. 有限元分析不同骨质下种植体长度及直径对初期稳定性的影响[J]. 中国组织工程研究, 2022, 26(6): 875-880.