上颌埋伏尖牙正畸引过程的仿真及其生物力学量化分析.docxVIP

上颌埋伏尖牙正畸引过程的仿真及其生物力学量化分析

摘要

本研究旨在通过建立精确的上颌埋伏尖牙正畸牵引三维有限元模型，模拟正畸牵引过程，对其生物力学进行量化分析，探究不同牵引方案下牙齿移动过程中牙周组织及颌骨的应力应变分布规律，为临床制定个性化、高效的上颌埋伏尖牙正畸牵引治疗方案提供理论依据和技术支持。

关键词

上颌埋伏尖牙；正畸牵引；仿真；生物力学；量化分析

一、引言

上颌埋伏尖牙是口腔正畸临床中常见的错颌畸形问题，其发生率约为1%-3%。埋伏尖牙不仅影响患者的咀嚼功能和面部美观，还可能导致邻牙牙根吸收、牙列拥挤等一系列并发症。正畸牵引是治疗上颌埋伏尖牙的主要方法，但该过程复杂，受到多种因素影响，如牵引方向、牵引力大小、牙齿及牙周组织的解剖结构等。传统的治疗方案主要依赖医生的临床经验，缺乏精准的生物力学指导，可能导致治疗周期延长、治疗效果不佳等问题。随着计算机仿真技术和生物力学理论的不断发展，通过对正畸牵引过程进行仿真及生物力学量化分析，能够深入了解牙齿移动过程中的力学机制，为优化治疗方案提供科学依据。

二、上颌埋伏尖牙正畸牵引三维有限元模型的建立

（一）数据采集

选取一名上颌埋伏尖牙患者，采用CBCT（锥形束CT）进行口腔颌面部扫描，获取高分辨率的断层图像数据，扫描层厚设置为0.25mm。同时，使用3D口扫仪对患者的牙列进行扫描，获取牙冠表面的三维形态数据。将CBCT数据和3D口扫数据以标准格式导入医学图像处理软件Mimics中。

（二）模型分割与重建

在Mimics软件中，根据不同组织的灰度值差异，对CBCT图像进行阈值分割，分别提取上颌骨、牙槽骨、埋伏尖牙、邻牙、牙周膜等组织结构。通过区域增长、编辑蒙版等操作，去除噪声和伪影，获得精确的各组织结构的三维模型。将3D口扫获取的牙冠数据与CBCT重建的牙根数据进行融合，构建完整的牙齿三维模型。

（三）模型优化与网格划分

将Mimics中重建好的三维模型以STL格式导入三维建模软件Geomagic中，进行模型的平滑、修补等优化处理，去除模型表面的缺陷和不连续面。然后将优化后的模型导入有限元分析软件ANSYS中，采用四面体单元对模型进行网格划分，通过调整网格尺寸参数，确保模型的网格质量满足计算精度要求。在牙齿与牙周膜、牙周膜与牙槽骨等接触部位进行网格细化，以更准确地模拟力学传递过程。

（四）材料属性设定

根据相关文献资料，为模型中的各组织结构赋予相应的材料属性。牙齿被视为均质、各向同性的弹性材料，弹性模量为18.6GPa，泊松比为0.31；牙周膜为超弹性材料，采用Mooney-Rivlin模型描述其力学行为，材料参数根据实验数据确定；牙槽骨和上颌骨为多孔弹性材料，弹性模量分别设定为13.7GPa和18.0GPa，泊松比均为0.3。同时，定义牙齿与牙周膜、牙周膜与牙槽骨之间的接触关系为摩擦接触，摩擦系数取0.2。

三、上颌埋伏尖牙正畸牵引过程的仿真

（一）牵引方案设计

设计多种不同的正畸牵引方案，包括不同的牵引方向（如垂直牵引、斜向牵引）和牵引力大小（如2N、3N、4N）。在有限元模型中，通过在埋伏尖牙牙冠上施加相应方向和大小的力来模拟正畸牵引过程。

（二）边界条件设定

将上颌骨的底部和后部进行固定约束，模拟颌骨在实际生理状态下的固定情况。同时，限制模型在X、Y、Z三个方向上的平移和旋转自由度，以确保计算的准确性。

（三）求解与计算

在ANSYS软件中，采用隐式动力学分析方法对正畸牵引过程进行求解计算。设定计算时间步长和总计算时间，根据牙齿移动的实际情况，逐步模拟牙齿在牵引力作用下的移动过程。在计算过程中，记录每个时间步长下牙齿、牙周组织及颌骨的应力应变分布情况。

四、生物力学量化分析

（一）应力分布分析

分析不同牵引方案下，埋伏尖牙、邻牙、牙周膜及牙槽骨等组织结构的应力分布情况。重点关注牙齿根尖区、牙周膜颈部和牙槽嵴顶等应力集中区域，比较不同方案下应力峰值的大小和位置变化。研究发现，垂直牵引时，牙周膜颈部和牙槽嵴顶的应力集中较为明显；斜向牵引时，应力分布相对较为均匀，但在牙齿移动方向的对侧牙周膜区域会出现较高的应力值。随着牵引力的增大，各组织结构的应力峰值均呈现上升趋势。

（二）应变分布分析

对应变分布进行量化分析，了解牙齿移动过程中牙周组织的变形情况。通过计算不同部位的应变值，评估牙周组织的受力变形程度。结果表明，牙周膜的应变值较大，尤其是在牙齿移动方向上，这与牙周膜在正畸牙齿移动中起到的缓冲和传递应力的作用相符。同时，牙槽骨的应变主要集中在与牙周膜接触的区域，且应变值随牵引力大小和牵引方向的变化而改变。

（三）牙齿移动位移分析

测量埋伏尖牙在不同牵引方案下的移动位移，分析牵引