不同牙槽骨有限元模型对牙周膜生物力学响应的影响
发布时间：2021-04-25
【摘要】口腔正畸研究中，牙颌组织模型对有限元分析结果有着重要影响。目前，常用的牙槽骨有限元模型主要有均质和非均质两种类型。均质型模型对牙槽骨进行整体材料赋值，每个网格单元的力学特性均相同;而非均质型模型则依据牙槽骨每个网格单元所处的计算机断层扫描(CT)图像灰度值对网格单元进行差别赋值。为了探究不同牙槽骨模型对牙周膜生物力学响应的影响，我们选择一例临床患者作为研究对象，构建了其下颌尖牙、牙周膜和两种牙槽骨模型，分别对尖牙施加1N的根向压入力和2Nmm根向扭矩，分析牙周膜在不同载荷形式下的静水压应力和最大对数主应变。研究结果表明无论尖牙平动或转动，牙周膜的力学响应均会受到牙槽骨模型的影响。相比于均质型模型，当牙槽骨模型为非均质模型时，尖牙沿根向平动时牙周膜的最大应力和应变分别降低了13.13%和35.57%;尖牙沿根向转动时牙周膜最大应力和应变分别降低了19.55%和35.64%。均质型牙槽骨模型会诱导正畸医生选择一个较小的矫治力，非均质牙槽骨模型更能体现真实牙槽骨骨质特性的差异性，更有利于获得准确的分析结果。
【关键词】口腔正畸;有限元模型;牙槽骨;牙周膜;生物力学
引言
临床正畸治疗中，错位牙齿的矫正主要依靠患者口腔佩戴的正畸矫治器。矫治器形变产生矫治力，矫治力首先作用在错位牙齿的牙冠上，并通过牙体传递给牙周膜，进一步经过牙周膜的缓冲、吸收作用于牙槽骨，进而激发牙槽骨改建，最终在矫治力的作用下，牙齿逐渐移动到理想位置[1-4]。在牙齿移动过程中，牙周膜作为介于牙齿和牙槽骨之间的一层结缔纤维组织，在激发骨改建、诱导牙齿移动中起着重要作用[5-6]。有研究表明牙齿矫治过程中的移动速度与牙周膜的应变密切相关[7-8]，牙根的吸收情况则与牙周膜的应力有关[9-10]。因此，研究牙周膜受载之后的生物力学响应对于明确正畸矫治机制、优化正畸治疗方案具有非常重要的意义。
对于牙周膜生物特性的研究方法主要包括实验测试和有限元模拟。牙周膜作为人体生物组织，一方面实验材料来源非常有限;另一方面，牙周膜一般厚度仅有0.2～0.3mm[11]，极易遭到破坏而失去生物活性，存放条件非常苛刻[12]。另外，对于牙周膜的临床测试也因患者意愿、装置结构复杂、跟踪记录周期长、数据处理繁琐等因素难以获得直接有效的试验结果[13]。因此，在实验条件有限的情况下，更多的人采用有限元法对其进行研究。目前有限元法已广泛用于口腔正畸医疗领域，具有准确度高、重复性好和可操作性强等优点，受到越来越多的研究人员青睐[9,14]。
在采用有限元法对正畸矫治机制进行研究时，构建牙颌组织有限元模型至关重要。正畸过程中涉及到的牙颌组织主要包括牙齿、牙周膜和牙槽骨。牙齿相比于其他牙颌组织具有很高的硬度，弹性模量可达20GPa[15-16]，因此在有限元分析中多把牙齿设定为刚体[10,17]。牙周膜因其厚度较薄、尺寸较小，则一般假定为厚度均匀、均质、各向同性的超-黏弹性体[12,18]。对于牙槽骨模型的构建方法则主要有两种：第一种是采用非均匀有理B样条(non-uniformrationalb-splines，NURBS)曲面拟合的方式将牙槽骨三角面片网格模型转化为实体模型，然后将其导入有限元软件进行整体网格划分、材料赋值得到牙槽骨有限元模型[3,19]。这种建模方法获得牙槽骨模型每个网格单元的材料特性都相同，属于均质型模型。虽然有些学者按照牙槽骨骨质特性将其人为地分为皮质骨和松质骨[7,19]，但对于整个皮质骨或者松质骨而言，模型的每个单元材料特性也都一样。另一种方法则直接利用牙槽骨三角面片网格模型生成体网格模型，并依据模型中每个网格单元所处的计算机断层扫描(computedtomography，CT)图像灰度值对单元进行差别材料赋值[18,20]，不再利用有限元软件进行网格划分和材料赋值，这种建模方法获得牙槽骨模型每个网格单元的材料特性均不相同，属于非均质模型。
目前，上述两种建模方法均有使用，如Liao等[9]为了研究牙周膜在正畸力优化中的作用构建了精细的均质型牙槽骨模型，Ashrafi等[15]、Cozzani等[16]和Benaissa等[19]也都在自己最新的研究中采用了第一种建模方法;而Wu等[18]、Liu等[20]和Bujtar等[21]则采用第二种方法构建了非均质牙槽骨模型用于有限元分析计算。Liu等[20]进一步对比分析了两种建模方法在构建下颌骨有限元模型中的优劣势，结果表明非均质模型更有利于构建具有复杂几何构造的颌骨模型。而针对不同类型的牙槽骨模型，牙周膜受载之后的生物力学响应是否存在差异却鲜有报道。多数学者习惯采用第一种建模方法构建均质型的牙槽骨模型，对于模型的材料赋值也通常机械式地借鉴他人的材料参数，然而有限元模型的精确度对其分析结果却有着不可忽视的影响。
为此，我们将基于一例临床患者的螺旋CT图像数据，选择其下颌尖牙作为研究对象，构建其牙颌组织有限元模型，分别将其牙槽骨设置为均质型模型和非均质型模型，对尖牙施加牙根向的压入力和扭矩，研究不同牙槽骨模型作用下的牙周膜应力和应变，进而明确不同牙槽骨仿真模型对牙周膜生物力学响应的影响，协助相关研究人员选择合适的牙槽骨建模方法，从而提高有限元分析的准确性。
1材料和方法
1.1牙颌组织有限元模型构建
参照医学伦理审查办法，在受试者知情同意的前提下选择浙江大学医学院第二附属医院一例临床病例作为研究对象，依据患者的螺旋CT图像数据构建其下颌尖牙牙颌组织有限元模型。首先，将患者CT数据导入医学影像处理软件Mimics15.0(Materialise公司，勒芬，比利时)，通过阈值设定、掩膜编辑、区域增长和三维建模等手段建立尖牙及其牙槽骨的三角面片网格模型。接着将尖牙和牙槽骨的三角面片网格模型导入网格模型处理软件3-Matic7.0(Materialise公司，勒芬，比利时)，将尖牙模型向外偏置0.25mm之后与牙槽骨模型求交，然后减去原来的尖牙模型，获得具有均匀厚度0.25mm的牙周膜三角面片网格模型[18,22]。进一步利用3-Matic7.0软件的网格编辑功能，对初步生成的三角面片网格模型进行网格细化、均匀等优化处理，然后基于模型表面的三角网格直接生成体网格模型，该体网格模型可直接被有限元分析软件Abaqus6.13(达索SIMULIA公司，波士顿，美国)调用，且不需要再进行网格划分。导入有限元Abaqus6.13之后牙颌组织模型如图1所示，主要包括尖牙、牙周膜和牙槽骨，其中牙周膜和牙槽骨的网格单元数量分别为3416和1915。为了保证模型精度和网格质量，我们进行了相应的模型收敛性分析。当牙周膜和牙槽骨模型进一步精细化，网格数量均增加1倍时，牙周膜的最大静水压应力变化量小于5%(模型收敛性测试时的标准误差)。因此，在这项研究中模型的精度和网格质量是可接受的。模型中网格单元类型均为四节点线性杂交四面体单元(C3D4H)。
2结果
当对尖牙施加1N大小的根向压入力时，牙周膜对①均质型(E=490MPa)和②非均质型牙槽骨模型的应力和应变响应如图3所示。牙周膜应力为正值表示牙周膜处于受压的状态，反之，牙周膜应力为负值表示牙周膜处于受拉的状态。牙周膜应变为正值表示牙周膜处于受拉的状态，而牙周膜应变为负值则表示牙周膜处于受压的状态。从图3中可以看出，尖牙受1N根向压入力时，对于不同类型的牙槽骨模型，牙周膜最大和最小应力应变出现的位置基本相同，但是具体数值却不一样。如表2所列，与①均质型牙槽骨模型(E=490MPa)相比，当牙槽骨模型为②非均质模型时，牙周膜的最大应力降低了13.13%，同时最小应力提高了49.29%;牙周膜的最大应变降低了35.57%，最小应变则提高了18.83%。可见，对于不同类型的牙槽骨模型，在牙齿平动时牙周膜的应力和应变响应均存在一定差异。
当对尖牙施加2Nmm的根向扭矩时，牙周膜对①均质型(E=490MPa)和②非均质型牙槽骨模型的应力和应变响应如图4所示。从图4中可以看出对于不同类型牙槽骨模型①和②，牙周膜最小应力、最大和最小应变出现的位置基本相同，但牙周膜最大应力出现的位置却不相同。此外，牙周膜最大和最小应力应变值也因牙槽骨模型的不同而存在一定差异。如表2所列，与①均质型牙槽骨模型相比，当牙槽骨模型为②非均质模型时，牙周膜的最大和最小应力分别降低了19.55%和40.05%，牙周膜最大和最小应变分别降低了35.64%和61.13%。可见，对于不同类型的牙槽骨模型，在牙齿转动时牙周膜的应力和应变响应也存在一定差异。
如图4所示，当对尖牙施加2Nmm的根向扭矩时，不同弹性模量的均质型牙槽骨模型①和③对牙周膜应力应变响应的影响也存在一定差异。当牙槽骨弹性模量提高至4900MPa之后，所引起的牙周膜最小应力、最大和最小应变出现的位置基本相同，但具体数值因牙槽骨弹性模量的不同而有所差异。如表2所列，当牙槽骨的弹性模量提高10倍之后，牙周膜的最大和最小应力值分别降低了19.92%和38.73%，牙周膜最大和最小应变值也分别降低了35.21%和59.59%。可见，矫治力作用下牙周膜的生物力学响应很大程度上受牙槽骨弹性模量的影响。
3讨论
通过上述的有限元分析可知，无论在牙齿平动或是转动时，牙槽骨模型均会对牙周膜的生物力学响应产生影响。当牙槽骨由①均质模型转为②非均质模型时，牙周膜的最大应力值在一定程度上降低了。有研究表明牙周膜的应力与牙齿移动过程中的牙根吸收有关[9-10]，当牙周膜应力过大时会造成牙周膜局部血流闭塞，诱发牙根吸收[1-2]。可见采用均质型牙槽骨模型研究矫治力作用下牙周膜的应力时更偏向于获得较高的应力值，进而选择一个更保守的矫治力。另一方面，牙周膜的应变与矫治过程中牙齿的移动速度相关，在一定范围内牙周膜的应变越大则牙齿的移动速度就越快[7-8]，当采用均质型牙槽骨模型研究牙周膜的应变分布时会获得一个较大的应变值，进而认为牙齿已具有较高的移动速度，同样会引导正畸医生选择一个较小的矫治力。矫治力过小会延长正畸治疗周期，降低患者配合治疗的意愿，难以实现理想的正畸治疗[4]。
另外，我们假设将①均质型牙槽骨模型弹性模量增大10倍，重点分析了在尖牙受到2Nmm的根向扭矩时，牙周膜对不同弹性模量的牙槽骨力学响应的差异性。从有限元分析结果(表2)可以看出在模型网格和节点完全相同的情况下，当牙槽骨弹性模量提高之后，牙周膜的最大应力和最大应变均有所降低。所以，采用高弹性模量的牙槽骨模型进行有限元分析时，分析结果会引导正畸医生选择较高的矫治力，相反，低弹性模量的牙槽骨模型则会诱导医生选择一个更低的矫治力。然而，更高或者更低的矫治力均不利于正畸治疗[4,8-9]。为了获得更加安全高效的正畸治疗，明确载荷作用下牙周膜的生物力学响应至关重要。因此，根据患者牙槽骨的实际情况，设定合适的牙槽骨力学参数，对于采用有限元法研究正畸治疗非常重要。
牙槽骨作为结构复杂的人体生物组织，内部含有神经管、骨小梁、骨髓小腔等组织结构，不同区域的力学特性差异很大[7,17]。位于牙槽骨表层的骨组织骨质致密、强度高，称为皮质骨;位于牙槽骨内部的骨组织则骨质疏松，呈海绵状，称为松质骨[20,27]。在认识到牙周膜生物力学响应受牙槽骨力学特性影响的前提下，构建符合人体牙槽骨生物力学特性的有限元模型对于提高有限元分析结果的精确度有着积极促进作用。相比于均质型牙槽骨模型，非均质型牙槽骨模型模糊了皮质骨和松质骨的界限，依据每个网格单元所处的CT图像灰度值对网格单元进行材料赋值，突显出牙槽骨不同区域骨质特性的差异性，更加贴近牙槽骨实际状况，有利于获得更加准确的分析结果。值得注意的是，非均质牙槽骨模型的整体材料属性在一定程度上受其网格单元尺寸的影响。模型中网格尺寸越小，每个单元越能依据其所在的微小局部CT图像进行个性化材料赋值，越容易获得更加细致的、差异性的牙槽骨模型。当网格尺寸变大时，其所覆盖的CT图像区域也更大，即便覆盖区域的CT图像有一定差异，也无法对该网格单元实现多种材料赋值。如图2a所示，从牙槽骨的横断面来看，皮质骨和松质骨的分界面并不明显，其中一部分原因就是该断面上牙槽骨网格单元的尺寸过大。不过，从模型收敛性测试结果来看，边缘区域过大的牙槽骨网格单元并没有对牙周膜的力学响应产生不可忽视的影响。
在非均质牙槽骨模型中，采用CT图像灰度值进行材料赋值的网格单元中最大的弹性模量是14420.28MPa，其量级与其他研究中报道的皮质骨弹性模量13700MPa非常接近[7,17]，间接证明了材料赋值公式的准确性。另外，如图2a所示，牙槽骨模型表面的部分区域被表示低弹性模量的蓝色所覆盖，这与处于真实牙槽骨表面的皮质骨力学特性不符。这主要是因为人体牙槽骨外还附着有骨膜和牙龈等组织结构，在牙槽骨模型构建过程中，有一部分口腔软组织被当作牙槽骨进行了重建，但这并不会对有限元分析结果产生影响，从中我们也可以发现非均质模型的另外一个优点。一般来说，对于牙槽骨模型边界的划定主要依靠人工进行手动分割，操作者的熟练程度和工作经验直接影响着牙槽骨模型的结构和大小，如果采用均质模型进行建模，那么模型构造上可能就与真实的牙槽骨存在结构差异，在一定程度上影响了分析结果。而当采用非均质建模方法时，因模型中的所有网格单元均依据其所在的CT图像灰度值进行材料赋值，即便有一部分外层软组织或者牙槽骨内部的骨髓小腔、神经管空腔等力学特性远异于骨组织的组织结构，被误以为牙槽骨进行有限元建模，也会因其CT图像灰度值偏低被赋予相应的材料参数，不会因模型精度影响有限元分析结果。
4结论
以一例临床病例作为研究对象，基于患者的CT图像构建了其尖牙、牙周膜和牙槽骨的有限元模型。其中，牙槽骨模型分别设定为均质模型和非均质模型，其他牙颌组织模型设定均相同。通过对尖牙施加根向矫治力和扭矩探究牙周膜在不同牙槽骨模型下的应力应变响应。结果表明：(1)无论牙齿平动或是转动，当牙槽骨为均质模型时，牙周膜在同等条件下会产生更大的应力和应变，会误导正畸医生选择一个较小的矫治力;(2)牙周膜的生物力学响应受牙槽骨力学参数影响较为显著，当牙槽骨弹性模量由490MPa增大10倍至4900MPa时，牙周膜的最大应力和应变值分别降低了19.92%和35.21%;(3)相比于均质牙槽骨模型，非均质牙槽骨模型依据每个网格单元所处的CT图像灰度值对模型进行材料赋值，更能体现真实牙槽骨组织骨质特性的差异性，有助于提高有限元分析的精确性。研究结果为构建更加精确的牙颌组织有限元模型提供了技术支持。——论文作者：仵健磊1，刘云峰2,3，李伯休4，王栋才2,3，董星涛2,3，周佳丽5
相关期刊推荐：《中国生物医学工程学报》是中国生物医学工程学会主办的专业性学术期刊是全国生物医学工程工作人才沟通学术思想、交流学术经验的园地。反映国内外该领域具有创造性的研究成果，促进国内外学术交流。报道有创造性的最新科研论文和研究简讯。包括学科有：人工器官和生物医用材料，生物效应，生物信息与控制;生物力学;生物医学仪器;中医工程;临床工程等。
