作者：王宗篪， 庄国铭， 谢培庆， 赖森财
【关键词】 牙种植；牙周膜；牙模型；有限元分析；牙应力
摘要: 目的 研究仿生种植牙、天然牙和骨性结合种植体3种模型骨界面上的应力分布，从理论上探讨种植体的新模型。方法 应用螺旋CT对上颌中切牙标本进行扫描，分别建立仿生种植牙、天然牙和骨性结合种植体模型，对3个模型施加同样的载荷，应用Super SAP93软件，计算出天然牙周膜和仿生牙周膜以及3种模型骨界面上的应力分布。结果 仿生牙周膜和天然牙周膜内外表面应力的变化趋势相同；骨性结合种植体模型骨界面上在颈部出现应力集中；仿生种植牙和天然牙2个模型骨界面上最大Von Mises应力从颈部到根端的变化比较均匀，没有出现颈部应力集中现象。结论 引入仿生牙周膜能起到天然牙周膜的生物力学功能，它能实现牙合力的分散和减缓作用，使仿生种植牙骨界面上的应力分布比较均匀，从生物力学相容性上说明了仿生种植牙模型的可行性。
关键词: 牙种植；牙周膜；牙模型；有限元分析；牙应力
ABSTRACT:Objective To research the stress distribution of the bone interface of the bionic dental implant, natural tooth and osseointegrated implant, and to study the new model of the implant in theory. Methods The shape data of the central maxillary incisor are obtained by helixCT, then the threedimensional finite element models of the bionic implant, natural tooth and osseointegrated implant can be set up, respectively.Stress distribution of the natural periodontal membrane, bionic periodontal membrane and the bone interface of the three models are calculated out by SuperSAP93, when the same load exerted to the three models. Results The stress change of the inside and outside surface of the natural periodontal membrane and bionic periodontal membrane is in the same tendency. The stress concentration appears in the cervix of the osseointegrated implant.The maximum Von Mises stress which is on the bone interface of the bionic implant and the natural tooth changes evenly from their cervix to their roots, the stress concentration doesn't appear in the cervix. Conclusion The bionic periodontal membrane has the same function as the natural periodontal membrane. It has the dispersing and cushioning action against stress, and it make the stress of the bone interface of the bionic implant distribute evenly.
KEY WORDS: dental implantation; periodontal lignment; dental models; finite element analysis; dental stress analysisbionic dental implant; periodontal membrane; finite element model; threedimensional finite element analysis; stress distribution
基于Branemark提出的“骨性结合种植牙理论”的人工种植牙技术[1]，现已在口腔修复科研与临床中得到广泛应用。采用三维有限元分析技术对种植体周围骨组织的应力分布进行分析与计算机模拟是牙种植体研究开发和临床应用的重要内容[2]，很多学者开展了种植体的材料、几何尺寸及形状对骨界面应力分布的影响的研究，推动了人工种植牙技术的发展 [36]。由于钛及其合金具有优异的生物和物理性能，把它作为种植体材料已较好地解决了种植体与牙周骨结合的生物相容性问题。但是，生物力学相容性问题还未在理论上和临床实践上得到彻底解决。种植牙与牙周骨刚性结合，无法解决牙合力的减缓与分散问题，过大或过小的应力均可引起周骨质的吸收或萎缩，从而导致人工种植牙的失败[7]。
为了研究种植牙的生物力学相容性问题，笔者从仿生学原理出发，引入仿生牙周膜概念。仿生牙周膜是模仿天然牙周膜的生物力学功能，希望它能起到对牙合力的传导和缓冲作用，使牙槽骨界面上的应力分布比较均匀，力图解决种植牙生物力学相容性问题。本工作把上颌中切牙制成标本，对标本进行螺旋CT横断层扫描，利用计算机图像处理技术和三维有限元软件Super SAP93，分别对天然牙、仿生种植牙和骨性结合种植牙建立起三维有限元模型[8]。对3个模型施加同样的载荷，应用三维有限元软件Super SAP93，计算出天然牙周膜和仿生牙周膜以及三种模型骨界面上应力分布，从理论上探讨种植牙的新模型。
为了研究种植牙的生物力学相容性问题，笔者从仿生学原理出发，引入仿生牙周膜概念。仿生牙周膜是模仿天然牙周膜的生物力学功能，希望它能起到对牙合力的传导和缓冲作用，使牙槽骨界面上的应力分布比较均匀，力图解决种植牙生物力学相容性问题。本工作把上颌中切牙制成标本，对标本进行螺旋CT横断层扫描，利用计算机图像处理技术和三维有限元软件Super SAP93，分别对天然牙、仿生种植牙和骨性结合种植牙建立起三维有限元模型[8]。对3个模型施加同样的载荷，应用三维有限元软件Super SAP93，计算出天然牙周膜和仿生牙周膜以及三种模型骨界面上应力分布，从理论上探讨种植牙的新模型。
1 材料和方法
1.1 建模
选择一颗近期拨除较为典型的成人上颌中切牙，标准参照文献[9]。所选的牙体完整无缺损，全长23 mm，冠长11 mm，冠宽8 mm，冠厚7 mm，颈宽7 mm，颈厚6 mm，根长12 mm。确定其中心点，采用复合树脂包埋形成25 mm×25 mm×25 mm的包埋块标本，且牙长轴与包埋块底面垂直。采用螺旋CT机（SOMATOM Sensation 4，德国Seimens公司）对标本进行断层扫描，断面与牙长轴垂直。扫描的层间距为0.5 mm，层厚0.5 mm（层间损失小，可忽略不计），得到46幅二维扫描断层图像，再按文献[8]建立起六面体为单元的天然牙、仿生种植体和骨性结合种植体三维有限元模型，3个模型分别有24656个节点、21390个单元。
1.2 材料参数与实验条件
材料参数见表1。将3个有限元模型中各种材料和组织假设为连续、均质、线性、各向同性的线弹性材料。材料受力为小变形。模型的边界条件简化为牙槽骨外周围固定约束。对3个有限元模型施加同样载荷，加载荷的部位为舌侧切1/3与中1/3交界处，加载方向与牙长轴呈30 °夹角，载荷为100 N。取Von Mises应力作为衡量应力水平的主要指标，对天然牙周膜和仿生牙周膜上的应力分布以及3种模型骨界面上应力分布进行比较分析。表1模型中各组成部分的材料特性（略）
2 结果
垂直牙长轴方向从颈部到根端部均分7个截面，截面1为颈部，截面4为牙根正中部，截面7为牙根端部。在加载荷作用下，通过Super SAP93有限元分析软件的计算，天然牙和仿生种植牙牙周膜内外表面上最大Von Mises应力分布、3种模型骨界面上最大Von Mises应力分布以及仿生膜取不同弹性模量下仿生种植牙模型骨界面上最大Von Mises应力分布见图1～4。
垂直牙长轴方向从颈部到根端部均分7个截面，截面1为颈部，截面4为牙根正中部，截面7为牙根端部。在加载荷作用下，通过Super SAP93有限元分析软件的计算，天然牙和仿生种植牙牙周膜内外表面上最大Von Mises应力分布、3种模型骨界面上最大Von Mises应力分布以及仿生膜取不同弹性模量下仿生种植牙模型骨界面上最大Von Mises应力分布见图1～4。
2.1天然牙周膜应力分布
天然牙周膜内表面上最大Von Mises应力沿牙长轴方向是递减变化的，高应力区出现在牙颈部区域，达到17.68 MPa，牙根中部和牙根端部区域应力都比较小了，牙根端部的应力大约是颈部的1/5，只有3.66 MPa，应力变化幅度达到Δσvom=14 MPa，说明牙周膜内表面应力沿牙长轴方向分布很不均匀（图1）。天然牙周膜外表面上的最大Von Mises应力沿牙长轴方向变化与内表面的情况有很大不同。沿牙长轴方向，从颈部到牙根中部的应力是缓慢减小的，到牙根正中部达到最小值后，应力随之往根端又缓慢增大。牙颈部应力为6.82 MPa，牙根正中部为2.9 MPa，牙根端部为3.2 MPa，应力变化幅度仅达到Δσvom=3.92 MPa，这说明了牙周膜外表面应力沿牙长轴方向分布比较均匀。此外还可看出，沿牙长轴方向，牙周膜外表面比内表面的应力要减小。颈部减小幅度最大，达61.43%。说明牙周膜对牙合力具有明显的分散和缓冲作用。
2.2 仿生牙周膜应力分布
当仿生膜弹性模量取天然牙周膜弹性模量的10倍值时，通过三维有限元法计算出仿生膜内外表面最大Von Mises应力沿牙长轴方向的分布（图2）。仿生膜内表面上高应力区也出现在牙颈部，牙根中部到牙根端部的应力也比较小，内表面上应力变化幅度达到Δσvom=13.38 MPa。仿生膜外表面上应力分布比较均匀，颈部最大Von Mises应力值为6.63 MPa，牙根中部应力最小为3.99 MPa，应力变化幅度仅为Δσvom=2.64 MPa。与天然牙周膜应力比较，这种弹性模量的仿生膜同样能起到应力分散和缓冲作用。
2.3 3种模型骨界面应力分布
仿生种植牙和天然牙两种模型骨界面上最大Von Mises应力沿牙长轴方向的变化都比较均匀，天然牙在骨界面上应力变化幅度为Δσvom=3.92 MPa，而仿生种植牙骨界面上的应力变化幅度仅为Δσvom=1.65 MPa。骨性结合种植牙模型骨界面上最大Von Mises应力沿牙长轴方向变化是递减的。从牙根中部到牙根端部的应力降低比较平缓，降低幅度只有10%左右，表明牙根中部到牙根端部的应力逐渐降至比较低的水平。从颈部至牙根中部的应力降低幅度比较大，达到63%，颈部附近区域出现应力集中现象，颈部处应力最大，达到23.11 MPa，它与天然牙相比该处应力增加幅度为238.86%（图3）。
2.4 仿生膜取不同弹性模量的仿生种植牙模型骨界面应力分布
仿生种植牙模型中仿生牙周膜取不同弹性模量其骨界面上最大Von Mises应力沿牙长轴方向的变化也不同。在颈部和根端部最大Von Mises应力随仿生膜弹性模量的增大而减小，当仿生膜弹性模量为689 MPa时，应力达到最小值，随后应力又随仿生膜弹性模量的增大而增大（图4）。可以看出，仿生膜弹性模量在600~689 MPa，仿生种植牙骨界面上的应力分布最均匀。
仿生种植牙模型中仿生牙周膜取不同弹性模量其骨界面上最大Von Mises应力沿牙长轴方向的变化也不同。在颈部和根端部最大Von Mises应力随仿生膜弹性模量的增大而减小，当仿生膜弹性模量为689 MPa时，应力达到最小值，随后应力又随仿生膜弹性模量的增大而增大（图4）。可以看出，仿生膜弹性模量在600~689 MPa，仿生种植牙骨界面上的应力分布最均匀。
3 讨论
天然牙模型与仿生种植牙模型骨界面上应力分布比较均匀，无应力集中现象。这是由于牙周膜或仿生膜对牙合力的分散与缓冲作用所致。骨性结合种植体颈部出现应力集中现象的结果，是由于骨性结合种植体无牙周膜对牙合力的分散与减缓作用，这可能引起骨性结合种植体牙槽骨破坏的重要原因之一，这一研究结果与其他学者的研究结果一致，也与临床所见的种植修复中骨吸收常由种植体颈部周围开始的现象相符[10，1113]。
为了避免骨性结合种植体应力集中现象，一些种植系统中采用了“应力吸收装置”(stress absorbing element,SAE)，希望能解决牙合力的减缓与分散问题，例如IMZ系统，在种植体与上部结构之间设置一层弹性模量与天然牙周韧带相似的塑料。但是，Rossen等学者研究结果表明，这样的“应力吸收装置”种植系统有应力减缓作用而无应力分散作用[14]。David等学者对采用弹性螺栓及钛螺栓的种植体进行了比较研究，结果表明弹性螺栓均未起到减缓牙合力向骨组织传导的作用[15]。对于仿生种植牙，由于仿生膜对牙合力的分散与减缓作用，使骨界面上的应力大大降低且分布比较均匀，从而可以避免应力集中现象。
取仿生牙周膜不同的弹性模量，仿生膜内外表面和骨界面上的应力分布都会发生改变。也就是说，应力分布与仿生膜的弹性模量取值有关。当仿生膜弹性模量在600～689 MPa，仿生种植牙骨界面上的应力分布最均匀。这结果可从天然牙模型与仿生种植牙模型的力学结构比较中得到解释。天然牙模型中的牙体弹性模量为18 600 MPa、牙周膜弹性模量为68.9 MPa、牙槽骨骨皮质的弹性模量为13 700 MPa，牙体和牙槽骨骨皮质的弹性模量在104 MPa数量级，牙周膜弹性模量在0.7×102 MPa数量级，因此，牙体与牙槽骨骨皮质之间有一层弹性模量小2个数量级的牙周膜能起到牙合力的分散和减缓作用。在仿生种植牙模型中，种植体和根套都由钛材料组成，弹性模量为105数量级，与天然牙模型的力学结构原理类比，连结种植体与钛根套之间的仿生牙周膜弹性模量应小2个数量级，则应在0.7×103 MPa的范围，即仿生牙周膜的弹性模量大约是天然牙周膜弹性模量的10倍左右才能对牙合力起到最佳的分散与减缓的效果。这也说明了弹性模量在600 MPa与689 MPa的仿生膜在骨界面上应力比较小且分布最均匀。
在仿生种植牙模型中，仿生牙周膜起到了天然牙周膜的生物力学功能，实现了牙合力的分散和减缓作用，使仿生种植牙骨界面上的应力分布比较均匀，从生物力学相容性上说明了仿生种植牙模型的可行性。
天然牙模型与仿生种植牙模型骨界面上应力分布比较均匀，无应力集中现象。这是由于牙周膜或仿生膜对牙合力的分散与缓冲作用所致。骨性结合种植体颈部出现应力集中现象的结果，是由于骨性结合种植体无牙周膜对牙合力的分散与减缓作用，这可能引起骨性结合种植体牙槽骨破坏的重要原因之一，这一研究结果与其他学者的研究结果一致，也与临床所见的种植修复中骨吸收常由种植体颈部周围开始的现象相符[10，1113]。
为了避免骨性结合种植体应力集中现象，一些种植系统中采用了“应力吸收装置”(stress absorbing element,SAE)，希望能解决牙合力的减缓与分散问题，例如IMZ系统，在种植体与上部结构之间设置一层弹性模量与天然牙周韧带相似的塑料。但是，Rossen等学者研究结果表明，这样的“应力吸收装置”种植系统有应力减缓作用而无应力分散作用[14]。David等学者对采用弹性螺栓及钛螺栓的种植体进行了比较研究，结果表明弹性螺栓均未起到减缓牙合力向骨组织传导的作用[15]。对于仿生种植牙，由于仿生膜对牙合力的分散与减缓作用，使骨界面上的应力大大降低且分布比较均匀，从而可以避免应力集中现象。
取仿生牙周膜不同的弹性模量，仿生膜内外表面和骨界面上的应力分布都会发生改变。也就是说，应力分布与仿生膜的弹性模量取值有关。当仿生膜弹性模量在600～689 MPa，仿生种植牙骨界面上的应力分布最均匀。这结果可从天然牙模型与仿生种植牙模型的力学结构比较中得到解释。天然牙模型中的牙体弹性模量为18 600 MPa、牙周膜弹性模量为68.9 MPa、牙槽骨骨皮质的弹性模量为13 700 MPa，牙体和牙槽骨骨皮质的弹性模量在104 MPa数量级，牙周膜弹性模量在0.7×102 MPa数量级，因此，牙体与牙槽骨骨皮质之间有一层弹性模量小2个数量级的牙周膜能起到牙合力的分散和减缓作用。在仿生种植牙模型中，种植体和根套都由钛材料组成，弹性模量为105数量级，与天然牙模型的力学结构原理类比，连结种植体与钛根套之间的仿生牙周膜弹性模量应小2个数量级，则应在0.7×103 MPa的范围，即仿生牙周膜的弹性模量大约是天然牙周膜弹性模量的10倍左右才能对牙合力起到最佳的分散与减缓的效果。这也说明了弹性模量在600 MPa与689 MPa的仿生膜在骨界面上应力比较小且分布最均匀。
在仿生种植牙模型中，仿生牙周膜起到了天然牙周膜的生物力学功能，实现了牙合力的分散和减缓作用，使仿生种植牙骨界面上的应力分布比较均匀，从生物力学相容性上说明了仿生种植牙模型的可行性。
参考文献：
[1]Branemark P I,Zarb G A,Albrektsson T. Tissueintegrated prostheses: osseointegration in clinical dentistry[M]. Chicago: Quintessence Publishing Co Inc, 1985.
[2]Geng J P,Keson B C Tan,Gui Rong Liu. Application of finite element analysis in implant dentistry: A review of the literature[J]. J Prosthet Dent, 2001,85:585.
[3]Yalcin C,Senay C. The effect of veneering materials on stress distribution in implantsupported fixed prosthetic restorations[J]. Int J Oral Maxillofac Implants, 2000,15:571.
[4]Plikcioglu, Haldun, Akca, Kivanc. Comparative evaluation of the effect of diameter, length and number of implants supporting threeunit fixed partial prostheses on stress distribution in the bone[J]. J Dent, 2002,30(1):41.
[5]Laurent P, Guy H,Michel B,et al. Two dental implants designed for immediate loading: a finite element analysis[J]. Int J Oral Maxillofac Implants, 2002,17(3):353.
[6]Chun H J,Cheong S Y,Han J H. Evaluation of design parameters of osseiontegrated dental implants using finite element analysis[J]. J Oral Rehabil, 2002,29(6):565.
[7]董玉英,董福生. 人工种植牙生物力学研究进展[J].中国口腔种植学杂志, 1999,4(1):39.
[8]王宗篪，谢培庆，庄国铭,等. 仿生种植牙三维有限元模型的建立[J]. 山东大学学报（医学版）, 2006,44(1):7781.
[9]樊瑜波,张晓峰,唐高妍. 生理载荷作用下上颌中切牙牙周膜应力分布的三维有限元研究[J]. 生物医学工程学杂志, 1999,16(1):22.
[10]Meijer H J A,Starmans F J M,Bosman F,et al. A comparison of three finite element model of edentulous mandible provided with implants[J]. J Oral Rehabil, 1993,20:147.
[11]李湘霞,韩科,卜奎晨,等. 下颌骨形态对种植体骨界面应力分布影响的研究[J]. 中华口腔医学杂志, 2002,37(6):13.
[12]Melo C,Matsushita Y. Comparative stress analyses of fixed freeend osseointegrated prostheses using the finite method[J]. J Oral Implantol, 1995,21(4):290294.
[13]Stegaroiu R,Sato T,Kusakari H,et al. Influence of restoration type on stress distribution in bone around implants: A threedimensional finite element analysis[J]. Int J Oral Maxillofac Implants, 1998,13:8290.
[14]Van Rossen I P,Broak L H,de Putter C,et al. Stressabsorbing elements in dental implants[J]. J Prosthet Dent, 1990,64:198.
[15]David C,Willian R. Comparison of stress transmission in IMZ implant system with polyoxymethylene or titanium intramobile element:a finite element stress analysis[J]. Int J Oral Maxillofac Implants, 1992,7:450.
