需要JoVE订阅才能观看此内容。 请登录或开始免费试用

方法文章

一种用于定位上颌牙齿阻抗中心的有限元方法

9.3K 次观看

DOI:

10.3791/60746

2020年4月8日

本文内容

摘要

本研究概述了利用上颌及上颌牙齿的低剂量三维锥形束患者影像建立有限元模型所需的工具。随后,利用这些患者模型精确定位所有上颌牙齿的CRES

摘要

阻力中心(CRES) 被视为可预测牙齿移动的基本参考点。用于估算 C 的方法RES 牙齿的研究方法从传统的放射学和物理测量,到对模型或尸体标本进行的体外分析不等。基于模型和单颗牙齿的高剂量微计算机断层扫描(micro-CT)图像进行有限元分析的技术已展现出广阔的应用前景,但针对新型低剂量、低分辨率锥形束计算机断层扫描(CBCT)图像的研究仍较为有限。此外,CRES 仅有少数特定牙齿(如上颌中切牙、尖牙和第一磨牙)的特征被描述过,其余牙齿则在很大程度上被忽视。此外,还有必要描述确定C点的方法学RES 详细说明,以便于复制和进一步拓展。

本研究利用常规的锥形束CT(CBCT)患者影像,开发用于获取上颌牙齿阻抗中心(CRES)有限元模型的工具与工作流程。通过对CBCT容积图像进行分割操作,提取与确定上颌牙齿CRES相关的三维(3D)生物结构。使用3matic软件对分割出的结构进行清理,并将其转换为由四面体(tet4)三角面构成的虚拟网格,最大边长为1 mm。随后,将这些模型进一步转化为最大边长为1 mm的实体四面体体积网格,用于有限元分析。采用工程软件Abaqus对模型进行预处理,建立装配体,并设定材料属性、相互作用条件、边界条件及载荷施加方式。在分析过程中,施加的载荷可模拟系统中的应力与应变,从而辅助定位CRES。本研究是实现牙齿移动精确预测的第一步。

引言

牙齿或牙段的阻抗中心(CRES)类似于自由物体的质心,这一术语借用于刚体力学领域。当单一力作用于阻抗中心(CRES)时,牙齿将沿力的作用线方向发生平移运动1,2。阻抗中心(CRES)的位置不仅取决于牙齿的解剖结构和物理特性,还受其周围环境(例如牙周膜、邻近骨组织、邻牙)的影响。牙齿属于受约束的物体,因此其阻抗中心(CRES)与自由物体的质心类似。在矫治器的操作过程中,大多数正畸医生会考虑力矢量与单个牙齿或牙组阻抗中心(CRES)之间的关系。实际上,当物体受到单一力作用时,其是否发生倾斜移动或整体移动,主要取决于该物体阻抗中心(CRES)的位置以及力矢量与阻抗中心(CRES)之间的距离。若能准确预测这一关系,治疗效果将显著提升。因此,对阻抗中心(CRES)进行精确估算可极大提高正畸牙齿移动的效率。

几十年来,正畸领域一直在重新审视关于特定牙齿、牙段或牙弓的CRES位置的研究1,2,3,4,5,6,7,8,9,10,11,12。然而,这些研究在方法学上存在诸多局限性。大多数研究仅针对少数牙齿确定了CRES位置,而忽略了大多数牙齿。例如,上颌中切牙和上颌切牙区段已被广泛研究。相比之下,关于上颌尖牙和第一磨牙的研究较少,而对其余牙齿则尚无相关研究。此外,许多研究基于牙齿的通用解剖数据、二维(2D)X线片测量结果以及2D绘图的计算来确定CRES的位置8。另外,部分现有文献使用通用模型或牙模的三维(3D)扫描数据,而非真实的人体数据4,8。随着正畸学逐步转向采用3D技术进行牙齿移动规划,重新审视这一概念至关重要,以建立对牙齿移动的三维、科学理解。

随着技术进步带来计算能力和建模能力的提升,创建和研究更复杂模型的能力也随之增强。计算机断层扫描和锥形束计算机断层扫描(CBCT)的引入,使模型和计算从二维世界迈向三维。计算能力与软件复杂性的同步提升,使研究人员能够利用三维放射影像提取精确的解剖模型,并在高级软件中用于分割牙齿、骨骼、牙周膜(PDL)以及各种其他结构7,8,9,10,13,14,15。这些分割出的结构可转化为虚拟网格,供工程软件使用,以计算在施加特定力或位移时系统的响应。

本研究提出了一种特定且可重复的方法学,可用于分析作用于从活体患者CBCT图像提取的模型上的假设性正畸力系统。采用该方法学,研究人员可估算不同牙齿的阻力中心(CRES),并综合考虑牙体结构的生物学形态,例如牙齿解剖结构、根数及其在三维空间中的方向、质量分布以及牙周附着结构。该过程的总体流程如图1所示,旨在帮助读者理解生成三维牙体模型以定位阻力中心(CRES)所涉及的逻辑步骤。

访问受限。请登录或开始试用以查看此内容。

方案

已获得机构审查委员会豁免,用于评估口腔颌面放射学部存档的锥形束CT(CBCT)数据集(IRB编号:17-071S-2)。

1. 体积选择与标准

  1. 获取头面部的锥形束CT(CBCT)图像16
  2. 检查图像的牙齿排列情况、缺牙状况、体素大小、视野范围以及图像的整体质量。
  3. 确保体素大小不超过 350 µm(0.35 mm)。

2. 牙齿和骨骼的分割

  1. 将CBCT图像的原始DICOM文件导入Mimics软件进行分割(图2)。点击图像 > 裁剪项目,裁剪图像以仅包含上颌骨及上颌牙齿。
    注意:视野应足够大,以完整包含上颌骨和上颌牙齿。确保图像包括牙冠、硬腭至鼻底、上颌窦、上颌牙齿的面部表面,以及硬腭后缘和上颌结节。
  2. 右键单击掩膜标签,为图像创建一个新掩膜。根据目标牙齿,将掩膜重命名为左侧UL1、UL2、…、UL7,右侧UR1、UR2、…、UR7。
  3. 在掩膜后的CBCT图像上识别目标牙齿(参见各视图)。使用清除掩膜工具擦除掩膜。由于牙齿与骨组织的灰度值相近,软件可能无法区分两者。
    注意:Mimics中的阈值工具无法将牙齿和骨组织分别分割,因此需要采用其他分割方法。
  4. 点击多切片编辑工具(Ctrl + M),选择视图(轴向冠状矢状),根据需要手动高亮(即绘制)部分切片。
    注意:高亮更多切片可为结构增加更多细节。
  5. 点击插值工具,填充被跳过的切片之间的体积,并应用。
  6. 右键单击掩膜,选择计算3D体积的选项,生成该牙齿的3D体积。
  7. 对上颌牙弓中的每一颗牙齿重复步骤2.2–2.6。
  8. 选择所有上颌3D牙齿(UL7–UR7),右键选择平滑。将平滑因子设为0.4,迭代次数设为4。
  9. 为分割上颌骨,右键单击掩膜标签,为图像创建一个新的掩膜。
  10. 从预设阈值集的下拉菜单中选择自定义,调整阈值以包含完整的上颌骨。应用阈值前,请务必勾选填充孔洞选项。
    注意:皮质骨中≤1 mm的小孔可接受,因其可在后续阶段轻松去除。
  11. 点击动态区域生长工具,以填充掩膜中可见的大孔洞。除勾选多层选项外,还需选择上颌骨掩膜作为该工具的目标。将最小值设为50,最大值设为150。按住Control键的同时,点击掩膜中未被高亮的皮质骨区域。
  12. 右键单击上颌骨掩膜,使用平滑掩膜功能。为获得最佳效果,重复此步骤3次。
  13. 右键单击掩膜,选择计算3D体积的选项,生成上颌骨的3D体积。
  14. 选择上颌3D骨结构,右键选择平滑。将平滑因子设为约0.4,迭代次数设为4。
  15. 选择上颌3D骨结构,右键选择包裹。将最小细节设为0.2 mm,间隙闭合距离设为1 mm,并勾选保护薄壁选项。点击确定
  16. 将上颌3D骨结构重命名为"Maxilla"。

3. 清理与网格划分

  1. 选择3D对象并复制(Ctrl + C).
  2. 打开 3matic 软件,然后粘贴(Ctrl + V选定的3D对象将作为3D结构出现在3matic的对象树和工作区中图3).
  3. 点击 固定 从工具栏中选择标签并使用 平滑 选项。在 操作步骤 框选所需的 3D 对象或实体 应用 默认参数。
  4. 点击 完成 从工具栏中选择标签并使用 局部平滑 选项。在 操作步骤 框选所需的3D对象或实体。使用光标手动平滑所需区域。
  5. 复制牙齿。在对象树中选择所有牙齿,右键单击,然后选择 重复.
  6. 选择 所有重复牙齿,分组并命名文件夹 "第1组"原始样本将作为最终用于分析的牙齿。
  7. 对于第1组中的重复牙齿,单击 曲线模块创建曲线 选项。对所有重复的牙齿,手动绘制一条围绕釉牙骨质界(CEJ)的曲线。
  8. 选择 曲线, 轮廓,以及 边界 实体在以下 平滑曲线 选项
  9. 将冠部和根部表面分别选择并分离为独立的部分 按曲线分割曲面 选项并左键单击三维对象以选择。
  10. 通过在釉牙骨质界(CEJ)处将牙齿分为牙根和牙冠,从牙齿的牙根结构中获取牙周膜(PDL)。
    1. 将第1组的3D对象(在步骤3.6中生成)复制为第2组。对于第2组,在对象树框中,单击 对象从表面列表中删除牙冠表面。对第2组中的所有对象执行此步骤。
    2. 对于第2组,点击 设计模块 > 中空. 应用所需的参数(表1).
    3. 点击 固定模块 > Fix Wizard单击各个部件,更新并按照给出的指示进行操作。
    4. 对所有部件重复步骤 3.10.3。将组 2 中的所有部件重命名 "UL1_PDL" 到 "UL7_PDL" 和 "UR1_PDL" 到 "UR7_PDL".
  11. 在第1组中,从对象树框中单击 对象从表面列表中删除根面。
  12. 选择 填充孔洞(常规) 选项并选择轮廓。单击 不良轮廓应用整个空间将被填满。
  13. 选择 设计模块 > 局部偏移 并选择整个牙冠表面。勾选以下选项: 方向 (选择外部) 偏移距离 (选择 0.5),以及 缩短距离 (选择 2.0) 应用.
  14. 重复步骤 3.13。
  15. 对上颌牙弓的每一颗牙齿重复步骤 3.11–3.14。
  16. 重新网格化图3)
    1. 点击 重绘模块 > 创建非流形装配体 > 主要实体 > 上颌骨 从对象树中选择 3.4 节(原始牙齿)中所有对象的相交实体并选择 应用.
    2. 点击 Remesh 模块. 拆分非流形装配体。
    3. 使用来自第1组的所有对象作为相交实体,重复步骤3.16.1-3.16.2 应用.
    4. 作为可选步骤,仅在需要时选择 完成模块 > 修剪 > 实体 > 上颌骨. 选择多余的结构(即噪声)并 应用.
    5. 点击 固定模块 > 修复向导 > 上颌骨 > 更新按照给出的说明进行操作。
    6. 使用第2组中的所有对象作为相交实体,重复步骤3.16.1 应用.
    7. 点击 Remesh 模块 > 自适应重网格划分. 从 3.16.6 中选择所有相交的实体 应用.
    8. 点击 Remesh 模块 > 非流形分裂装配.
    9. 点击 重绘模块 > 创建非流形装配体 > 主要实体 > 单个物体 (PDL)从对象树的第2组中选择 相交实体 > 选择相应对象 来自步骤 3.4(对应于该牙齿类型)以及 应用.
    10. 点击 重网格模块 > 自适应重网格划分. 从 3.16.9 中选择相交的实体 应用.
    11. 点击 重绘模块 > 非流形分裂装配.
    12. 对每颗牙齿重复步骤 3.16.9–3.16.11。
  17. 点击 重绘模块 > 保持质量的三角面片简化在对象树中选择所有实体(即牙齿、牙周膜和上颌骨)并 应用.
  18. 点击 重绘模块 > 创建体积网格 > 选择实体选择 网格参数.
  19. 对所有实体(即牙齿、牙周膜和上颌骨)重复步骤 3.18。
  20. 手动将输入文件(.inp)从 3Matic 导出至 Abaqus图4).

4. 有限元分析

注意:所有自定义 Python 脚本均可在补充附件中找到。这些脚本是使用 Abaqus 中的宏管理器功能生成的。

  1. 预处理设置
    1. 打开 Abaqus 并选择 Standard Model。点击 File > Set the Work Directory > Select Location for File Storage
    2. 点击 File > Run Script,然后选择 Model_setup_Part1.py
    3. 在模型目录中指定加载 .inp 文件至 Abaqus 的文件路径。
    4. 点击 Models > Simulation > Parts > Maxilla > Surfaces
    5. 在对话框中将表面命名为 "UL1 _socket"。
    6. Select the Region of the Surface 下选择 By Angle,并输入角度值 "15"。
    7. 确保牙槽的所有区域均已被选中,完成后点击 Done
    8. 对各个牙槽重复步骤 4.1.4–4.1.7。
    9. 点击 Models > Simulation > Parts,然后选择 UL1 > Surfaces。将表面命名为 "UL1"。
    10. Select the Region of the Surface 下选择 "Individually",从屏幕上选择牙齿,然后点击 Done
    11. 对所有牙齿重复步骤 4.1.9–4.1.10。
    12. 点击 Models > Simulation > Parts,然后选择 UL1_PDL > Surfaces。将表面命名为 "UL1_PDL_inner"。
    13. Select the Region of the Surface 下选择 By Angle,并输入角度值 "15"。
      注意:若在最终模拟过程中发现错误,应减小角度并重新选择表面。
    14. 确保牙周膜(PDL)的整个内表面区域均被选中,完成后点击 Done
    15. 选择 UL1_PDL > Surfaces,将表面命名为 "UL1_PDL_outer"。
    16. Select the Region of the Surface 下选择 By Angle,并输入角度值 "15"。
      注意:若在最终模拟过程中发现错误,应减小角度并重新选择表面。
    17. 确保牙周膜(PDL)的整个外表面区域均被选中,完成后点击 Done
    18. 对所有牙周膜(PDL)重复步骤 4.1.13–4.1.19。
    19. 点击 File > Run Script,然后选择 Model_setup_Part2.py
    20. 点击 Models > Simulation > BCs。将边界条件命名为 BC_all,然后选择 StepInitial。在类别下选择 "Mechanical",在 "Types of Selected Step" 下选择 "Displacement/Rotation"。点击 Continue
    21. Select Regions for the Boundary Condition 下选择 By Angle,并输入角度值 "15"。勾选 Create Set,并为14颗牙齿分别选择各自的牙槽。完成后点击 Done
      注意:此操作有助于模拟牙齿的瞬时移动。
    22. 点击 Models > Simulation > Assembly > Sets > Create Set。将集合命名为 "U1_y_force"。
    23. Select the Nodes for the Set 下选择 Individually
      注意:在随机选择的牙冠节点上施加1牛顿的集中力,方向为正Y方向(模拟远中移动力)或正Z方向(模拟压低力)。
    24. 在上颌中切牙(U1)唇侧牙冠中心选择一个节点,然后点击 Done
    25. 点击 Sets > Create Set。将集合命名为 "U1_z_force"。
    26. 重复步骤 4.1.23–4.1.24。
    27. 对所有牙齿重复步骤 4.1.22–4.1.26。
      注意:在为特定牙齿生成集合(如步骤 4.1.25)之前,需先对该牙齿执行 Instance > Resume 操作。
  2. 模型设置
    1. 点击 Models > Simulation > Assembly > Instances。选择 All Instances,然后点击 Resume
    2. 点击 Tools > Query > Point/Node。在随机选择的中切牙中心选择一个节点,然后点击 Done
    3. 在页面底部的命令中心中,复制步骤 4.2.2 所选节点的 X、Y 和 Z 坐标。
    4. 在垂直工具栏中选择 Translate Instance,并在屏幕上选择整个装配体(即所有实例),然后点击 Done
    5. Select a Start Point for the Translation Vector 框中,粘贴步骤 4.2.3 复制的坐标,或手动输入 X、Y 和 Z 值,然后点击 Enter
    6. Select an End Point for the Translation Vector or enter X,Y,Z 中输入坐标 "0.0"、"0.0" 和 "0.0",然后点击 Enter
    7. 对于 Position of Instance,点击 Ok
    8. 点击 Tools > Query > Point/Node,选择位于中切牙中线正上方的一个节点,然后点击 Done
    9. 在页面底部的命令中心中,复制步骤 4.2.8 所选节点的 X、Y 和 Z 坐标。
    10. 在垂直工具栏中选择 Translate Instance,并在屏幕上选择整个装配体(即所有实例),然后点击 Done
    11. 将复制的坐标粘贴至 Select a Start Point for the Translation Vector - or Enter X,Y,Z 框中,然后点击 Enter
    12. Select an End Point for the Translation Vector - or enter X,Y,Z 中输入步骤 4.2.9 复制的坐标,并将 X 坐标修改为 0.0,然后点击 Enter
    13. 对于 Position of Instance,点击 Ok
    14. 点击 File > Run Script,然后选择 Model_setup_Part3.py。插入或修改材料属性。
    15. 点击 Models > Simulation > Materials,然后点击 Bone/PDL/Tooth。输入组织特异性属性。
    16. 点击 File > Run Script,然后选择 Functions.py
  3. 模型处理
    1. 点击 File > Run Script,然后选择 Job_submission.py
      注意:作业模块用于在模型上设置一个或多个操作,作业管理器则用于启动模型分析、显示进度并提示完成状态。
    2. 在标题为 Suppress All 的对话框中,根据约束条件输入牙齿的侧别(L 或 R)(在 Models > Simulation > Constraints 下)。点击 Ok
    3. 在标题为 Job Submission 的对话框中输入 "Y",以运行指定牙齿的分析。点击 Ok
    4. 在标题为 Directions for Analysis 的对话框中输入 "Y",以指定力的施加方向。点击 Ok
  4. 用于 CRES 估算的后处理
    1. 选择 File > Run Script > Bulk_process.py
    2. 在标题为 Analyze Multiple Jobs 的对话框中,为指定的牙齿输入 "Y"。点击 Ok
    3. 在标题为 Directions for Analysis 的对话框中,为指定力的施加方向输入 "Y"。点击 Ok
    4. 在标题为 Get Input 的对话框中,输入如命名实例所示的 Specific Tooth Number(例如 UL1 或 UL5 等)。点击 Ok
    5. 在命令框中核对 Force About PointEstimated Location 的坐标。若两者不一致,则重复步骤 4.3.1–4.4.4。
      注意:在运行每一步的作业后,将在 Abaqus 界面内运行一个用 Python 编写的用户自定义算法,用于分析反作用力系统及加载后产生的力矩。该算法会自动建议一个新的节点位置以施加载荷,使得力系统中的力矩接近零。该过程为迭代进行,直至找到或估算出施加力时能产生最接近零力矩的节点位置。该算法的详细说明见讨论部分。

访问受限。请登录或开始试用以查看此内容。

结果

为了验证操作步骤部分(步骤2)中所述的分割和手动勾画方法,从一个干燥颅骨中取出上颌第一磨牙,并拍摄了锥形束CT(CBCT)图像。按照步骤2所述,使用图像处理与编辑软件Mimics对该牙齿进行手动勾画。随后进行网格划分,利用3matic软件对分割后的模型进行清理,并将其导入Abaqus进行分析。我们未发现牙齿的有限元模型与实验室中实际测量的牙齿在线性及体积测量值之间存在显著差异(补充文件4)。

为验证用户自定义算法在确定物体CRES时的有效性,在脚本创建初期使用了一个梁被包裹于鞘内的简化模型(图5A)。钢制包层被约束为三个方向的位移自由度,梁与鞘界面处的节点被绑定在一起。随机选择施力节点,并以迭代方式应用子程序,直至解收敛。在该简化模型中,长度为30单位、宽度为10单位的梁被包裹于鞘内。通过遵循既定算法及其...

访问受限。请登录或开始试用以查看此内容。

讨论

本研究展示了一套工具,用于建立基于患者锥形束CT(CBCT)图像构建的上颌牙齿模型进行有限元分析(FEA)以确定其牙根表面中心(CRES)的一致性工作流程。对于临床医生而言,清晰且直观的上颌牙齿CRES分布图将是一种极为宝贵的临床工具,可用于规划牙齿移动并预测副作用。有限元法(FEM)于1973年引入牙科生物力学研究17,此后已被广泛应用于分析牙槽支持结构中的应力与应变场6,7,8,9,10,11,12。从工作流程中列出的步骤数量可以看出(图1<...

访问受限。请登录或开始试用以查看此内容。

披露

作者无任何利益冲突需要披露。

致谢

作者谨此感谢 Charles Burstone 基金会奖项对本项目的支持。

访问受限。请登录或开始试用以查看此内容。

材料

本文使用的材料清单
姓名公司目录编号评论
3-matic 软件Materialise, Leuven, Belgium.清洁与网格化处理
Abaqus/CAE 软件,版本 2017Dassault Systèmes Simulia Corp., Johnston, RI, USA.有限元分析
Mimics 软件,版本 17.0Materialise, Leuven, Belgium.牙齿与骨骼的分割

参考文献

  1. Smith, R. J., Burstone, C. J. Mechanics of tooth movement. American Journal of Orthodontics. 85 (4), 294-307 (1984).
  2. Christiansen, R. L., Burstone, C. J. Centers of rotation within the periodontal space. American Journal of Orthodontics. 55 (4), 353-369 (1969).
  3. Tanne, K., Nagataki, T., Inoue, Y., Sakuda, M., Burstone, C. J. Patterns of initial tooth displacements associated with various root lengths and alveolar bone heights. American Journal of Orthodontics and Dentofacial Orthopedics. 100 (1), 66-71 (1991).
  4. Burstone, C. J., Pryputniewicz, R. J. Holographic determination of centers of rotation produced by orthodontic forces. American Journal of Orthodontics. 77 (4), 396-409 (1980).
  5. Dermaut, L. R., Kleutghen, J. P., De Clerck, H. J. Experimental determination of the Cres of the upper first molar in a macerated, dry human skull submitted to horizontal headgear traction. American Journal of Orthodontics and Dentofacial Orthopedics. 90 (1), 29-36 (1986).
  6. Tanne, K., Sakuda, M., Burstone, C. J. Three-dimensional finite element analysis for stress in the periodontal tissue by orthodontic forces. American Journal of Orthodontics and Dentofacial Orthopedics. 92 (6), 499-505 (1987).
  7. Meyer, B. N., Chen, J., Katona, T. R. Does the Cres depend on the direction of tooth movement? American Journal of Orthodontics and Dentofacial Orthopedics. 137 (3), 354-361 (2010).
  8. Kojima, Y., Fukui, H. A finite element simulation of initial movement, orthodontic movement, and the centre of resistance of the maxillary teeth connected with an archwire. European Journal of Orthodontics. 36 (3), 255-261 (2014).
  9. Reimann, S., Keilig, L., Jäger, A., Bourauel, C. Biomechanical finite-element investigation of the position of the centre of resistance of the upper incisors. European Journal of Orthodontics. 29 (3), 219-224 (2007).
  10. Viecilli, R. F., Budiman, A., Burstone, C. J. Axes of resistance for tooth movement: Does the Cres exist in 3-dimensional space? American Journal of Orthodontics and Dentofacial Orthopedics. 143 (2), 163-172 (2013).
  11. Ammar, H. H., Ngan, P., Crout, R. J., Mucino, V. H., Mukdadi, O. M. Three-dimensional modeling and finite element analysis in treatment planning for orthodontic tooth movement. American Journal of Orthodontics and Dentofacial Orthopedics. 139 (1), 59-71 (2011).
  12. Sia, S., Koga, Y., Yoshida, N. Determining the center of resistance of maxillary anterior teeth subjected to retraction forces in sliding mechanics. An in vivo study. Angle Orthodontics. 77 (6), 999-1003 (2007).
  13. Cattaneo, P. M., Dalstra, M., Melsen, B. Moment-to-force ratio, center of rotation, and force level: a finite element study predicting their interdependency for simulated orthodontic loading regimens. American Journal of Orthodontics and Dentofacial Orthopedics. 133 (5), 681-689 (2008).
  14. Tominaga, J. Y., et al. Effect of play between bracket and archwire on anterior tooth movement in sliding mechanics: A three-dimensional finite element study. Journal of Dental Biomechanics. 3, 1758736012461269(2012).
  15. Cai, Y., Yang, X., He, B., Yao, J. Finite element method analysis of the periodontal ligament in mandibular canine movement with transparent tooth correction treatment. BMC Oral Health. 15 (106), (2015).
  16. Pauwels, R., Araki, K., Siewerdsen, J. H., Thongvigitmanee, S. S. Technical aspects of dental CBCT: state of the art. Dentomaxillofacial Radiology. 44 (1), 20140224(2015).
  17. Farah, J. W., Craig, R. G., Sikarskie, D. L. Photoelastic and finite element stress analysis of a restored axisymmetric first molar. Journal of Biomechanics. 6 (5), 511-520 (1973).
  18. van Driel, W. D., van Leeuwen, E. J., Von den Hoff, J. W., Maltha, J. C., Kuijpers-Jagtman, A. M. Time-dependent mechanical behavior of the periodontal ligament. Proceedings of the Institution of Mechanical Engineers, Part H: Journal of Engineering in Medicine. 214 (5), 497-504 (2000).
  19. Bourauel, C., et al. Simulation of orthodontic tooth movements. A comparison of numerical models. Journal of Orofacial Orthopedics. 60 (2), 136-151 (1999).
  20. Schneider, J., Geiger, M., Sander, F. G. Numerical experiments on longtime orthodontic tooth movement. American Journal of Orthodontics and Dentofacial Orthopedics. 121 (3), 257-265 (2002).
  21. Ten Cate, A. R. Oral histology, development, structure and function (5th ed). , St. Louis Mosby. (1998).
  22. McCormack, S. W., Witzel, U., Watson, P. J., Fagan, M. J., Gröning, F. The Biomechanical Function of Periodontal Ligament Fibres in Orthodontic Tooth Movement. PLoS One. 9 (7), e102387(2014).
  23. Huang, H., Tang, W., Yan, B., Wu, B., Cao, D. Mechanical responses of the periodontal ligament based on an exponential hyperelastic model: a combined experimental and finite element method. Computer Methods in Biomechanics and Biomedical Engineering. 19 (2), 188-198 (2016).
  24. Yang, J. A new device for measuring density of jaw bones. Dentomaxillofacial Radiology. 31 (5), 313-316 (2002).
  25. Gradl, R., et al. Mass density measurement of mineralized tissue with grating-based X-ray phase tomography. PLoS One. 11 (12), e01677979(2016).
  26. Jiang, F., Kula, K., Chen, J. Estimating the location of the center of resistance of canines. Angle Orthodontics. 86 (3), 365-371 (2016).
  27. Nyashin, Y., et al. Center of resistance and center of rotation of a tooth: experimental determination, computer simulation and the effect of tissue nonlinearity. Computer Methods in Biomechanics and Biomedical Engineering. 19, 229-239 (2016).
  28. Toms, S. R., Eberhardt, A. W. A nonlinear finite element analysis of the periodontal ligament under orthodontic tooth loading. American Journal of Orthodontics and Dentofacial Orthopedics. 123 (6), 657-665 (2003).
  29. Osipenko, M. A., Nyashin, M. Y., Nyashin, Y. I. Centre of resistance and centre of rotation of a tooth: the definitions, conditions of existence, properties. Russian Journal of Biomechanics. 3 (1), 5-15 (1999).
  30. Dathe, H., Nägerl, H., Dietmar, K. M. A caveat concerning center of resistance. Journal of Dental Biomechanics. 4, 1758736013499770(2013).
  31. Hohmann, A., et al. Influence of different modeling strategies for the periodontal ligament on finite element simulation results. American Journal of Orthodontics and Dentofacial Orthopedics. 139 (6), 775-783 (2011).

访问受限。请登录或开始试用以查看此内容。

重印与许可

标签

Abaqus