正畸弓丝与托槽间摩擦力建模及试验研究

姜金刚, 吴殿昊, 王 磊, 张永德, 刘 怡, 姚 亮

(1.哈尔滨理工大学 先进制造智能化技术教育部重点实验室, 黑龙江 哈尔滨 150080;2.北京大学 口腔医学院 北京 100081)

根据世界卫生组织的研究，错颌畸形是口腔三大疾病之一.全球儿童和青少年错颌畸形患病率为56%，主要是由牙位不正确引起的[1-2]，它的发生将影响咀嚼和发音，容易引发龋齿、牙周炎和呼吸道疾病等[3].最有效和成熟的治疗方法是正畸矫治，如图1所示，金属托槽预先粘贴于牙齿表面，托槽相对于牙齿是固定的，利用结扎丝将托槽与弓丝捆扎，正畸矫治是通过正畸弓丝形变产生有效矫治力带动托槽及牙齿产生生理性移动[4-5]，但结扎丝不能完全保证二者的相对固定，托槽将沿着弓丝形变的方向产生相对滑动或滑动趋势，进而在接触面产生正畸摩擦力[6-8].研究显示临床上50%～60%的正畸力要用来克服正畸过程中产生的摩擦力，矫治力必须克服摩擦力，并且不损伤牙齿和周围组织的健康，进而获得牙齿与周围组织的力学平衡，才能使牙齿产生生理性的移动[8-9].因此，了解正畸治疗过程中产生的摩擦力的强度是促进牙齿最佳生物运动的关键[10].

Fig.1 Orthodontic treatment technology图1 正畸矫治技术

近年来，轻力矫治和数字化正畸的理念越来越受到正畸医生的关注[11-12]，大量学者对正畸摩擦力展开研究.Razali和Gómez-Gómez等[13-15]采用有限元法探究了正畸摩擦力在正畸过程中的影响.Tiziano Baccetti等[16]通过试验表明，低摩擦系统有助于释放较大的正畸力，保证正畸过程中牙齿移动效率.许多研究学者利用特定测量装置及测量方法，通过试验分别测得不同干湿条件、不同正畸弓丝-托槽组合及不同结扎程度等情况下弓丝与托槽间的正畸摩擦力，致力于探究正畸摩擦力的主要影响因素及变化规律，并寻求使得正畸摩擦力最小的弓丝-托槽组合及接触方式[17-24].但这种方式大都需要多组且大量重复试验，对比分析后才能得到最佳的组合方案，且在不同的环境中所得到的数据具有差异性，不能为正畸医师提供统一化的标准方法，很难确切地控制摩擦力大小.因此亟需建立正畸摩擦力的预测模型，通过计算来精准预测正畸摩擦力.周珊等[25]根据摩擦力测量试验结果推算出正畸摩擦力预测公式，正畸摩擦力预测公式能够粗略推测出正畸矫治力需要克服的正畸摩擦力大小，但未能考虑正畸弓丝与托槽间的几何关系、力学关系及物理参数对预测模型的影响；Yukio Kojima等[26]建立了牙齿正畸滑动力学模型，根据弓丝与托槽间的力学关系，以分力叠加的方式量化正畸摩擦力，但一般弓丝与托槽之间存在微小的接触角，即牙齿在矫正滑动时引起的倾斜槽沟与参与矫正的弓丝之间接触而形成角度，该研究未能考虑接触角度对摩擦力的影响.

综上所述，目前现有的正畸摩擦力预测模型未能全面考虑弓丝托槽间受力关系和弓丝托槽的物理参数等因素，缺乏标准的参数化表达机制，因此不能保证模型预测的精确度和通用性.针对以上问题，依据正畸弓丝与托槽间的几何关系及力学关系，考虑正畸弓丝-托槽间接触特点及多物理参数，提出一种基于分力叠加法的正畸摩擦力预测模型建立方法，根据弓丝在托槽内的滑动机制，搭建基于六维力传感器的正畸摩擦力测量系统，进行不同弓丝-托槽组合和不同接触角度下的摩擦力测量及分析，以证明所建立的正畸摩擦力预测模型的有效性.

1 正畸摩擦力分析及预测模型建立

1.1 正畸摩擦力影响因素及力学基础

正畸矫治器由托槽和弓丝组成，托槽粘贴于待矫治的牙齿上，正畸医师根据临床矫治需求将弓丝弯制成特定形状，用以产生矫治畸形牙齿所需要的正畸矫治力，然后将弯制好的弓丝放入托槽内并用结扎丝进行固定，以便对弓丝位置及姿态进行限定[5,8].弓丝需要通过托槽才能将矫治力作用于牙齿本身，弓丝与托槽必然会接触，当弓丝与托槽间发生相对位移或者有相对位移的趋势时，就会产生正畸摩擦力，其大小与弓丝运动方向相反，在正畸矫治过程中，正畸力需要克服12%～60%的正畸摩擦力才能作用到牙齿上[8]，若摩擦力增大，矫治器的有效性将降低，治疗时间延长.另外，摩擦阻力对作用于牙齿的力矩/力也产生影响，继而影响牙齿的旋转中心，出现牙齿无法移动和牙根吸收等并发症，临床矫治的周期会被严重延长.

影响正畸摩擦力的因素分为两大类：机械因素和生物因素[8,19].由于生物因素针对不同人群的个性化差异大，且生物参数很难定量评估，因此在本研究中致力于探究机械因素对正畸摩擦力的影响.如图2所示，机械因素主要为弓丝材料、截面形状和尺寸、托槽材料及宽度、槽沟宽度、弓丝与托槽的表面粗糙度、弓丝与托槽的接触角度以及结扎方法和松紧度等.

Fig.2 Factors affecting orthodontic friction图2 影响正畸摩擦力的因素

弓丝、托槽和结扎丝是影响正畸摩擦力的3个主要装置，利用结扎丝将弓丝固定在托槽槽沟内，在弓丝与托槽槽沟底面间产生垂直于托槽基底方向的压力，当弓丝形变时生成正畸力F，同时在弓丝与托槽间产生摩擦力Ff.通常，根据弓丝与托槽的相对接触情况可知正畸摩擦力由三种接触分力组成，作用原理如图3所示.弓丝与托槽表面接触摩擦力称为经典摩擦力Ff；由于牙齿倾斜或弓丝弯曲导致弓丝与托槽槽沟之间的倾斜达到一定的角度时，弓丝与托槽壁接触而产生阻力称为约束力FB；当托槽与弓丝之间的成角增大使得弓丝发生永久变形时，弓丝与槽沟之间产生的阻力称为刻痕阻力FS[27].

Fig.3 Force analysis between orthodontic archwire and bracket图3 正畸弓丝与托槽间受力分析

1.2 正畸摩擦力预测模型建立

通过对正畸摩擦力原理分析可知，在正畸矫治过程中，依据正畸弓丝与托槽间接触角大小，正畸摩擦力三种接触分力可能阶段性地存在或同时存在，因此提出基于分力叠加原理的计及接触角 θc的正畸摩擦力预测模型建立方法，避免了传统正畸弓丝与托槽间摩擦力预测模型中对于正畸弓丝与第i个托槽间的接触角 θc的忽略而造成预测模型的精度损失.在摩擦力建模过程中先对接触角 θc进行分析，然后依次对约束力FB、经典摩擦力Ff和刻痕阻力FS进行建模，最终总摩擦力FR为三者之和.

如图4所示，以上牙列为例，首先考虑当弓丝与托槽槽沟上下壁两侧同时接触时，通过对弓丝及托槽参数分析及受力分析可以得到，正畸弓丝与第i个托槽间的接触角为

Fig.4 Position relation between adjacent brackets and friction analysis between the archwire and brackets图4 相邻托槽间位置关系及弓丝-托槽间摩擦力分析

式中i表示托槽的个数，每个托槽对应所在牙列上的1颗牙齿(正常人牙列为14颗牙齿)，i=1、2、....14，bi为第i个托槽槽沟间隙宽度，Wi为第i个托槽作用在接触段弓丝的宽度，D为弓丝横截面直径.

对于处于正畸移动过程中的正畸弓丝与托槽之间的接触作用，牙齿倾斜和正畸弓丝弯曲导致正畸弓丝与托槽槽沟之间形成倾斜角，则正畸弓丝与第i个托槽上下壁的两边同时接触产生约束力，正畸弓丝与第i个托槽间约束力是由正畸弓丝与第i个托槽下壁直接接触产生的正压力和第i个托槽内正畸弓丝弯曲变形产生的力矢量和相加求得：

式中 θi是 在约束力方向上与第i个托槽间的夹角，cosθi是平行于正畸弓丝表面的分量，μi是正畸弓丝与第i个托槽间的动摩擦系数，并且

任意两个相邻托槽间的正畸弓丝的挠曲线方程表达为

式中M(x)为任意两个相邻托槽间的弓丝在点x处的弯矩，E为杨氏模量，I为面积惯性矩，对于圆丝D为圆丝直径，对于方丝c2为矩形丝截面上平行于沟槽表面一边的长度，c1为矩形丝截面上垂直于沟槽表面一边的长度.

如图5所示，A点为第i个托槽上壁右端与正畸弓丝的接触点，即力矩分析过程中，A点所在托槽位置作为固定铰支座，B点为第i+1个托槽左端边线与正畸弓丝的交点，即B点所在托槽位置为滑移支座；C点为A点到弓丝中点间的任意一点，D点为弓丝中点到B点间的任意一点.A点所在托槽上下两接触点将朝着接触处弓丝的切线方向移动，托槽中心将沿着弓丝的施力方向移动，除此之外，中心还将产生逆时针的转动；同样B点所在托槽中心及上下壁端点将朝着弓丝的施力方向移动.第i个托槽与第i+1个托槽之间正畸弓丝的长度为Si，即AB间的弓丝段长度，i=14表示图4中上牙列最后1个托槽，根据临床经验，此托槽主要起到固定弓丝末端和颌内及颌间牵引作用，其受力情况不定，因此暂不对该托槽进行摩擦力的研究，为保证定义的完整性，特殊地，S14=0.两托槽间水平间距近似为是第i个牙齿的托槽底边中点到第i+1个相邻牙齿的托槽底边中点的距离.

Fig.5 Motion analysis of special points on brackets and relative position of adjacent brackets图5 托槽上特殊点运动分析及相邻托槽相对位置

如图6所示，第i个托槽和第i+1个托槽之间的正畸弓丝中的里任意一点C的力矩平衡表达式为

式中：Ay是 点A处y方向的反作用力，MA是点A处的反作用力矩，M1(x)是第i个托槽和第i+1个托槽之间的正畸弓丝中的里任意一点C的力矩，点C的横坐标为x.

Fig.6 Moment analysis of orthodontic archwire图6 正畸弓丝力矩分析

对于第i个托槽和第i+1个托槽之间的正畸弓丝的其余部分，即部分的正畸弓丝，任意一点D的力矩的平衡表达式为

式中，M2(x)是第i个托槽和第i+1个托槽之间的正畸弓丝上的之间任意一点D处的力矩，Mθi为第i个托槽和第i+1个托槽之间的正畸弓丝发生偏转时所产生力矩.

x处弯矩表示为

使用奇异函数，得到方程：

第i个托槽和第i+1个托槽之间的正畸弓丝的刚度EI是恒定的，将(7)代入方程(4)表示为

对x做两次积分，得到转角方程和挠度方程u(x)的关系，分别为

通过使用边界条件 θ (x)|x=0=0，u(x)|x=0=0来确定方程(11)中的积分常数，计算得出C1=C2=0；通过使用θ(x)|x=Si=0，u(x)|x=Si=0边界条件确定方程(11)中的第i个托槽和第i+1个托槽之间的正畸弓丝在x=0处的支反力Ay：

将(12)代入(8)，通过使用边界条件 θ (x)|x=Si=0，u(x)|x=Si=0进行求解得到：

因此，第i个托槽和第i+1个托槽之间的正畸弓丝的转角方程和挠度方程表达式为

由于第i个托槽和第i+1个托槽之间的正畸弓丝发生偏转时所产生力矩Mθi是由第i个托槽内正畸弓丝弯曲变形产生的力与Mθ的力矩臂L相乘所得，力矩臂L的长度为并且在第i个托槽和第i+1个托槽之间的正畸弓丝的中点处的转角 θ (x)为角度 θi-的正切值，通过(14)可以求解第i个托槽内正畸弓丝弯曲变形产生的力：

则

可以得到：

因此，正畸弓丝与第i个托槽下壁直接接触产生的正压力和第i个托槽内正畸弓丝弯曲变形产生的力分别相加得到正畸弓丝与第i个托槽间约束力作用于托槽上的反力

将(19)代入(2)，可以得出正畸弓丝与第i个托槽间约束力的表达式为

首先，正畸弓丝通过结扎丝固定在第i个托槽槽沟内，弓丝与托槽槽沟表面接触，并且受到结扎丝施加的法向压力作为已知量可利用结扎丝与弓丝间加力方式确定.此时会产生法向摩擦力，根据第一摩擦理论，可知：

其次，弓丝置于托槽内时与托槽上下壁接触，弓丝与第i个托槽下壁直接接触产生的正压力由滑动摩擦力计算公式得出正畸弓丝与第i个托槽间滑动摩擦力预测模型的表达式为：

第i个托槽间经典摩擦力应为法向摩擦力与滑动摩擦力之和

当托槽与弓丝之间形成的倾斜角增大使得弓丝发生永久变形时，弓丝与第i个槽沟之间产生的阻力称为刻痕阻力刻痕阻力由与第i个托槽接触的正畸弓丝弯曲变形产生的力求得：

根据式(24)和(18)可知弓丝与第i个槽沟之间产生的阻力称为刻痕阻力的表达模型为

正畸弓丝同第i个托槽间的总摩擦力由正畸弓丝与第i个托槽间滑动摩擦力、正畸弓丝与第i个托槽间约束力和正畸弓丝与第i个托槽间刻痕阻力叠 加而成，因此正畸弓丝与托槽间总摩擦力表示为

即

2 正畸摩擦力测量试验

2.1 试验材料

本研究中主要探究弓丝截面形状尺寸及材料、托槽尺寸及槽沟宽度和弓丝与托槽间接触角度为影响正畸摩擦力的主要因素，所以应采用不同的托槽与弓丝组合和不同的接触角度进行测量.试验中使用的不锈钢正畸弓丝由中国仁云医疗器械公司制造，澳丝正畸弓丝由美国A J Wilcock公司制造，结扎丝采用杭州奥索医疗器械有限公司生产的直径为0.25 mm的不锈钢结扎丝.为了探究不同弓丝与托槽组合对正畸摩擦力的影响，摩擦力测量试验采用包括两种材料和两种截面形状的10组不同的弓丝与三种不同槽沟尺寸的托槽组合测量，弓丝代号分别为I-S1616、Ⅱ-S1622、Ⅲ-S1625、Ⅳ-S0014、Ⅴ-S0016、Ⅵ-S0018、Ⅶ-S0020、Ⅷ-A0016、Ⅸ-A0018和Ⅹ-A0020.托槽代号为：T35、T40和T45.上述的弓丝代号命名规则为S表示不锈钢弓丝，A表示澳丝，后面四位数字为弓丝截面尺寸，代号00加弓丝直径表示的是圆丝截面尺寸，例如S1622表示不锈钢方丝，截面尺寸为0.406 4 mm×0.558 8 mm；A0016表示澳丝，截面直径为0.558 8 mm[5-6].托槽代号规则如下：T表示传统托槽，后面两位数字表示托槽槽沟宽度，单位为英寸，例如：T35表示槽沟宽度为b=0.889 mm的国产传统托槽.作用在接触段弓丝的托槽宽度Wi对于不同牙齿的托槽略微差异，考虑到制造误差，因此可通过游标卡尺实测.针对本研究中采用的国产不锈钢弓丝和澳丝，同种材料已采用哈尔滨工业大学5569电子万能材料试验机[28]通过单轴拉伸试验获得弹性模量E，国产不锈钢弓丝：E=8.52×104MPa；澳丝弓丝：E=7.49×104MPa，与不锈钢托槽间摩擦系数分别为0.15和0.12.对于圆丝对于方丝

2.2 测量系统搭建

通过分析正畸摩擦力的产生机理可知，由于弓丝与托槽之间相互接触且弓丝不断形变，在矫治过程中出现相对滑动现象或产生相对滑动的趋势，因此在二者之间会产生摩擦力.所以基于弓丝与托槽间的滑动机制，采用有限滑动法对正畸摩擦力进行测量.如图7所示，测量系统主要包括：正畸矫治器(弓丝、结扎丝及上牙列相邻3个托槽)、万向夹持器、模拟牙列、牙列支撑平台、电机控制器、弓丝拉动滑台、六维力传感器、力采集器、游标卡尺和电子角度测量仪.试验中使用的传感器为合肥旭宁科技公司生产的六维力传感器，其量程为20 N，单通道最大采样频率为5 000 Hz，分辨率为0.01 N.模拟牙列以铜柱作为模拟牙齿，为了实现弓丝良好的拉动顺畅性且又能合理安置力传感器，模拟牙齿安装在模拟牙列槽沟曲率较小的端部位置，并通过100 g砝码施加约为1 N的结扎压力.力传感器与活动模拟牙齿固接，活动模拟牙齿底端悬空，中间托槽粘接于活动模拟牙齿上，相邻两侧托槽粘接于固定模拟牙齿上且上下边与模拟牙列安装面平行，从而可以使得弓丝段的变形较小，易于测量所需的 θi值.每次试验前通过游标卡尺测量，调整牙齿位置使得待测托槽

Fig.7 Orthodontic friction measuring device图7 正畸摩擦力测量装置

在进行不同弓丝和托槽组合正畸摩擦力试验时，经过多次测试观察中间托槽底边与水平线间倾斜角度为10°(以划线方式预先标记)，微调中间托槽位置可确保弓丝与托槽槽沟上下壁面同时接触，然后验证公式(27)，测量十组弓丝与不同托槽组合情况下正畸摩擦力的数值.在进行不同接触角度下正畸摩擦力测量试验时，统一采用T40托槽，粘贴时中间托槽与相邻托槽上下边平行，经过计算该托槽与试验所用弓丝的理论最大接触角不超过9°，因此设定四组接触角度值，即0°、3°、6°和9°，可以仅考虑弓丝与托槽槽沟单侧壁面接触，验证公式(28).通过改变两侧模拟牙齿上螺柱高度并微调万向夹持器，可以改变中间托槽与弓丝的接触角，测量十组弓丝在该四组接触角度情况下正畸摩擦力的数值.

在测量时将模拟牙列固定于牙列支撑平台，根据试验要求可上下调节牙列支撑平台，保证弓丝在抽丝端水平附近，用95%乙醇棉球对托槽和弓丝进行脱脂处理，以消除各类附着物对试验结果造成的影响.由于所建立的正畸摩擦力预测模型未考虑生物因素，因此试验在干态环境室温(24±2 ℃)下进行.采用电子角度测量仪确定或保证每次试验涉及到的角度值，将弓丝抽丝端用夹紧枪头夹紧，启动电机驱动滑台使弓丝沿抽丝端延长线方向平稳拉动，考虑到临床正畸过程中牙齿移动速度，为了更加真实模拟弓丝在托槽间滑动情况，弓丝的拉动速度应尽量降低，根据试验研究[29]可知，正畸弓丝的滑动速度处于5.0×10-7～1.0×10-1mm/s时，由于滑动速度变化对正畸摩擦力影响程度较小，因此设置滑台移动速度为0.1 mm/s，并在有限的滑动行程3 mm内进行测定，即可模拟正畸治疗的条件[30].通过力采集器获得Z轴力数据即为摩擦力数值.

2.3 结果及讨论

利用正畸摩擦力测量装置对正畸摩擦力进行测量，为了保证所获得的测试结果具有可靠性和通用性，对同种代号弓丝的测量试验取多组样本，接触角度的测量误差在0.1°以内认为是有效样本，每组试验使用同一根弓丝截取相同长度的3段进行3次重复试验，取数值稳定状态下力的最大值，对测得的3个正畸摩擦力取算数平均值作为该次试验对应的有效测量值.

2.3.1 不同弓丝和托槽组合试验

不同弓丝和托槽组合下正畸摩擦力的试验结果及理论结果如图8所示.在测定不同托槽和弓丝组合摩擦力时，始终存在一定接触角度，代号T35、T40和T45所代表的国内传统托槽槽沟宽度依次增大.图8(a)、8(b)和8(c)中对正畸摩擦力试验值和预测模型计算所得理论值进行了对比.代号Ⅰ、Ⅱ和Ⅲ的不锈钢方丝，如图8(a)所示，根据试验结果对比曲线可知，相同截面尺寸的方丝与托槽间的正畸摩擦力随托槽槽沟宽度增大而减小，相同托槽配合的方丝与托槽间正畸摩擦力随截面尺寸(截面积)的增大而增大，如图8(d)所示，试验值与理论值的误差率处于1.35%～7.35%以内；代号Ⅳ、Ⅴ、Ⅵ和Ⅶ的不锈钢圆丝，代号Ⅷ、Ⅸ和Ⅹ的澳丝圆丝，如图8(b)和图8(c)所示，根据试验结果对比曲线可知，相同截面直径的圆丝与托槽间的正畸摩擦力随托槽槽沟宽度增大而减小，相同托槽配合的圆丝与托槽间正畸摩擦力随截面直径的增大而增大，如图8(d)所示，对于不锈钢圆丝，试验值与理论值的误差率处于1.70%～9.75%以内，对于澳丝圆丝，试验值与理论值的误差率均处于2.10%～8.65%以内.

综合分析代号Ⅴ、Ⅵ、Ⅶ的不锈钢圆丝和代号Ⅷ、Ⅸ、Ⅹ的澳丝圆丝，如图8(e)所示，在与相同托槽组合时，不锈钢圆丝与托槽间的正畸摩擦力高于相同截面尺寸澳丝圆丝与托槽间的正畸摩擦力.对比分析代号Ⅴ、Ⅵ、Ⅶ的不锈钢圆丝和代号Ⅰ、Ⅱ、Ⅲ的不锈钢方丝，在受到托槽约束的情况下，截面积较小的圆丝所产生的摩擦力可以达到甚至超过截面积较大的方丝，考虑其原因为随着矩形弓丝长边长度的增加而短边不变，矩形弓丝与托槽槽沟接触面积基本不变，且随着截面积增大，约束力与托槽间底边夹角变化幅度减小，导致正畸摩擦力的增加幅度较小.

结合以上三类弓丝试验结果进行分析，在与托槽槽沟双侧接触的情况下，正畸摩擦力与托槽槽沟宽度成负相关，托槽槽沟宽度越小，θi越大，从而越大，根据式(27)可知正畸弓丝所受到的三种分力都增大，因此摩擦力就会越大；正畸摩擦力与弓丝截面尺寸成正相关，弓丝截面尺寸越大，即在式(27)中I值逐渐增加，约束力与刻痕阻力越大，因此摩擦力也会越大；但正畸摩擦力受托槽槽沟宽度的影响程度明显高于弓丝截面尺寸的影响，即 θi变化幅度对结果的影响超过I变化的影响.另外，与方丝相比，圆丝的截面尺寸的变化引起的正畸摩擦力的变化较为敏感；与圆丝相比，托槽槽沟宽度的变化引起方丝与托槽间的正畸摩擦力变化较为敏感.

Fig.8 Experimental results of different orthodontic arch wire-bracket combination orthodontic friction force图8 不同弓丝-托槽组合正畸摩擦力测量试验结果

2.3.2 不同接触角度试验

不同接触角度下正畸摩擦力的试验结果及理论结果如表1所示，在测定不同接触角度情况下正畸摩擦力时，始终采用T40托槽与正畸弓丝组合，设置了四种接触角度对十种弓丝分别进行摩擦力测量，代号Ⅰ、Ⅱ和Ⅲ的不锈钢方丝，根据表1中试验结果可知，相同截面尺寸的方丝与托槽间的正畸摩擦力随接触角度增大而增大，试验值与理论值的误差率处于1.90%～8.30%之间；代号Ⅳ、Ⅴ、Ⅵ和Ⅶ的不锈钢圆丝，代号Ⅷ、Ⅸ和Ⅹ的澳丝圆丝，根据表1中试验结果可知，相同截面直径的圆丝与托槽间的正畸摩擦力随接触角度增大而增大，对于不锈钢圆丝，试验值与理论值的误差率处于0.55%～8.55%范围，对于澳丝圆丝，试验值与理论值的误差率均处于0.70%～9.65%之间.对于以上两种试验结果，理论模型与实际模型存在的误差，可能是由结扎丝施力不均和弓丝拉动过程中接触角与所设定角度的微小改变造成的，但在实际矫治中，牙齿移动速度极低，以上问题将得以缓解，因此本文中误差率的大小在临床正畸中是允许的.

表1 不同接触角度正畸摩擦力试验数值、理论数值及偏差率Table 1 Experimental values, theoretical values and deviation rate of orthodontic friction at different contact angles

综合分析代号Ⅴ、Ⅵ、Ⅶ的不锈钢圆丝和代号Ⅷ、Ⅸ、Ⅹ的澳丝圆丝的正畸摩擦力，在弓丝组合的托槽参数不变的情况下，对于同一接触角度，不锈钢圆丝与托槽间的正畸摩擦力高于相同截面尺寸澳丝圆丝与托槽间的正畸摩擦力.分析代号Ⅰ、Ⅱ、Ⅲ的不锈钢方丝和代号Ⅴ、Ⅵ、Ⅶ的不锈钢圆丝的正畸摩擦力，在接触角为0°和3°时，矩形弓丝较圆形弓丝产生的摩擦力大，考虑其原因是矩形弓丝的有效转矩更大，对托槽施加的有效压力也大，但随着接触角的增加，不锈钢圆丝比不锈钢方丝摩擦力增加得更迅速，在逐渐趋近托槽约束的过程中，接触角在6°和9°时，圆丝出现了比方丝更大的摩擦力.

结合以上三组试验结果进行分析，在弓丝组合的托槽参数不变的情况下，接触角的增大将引起 θi的大幅度增大，因此正畸摩擦力与接触角度成正相关.另外，与方丝相比，圆丝的截面尺寸的变化对正畸摩擦力较为敏感.

综上所述，在满足正畸治疗效果的前提下，为实现轻力矫治的效果，选用摩擦系数较小的矫治器材料，确保牢固的弓丝结扎及均匀的施力方式，并尽可能采用较细的正畸弓丝或槽沟尺寸较宽的托槽，将托槽与正畸弓丝间接触角度调整至最小，更易于避免摩擦力的产生，另外，使用澳丝进行矫治更能减小正畸摩擦力的影响进而提高矫治效率.在托槽接触约束较大或需要特殊弓丝形态的情况下，同种材料中可首先考虑方丝矫治.

3 结论

正畸摩擦力的有效控制是促进牙齿最佳生物运动的关键，本文中针对影响正畸摩擦力的机械因素，根据弓丝与托槽的接触情况将正畸摩擦力分为三种接触分力，依据正畸弓丝与托槽间的几何关系、力学关系及物理参数，提出一种基于分力叠加原理的计及接触角度的正畸摩擦力预测模型建立方法，将弓丝与托槽间接触角度作为正畸摩擦力预测模型的影响因素，避免了因考虑因素不全而导致的精度丢失问题，从而有效提高了正畸摩擦力数学模型的精确度.根据弓丝和托槽间滑动机制提出了以有限滑动法进行正畸摩擦力测量，模拟弓丝在托槽间的相对位移搭建了正畸摩擦力测量系统，以影响正畸摩擦力大小的弓丝截面尺寸、槽沟宽度、弓丝与托槽接触角为变量进行摩擦力的试验测量，探究了正畸摩擦力随着关键参数的变化规律，将测量得到的试验数据与预测模型得出的理论数据进行分析，正畸摩擦力模型得出的理论数据与试验数据的误差率处于0.55%～9.65%之间，将来可利用该模型建立正畸摩擦力预测系统，通过理论计算与仿真的方式对个性化正畸摩擦学行为进行精确预测，从而辅助数字化正畸治疗，达到轻力矫治的效果.未来将进一步将生物因素考虑到正畸摩擦力预测模型中，尽可能模拟在口内的真实环境进行建模和试验.