
图1 物体变形前后图像
Fig.1 Image of object before and after deformation
如果取P为研究点,以P为中心取子集A,假设为m×m个像素,则A就记录了P点周围随机分布的斑点灰度值的信息,移动后,原来子集A处的斑点,就位于子集B处相应的位置,斑点间一一对应,这是另一样本空间。由概率和统计理论知道,两样本空间A,B完全匹配即相关,则相关系数为1,若有变形,相关系数会下降,则相关系数就小于1。若取一个非相应位置处的子集B’,则B’和A就不完全相关,其相关系数更小或为0[3]。设f(xi,yj)为变形前图像子集A中(i,j)处的光强,g(x*i;y*j)为变形后图像子集B中对应点坐标(x*i;y*j)的光强,f,g-D, 为计算子集区域的平均灰度值。则子集A和子集B的相关系数可用标准化协方差相关函数来估计:

可以由A、B的相关来确定某点的P的位移及其导数,只要两者相关,对于相关系数分布的C(u,v)曲面,根据相关系数取最大值时的位移就是真实的散斑位移[4]。因此由上述原理可以得到物体变形的位移信息,从而用数值方法得到相应的应变。
2 实验过程
2.1 试件的制备 试件主要分为两部分:一部分为中大医院提供的一组新鲜离体无龋、无损伤的正常后磨牙(浸泡于生理盐水中),包括上磨牙和下磨牙(男性的);另一部分是制备的牙周组织模拟体,由上海第二医科大学医药材料厂生产的自凝透明粉(20g)+(80ml)义齿基托树脂液剂Ⅱ型(又称:自凝牙托水,上海珊瑚化工厂生产),配置比例为1:4,放置于空气中自然固化,待液体呈粘稠状,稍微具有一定强度时,将后磨牙埋置于其中,模拟牙齿在牙周组织中的形态。待材料完全固化好后即可进行实验。实验后对这些正常的后磨牙使用桩核修复技术进行加工即可得到修复后的试件模型。
2.2 实验装置 实验系统由安置在固定结构上的白光光源、被测试件、CCD摄像头、砝码-杠杆加载装置、图像采集与处理设备(计算机)等,实验在防振台上进行。实验装置简图如图2。

图2 实验装置简图
Fig.2 Sketch of experiment installation
2.3 实验方法 先将试件定位在砝码-杠杆加载装置的测试台上,试件下缘与测试台台面均匀接触,防止下缘呈点接触,影响测量的位移和应变数值。然后对试件进行分级加载,每级荷载为25KN,加载等级为20级,加力方向拟正中咬时垂直受力状态。为了突出牙齿在牙周模型中的力学特性,提高测试精度,加载时加大了加荷力度。因为对于男性,每个磨牙在正常情况下,磨牙牙周膜所能耐受的咀嚼压力为350KN,但一般通过训练咀嚼压力可以有很大提高,且不会形成牙周组织的损伤,而这些强大的咬力迅速向临近牙槽骨扩散,以增加单位面积的承受能力,而不形成对牙周组织的创伤性力[5]。所以这里荷载加到500KN是可行的。
在定位好试件后,根据需要选择CCD的镜头、光圈和调整焦距,采用稳定、均匀光源,以使得到散斑图中的散斑颗粒清晰、大小合适。为了减少离面位移的影响、消除间隙,对试件预加微小初始力(杠杆及立柱的重量),稳定后采集第一幅图像,作为样本图像,然后分级加载,采集每级荷载下的图像各5幅,作为目标图像,从而得到整个加载过程的散斑图,图幅大小为768*576pixel。这些散斑图就是数字图像相关方法进行变形分析的原始数据,利用数字相关分析软件就可以得到整个加载过程的位移与应变演化情况。因计算结果是以像素为单位,为了将像素换算成长度单位,需要对记录系统进行标定,本实验采用简单的标定方法,在实验对象的同一位置处放置一把刻度尺,采集图像,找出像素数与单位长度之间的对应关系,就完成了标定。应该注意的是标定后对系统不能再做任何调整,否则需要重新标定。这里对修复前后的实验应分别进行标定。分析计算后得到的是各个状态的位移场和应变场,结果保存在数据文件里,然后用后处理软件(origin、matlab等 )处理和显示出来。
对正常后磨牙试件和桩核修复后的后磨牙分别进行上述实验,得到各自的处理结果以用于分析比较。图3为对应的桩核修复前后的上磨牙散斑图像:

(a)修复前 (b) 修复后
图3 桩核修复前后的上磨牙散斑图像
(a)Before restored (b)After restored
Fig.3 The speckle image of teeth before and after restored
3 结果与分析
由于实验均是对同一牙齿进行桩核修复前后的位移及应变进行测试,其试验的模型、条件、环境、方式等均一致,从而使得修复前后的测试结果具有最大的可比性。本研究是对修复前后的后磨牙在各级荷载下的散斑图进行处理,得到其位移场与应变场,再利用origin处理软件得出两者的位移对比曲线及变形对比曲线。这些曲线即可用来分析修复技术的可行性及修复的好坏。对于上磨牙,它有三个牙根,为了分析修复前后各牙根的受力情况,下面分别给出了修复前后各牙根相对于牙周组织移动的u、v场对比曲线。图4分别为正常后磨牙和桩核修复后的后磨牙各个牙根的u场曲线,图5分别为对应的修复前后的后磨牙各个牙根的v场曲线。

(a)修复前 (b)修复后
图4 桩核修复前后后磨牙各牙根u场对比曲线
(a)Before restored (b)After restored
Fig.4 Displacement of the teeth before and after restored in the x direction

(a)修复前 (b)修复后
图5 桩核修复前后后磨牙各牙根的v场对比曲线
(a)Before restored (b) After restored
Fig.5 Displacement of the teeth before and after restored in the y direction
由上面的曲线可以看出,修复前后的后磨牙各牙根随着荷载的增大在x、y方向的位移均呈增大趋势。注意的是为了评价修复技术的好坏,本文在分析结果时侧重于对曲线进行对比分析。
从上面曲线的对比中,可以清楚地看到修复后的后磨牙在受咬力下各牙根相对于牙周膜的移动情况与正常后磨牙基本保持一致,进而从另一角度说明修复前后的后磨牙对牙槽骨的作用力也保持一致,没有畸变。同时也进一步说明咬力在修复后的后磨牙各牙根的分配、传递情况均与正常后磨牙保持一致。
为了更有力的说明修复前后的后磨牙在荷载下的状态,下面给出了修复前后的后磨牙牙根的整体平均应变曲线对比图。图6为桩核修复前后的后磨牙在x、y方向的应变。

(a)x方向的应变 (b)y方向的应变
图6 桩核修复前后的后磨牙分别在x、y方向的应变
(a)Strain in x direction (b)Strain in y direction
Fig.6 Stain of the teeth before and after restored in x and y directions
由上面的曲线可以看出,修复前后的后磨牙在荷载下的应变曲线形状基本相同、大小偏差很小。这些曲线表明在初始的几级荷载下应变几乎为零,随着荷载的继续加大,开始出现弹性压应变,与实际情况较符合。
以上结论都说明桩核修复后的后磨牙其性能与正常后磨牙基本保持一致,不会影响正常咀嚼,且不会给患者带来不适。
4 小结
桩核修复残冠残根已被普遍使用,但桩核冠修复体在载荷作用下对牙周组织的影响机理仍然不是十分清楚。因此,本文利用图像采集系统采集修复前后的后磨牙在其支持组织模拟结构中受载时的图像,运用数字图像相关技术对其进行了分析,为桩核修复技术提供实验依据。
实验结果表明,对于桩核修复后的后磨牙,其力学特性与正常后磨牙保持一致,不会影响咀嚼,也不会给患者带来不适感,从而说明该桩核修复技术进行临床应用的可行性。
从该实验我们也可以看到数字图像处理技术运用于口腔修复力学研究的可行性及其优点,迈出了数字图像相关技术在口腔医学领域的应用第一步,同时也丰富了医学界关于核桩修复技术的现有成果。
【参考文献】
1 马少鹏,潘一山.岩石材料基于天然散斑场的变形观测方法研究.岩土力学与工程学报,2002,21(6):792-796.
2 单保华,欧进萍.散斑图像相关数字技术原理及应用.实验力学,2003,18:409-418.
3 王怀文,亢一澜,富东慧. 应用数字散斑相关技术进行薄膜材料断裂问题研究.烟台大学学报(自然科学与工程版).2001,14(2):100-104.
4 芮嘉白,金观昌,徐秉业. 一种新的数字散斑相关方法及其应用.力学学报.1994,26(5):599-607.
5 张富强,杨宠莹,苏炳华.力在牙槽骨内传递的力学分析.医用生物力学,1994,9:202-206.
中华现代中医学杂志
| • | 股骨头坏死早期和晚期疼痛的表现 | • | 新概念治疗股骨头坏死! |
| • | 吸鼻涕当心患中耳炎 | • | 慢性化脓性中耳炎中耳疾病 |
| • | 鼻窦内窥镜下治脑脊液鼻漏 | • | 长期的红霉素可治疗慢性鼻窦炎 |
| • | 慢性扁桃体炎需要做哪些检查 | • | 急性扁桃体炎 |
| • | 化妆品与青春痘 | • | 春季皮肤护理常识 |
| • | 精液是不是越多越好? | • | 男人一生精液有多少 |
| ........................更多>>> | |||