第二军医大学学报  2016, Vol. 37 Issue (7): 821-826   PDF    
两种种植体支持式下颌半口固定义齿种植体应力有限元分析
陆永健1,2, 唐震1, 龚璐璐3, 李笑梅1, 汪大林1     
1. 第二军医大学长海医院口腔科, 上海 200433;
2. 上海交通大学医学院附属新华医院口腔科, 上海 200092;
3. 同济大学生命科学与技术学院, 上海 200092
摘要: 目的 设计两种种植体支持式下颌半口全牙列固定义齿,研究在不同下颌骨垂直高度降低情况下种植体应力分布情况。 方法 利用锥束计算机断层扫描成像(CBCT)扫描下颌骨及上部义齿,形成下颌骨基准模型,将基准模型中下颌骨垂直高度分别下降0、5、10和15 mm,组装种植体后,建立颏孔前6枚平行种植体支持和“All-on-four”支持的不同下颌骨垂直高度的种植体有限元模型,分别为6-种植体模型和4-种植体模型,共两组(8个)。在右下颌第一磨牙中央窝静态加载250 N垂直力,用Ansys 15.0软件进行模型各部分的应力分析。 结果 在相同的加载条件下,种植体最大应力值位于颈部;两组8个模型最大应力值分布显示,不同下颌骨垂直高度下4-种植体组的最大应力值在40.12~49.06 MPa之间,6-种植体组的最大应力值在80.62~109.64 MPa之间;6-种植体组的最大应力值是4-种植体组的2倍。两组模型中种植体最大应力在下颌骨垂直高度降低5 mm时均表现为最低,随着下颌骨垂直高度降低,应力逐渐增大,而垂直高度下降0 mm的模型最大应力要略大于下降5 mm模型。 结论 4-种植体和6-种植体支持的无牙颌固定义齿在相对极端的载荷下种植体都没有出现破坏性应力,4-种植体有着更加合理的应力分布;下颌骨垂直高度下降5 mm组应力最小,说明适当降低垂直高度有利于种植体应力更加合理的分布。
关键词: 牙种植     三维有限元分析     无牙颌模型    
Three-dimensional finite element analysis of implant stress of fixed full-arch implant-supported prosthesis on edentulous mandible in two different designs
LU Yong-jian1,2, TANG Zhen1, GONG Lu-lu3, LI Xiao-mei1,     
1. Department of Stomatology, Changhai Hospital, Second Military Medical University, Shanghai 200433, China;
2. Department of Stomatology, Xinhua Hospital, Shanghai Jiaotong University School of Medicine, Shanghai 200092, China;
3. School of Life Sciences and Technology, Tongji University, Shanghai 200092, China
Supported by Scientific Research Project of Science and Technology Committee of Shanghai (12ZR1438400).
Abstract: Objective To design two kinds of fixed full-arch implant-supported prosthesis and to investigate their stress distribution in the different vertical heights of the edentulous mandible. Methods The cone-beam computed tomography(CBCT) scan data of the upper denture and mandible were collected to establish the base models of the mandible. The vertical heights of the mandible models were lowered by 0, 5, 10 and 15 mm, and the three-dimensional finite element models of implants with different vertical heights supported by mental foramen 6 parallel implants and supported by "All-on-Four" were established, including 2 groups (8 models):6-implant group and 4-implant group. A 250 N vertical force was applied to the central fossa of the right mandibular first molar. The stress of each part of the models was analyzed by Ansys 15.0 software. Results Under the same loading condition, the maximum stress value of the implants was in the neck of the implants. The maximum stress distribution of 8 models showed that the maximum stress value was 40.12-49.06 MPa for the 4-implant group and was 80.62-109.64 MPa for the 6-implant group, with the latter being two folds that of the former. In two groups, the maximum stress was the lowest when the vertical height reduction was 5 mm. With the decrease of the mandible vertical height, the stress increased gradually. However, the maximum stress of the 0 mm-height-reduced models was larger than that of 5 mm-height-reduced models. Conclusion There is no destructive stress under certain extreme load in both 4-implant group and 6-implant group, with the 4-implant group having a more reasonable stress distribution. The 5 mm-hight-reduced models have the minimum stress, indicating that appropriate vertical height reduction can lead to more reasonable implant stress distribution.
Key words: dental implantation     three-dimensional finite element analysis     edentulous mandible model    

下颌牙列缺失的患者在进行传统下半口义齿修复时常存在固位困难的问题,为解决这一问题,我们通常使用种植体支持式覆盖义齿或固定义齿修复。在相同条件下,无论是在咀嚼效率还是美观舒适方面,种植体支持式固定义齿都具有较大优势。

影响种植体支持式固定义齿设计及效果的因素很多,包括牙槽骨条件、咬合力大小、全身疾病、费用等因素,其中牙槽骨因素最为重要。全口牙列缺失的患者通常伴有不同程度的牙槽骨萎缩,而牙槽骨的萎缩程度对修复体设计具有至关重要的影响。牙槽骨萎缩分为垂直向萎缩和颊舌向萎缩。当牙槽骨颊舌向过度吸收时,如要采用种植体支持式固定义齿修复,则需进行骨移植手术,考虑到无牙患者的年龄及身体情况,一般改用种植体支持式覆盖义齿修复;而对于以垂直向吸收为主的无牙患者,可考虑采用种植体支持式固定义齿修复[1-2]。为了支持下颌固定义齿,一般采用颏孔前4~6颗种植体的设计方案[3-4]。从生物力学角度分析,在垂直距离不变的情况下,当下颌骨垂直高度逐渐降低时,种植体上部结构与种植体之间的冠根比会随之增加,种植体周围应力逐渐增加。当下颌骨垂直高度降低至一定程度时,种植体周围是否会出现破坏性的应力而导致种植体失败,目前还没有相关的文献报道。本实验设计了两种种植体支持式下半口全牙列固定义齿,研究在不同下颌骨垂直高度降低下种植体的应力分布情况。

1 材料和方法 1.1 下颌骨基准模型及义齿模型

下颌牙列缺失患者1例,男,60岁,下颌骨丰满,前牙区牙槽骨宽度大于6 mm,制作传统下半口义齿,义齿中加入阻射剂(硫酸钡与树脂基托材料的混合比例为1∶3)。将制作完成的下半口义齿戴入患者口中,在锥束计算机断层扫描成像(cone-beam computed tomography, CBCT)下,得到完整的下颌骨及上部义齿的影像图(图 1),包括皮质骨、松质骨、下颌神经管、颏孔以及上部义齿结构;义齿组织面与下颌骨之间的距离为黏膜厚度。将影像数据以DICOM格式导出[5]

图 1 戴有下半口义齿的下颌骨CBCT影像 Fig 1 Cone-beam computed tomography (CBCT) image of edentulous mandible and prosthesis

将DICOM格式数据导入ITK-SNAP软件(vesion 2.2, 宾夕法尼亚大学, 美国)进行影像处理, 形成三维标准三角语言(STL)表面.修剪STL表面构成骨小梁网格, 经过布尔运算构建并优化颌面部骨骼模型、无牙颌模型、下半口义齿模型, 将输出数据传送到HyperMesh 10.0软件中, 最后形成下颌骨基准模型及上部义齿结构模型(图 2).下颌骨基准模型包括皮质骨和松质骨; 下半口模型进行颊舌侧的基托边缘修整, 保留义齿及少量基托, 作为上部固定义齿模型.

图 2 下颌骨(A)以及上部义齿(B)的三维有限元模型 Fig 2 Three-dimensional finite element model of edentulous mandible (A) and prosthesis (B)

1.2 种植体模型

种植体采用柱形种植体外观,直径4.1 mm,长10 mm (ITI, SW)。根据种植体外形数据应用HyperMesh 10.0软件建立种植体模型。在本实验中忽略种植体螺纹。

1.3 建立不同垂直高度的无牙下颌三维模型

在下颌骨基准模型上,通过牙槽嵴顶下方间隔取像的方法,将牙槽嵴顶平面下移0、5、10、15 mm,形成4个不同程度下颌骨垂直高度的无牙模型,保留最上方牙槽嵴弧度。

1.4 建立包含种植体及上部结构的无牙

下颌三维模型分别采用以下两种设计模式,应用HyperMesh 10.0软件,在4个不同垂直高度无牙下颌模型上按以下两种设计模式装配种植体和上部义齿。模式1(6-种植体模型):颏孔前6枚相互平行的种植体,支持12颗人造牙的固定桥修复体,表示为模型6-0、6-5、6-10、6-15;模式2(4-种植体模型):颏孔前4枚种植体,两末端种植体向远中倾斜30°,支持12颗人造牙的固定桥修复体,表示为模型4-0、4-5、4-10、4-15。共建立8个三维有限元分析模型(图 3)。

1.5 无牙下颌三维模型材料性能及加载方式的设定

将模型中下颌骨、种植体和上部结构均考虑为连续、均质、各向同性的线弹性材料。弹性模量及泊松比[6-8]表 1,分别对下颌骨骨皮质与骨松质进行材料赋值,材料变形为小变形,种植体与骨界面为100%骨性结合,加载不发生相对滑动。边界条件:在下颌骨双侧升支后缘及下颌骨下缘咀嚼肌附丽区给予刚性约束,阻止下颌骨的刚体移动[9-11]。加载力参考现有的研究报道[5-7],设定为在右下颌第一磨牙中央窝加载250 N的垂直力,方向平行于种植体长轴并指向下颌骨。对模型进行加载设置后,用Ansys 15.0软件进行应力分析。

图 3 无牙下颌三维模型设计分组 Fig 3 Different designs of the edentulous mandible models

表 1 下颌骨有限元模型各结构的弹性模量及泊松比 Tab 1 Young's modulus and Poisson's ratio assigned to the mandible model components

2 结果

在相同的加载条件下,种植体有限元模型应力云图(图 4)显示,8个模型中越靠近加载侧的种植体(R5)承受的应力越大;种植体的最大应力值位于种植体颈部。

两组8个模型靠近加载侧种植体的最大应力值分布(表 2)显示,4-种植体组,模型4-0、4-5、4-10、4-15最大应力值分别为43.76、40.12、44.26、49.06 MPa;6-种植体组,模型6-0、6-5、6-10、6-15的最大应力值在83.13、80.62、88.88、109.64 MPa;6-种植体组的最大应力值是4-种植体组的2倍,显著大于4-种植体组。

两组模型中,种植体最大应力值在下颌骨垂直高度降低5 mm时表现为最低,随着下颌骨垂直高度的降低,应力逐渐增大;降低0 mm时的最大应力值要大于降低5 mm时。

图 4 8个有限元模型的种植体应力云图 Fig 4 The von Mises stress distribution (in MPa) within implants (8 models)

表 2 R5种植体颈部最大应力值 Tab 2 Maximum stress of implant neck of R5

3 讨论

本实验模拟了牙槽骨不同程度萎缩的下颌骨模型,比较了两种修复设计方案对种植体应力分布的影响。

对于下颌无牙的种植体支持式固定修复体,支持的种植体数量可由3~12枚不等,通常设计4~6枚种植体。在保证对义齿有效支持和固位、确保义齿功能有效发挥、无损种植义齿周围组织健康的前提下,种植体的数量越少越好。减少种植体数量,可以减轻患者经济负担,减少手术创伤和风险,降低上部结构制作和戴入的难度[12]。近年来,有学者提出在4个种植体上完成重建的理念(All-on-four),在无牙下颌前部颏孔间植入4枚种植体,用固定义齿修复,临床成功率高[13-16]。All-on-four要求双侧远中种植体平台向远中倾斜30°,以降低远中悬臂长度。本实验比较了两种设计方案,其中4-种植体模型按“All-on-four”要求设计。从实验结果来看,无论下颌骨的萎缩程度如何,4-种植体组承载的应力明显低于6-种植体组。虽然种植体的数量减少,但两端的种植体向远中倾斜后,使上部义齿的游离悬臂缩短,从而使义齿的支持结构更加合理,相应的种植体及周围骨应力分布也相对合理。由此可见,游离悬臂的长度对种植体周围应力分布有较大的影响[17]。因此,在牙槽骨严重萎缩的情况下,采用4-种植体的修复方式相对合理。

对于种植体冠根比的研究目前存有较大分歧。文献报道,随着冠根比的增加,种植体周围特别是种植体颈部的应力增加;但也有报道短种植体支持的固定义齿修复也能获得良好的临床效果[10-11]。本实验设定的条件是,当下颌骨垂直高度下降时,垂直距离保持不变,冠根比例随之增加。在下颌骨垂直高度下降5~15 mm范围内,4-种植体组和6-种植体组的种植体颈部应力均随着下颌骨垂直高度下降而增加,表明随着冠根比的增加,从加载点到种植体颈部的力臂增加,导致种植体颈部应力增加。而对于基准模型种植体颈部应力略高于下颌骨垂直高度降低5 mm组的现象,分析其原因可能是加载负荷通过上部义齿结构直接作用于种植体的同时,上部结构材料对负载也有一定的吸收和缓冲的作用,当在力臂较短的情况下,此结构上部结构对于负载的吸收或缓冲作用最小,因而加载在种植体上的负载就相对较大,所以基准模型虽然冠根比较小,却表现出较大的应力。随着基台高度的增加,上部结构吸收负荷的影响小于力臂长度增加的影响,冠根比因素起了主导作用。因此,种植体颈部应力表现出随着冠根比增大而逐渐增大的规律,但其变化是非线性的,这也从另一角度证实了上部结构吸收负荷这一因素的存在[18]。综上所述,本实验显示下颌骨垂直高度下降5 mm时,种植体应力分布更加合理。从临床角度看,在下颌骨有适当萎缩的情况下,一方面可以增加口腔内的垂直修复空间,另一方面也可以使种植体应力分布更加合理,延长种植体的寿命。

咬合力与种植义齿支持结构的稳定性密切相关,人的咬合力存在相当大的差异,咬合力大的人需要更加坚强的支持结构。Correa等[19]测定种植体支持的无牙颌修复体平均最大咬合力前牙区为166 N,磨牙区为240~248 N,前磨牙区为178~182 N。本实验采用单侧施加负荷的方式,加载力为250 N,旨在研究在异常极端咬合情况下,种植体是否会出现破坏性的应力。实验结果显示,4-种植体组的最大应力值为49.06 MPa,6-种植体组的最大应力值为109.64 MPa。而种植体钛金属的抗压强度与抗拉强度均为616 MPa,疲劳强度为364 MPa[20]。因此无论是4-种植体组还是6-种植体组,随着下颌骨垂直高度的下降,虽然种植体颈部应力增加,但并未出现破坏性应力。上述结果说明两种下颌骨种植体支持的固定义齿修复方式都是可行的,但从两组种植体的应力值来看,4-种植体的设计方案更加合理。

实验中,在相同加载条件下,各模型最大的应力均集中在距离加载点最近的种植体上,而且最大应力值均集中在舌侧区。从理论上分析,由于牙槽骨的结构在颊舌向具有非对称性,颊侧皮质骨密度和厚度都要大于舌侧,所以种植体周围骨应力分布在颊侧和舌侧也有明显不同。这种差异随骨密度的降低被明显放大,加剧了一侧的应力集中度。在咬合力作用下,种植体必然向刚度相对低的一侧倾斜,轴向加载就转化为力矩作用。一侧骨密度越低,倾斜程度加大,力矩作用就更加明显。因此,应力在种植体颈部集中的程度就会加剧[21-22]

在临床种植体支持的下颌半口义齿病例中,松动或者中央螺栓断裂的情况往往出现在远端的种植体。本实验结果发现在加载单一垂直力的情况下,所有的模型组的最大应力值都表现在远端的种植体上,而且集中于颈部皮质骨区域。皮质骨弹性模量远大于松质骨。在相同的形变位移时,弹性模量越大,应力就越大,颈部皮质骨就可能发生吸收,从而导致临床种植体松动[23]。但本实验中加载力只有垂直向,没有考虑侧向力,且实验只能有条件地模拟临床情况,在现实条件下咬合力情况更加复杂。牙槽骨吸收的原因也很复杂,其中的力学机制还有待进一步深入研究。

本实验采用逆向工程法建立了完整的下颌骨、固定义齿、种植体模型,相比传统方法得到更多的细节以及复杂的几何结构,该模型可模拟临床情况,采用三维有限元方法进行应力分析研究。三维有限元方法具有很多优点,有限元是目前最成功建立仿生数学模型的方法之一,可利用其对结构、形状、载荷和材料力学性能等极其复杂的构建进行应力分析,并已广泛应用于口腔生物力学研究领域[20]。但三维有限元法也有其局限性,比如限定所有的修复材料以及下颌骨均为各向同性的线弹性材料。这与临床实际情况有偏差,而且在实际生活中,义齿的受力情况复杂。因此,本实验中应力的绝对数值不是研究关注的重点,只是一个参考。本研究所关注的是各个分析模型应力间的相对数值,以此作为选择两种修复方案的参考[24]

综上所述,4-种植体和6-种植体支持的无牙颌固定义齿在相对极端的载荷下均未对种植体产生破坏性应力,4-种植体具有更加合理的应力分布;随着下颌骨垂直高度的下降,种植体周围的应力改变有一定规律,其中下颌骨垂直高度降低5 mm组应力最小,说明适当降低牙槽骨高度有利于种植体应力更加合理的分布。

参考文献
[1] BATENBURG R H, MEIJER H J, RAGHOEBAR G M, VISSINK A. Treatment concept for mandibular overdentures supported by endosseous implants:a literature review[J]. Int J Oral Maxillofac Implants , 1998, 13 :539–545.
[2] VISSER A, RAGHOEBAR G M, MEIJER H J, BATENBURG R H, VISSINK A. Mandibular overdentures supported by two or four endosseous implants. A 5-year prospective study[J]. Clin Oral Implants Res , 2005, 16 :19–25.
[3] CHEE W, JIVRAJ S. Treatment planning of the edentulous mandible[J]. Br Dent J , 2006, 201 :337–347. DOI:10.1038/sj.bdj.4814041
[4] SALEH SABER F, GHASEMI S, KOODARYAN R, BABALOO A, ABOLFAZLI N. The comparison of stress distribution with different implant numbers and inclination angles in all-on-four and conventional methods in maxilla:a finite element analysis[J]. J Dent Res Dent Clin Dent Prospects , 2015, 9 :246–253. DOI:10.15171/joddd.2015.044
[5] DEVOCHT J W, GOEL V K, ZEITLER D L, LEW D. Experimental validation of a finite element model of the temporomandibular joint[J]. J Oral Maxillofac Surg , 2001, 59 :775–778. DOI:10.1053/joms.2001.24292
[6] HUANG H L, CHANG C H, HSU J T, FALLGATTER A M, KO C C. Comparison of implant body designs and threaded designs of dental implants:a 3-dimensional finite element analysis[J]. Int J Oral Maxillofac Implants , 2007, 22 :551–562.
[7] O'BRIEN W J. Dental materials and their selection. 2nd ed. Chicago:Quintessence[M]. 2002 .
[8] YANG H S, LANG L A, MOLINA A, FELTON D A. The effects of dowel design and load direction on dowel-and-core restorations[J]. J Prosthet Dent , 2001, 85 :558–567. DOI:10.1067/mpr.2001.115504
[9] CRUZ M, LOURENÇO A F, TOLEDO E M, DA SILVA BARRA L P, DE CASTRO LEMONGE A C, WASSALL T. Finite element stress analysis of cuneiform and cylindrical threaded implant geometries[J]. Technol Health Care , 2006, 14 :421–438.
[10] CRUZ M, WASSALL T, TOLEDO E M, BARRA L P, LEMONGE A C. Three-dimensional finite element stress analysis of a cuneiform-geometry implant[J]. Int J Oral Maxillofac Implants , 2003, 18 :675–684.
[11] CRUZ M, WASSALL T, TOLEDO E M, DA SILVA BARRA L P, CRUZ S. Finite element stress analysis of dental prostheses supported by straight and angled implants[J]. Int J Oral Maxillofac Implants , 2009, 24 :391–403.
[12] SANTOS L S, ROSSI A C, FREIRE A R, MATOSO R I, CARIA P H, PRADO F B. Finite-element analysis of 3 situations of trauma in the human edentulous mandible[J]. J Oral Maxillofac Surg , 2015, 73 :683–691. DOI:10.1016/j.joms.2014.10.014
[13] MALHOTRA A O, PADMANABHAN T V, MOHAMED K, NATARAJAN S, ELAVIA U. Load transfer in tilted implants with varying cantilever lengths in an all-on-four situation[J]. Aust Dent J , 2012, 57 :440–445. DOI:10.1111/adj.12002
[14] TAKAHASHI T, SHIMAMURA I, SAKURAI K. Influence of number and inclination angle of implants on stress distribution in mandibular cortical bone with All-on-four Concept[J]. J Prosthodont Res , 2010, 54 :179–184. DOI:10.1016/j.jpor.2010.04.004
[15] BEVILACQUA M, TEALDO T, MENINI M, PERA F, MOSSOLOV A, DRAGO C, et al. The influence of cantilever length and implant inclination on stress distribution in maxillary implant-supported fixed dentures[J]. J Prosthet Dent , 2010, 105 :5–13.
[16] BONNET A S, POSTAIRE M, LIPINSKI P. Biomechanical study of mandible bone supporting a four-implant retained bridge finite element analysis of the influence of bone anisotropy and foodstuff position[J]. Med Eng Phys , 2009, 31 :806–815. DOI:10.1016/j.medengphy.2009.03.004
[17] BEVILACQUA M, TEALDO T, MENINI M, PERA F, MOSSOLOV A, DRAGO C, et al. The influence of cantilever length and implant inclination on stress distribution in maxillary implant supported fixed dentures[J]. J Prosthet Dent , 2010, 105 :5–13.
[18] AKÇA K, AKKOCAOGLU M, CÖMERT A, TEKDEMIR I, CEHRELI M C. Bone strains around immediately loaded implants supporting mandibular overdentures in human cadavers[J]. Int J Oral Maxillofac Implants , 2007, 22 :101–109.
[19] CORREA S, IVANCIK J, ISAZA J F, NARANJO M. Evaluation of the structural behavior of three and four implant-supported fixed prosthetic restorations by finite element analysis[J]. J Prosthodont Res , 2012, 56 :110–119. DOI:10.1016/j.jpor.2011.07.001
[20] VAN STADEN R C, GUAN H, LOO Y C. Application of the finite element method in dental implant research[J]. Computer Methods Biomech Biomed Engin , 2006, 9 :257–270. DOI:10.1080/10255840600837074
[21] TEPPER G, HAAS R, ZECHNER W, KRACH W, WATZEK G. Three dimensional finite element analysis of implant stability in the atrophic posterior maxilla:a mathematical study of the sinus floor augmentation[J]. Clin Oral Implants Res , 2002, 13 :657–665. DOI:10.1034/j.1600-0501.2002.130613.x
[22] SERTGÖZ A, GVVENER S. Finite element analysis of the effect of cantilever and implant length on stress distribution in an implant-supported prosthesis[J]. J Prosthet Dent , 1996, 76 :165–169. DOI:10.1016/S0022-3913(96)90301-7
[23] GENG J P, TAN K B, LIU G R. Application of finite analysis in implant dentistry:a review of literature[J]. J Prosthet Dent , 2001, 85 :585–598. DOI:10.1067/mpr.2001.115251
[24] LIU J, PAN S, DONG J, MO Z, FAN Y, FENG H. Influence of implant number on the biomechanical behavior of mandibular implant-retained/supported overdentures:a three-dimensional finite element analysis[J]. J Dent , 2013, 41 :241–249. DOI:10.1016/j.jdent.2012.11.008