医学毕业论文--Amira与UG及LGOR软件联合建立下颌磨牙三维有限元模型.doc

Amira与UG及ALGOR软件联合建立下颌磨牙三维有限元模型 作者：杨萍，张少锋，刘继鹏，周冰，孟庆岩，高婧 【摘要】 目的：建立人下颌Ｄ区５６７三维有限元模型. 方法：采用薄层CT扫描技术、医学影像三维重建软件Amira和Unigraphics NX造型软件以及有限元分析软件ALGOR相结合建立下颌567及其支持组织的三维有限元模型. 结果：建立了下颌Ｄ区５６７三维有限元模型，总节点为93 260个，六面单元数为167 111个. 结论：所建模型结构完整，空间结构测量准确度高，单元划分精细、能够较精确地模拟实体状态，为进一步的生物力学研究提供了基础. 【关键词】 有限元分析法 　　0引言 　　有限元分析法(finite element analysis method，FEAM)是一种理论力学的分析计算方法，近年来,随着电子计算机技术的飞速发展和各种功能齐全的软件开发, 已被广泛运用到口腔医学领域中用于各种复杂问题的力学研究.该方法的基础是模型的建立，如何提高模型的相似性是研究的重点［1］ ，只有提高模型的相似性，才能保证实验结果的准确性. 由于牙颌组织外形结构复杂，以往建模难度大，花费高，建模时问长. 我们将薄层CT扫描技术、医学影像三维重建软件Amira，Unigraphics NX以及有限元分析软件ALGOR相结合建模，对如何提高有限元模型精确性、加快建模的速度、简化建模过程进行探讨. 　　1材料和方法 　　1.1材料GX26O商用台式机(处理器Intel Pentium 43.0 G，内存2 GB，硬盘250 GB，操作系统Windows 2000 Professiona1)；PQ6000螺旋CT扫描机(美国Picker公司)；应用软件Amira4.0（美国TGS公司）；UGNX4(美国EDS公司)；有限元分析软件ALGORv19（美国ALGOR公司）；根据参考文献［2］，选择一牙列完整、咬合关系正常、无明显牙周疾患及牙槽骨吸收成年男性志愿者，作为下颌牙列及下颌骨三维形态数据测量的样本. 　　1.2方法 　　1.2.1CT扫描及原始数据的获取测试者取仰卧位，颏部抬高，头部固定，戴用预制咬合板，避免上下牙列接触. 用螺旋CT扫描机进行扫描，层厚0.２ mm，选择得到断层影像100张，将扫描图像在CT工作站转换为bmp格式文件，记录并存入计算机. 　　1.2.2ＣＴ图像的三维重建将得到的bmp文件输入医学影像三维重建软件Amira中，读取二维图像数据生成牙齿的三维点云图（图1），转换为可以在UG中输入的tiff格式,再将模型导入UG造型软件实体网格构建. 将多面体网格填充方式由点云形成粗略实体，以牙冠为特征形成XY轴，使用UG中的Imageware模块对齐坐标轴，进行点云拟合和形成髓腔. 平行切割点云形成相交线. 利用已得到的相交线，通过曲线命令形成牙体根部，用缝合命令封闭根尖端完成髓腔实体. 用布尔运算形成牙的实体，将模型保存为igs格式（图2）. 　　图1－图2　略 　　1.2.3三维有限元模型的建立将建立的三维空间模型数据导入有限元分析软件中，利用ALGOR中分析模块对模型进行自动网格划分，采用8节点六面体单元自动、依次划分各部分的单元，获得三维有限元模型. 根据要求设定材料特性、边界条件，假设模型中的各材料和组织为连续、均质、各向同性的线弹性材料；其近远中径及底部完全固定，受力时模型各界面均不产生相互滑动. 参考文献［3-5］ ，所用主要材料参数如表1所示. 　　2结果 　　建立了下颌567实体模型， 将模型导入有限元分析软件ALGOR后， 采用8节点六面体单元自动、 依次划分各部分的单元， 获得三维有限元模型(图3). 建成的模型总单元数167 111个， 节点数93 260个， 其中下颌5单元数为36 093个， 节点21 830个， 下颌6单元数为73 984个， 节点40 694个， 下颌7单元数为57 034个， 节点30 736个. 所构建的有限元模型能够以单独或组合的方式观察各组成部分的情况， 还可以按照不同的研究目的和要求进行删除、 旋转、 切割和添加， 方便个别修改和调用. 　　表1有关材料的力学参数　略 　　3讨论 　　3.1三维有限元模型建立的基础FEAM分析研究的基础是建立有限元模型，理想模型的建立是准确计算和分析的前提. 而二维截面的生成是模型与实体相似度高低的关键，CT扫描技术具有定位准确，准确再现复杂细微结构，数据精确及可重复使用等优势. 故我们选择此方法较好地再现了牙颌组织结构. 建模过程中，由于牙齿的表面形态及内部结构细微复杂，在建模过程中反映其真实形态难度很大，尤其是磨牙的建模更为复杂，以往建模过程中部分方法不能获得精细的解剖外形及内部结构. 面对模型的“简化”意味着几何