隐形矫治系统是三维测量技术、计算机技术以及快速成型技术完美的产物，使用隐形矫治系统可以为患者制作个性化透明牙套，实现舒适、隐形的矫治过程，并能模拟整个矫治过程，提前预测矫治结果。在矫治过程中，牙齿在矫治力的作用下位置发生改变，牙龈也随之产生变形，矫治结果的预测需要模拟牙齿移动的变化，包含牙齿位置变换和牙龈变形两个部分。

牙齿位置的移动属于刚体运动，通过坐标线性转换完成，而牙龈变形是牙齿移动仿真的难点之一，属于人体变形的范畴，人体变形是指人体在自身运动或受到外力作用时，骨骼（刚体）改变位置而导致的肌肉或皮肤等软组织发生形状改变。在隐形矫治系统中，牙龈变形仿真要在满足计算实时性的前提下保证良好的真实感。目前的口腔正畸系统，如Align公司的Invisalign系统、Cadent公司的OrthoCAD系统等都具有牙齿移动过程中牙龈跟随变形的仿真模拟。国内仅有首都医科大学的白玉兴等开发了国内首个无托槽矫治系统，并开展了临床试验，但从文献中还未有相关牙齿移动仿真相关技术研究的报道医`学教育网搜集整理。

牙齿移动过程中牙龈变形是牙齿移动仿真的关键部分之一，属于软组织变形模拟，主要有两种仿真方法：质点-弹簧建模法和有限元建模法，这两种方法都是将连续的几何体在空间上离散化，把无限问题转化为有限问题，通过求解相应的微分方程得到问题的近似解。质点-弹簧模型由于结构简易用、算法容易实现并且计算复杂度较低，己被广泛应用于许多领域。而有限元模型可伸缩性好，可以很方便地用相同的网格结构实现对不同复杂程度和精度的计算，而且模型的参数易于调节，能够方便地实现材质的各种属性。但是存在的主要问题是涉及大量复杂的计算，计算复杂度高。有许多学者对相似的软组织变形进行了研究。Nedel等使用质点-弹簧模型构造肌肉体，将拉格朗日动力学方程应用于该模型，然后使用龙格-库塔法来求解动力学方程，模拟了肌肉体的变形效果。Zhu等将8结点6面体为单元的线性有限元法应用于拉格朗日动力学方程，同样模拟了肌肉的变形。Wilhelms使用质点-弹簧模型模拟了皮肤的变形。Keeve使用有限元方法模拟了面部组织的变形。本文在分析牙颌模型数据表示方式的基础上，考虑到计算效率，采用弹簧质点模型仿真模拟牙齿的刚性移动和牙龈弹性变形过程，显示牙齿的虚拟移动。

1 材料与方法

1.1 材料与设备

本文使用国产3D CaMega 光学三维扫描系统采集牙颌石膏模型三维数据。国产3D CaMega 光学三维扫描系统利用精密机械系统精确移动定位实现对物体进行多方位、多角度的拍摄，生成全面、统一的三维型面点云数据。测量主要参数：单次拍摄范围30×24 mm；拍摄距离175 mm；X、Y分辨率0.023 mm；图像分辨率1280×1024；测量精度0.010 mm；测量密度130万点；快门0.4秒。本文使用的微型计算机配置：PIV 2.8G；512M内存；120G硬盘；GeForce2 MX440显卡；32M显存。选取翻制硬石膏模型开展试验。

1.2 模型数字化

全口牙列的三维测量提供了一种对患牙模型进行数字化的手段，为后续探测牙弓线、牙齿移动变形等提供了基础数据。本文采用光学三维扫描系统对牙齿石膏模型进行测量。测量头主要由一个投影装置和一个数码摄像机组成。测量时，将石膏备牙体装夹在工作平台上，投影装置投射一组光栅到石膏牙体表面，光栅的条纹会随着牙体表面的高度起伏而发生弯曲、变形，弯曲变形的程度包含了牙体的高度信息，通过反射由摄影机采集条纹和图像，经三维图像处理软件，对条纹图像进行处理，计算生成三维数据。

三维光学测量系统扫描获取的数据是石膏模型表面的点云数据，由于在测量过程中存在牙齿的重叠、遮挡等，需要进行多视角测量并对多次测量的数据进行拼接，对扫描得到的点云数据进行处理并三角化之后得到三角网格模型，由于在测量中不可避免的遇到噪声、扫描盲区等缺陷，需要对测量获取的模型进行去噪音、孔洞修补等处理，最终得到的点云数据和三角网格模型。

1.3 数字牙齿移动仿真

牙齿刚性移动可以采用线性坐标变换的方式进行，本文重点介绍由牙齿刚性移动产生的牙龈变形仿真模拟。

1.3.1 牙龈质点-弹簧变形建模

使用质点-弹簧可以为软组织建立面模型和体模型，分别由一系列平面片或多面体组成。在本文中，数字化后的牙颌模型是三角网格曲面模型，可以自然地把三角网格看作弹簧-质点系统网格，三角面片的顶点作为弹簧质点系统的质点，三角网格的边则作为连接质点的弹簧。这样，质点弹簧系统模型就将变形的牙龈简化成由线弹簧连接的线弹性质点系统，可以利用质点弹簧的运动规律来模拟牙龈的弹性变形过程。

在质点-弹簧系统中，每一个质点与周围相邻的若干个质点由遵守胡克定律的弹簧相连，质点运动则受到与其速度成正比的阻尼力的约束，质点运动而产生的弹簧变形力通过胡克定律计算。在典型的质点弹簧模型中，连接质点的弹簧包括拉力弹簧和交叉弹簧。拉力弹簧将一个质点和与它直接相连的质点连接，它产生的作用力抵抗压缩或拉伸变形；交叉弹簧跨过相邻三角形的公共，连接两个有公共边的三角形上的质点，交叉弹簧产生的力主要抵抗弯曲和剪切变形。为了提高计算的速度，本文使用简化的质点-弹簧模型，不考虑牙龈变形时产生的弯曲和剪切变形，简化的质点-弹簧模型。

1.3.2 系统动力学方程

在质点-弹簧模型中，质点所受的内力是弹簧产生的弹性变形力和引入的阻尼力。

本文使用Verlet积分法来求解质点-弹簧系统的动力学方程。Verlet算法并不显式的计算质点的运动速度，而是通过质点的位置隐式的计算速度，因此Verlet积分法相对比较稳定。

2 试验结果

2.1 试验工具

本文在WinXP操作系统平台上，利用Hoops图形开发包和VC++6.0软件开发平台自主研发了牙齿移动仿真模拟程序。

2.2 试验结果

牙齿移动仿真模拟计算的结果。模拟了牙齿平移2mm和6mm时以及相应的牙龈连接部分变形的效果。模拟了牙齿绕自身转轴旋转10和30时以及相应的牙龈连接部分变形的效果。

3 讨论

（1）透明隐形矫治技术涉及三维测量技术、数字化设计制造技术、口腔正畸修复学、口腔解剖学等多学科的交叉。由于透明隐形矫治器具有美观、舒适等优越性，在欧美等国已经广泛使用，而国内仅有少数口腔医学院校在开展研究，技术落后很多。迫切需要研发具有自主知识产权的数字化口腔正畸矫治技术。本文初步实现了牙齿移动、牙龈变形仿真模拟等关键技术，为透明隐形矫治器一系列母模的快速制作提供了方法和基础数据，具有重要应用价值。

（2）本文使用简化的质点弹簧方法模拟了牙龈随牙齿移动时的变形过程，使用Verlet积分法求解了系统的动力学方程，实现了由牙齿移动产生的牙龈变形的仿真。该方法在牙齿位置变动较小的情况下可以给出比较好的模拟结果，考虑到牙齿矫正过程中牙齿的变动位置一般比较小，因此该方法可以满足应用需求。

（3）该方法在求解质点运动方程需要耗费大量的时间，实时性受到一定影响，同时模拟的真实感有待改进。如何加快求解速度、如何实现实时性较好的变形仿真模拟，提高模拟的精度都是今后研究的重点。