矫形力作用下下颌骨重建过程的数值模拟.PDFVIP

2004年8月 重庆大学 学报 Aug．2004 第27卷第8期 Journal of Chongqing University Vn1．27 No．8 文章编号：1000—582X(2004)08—0020—04 矫形力作用下下颌骨重建过程的数值模拟 严 波，刘海京，张培 源 (重庆大学工程力学系，重庆 400030) 摘 要：下颌骨在矫形力作用下重建过程研究，对矫治方案确定等具有重要的指导意义。作者首先 建立下颌骨与颞下颌关节的三维有限元模型，并以此为基础，依据Cowin重建规律，对在矫形外力作用 下下颌骨矫治过程中骨的法向重建进行了数值模拟。对重建前后数值模拟结果的比较表明，法向重建 过程对于下颌骨的表面应力分布和应力梯度的分布具有较大的影响，因而在下颌骨矫治方案的制定过 程中这一因素不可忽视。 关键词：下颌骨；颞下颌关节；重建；数值模拟 中图分类号：Q66 文献标识码：A 有限元方法在口腔生物力学的研究中已经得到了 者建立了下颌骨与颞下颌关节的三维有限元模型，进 广泛应用。一些研究者建立了三维下颌骨与颞下颌关 而分析了髁突的应力分布，忽略重建对密度分布的影 节有限元模型，并研究了在颏兜力作用下髁突表面的 响，应用Cowin重建方程对矫形力作用下骨的重建过 应力分布，并建立了包含肌肉的有限元模型 j，有的 程进行了数值模拟。 不仅建立了类似的模型，还给出了肌肉的坐标与横截 1 三维有限元模型的建立 面积，并分析了颞下颌关节在不同的咬合情况下的应 力分布 j。有的还建立了相应的各向异性的有限元 选择一名l5岁下颌发育正常的男孩，对其下颌骨 模型 j，还研究了在不同方向的颏兜力作用下髁突 进行水平断层CT扫描，共扫描38张。将扫描图片转 表面和关节盘的应力分布，得到应力分布与颏兜力的 入计算机，然后转化成下颌边缘点的坐标，建立了包含 大小和颏兜力与 Frankfort平面问的夹角有关的 下颌关节和肌肉在内的三维实体模型，如图1所示。 结论[71。 计算中松质骨、关节盘和肌肉被视为各向同性弹性介 下颌骨在矫治过程中，在颏兜力的作用下会发生 质，其力学性能参数取值见表1所列 j。在进行有 骨的重建。骨组织的重建与应力应变的关系研究已经 限元网格划分时，由于模型的几何形状不规则，采用四 取得较大进展。早在1892年Wolff就提出了骨是依照 面体单元对下颌骨、关节盘和关节窝进行了划分，肌肉 作用于它们的应力和应变而变化的观点，并认为骨小 则用只受拉不受压的杆单元模拟。模型总共划分 梁按最大主应力的方向排列。将应变作为重建刺激量 20 019个单元，4 438个节点。矫形力的大小为7．ON， 建立了一个实用的可供计算的活体骨组织应力与表面 与Frankfort平面夹角为40。，方向指向冠状面。 重建适应模型 表1 骨和关节盘的力学性能参数 v(x0，几，t)=C；j(几， 0)[ ij( 0，t)一 ( 0)](1) 这里V( ，n，t)为在位置 处表面外法线几方向上的 生长速度，Cii(凡， )定义为生长速率系数， ij( ，t)为 直角坐标系中的应变分量，而 ( 。)为生长平衡应 变，t为时间。方程(1)即为Cowin重建方程。 为了研究下颌骨在矫形力作用下的重建行为，作 · 收稿日期：2004—04—14 作者简介：严波(1965一)，男，重庆永川人，重庆大学副教授，从事工程力学与