1.1 设计 胫骨假体设计及改进的有限元分析。
1.2 时间及地点 实验于2022年2月至2023年4月在上海工程技术大学材料科学与工程学院完成。
1.3 材料
1.3.1 患者资料 采集一名女性双侧膝骨关节炎志愿者的CT以及MRI影像数据,58岁,体质量62 kg,157 cm,双膝均无创伤史。扫描前受试者左右两侧膝关节的K-L等级基于门诊评估分别为Ⅰ级和Ⅲ级,右侧膝关节患有重度膝骨关节炎。
志愿者自愿参加,并且对实验过程完全知情,已签署“知情同意书”。此次研究经上海工程技术大学材料科学与工程学院医学伦理委员会批准。
1.3.2 实验仪器和软件 CT(Light Speed VCT;GE,USA);MRI(Signa 3.0T,GE,USA);Hypermesh 14.0(Altair Hyper-Works,USA);ABAQUS(Dassault,France);Mimics 19.0(Materialise,Belgium);3-Matic 11.0(Materialise,Belgium)。
1.4 方法
1.4.1 重度骨关节炎膝关节有限元模型的建立 该研究延续了张震等[13]的工作,采集一名女性志愿者的CT以及MRI影像数据,扫描时受试者身体保持仰卧位双腿呈伸直状。将扫描完成的双侧膝关节DICOM格式文件导入Mimics 19.0(Materialise,Belgium)中。通过识别影像数据中序列的坐标参数生成冠状位、轴状位和矢状位3个方向的视图。重度骨关节炎膝关节的影像图如图1所示。
对已提取出的骨性组织蒙版进行逐层重新分割、孔洞填补以及光顺处理,修复后的蒙版在兼顾骨关节炎膝关节骨赘特征的同时确保了骨性组织外形轮廓的完整性以及层与层之间的连续性。所建立三维解剖模型,如图2A所示。
根据受试者右侧重度骨关节炎膝关节股骨远端和胫骨近端在CT影像中反映的外形大小以及关节间隙,选用LEGION后稳定型骨水泥假体进行模拟置入,其中股骨髁大小选用3号,胫骨平台大小选用2号,垫片厚度选用9 mm。用计算机断层扫描仪对置入假体的三维几何结构进行扫描,利用Geomagic Studio进行模型的锐化处理和3-Matic构建假体的三维数值模型,如图2B所示,假体由徐州医科大学附属连云港医院提供。
全膝关节置换术中截骨的关键点在于保证膝关节0°和90°屈伸时关节间隙的平衡,此文基于右侧重度骨关节炎膝关节三维解剖模型,在Mimics中以胫骨侧为基准依次完成胫骨近端、股骨远端、股骨髁前、髁后以及髁中间部分的截骨,再次通过3-Matic实现骨与关节假体的全膝关节置换装配,如图2C所示。
1.4.2 假体改进 全膝关节置换的主要失效形式表现为感染、松动、髌骨脱位、胫骨平台断裂和疲劳破裂几种类型,其中无菌性松动导致失效约占总失效案例的25%;其次为感染,约占10%;而发生无菌性松动最主要的原因是骨吸收造成的骨量不足[4]。
全膝关节置换后膝关节运动承载时,弹性模量更高的人工膝关节假体负担较多的载荷,而弹性模量较低的骨骼负载降低,相当于假体遮挡了骨骼部分应力。根据Wolf定律,由于一个弹性模量更高的结构置入,没有足够的负荷转移到骨上时,应力屏蔽和骨吸收就会发生[14]。
通过对胫骨假体基座进行改进,使更多的负载能够传到骨骼上,减少应力遮挡效应。参照文献[4]四桩设计的金属基座,适当增加桩的长度、增加桩的直径和将桩外移都有利于增加固定的旋转稳定性。所以对单柄髓内柄假体(图3A)做出如下改进:①将单柄髓内柄假体的单柄和侧翼部分通过拉伸、切除命令去除,只留下胫骨平台;在胫骨平台基础上通过编辑草图、拉伸凸台命令拉伸出长为10 mm,直径为9 mm的四桩胫骨髓内柄,对胫骨柄进行2 mm的倒圆角处理,如图3B所示,即为四桩胫骨髓内柄假体。②在四桩胫骨髓内柄假体基础上,通过绘制草图、拉伸切除等命令对胫骨平台中部进行改进,拉伸切除出一个矩形孔洞,如图3C所示,即为有孔四桩胫骨髓内柄假体。
1.4.3 膝关节数值模拟计算
(1)网格划分:由于主要研究的是胫骨基座在膝关节屈曲时的应力分布与位移变化情况,在计算时只考虑股骨和股骨假体、胫骨和胫骨假体以及聚乙烯衬垫的约束固定关系,不考虑腓骨和髌骨中力的作用关系,也为了使有限元计算结果更准确,故可以去除髌骨和腓骨。将股骨、胫骨、股骨假体、胫骨假体、聚乙烯衬垫按照四面体实体单元(C3D4)网格划分,其中股骨划分为81 570个单元,股骨假体划分为49 561个单元,聚乙烯衬垫划分为77 066个单元,胫骨假体1划分为43 947个单元,胫骨假体2划分为56 374个单元,胫骨假体3划分为53 456个单元,胫骨1划分为90 698个单元,胫骨2划分为72 223个单元,胫骨3划分为71 361个单元。
(2)材料赋值:由于人工膝关节中的构成组件较多,且每个组件的材料特性均不同。在人工膝关节中,股骨假体和胫骨假体则采用线弹性假设,材料选用CoCrMo并赋值弹性模量和泊松比为220 000 MPa和0.3[3]。胫骨衬垫的材料通常选用高分子聚乙烯 [15]。由于日常行走、上下楼梯等活动时负荷频率大于0.1 Hz,骨骼变形仍处于弹性反应阶段,因此目前普遍认可的是对骨骼采用各向同性线弹性材料模拟。不同骨骼的材料参数不尽相同,参考文献[16]中骨骼的材料参数对模型进行材料定义。相关参数如表1所示。
(3)边界约束条件:膝关节假体模型存在接触,需要按照真实情况下关节间的相互作用进行约束,因为在此次模拟中忽略了关节软骨的作用,由于主要考虑安装了假体后的股骨和胫骨的应力分布和变形情况,对股骨假体进行约束,力的作用方向朝下,保证股骨假体在竖直方向移动,作用在衬垫上;对胫骨进行完全约束,限制X、Y、Z轴方向的位移约束,膝关节位移约束及载荷施加如图4所示。
在临床应用中,胫骨平台假体是与胫骨平台通过骨水泥固定,为了简化有限元模型,此次研究中将胫骨假体以及垫片设定为约束固定,限制其假体平移、旋转运动[17]。衬垫与股骨假体采用表面与表面接触,摩擦系数为0.04[18];衬垫与胫骨假体采用表面与表面接触,摩擦系数为0.14,对其他部件进行绑定接触设置,见表2,从而完成膝关节有限元模型的接触类型设置。
(4)载荷条件:对模型施加载荷,按照体质量62 kg进行加载,根据文献[18]按照胫骨轴线和股骨髁轴线方向施加不同倍数的人体体质量。当人体处于站姿时,力加载于股骨上,大小为人体质量;当假体处于屈膝30°时,股骨髁轴线方向上,大小为0.8倍人体质量;胫骨轴线方向上,大小为0.6倍人体质量;当假体屈膝60°时,股骨髁轴线方向上,大小为0.85倍人体质量;胫骨轴线方向上,大小为0.55倍人体质量;当假体90°屈膝时,处于跪姿,股骨髁轴线方向上,大小为0.9倍人体质量;胫骨轴线方向上,大小为0.5倍人体质量。
1.5 主要观察指标 在0°,30°,60°,90°角度下,观察胫股关节、胫骨、股骨、股骨假体、胫骨托和胫骨衬垫的Mises应力和变形云图。