低弹性模量骨水泥对骨质疏松压缩性骨折椎体及邻近椎体应力的影响：三维有限元分析

1、www.CRTER.org包拥政，等. 低弹性模量骨水泥对骨质疏松压缩性骨折椎体及邻近椎体应力的影响：三维有限元分析低弹性模量骨水泥对骨质疏松压缩性骨折椎体及邻近椎体应力的影响：三维有限元分析包拥政1，祝周兴2，冯云升3，吴 强1，胡孔和1，何小龙1，朱文刚1，席新华1，钟学仁1，周龙泽1，廖军健1，戴相恒1 (1汕头大学医学院附属粤北人民医院骨科四区，广东省韶关市 512026，2广东韶关乐昌市人民医院骨科，广东省韶关市 512200，3广东韶关南雄市中医院，广东省韶关市 512400)引用本文：包拥政，祝周兴，冯云升，吴强，胡孔和，何小龙，朱文刚，席新华，钟学仁，周龙泽，廖军健，戴相恒.

2、低弹性模量骨水泥对骨质疏松压缩性骨折椎体及邻近椎体应力的影响：三维有限元分析J.中国组织工程研究，2016，20(16): 2285-2293.DOI: 10.3969/j.issn.2095-4344.2016.16.001 ORCID: 0000-0001-9380-5438(包拥政)文章快速阅读：低弹性模量骨水泥对骨质疏松压缩性骨折椎体及邻近椎体应力的影响结局意义：骨水泥注射后骨折椎体及邻近椎体应力均较治疗前明显增加，骨折椎体随骨水泥弹性模量的增加应力也随之增大，而邻近椎体应力无明显变化。包拥政，男，1985年生，广东省韶关市人，汉族，南方医科大学在读博士，主治医师，主要从事脊柱外科方面

3、的研究。 中图分类号:R318文献标识码:A文章编号:2095-4344(2016)16-02285-09稿件接受：2016-02-12http:/WWW.方法：首次模拟椎体成形注射不同弹性模量骨水泥后各运动状态下骨折椎体及邻近椎体的应力改变。造模：建立腰椎骨质疏松椎体压缩性骨折模型。文题释义：弹性模量：指材料在弹性变形阶段内正应力和对应的正应变的比值。骨水泥的弹性模量明显大于骨质疏松的骨皮质及松质骨，注入骨水泥后使椎体的刚度及强度明显增加，进而可能导致再发骨折。低弹性模量骨水泥：降低骨水泥弹性模量能是减少修复后引起骨折椎体及邻近椎体应力变化而发生骨折的一种方法。临床常用的聚

4、甲基丙烯酸甲酯骨水泥可通过加入某些物质降低凝固后的弹性模量，但目前尚未发现较为理想的物质或本身较低弹性模量的骨水泥。摘要背景：椎体成形及椎体后凸成形可有效修复骨质疏松性椎体压缩性骨折，但注入骨水泥强化后可使骨折椎体及邻近椎体应力改变易导致新发骨折。目的：应用三维有限元法分析不同弹性模量骨水泥注射椎体成形后骨折椎体及邻近椎体的应力变化。方法：选取1个有代表性的健康成年男性志愿者行腰椎CT扫描，将图像处理后导出应用Mimics进行三维重建，应用Geomagic对三维模型进行光滑、打磨、去噪，Ansys Workbench下装配实体模型，赋值后建立L2-L4段骨质疏松椎体压缩性骨折模型。设定在L3椎

5、体注入4 mL弹性模量分别为8 000，4 000，2 000及1 000 MPa的骨水泥，骨水泥在椎体中央呈球形分布。在L2上表面施加500 N预载荷，附加弯矩为50 Nm，约束L4下表面自由度。模拟L2-L4前屈、后伸、右屈及右侧轴向旋转4种运动状态，比较不同弹性模量骨水泥注射前后骨折椎体及上下邻近椎体的应力变化。结果与结论：骨水泥注射后骨折椎体及邻近椎体应力均较治疗前明显增加，随骨水泥弹性模量的增加骨折椎体应力也随之增大，而邻近椎体应力无明显变化。提示降低骨水泥弹性模量可能是减少椎体成形后因骨折椎体及邻近椎体应力变化而发生骨折的一种方法。关键词：生物材料；骨生物材料；骨质疏松椎体压缩性骨

6、折；椎体成形术；骨水泥；弹性模量；骨折椎体；邻近椎体；生物力学；三维有限元分析3 P.O.Box 1200,Shenyang 110004 kf23385083Bao Yong-zheng, Studying for doctorate, Attending physician, Department of Orthopedics, Yue Bei Peoples Hospital, Shantou University Medical School, Shaoguan 512026, Guangdong Province, China主题词：脊柱骨折；骨质疏松性骨折；椎体成形术；生物力学；组

7、织工程基金资助：2013年韶关市医药卫生科研计划项目(Y13175)Effect of bone cement with a low elastic modulus on the fractured and adjacent vertebrae in patients with osteoporotic vertebral compression fractures: a three-dimensional finite analysisBao Yong-zheng1, Zhu Zhou-xing2, Feng Yun-sheng3, Wu Qiang1, Hu Kong-he1, He Xi

8、ao-long1, Zhu Wen-gang1, Xi Xin-hua1, Zhong Xue-ren1, Zhou Long-ze1, Liao Jun-jian1, Dai Xiang-heng1 (1Department of Orthopedics, Yue Bei Peoples Hospital, Shantou University Medical School, Shaoguan 512026, Guangdong Province, China; 2Department of Orthopedics, Lechang City Peoples Hospital, Shaogu

9、an 512200, Guangdong Province, China; 3Nanxiong Hospital of Traditional Chinese Medicine, Shaoguan 512400, Guangdong Province, China)AbstractBACKGROUND: Vertebroplasty and kyphoplasty can effectively repair osteoporotic vertebral compression fractures, but bone cement injection can cause the change

10、of stress in the fractured vertebrae and adjacent vertebrae after surgery, leading to new fractures.OBJECTIVE: To analyze the stress changes of the fractured vertebrae and adjacent vertebrae after vertebroplasty with different elastic modulus bone cement by a three-dimensional finite element method.

11、METHODS: One healthy adult male volunteer was selected for lumbar spine CT scan. The acquired images were imported for three-dimensional reconstruction using Mimics. The three-dimensional model was smoothed, polished and denoised by Geomagic software, and then the solid mode was built under Workbenc

12、h Ansys. An osteoporotic vertebral compression fracture model in L2-4 segments was established after assignment. Bone cement (4 mL) with different elastic moduli (8 000, 4 000, 2 000 and 1 000 MPa) injected into the L3 segment distributed in the middle of the vertebrae as spherical shape. 500 N pre-

13、load was applied on the L2 surface with an additional bending moment of 50 Nm. The lower surface free degree of L4 was restrained. The L2-4 forward flexion, posterior extension, right flexion and axial rotation on the right side were stimulated. The stress changes of the fractured vertebrae and the

14、upper and lower adjacent vertebrae before and after bone cement injection with different elastic moduli were compared. RESULTS AND CONCLUSION: The stress of the fractured vertebrae and adjacent vertebrae were significantly increased compared with that before operation. With the increase of elastic m

15、odulus, the stress of the fractured vertebrae increased, but there were no changes in the stress of adjacent vertebrae. These findings indicate that the elastic modulus of bone cement may be a method to reduce new fractures of the fractured and adjacent vertebrae after bone cement injection.Subject

16、headings: Spinal Fractures; Osteoporotic Fractures; Vertebroplasty; Biomechanics; Tissue EngineeringFunding: the Shaoguan Medical and Health Research Project, No. Y13175Cite this article: Bao YZ, Zhu ZX, Feng YS, Wu Q, Hu KH, He XL, Zhu WG, Xi XH, Zhong XR, Zhou LZ, Liao JJ, Dai XH. Effect of bone c

17、ement with a low elastic modulus on the fractured and adjacent vertebrae in patients with osteoporotic vertebral compression fractures: a three-dimensional finite analysis. Zhongguo Zuzhi Gongcheng Yanjiu. 2016;20(16):2285-2293.2289ISSN 2095-4344 CN 21-1581/R CODEN: ZLKHAH0 引言 Introduction骨质疏松症是以骨量减

18、少和微结构的退化(指由于骨组织的吸收和形成失衡等原因所致，表现为骨小梁结构破坏、变细和断裂)为特征伴骨脆性增加和易发骨折(指对负荷承受能力降低而易于发生微细骨折或完全骨折)的一种全身性骨骼疾病。骨质疏松是一种进行性并且难于逆转的病理生理过程，发病后很难再完全恢复至其正常骨组织的结构和质量。目前全世界约有2亿骨质疏松患者，骨质疏松椎体压缩性骨折是老年骨质疏松患者最常见的骨折之一，严重影响患者日常活动及生活质量，并可导致身高降低、疼痛、后凸畸形、卧床、焦虑、抑郁甚至死亡等并发症。椎体成形术是目前临床上微创修复骨质疏松椎体压缩性骨折的常规方式，注入的骨水泥凝固后即可达到止痛、稳定骨折椎体的目的，第2

19、天即可下地活动，极大改善了患者生活质量并减轻了社会、家庭的经济负担。但研究发现骨水泥的弹性模量明显大于骨质疏松骨折的骨皮质及松质骨，注入骨水泥后使椎体的刚度及强度明显增加，易导致骨折椎体再发骨折、邻近椎体骨折等并发症。降低骨水泥弹性模量可能是减少术后引起骨折椎体及邻近椎体应力变化而发生骨折的一种方法。文章拟建立腰椎骨质疏松椎体压缩性骨折模型，通过模拟治疗前及椎体成形后骨折椎体及上下邻近椎体各运动状态，分析治疗前与注射不同弹性模量骨水泥后骨折椎体及上下邻近椎体的应力变化。1 对象和方法 Subjects and methods 1.1 设计 计算机模拟试验。1.2 时间及地点 于2015年5至1

20、0月在粤北人民医院骨科完成。1.3 材料 CT机：General Electrics 64层螺旋CT机。计算机工作站：Intel Core i7 2600 四核3.4 GHz，16 G DDR3内存，1 G GDDR5 显存，1 T硬盘，操作系统Windows 7 Home Premium(64bit)。医学图像处理软件：Mimics 10.01(医学三维重建软件，比利时Materialise公司)，Geomagic 12(逆向工程软件，美国Geomagic公司)，Ansys Workbench 14.0(通用有限元计算软件，美国Ansys公司)。1.4 对象 志愿者1名，男，33岁，身高17

21、4 cm，体质量75 kg，无先天性、后天性脊柱畸形；无脊柱肿瘤、结核、腰椎间盘突出、骨质增生等疾病；无腰椎手术史，志愿者对试验方案知情同意。排除骨质疏松者。1.5 方法螺旋CT扫描：采用General 64层螺旋CT对已经选定的对象进行腰椎螺旋扫描及断层图像处理，扫描参数：层厚0.625 mm，球管电流200 mA、电压120 kV。扫描时志愿者采取仰卧位静止不动，保持扫描断面与身体长轴垂直。CT图像处理及保存：在CT工作站中调整图灰度、增加对比度，对图像观察细节进行处理，得到清晰的骨窗断层图像，保存为Dicom格式刻录光盘保保存。CT图像处理及腰段三维图像的重建：将Dicom格式的图像数据

22、导入三维重建软件Mimics。在Mimics中逐层分割提取已选取的CT图像，去除骨骼周围软组织图像，把腰椎L2-L4段从背景中分割。得到处理后每一个断层的CT图像，然后重建出腰段的三维图像，对其进行Mask蒙皮处理之后，以STL格式文件导出。腰段椎体三维实体模型的建立和光滑处理：将前面得到的STL文件读入Geomagic软件，对三维模型进行光滑、打磨、去噪等一系列图像处理，生成三维图形IGES文件格式。Workbench前处理：得到的IGES格式文件导入到Workbench中进行布尔运算操作，得到皮质骨、松质骨、后部结构、纤维环、髓核等组织结构。网格划分：将装配好的实体模型，赋好材料参数(表1

23、)，设置好接触，然后进行合适的单元设定及划分网格，单元尺寸2.0-3.0 mm，韧带单元采用杆单元(Link10)，设定其为只受拉不受压，其余部件均采用Solid 187实体单元。接触设置：由于各部件之间不会发生分离，故接触均设置为绑定。模拟运算：设定在L3为骨折椎体，注入4 mL弹性模量分别为8 000，4 000，2 000及1 000 MPa的骨水泥，分为设为正常弹性模量组及低弹性模量组1，2，3，骨水泥在椎体中央呈球形分布。在L2上表面施加500 N预载荷，附加弯矩为50 Nm，约束L4下表面自由度。模拟前屈、后伸、右屈及右轴向旋转(因左右对称，为减少运算量，本研究仅模拟右屈、右旋状态

24、)4个方向L2-L4椎体应力变化情况。1.6 主要观察指标 观察三维有限元模型中骨折椎体及邻近椎体在椎体成形注射不同弹性模量骨水泥前后各运动状态下的应力变化。1.7 统计学分析 采用SPSS 18.0统计软件进行数据分析，对各组数据行正态性检测，然后对治疗前与不同弹性模量骨水泥注射后相同运动状态数据行配对t 检验，P 0.05为差异有显著性意义。2 结果 Results 2.1 模型生成结果 组装生成的L2-L4腰椎骨质疏松椎体压缩性骨折模型见图1，椎体成形后骨水泥分布示意图见图2(骨水泥体积为4 mL)。2.2 治疗前及不同弹性模量骨水泥注射后骨折椎体及邻近椎体应力变化 见表2。2.2.1

25、治疗前及不同弹性模量骨水泥注射后L2应力 见图3。治疗前：最大前屈应力0.747 MPa，最大后伸应力0.714 MPa，最大右屈应力0.817 MPa，最大右旋应力0.600 MPa。表1 有限元模型中各材料属性1-3Table 1 Material properties of finite element models图2 注射4 mL骨水泥椎体成形术后示意图Figure 2 Model of osteoporotic fractures after treated by vertebroplasty with 4 mL bone cement BA图1 组装完成的L2-L4骨质疏松模型F

26、igure 1 The osteoporotic model of L2-4 segments图注：图中A为正面观，B为侧面观。表2 骨水泥注入前后各椎体不同运动状态下的应力比较 (s, MPa)Table 2 Comparison of mean stress of L2-4 vertebrae before and after bone cement injection with different elastic moduli under different movements椎体状态治疗前应力治疗后应力正常弹性模量组(8 000 MPa)低弹性模量组1(4 000 MPa)低弹性模量组

27、2(2 000 MPa)低弹性模量组3(1 000 MPa)L2右屈0.6050.1673.3250.001a2.3260.002 a2.3260.003 a2.3270.002 a旋转0.5890.0093.8310.004 a3.8300.004 a3.8270.004 a3.8260.003 a前屈0.7320.0143.4410.001 a3.4410.002 a3.4410.002 a0.4420.002 a后伸0.6870.0323.6270.005 a3.6300.003 a3.6310.003 a3.6330.003 aL3右屈0.6220.0123.6370.106 a2.9

28、620.012 a2.4010.007 a1.6260.011 a旋转0.8980.0198.2080.015 a5.9320.012 a4.0300.013 a2.5140.009 a前屈0.3350.0273.1460.028 a2.1540.067 a1.3660.013 a0.7830.009 a后伸0.2350.0253.7140.060 a2.6110.014 a1.6210.015 a1.0190.012 aL4右屈0.4020.0103.5010.010 a3.4940.016 a3.5020.015 a3.4990.022 a旋转0.8680.0413.6030.014 a3

29、.6040.016a3.6090.015 a3.6040.018 a前屈0.3090.0194.8520.008 a4.8510.008 a4.8510.007 a4.8510.007 a后伸0.3190.0255.7090.013 a5.7100.011 a5.7110.009 a5.7090.011 a材料名称弹性模量(MPa)泊松比横截面积(mm2)骨折后皮质骨3 3600.3骨折后松质骨280.2骨水泥正常弹性模量组8 0000.41低弹性模量组14 0000.41低弹性模量组22 0000.41低弹性模量组31 0000.41后部结构2 3450.25终板6700.4髓核8.40.4

30、5纤维环8.40.45前纵韧带200.338后纵韧带700.320黄韧带500.360棘间韧带280.335棘上韧带280.335关节囊韧带1000.340横突间韧带500.310表注：与治疗前比较，aP 0.05。各弹性模量组椎体成形后各椎体运动状态下应力均较治疗前明显增高(P 0.05)，L3椎体应力在各运动状态下随骨水泥弹性模量的增加逐渐增大(P 0.05)。正常性模量组(8 000 MPa)：最大前屈应力3.442 MPa，最大后伸应力3.632 MPa，最大右屈应力2.329 MPa，最大右旋应力3.836 MPa。低弹性模量组1(4 000 MPa)：最大前屈应力3.443 MPa

31、，最大后伸应力3.633 MPa，最大右屈应力2.329 MPa，最大右旋应力3.834 MPa。治疗前前屈 后伸 右屈 右旋正常性模量组(8 000 MPa)低弹性模量组1(4 000 MPa)低弹性模量组2(2 000 MPa)低弹性模量组3(1 000 MPa)图3 治疗前及不同弹性模量骨水泥注射后L2应力云图Figure 3 Mises stress cloud of L2 vertebral body in flexion, extension, right lateral flexion and right axial rotation before and after bone

32、cement injection with different elastic moduli前屈 后伸 右屈 右旋治疗前正常性模量组(8 000 MPa)低弹性模量组1(4 000 MPa)低弹性模量组2(2 000 MPa)低弹性模量组3(1 000 MPa)图4 治疗前及不同弹性模量骨水泥注射后L3应力云图Figure 4 Mises stress cloud of L3 vertebral body in flexion, extension, right lateral flexion and right axial rotation before and after bone cem

33、ent injection with different elastic moduli前屈 后伸 右屈 右旋正常性模量组(8 000 MPa)低弹性模量组1(4 000 MPa)低弹性模量组2(2 000 MPa)低弹性模量组3(1 000 MPa)治疗前图5 治疗前及不同弹性模量骨水泥注射后L4应力云图Figure 5 Mises stress cloud of L4 vertebral body in flexion, extension, right lateral flexion and right axial rotation before and after bone cement

34、 injection with different elastic moduli低弹性模量组2(2 000 MPa)：最大前屈应力3.443 MPa，最大后伸应力3.634 MPa，最大右屈应力2.329 MPa，最大右旋应力3.832 MPa。低弹性模量组3(1 000 MPa)：最大前屈应力3.444 MPa，最大后伸应力3.636 MPa，最大右屈应力2.329 MPa，最大右旋应力3.828 MPa。2.2.2 治疗前及不同弹性模量骨水泥注射后L3应力 见图4。治疗前：最大前屈应力0.364 MPa，最大后伸应力0.262 MPa，最大右屈应力0.636 MPa，最大右旋应力0.918

35、 MPa。正常性模量组(8 000 MPa)：最大前屈应力3.167 MPa，最大后伸应力3.790 MPa，最大右屈应力1.796 MPa，最大右旋应力8.226 MPa。低弹性模量组1(4 000 MPa)：最大前屈应力2.163 MPa，最大后伸应力2.628 MPa，最大右屈应力2.974 MPa，最大右旋应力5.942 MPa。低弹性模量组2(2 000 MPa)：最大前屈应力1.379 MPa，最大后伸应力1.639 MPa，最大右屈应力2.410 MPa，最大右旋应力4.042 MPa。低弹性模量组3(1 000 MPa)：最大前屈应力0.792 MPa，最大后伸应力1.032

36、MPa，最大右屈应力1.638 MPa，最大右旋应力2.526 MPa。2.2.3 治疗前及不同弹性模量骨水泥注射后L4应力 见图5。治疗前：最大前屈应力0.326 MPa，最大后伸应力0.344 MPa，最大右屈应力0.416 MPa，最大右旋应力0.909 MPa。正常性模量组(8 000 MPa)：最大前屈应力4.862 MPa，最大后伸应力5.722 MPa，最大右屈应力3.512 MPa，最大右旋应力3.619 MPa。低弹性模量组1(4 000 MPa)：最大前屈应力4.861 MPa，最大后伸应力5.722 MPa，最大右屈应力3.512 MPa，最大右旋应力3.620 MPa。

37、低弹性模量组2(2 000 MPa)：最大前屈应力4.860 MPa，最大后伸应力5.722 MPa，最大右屈应力3.520 MPa，最大右旋应力3.622 MPa。低弹性模量组3(1 000 MPa)：最大前屈应力4.860 MPa，最大后伸应力5.723 MPa，最大右屈应力3.527 MPa，最大右旋应力3.624 MPa。3 讨论 Discussion自2002年经皮椎体成形术引进国内治疗骨质疏松椎体压缩性骨折以来，目前已广泛应用于治疗畸形不严重及无脊髓、神经受压症状的骨质疏松椎体压缩性骨折。与卧床保守治疗相比具有见效快、卧床时间短等优点，与椎弓根螺钉内固定手术相比具有创伤小、并发症少

38、等优点，减轻了患者痛苦及减少了卧床相关并发症。虽然目前研究也尚未明确是否为椎体成形术的应力改变是否会引起再骨折还只是骨质疏松自然病程一部分。但有学者随访发现部分患者术后易再出现再骨折与多种因素相关4-10，且椎体成形术后再骨折时间明显较保守治疗组节段再骨折时间短11-12，Tanigawa等13对194例患者随访发现椎体成形术后大约有1/3患者会发生再骨折，其中约2/3发生在邻近椎体。黄天霁等14对148例因为单节段椎体压缩骨折行椎体成形术或椎体后凸成形术的患者进行平均34.4个月的随访，其中再发椎体骨折29例(共42个椎体)，邻近椎体骨折占55.17%，非邻近椎体骨折占44.83%。Chen

39、等15对1 800例患者共1 820个骨折椎体进行长达2年的随访，发现共有10例骨折椎体行椎体成形术后出现再骨折，均位于胸腰段，过于恢复高度以及填充的骨水泥呈块板分布是再骨折的原因。侯文根等16对180例患者行椎体后凸成形术后通过X射线计算骨水泥体积分数，发现单节段骨质疏松性椎体压缩性骨折行椎体后凸成形术患者中邻近椎体骨折发生率及骨水泥渗透率随骨水泥分数的增加而增加。弹性模量是指材料在弹性变形阶段内，正应力和对应的正应变的比值。骨水泥的弹性模量明显大于骨质疏松的骨皮质及松质骨，注入骨水泥后使椎体的刚度及强度明显增加，进而可能导致再发骨折。目前三维有限元分析已广泛应用于骨质疏松椎体压缩性骨折的研

40、究17-27。Wilcox 28通过三维有限元分析对L2-L3脊柱单位注入不等量的骨水泥，发现邻近节段再骨折的风险与注入骨水泥的量成正比。Yan等29建立有限元模型，按注入骨水泥量占椎体容积分为两组(8.6%-11.4%，13.1%-22.8%)。结果椎体的强度和刚度随着骨水泥量的增加而增加，且强度的增加较刚度增加明显。有限元法模拟研究发现，椎体成形术后终板下的骨水泥具有相当于垂直柱的作用，由此减少了椎体终板的生理凹陷，邻近椎体终板向内凸约17%，使上下终板向外变形，髓核内压力升高(19%)，增加相邻椎体的应力，易导致相邻椎体骨折30-31。贾小林等32对191例单椎体骨折行椎体后凸成形术的患

41、者进行回顾性分析，发现局部后凸角度、椎体压缩程度和术后骨折高度恢复情况是相邻及间隔椎体继发骨折的共同危险因素，低体质量指数和骨水泥椎间隙的渗漏是邻近椎体的单独危险因素。更但也有学者研究得出相反的结论，王涛等33对1例L2椎体骨质疏松性骨折行椎体后凸成形术治疗并随访1年半的1例患者，根据治疗前后患者的腰椎CT资料，构建椎体生物力学的有限元模型，对第L1-L3椎体在前后屈曲、左右侧屈及旋转状态下进行应力载荷实验，观察L2及邻近L1、L3椎体终板、椎体骨小梁、椎体骨水泥所承受的应力改变情况。研究发现术后L2椎体终板、椎体骨水泥在前屈、后伸、左侧屈、右侧屈状态下的应力值均高于术前(P 0.05)，椎体

42、后凸成形术前后椎体骨小梁各种运动状态下的应力值差异无显著性意义，L1、L3椎体后凸成形术前后各种模拟状态下应力变化差异无显著性意义，提示椎体后凸成形术后邻近椎体不会出现应力变化。 此次研究发现骨折椎体应力较术前明显增加，且随着注入骨水泥弹性模量的减小逐渐降低。提示降低骨水泥弹性模量可能是减少术后因骨折椎体及邻近椎体应力变化而发生骨折的一种方法。Persson等34把不同类型的脂肪酸和三酰甘油加入到加入到丙烯酸骨水泥中，发现均可改变骨水泥的刚度，加入极少量的亚油酸( 0.05)，结合本课题组的另一部分研究及既往研究提示邻近椎体应力主要与骨水泥注入体积相关，随骨水泥注入量的增加而增大(P 0.05

43、)，骨水泥注入体积是引发再骨折的主要因素39-43。因此在临床中不应追求注入高剂量的骨水泥，小剂量的骨水泥也可获得良好的临床疗效，并可减少骨水泥渗漏及继发骨折的风险。 研究通过CT扫描建立腰段骨质疏松椎体压缩性骨折椎体成形术后三维有限元模型，发现骨折椎体及邻近椎体应力治疗后均较治疗前明显增加，骨折椎体随骨水泥弹性模量的增加承受的应力也随之增大，而邻近椎体应力无明显变化。实验结果提示使用低弹性模量可水泥可降低骨折椎体应力从而减少再骨折的发生，为研发新型骨水泥提供理论依据。虽然目前尚未有成熟的低弹性模量水泥供临床使用，但将来必定有更也的骨水泥材料服务于临床。致谢：感谢粤北人民医院CT室同事的帮助。

44、作者贡献：包拥政进行试验设计，祝周兴、冯云升、胡孔和选择志愿者及CT扫描，包拥政、何小龙、朱文刚、席新华建立三维有限元模型及有限元分析，包拥政、钟学仁、周龙泽进行统计学处理，包拥政、廖军健、戴相恒成文，吴强审校。利益冲突：所有作者共同认可文章无相关利益冲突。伦理问题：志愿者者对研究项目知情并签署知情同意书，并同意影像学资料用于研究用途，试验获得粤北人民医院伦理委员会批准。文章查重：文章出版前已经过CNKI反剽窃文献检测系统进行3次查重。文章外审：本刊实行双盲外审制度，文章经国内小同行外审专家审核，符合本刊发稿宗旨。作者声明：文章第一作者对研究和撰写的论文中出现的不端行为承担责任。论文中涉及的原

45、始图片、数据(包括计算机数据库)记录及样本已按照有关规定保存、分享和销毁，可接受核查。文章版权：文章出版前杂志已与全体作者授权人签署了版权相关协议。4 参考文献 References1 邱兴,杨圣,芦健民,等.骨质疏松骨折椎体及相邻椎体骨水泥注入前后的有限元分析J.中国组织工程研究,2013, 17(17):3041-3048.2 严伟洪,仇胥斌,杨惠林,等. 胸腰段骨质疏松三维有限元模型的建立及临床应用J. 吉林医学,2011,2(32): 1053-1056.3 李秋军. 骨质疏松脊椎有限元模型的建立及临床应用D.北京协和医学院,2009.4 张鹭,吴军,路锟,等. 经皮椎体成形术治疗老年

46、骨质疏松椎体压缩性骨折术后邻近椎体骨折的相关危险因素J. 中国老年学杂志,2015,35(16):4612-4614.5 张斌,尚咏. 经皮椎体后凸成形术后继发邻近椎体骨折的危险因素分析J. 临床骨科杂志,2014,17(3):249-252.6 李智斐,付拴虎,张家立,等. 椎体成形术后再骨折相关因素分析J. 中国脊柱脊髓杂志,2014,24(9):790-794.7 易志坚,曹家树,王茂林,等. 经皮椎体成形术后椎体再发骨折的危险因素分析J. 颈腰痛杂志,2013,6:474-476.8 罗斌,吴星火,邸方. 椎体后凸成形术后邻近椎体再发骨折发生率及相关危险因素的分析研究J. 创伤外科杂志

47、, 2013,15(1):23-27.9 刘建,俞雷钧,宋红浦,等. 经皮椎体成形术后非手术椎体骨折的危险因素分析J. 中国骨伤,2013,27(3):190-193.10 闫亮,昌震,贺宝荣,等. 椎体成形术后邻近椎体骨折的相关危险因素分析J. 陕西医学杂志, 2013,42(6):676-678.11 Yi X, Lu H, Tian F, et al. Recompression in new levels after percutaneous vertebroplasty and kyphoplasty compared with conservative treatment. Arc

48、h Orthop Trauma Surg. 2014;134(1):21-30. 12 徐德兴,陈昌礼,陈硕. 经皮椎体成形术后非手术椎体新发骨折的临床分析J. 华西医学, 2015,30(7):1217-1221.13 Tanigawa N, Kariya S, Komemushi A, et al. Percutaneous vertebroplasty for osteoporotic compression fractures: long-term evaluation of the technical and clinical outcomes. AJR Am J Roentgeno

49、l. 2011;196(6):1415-1418. 14 黄天霁,寇玉辉,殷晓峰,等. 椎体强化术后再发椎体骨折的临床特点和危险因素J. 北京大学学报(医学版), 2015, 2(47):237-241.15 Chen LH, Hsieh MK, Liao JC, et al. Repeated percutaneous vertebroplasty for refracture of cemented vertebrae. Arch Orthop Trauma Surg. 2011;131(7):927-933. 16 侯文根,孙晓辉,张超,等. 骨水泥体积分数对单节段骨质疏松性胸腰椎体压缩

50、骨折行经皮椎体成形术患者邻近椎体骨折的影响J. 新乡医学院学报, 2015,6(12):1119-1121.17 Bleiler C, Wagner A, Stadelmann VA, et al. Multiphasic modelling of bone-cement injection into vertebral cancellous bone. Int J Numer Method Biomed Eng. 2015;31(1):e02696. 18 Matsuura Y, Giambini H, Ogawa Y, et al. Specimen-specific nonlinear

51、finite element modeling to predict vertebrae fracture loads after vertebroplasty. Spine (Phila Pa 1976). 2014;39(22):E1291-1296. 19 Cho AR, Cho SB, Lee JH, et al. Effect of Augmentation Material Stiffness on Adjacent Vertebrae after Osteoporotic Vertebroplasty Using Finite Element Analysis with Diff

52、erent Loading Methods. Pain Physician. 2015;18(6):E1101-1110. 20 Imai K. Computed tomography-based finite element analysis to assess fracture risk and osteoporosis treatment. World J Exp Med. 2015;5(3):182-187. 21 Giambini H, Qin X, Dragomir-Daescu D, et al. Specimen-specific vertebral fracture modeling: a feas