建立上颌第一磨牙非线性三维有限元模型.docx

《建立上颌第一磨牙非线性三维有限元模型.docx》由会员分享，可在线阅读，更多相关《建立上颌第一磨牙非线性三维有限元模型.docx（7页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、建立上颌第一磨牙非线性三维有限元模型 目的：建立上颌第一磨牙及其牙周支持组织非线性三维有限元模型.方法：运用3DSS对上颌第一磨牙标准模型进行扫描；运用逆向工程软件GeomagicStudio8生成三维实体模型；有限元软件ANSYS10.0中生成三维有限元模型.结果：建立了包含上颌第一磨牙、牙周膜、硬骨板、松质骨及密质骨的非线性三维有限元模型，共96875个10节点四面体单元，132838个节点.结论：建立的三维有限元模型有很高的几何相像性，并且结构完整，网格质量较好，为进一步的生物力学探讨奠定了基础. 0引言上颌第一磨牙是正畸传统支抗的重要组成部分，其在正畸临床中有着极其重要的地位.临床上矫

2、治安排的制定常涉及扩大或缩小上颌第一磨牙宽度及压低、扶正、远中移动上颌第一磨牙，所以对上颌第一磨牙生物力学性质的探讨至关重要.三维有限元分析是对生物力学进行探讨的重要手段，其在口腔生物力学探讨中已经得到了广泛的应用与发展1.我们运用3DSS结合逆向工程技术建立上颌第一磨牙及其牙周支持组织三维有限元模型，并对牙周膜参数进行非线性的超弹性设定，为进一步的生物力学探讨奠定基础.1材料和方法1.1材料采纳王惠芸等测量和统计的中国人恒牙大小形态数据扩大15倍制成的上颌第一磨牙标准模具（第四军医高校口腔医学院解剖生理教研室）；PC机：IntelPentiumD3.0GHzCPU，2G内存，WindowsX

3、P操作系统；三维扫描仪（ThreeDimentionalSensingSystem，3DSS）：彩色标准型（3DSS?STDC?0601，上海数造机电科技有限公司，第四军医高校口腔医学院修复学教研室）；GeomagicStudio8逆向工程软件（RaindropGeomagic,Inc）；ANSYS10.0有限元分析软件（ANSYS,Inc）1.2方法1.2.13DSS扫描调整3DSS摄像头的高度与焦距，将15倍大小的上颌第一磨牙标准模具置于摄像头前方约2m处，并与摄像头等高，以黑色非反光幕布为背景，运用3DSS自带的扫描软件3dsscolor分别对上颌第一磨牙标准模具的近远中面、颊舌面、牙合

4、面及根面进行扫描，扫描得到的点云数据以.asc格式保存.1.2.2三维实体模型的建立将扫描得到的6部分点云数据输入GeomagicStudio8逆向工程软件中，去除多余的杂点与坏点，运用Merge吩咐进行拼接，得到完整的上颌第一磨牙点云模型.利用点云通过逆向工程的方式重构模型，并对模型进行适当的修补及曲面优化，得到NURBS曲面，进而获得上颌第一磨牙CAD三维实体模型.运用scale吩咐将模型缩小15倍并进行正交化，得到位于坐标系中心、正常大小的上颌第一磨牙三维实体模型，将此模型转变为片面模型，利用offset吩咐沿片面的法线方向向外匀称扩展0.25mm得到牙周膜片面模型，通过形成NURBS曲

5、面最终得到牙周膜CAD三维实体模型.按上述方法将牙周膜模型向外匀称扩展0.25mm得到硬骨板CAD三维实体模型.将上颌第一磨牙、牙周膜及硬骨板的三维实体模型输出为.iges格式文件分别进行保存.1.2.3三维有限元模型的建立利用ANSYS10.0中的建模工具形成松质骨及密质骨（厚2mm）模型，将得到的实体模型输入到ANSYS10.0有限元分析软件中，以腭根长12mm的位置作为冠根分界线，利用divide吩咐以工作平面对模型进行切割，再运用布尔运算的overlap吩咐与glue吩咐对上述实体模型进行运算，得到包含上颌第一磨牙、牙周膜、硬骨板、松质骨及密质骨的完整三维模型.假设牙齿、硬骨板、松质骨

6、及密质骨为连续、匀称、各向同性的线弹性材料；牙周膜假设为非线性的超弹性材料.受力时模型各界面之间不产生相对滑动.根据表1设定单元类型、尺寸及材料参数，进行网格划分，最终得到包含5个部分、4种材料的上颌第一磨牙及其牙周支持组织的非线性三维有限元模型.1.2.4材料参数及牙周膜的非线性设定模型各部分设定的单元类型、尺寸及材料参数如表1.牙周膜设定为非线性的超弹性模型.Vollmer等2探讨得到牙周膜双线性应力应变曲线（图1）：牙周膜在应变量达到=7.5%之前，其应力应变关系体现在第一个应力应变区，此时弹性模量E1=0.05MPa；当应变量达到=7.5%之后，应力应变关系体现为其次个应力应变区，其弹

7、性模量变为E2=0.22MPa.即牙周膜的应力应变不成线性关系，当应变量达到肯定限度时（=7.5%），其力学性质发生改变，弹性模量增加，较小的应变量即会产生较大的应力.模型柏松比=0.30.根据Mooney?Rivlin形式本构模型3进行拟合，得到牙周膜的双参数超弹性模型：依据公式：G=E/3=2(C10+C01)C01/C10=E1K=d/2d=(1?2)/(C10+C01)得到牙周膜双参数超弹性Mooney?Rivlin模型：C10=0.0079，C01=0.0004，d=48.2表1牙齿、牙周支持组织的三维有限元单元及力学参数（略）图1牙周膜双线性应力应变曲线（略）2结果建立的上颌第一磨

8、牙及其支持组织的非线性三维有限元模型具有很高的几何相像性，结构完整，网格质量较好.建立的三维有限元模型共96875个10节点四面体单元，132838个节点：上颌第一磨牙40415个单元，58605个节点；牙周膜10258个单元，20602个节点；硬骨板12202个单元，24413个节点；密质骨4132个单元，8053个节点；松质骨29868个单元，44977个节点.牙齿全长20.1mm，腭根长12.0mm，近中根11.6mm，远中根10.9mm，近远中宽9.2mm,颊舌向宽8.4mm（图2）.a：整体模型；b：牙齿模型；c：牙周膜模型；d：硬骨板模型；e：松质骨与密质骨模型.图2上颌第一磨牙及

9、其牙周支持组织三维有限元模型（略）3探讨上颌第一磨牙结构困难，形态不规则，建立完整、精确、相像性高的三维有限元模型对计算结果的精确性至关重要.3DSS是高速高精度的工业级三维扫描测量设备，采纳的是目前国际上最先进的结合结构光技术、相位测量技术、计算机视觉技术的复合三维非接触式测量技术，可以达到0.01mm的扫描精度，在数秒内即可完成对模型的高速高密度测量，输出三维点云供进一步后期处理.逆向工程是依据已经存在的产品模型，反向推出产品设计数据(包括设计图纸或数字模型)的过程是将产品样件转化为三维模型的相关数字化技术和几何建模技术的总称.GeomagicStudio8是四大逆向工程专业软件之一，利用

10、它可以在可视化的界面下对模型进行修改，大大缩短了建模时间，提高了建模的效率和可操作性.将3DSS与先进的逆向工程技术相结合，可以避开人为干扰，将扫描的点云数据干脆重建为三维实体模型，并且精确度极高，细微环节表达完整，相像性好，为计算结果的精确性供应了保证.但是由于其仅能得到模型的表面数据，应用范围受到了肯定程度的限制.试验中我们得到了以下阅历：运用Geomagic的offset吩咐与ANSYS的overlap吩咐，找到了一条建立困难牙齿牙周膜及硬骨板的有效便捷之路；在GeomagicStudio8中生成三维实体模型时，避开产生一些相交成锐角的线，以提高后期的网格划分质量；对模型不同的部分采纳不

11、同的网格尺寸进行网格划分，在确保计算精度的前提下可以提高运算效率；虽然.iges格式文件具有良好的软件接口兼容性，但在导入导出过程中仍会造成数据的丢失，所以要尽量避开模型以.iges格式反复的导入与导出.在众多探讨中，大多数学者常把牙周膜假设为均质、各向同性的线弹性材料6-7.但事实上，牙周膜是非均质、完全各向异性的非线性材料，这在少数学者的探讨中得到了部分的考虑.牙周膜的非线性主要表现为粘弹性与超弹性，粘弹性主要是探讨肯定加载作用下，在牙周膜达到稳定状态前的这段时间内，其应力应变随时间的改变关系8-10；当达到稳定状态后，应对牙周膜进行超弹性设定，以进行非线性应力分析.本模型的建立主要是为了探讨牙周膜达到稳定后的应力分布状况，所以对牙周膜进行了超弹性假设.我们在前人探讨的基础上尝试运用3DSS进行扫描获得建模原始数据，并对模型进行细化，建立了包括牙齿、牙周膜、硬骨板、松质骨及密质骨的三维有限元模型；同时，在肯定探讨假设的状况下对牙周膜进行了非线性参数设定，在肯定程度上提高了模型的仿生性，为进一步的生物力学探讨奠定了良好的基础.