可灌注椎弓根螺钉优化设计的有限元分析

可灌注椎弓根螺钉优化设计的有限元分析贾宗海1 李景龙2 路青林作者简介：贾宗海（1983），男，研究方向：创伤及脊柱外科通讯作者：路青林，男，教授，硕士研究生导师。Tel：（0531）89268696；13791120838；Email:luqinglin 张文强1 彭大勇1 单仕营1（1.山东大学附属千佛山医院 骨科，山东济南250014；2.山东大学岩土与结构工程中心，山东济南 250061）摘要 目的 探讨可灌注椎弓根螺钉的最优化中孔直径和侧孔分布。方法 根据GSS椎弓根螺钉的几何数据和可灌注椎弓根螺钉的设计思路，建立普通椎弓根螺钉和可灌注椎弓根螺钉的有限元模型，模拟悬臂梁弯曲、拉伸和扭转试验对可灌注椎弓根螺钉进行优化分析。结果 1.普通椎弓根螺钉和中孔直径为1.0mm、1.5mm、2.0mm、2.5mm可灌注椎弓根螺钉模型的屈服载荷分别为994N、1026N、1120N、1083N、1064N。中孔直径为1.5mm时模型的屈服载荷最大，为普通椎弓根螺钉的113%。2.中孔直径为1.5mm，侧孔为3、4、6孔时螺钉模型的屈服载荷分别为1046N、997N、1013N。三个螺钉模型所能承受的最大拉力为4950N、4510N、4675N；最大扭矩为10.3Nm、8.7Nm、9.9Nm。结论 中孔直径为1.5mm，侧孔为3、4、6孔的可灌注椎弓根螺钉的机械强度与普通椎弓根螺钉相比没有降低，可满足临床的需要。关键词：椎弓根；螺钉；优化；设计；有限元分析The optimum design of cannulated pedicle screw: three-dimensional finite element analysisJIA Zong-hai1, LI Jing-long2, LU Qing-lin1, ZHANG Wen-qiang1 , PENG Da-yong1, SHAN Shi-ying 1 (1.Department of Orthopaedic,Qianfoshan Hospital, ShanDong University, Jinan 250014, China; 2.Geotechnical and Structure Engineering Technique Research Center, ShanDong University, Jinan 250061, China)Abstract: Objective To ascertain the optimum diameter of central hole and the distribution of radial holes in cannulated pedicle screw. Method Set up the three-dimensional finite element models of simplified solid pedicle screw and cannulated pedicle screws according to the geometric datas of general spine system pedicle screw and the design idea of cannulated pedicle screw in computer. Then optimize the design of cannulated pedicle screw through simulating cantilever bending test, tensile test and torsion test. Results 1. In cantilever bending test, the yielding load of the finite element models of solid pedicle screw and cannulated pedicle screws with the diameter of central hole 1.0 mm, 1.5 mm, 2.0 mm and 2.5 mm was 994 N, 1026 N, 1120 N, 1083 N, 1064 N respectively. The yielding load of cannulated pedicle screw with central hole 1.5 mm was 113% of that of solid pedicle screw. 2. The yielding load of the finite element models of cannulated pedicle screws with the diameter of central hole 1.5mm and radial holes 3, 4 or 6 was 1046 N, 997 N, 1013 N respectively. The maximum tensile force of three models was 4950 N, 4510 N, 4675 N respectively.The maximum torque was 10.3 Nm, 8.7 Nm, 9.9 Nm respectively. Conclusion The mechanical strength of cannulated pedicle screws with the diameter of central hole 1.5 mm and radial holes 3, 4 or 6 is equivalent to that of solid pedicle screw.Key word：Pedicle; Screw; Optimized; Design; Finite element analysis椎弓根螺钉内固定技术于1959年由Boucher【1】首次提出后，经过不断的改进与演变，目前该技术已经在临床上获得广泛应用。但对伴有骨质疏松的患者，因骨量减少,螺钉把持力降低而使螺钉松动、切割的发生率增加，限制了椎弓根螺钉（Solid Pedicle Screw, SPS）的应用。最近国内外有文献报道通过可灌注椎弓根螺钉(Cannulated Pedicle Screw, CPS)灌注骨水泥的方法增加螺钉在骨质疏松椎体中把持力【2-6】，但螺钉的规格尚没有统一的标准。本研究应用三维有限元法，从机械强度方面探讨可灌注椎弓根螺钉的最优化中孔直径和侧孔分布。1 可灌注椎弓根螺钉中孔直径的最优化1.1 椎弓根螺钉有限元模型的建立建立简化的通用型脊柱内固定系统（general spine system, GSS）椎弓根螺钉和中孔直径为1.0mm、1.5mm、2.0mm、2.5mm可灌注椎弓根螺钉的有限元模型(规格：6.5mm45mm)。根据GSS椎弓根螺钉的几何数据（螺钉尾端为14mm14mm11mm的“U”型结构，螺纹为自上而下外径一致的圆柱形，钉杆部分为锥形，尾端直径6.5mm，头端直径4.5mm），应用ANSYS 9.0建立两种五个螺钉的有限元模型，建模时忽略钉体螺纹结构(见图1)。定义模型为连续的钉体，采用4节点体Solid45单元划分网格，依中孔直径不同，每个模型约5 -15万个单元。假定各单元为均质、连续和各向同性的完全弹性体，且保持独自的介质特征。模型的材料力学属性见表1，赋值时将可灌注椎弓根螺钉中孔部分定义为骨水泥（聚甲基丙烯酸甲脂 Polylnethylmethacrylate，PMMA）的材料属性。两种材料均按弹性材料模拟，不考虑材料间的摩擦接触问题。图1 简化的GSS椎弓根螺钉有限元模型图及悬臂梁弯曲试验示意图Fig.1 The finite element model of simplified solid pedicle screw or schematic diagram of cantilever bending test表1 模型的材料力学属性Tab.1 Material properties for models of pedicle screw性能钛合金（Ti-6AL-4V）PMMA 弹性模量（Gpa）1104.0泊松比0.30.3屈服强度(Mpa)98038.70抗拉强度(Mpa)108044.90引自骨科基础科学【7】1.2 模型的边界条件和载荷模型左端长方体结构八个角点为全约束，右端锥形柱状体为自由端，其轴向方向为Z方向，指向头端为正；竖直方向为Y方向，向上为正。模拟两种五个螺钉模型的悬臂梁弯曲试验，在螺钉体部距离钉尖0-40mm范围内最上端的节点处施加荷载，方向竖直向下(如图1)。逐渐增大载荷值，直至模型局部达屈服强度。观察载荷作用下，模型的应力分布情况及局部达屈服强度时模型所能承受的最大总载荷，以此确定可灌注椎弓根螺钉中孔的最优化直径。2 可灌注椎弓根螺钉侧孔分布的分析2.1 建立不同侧孔位置的螺钉有限元模型在中孔直径最优化分析结果的基础上，建立侧孔夹角为120（3孔）、90（4孔）和120（3孔2组）三种螺钉的有限元模型。侧孔直径为1.5mm，侧孔位置见图2。其余条件同上，赋值时中孔、侧孔部分定义为骨水泥的材料属性。 图2 侧孔位置轴向示意图 孔(16)：第16侧孔的位置；数值(mm)：侧孔位置至螺钉头端的距离Fig.2 Axial schematic diagram of the location of radial holes hole (16):the location of the first to sixth radial hole; Num(mm):the distance from radial hole to the top of screw2.2 模型的边界条件和载荷约束条件同上，模拟3、4、6孔螺钉模型的悬臂梁弯曲、拉伸和扭转试验。悬臂梁弯曲模拟试验中加载方式同前。在拉伸模拟试验中为避免受力处局部应力集中，轴向的拉力载荷施加于模型头端平面的55个节点上（如图7a）。扭转模拟试验的加载是通过在模型端部的轴截面上、下两侧施加等值、反向的线荷载来模拟扭矩（如图8a）。逐渐增大作用于模型上的载荷，直至模型局部达到极限强度。观察在载荷作用下，模型的应力分布情况及局部达极限强度时模型所能承受的最大总载荷，以此分析侧孔分布对螺钉强度的影响。3 结果3.1 可灌注椎弓根螺钉中孔直径的最优化普通椎弓根螺钉和可灌注椎弓根螺钉的有限元模型及悬臂梁弯曲试验中螺钉模型的应力分布情况见图3。垂直载荷作用下，不同中孔直径螺钉模型的应力分布情况与普通椎弓根螺钉一致，最大应力均发生在模型根部，向头端方向应力逐渐减小，同时钉体还发生向下的位移。五个模型局部达屈服强度时所施加的最大总荷载值见图4，当中孔直径为1.5mm时，螺钉所能承受的总荷载最大，为普通椎弓根螺钉的113%。3a3b 3c 3d 3e3f图3 螺钉有限元模型及悬臂梁弯曲试验中模型的应力分布图(a、b为模型的x-平面图，cf为y-平面图)a:SPS位移等值线示意图；b:SPS的应力分布等值线图；cf:中孔直径1.0mm、1.5mm、2.0mm、2.5mm CPS等效应力达屈服强度时应力分布等值线图Fig.3 The finite element models of pedicle screws; von Mises stress distribution of pedicle screws in cantilever bending test (a and b:the x-plan of model of pedicle screw; cf:the y-plan) a:the isogram of displacement in solid pedicle screw; b:the isogram of von Mises stress distribution in solid pedicle screws; cf:the isogram of von Mises stress distribution in cannulated pedicle screws with the diameter of central hole 1.0mm, 1.5mm, 2.0mm, 2.5mm图4 螺钉模型悬臂梁弯曲试验最大总载荷（N）Fig.4 The yielding load of models of pedicle screws in cantilever bending test (N)3.2 可灌注椎弓根螺钉侧孔的分布通过最优化中孔直径的分析，建立中孔直径1.5mm，侧孔分别为3、4、6孔的螺钉有限元模型，螺钉内骨水泥模型见图5。图5 3、4、6孔CPS中孔及侧孔部位的骨水泥模型Fig.5 The finite element models of bone cement in cannulated pedicle screws悬臂梁弯曲模拟试验中，3、4、6孔可灌注椎弓根螺钉的有限元模型及载荷下模型的应力分布见图6。三种情况下，模型的最大应力均出现在模型根部，向头端方向应力逐渐减小。模型局部达屈服强度时所施加的总载荷分别为1046N、997N、1013N，见图9。6a6b6c图6 ac分别为3、4、6孔CPS等效应力达到屈服强度时应力等值线分布图Fig.6 The finite element models of cannulated pedicle screws with the diameter of central hole 1.5mm and radial holes 3, 4 or 6; ac:the isogram of von Mises stress distribution for cannulated pedicle screws in cantilever bending test拉伸模拟试验中,模型的受力模式及应力分布情况见图7。钉体的拉应力自左向右逐渐增大，最大应力出现在最右端侧孔的两侧。当模型局部达抗拉强度时，模型上施加的总拉力分别为4950N、4510N、4675N，见图9。7a7b7c7d图7 a:螺钉模型拉伸试验示意图；bd:3、4、6孔CPS拉应力达到抗拉强度时应力等值线分布图Fig.7 a:Schematic diagram of tensile test; bd:The isogram of von Mises stress distribution for cannulated pedicle screws in tensile test扭转模拟试验中，模型的受力模式及应力分布情况见图8。3孔和6孔螺钉模型的最大应力出现在距离头端最近的侧孔边缘处，4孔螺钉模型的最大应力出现在施加扭矩处。当局部达屈服强度时，3、4、6孔螺钉模型上施加的总扭矩分别为10.3Nm、8.7Nm和9.9Nm。8a8b8c8d图8 a:螺钉模型扭转试验示意图；bd:3、4、6孔CPS等效应力达到屈服强度时应力等值线分布图Fig.8 a:Schematic diagram of torsion test; bd:The isogram of von Mises stress distribution for cannulated pedicle screws in torsion test图9 3、4、6孔CPS模型的最大总载荷Fig.9 The maximum loads of cannulated pedicle screws in three simulated tests4 讨论老年骨质疏松患者合并腰椎退行性疾病、椎管狭窄、肿瘤、外伤等疾患需行后路椎弓根螺钉系统内固定时，因骨量减少螺钉把持力降低往往导致脊柱内固定失败。应用可灌注椎弓根螺钉可在置钉后经螺钉中孔、侧孔向椎体内注入骨水泥，既避免了普通椎弓根螺钉强化技术所导致的并发症，又可将骨-螺钉界面转化为骨水泥-骨界面，通过增大接触面积提高螺钉的把持力【3-5】。李书纲等【8】研究证实GSS螺钉较其他椎弓根螺钉固定系统有较好的椎体锚固作用。螺钉直径和长度的增加可以提高螺钉在骨质疏松椎体中的固定强度，但增加的程度受椎体和椎弓根解剖参数的限制，Strempel等【9】认为椎弓根螺钉的直径最大不应超过椎弓根横径的80%，而拧入深度应达进钉点至椎体前缘骨皮质距离的80%90%，席伟等【10】研究认为临床选择椎弓根螺钉的长度应不少于45mm。依据上述标准及国人腰椎椎弓根的形态学数据11 ，本研究在国产GSS椎弓根螺钉的基础上，选择直径6.5mm，长度45mm的椎弓根螺钉进行改良设计。中孔直径的大小直接关系螺钉的机械强度，杨述华等【6】对其设计的中孔直径为2.5mm的空心侧孔椎弓根螺钉与普通椎弓根螺钉进行三点折弯试验，结果显示空心侧孔椎弓根螺钉的三点抗折弯强度较普通椎弓根螺钉降低，但统计学比较无差异，可满足临床的需要。杨惠林等【4】对其设计的中孔直径2.0mm的可灌注椎弓根螺钉进行机械强度测试，得到同样的结论。侧孔的分布决定了骨水泥的弥散方向及骨水泥-骨的接触面积，对螺钉的固定强度其决定性作用【5】。目前文献中螺钉侧孔主要以90分布为主【2-5】，本研究假定中孔直径为1.0mm、1.5mm、2.0mm、2.5mm，并提出120（3孔）和120（3孔2组）两种侧孔分布方式，应用有限元法从机械强度方面探讨了可灌注椎弓根螺钉的最优化中孔直径和侧孔分布。三维有限元法为骨科临床内固定设计提供了一种有效的计算机模拟方法，其精确度能够满足临床的需要【12】。本研究应用三维有限元法，根据GSS椎弓根螺钉的几何数据及可灌注椎弓根螺钉的设计思路建立简化的螺钉模型。赋值时同时将中孔和侧孔部分定义为骨水泥的材料属性，这样更符合可灌注椎弓根螺钉应用时的实际情况。在螺钉中孔直径的优化分析中，模拟悬臂梁弯曲试验在螺钉体部施加垂直载荷，这一受力模式与椎弓根螺钉在人体中的受力情况相似，可以更合理的评价椎弓根螺钉的机械强度【13】。本研究结果显示垂直载荷下模型局部达屈服强度时，最大应力均发生在模型的根部，这与临床上螺钉断裂多发生在螺钉颈部相一致。通过本研究表明，中孔直径为1.0mm、1.5mm、2.0mm、2.5mm的可灌注椎弓根螺钉所能承受的最大弯曲载荷均较普通椎弓根螺钉增大，其中中孔直径为1.5mm时螺钉所能承受的荷载值最大，为普通椎弓根螺钉的113%。在侧孔分布的分析中，本研究建立了中孔直径为1.5mm，三种侧孔分布的可灌注椎弓根螺钉模型。悬臂梁弯曲模拟试验结果显示中孔直径1.5mm，三种侧孔分布的螺钉模型较不建立侧孔时模型的抗弯曲强度有所降低，但仍高于普通椎弓根螺钉，其中3孔和6孔螺钉模型的抗弯曲强度要优于4孔。在模拟拉伸试验中，3、4、6孔螺钉模型所能承受的最大拉力为4950N、4510N、4675N，足以满足临床的需要。在模拟扭转试验中，3孔和6孔螺钉模型的最大应力出现在距离头端最近的侧孔边缘处，4孔螺钉模型的最大应力出现在施加扭矩处。3、4、6孔螺钉模型所能承受的最大扭矩分别为10.3Nm、8.7Nm和9.9Nm，分析其原因可能因为侧孔之间的夹角不同造成最大应力出现的位置不同，但三种模型所承受的最大扭矩足以对抗螺钉在使用时承受的扭矩【4、8】。综上所述，本研究认为：1. 120（3孔）和120（3孔2组）两种侧孔分布螺钉模型的抗折弯、抗拉、抗扭转强度优于90（4孔）的螺钉模型；2. 中孔直径为1.5mm，侧孔为3、4、6孔的可灌注椎弓根螺钉在机械强度方面较普通椎弓根螺钉没有降低，足以承受上半身的重量，满足临床需要。本研究应用三维有限元法分析了可灌注椎弓根螺钉的最优化中孔直径及侧孔分布，但基于有限元法理论分析的局限性，该结论还有待于实体螺钉试验的进一步验证。参考文献1 Boucher HH. A method of spinal fusionJ. J Bone Joint Surg, 1959, 41B:248-259.2 Becker S, Chavanne A, Spitaler R, et al. Assessment of different screw augmentation techniques and screw designs in osteoporotic spinesJ. Eur Spine J, 2008,17:1462-1469.3 Yazu M, Kin A, Kosaka R, et al. Efficacy of novel-concept pedicle screw fixation augmented with calcium phosphate cement in the osteoporotic spineJ. J Orthop Sci, 2005,10:56-61.4 杨惠林, 王志荣, 王根林, 等. 可灌注骨水泥椎弓根螺钉的生物力学研究J. 中华骨科杂志, 2009, 29(3):241-247.5 Lih-Huei Chen, Ching-Lung Tai, Po-Liang Lai, et al. Pullout strength for cannulated pedicle screws with bone cement augmentation in severely osteoporotic bone: Influences of radial hole and pilot hole tappingJ. Clinical Biomechanics, 2009, 24:613-618.6 杨述华, 胡勇, 陈中海, 等. 空心侧孔椎弓根螺钉添加聚甲基丙烯酸甲酯骨水泥的生物力学研究J. 中华创伤杂志, 2002, 18 (1):17-22.7 Joseph AB, Thomas AE, Sheldon RS. 骨科基础科学M. 陈启明等, 译. 2版. 北京: 人民卫生出版社, 2001:172-180.8 李书纲, 邱贵兴, 翁习生, 等. 通用型脊柱内固定系统椎弓根螺钉翻修作用的生物力学研究J. 中华骨科杂志, 2002, 22(11):648-652.9 Strempel VA, Kuehle J, Plitz W. Stablilitaet von pedikelschrauben. Teil 2: maximale auszugskraefte unter beruecksichtigungder knochendichteJ. J Orthop, 1994,132(1):82-86.10 席伟, 雷伟, 严亚波. 椎弓根螺钉长度变化对螺钉-骨复合体模型应力影响的三维有限元分析研究J. 医用生物力学, 2010, 25(3):206-211.11 侯树勋, 史亚民. 国人下胸椎及腰椎椎弓根形态学特点及其临床意义J. 中华骨科杂志, 1994, 14:222-225.12 Lafage V, Gangnet N, Senegas J, et al. New interspinous implant evaluation using an in vitro biomechanical study combined with a finite-element analysisJ. Spine, 2007, 31(16): 1706-1713.13 Ching-Kong Chao, Ching-Chi Hsu, Jaw-Lin Wang, et al. Increasing bending strength and pullout strength in conical pedicle screw: Biomechanical tests and finite element analysesJ. J Spinal Disord Tech, 2008, 21:130-138.