3D打印假体与单纯钉棒系统用于全骶骨切除重建的有限元分析比较.docx

? ? 3D打印假体与单纯钉棒系统用于全骶骨切除重建的有限元分析比较 ? ? 张亚 马晓孟 孙允龙 李锋 全骶骨切除术是根治高位骶骨肿瘤或骶骨恶性肿瘤的重要手段，骶骨全切造成脊柱—骨盆完全分离，Galveston及其改良技术的出现使腰椎—骨盆的连续性得以重建，但仍然充满挑战性[1-2]。研究表明，16.1%～25%的全骶骨切除术患者发生内固定失败[3-5]，单纯钉棒固定存在着载荷完全由后路钉棒系统转移，无法维持结构平衡的缺陷。3D打印假体为腰椎—骨盆重建方式带来了新的可能，假体能够个性化定制，完美匹配骨缺损，一次性重建腰椎前柱和骨盆后环，降低操作难度。假体—骨界面的类骨小梁多孔结构也有助于促进骨长入，提高假体的稳定性。由于骶骨肿瘤发病率低，3D打印假体在全骶骨切除术中的临床应用案例非常少[6-9]，且多为个案报道，缺乏长期随访结果和生物力学分析，短期内无法为医生提供可靠的临床预后信息。本研究拟借助有限元分析方法，对骶骨全切后单侧双棒双髂骨钉重建方式与联合应用3D打印假体的两种重建方式进行评估，比较不同重建方式的生物力学性能，探索3D打印假体重建腰椎—骨盆稳定性的可行性。 1 材料与方法 1.1 模型设计 选取45岁成年女性的CT数据，经Mimics 20.0（比利时Materialise公司）软件处理后用于模型重建，保留L3、L4、L5椎体，L3/4、L4/5椎间盘（纤维环和髓核分别重建）及截骨后的双侧髋骨。使用3-Matic 11.0（比利时Materialise公司）软件设计三种腰椎－骨盆重建模型。模型1使用单侧双髂骨钉双棒重建；模型2使用单侧双髂骨钉双棒＋3D打印假体重建，其中髂骨钉穿过假体两翼，使用松质骨螺钉将假体与L5椎体及两侧髂骨连接；模型3也为单侧双髂骨钉双棒＋3D打印假体重建，但髂骨钉不穿过假体。腰椎椎弓根钉、松质骨螺钉及连接棒直径均为6.5 mm，髂骨钉直径7.5 mm，内固定钉棒及假体均为钛合金（Ti6Al4V）材料。将各模型部件保存为STL点云格式备用。 1.2 有限元模型的创建 使用Geomagic Studio 2014（美国Geomagic公司）软件将STL点云格式的模型部件逆向转化为曲面片文件，保存为IGES格式。使用SolidWorks 2020（法国Dassault公司）软件将IGES格式文件转化为STEP格式的实体文件（见图1A），导入Abaqus CAE（法国Dassault公司）有限元分析软件进一步处理。 1.3 材料属性赋值 模型各部分材料属性根据已有文献报道设定（见表1）。其中假体为3D打印钛合金（Ti6Al4V）多孔结构，孔隙率为70%～80%，孔径为600～800μm，与已获批上市的3D打印的钛合金内植物产品类似（孔隙率80%，孔径800μm），故参考相关公开的文献数据设定假体材料属性[10]。 表1 各模型部件材料属性 1.4 约束条件与载荷 内固定钉棒与骨骼和假体、椎体与椎间盘及纤维环与髓核间设定为绑定连接，忽略微动。假体两翼与髂骨截骨面之间定义为表面—表面接触，摩擦系数为0.3。限制双侧髋臼在三维坐标系上任意方向的移动，模拟双足站立姿势。在L3椎体上终板施加960 N的垂直载荷，模拟上半身的重力（见图1B）。 1.5 网格划分、装配、分析步与作业 使用四节点线性四面体单元（C3D4）为各部件划分网格并装配模型（见图1C），网格单元及节点数量如下（见表2），创建分析步并提交作业，应力分布情况以冯米斯应力云图展示，位移情况以位移云图展示。 图1 A.模型曲面实体；B.约束条件与载荷；C.网格划分 表2 各部件单元数及节点数 1.6 观察指标 （1）L5下沉距离：选取各模型L5下终板前缘10个点，查询其在垂直轴上的未缩放位移，取平均值作为L5下沉距离。 （2）髂骨钉平均最大应力：与L3椎弓根钉相连的两颗髂骨钉记为上位髂骨钉，与L4和L5椎弓根钉相连的两颗髂骨钉记为下位髂骨钉。在对称的两颗髂骨钉上选择应力集中区域最大应力的5个点，取平均值作为该髂骨钉平均最大应力。 （3）椎弓根钉平均最大应力：包括L3椎弓根钉、L4椎弓根钉和L5椎弓根钉。在对称的两颗椎弓根钉上选择应力集中区域最大应力的5个点，取平均值作为该椎弓根钉平均最大应力。 （4）锁定钉平均最大应力：在模型2与模型3中将连接假体和髂骨的松质骨螺钉记为锁定钉。选取锁定钉上应力集中区域最大应力的5个点，取平均值作为锁定钉平均最大应力。 （5）连接棒平均最大应力：与L3椎弓根钉相连的两根连接棒记为外侧连接棒，与L4和L5椎弓根钉相连的两根连接棒记为内侧连接棒。在对称的两根连接棒上选择应力集中区域最大应力的5个点，取平均值作为该连接棒平均最大应力。 比较三个模型L5下沉距离，髂骨钉、椎弓根钉和连接棒平均最大应力差异，同一模型内