背景:随着检查技术的革新,不同分期脊柱转移瘤患者逐年增多,经皮椎体成形是脊柱转移瘤重要的治疗手段,但其治疗不同分期、不同活动方式肿瘤的术后生物力学疗效尚未见报道. 目的:基于三维有限元模型模拟胸椎T10骨肿瘤转移不同位置的应力及位移情况. 方法:根据1例30岁健康男性胸椎三维CT扫描数据,采用Mimics软件构建胸椎(T9-T11)的三维几何模型,包括肋骨、韧带及椎间盘.模拟被胸椎转移瘤侵袭的T9-T11椎体及其后胸椎不同部位的三维模型,包括椎体结构完整的对照组,单侧转移累及椎体区(实验组1),单侧转移累及椎体和椎弓根区(实验组2),单侧转移累及椎体、椎弓根和横突区(实验组3),双侧转移累及椎体、椎弓根和横突区(实验组4).利用Abaqus软件创建三维有限元模型,对负重、屈曲、伸展、旋转工况下Von Mises应力分布及模型的位移情况进行分析. 结果与结论:①在对不同实验组加载点在负重、屈曲、伸展、旋转工况下最大总位移的研究中,随着转移瘤侵袭部位及侵袭面的增多,加载点的总位移在增大,整体刚度降低,尤其实验组4加载点的总位移最大;②在屈曲工况下椎体、椎弓根破坏后明显增加最大Von Mises应力值,而在此基础上增加胸肋关节破坏,其最大Von Mises应力值基本不变;③有限元分析并模拟骨肿瘤模型的基础上,将骨水泥区域的单元设置为一个单独的set集合,后将骨水泥区域设置为对应的材料属性以模拟骨水泥填充,结果发现,负重、屈曲、伸展、旋转工况下最大总位移均小于各实验组;④模拟经皮椎体成形治疗后患者在负重、屈曲、伸展、旋转工况下的最大应力值相较于骨肿瘤模型大幅度降低;⑤提示基于胸椎T9-T11的三维有限元模型构建有利于对胸椎骨肿瘤转移生物力学特性的研究,并能够在胸椎骨肿瘤转移模型的基础上准确模拟经皮椎体成形后加载点在不同工况下的总位移及最大Von Mises应力情况.
BACKGROUND:With the innovation of examination technique,the number of patients with spinal metastases in different stages is increasing year by year.Percutaneous vertebroplasty is an important treatment for spinal metastases;however,there is no report on the biomechanical effect in different stages and different activities after operation. OBJECTIVE:To simulate thoracic T10 bone stress and displacement of the different locations of the tumor metastasis based on the three-dimensional finite element model. METHODS:According to thoracic three-dimensional CT images of a 30-year-old healthy male,Mimics software was used to construct a three-dimensional geometric model of thoracic vertebrae(T9-T11),including ribs,ligaments and intervertebral discs.Three-dimensional models of T9-T11 vertebral bodies and different parts of the posterior thoracic vertebrae invaded by thoracic metastatic tumors were simulated,including the control group with intact vertebral structure,unilateral metastasis involving the vertebral body area(experimental group 1),unilateral metastasis involving the vertebral body and pedicle area(experimental group 2),unilateral metastasis involving the vertebral body,pedicle and transverse process area(experimental group 3),and bilateral metastasis involving the vertebral body,pedicle and transverse process area(experimental group 4).Abaqus software was used to create a three-dimensional finite element model.The von Mises stress distribution and the displacement of the model were analyzed under the loading condition,buckling condition,extension condition,and rotation condition. RESULTS AND CONCLUSION:(1)In the study of the maximum total displacement of loading points in different experimental groups under loading,flexion,extension,and rotation conditions,with the increase of metastatic tumor invasion site and invasion surface,the total displacement of loading points increased,and the overall stiffness decreased,especially the total displacement of loading points in experimental group 4 was the largest.(2)Under flexion condition,the maximum Von Mises stress value increased significantly after vertebral body and pedicle destruction,while the maximum Von Mises stress value was almost unchanged when the thoracocostal joint destruction was added.(3)On the basis of finite element analysis and simulation of bone tumor model,the elements in the bone cement region were set as a single set,and the bone cement region was set as the corresponding material properties to simulate bone cement filling.The results showed that the maximum total displacement under loading,flexion,extension,and rotation conditions was less than that of each experimental group.(4)The maximum stress values of the simulated percutaneous vertebroplasty patients in the loading,flexion,extension and rotation conditions were significantly lower than those of the femoral model.(5)It is concluded that the three-dimensional finite element model based on thoracic T9-T11 conducive to the biomechanics characteristics of thoracic vertebrae tumor metastasis,and on the basis of the thoracic vertebrae tumor metastasis model can accurately simulate load point after percutaneous vertebral body under different conditions of total displacement and the maximum Von Mises stress situation.