下颌角截骨整形手术三维有限元模型的建立及生物力学研究

摘要

1.研究背景
　　 目前下颌角截骨整形术主要以口内入路为主,传统手术方法是用摆动锯或来复锯进行截骨。由于口内入路空间小,手术暴露和用锯截骨时操作困难,锋利的锯片容易损伤血管和周围组织。因此柳大烈等提出应用钻凿法进行下颌角截骨,认为此法相比传统方法有更多优点,值得推广应用。钻凿法下颌角截骨整形术作为矫治下面部轮廓整形的重要术式受到很多研究者和临床医生的关注。
　　 2.目的和意义
　　 通过本课题研究,探索研究下颌角生物力学原理的方法和手段,建立高度仿真的下颌骨和颞下颌关节的三维有限元力学模型；通过对钻凿法下颌角截骨整形手术的仿真模拟,了解和阐述其生物力学原理:通过钻凿法下颌角截骨整形手术仿真模型,提出和分析对其的改进和优化的方法。
　　 3.方法
　　 3.1不同下颌骨三维有限元模型建模方法的比较
　　 3.1.1 Mimics直接建模法建立下颌骨三维有限元模型
　　 直接法为直接按照物体结构系统的几何外形建立节点和单元。基于女性青年志愿者头颅CT扫描图像,利用轮廓提取和阈值分割的方法提取下颌骨,删除附近其他骨骼和软组织。通过软件的三维重建和Remesh模块得到光滑连续的骨骼表面单元模型,导入Ansys软件中建立了下颌骨三维有限元模型。通过Mimics内部的FEA模块利用CT图像中骨组织内灰度信息,为每个体单元赋材料性质参数,模拟了骨组织材料性质非均匀的特性。建立了下颌骨三维有限元模型(模型-1)。
　　 3.1.2 Solidworks间接建模法建立下颌骨三维有限元模型
　　 间接建模法,需要建立物体的实体模型,再对实体模型进行自动网格划分,形成有限元模型。根据与上方法同源的原始CT图片,先在图像处理软件中进行图像分割把下颌骨分离出来,以点云形式输入到逆向工程软件Solidworks。软件对曲面域基于曲率划分,生成曲面修补,调整曲面修补的特征线消除扭曲曲面,导出曲面线框模型在Ansys软件中生成曲面,曲面围成的部分即自动生成所需的实体模型。有限元分析软件ANSYs对实体进行网格划分,生成下颌骨三维有限元模型(模型-2)。骨组织的材料参数与Mimics模型取相同值。
　　 3.1.3直接建模法和间接建模法比较
　　 为对两种建模方法进行比较分析,对Mimics直接建模法下颌骨三维有限元模型(模型-1)和Solidworks间接建模法下颌骨三维有限元模型(模型-2)采用相同的边界约束条件,在颏联合水平加载同等大小前后向冲击载荷,在同一求解器进行计算。比较两模型的Von Mises分布,对应部位应力大小进行统计学分析。
　　 3.2钻凿法下颌角截骨整形手术的生物力学研究
　　 3.2.1 带颞下颌关节的下颌骨三维有限元模型的建立
　　 采用女性青年自愿者头颅CT扫描图像,间接法建立包括颅骨,上颌骨,颞下颌关节盘和下颌骨三维有限元模型。杆单元模拟咀嚼肌约束,设定下颌骨和颞下颌关节窝与关节盘为3D面-面接触关系,约束关节窝的前、后和上表面个节点和杆单元上端所有方向的位移。建立包括颞下颌关节和咀嚼肌及韧带在内的下颌骨三维有限元模型(模型-3)。
　　 3.2.2钻凿法下颌角截骨整形手术模型的建立和生物力学分析
　　 根据临床手术设计和操作方法,通过施加布尔运算,对模型 3进行简化,钻孔切割等操作,模拟手术方案,建立钻凿法下颌角截骨整形手术模型(模型-4)。模拟手术操作在截骨线下缘施加载荷模拟截骨过程,根据失效准则判定失效载荷,观察分析截骨断裂过程中Von Mises应力的在分布下颌骨,颞下颌关节的分布。
　　 3.3钻凿法下颌角截骨整形手术的生物力学优化
　　 3.3.1 非均匀钻孔手术设计
　　 根据截骨线局部解剖测量和失效载荷最小化的原则,在模型-4上进行布尔运算,设计了下颌骨截骨线两端钻孔加密的新截骨方式,边界条件同模型-4,建立非均匀截骨线下颌角截骨整形手术模型(模型-5)。载荷条件和观察指标同模型-4。并将计算结果与模型-4进行对比。
　　 3.3.2下颌角颊侧皮质骨去除手术设计
　　 根据截骨线局部解剖测量和失效载荷最小化的原则,结合手术实际操作过程,设计了将截骨线上内侧下颌骨外板皮质骨去除的手术方法,对均匀截骨法下颌角截骨整形手术模型(模型-4)和非均匀截骨法下颌角截骨整形手术模型(模型-5)进行优化。边界条件、载荷条件和失效准则同模型 4 生成去除下颌骨外板皮质骨均匀截骨手术模型(模型-6)和去除下颌骨外板皮质骨非均匀截骨手术模型(模型-7)。观察指标同模型-4,并将模型4、5、6、7的计算结果进行对比,筛选出优化下颌角截骨手术的模型。
　　 3.3.3约束条件变化对下颌角截骨手术的影响
　　 在模型-7基础上模拟肌肉的杆单元上加载最大肌力的轴向预张力建立模型-8；约束手术侧截骨线上以上下颌支后缘部分节点建立模型-9;同时采用以上两种附加约束方法模型-10。并将模型7,8,9,10的计算结果进行比较。
　　 4.结果
　　 4.1 Mimics直接建模法建立下颌骨三维有限元模型
　　 建立了下颌骨三维有限元模型(模型-1),包含皮质骨和松质骨,共214354个节点,147219个四面体单元。
　　 4.2 Solidworks间接建模法建立下颌骨三维有限元模型
　　 建立了下颌骨三维有限元模型(模型-2),包含皮质骨和松质骨,共79538个节点,73999个四面体单元。
　　 4.3 Mimics建模法和Solidworks建模法比较
　　 两模型在冲击载荷下的响应、应力传导和分布趋势是一致的；在局部的最大应力值有差异；模型-1和模型-2在冲击载荷峰值时体部、磨牙区、髁突颈采样节点von Mises应力值数据,通过SPSS13.0软件进行统计学分析,配对t检验结果,t分别为-1.496,0.450,-0.286；P分别为0.23l,0.684,0.794；模型-1和模型-2之间的差异不显著(p>0.05)。
　　 4.4带颞下颌关节的下颌骨三维有限元模型的建立
　　 建立了包括颞下颌关节,颞下颌关节盘、咀嚼肌和韧带的下颌骨三维有限元模型(模型-3),包含1895460个单元,2863257个节点,24个杆单元,4个接触单元对。
　　 4.5钻凿法下颌角截骨整形手术的生物力学研究
　　 建立了均匀钻孔钻凿法下颌角截骨整形手术三维有限元模型(模型-4)；计算结果显示,下颌角断裂失效载荷为2426N；von Mises应力主要分布在截骨线下半段,同时向颏联合部传导；骨折线下半段和颏联合部都有明显的应力集中。手术侧关节盘应力较非手术侧大,大部分区域的应力值在0.56-7.54 MPa之间。
　　 4.6对钻凿法下颌角截骨整形手术的改进
　　 4.6.1 建立了下颌角非均匀钻孔的钻凿法下颌角截骨整形手术三维有限元模型(模型-5)；计算结果,模型-5失效断裂载荷为2018N。下颌骨von Mises应力分布显示,应力集中于截骨线下端并向颏部传导,最大应力点位于颏部,颏部红色显示的应力区域明显小于模型-4。
　　 4.6.2建立去除下颌角颊侧皮质骨均匀钻孔的钻凿法下颌角截骨整形手术三维有限元模型(模型-6)和去除下颌角颊侧皮质骨的非均匀钻孔的钻凿法下颌角截骨整形手术三维有限元模型(模型-7)；模型-6下颌角断裂失效载荷为1854N,模型-7为1447N；von Mises应力分布,模型-6和模型-7应力分布均较模型-4和模型-5明显集中,模型-7比模型-6应力分布更加集中,最大应力分布位于截骨线下半段,颏部无明显应力集中区域。
　　 4.6.3 计算结果模型-8失效载荷1183N；模型-9失效载荷1028N；模型-10失效载荷876N；模型-8 von Mises应力分布颏部的应力分布较模型-7明显减小；模型-9应力集中在截骨线及其上方区域内,而且较模型-7更加集中,方向与截骨线相同,下颌骨其他部位应力明显减小；模型-10应力集中更加明显,应力分布范围较模型-9更小。模型-8关节盘von Mises应力大部分在0.04-0.95 MPa之间,手术侧稍大,少数区域为0.95-8.2MPa；模型-9关节盘von Mises应力手术侧较非手术侧明显增大,关节盘上面局部地方von Mises应力在32.6-45.1MPa之间,后外侧边缘甚至超过了45.1MPa；模型-10关节盘von Mises应力手术侧大于非手术侧,大部分区域为0.52-3.25MPa之间,最大不超过9.86MPa。在模型-8中肌肉韧带的拉应力最大,模型-9最小,模型-10居中但仍然小于模型-7。
　　 4.7 存在的问题和今后研究展望由于人体组织结构的不规则、材料的非线性,而且大部分受力为动态,特别是颌骨组织,这对计算机提出了高性能的计算和制图能力要求。
　　 5.结论
　　 根据生物力学原理结合临床实践经验,利用仿真力学实验方法,对模型进行实验条件的仿真手术,对模型在不同条件下的截骨线局部破坏极限进行分析,以失效断裂载荷最小化为标准,寻求优化的截骨方法。通过对模型的计算结果分析,我们认为非均匀钻孔的方法优于均匀钻孔的方法,去除下颌角颊侧皮质骨的方法优于不去除下颌角颊侧皮质骨的方法；在非均匀钻孔的去下颌角颊侧皮质骨手术模型上,单独增加最大咬合力约束或下颌支后缘约束与同时增加最大咬合力和下颌支后缘约束都能明显减少失效载荷,但同时施加最大咬合力和下颌支后缘约束的方法更优。

目录

文摘

英文文摘

前言

第一章 下颌骨三维有限元模型的建立方法研究

第一节 Mimics建立下颌骨有限元模型

1 材料

2 方法

3 结果

4 讨论

第二节 solidworks建立下颌骨有限元模型

1 材料

2 方法

3 结果

4 讨论

第三节 两种建模方法的对比研究

1 材料

2 方法

3 结果

4 讨论

第四节 本章小结

参考文献

第二章 下颌角肥大钻凿法截骨整形手术的生物力学研究

第一节 包括颞下颌关节的下颌骨三维有限元模型的建立

1 颞下颌关节和咀嚼肌解剖及功能

2 材料

3 方法

4 结果

5 讨论

第二节 钻凿法下颌骨截骨整形手术的生物力学研究

1 材料

2 方法

3 结果

4 讨论

第三节 本章小结

参考文献

第三章 基于生物力学分析的下颌角钻凿法截骨手术的改进

第一节 非均匀钻孔设计的下颌角截骨整形手术模型

1 材料

2 方法

3 结果

4 讨论

第二节 去除下颌角颊侧皮质骨模型

1 材料

2 方法

3 结果

4 讨论

第三节 约束条件变化对下颌角截骨手术的影响

1 材料

2 方法

3 结果

4 讨论

第四节 本章小结

参考文献

全文总结

文献综述

参考文献

中英文缩略语对照表

攻读学位期间成果

致谢

统计学审稿证明