空心螺钉联合同种异体腓骨治疗股骨颈骨折的三维有限元力学研究*

葛辉，苏丽霞，陈鹏，汤永南，罗东斌，焦锋，何伟

（1.广州市中西医结合医院 骨科，广东 广州 510800；2.广州中医药大学第一附属医院 骨科，广东 广州 510405；3.广州中医药大学第三附属医院 骨科，广东 广州 510378）

股骨颈骨折是临床上常见的骨折之一，约占全身骨折的3.6%，占股骨近段骨折的51.97%［1］。三枚空心螺钉内固定术是目前临床最常用的股骨颈骨折内固定方式［2-3］。尽管三枚空心螺钉内固定构型具有较强的抗剪力作用，使骨折断端能获得较好的力学稳定性，但股骨头坏死仍是其主要的并发症之一，根据研究显示：股骨颈骨折术后股骨头坏死发生率高达14.3%［4］。一旦发生股骨头坏死，股骨头将面临塌陷的风险，一旦塌陷将快速诱发髋关节骨性关节炎，最终不得不进行人工关节置换，这对于中青年股骨颈骨折患者将是灾难性的［5-11］。股骨颈骨折内固定术后如何预防股骨头坏死塌陷成为人们关注的重点。理想的股骨颈骨折内固定术式，是在股骨颈骨折复位后能提供有效的力学稳定性，同时在股骨头坏死发生后仍能形成稳定的力学支撑，预防股骨头塌陷，并促进股骨头坏死的修复。

相关研究证实腓骨在股骨头坏死后提供的有效力学支撑可以预防股骨头塌陷，同时腓骨的生物相容性还能促进坏死股骨头的修复［12］，因此腓骨移植支撑术常用于治疗非创伤性股骨头坏死。然而，股骨颈骨折术后继发股骨头坏死患者，如采用腓骨移植支撑术治疗股骨头坏死患者需面临二次手术，同时很多继发股骨头坏死患者诊断明确时股骨头已经塌陷，错过了最佳手术时机。如果在应用空心螺钉固定股骨颈骨折的同时联合腓骨植入，这样既解决了骨折固定的问题，还解决了术后股骨头坏死的防治问题，也避免了二次手术。腓骨没有螺纹，其力学强度也不及空心螺钉，空心螺钉联合同种异体腓骨内固定构型能否为股骨颈骨折的固定提供良好的力学稳定性，目前尚缺乏有效的力学依据。

为此本研究拟应用受试对象的CT数据，应用三维重建、逆向工程、有限元生物力学分析等技术，建立空心螺钉联合同种异体腓骨治疗股骨颈骨折的三维有限元力学模型，进行力学分析，评估其力学稳定性，为解决股骨颈骨折术后继发股骨头的治疗难题做出新的探索。

1 材料与方法

1.1 几何模型数据来源

1名45岁骨盆及股骨符合正常生理解剖的男性志愿者，体重70 kg，身高175 cm。应用Somatom Balance 64 排螺旋CT（德国西门子公司）获取研究对象的骨盆及股骨CT数据。

1.2 软件设备

Mimics 14.1软件（比利时Materialise公司），Geomagic Studio 2013（美国Geomagic公司），Geomagic Design X软件（美国Geomagic公司），Solidworks 2014（美国Dassault Systemes S.A），Abaqus 6.10（美国Dassault Systemes S.A），Rexcan CS+3D激光扫描仪（韩国首尔Solutionix公司），6.5 mm空心螺钉（Synthes公司）。

1.3 模型构建

1.3.1 股骨颈骨折三维模型的建立 将研究者的CT数据导入Mimics软件，通过Thresholding，Region Growing，Edit Masks和Calculating 3D等功能，重建出包括髂骨松质骨，髂骨皮质骨，股骨松质骨和股骨皮质骨的三维模型，然后利用Smoothing和Wrap等功能对三维模型进行初步优化处理，并导出STL格式文件。导入Geomagic Studio软件进一步行光滑、网格划分、轮廓探测及拟合曲面等处理，导出IGS格式文件。将文件导入Soildworks软件，应用拉伸切除功能，进行股骨颈骨折的模拟，建立股骨颈骨折的三维模型，本研究对骨折线垂直股骨颈的骨折类型进行建模，见图1。

图1 股骨颈骨折模拟示意图

1.3.2 空心螺钉固定股骨颈骨折的三维模型建立应用3D激光扫描仪扫描空心螺钉实体，建立空心螺钉的点云数据，保存为STL格式。将点云数据导入Geomagic Design X软件，进入“点云模块”，应用 “过滤杂点”、“平滑”等功能对点云数据进行优化，然后应用“构造曲面”功能将空心螺钉的点云数据制作成曲面模型。再应用Geomagic Design X软件的曲面片模块的“平滑”、“穴填补”、“面片优化”、“包裹”等功能对空心螺钉的曲面数据模型进行优化，制作成精细的曲面数据模型。然后应用“变换实体”把曲面数据模型转换成三维实体数据模型，见图2，保存为IGS格式。将建立的股骨颈骨折与空心螺钉三维模型导入到Solidworks软件中，进入装配模块，根据手术中空心螺钉的排列位置，进行空心螺钉固定股骨颈骨折的手术模拟，建立空心螺钉固定股骨颈骨折的三维实体模型，螺钉的排列方式为倒“品”字型，见图3，螺钉的植入深度为股骨头皮质骨下5 mm。

图2 空心螺钉三维实体模型

图3 三枚空心螺钉倒“品”字植入示意图

1.3.3 空心螺钉联合同种异体腓骨固定股骨颈骨折的三维模型建立 进入Soildworks软件零件模块，应用草图、凸台拉伸、拉伸切除等功能，建立腓骨的三维模型。根据腓骨的生理解剖，设置同种异腓骨的参数为长95 mm、外直径12 mm、内直径8 mm，见图4。将建立的股骨颈骨折、空心螺钉、腓骨三维模型导入到Solidworks软件中，进入装配模块，根据手术中空心螺钉和腓骨的排列位置，进行空心螺钉联合同种异体腓骨固定股骨颈骨折的手术模拟，建立空心螺钉联合同种异体腓骨固定股骨颈骨折的三维模型。螺钉和腓骨的排列方式为倒“品”字型，螺钉和腓骨的植入深度为股骨头皮质骨下5 mm，股骨头坏死常发生于股骨头的前上方，腓骨的植入位置位为倒“品”字型的左上方，见图5。

图4 同种异体腓骨三维实体模型

图5 空心螺钉联合同种异体腓骨倒“品”字植入示意图，腓骨位于前上方

1.3.4 关节软骨、韧带、肌肉三维模型的建立根据髋臼、股骨头软骨的厚度和关节间隙，应用Geomagic Design X软件建立髋臼、股骨头软骨三维模型。根据髋关节韧带的解剖学结构和生物力学特性，应ABAQUS软件桁架（TRUSS）模型模拟关节韧带。基于髋关节肌肉的解剖结构与生物力学特点，应ABAQUS软件的创建连接模块的轴（Axial）模拟髋关节7条主要肌肉（大收肌、梨状肌、臀小肌、臀中肌、长收肌、阔筋膜张肌、臀大肌）。根据SVERDLOVA等［13］报道的肌力和肌肉横断面面积关系，估算站立位关节周围7条肌肉的肌力大小见表1，在ABAQUS软件进行肌力赋值。

表1 髋关节7条主要肌肉的肌力 （N）

1.3.5 三维有限元模型的网格划分与材料属性赋值 进入ABAQUS有限元分析软件的“Mesh”模块，应 用“Seed Part Instance”、“Mesh Part Instance”等功能对上述建立的骨盆、股骨、韧带、关节软骨、腓骨、螺钉等组件进行网格划分。空心螺钉联合同种异体腓骨术式模型的单元总数为627 089，节点总数为161 869。多枚空心螺钉术式模型的单元总数为592 863，节点总数为160 798。

进入ABAQUS软件的“Property”模块，应用“Create Material”等功能，给上述的骨盆、股骨、韧带、软骨的三维模型，赋予各向同性匀质材料属性。具体赋值如表2所示［14-18］。

表2 有限元模型的材料属性

1.3.6 三维有限元模型的边界条件设置与力学加载 进入到ABAQUS软件的“Interaction”模块，“Create Interaction”、“Create Interaction Property”、“Create Constraint”等功能，将髂骨的骶髂关节和耻骨联合定义为完全约束；关节定义为无摩擦接触（Smooth）。韧带、肌肉与骨表面定义为捆绑（Tie）约束，植入物与宿主骨之间的接触定义为摩擦系数为0.6的摩擦接触，骨折断端的接触关系定义为摩擦系数0.46的摩擦接触，各接触关系见表3［19-22］。进入到“Load”模块，应用“Create Load”等功能，等效于体重的地面反作用力（700 N）施加于与股骨远端捆绑（Tie）在一起的刚体平面上，模拟人体单腿站立。已赋予边界条件的股骨颈骨折三维有限元模型，见图6。

表3 髋关节模型各个接触关系

图6 股骨颈骨折三维有限元模型

2 结果

2.1 股骨头颈部皮质骨的应力和位移情况

空心钉联合腓骨术式和多枚空心钉术式固定股骨颈骨折皮质骨最大应力均出现在股骨头皮质的头顶部，空心钉联合腓骨术式最大应力为81.38 MPa，多枚空心钉术式的最大应力为98.00 MPa，两者的差值为16.62 MPa，应力云图见图7。空心钉联合腓骨术式和多枚空心钉固定术式的股骨头颈位移均从股骨头颈上缘向下缘逐渐减小，最大位移均发生在股骨头颈张力最大的外上缘，空心钉联合腓骨术式最大位移为2.759 mm，多枚空心钉术式的最大位移为3.066 mm，两者的差值为0.307 mm，位移云图见图8。空心螺钉联合同种异体腓骨术式的股骨头颈部皮质骨的最大应力和最大位移分布位置与多枚空心螺钉术式差异无统计学意义，但股骨头颈部皮质骨的最大应力和最大位移值均比多枚空心螺钉术式小。

图7 股骨头颈部皮质骨应力云图

图8 股骨头颈部位移云图

2.2 骨组织内部和内固定的应力情况

应用剖面视图，观察股骨头内部骨组织的力学变化。空心螺钉联合同种异体腓骨术式和多枚空心螺钉术式固定股骨颈骨折的骨组织最大应力均出现在骨组织的皮质骨部位，空心螺钉联合同种异体腓骨术式最大应力52.76 MPa，多枚空心螺钉术式的最大应力65.09 MPa，两者的差值为12.33 MPa，应力云图见图9。空心螺钉联合同种异体腓骨术式内固定的最大应力为538.9 MPa，多枚空心螺钉术式内固定的最大应力为650.0 MPa，两者的差值为111.1 MPa，两者的应力较大区域均出现在骨折断端及内固定的尾部，应力云图见图10。空心螺钉联合同种异体腓骨术式的骨组织内部和内固定的最大应力分布位置与多枚空心螺钉术式无明显差异，但骨组织内部和内固定的最大应力值均比多枚空心螺钉术式小。

图9 骨组织内部应力云图

图10 内固定应力云图

3 讨论

稳定的股骨颈骨折内固定构型能有效的重建股骨颈内部的应力传导，避免应力集中，降低骨组织及内固定需要承受的最大应力。一旦发生应力集中，骨组织和螺钉承受的应力超过其屈服强度，将发生骨组织被切割和断钉等现象。

虽然空心螺钉联合同种异体腓骨术式与多枚空心螺钉术式的股骨头颈部皮质的最大应力和位移及骨组织内部和内固定的最大应力分布位置无明显差异，但空心螺钉联合同种异体腓骨术式的股骨头颈部皮质的最大应力和位移及骨组织内部和内固定的最大应力均比多枚空心螺钉术式小，说明空心螺钉联合同种异体腓骨内固定构型能更好地重建股骨近端的应力传导，使骨组织和内固定整体结构的应力相对分散，骨组织和内固定需要承受的最大应力更小，失效可能性低，具有更好的力学稳定性。

结合以上研究结果，本研究认为空心螺钉联合同种异体腓骨术式固定股骨颈骨折的力学稳定性优于多枚空心螺钉术式，能够为股骨颈骨折愈合提供良好的力学稳定性。同时同种异体腓骨具有力学支撑作用及生物相容性，股骨头坏死发生后能提供有效的力学支撑，预防股骨头的塌陷，腓骨的生物相容性能促进坏死的修复，因此空心螺钉联合同种异体腓骨术式为较理想的股骨颈骨折内固定术式。

本研究的优势：既往的髋部生物力学研究，进行髋关节三维有限元模型建立时，常常忽略关节软骨、韧带、肌肉对关节应力的影响；也未对骨盆进行建模，直接将力加载在股骨头上，无法真实的反应髋部的应力及股骨头的受力点［23］。本研究对上述问题进行了改良，使力学结果更加接近实际，对临床应用更具有指导意义。其次，以往对股骨颈骨折的力学研究着重于寻求更加稳定的内固定构型，避免内固定失效，忽略股骨颈骨折术后继发股骨头坏死的防治［24］。本研究应用腓骨替代一枚空心螺钉，建立空心螺钉联合腓骨的全新内固定构型，既保证了股骨颈骨折固定的力学稳定性，股骨头坏死发生后又能提供有效的力学支撑预防坏死股骨头塌陷，同时解决了股骨颈骨折固定和股骨头坏死的防治问题，本术式具有较好的临床推广意义。

本研究的局限性：本研究只对股骨颈骨折的一种常见类型进行了分析，未对股骨颈骨折的各种分型进行系统性的分析；同时，本研究在进行空心螺钉联合同种异体腓骨术式及多枚空心螺钉术式固定股骨颈骨折手术模拟时，假设股骨颈骨折已获得理想的解剖复位，尚缺乏对未获得解剖复位情况下生物力学稳定性的研究，对临床的参考尚存在一定的局限性，在进一步的研究中，笔者将对上述内容深入探讨。