本发明公开了一种T形曲动态矫治力作用下牙齿移动模型建立方法，涉及正畸治疗技术领域，它包含如下步骤：1)分析T形曲的结构特征和加载特征；2)建立T形曲竖直臂矫治力预测模型；3)建立T形曲水平臂矫治力预测模型；4)建立T形曲静态矫治力预测模型；5)建立蜡制颌堤模拟牙齿移动过程中动态阻力模型；6)建立T形曲动态矫治力预测模型；7)建立T形曲动态矫治力作用下牙齿移动模型。本发明能够有效的对医师所施加的T形曲矫治力值进行参数化表达，准确地预测医师所施加的T形曲动态矫治力值作用下牙齿的移动情况，辅助医生提高正畸治疗的安全性和预见性，提高口腔正畸诊疗的数字化程度。
一种T形曲动态矫治力作用下牙齿移动模型建立方法
技术领域
[0001] 本发明涉及一种T形曲动态矫治力作用下牙齿移动模型建立方法，属于正畸治疗技术领域。
背景技术
[0002] 固定矫治技术是目前最为有效的正畸治疗方法，T形曲是目前最为常用的关闭牙齿间隙的曲，在传统固定矫治技术的诊断过程中，各矫治阶段所使用正畸弓丝产生的矫治力大小、矫治效果的预测多根据正畸医生的经验进行判断，虽然这种依靠正畸医师经验的传统正畸治疗手段在大多数患者的治疗中能发挥一定的作用，但正畸力缺乏量化标准，治疗结果完全依赖于医生水平，易对患者造成伤害并导致治疗效率的降低。
[0003] 基托蜡具有质地较软，韧性好，加热变软不粘手等物理特性，在口腔临床应用中可基于患者牙齿数据制成蜡质颌堤，将蜡制颌堤浸入水浴箱可以将正畸矫治过程大大缩短，区别与其他刚性牙齿模型，将蜡制颌堤浸入恒温水浴环境中能够观察到牙齿随正畸弓丝的移动情况，实现正畸矫治过程中牙齿与牙周组织间相互作用这一动态过程的模拟，从而对矫治方案进行修正。
[0004] 虽然目前蜡制颌堤的应用能够在一定程度上提升矫治方案的合理性，但在应用基托蜡模拟正畸治疗的过程中，正畸弓丝形状与牙齿所受动态矫治力、动态矫治力作用下牙齿移动情况三者之间的关系并未确立，缺乏相应的量化标准，正畸医生仍旧无法通过该模拟过程预测在正畸弓丝影响下患者牙齿移动的情况，因此建立基于蜡制颌堤的T形曲动态矫治力作用下牙齿移动模型，对T形曲动态矫治力作用下牙齿移动情况进行预测，对于开展口腔数字化诊疗，辅助医生提高正畸治疗的安全性和预见性具有十分重要的意义。
发明内容
[0005] 针对上述问题，本发明要解决的技术问题是提供一种T形曲动态矫治力作用下牙齿移动模型建立方法，确立应用蜡制颌堤模拟正畸治疗的过程中，T形曲形状、牙齿所受动态矫治力、动态矫治力作用下牙齿移动情况三者之间的关系，辅助正畸医师设计正畸弓丝参数，从而提高正畸治疗的安全性和预见性。
[0006] 上述目的主要通过以下方案实现：
[0007] 本发明的一种T形曲动态矫治力作用下牙齿移动模型建立方法，其特征在于：所述方法的具体实现过程为：
[0008] 1)分析T形曲的结构特征和加载特征；
[0009] 2)建立T形曲竖直臂矫治力预测模型；
[0010] 3)建立T形曲圆弧部分水平臂矫治力预测模型；
[0011] 4)建立T形曲静态矫治力预测模型；
[0012] 5)建立蜡制颌堤模拟牙齿移动过程中动态阻力模型；
[0013] 6)建立T形曲动态矫治力预测模型；
[0014] 7)建立T形曲动态矫治力作用下牙齿移动模型。
[0015] 作为优选，所述的步骤1)中，由T形曲的结构特性可知，在T形曲对牙齿进行矫治时，矫治力由T形曲水平臂释放，T形曲发生变形的竖直臂及圆弧部分在T形曲的两侧相互对称，因此，在进行T形曲矫治力模型建立时，只需对对称的一部分进行分析；T形曲圆弧部分弯曲半径为R，整体高度为h，关闭间隙为b，关闭间隙由临床上后抽弓丝加力时产生；
[0016] 对矫治过程中T形曲的变形情况进行分析后可知，T形曲在竖直臂及T形曲圆弧部分水平臂发生变形，T形曲产生的正畸力由两变形区域产生的回复力组成，因此需对T形曲竖直臂及T形曲圆弧部分水平臂分别进行力学分析，并将二者叠加，建立T形曲矫治力预测模型；
[0017] 作为优选，所述的步骤2)中，T形曲竖直臂的转角方程θ(x)和挠度方程v(x)可表达为：
[0018]
[0019] 式中，M(x)是竖直部分上x距离处所受弯矩，E为材料的弹性模量，Iz为弓丝截面对z轴的惯性矩，对于圆丝Iz＝πd4/64，d为圆丝直径，对于矩形丝Iz＝c1c23/12，c2为矩形丝截面上与z轴平行边的长度，c1为矩形丝截面上与z轴垂直边的长度，C0和D0是积分常数，C0和D0由边界条件确定，T形曲竖直臂的弯矩方程为：
[0020] M(x)＝-P(y-x) (2)
[0021] 式中，P是产生竖直臂转角θ(x)变形所需的力，y为竖直臂变形前长度；
[0022] 将式(2)带入到式(1)中进行积分，可得：
[0023]
[0024] 为确定式(3)中的积分常数C0、D0，需要确定T形曲竖直臂的边界条件，基于支撑条件，其挠度或转角常为零或已知，对x＝0，即T形曲竖直臂及圆弧部分水平臂的连接处的变形量进行求解，由于T形曲竖直臂及圆弧部分水平臂的连接处有一纵向对称面，且外力均作用于这一对称面上，因此，变形后的曲梁轴线仍位于该纵向对称面内，该变形属于曲梁的平面弯曲变形问题，故可将连接处的圆弧等效为弧度为π/4的弯曲梁，对其取一段弧度为dα微元；
[0025] 在曲梁平面弯曲的情况下，外力都在曲梁的纵向对称面内，变形后的轴线仍为这一对称面内的曲线，曲梁并无扭转变形，此时仍可使用平面假设，可得到竖直臂及圆弧部分水平臂连接处曲梁变形后挠曲线的微分方程如式(4)，在给定的边界条件下，积分该方程便可确定曲梁的变形；
[0026]
[0027] 式中，u为竖直臂及圆弧部分水平臂连接处曲梁横截面在x方向上的位移，连接处所受转矩M0＝M|x＝0＝-Py，Iω为连接处曲梁横截面对ω轴的惯性矩，由于T形曲竖直臂与连接处曲梁的弯曲类型一致，有Iω＝Iz。
[0028] 由弧长公式可知ds＝Rdα，因此竖直臂及圆弧部分水平臂连接处曲梁变形后挠曲线的微分方程可变为：
[0029]
[0030] 解连接处曲梁变形后挠曲线的常系数非齐次微分方程求得：
[0031]
[0032] 由于竖直臂及圆弧部分水平臂连接处曲梁沿纵向对称面对称，竖直臂及圆弧部分水平臂曲梁的边界条件为 解得：
[0033] 连接处曲梁的挠度方程可表示为：
[0034]
[0035] 弯梁的转角方程为：
[0036]
[0037] 因此，有边界条件 解得：
[0038] 将C0，D0带入(3)中可得：
[0039]
[0040] 因为最大转角及最大挠度均产生在抽丝端处，即x＝y处，最大挠度即为抽丝端T形曲水平臂移动距离m，可知：
[0041]
[0042] 由作用力反作用力原则，T形曲竖直臂形变产生的矫治力F1即为产生竖直臂转角θ(x)变形所需力的反力：
[0043]
[0044] 作为优选，所述的步骤3)中，T形曲圆弧部分水平臂与T形曲竖直臂变形情况相同，因此可以使用相同的力学建模方法，在正畸治疗过程中，正畸医师会预先拉动T形曲抽丝端引起T形曲形变，并将形变后的T形曲安装在病人的牙齿上，由T形曲形变后产生的回复力拉动患者牙齿移动，达到矫治患者畸形牙齿的目的，T形曲圆弧部分水平臂的弯曲挠度，即T形曲圆弧部分水平臂沿Z轴方向的位移，在T形曲形变过程中，T形曲圆弧部分水平臂与T形曲竖直臂连接处的对称中心在空间中的位置是不断变化的，因此，为了计算T形曲圆弧部分水平臂弯曲挠度s，需对T形曲形变过程中连接处对称中心的运动情况进行简化，将变形前连接处的对称中心重合，计算变形后竖直臂理论长度与竖直臂变形前长度y的差值，即T形曲圆弧部分水平臂弯曲挠度s为：
[0045]
[0046] 可得到T形曲圆弧部分水平臂挠曲线近似微分方程为：
[0047]
[0048] 式中，M(l)为T形曲圆弧部分水平臂的l距离处的弯矩；
[0049] T形曲圆弧部分水平臂的转角方程θ(l)和挠度方程v(l)可通过对式(13)进行积分得到：
[0050]
[0051]
[0052] 式中，C1和D1是由边界条件确定的积分常数，而T形曲圆弧部分水平臂的弯矩方程可由式(16)表达：
[0053] M(l)＝-G(w-R-l) (16)
[0054] 式中，G为圆弧部分水平臂产生的矫治力，w为T形曲圆弧部分水平臂的长度；
[0055] 将式(13)和(16)代入到式(14)和(15)中进行积分，可得：
[0056]
[0057]
[0058] T形曲圆弧部分水平臂的边界条件求解与T形曲竖直臂相同，将一侧的侧面圆弧弧度为π/4，对其取一段弧度为dβ微元，由于对竖直臂边界条件的求解给出了详细的过程，因此，在这里直接给出解圆弧段挠度常系数非齐次方程得到的边界条件式；
[0059]
[0060] 由于T形曲侧面圆弧的对称性，侧面圆弧的边界条件为 解得：
[0061] 则弯梁的挠度方程为：
[0062]
[0063] 弯梁的转角方程为：
[0064]
[0065] 因此，有边界条件v|l＝0＝u|β＝0＝0，θ|l＝0＝ε|β＝0＝0，解得：C1＝0，D1＝0；
[0066] 将C1和D1代入(17)、(18)中可得：
[0067]
[0068]
[0069] 式中，T形曲圆弧部分水平臂最大转角及最大挠度产生在l＝w-R处，最大挠度即为T形曲圆弧部分水平臂(4-3)移动距离s，可知：
[0070]
[0071] 沿抽丝端方向的正畸力F2即为由T形曲圆弧部分水平臂形变产生圆弧部分水平臂转角θ(l)变形所需力沿抽丝端方向分力的反力，即：
[0072]
[0073] 作为优选，所述的步骤4)中，T形曲(4)产生的静态矫治力即为竖直臂(4-1)及圆弧部分(4-2)形变后产生的回复力的合力，即：
[0074]
[0075] 作为优选，本方法应用于一种模拟牙齿移动用蜡制颌堤。
[0076] 作为优选，本方法适用的一种模拟牙齿移动用蜡制颌堤由基托蜡颌堤、树脂牙齿模型、正畸托槽、T形曲组成。
[0077] 作为优选，所述的步骤5)中，被测牙齿与测量原件以树脂圆柱体连接，牙齿在蜡制颌堤中的移动实际为圆柱连接体在蜡制颌堤中的运动，因此，以圆柱体作为基本构件进行分析；当牙齿在蜡制颌堤中移动的速度为vt时，vt为蜡制颌堤在t时刻的流动速度，圆柱体上受到沿移动方向的作用力为绕流拖曳力；摩擦拖曳力和压差拖曳力共同组成绕流拖曳力；摩擦拖曳力是由于流体的粘滞性在柱体表面形成边界层，在此边界层范围内，流体产生速度梯度，摩擦效应显著，产生了摩擦切应力；压差拖曳力是由于边界层在圆柱体表面某点处分离，在分离点下游即在圆柱体后部形成很强的漩涡尾流，使得圆柱前后产生压力差，进而在流动方向产生了一个力，而在流体流动中，圆柱体旋涡尾流是随雷诺数的Re的变化发展的，牙齿在蜡制颌堤中的移动雷诺数Re＜5，因此无旋涡尾流的产生，无压差拖曳力产生；
[0078] 对单位长度圆柱体上的拖曳力fD可用式(27)计算：
[0079]
[0080] 式中，v0为未受绕流影响垂直于圆柱体轴线的牙齿移动速度分量，ρ(t)为在t时刻实验温度下蜡制颌堤的密度，A为单位长度圆柱体垂直于移动方向的投影面积，对于圆柱体，A＝1×D，D为圆柱体的直径，CD为拖曳力系数，它集中反映了流体的粘滞性而引起的粘滞效应，与雷诺数Re和圆柱面粗糙度Ra有关；
[0081] 假设本研究中的蜡制颌堤流体是不可压缩的理想流体，排蜡体积为 的圆柱体在移动速度vt＝v(x,y,z,t)的蜡制颌堤流场中移动；暂不考虑圆柱体对蜡制颌堤流场的影响，即假定蜡制颌堤流场内的压强分布不因圆柱体的存在而改变，那么圆柱体的边界作为加速流体边界的一部分，也就是被圆柱体置换的那部分体积内的蜡制颌堤流体，它本应以静止的状态存在于蜡制颌堤流场中，但实际上由于圆柱体移动的存在，这部分静止的蜡制颌堤流体将被加速至与圆柱体边界移速相同的状态；因此加速的蜡制颌堤流体将会对排蜡体积为 的圆柱体沿流动方向作用一个惯性力Fk，惯性力Fk的数值等于圆柱体的排蜡质量M0与体积 内蜡制颌堤流体的平均加速度 的乘积，即：
[0082]
[0083] 对于研究中的圆柱体来说， 可以取圆柱体轴中心位置处的流体加速度 来表示，此时：
[0084]
[0085] 但由于圆柱体存在于蜡制颌堤流场中，必将使圆柱体周围的流体质点受到扰动引起速度变化，从而改变蜡制颌堤流场内的压强分布，所以，圆柱体的扰动是圆柱体周围改变了原来运动状态的那部分附加流体的质量Mw沿流体流动方向也将对圆柱体产生一个附加惯性力，即附加质量力；因此加速的流体沿流动方向实际作用在圆柱体上的绕流惯性力fL可表示为：
[0086]
[0087] 令Mw＝CmM0，则式(30)可表示为：
[0088]
[0089] 式中，Cm为附加质量系数，CM为质量系数，亦称为惯性力系数，集中反映了由于流体的惯性以及圆柱体的存在，使圆柱体周围蜡制颌堤流场的速度改变而引起的附加质量效应；
[0090] 经过以上分析，可获得牙齿在蜡制颌堤移动过程中所受阻力情况，牙齿在正畸弓丝变形产生的正畸力影响下在蜡制颌堤中移动，移动过程中，由于流场的绕流特性，牙齿受到绕流惯性力fL和拖曳力fD的影响；
[0091] 受热量交换影响，热场中蜡制颌堤模型内部温度是随时间变化的，内部温度的变化将引起蜡制颌堤模型密度的变化，进而影响牙齿在蜡制颌堤内部移动受到的阻力。牙齿模型在蜡制颌堤内部移动时遵循粘性流体能量方程中的规律，令e代表单位质量流体所具有的内能，则ρ(t)e为单位体积流体具有的内能，ρ(t)vt2/2代表单位体积的动能，从而单位体积流体包含的总能量E＝ρ(t)e+ρ(t)vt2/2；
[0092] 经过简化整理，能量守恒原理可近似地表示为：
[0093]
[0094] 式中，cp为无量纲压强系数，Φ为牙齿模型在蜡制颌堤流体中移动时消耗的机械功，k为计算系数，▽Τ为基托蜡流体热场的温度梯度，q为热流密度；
[0095] 对▽Τ进行求解，设蜡制颌堤的厚度为2δ，初始温度为t0；在初始瞬间将它放置于温度为t∞的流体中，流体与蜡制颌堤间的表面传热系数h为常数，蜡制颌堤两边对称受热，蜡制颌堤内部温度分布必以其中心截面为对称面，因此，只需要研究厚度为δ的半块蜡制颌堤，把x轴的原点置于蜡制颌堤的中心截面上，对于x≥0的半块蜡制颌堤，可列出如下导热微分方程：
[0096]
[0097] 式中，a为热扩散率，式(33)两边对x积分，可得：
[0098]
[0099] 蜡制颌堤在水浴环境下均匀受热可以简化为一维热场问题，因此有:
[0100]
[0101] 将蜡制颌堤流体热场的温度梯度 代入式(32)中整理得到：
[0102]
[0103] 对式(36)等式两边对t积分整理得到：
[0104]
[0105] 式中，T为蜡制颌堤流体热场的温度；
[0106] 傅里叶定律用热流密度q表示时有如下形式：
[0107]
[0108] 式中，λ为导热系数；
[0109] 将式(38)代入式(37)中可得蜡质颌堤密度ρ随时间t变化的表达式：
[0110]
[0111] 蜡制颌堤模拟牙齿移动过程中动态阻力模型可由式(40)表达：
[0112]
[0113] 式中，f为蜡制颌堤模拟牙齿移动过程中的动态阻力。
[0114] 作为优选，所述的步骤6)中，基于T形曲静态矫治力预测模型及蜡制颌堤模拟牙齿移动过程中动态阻力模型，T形曲动态矫治力预测模型可表达为：
[0115]
[0116] 式中，F即为应用T形曲动态矫治力预测模型得出的T形曲动态矫治力。
[0117] 作为优选，所述的步骤7)中，物理学中基于加速度通过积分可以得到位移和速度矢量关系：
[0118]
[0119] 牙齿移动的加速度a0、牙齿的质量m0、作用力间的关系可由式(43)表达：
[0120] F＝m0a0 (43)
[0121] 则受到正畸力F0、绕流惯性力fL和拖曳力fD作用在蜡质颌堤中移动的牙齿的加速度a0由式(44)表达：
[0122]
[0123] 将式(44)代入式(42)中可得到基于蜡制颌堤的T形曲作用下牙齿移动预测模型：
[0124]
[0125] 式中，Sm为基于蜡制颌堤的T形曲作用下牙齿移动，m0为牙齿的质量，t1为牙齿移动的时间。
[0126] 本发明的有益效果为：
[0127] 1、采用参数化的建模方法，能够更直观地反映出各因素对T形曲作用下牙齿移动情况的影响效果，便于医师对弯制的弓丝进行调整，进而获得最佳的牙齿矫治效果。
[0128] 2、本方法适用于基托蜡质颌堤，相较于基于传统刚性颌堤得出的预测模型，本方法能够体现真实正畸过程中牙齿相对移动及正畸弓丝力衰减现象对矫治力的影响，所以本方法算得的矫治力能够体现真实正畸过程中的动态特性，在矫治力及牙齿移动预测方面本方法的精确度更高。
附图说明
[0129] 为了易于说明，本发明由下述的具体实施及附图作以详细描述。
[0130] 图1为本发明方法T形曲动态矫治力作用下牙齿移动模型建立流程图；
[0131] 图2为本发明T形曲曲竖直臂力学分析示意图；
[0132] 图3为本发明T形曲弯曲部分力学分析示意图；
[0133] 图4为本发明T形曲圆弧部分变形分析示意图；
[0134] 图5为本发明所使用的一种模拟牙齿移动用蜡制颌堤示意图；
[0135] 图6为本发明方法蜡制颌堤模拟牙齿移动过程中动态阻力模型建立流程图；
[0136] 图7为本发明蜡制颌堤中牙齿移动示意图；
[0137] 图8为本发明蜡制颌堤受热示意图；
[0138] 图中：1-基托蜡颌堤；2-树脂牙齿模型；3-正畸托槽；4-T形曲；4-1-竖直臂；4-2-圆弧部分；4-3-圆弧部分水平臂；4-4-T形曲水平臂。
具体实施方式
[0139] 为使本发明的目的、技术方案和优点更加清楚明了，下面通过附图中示出的具体实施例来描述本发明。但是应该理解，这些描述只是示例性的，而并非要限制本发明的范围。此外，在以下说明中，省略了对公知结构和技术的描述，以避免不必要地混淆本发明的概念。
[0140] 如图1、图2、图3、图4、图5、图6、图7、图8所示，本具体实施方式采用以下技术方案：本发明的一种T形曲动态矫治力作用下牙齿移动模型建立方法，其特征在于：所述方法的具体实现过程包括以下步骤：
[0141] 1)分析T形曲的结构特征和加载特征；
[0142] 2)建立T形曲竖直臂矫治力预测模型；
[0143] 3)建立T形曲圆弧部分水平臂矫治力预测模型；
[0144] 4)建立T形曲静态矫治力预测模型；
[0145] 5)建立蜡制颌堤模拟牙齿移动过程中动态阻力模型；
[0146] 6)建立T形曲动态矫治力预测模型；
[0147] 7)建立T形曲动态矫治力作用下牙齿移动模型。
[0148] 进一步地，所述的步骤1)中，由T形曲的结构特性可知，在T形曲对牙齿进行矫治时，矫治力由T形曲水平臂4-4释放，T形曲发生变形的竖直臂4-1及圆弧部分4-2在T形曲的两侧相互对称，因此，在进行T形曲矫治力模型建立时，只需对对称的一部分进行分析；T形曲圆弧部分4-2弯曲半径为R，整体高度为h，关闭间隙为b，关闭间隙由临床上后抽弓丝加力时产生；
[0149] 对矫治过程中T形曲的变形情况进行分析后可知，T形曲在竖直臂4-1及T形曲圆弧部分水平臂4-3发生变形，T形曲产生的正畸力由两变形区域产生的回复力组成，因此需对T形曲竖直臂4-1及T形曲圆弧部分水平臂4-3分别进行力学分析，并将二者叠加，建立T形曲矫治力预测模型。
[0150] 进一步地，所述的步骤2)中，T形曲竖直臂4-1的转角方程θ(x)和挠度方程v(x)表达为：
[0151]
[0152] 式中，M(x)是竖直部分上x距离处所受弯矩，E为材料的弹性模量，Iz为弓丝截面对z轴的惯性矩，对于圆丝Iz＝πd4/64，d为圆丝直径，对于矩形丝Iz＝c1c23/12，c2为矩形丝截面上与z轴平行边的长度，c1为矩形丝截面上与z轴垂直边的长度，C0和D0是积分常数，C0和D0由边界条件确定，T形曲竖直臂4-1的弯矩方程为：
[0153] M(x)＝-P(y-x) (2)
[0154] 式中，P是产生竖直臂转角θ(x)变形所需的力，y为竖直臂变形前长度；
[0155] 将式(2)带入到式(1)中进行积分，可得：
[0156]
[0157] 为确定式(3)中的积分常数C0、D0，需要确定T形曲竖直臂4-1的边界条件，基于支撑条件，其挠度或转角常为零或已知，对x＝0，即T形曲竖直臂4-1及圆弧部分水平臂4-3的连接处的变形量进行求解，由于T形曲竖直臂4-1及圆弧部分水平臂4-3的连接处有一纵向对称面，且外力均作用于这一对称面上，因此，变形后的曲梁轴线仍位于该纵向对称面内，该变形属于曲梁的平面弯曲变形问题，故可将连接处的圆弧等效为弧度为π/4的弯曲梁，对其取一段弧度为dα微元；
[0158] 在曲梁平面弯曲的情况下，外力都在曲梁的纵向对称面内，变形后的轴线仍为这一对称面内的曲线，曲梁并无扭转变形，此时仍可使用平面假设，可得到竖直臂4-1及圆弧部分水平臂4-3连接处曲梁变形后挠曲线的微分方程如式(4)，在给定的边界条件下，积分该方程便可确定曲梁的变形；
[0159]
[0160] 式中，u为竖直臂4-1及圆弧部分水平臂4-3连接处曲梁横截面在x方向上的位移，连接处所受转矩M0＝M|x＝0＝-Py，Iω为连接处曲梁横截面对ω轴的惯性矩，由于T形曲竖直臂4-1与连接处曲梁的弯曲类型一致，有Iω＝Iz。
[0161] 由弧长公式可知ds＝Rdα，因此竖直臂4-1及圆弧部分水平臂4-3连接处曲梁变形后挠曲线的微分方程可变为：
[0162]
[0163] 解连接处曲梁变形后挠曲线的常系数非齐次微分方程求得：
[0164]
[0165] 由于竖直臂4-1及圆弧部分水平臂4-3连接处曲梁沿纵向对称面对称，竖直臂4-1及圆弧部分水平臂4-3曲梁的边界条件为 解得：
[0166] 连接处曲梁的挠度方程可表示为：
[0167]
[0168] 弯梁的转角方程为：
[0169]
[0170] 因此，有边界条件 解得：
[0171] 将C0，D0带入(3)中可得：
[0172]
[0173] 因为最大转角及最大挠度均产生在抽丝端处，即x＝y处，最大挠度即为抽丝端T形曲水平臂4-4移动距离m，可知：
[0174]
[0175] 由作用力反作用力原则，T形曲竖直臂4-1形变产生的矫治力F1即为产生该变形量所需力的反力：
[0176]
[0177] 进一步地，所述的步骤3)中，T形曲圆弧部分水平臂4-3与T形曲竖直臂4-1变形情况相同，因此可以使用相同的力学建模方法，在正畸治疗过程中，正畸医师会预先拉动T形曲4抽丝端引起T形曲4形变，并将形变后的T形曲安装在病人的牙齿上，由T形曲4形变后产生的回复力拉动患者牙齿移动，达到矫治患者畸形牙齿的目的，T形曲圆弧部分水平臂4-3的弯曲挠度，即T形曲圆弧部分水平臂4-3沿Z轴方向的位移，在T形曲形变过程中，T形曲圆弧部分水平臂4-3与T形曲竖直臂4-1连接处的对称中心在空间中的位置是不断变化的，因此，为了计算T形曲圆弧部分水平臂4-3弯曲挠度s，需对T形曲4形变过程中连接处对称中心的运动情况进行简化，将变形前连接处的对称中心重合，计算变形后竖直臂理论长度与竖直臂变形前长度y的差值，即T形曲圆弧部分水平臂弯曲挠度s为：
[0178]
[0179] 可得到T形曲圆弧部分水平臂4-3挠曲线近似微分方程为：
[0180]
[0181] 式中，M(l)为T形曲圆弧部分水平臂4-3的l距离处的弯矩；
[0182] T形曲圆弧部分水平臂4-3的转角方程θ(l)和挠度方程v(l)可通过对式(13)进行积分得到：
[0183]
[0184]
[0185] 式中，C1和D1是由边界条件确定的积分常数，而T形曲圆弧部分水平臂4-3的弯矩方程可由式(16)表达：
[0186] M(l)＝-G(w-R-l) (16)
[0187] 式中，G为圆弧部分水平臂产生的矫治力，w为T形曲圆弧部分水平臂4-3的长度；
[0188] 将式(13)和(16)代入到式(14)和(15)中进行积分，可得：
[0189]
[0190]
[0191] T形曲圆弧部分水平臂4-3的边界条件求解与T形曲竖直臂4-1相同，将一侧的侧面圆弧弧度为π/4，对其取一段弧度为dβ微元，由于对竖直臂4-1边界条件的求解给出了详细的过程，因此，在这里直接给出解圆弧段挠度常系数非齐次方程得到的边界条件式；
[0192]
[0193] 由于T形曲4侧面圆弧的对称性，侧面圆弧的边界条件为 解得：
[0194] 则弯梁的挠度方程为：
[0195]
[0196] 弯梁的转角方程为：
[0197]
[0198] 因此，有边界条件v|l＝0＝u|β＝0＝0，θ|l＝0＝ε|β＝0＝0，解得：C1＝0，D1＝0；
[0199] 将C1和D1代入(17)、(18)中可得：
[0200]
[0201]
[0202] 式中，T形曲圆弧部分水平臂4-3最大转角及最大挠度产生在l＝w-R处，最大挠度即为T形曲圆弧部分水平臂4-3移动距离s，可知：
[0203]
[0204] 沿抽丝端方向的正畸力F2即为由T形曲圆弧部分水平臂4-3形变产生圆弧部分水平臂转角θ(l)变形所需力沿抽丝端方向分力的反力，即：
[0205]
[0206] 进一步地，所述的步骤4)中，T形曲4产生的静态矫治力即为竖直臂4-1及T形曲圆弧部分4-2形变后产生的回复力的合力，即：
[0207]
[0208] 进一步地，本方法应用于一种模拟牙齿移动用蜡质颌堤。
[0209] 进一步地，本方法适用的一种模拟牙齿移动用蜡质颌堤由基托蜡颌堤1、树脂牙齿模型2、正畸托槽3、T形曲4组成。
[0210] 进一步地，所述的步骤5)中，被测牙齿与测量原件以树脂圆柱体连接，牙齿在蜡制颌堤中的移动实际为圆柱连接体在蜡制颌堤中的运动，因此，以圆柱体作为基本构件进行分析；当牙齿在蜡制颌堤中移动的速度为vt时，vt为蜡制颌堤在t时刻的流动速度，圆柱体上受到沿移动方向的作用力为绕流拖曳力；摩擦拖曳力和压差拖曳力共同组成绕流拖曳力；摩擦拖曳力是由于流体的粘滞性在柱体表面形成边界层，在此边界层范围内，流体产生速度梯度，摩擦效应显著，产生了摩擦切应力；压差拖曳力是由于边界层在圆柱体表面某点处分离，在分离点下游即在圆柱体后部形成很强的漩涡尾流，使得圆柱前后产生压力差，进而在流动方向产生了一个力，而在流体流动中，圆柱体旋涡尾流是随雷诺数的Re的变化发展的，牙齿在蜡制颌堤中的移动雷诺数Re＜5，因此无旋涡尾流的产生，无压差拖曳力产生；
[0211] 对单位长度圆柱体上的拖曳力fD可用式(27)计算：
[0212]
[0213] 式中，v0为未受绕流影响垂直于圆柱体轴线的牙齿移动速度分量，ρ(t)为在t时刻实验温度下蜡制颌堤的密度，A为单位长度圆柱体垂直于移动方向的投影面积，对于圆柱体，A＝1×D，D为圆柱体的直径，CD为拖曳力系数，它集中反映了流体的粘滞性而引起的粘滞效应，与雷诺数Re和圆柱面粗糙度Ra有关；
[0214] 假设本研究中的蜡制颌堤流体是不可压缩的理想流体，排蜡体积为 的圆柱体在移动速度vt＝v(x,y,z,t)的蜡制颌堤流场中移动；暂不考虑圆柱体对蜡制颌堤流场的影响，即假定蜡制颌堤流场内的压强分布不因圆柱体的存在而改变，那么圆柱体的边界作为加速流体边界的一部分，也就是被圆柱体置换的那部分体积内的蜡制颌堤流体，它本应以静止的状态存在于蜡制颌堤流场中，但实际上由于圆柱体移动的存在，这部分静止的蜡制颌堤流体将被加速至与圆柱体边界移速相同的状态；因此加速的蜡制颌堤流体将会对排蜡体积为 的圆柱体沿流动方向作用一个惯性力Fk，惯性力Fk的数值等于圆柱体的排蜡质量M0与体积 内蜡制颌堤流体的平均加速度 的乘积，即：
[0215]
[0216] 对于研究中的圆柱体来说， 可以取圆柱体轴中心位置处的流体加速度 来表示，此时：
[0217]
[0218] 但由于圆柱体存在于蜡制颌堤流场中，必将使圆柱体周围的流体质点受到扰动引起速度变化，从而改变蜡制颌堤流场内的压强分布，所以，圆柱体的扰动是圆柱体周围改变了原来运动状态的那部分附加流体的质量Mw沿流体流动方向也将对圆柱体产生一个附加惯性力，即附加质量力；因此加速的流体沿流动方向实际作用在圆柱体上的绕流惯性力fL可表示为：
[0219]
[0220] 令Mw＝CmM0，则式(30)可表示为：
[0221]
[0222] 式中，Cm为附加质量系数，CM为质量系数，亦称为惯性力系数，集中反映了由于流体的惯性以及圆柱体的存在，使圆柱体周围蜡制颌堤流场的速度改变而引起的附加质量效应；
[0223] 经过以上分析，可获得牙齿在蜡制颌堤移动过程中所受阻力情况，牙齿在正畸弓丝变形产生的正畸力影响下在蜡制颌堤中移动，移动过程中，由于流场的绕流特性，牙齿受到绕流惯性力fL和拖曳力fD的影响；
[0224] 受热量交换影响，热场中蜡制颌堤模型内部温度是随时间变化的，内部温度的变化将引起蜡制颌堤模型密度的变化，进而影响牙齿在蜡制颌堤内部移动受到的阻力。牙齿模型在蜡制颌堤内部移动时遵循粘性流体能量方程中的规律，令e代表单位质量流体所具有的内能，则ρ(t)e为单位体积流体具有的内能，ρ(t)vt2/2代表单位体积的动能，从而单位2
体积流体包含的总能量E＝ρ(t)e+ρ(t)vt/2；
[0225] 经过简化整理，能量守恒原理可近似地表示为：
[0226]
[0227] 式中，cp为无量纲压强系数，Φ为牙齿模型在蜡制颌堤流体中移动时消耗的机械功，k为计算系数，▽Τ为基托蜡流体热场的温度梯度，q为热流密度；
[0228] 对▽Τ进行求解，设蜡制颌堤的厚度为2δ，初始温度为t0；在初始瞬间将它放置于温度为t∞的流体中，流体与蜡制颌堤间的表面传热系数h为常数，蜡制颌堤两边对称受热，蜡制颌堤内部温度分布必以其中心截面为对称面，因此，只需要研究厚度为δ的半块蜡制颌堤，把x轴的原点置于蜡制颌堤的中心截面上，对于x≥0的半块蜡制颌堤，可列出如下导热微分方程：
[0229]
[0230] 式中，a为热扩散率，式(33)两边对x积分，可得：
[0231]
[0232] 蜡制颌堤在水浴环境下均匀受热可以简化为一维热场问题，因此有:
[0233]
[0234] 将蜡制颌堤流体热场的温度梯度 代入式(32)中整理得到：
[0235]
[0236] 对式(36)等式两边对t积分整理得到：
[0237]
[0238] 式中，T为蜡制颌堤流体热场的温度；
[0239] 傅里叶定律用热流密度q表示时有如下形式：
[0240]
[0241] 式中，λ为导热系数；
[0242] 将式(38)代入式(37)中可得蜡质颌堤密度ρ随时间t变化的表达式：
[0243]
[0244] 蜡制颌堤模拟牙齿移动过程中动态阻力模型可由式(40)表达：
[0245]
[0246] 式中，f为蜡制颌堤模拟牙齿移动过程中的动态阻力。
[0247] 进一步地，所述的步骤6)中，基于T形曲静态矫治力预测模型及蜡制颌堤模拟牙齿移动过程中动态阻力模型，T形曲动态矫治力预测模型可表达为：
[0248]
[0249] 式中，F即为应用T形曲动态矫治力预测模型得出的T形曲动态矫治力。
[0250] 进一步地，所述的步骤7)中，物理学中基于加速度通过积分可以得到位移和速度矢量关系：
[0251]
[0252] 牙齿移动的加速度a0、牙齿的质量m0、作用力间的关系可由式(43)表达：
[0253] F＝m0a0(43)
[0254] 则受到正畸力F0、绕流惯性力fL和拖曳力fD作用在蜡质颌堤中移动的牙齿的加速度a0由式(44)表达：
[0255]
[0256] 将式(44)代入式(42)中可得到基于蜡制颌堤的T形曲作用下牙齿移动预测模型：
[0257]
[0258] 式中，Sm为基于蜡制颌堤的T形曲作用下牙齿移动，m0为牙齿的质量，t1为牙齿移动的时间。
[0259] 进一步地，在应用所述一种模拟牙齿移动用蜡制颌堤预测T形曲矫治效果时，首先将正畸托槽粘贴在树脂牙齿模型的外表面上，将T形曲固定在正畸托槽上，此时，将蜡制颌堤与T形曲同时浸入75℃恒温水浴的环境中，2min后取出，观察水浴前后牙齿位置即可清楚地了解T形曲作用下牙齿移动的情况，应用本方法提出的一种T形曲动态矫治力预测模型可以通过T形曲的形状参数计算产生动态矫治力的大小，依据牙齿移动的情况，通过本模型计算并调整T形曲的形状参数获得最佳的治疗效果，进而辅助正畸医生制定更合理的矫治方案。
[0260] 以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是说明本发明的原理，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
