用于评估手术辅助快速腭扩张模式的有限元分析模型

Summary

创建了一套新的手术辅助快速腭扩张 （SARPE） 有限元模型，该模型可以在不同角度的颊骨截骨术中执行临床所需的扩张器激活量，以进一步分析所有三个维度的半上颌骨扩张模式。

Abstract

引入手术辅助快速腭扩张 （SARPE） 以释放骨阻力，以促进骨骼成熟患者的骨骼扩张。然而，据报道，在所有 SARPE 患者中，有 7.52% 的患者出现左右两侧不对称扩张，其中 12.90% 的患者必须接受第二次手术进行矫正。导致不对称扩张的病因尚不清楚。有限元分析已被用于评估颌面结构中与SARPE相关的应力。然而，由于LeFort I截骨部位的骨骼碰撞仅在一定程度的扩张后发生，因此大多数现有模型并不能真正代表力分布，因为这些现有模型的扩张量很少超过1 mm。因此，有必要创建一个新的SARPE有限元模型，该模型可以执行临床所需的扩张激活量，以进一步分析所有三个维度的半上颌骨扩张模式。将锥形束计算机断层扫描 （CBCT） 的三维 （3D） 颅骨模型导入 Mimics 并转换为数学实体，以分割上颌复合体、上颌第一前磨牙和上颌第一磨牙。这些结构被转移到Geomagic中，用于表面平滑和松质骨和牙周韧带的创建。然后保留上颌骨的右半部分并镜像，以在 SolidWorks 中创建完美对称的模型。构建了 Haas 扩张器并将其绑扎到上颌第一前磨牙和第一磨牙上。在Ansys中对不同角度、间隙为1 mm的各种颊切骨术组合进行了有限元分析。进行收敛测试，直到达到所需的两侧膨胀量（总共至少 6 毫米）。本研究为评估颊骨切开成角如何影响SARPE的扩张模式奠定了基础。

Introduction

手术辅助快速腭扩张 （SARPE） 是一种常用的技术，用于骨骼成熟患者的上颌骨结构和牙弓横向扩张¹.该手术包括 LeFort I 截骨术、腭中皮质切开术，以及翼上颌裂² 的松解术。然而，已有报道称，SARPE 不希望的扩张模式，例如左半上颌骨和右半上颌骨³ 之间的不均匀扩张和牙槽突颊侧倾/旋转⁴，这可能导致 SARPE 失败，有时甚至需要额外的手术进行矫正⁵。先前的研究表明，上颌周截骨术的变化可能在 SARPE 扩张后模式中起重要作用^2,3，因为 Le Fort I 截骨部位的骨块之间的碰撞会导致半上颌骨横向扩张的不均匀阻力和半上颌骨的旋转，切口下方的牙槽边缘向内移动，而牙槽突扩张 ^3，4.因此，有必要研究不同截骨方向，尤其是颊骨截骨术对 SARPE 后扩张模式的影响。

已经建立了几个有限元分析（FEA）模型来评估SARPE过程中的力分布。然而，这些模型中的膨胀量限制在最多 1 mm，远低于所需的临床量 6,7,8,9,10,11,12。FEA 模型的扩展不足可能导致对 SARPE 后结果的错误预测。更具体地说，正如 Chamberland 和 Proffit⁴ 所报道的那样，如果扩张器没有充分转动，可能无法证明截骨部位的骨骼之间的碰撞，这可能不能反映真实的临床现实。由于先前模型中内置的扩展量有限，因此这些模型的结果评估集中在应力分析上。然而，牙科中有限元分析的应力分析通常是在静态载荷下进行的，材料的力学性能设置为各向同性和线性弹性，这进一步限制了有限元分析研究的临床相关性¹³。

此外，这些研究中的大多数没有考虑截骨部位⁶、^7、8、¹⁰、^11、12 处的手术器械厚度，通常将切口处的摩擦设置为零作为边界条件的一部分。然而，这种设置过度简化了硬组织和软组织之间的接触。它可能会显着影响力的分布和由此产生的半上颌骨扩张模式。

然而，没有可用的文献使用有限元分析 （FEA） 模型研究截骨术对 SARPE 后不对称性的影响。目前所有的研究都采用了对称截骨模式6,7,8,9,10,11,12,14的模型^，这并不能反映临床实践的现实，即截骨术在颅骨的每一侧可能不同。缺乏研究不对称截骨术对 SARPE 后不对称影响的文献，这是一个必须解决的重大知识差距。

因此，本研究的目的是开发一种新的 SARPE FEA 模型，该模型可以真正模拟临床状况，包括扩张量和截骨间隙，并研究各种截骨设计下半上颌骨在三个维度的扩张模式。这种方法将为 SARPE 后扩增模式背后的机制提供有价值的见解，并作为临床医生规划和执行 SARPE 程序的有用工具。

Protocol

本研究利用了一名患有 SARPE 的患者的预先存在的、去识别化的治疗前 CBCT 图像作为治疗计划的一部分。该研究是根据《赫尔辛基宣言》进行的，并得到了机构审查委员会（协议#853608）的批准。

1. 样本采集和牙齿分割

1. 在自然头部位置获取头部的人类 CBCT 图像，包括患者的上颌骨复合体，包括上颌基底骨、上颌牙槽骨和上颌牙列。
2. 将 CBCT Digital Imaging and Communications in Medicine （DICOM） 文件导入 Mimics 软件。
   1. 创建 新项目 （Ctrl + N），选择所有 DICOM 图像，然后单击 "下一步 "和 "转换"。
   2. 定义模型的方向（A:前方，P:后方，T:顶部，B:下部，L:左，R:右），然后单击 确定。
3. 将锉刀分割为上颌复合体、上颌第一前磨牙和上颌第一磨牙。
   1. 单击阈值，选择适当的 阈值来分割骨骼，然后单击 应用。
   2. 创建新面罩并单击" 编辑面罩"，使用 "绘制 和 擦除 "来分割患者的上颌骨复合体、上颌第一前磨牙和上颌第一磨牙。
4. 将目标导出为立体光刻 （STL） 文件。
   1. 右键单击蒙版，然后选择 "计算 3D "以生成 3D 对象。
   2. 右键单击 3D 对象，选择 STL+，选择所需的对象，然后按 "添加 "和 "完成 "以创建 STL 文件。

2. 表面平滑和形成松质骨和牙周韧带间隙

1. 将 STL 文件导入 Geomagic 软件。
   1. 单击 "文件">"打开"，选择 STL 文件，然后按 "打开"。
   2. 为单位中的数据弹出窗口选择"毫米"，然后单击"确定"。
2. 平滑上颌复合体、上颌第一前磨牙和上颌第一磨牙的表面。
   1. 单击 "多边形">"移除尖峰"，单击" 低"附近的平滑度级别并将其拖动，单击 "应用"，然后单击" 确定"。
   2. 单击 "多边形">"松弛多边形"，单击并拖动" 最小"附近的平滑度级别，单击 "应用"，然后单击" 确定"。
   3. 单击"多边形">"修复交叉点"，在"模式"窗口中选择"放松/清洁"，单击"应用"，然后单击"确定"。
3. 将模型的曲面修改为连续的闭合区域。
   1. 单击并拖动锐曲面，然后按 Delete 键以创建孔。
   2. 单击"填充孔">"多边形"，在"填充方法"窗口中使用"填充"、"填充部分"、"创建桥梁"来填充孔，然后单击"应用"，然后单击"确定"。
4. 将 2D 曲面转换为 3D 实体模型，并将其导出为计算机辅助设计 （CAD） 文件。
   1. 单击" 编辑">"相位">"形状相位"，选择 "编辑轮廓 "以草绘曲面的轮廓，然后单击" 确定"。
   2. 单击"绘制面片布局"（Draw Patch Layout） 并绘制四边形网格以覆盖所有曲面，然后单击"确定"（OK）。
   3. 单击 "构造网格"（Construct Grids），定义适当的"分辨率"（ Resolution），然后单击 "确定 "（OK） 以生成更精细的网格。
   4. 单击 "拟合曲面"（Fit Surfaces），单击"应用"（Apply），然后单击" 确定"（ OK ） 以构建 3D 实体模型。
   5. 单击 "文件"（File） >"另存为"（Save as ） 以输出 3D 模型并将其保存在 IGES 文件（名为 Maxilla）中。
5. 通过将上颌复合体的体积从颊牙槽表面减少 1 毫米来创建松质骨。通过将牙根的轮廓扩大 0.2 毫米来创造牙周韧带空间。
   1. 单击"面相位"（Polygon Phase），在"等值线"（Contour Lines） 窗口中选择"删除"（Delete），在"Patch 布局"窗口中选择"保留"（Preserve），然后按"确定"（OK） 将 3D 实体模型转换为 2D 曲面。
   2. 单击 多边形>偏移，在松质骨和牙周韧带的距离面板中输入 -1 mm 和 0.2 mm，然后单击 应用 和 确定。
   3. 单击 Edit > Phase > Shape Phase，选择 Restore Patch Layout，然后按 OK。
   4. 单击 "构造网格"（Construct Grids），定义适当的"分辨率"（ Resolution），然后单击 "确定 "（OK） 以生成更精细的网格。
   5. 单击 "拟合曲面"（Fit Surfaces），单击"应用"（Apply），然后单击" 确定"（ OK ） 以构建 3D 实体模型。
   6. 单击" 文件"（File） >"另存为"（Save as ） 以输出 3D 模型并将其保存在 IGES 文件（名为 CB 和 PL）中。

3. 构建解剖对称上颌骨模型

1. 将 CAD 文件输入到 SolidWorks 中。
   1. 单击 "文件">"打开"，选择上颌骨文件，然后按 "打开 "导入 CAD 文件。
   2. 单击 "文件">"保存 "将文件保存为 "部件 "格式。
2. 在腭平面 （PP） 下方构建松质骨。
   1. 单击" 插入>零件"，选择 CB 文件，然后按 "打开 "导入 CAD 文件。
   2. 单击" 插入>参照几何>平面"（Plane ），在腭面上选取三个特征点，然后单击 "确定 "（OK） 以创建切割平面。
   3. 单击 "插入>特征">"分割"，在 "修剪工具"中选择腭平面，然后单击"切割零件"以创建 切割 预览。
   4. 勾选 "生成的实体"中的复选框，然后单击" 确定 "以分离松质骨。
   5. 单击腭平面上方的松质骨，右键单击并按"正文"部分中的"删除"。
3. 构建上颌第一前磨牙和上颌第一磨牙的牙周韧带。
   1. 单击 "插入零件>"（Insert Part），选择 PL 文件，然后按" 打开 "（Open） 导入 CAD 文件。
   2. 单击" 插入">"特征">"相交"，然后在 "选择 "窗口中选择"上颌和PL"。
   3. 在"选择"窗口中选择"同时创建"，在"区域列表"中选择牙周韧带部分，然后单击"确定"以生成韧带。
4. 执行从前鼻棘 （ANS） 到后鼻棘 （PNS） 的腭中部切割平面，并保留上颌复合体的右半部分。
   1. 单击" 插入>参照几何">平面"（Insert Reference Geometry），在腭中平面上选取三个特征点，然后单击"确定"（ OK ） 创建切割平面。
   2. 单击" 插入">"特征">"分割"，在 "修剪工具"中选择腭平面，然后单击 "切割零件 "以创建切割预览。
   3. 勾选" 生成的实体"中的复选框，然后单击 "确定 "以分离上颌复合体。
   4. 单击上颌复合体的左半部分，单击鼠标右键，然后按"正文"部分中的"删除"。
5. 镜像上颌复合体的右半部分并创建一个相同的左半部分。
   1. 单击" > Mirror > Mirror"中插入图案/镜像，然后在 "镜像面/面"中选择中腭平面。
   2. 在 要镜像的身体中选择所有右半上颌复合体，然后单击 确定 以生成上颌复合体的左半部分。

4. 在上颌第一前磨牙和第一磨牙上创建一个 Haas 扩张器和束带

1. 构建前磨牙带和磨牙带。
   1. 单击 "插入零件"（Insert > Part），选择 PL 文件，然后按" 打开 "（Open） 导入 CAD 文件。
   2. 单击" 插入>特征">"分割"，选择 PL 文件中的齿，并将 "均匀缩放 "设置为 1.05。单击 "确定 "（OK） 生成厚度为 0.5 mm 的条带。
   3. 单击" 插入>参考几何>平面"，在遮挡平面上选择三个特征点，然后单击 "确定 "以创建参考平面。
   4. 单击" 插入>参照几何>平面"（Reference Geometry Plane），选择遮挡平面，并将偏移距离设置为 1.5 mm。单击 "确定 "（OK） 创建第一个切割平面。
   5. 单击 "插入>参照几何>平面"（Plane ），选择遮挡平面，并将偏移距离设置为 4.0 mm。单击 "确定 "（OK） 创建第二个切割平面。
   6. 单击" 插入">"特征">"分割"，然后在 "修剪工具"（Trim Tools ） 中选取第一个和第二个平面，并在 "目标主体"（Target Bodies） 中选取齿。单击 "切割主体"（Cut Bodies ） 以创建切割预览。
   7. 勾选" 生成的实体"中的复选框，然后单击 "确定 "以分离齿。
   8. 单击第一个平面上方和第二个平面下方的波段，单击鼠标右键，然后按"正文"部分中的"删除"。
2. 构建亚克力板。
   1. 单击" 插入>参照几何>平面"（Insert Reference Geometry），在硬腭平面上选取三个特征点，然后单击"确定"（ OK ） 创建草图平面。
   2. 单击 "插入">"草图"，绘制参考 Haas 扩展器的亚克力板，然后单击 "退出草图"。
   3. 单击" 插入">凸台/底座>"拉伸"，选择亚克力板的草图，在 "深度"中设置 5 mm，然后单击 "确定"。
   4. 单击 "插入>特征">"弯曲"，然后弯曲亚克力板以适应上颚的解剖结构。
   5. 单击 "插入>特征"（Insert Features > Fillet/Round"（圆角/圆形），然后在半径为 1 mm 的范围内对亚克力板的锋利边缘进行圆角。
3. 构造扩张臂。
   1. 单击" 插入>参照几何">平面"（Plane ），在波段上选取三个特征点，然后单击 "确定 "（OK） 以创建草图平面（命名为 P1）。
   2. 单击" 插入草绘">，绘制一个直径为 2 mm 的圆，然后单击" 退出草绘 "（名为 C1）。
   3. 单击" 插入>参照几何体>平面"，在亚克力板上选择三个特征点，然后单击 "确定 "以创建草图平面（命名为 P2）。
   4. 单击 "插入">"草绘"，绘制一个直径为 2 mm 的圆，然后单击" 退出草绘 "（名为 C2）。
   5. 单击" 插入>参照几何>平面"（Insert Reference Geometry），选取 P2 平面，并将偏移距离设置为 6 mm。单击 "确定 "（OK） 进入草绘平面。
   6. 单击" 插入草绘"（Sketch>），绘制一个直径为 2 mm 的圆，然后单击" 退出草绘"（Exit Sketch ）（名为 C3）。
   7. 单击" 插入">"凸台/底座">放样"，然后在 "截面梁"（Profiles ） 窗口中选取 C1、C2 和 C3 草图。
   8. 在"特征范围"窗口中选择波段和亚克力板，勾选"选项"窗口中的"合并结果"，然后单击"确定"。

5. 设计截骨术

1. 创建一个 1 毫米厚的平面，相当于外科医生通常使用的车针直径，从梨状孔 （Alar） 的角到颧下嵴 （IZC） 从水平面不同角度。
   1. 单击" 插入>参考几何>平面"，在截骨平面上选择三个特征点（与水平面成 0°、10°、20° 或 30°），然后单击" 确定 "以创建平面（命名为 O1）。
   2. 单击" 插入>参照几何>平面"（Insert Reference Geometry），选择截骨平面，并将偏移距离设置为 1.0 mm。单击 "确定 "（OK） 创建下部切割平面（命名为 O2）。
   3. 单击" 插入">"特征">"分割"，在 "修剪工具"中选择"O1"和"O2"平面，然后单击 "切割零件 "以创建切割预览。
   4. 勾选" 生成的实体"中的复选框，然后单击 "确定 "以分离上颌复合体。
   5. 单击 O1 和 O2 平面之间的主体，单击鼠标右键，然后在主体部分中按 Delete 键。
2. 在 Parasolid 模型零件文件 （X_T） 中导出具有不同颊骨截骨角度的模型以进行分析。
   1. 单击"文件">"另存为"，然后在"文件类型"列表中选择"Parasolid （x_t）"。
   2. 单击 "保存 "以导出有限元分析软件的模型。

6. 有限元分析

1. 将上颌骨复杂模型的材料参数导入并设置到Ansys软件中。
   1. 单击并拖动"工具箱中的静态结构"以创建分析工作空间。
   2. 双击工程数据，然后在属性中设置所有材料的杨氏模量和泊松比。不同结构¹²、^15、16的材料性能列于表1中。
   3. 双击 "几何"，单击 "文件">"导入外部几何文件"，然后单击 "生成 "以导入上颌复杂模型。
   4. 单击" 创建>布尔值"，然后用松质骨和牙齿生成皮质骨和牙周韧带。
2. 设置有限元分析模型。
   1. 双击 "模型"（Model），然后单击 "几何"（Geometry ） 以选择每个零件的材料属性。
   2. 右键单击"网格"（ Mesh ），然后单击 "生成网格"（Generate Mesh ） 以在模型上构建单元。
   3. 单击 "连接"（Connections），在" 接触主体 "（Contact Bodies） 中指定软/小零件，在 "目标实体"（Target Bodies） 中指定刚性/大零件。
   4. 在"定义"中指定接触类型和摩擦系数。表2列出了不同部件¹⁷的连接特性。
   5. 右键单击" 连接"，单击" 插入>弹簧 "以连接截骨平面的上部和下部。将弹簧设置为 1 mm 长，弹簧常数 k = 60 N/mm，并在每个网格节点处放置一个弹簧。
3. 在各种截骨组合上沿亚克力板上的x轴（垂直于中线）设置临床上可接受的力。
   1. 右键单击 "静态结构"（Static Structural），单击 "插入>固定支撑"（Fixed Support ），并将腭面上的结构设置为不可移动。
   2. 右键单击 "静态结构"（Static Structural），单击 "插入> 力"（Insert Force），然后设置一个 150 N 的力以施加在远离内线的方向的亚克力板上。
   3. 右键单击 "求解"（Solution），然后单击 "插入>变形"（Deformation > Total ） 以监视展开的变形。
4. 进行收敛测试，直到实现两端的扩展。
   1. 单击工具栏上的 "求解 "（Solve），然后等待"力收敛"（Force Convergence） 级别达到 "力标准"（Force Criterion）。
   2. 单击 "总变形"（Total Deformation ） 以显示展开结果。
5. 测量解剖标志在所有三个维度上的位移，作为扩张的结果。建议使用以下地标来评估扩展模式:
   上颌中切牙 （U1） 的近切线角。
   上颌第一前磨牙（U4）的颊尖。
   上颌第一磨牙（U6）的近中颊尖。
   梨状孔（Alar）的外侧下角。
   颧下嵴 （IZC）。
   扩张器的中点。

Results

该演示模型利用了一名 47 岁患有上颌骨缺陷的女性的 CBCT 图像。在生成的模型中，保留了鼻腔、上颌窦和扩张锚定牙齿（第一前磨牙和第一磨牙）的牙周韧带空间的解剖结构（图 1）。

为了准确模拟手术过程，在所有模拟中，鼻中隔、鼻腔侧壁和翼上颌裂与上颌体分离。此外，还创建了一个代表手术期间颊口截骨术的平面，厚度为 1 毫米。该平面从梨状孔（Alar）的角开始，向后延伸到翼状上颌裂（PMF）（图2A-D）。

对左右两侧均有对称零度切口的模型进行了初步测试（图2E），结果显示，150 N的力导致膨胀器的膨胀超过8 mm（图2F），超过了大多数文献中看到的膨胀量。这一结果被认为是合适的，因为它属于 SARPE 患者最常需要的扩张范围。...

Discussion

SARPE 中颊骨切开术的方向可以是上颌骨支撑区域向下踩踏之前从鼻孔水平切开，也可以是从梨状边缘向上颌第一磨牙相对应的支撑的斜坡切口，如 Betts² 所述。无论哪种方式，截骨术都延伸到上颌骨的颧突下方。然而，目前大多数关于SARPE的有限元分析研究都使用水平切口，向后延伸，与梨状边缘6,7,12,14相同。这偏离了临床上通常进行的操作，并改变了 FEA 中的条件，例如半上颌骨的质心以及截骨术的方向和接触区域。由于膨胀力并不总是穿过质心，因此在有限元分析期间，半上颌骨必然会发生旋转。然而，在临床情况下，截骨线处可能会发生碰撞，并且由此产生的旋转中心随后可能会发生变化。因此，为了产生临床上适用的结果，FEA 中的截骨术必须模仿现实生活中的手术模式。当前研究中引入的模型允许研究人员从不同角度构建截骨术（

Materials

List of materials used in this article

Name	Company	Catalog Number	Comments
Ansys	Ansys	版本 2019Ansys	是一款用于有限元分析的软件，可以基于微分方程求解复杂的模型。通过该软件分析不同颊侧截骨角度的扩展结果。
Geomagic Studio	3D Systems	版本 10Geomagic	Studio 是一款用于逆向工程的软件，可以根据物理扫描点生成数字模型。这项研究通过该软件构建了松质骨和牙周韧带。
Mimics	Materialise	版本 16Mimics	是一款基于医疗 3D 图像的工程软件，可高效地将 CT 图像转换为 3D 模型。本研究通过患者的 DICOM 图像重建了上颌骨复合体。
SolidWorks	Dassault Systèmes	2018 版	SolidWorks 是一款计算机辅助设计软件，供设计师和工程师创建 3D 模型。在本研究中，通过该软件设计和绘制了一台 Haas 膨胀机。