Summary
本文提出的方法是建立和验证体外3D 模型, 能够测量由两个支架之间放置的 V 形弯曲的不同矫正线所产生的力系统。另外的目标是比较这个力系统与不同类型的矫正线和以前的模型。
Abstract
正确理解各种正畸器械所产生的力系统, 可以使病人的治疗效率和可预测性。将复杂的多支架装置减少到简单的双支架系统, 以进行力系统评估, 将是朝这个方向迈出的第一步。然而, 这方面的正畸生物力学大部分限于2D 实验研究, 计算机建模/分析或理论外推现有模型。本协议的目的是设计、构造和验证一种体外3D 模型, 能够测量在两个支架之间放置 V 形弯曲的弓产生的力和力矩。另外的目标是比较不同类型的矫正线在他们自己和以前的模型中产生的力系统。为此, 我们模拟了一个 2 x 4 装置, 代表一个臼齿和一个门牙。正畸导线测试仪 (OWT) 是由两个多轴力传感器或负载细胞 (纳米传感器) 构成的, 正畸支架连接。负载单元能够测量所有三平面空间中的力系统。测试了两种类型的矫正线, 不锈钢和β钛三不同大小 (0.016 x 0.022 英寸, 0.017 x 0.025 英寸和 0.019 x 0.025 英寸)。每根导线接受一个垂直的 V 弯被系统地安置在一个特定位置以预定义的角度。类似的 V 形弯曲在不同的矫正线上复制, 在11个不同的位置之间的磨牙和切牙附件。这是第一次在体外尝试使用 V 形弯曲对不同矫正线的正畸装置进行模拟。
Introduction
临床正畸治疗的一个重要方面是了解 multibracket 装置产生的力系统。明确了解基本的生物力学原理有助于提供可预测的结果, 并尽量减少潜在的副作用1。近年来, 随着支架位置和设计的增加, 在矫正线中放置弯曲的趋势越来越远;然而, 综合正畸治疗仍需要在矫正线中放置弯曲。弯曲, 当放置在不同的类型和大小的矫正线, 可以创建各种力量系统, 适合不同类型的牙齿运动。虽然在考虑多齿时, 力系统会变得相当复杂, 但一个有用的起点可能涉及简单的双支架系统。
迄今为止, V 型弯曲力学主要是在第二阶分析, 利用数学模型1,2,3,4,5和/或基于计算机的分析/模拟6. 这已经对涉及相邻支架的拱导线的二阶相互作用所涉及到的力系统有了基本的理解 (图 1)。然而, 这些方法施加一定的边界条件, 以运行模拟, 可能不会在实际的临床情况下, 可能会发生偏差。最近, 提出了一种新的体外模型, 涉及力传感器, 以测量三维 (3D) 的力量和创造的时刻, 不仅评估第二级弓支架相互作用, 但也在第三级7。然而, 不同类型的矫正线对沿切牙磨牙弓的不同弯曲位置的力系统的影响没有得到评价。此外, 这项研究仅涉及弹性正畸矫正线的评价, 这不是牙齿运动发生的主要矫正线。因此, 本研究的目的是评估在不同位置的方形不锈钢和β钛矫正线在一个3D 设置涉及磨牙和切牙支架所产生的力系统。临床医生需要知道的力量系统应用于牙列时, 具体组合的弓支架组合被用来修复错用。
本文研究了三平面空间的正畸力系统, 模仿临床实际。据了解, 临床上测量力系统是极其困难的;因此, 这种测量必须在体外进行。假设在实验室中由 V 形弯曲产生的力系统在病人口中被复制是相似的。创建了一个工作流来评估如何配置实验设置 (图 2)。
正畸钢丝测试仪 (OWT) 是一种创新产品, 由正畸的分工与生物工程 & 北京大学生物动态实验室, 康涅狄格健康, 法明顿, CT, 美国 (图 3) 合作开发。它的目的是准确模仿上颌牙齿的安排, 在口腔内和一些内部的条件, 同时提供测量的力量系统创建的所有三平面空间。OWT 的主要机械部件是数据采集装置 (数据采集器)、纳米力/扭矩传感器、湿度传感器、温度传感器和个人电脑。测试装置放置在具有温度/湿度控制的玻璃外壳内。这允许部分模拟口腔内环境。该数据采集器作为三传感器的接口: 湿度传感器、力/力矩传感器、热敏电阻以及位于平台上的传感器的测试设备 (图 3)。这些链接与一个软件程序。软件是可视化编程的平台和开发环境, 用于控制不同类型的硬件。它被选择自动化的正畸线测试仪。
在检测仪器上设置了一系列铝钉, 以表示上颌牙弓的牙齿。两个代表右中央切牙和右第一磨牙的钉子连接到传感器/负载单元 (S1 和 S2)。负载单元是一种机械装置, 可以测量所有三个平面 (x-y) 中应用到它的力和力矩: fx、fy和 fz;mx, my, 和 mz。该钉是有系统地定位, 以创建牙科弓形式。每一个 peg 是由一个精确记录的测量分离, 用平均牙齿宽度计算, 在患者正畸治疗。为实验选择的形状是由标准模板创建的 "卵球形" 拱形式。
Subscription Required. Please recommend JoVE to your librarian.
Protocol
1. 实验设置
- 使用定制的 "夹具", 标记在 OWT 铝钉上放置磨牙管和门牙托架的精确位置。
- 粘结标准自结扎支架与复合材料。光治疗40秒。
- 将 0.021 x 0.025 英寸不锈钢 (SS) ' 卵球形 ' 上颌弓插入支架插槽中。
- 将测试仪器放在玻璃腔内。
- 检查是否有任何意外的弓激活。任何激活的弓将自动创建一个力系统, 将显示在计算机屏幕上。
- 如果观察到任何弓激活, 则重新定位托架。重复步骤 1.2-1.5。
2. 制作模板弓 (图 4)
- 在测试仪器中放置一个弓 (0.021 x 0.025 SS)。
- 使用永久性标记来指示以下内容: 1) 中线, 2) 一点, 紧靠门牙支架 (I) 和 3), 立即内侧到磨牙管 (M)。对弓的对侧也做同样的事。这是模板拱线。
- 将弓与标记的点转移到图表纸上。
- 在图形纸上制作弓的精确副本。
注: 此图纸可用于确定样品的所有矫正线的 V 形弯的位置。 - 计算拱线段 (L) 从 I 到 m 的周长。
- 现在, 马克11点从我到 m。每一点都是未来的 V 折弯位置。
- 将每个点标记为0到10。
- 请确保每个折弯位置与另一项之间的大小相等。
- 通过计算每个位置的/L 为每个折弯位置获取唯一的数字/比率。
3. V 型弯的放置
- 从样品中取一个新的弓。
- 把它放在模板弓/图形纸上, 并将其中一个十一折弯位置双边地转移到弓。
- 使用矩形弓钳或轻钢丝钳, 使对称的 V 弯在两个位置。
- 将弓放在玻璃板/平面平台上, 并检查弓两端与量角器的角度测量。
- 如果需要, 调整两端, 以便创建150°的角度。
- 对示例的所有矫正线重复步骤3.1 到3.5。
4. 测量力系统 (图5和 6)
- 打开用于数据记录的软件程序 (参见材料表)。
- 为要保存的数据创建一个新文件夹。
- 单击 "运行" 启动软件。该程序将在每个传感器上实时显示三个力和三力矩值。
- 等待大约10-15 秒的数据记录软件的波动停止。确保该软件上的图形线为部队系统的所有组件显示 "平坦" 线。
注: 所有六测量在每个传感器将显示微不足道的价值 (力量 < 1 g 和时刻 < 10 g 毫米)。 - 从平台上轻轻地移除 "测试设备"。使用 Weingart 钳将弓插入到摩尔管中。
- 用牙周洁打开门牙支架的门。
- 提起弓的前部, 并将其插入托架插槽中。确保弓中线与测试仪器的中线重合。
- 将测试仪器返回平台并关闭玻璃室的门。
- 将温度设置为摄氏37摄氏度。等待一分钟的温度, 玻璃室调整。
- 单击软件上的 "开始保存" 按钮, 并允许软件保存/传输数据至少10秒。再次单击 "开始保存" 按钮以结束数据传输, 然后单击 "停止"。
注: 每个测量周期在10秒内为每个组件 (fx、fy、fz、mx、my和 mz) 生成100读数。 - 转到包含已保存数据的文档, 并将数据集复制/导出到自定义设计的数据分析电子表格 (请参阅辅助表)。选择正确的 V 形折弯位置编号和特定导线样本以插入数据。
- 对该特定折弯位置的10矫正线重复步骤4.3 到4.11。
- 现在, 将矫正线的计算方法和标准偏差复制到单独的电子表格中, 以创建数据的图形表示形式。
- 对所有折弯位置和矫正线类型重复步骤4.2 到4.13。
注: 矫正线包括不锈钢 (SS) 和β钛 (ß钛), 具有以下尺寸: 0.016 x 0.022 英寸, 0.017 x 0.025 英寸, 0.019 x 0.025 英寸。
5. 错误评估
- 按照步骤 4.1-4.4 中的说明运行计算机/软件
- 从平台上取出 "测试设备"。
- 获得一个直长 0.021 x 0.025 英寸 SS 线。使用轻线钳, 弯曲线的一端成一个小钩。将弓的游离端插入从远端的摩尔管中。
- 将测试仪器放回平台上。
- 将已知重量 (50 克) 附加到挂钩上。让它在垂直平面上自由悬挂, 消除任何类型的干扰。关上玻璃室的门。
- 执行步骤 4.10-4. 11。
- 对切牙支架重复步骤 5.1-5.6。
- 输入作为 "测量值" 的摩尔管的托架和 Mx 的熔体值。
- 现在应用平衡等式 (参见补充文本) 计算 "预期值"。
Subscription Required. Please recommend JoVE to your librarian.
Representative Results
传感器板中心的每个传感器所经历的总力和总力矩由其三个正交分量表示: fx、fy和 z 轴, 表示沿 x 轴、y 轴和 z 轴线的作用力,分别mx, my, mz表示在同一轴周围的时刻。传感器的初始测量将数学转换为支架所经历的力和力矩值 (图 7)。
一系列图表显示垂直力在摩尔 (熔米) 和切牙托架 (熔i), 力矩 (牙冠小费) 在摩尔支架 (mxm), 力矩/扭矩 (labio 舌小费) 在门牙支架 (mxi) 相对于从原始数据中创建的与单个齿坐标系相关的一/升比率。牙冠比表示每个 V 形弯曲的位置, 其中 ' a ' 是切牙支架的远端边缘和 V 形弯的顶端之间的距离, 和 ' L ' 的摩尔管的内侧边缘和切牙托架的远端边缘之间的距离easured 沿弓 (37 毫米)。a/升比率为 0.0 (0 毫米/37 毫米) 表示与门牙支架相邻的折弯, 每一个连续的弯曲 (a/l = 0.1, 0.2等) 的间距3.7 毫米远离上一个折弯结束与/升 = 1.0 (37 mm/37 毫米), 代表一个弯曲相邻的摩尔括号。力系统的方向由一个负/正符号表示。图表按导线类型和大小分组 (图 9和 10)。图上的每个点都表示十相似矫正线的平均值, 误差线表示上述平均值的一个标准偏差。靠近水平轴的点 (上面或下方) 表示具有低震级的力或力矩, 而从水平轴 (上方或下方) 远一点的点表示具有较高震级的力或力矩。
垂直力 (FZ) 显示六线类型的对称和线性模式 (图 8)。接近 V 弯到任一支架, 更高是垂直力量。当折弯从支架上移开, 向中心移动时, FZ的大小减小, 直到到达某一点时, 两个力都近似为零 (中性区)。当弯向更远的位置移动, FZ逐渐增加。然而, 个别力 (f 和 f子)的方向是颠倒的。定量地, SS 矫正线创造了一个明显地更大的力量系统比ß-Ti 矫正线。此外, 更高维度矫正线创建更大的力系统。令人惊讶的是, 矫正线在两个支架上创建的相对力系统在弓的大小和类型方面都相当相似。
相反, 矩 (MX) 显示了非线性和非对称模式 (图 9)。当 V 弯被放置在接近的摩尔管 (x/升比率 > 0.6), 以及在摩尔管 (红色) 的力矩方向逆转从x/升的0.0 到 0.2, 是类似于所有矫正线, 也许代表弓支架相互作用和支架方向 (二阶与第三阶) 的一个更基本的性质。在两个托架上的矩比率显示了在所有矫正线测试中观察到的特定模式 (图 10)。弯曲, 被放置在接近门牙 (0.0-0.3 为ß和 0.0-0.2 为 SS) 有两个时刻在同一方向 (mxi/mxm > 0)。从 0.3-0.6 的ß和 0.3-0.4 的 a/升为 SS, 瞬间是相反的方向 (mxi/mxm < 0) (中性区)。弯曲在0.6 或更大没有创造一个重要时刻在门牙 (≈0 g 毫米), 但一个巨大的片刻在摩尔管 (mxi/mxm≈0) 产生。
定量地, 再和垂直力一样, SS 弓所产生的力矩的大小在统计学上和临床上都大于ß-Ti 矫正线产生的, 无论是升比还是拱线的大小。
百分比错误由以下等式计算:
小于50克的权重的% 误差被发现为 5%, 从50到500克的权重计算为0.5%。
中性区 (相等和相反弯矩) 被发现在 0.3-0.4 的ß比率为 SS 矫正线和 0.4-0.5。在这些特定的弯曲位置, 垂直力是最小的时刻作用在切牙和摩尔括号对面的方向。根据 a/d 比率, 在磨牙和切牙支架之间由 V 形弯曲产生的力系统可分为三种不同的类别 (图 11)。
图 1: 在第二个顺序中由两个共线托架创建的强制系统.L 是两个括号之间的距离;a 是 V 形弯从支架 a 的位置;fa和 fb是分别在支架 A 和 b 上产生的垂直力;Ma是片刻在 a;Mb是括号 b 的时刻.请点击这里查看这个数字的大版本。
图 2: 工作流.请单击此处查看此图的较大版本.
图 3: 正畸钢丝测试仪 (OWT).A: 检测仪器, B: 测量平台, C: 温度监测器。请单击此处查看此图的较大版本.
图 4: 两个附件之间的折弯位置的示意图表示形式.每一个蓝色点是一个折弯位置, 代表距离 ' a ' 测量从切牙支架沿弓。将有11个不同的值为 ' a ' 的增量为3.7 毫米. (即蓝点与相邻的蓝点分开3.7 毫米)。L 是从切牙支架远端表面到沿弓的摩尔管远端表面测量的周长长度。请单击此处查看此图的较大版本.
图 5: 弓插入和持有的支架上的铝钉连接到传感器.请单击此处查看此图的较大版本.
图 6: 显示原始数据 (蓝色和红色盒) 的软件程序, 从连接到门牙和臼齿托架的两个传感器 (S1 和 S2) 获得.请单击此处查看此图的较大版本.
图 7: X Y Z 坐标及其与 OWT 的方向关系.X: 横向平面;水平平面;Z: 垂直平面。请单击此处查看此图的较大版本.
图 8: 垂直力 (熔) 在两个托架上的图形表示.请单击此处查看此图的较大版本.
图 9: 在两个托架上的横向平面 (Mx) 中的矩的图形表示.请单击此处查看此图的较大版本.
图 10: 相对力系统横跨不同的弓类型和大小通过时刻的比率描述 [Mx (i)/mx (m)].请单击此处查看此图的较大版本.
图 11: 三不同的力量系统从 V 弯.每个区域代表一个唯一的 f# 系统。"蓝色" 阴影区域描述了与类似的相对力系统的 a/d 比率。请单击此处查看此图的较大版本.
| 力系统 | 摩尔 (m) | 切牙 (i) |
| 熔 (+) | 侵入 | 侵入 |
| 熔 (-) | 突出 | 突出 |
| Mx (+) * | 内侧提示 | 面部/唇尖 |
| Mx (-) * | 远端尖端 | 腭/舌尖 |
| * 所有测量均在支架上进行 | ||
表 1: 签署了部队系统的公约和方向。
补充图 1: x 轴周围力矩的平衡图 (Mx).注意: 图只是比较矩的大小。mx(m) + Mx(i)和 fz(m)或 fz(i) x D 的方向将始终是相反的。因此, ΣMx= 0 (参见补充文本)。请单击此处查看此图的较大版本.
补充案文. 请单击此处下载此文件.
Subscription Required. Please recommend JoVE to your librarian.
Discussion
正畸矫正线已研究了不同的方式8,9,10,11。他们也被评估了各种机械性能, 但他们很少被分析, 以确定的力量系统, 他们将创建12,13,14,15。三点弯曲试验普遍用于评价正畸矫正线;然而, 它们通常是在没有任何弯曲的直导线上进行的。体外评估通常是优化, 以观察只有1或2个变量的时候, 不允许的结果是容易适应的临床情况。本研究的重点是试验性地确定由垂直 V 形弯曲所产生的3D 力系统, 在 interbracket 距离上放置在不同位置的矩形矫正线作为 2 x 4 装置。此协议与以前分析 V 形弯曲力学的方法有很大区别。这是第一次建立一个实际的体外设置是利用纳米传感器模仿工作的两个支架-弓几何, 而不是依赖于计算机模型或有限元方法。这个机械模型不仅测量弯矩 (第二阶导线支架相互作用), 而且扭矩 (三阶导线支架相互作用)。不施加边界条件。换句话说, 以前的研究从来没有考虑到弓的曲率, 因为它从磨牙到门牙的支架。由于这一曲线, 切牙和摩尔括号没有定位在同一平面上, 也不是平行于彼此。这种安排可以增加对力系统的分析的复杂性, 这使得他们在临床上比仅涉及两个相同的括号排列在一条直线和平行3,4更相关。
传感器的功能和数据输出可以很容易地影响的因素, 如设备的误差, 传感器灵敏度, 过热的 OWT, 人为错误的电线活化, 弯曲, 结扎, 形状, 不当的导线定位, 失活的导线在最后插入之前, 变形弓等, 因此, 用新的矫正线进行重复测量, 并运用平衡定律验证数据是很重要的。此外, 只有少数矫正线应插入测量, 以避免过热的 OWT。
每个折弯位置仅由3.7 毫米分开。因此, 准确地放置 V 型弯沿弓也是重要的。小偏差从期望的位置可能从根本上改变被记录的力量系统。一种自定义的图形纸, 包含带有 V 形弯曲位置的弓模板, 有助于达到所需的精度。铝钉上的托架定位不当也可以做同样的事。因此, 如果出现粘结故障, 则使用定制的精密夹具来获得支架的位置。
在实验过程中, 如果支架得到脱胶, 一个新的支架必须精确地放回同一位置。定制设计的夹具可以帮助找到所需的位置。无任何弯曲的被动矫正线必须用于确保支架的位置正确。如果没有, 它将不得不 rebracketed。重要的是不要重用脱胶支架, 因为有增加的可能性的支架变形。
当前方法的一个缺点是, 仅使用了两个传感器。增加更多的传感器将允许研究更复杂的力系统, 如那些包括三或更多的支架设置在拱。另一个潜在的缺点是无法模拟口腔环境。温度、唾液、闭塞等因素可能会影响产生的力系统。然而, 在这一点上是不可能同时测量的力量系统和观察到的牙齿运动在临床水平。
使用有限元分析 (FEM) 的计算机建模和模拟是一个快速出现的区域, 用于解码各种正畸器具的生物力学16、17、18、 19. 然而, 验证这些方法的一个先决条件是精确地将复杂的弓支架相互作用结合在一起, 并将假设保持在最低限度。在第二阶和第三阶中, 弓支架的相互作用很大程度上是未知的, 这可能会限制这些程序的准确性。为了更好地进行计算机模拟, 首先要找出在各种临床情况下存在的力系统, 生成一个相当可观的生物力学数据库, 然后建立一个基于该数据集的计算机模型。换言之, 更好的建模和预测将需要本协议提供的实际实验。
Subscription Required. Please recommend JoVE to your librarian.
Disclosures
作者没有什么可透露的。
Acknowledgments
作者希望感谢所有使这项工作成为可能的同事, 特别是 Drs. 阿迪亚齐柏和拉文答腊。作者感谢北京大学生物动态 & 生物工程实验室在康涅狄格健康方面所提供的设施在这个项目的发展。
Materials
| Name | Company | Catalog Number | Comments |
| Force/Torque Sensors/Transducers | Nano17 F/T Sensors, ATI Industrial Automation, Apex, NC, USA | Part of the OWT | |
| CHS Series Humidity Sensor Units | TDK Corporation | Part of the OWT | |
| Temperature sensors | (Murata NTSDXH103FPB30 thermistor) Murata Manufacturing Co., Ltd | Part of the OWT | |
| LabVIEW 7.1. | Laboratory Virtual Instrumentation Engineering Workbench, Version 7.1 | Software Program | |
| Self-Ligating brackets | Empower Series, American Orthodontics. | Orthodontic Brackets | |
| Stainless steel archwires | Ultimate Wireforms, Inc. in Bristol, CT | Archwires | |
| Beta-Titanium Archwires | Ultimate Wireforms, Inc. in Bristol, CT | Archwires | |
| Data acquisition device (DAQ) | National Instruments (NI) USB 6210 | Part of the OWT | |
| Ortho Form III (Archform template) | 3M Oral Care, St. Paul, MN, USA | Ovoid arch form | |
| Weingart Plier | Hu-Friedy Mfg. Co., LLC Chicago, IL | Orthodontic Plier | |
| Light wire Plier | Hu-Friedy Mfg. Co., LLC Chicago, IL | Orthodontic Plier |
References
- Burstone, C. J., Koenig, H. A. Force systems from an ideal arch. Am J Orthod. 65 (3), 270-289 (1974).
- Koenig, H. A., Burstone, C. J. Force systems from an ideal arch: Large deflection considerations. Angle Orthod. 59 (1), 11-16 (1989).
- Burstone, C. J., Koenig, H. A. Creative wire bending: The force system from step and V bends. Am J Orthod and Dentofac Orthop. 93 (1), 59-67 (1988).
- Ronay, F., Kleinert, W., Melsen, B., Burstone, C. J. Force system developed by V bends in an elastic orthodontic wire. Am J Orthod and Dentofac Orthop. 96 (4), 295-301 (1989).
- Demange, C. Equilibrium situations in bend force systems. Am J Orthod and Dentofac Orthop. 98 (4), 333-339 (1990).
- Isaacson, R. J., Lindauer, S. J., Conley, P. Responses of 3-dimensional arch wires to vertical V bends: Comparisons with existing 2-dimensional data in the lateral view. Semin Orthod. 1 (1), 57-63 (1995).
- Upadhyay, M., Shah, R., Peterson, D., Takafumi, A., Yadav, S., Agarwal, S. Force system generated by elastic archwires with vertical V bends: A three-dimensional analysis. Eur J Orthod. 39 (2), 202-208 (2017).
- Gurgel, J. A., Kerr, S., Powers, J. M., LeCrone, V. Force-deflection properties of superelastic nickel-titanium archwires. Am J Orthod Dentofacial Orthop. 120 (4), 378-382 (2001).
- Gurgel, J. A., Kerr, S., Powers, J. M., Pinzan, A. Torsional properties of commercial nickel-titanium wires during activation and deactivation. Am J Orthod Dentofacial Orthop. 120 (1), 76-79 (2001).
- Hazel, R. J., Rohan, G. J., West, V. C. Force relaxation in orthodontic arch wires. Am J Orthod. 86 (5), 396-402 (1984).
- Lundgren, D., Owman-Moll, P., Kurol, J., Martensson, B. Accuracy of orthodontic force and tooth movement measurements. Br J Orthod. 23 (3), 241-248 (1996).
- Goldberg, A. J., Burstone, C. J. An evaluation of beta titanium alloys for use in orthodontic appliances. J Dent Res. 58 (2), 593-600 (1979).
- Kusy, R. P., Whitley, J. Q. Thermal and mechanical characteristics of stainless steel, titanium-molybdenum, and nickel titanium archwires. Am J Orthod Dentofacial Orthop. 131 (2), 229-237 (2007).
- Kapila, S., Sachdeva, R. Mechanical properties and clinical applications of orthodontic wires. Am J Orthod Dentofacial Orthop. 96 (2), 100-109 (1989).
- Verstrynge, A., Humbeeck, J. V., Willems, G. In-vitro evaluation of the material characteristics of stainless steel and beta-titanium orthodontic wires. Am J Orthod Dentofacial Orthop. 130 (4), 460-470 (2006).
- Tominaga, J. Y., Tanaka, M., Koga, Y., Gonzales, C., Kobayashi, M., Yoshida, N. Optimal loading conditions for controlled movement of anterior teeth in sliding mechanics. Angle. 79 (6), 1102-1107 (2009).
- Cattaneo, P. M., Dalstra, M., Melsen, B. The finite element method: A tool to study orthodontic tooth movement. J Dent Res. 84 (5), 428-433 (2005).
- Fotos, P. G., Spyrakos, C. C., Bernard, D. O. Orthodontic forces generated by a simulated archwire appliance evaluated by the finite element method. Angle Orthod. 60 (4), 277-282 (1990).
- Geramy, A. Alveolar bone resorption and the center of resistance modification (3-D analysis by means of the finite element method. Am J Orthod Dentofacial Orthop. 117 (4), 399-405 (2000).
