CN111588497A - 一种基于正畸弓丝弯制点角距比和的平面等半径圆域划分半径确定方法 - Google Patents

一种基于正畸弓丝弯制点角距比和的平面等半径圆域划分半径确定方法 Download PDF

Info

Publication number
CN111588497A
CN111588497A CN202010484782.6A CN202010484782A CN111588497A CN 111588497 A CN111588497 A CN 111588497A CN 202010484782 A CN202010484782 A CN 202010484782A CN 111588497 A CN111588497 A CN 111588497A
Authority
CN
China
Prior art keywords
bending
radius
equal
circular
circle
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202010484782.6A
Other languages
English (en)
Other versions
CN111588497B (zh
Inventor
姜金刚
吴殿昊
姚亮
郭亚峰
张永德
曾阳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin University of Science and Technology
Original Assignee
Harbin University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Harbin University of Science and Technology filed Critical Harbin University of Science and Technology
Priority to CN202010484782.6A priority Critical patent/CN111588497B/zh
Publication of CN111588497A publication Critical patent/CN111588497A/zh
Application granted granted Critical
Publication of CN111588497B publication Critical patent/CN111588497B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61CDENTISTRY; APPARATUS OR METHODS FOR ORAL OR DENTAL HYGIENE
    • A61C7/00Orthodontics, i.e. obtaining or maintaining the desired position of teeth, e.g. by straightening, evening, regulating, separating, or by correcting malocclusions
    • A61C7/12Brackets; Arch wires; Combinations thereof; Accessories therefor
    • A61C7/20Arch wires
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61CDENTISTRY; APPARATUS OR METHODS FOR ORAL OR DENTAL HYGIENE
    • A61C7/00Orthodontics, i.e. obtaining or maintaining the desired position of teeth, e.g. by straightening, evening, regulating, separating, or by correcting malocclusions
    • A61C7/002Orthodontic computer assisted systems
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61CDENTISTRY; APPARATUS OR METHODS FOR ORAL OR DENTAL HYGIENE
    • A61C7/00Orthodontics, i.e. obtaining or maintaining the desired position of teeth, e.g. by straightening, evening, regulating, separating, or by correcting malocclusions
    • A61C7/02Tools for manipulating or working with an orthodontic appliance
    • A61C7/026Tools for manipulating or working with an orthodontic appliance for twisting orthodontic ligature wires

Landscapes

  • Health & Medical Sciences (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Dentistry (AREA)
  • Epidemiology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Engineering & Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Dental Tools And Instruments Or Auxiliary Dental Instruments (AREA)

Abstract

本发明公开了一种基于正畸弓丝弯制点角距比和的平面等半径圆域划分半径确定方法,它涉及正畸弓丝弯制技术领域,本发明针对弯制点密集程度较小的一类正畸弓丝曲线,基于正畸弓丝曲线弯制点信息集、弯制点的机器人弯制信息集,结合机器人弯制正畸弓丝的运动特点,建立一种基于正畸弓丝弯制点角距比和的平面等半径圆域划分半径确定方法。技术要点为:导入圆域划分数据并将正畸弓丝曲线T转换为平面曲线T′;根据弯制点角距比信息确定等半径圆域划分个数初始值n;划分等半径确定圆域;寻找最佳划分个数;输出合理角距比等半径圆域划分半径requal。本发明采用角距比和判定调整试划分圆域个数,确定合理角距比等半径圆域的划分半径值,提高等半径圆域划分的效率,进而提高正畸弓丝弯制规划的效率,避免了机器人弯制正畸弓丝过程中出现干涉及弯制复杂的问题。

Description

一种基于正畸弓丝弯制点角距比和的平面等半径圆域划分半 径确定方法
技术领域
本发明专利涉及一种基于正畸弓丝弯制点角距比和的平面等半径圆域划分半径确定方法,属于正畸弓丝弯制技术领域。
背景技术
错颌畸形是危及人体健康的第三大口腔疾病,呈现较高的发病率,现代口腔医学中,固定矫治是一种常用且有效的正畸治疗手段,而正畸弓丝的弯制是固定矫治技术的关键,近年来,深受数字化制造技术的影响,传统的口腔制造加工工艺正发生革命性变化,口腔正畸领域也受益于数字化技术,正畸矫治器中弓丝的加工正在向数字化发展。
在机器人弯制个性化正畸弓丝的过程中,个性化正畸弓丝与机器人弯制手爪间可能发生干涉,干涉即个性化正畸弓丝与机器人弯制手爪发生碰撞,干涉发生后会大大影响个性化正畸弓丝的弯制精度,进而影响矫治效果,致使所弯制的个性化正畸弓丝无法应用于临床治疗;研究发现,在正向弯制个性化正畸弓丝的过程中,正向弯制即由未弯制的正畸弓丝弯制成复杂的成形弓丝,干涉往往是由于不合理的成形控制点弯制顺序引起的,而合理的正向弯制成形控制点弯制顺序可以有效地避免干涉的发生,获得合理的正向弯制成形控制点弯制顺序是实现正畸弓丝数字化弯制的必要前提,但目前正畸弓丝弯制技术领域缺少规划正向弯制成形控制点弯制顺序的方法,难以实现正畸弓丝数字化弯制。
对于正畸弓丝弯制规划的划分方法的研究,在本发明人已经授权的授权公告号为CN107647925B的发明专利《一种用于正畸弓丝弯制规划的圆域划分方法》中提出了一种等半径圆域划分方法,在正畸弓丝曲线上进行区域划分,最终对每个区域进行排序,从而得出最终弯制点的弯制顺序,虽然这种方法对正畸弓丝弯制规划具有一定应用价值,但该方法中提及的等半径圆域划分方法仅是部分弧段对平面正畸弓丝曲线进行区域切割,无法形成完整的规则区域,因此难以在规定区域内定义合理的划分依据,仅以一种无依据的均匀化标准对正畸弓丝曲线进行划分,所划分的圆域区间内通常会存在弯制点密集程度过大、弯制复杂程度过大的情况,即划分后的区间未充分考虑到正畸弓丝曲线上弯制点分布信息的个性化特点,比如患者个性化正畸弓丝上的弯制点常常存在整体弯制点密集程度相对较小的情况,此时通常可采用一种衡量正畸弓丝曲线上单个弯制点与相邻弯制点间紧密程度的圆单位圆域弯制点密度
Figure BDA0002518648960000012
来验证,表达式为
Figure BDA0002518648960000011
其中公式中的数值1表示单位圆域中的一个弯制点,lj表示弯制点tj与其距离最近的弯制点之间的直线距离,j的取值范围为1≤j≤i,i表示正畸弓丝曲线上的弯制点个数,若通过计算可知所有单位圆域弯制点密度中的最大值
Figure BDA0002518648960000013
不大于所规定的单位圆域弯制点密度的上限值
Figure BDA0002518648960000014
可知在所取的个性化正畸弓丝曲线上,无论划分后的圆域上包含符合弯制点个数限制的任意数量的弯制点,划分后的圆域依然能保证圆域弯制点的密集程度在可接受的范围内,能够满足弯制系统对弯制点密度的要求,因此在满足上述要求的一类正畸弓丝曲线的圆域划分的过程中无需考虑弯制点密度因素的影响,即该弓丝弯制点分布具有特殊属性,在对此类个性化正畸弓丝进行圆域划分时,现有的正畸弓丝成形控制点弯制顺序规划方法合理性差,无法对此类个性化正畸弓丝实现高效数字化弯制,从而不能有效避免由于弯制顺序规划不合理导致的弯制机器人产生空程无效动作、弯制过程相互干涉动作和弯制运动复杂动作,不利于发挥弯制机器人优势的最大化,也不能明显提高弯制效率。
发明内容
针对上述问题,本发明提出一种基于正畸弓丝弯制点角距比和的平面等半径圆域划分半径确定方法,解决目前正畸弓丝弯制技术领域缺少用于规划一类整体单位弯制点密度相对较小的正畸弓丝弯制顺序的等半径圆域的划分半径确定方法,在等半径圆域划分半径确定的过程中提出合理的圆域限制参数,仅对等半径圆域所划分的弯制点弯制复杂程度定量约束,获得一系列区域内弯制点弯制复杂程度符合要求的合理角距比等半径圆域,最终得到适用于规划此类正畸弓丝曲线弯制点的等半径圆域的划分半径,为具有特殊属性的一类正畸弓丝弯制规划时所用的等半径圆域划分方法提供便利,从而提高了正畸弓丝弯制规划的效率,发挥弯制机器人优势的最大化,保证正畸弓丝弯制过程正常运作,避免了机器人弯制此类正畸弓丝过程中出现干涉的问题。
上述目的主要通过以下方案实现:一种基于正畸弓丝弯制点角距比和的平面等半径圆域划分半径确定方法,所述方法的具体实现过程为:
步骤一、等半径确定圆域划分数据导入及正畸弓丝曲线转换:
根据患者有i个弯制点的正畸弓丝曲线,计算并输入正畸弓丝曲线弯制点信息集T={t1,t2,t3,...,ti},ti=(xi,yi,zi)′为每个正畸弓丝曲线弯制点的坐标,在每个弯制点ti上机器人执行不同的弯制运动,每一个正畸弓丝曲线弯制点ti均对应一个弯制点机器人弯制信息单元ri,输入弯制点的机器人弯制信息集为R={r1,r2,r3,...,ri},ri=(xi,yi,zi,αi)′表示机器人在弯制该点时的弯制点坐标及弯制角度,αi为机器人作用在弯制点ti上的弯制角度;
Figure BDA0002518648960000021
Figure BDA0002518648960000022
将个性化正畸弓丝曲线成形控制点信息集中各弯制点的坐标ti=(xi,yi,zi)′中的zi赋值为0,即令zi=0,获得正畸弓丝曲线转换平面正畸弓丝曲线T′;
Figure BDA0002518648960000023
步骤二、计算等半径确定圆域初始试划分个数:
按照
Figure BDA0002518648960000024
预先计算正畸弓丝曲线上i个弯制点的弯制点角距比,其中第j个弯制点的弯制点角距比规定
Figure BDA0002518648960000025
Ej是对第j个弯制点弯制复杂程度的量化描述,αj为作用在弯制点tj处的弯制角度,
Figure BDA0002518648960000026
表示作用在弯制点tj处弯制距离,即弯制点tj-1与tj之间曲线段的长度,特别地,由于第一个弯制点t1无需弯制,规定弯制点t1的弯制点角距比E1=0;首次将试划分n个等半径确定圆域,n的初始值为
Figure BDA0002518648960000031
其中
Figure BDA0002518648960000032
表示对式∑E/(∑E)max所计算结果的取整,(∑E)max表示在正畸弓丝曲线上将要划分的任意一个等半径确定圆域an内的圆域弯制点角距比和
Figure BDA0002518648960000033
所要求的上限值,
Figure BDA0002518648960000034
表示正畸弓丝曲线上第n个等半径确定圆域an内的所有弯制点的弯制点角距比之和,即对等半径确定圆域an所划分的弯制点整体弯制复杂程度的量化描述,当变半径划分圆域an内的弯制点分别为tq+1、tq+2、tq+3、…、
Figure BDA0002518648960000035
时,规定
Figure BDA0002518648960000036
圆域弯制点个数
Figure BDA0002518648960000037
为半径为
Figure BDA0002518648960000038
的等半径确定圆域an内的弯制点个数,q表示正畸弓丝曲线上在等半径确定圆域an之前生成的n-1个圆域内所有弯制点的个数,即
Figure BDA0002518648960000039
以上依次提及的弯制点角距比E、圆域弯制点角距比和
Figure BDA00025186489600000310
圆域弯制点个数
Figure BDA00025186489600000311
三个参数统称为圆域限制参数,跳转至步骤三;
步骤三、试划分等半径确定圆域:
在以第一个弯制点t1为起点、以最后一个弯制点ti为终点的正畸弓丝曲线段上选取n+1个点作为圆域形成点,首个圆域形成点为弯制点t1所在的点,末尾圆域形成点为弯制点ti所在的点,使得每个圆域形成点与相邻圆域形成点间连接所得的n条直线段长度都相等,规定水平向右的向量按顺时针方向扫过的直线段依次用
Figure BDA00025186489600000312
表示,且存在
Figure BDA00025186489600000313
试划分等半径确定圆域以弯制点t1所在的圆域形成点为起始点进行,依次以
Figure BDA00025186489600000314
的中点作为圆心,以
Figure BDA00025186489600000315
作为半径生成n个等半径确定圆域,每个等半径确定圆域的边界线通过两个圆域形成点,相邻的两个等半径确定圆域的边界线相交于一个共用的圆域形成点,即第n-1个等半径确定圆域an-1右端的圆域形成点恰好是第n个等半径确定圆域an左端的圆域形成点,规定圆域边界线所截的正畸弓丝曲线段所包含的弯制点被该圆域所划分,当两个等半径确定圆域的边界线共用的圆域形成点所在的点恰好为正畸弓丝曲线上的一个弯制点,则规定该交点所在的弯制点被前一个等半径确定圆域所划分,假如第n-1个等半径确定圆域an-1与第n个等半径确定圆域an共用的圆域形成点为弯制点tj,则弯制点tj被等半径确定圆域an-1所划分,跳转至步骤四;
步骤四、寻找最佳试划分个数:
分别计算步骤三中生成的n个等半径确定圆域的圆域弯制点个数
Figure BDA00025186489600000316
可得到圆域弯制点个数集
Figure BDA00025186489600000317
对圆域弯制点个数集Q中的n个圆域弯制点个数进行降序排列,取出其中最大的圆域弯制点个数,记为Qam,根据所要求的圆域弯制点个数上限值Qmax,Qmax表示在正畸弓丝曲线上将要划分的任意一个等半径确定圆域an内的圆域弯制点个数
Figure BDA00025186489600000318
所要求的上限值,特别地,Qmax=5,判断是否存在Qam≤5,
具体为:
如果Qam≤5不成立,说明所生成的n个等半径确定圆域中存在不符合圆域弯制点个数上限值Qmax要求的圆域,其中Qmax=5,可知此时n值不是最佳试划分个数,则需要通过改变圆域的个数,从而改变圆域半径的大小,进行重新试划分等半径确定圆域,令n=n+1,即下一次试划分等半径确定圆域时在此次划分个数的基础上增加一个,进而跳转至步骤三;
如果Qam≤5成立,说明所生成的n个等半径确定圆域都符合圆域弯制点个数上限值Qmax的要求,进而按照
Figure BDA0002518648960000041
计算步骤三中生成的n个等半径确定圆域的圆域弯制点角距比和
Figure BDA0002518648960000042
可得到圆域弯制点角距比和集
Figure BDA0002518648960000043
对圆域弯制点角距比和集M中的n个圆域弯制点角距比和进行降序排列,取出其中最大的圆域弯制点角距比和,记为(∑E)am,根据所要求的圆域弯制点角距比和上限值(∑E)max,在Qam≤5成立的情况下判断是否存在(∑E)am≤(∑E)max
具体为:
如果(∑E)am≤(∑E)max成立,说明步骤三中所生成的n个等半径确定圆域都符合圆域弯制点角距比和上限值(∑E)max的要求,即此时所有的等半径确定圆域均符合划分要求,可知此时n值恰好为最佳试划分个数,称n个圆域a1、a2、…、an都为合理角距比等半径圆域,跳转至步骤五;
如果(∑E)am≤(∑E)max不成立,则说明步骤三中所生成的n个等半径确定圆域存在不符合圆域弯制点角距比和上限值(∑E)max要求的圆域,可知此时n值不是最佳试划分个数,则需要通过改变圆域的个数,从而改变圆域半径的大小,进行重新试划分等半径确定圆域,令n=n+1,即下一次试划分等半径确定圆域时在此次划分个数的基础上增加一个,进而跳转至步骤三;
步骤五、输出合理角距比等半径圆域划分半径
得到步骤四中输出的n个合理角距比等半径圆域及每个合理角距比等半径圆域的划分半径,划分半径值依次为
Figure BDA0002518648960000044
Figure BDA0002518648960000045
则requal即为可以在正畸弓丝曲线上划分出n个合理角距比等半径圆域的通用划分半径,输出合理角距比等半径圆域划分半径requal,程序结束。
本发明的有益效果为:
1、本发明针对等半径圆域的划分半径的确定,采用了三个圆域限制参数作为等半径圆域划分半径的计算依据,提及圆域弯制点个数
Figure BDA0002518648960000046
弯制点角距比E和圆域弯制点角距比和
Figure BDA0002518648960000047
的概念,对单个弯制点弯制复杂程度、圆域内弯制点总体弯制复杂程度的量化描述,首先基于弯制点角距比均分原则,确定圆域试划分个数n的初始值,以n的初始值为起始值进行划分,而不是以无依据的n的初始值为起始进行等半径圆域试划分,有效提高了寻找最佳划分个数n的速度;在生成等半径确定圆域后首先利用Qmax进行弯制点个数约束,再利用圆域弯制点角距比和上限值(∑E)max进行条件约束,即不满足Qmax的条件限制时就可立即反馈,可充分提升算法的计算效率,通过二者的上限限制,可在一条正畸弓丝曲线上形成若干个符合设定要求的合理角距比等半径圆域,进而得到符合要求的等半径圆域的划分半径requal,采用本方法确定的划分半径对具有特殊属性一类正畸弓丝进行圆域划分时,可有效避免所划分的圆域出现弯制复杂程度过高的现象,从而最大程度避免了机器人在弯制过于复杂的问题,圆域限制参数作为一种规划指标应用于正畸弓丝弯制规划领域,为一类单位弯制点密度较小的正畸弓丝的圆域划分半径确定提供了理论依据。
2、同本发明人同日申报的发明专利《一种正畸弓丝弯制规划等半径圆域划分半径确定方法》相比,本方法以一类个性化正畸弓丝曲线上弯制点具有单位弯制点密度较小的特殊属性为前提,由于在确定等半径圆域划分半径前对单位圆域弯制点密度
Figure BDA0002518648960000051
进行整体验证,因此在试划分圆域过程中省略了对圆域的弯制点密度
Figure BDA0002518648960000052
的判定,从而圆域试划分过程完全基于圆域弯制点角距比和
Figure BDA0002518648960000053
进行,最终确定等半径圆域的划分半径requal,不仅满足了所确定的划分半径的合理性,而且简化了划分半径确定过程,使划分半径确定算法的复杂程度降低,提高了寻找合理划分半径的效率。
3、同本发明人同日申报的发明专利《一种基于正畸弓丝弯制点密度的平面等半径圆域划分半径确定方法》相比,虽然两方法都适用于具有特殊属性的一类个性化正畸弓丝曲线,但对于《一种基于正畸弓丝弯制点密度的平面等半径圆域划分半径确定方法》中提及的方法所侧重的前提是单个弯制点角距比都满足设定要求,进而仅以圆域弯制点个数
Figure BDA0002518648960000054
和圆域弯制点密度
Figure BDA0002518648960000055
作为确定划分半径的依据,本方法所侧重的前提是单位弯制点密度满足设定要求,进而仅以圆域弯制点个数
Figure BDA0002518648960000056
和圆域弯制点角距比和
Figure BDA0002518648960000057
作为确定划分半径的依据,两种方法在进行等半径圆域划分半径的确定时应用于不同类别的正畸弓丝,因此本方法的提出与另种方法互相补偿,进而完善了正畸弓丝弯制点等半径圆域划分半径确定的系列方法。
4、同本发明人已经授权的授权公告号为CN107647925B的发明专利《一种用于正畸弓丝弯制规划的圆域划分方法》相比,发明专利《一种用于正畸弓丝弯制规划的圆域划分方法》属于一种等半径圆域划分的方法,存在所划分的圆域区间内弯制点密度过大、弯制难度过大的情况,即所生成的圆域区间未充分考虑到正畸弓丝曲线上弯制点分布信息的个性化特点,仅以一种无依据的均匀化标准对正畸弓丝曲线进行划分,且所提出的圆域划分过程仅仅是以圆弧对正畸弓丝曲线进行分割进而得到区域,不是严格意义上的圆域划分,而本专利提供的方法属于一种等半径圆域划分半径的确定方法,针对一类单位圆域弯制点密度相对较小的正畸弓丝而提出,在确定过程中同样涉及到等半径圆域的划分,可形成规则完整的圆域,并基于规则的区域对所需的圆域限制参数进行定义,依据所提出的圆域弯制点个数及圆域弯制点角距比和对所划分等半径圆域的弯制点进行弯制复杂程度的定量约束,从而导致等半径圆域的划分半径为符合圆域限制参数的规定而发生变化,最终得到符合一类单位圆域弯制点密度较小的正畸弓丝曲线个性化特点的等半径圆域的合理划分半径,利用所确定的划分半径进行圆域划分,可有效避免完成划分后的各圆域内弯制点的弯制复杂程度相差悬殊的情况,提高了各区域的均匀性,弯制机器人将不会产生空程无效动作或弯制运动复杂动作动作,从而可以发挥弯制机器人优势的最大化,保证正畸弓丝弯制过程正常运作,提高正畸弓丝弯制规划的效率,避免了机器人弯制特殊属性正畸弓丝过程中出现干涉的问题。
附图说明
为了易于说明,本发明由下述的具体实施及附图作以详细描述。
图1为一种基于正畸弓丝弯制点角距比和的平面等半径圆域划分半径确定方法流程图;
图2为空间个性化正畸弓丝弯制点分布示意图;
图3为等半径圆域初始试划分平面个性化正畸弓丝曲线示意图;
图4为等半径圆域完成试划分平面个性化正畸弓丝曲线示意图;
具体实施方式
为使本发明专利的目的、技术方案和优点更加清楚明了,下面通过附图中示出的具体实施例来描述本发明专利,但是应该理解,这些描述只是示例性的,而并非要限制本发明专利的范围,此外,在以下说明中,省略了对公知结构和技术的描述,以避免不必要地混淆本发明专利的概念。
实施例1:如图1、图2、图3、图4所示,本具体实施方式采用以下技术方案:一种基于正畸弓丝弯制点角距比和的平面等半径圆域划分半径确定方法,所述方法的具体实现过程为:
步骤一、等半径确定圆域划分数据导入及正畸弓丝曲线转换:
根据患者有i个弯制点的正畸弓丝曲线,计算并输入正畸弓丝曲线弯制点信息集T={t1,t2,t3,...,ti},ti=(xi,yi,zi)′为每个正畸弓丝曲线弯制点的坐标,在每个弯制点ti上机器人执行不同的弯制运动,每一个正畸弓丝曲线弯制点ti均对应一个弯制点机器人弯制信息单元ri,输入弯制点的机器人弯制信息集为R={r1,r2,r3,...,ri},ri=(xi,yi,zi,αi)′表示机器人在弯制该点时的弯制点坐标及弯制角度,αi为机器人作用在弯制点ti上的弯制角度;
Figure BDA0002518648960000061
Figure BDA0002518648960000062
将个性化正畸弓丝曲线成形控制点信息集中各弯制点的坐标ti=(xi,yi,zi)′中的zi赋值为0,即令zi=0,获得正畸弓丝曲线转换平面正畸弓丝曲线T′;
Figure BDA0002518648960000063
步骤二、计算等半径确定圆域初始试划分个数:
按照
Figure BDA0002518648960000064
预先计算正畸弓丝曲线上i个弯制点的弯制点角距比,其中第j个弯制点的弯制点角距比规定
Figure BDA0002518648960000065
Ej是对第j个弯制点弯制复杂程度的量化描述,αj为作用在弯制点tj处的弯制角度,
Figure BDA0002518648960000066
表示作用在弯制点tj处弯制距离,即弯制点tj-1与tj之间曲线段的长度,特别地,由于第一个弯制点t1无需弯制,规定弯制点t1的弯制点角距比E1=0;首次将试划分n个等半径确定圆域,n的初始值为
Figure BDA0002518648960000071
其中
Figure BDA0002518648960000072
表示对式∑E/(∑E)max所计算结果的取整,(∑E)max表示在正畸弓丝曲线上将要划分的任意一个等半径确定圆域an内的圆域弯制点角距比和
Figure BDA0002518648960000073
所要求的上限值,
Figure BDA0002518648960000074
表示正畸弓丝曲线上第n个等半径确定圆域an内的所有弯制点的弯制点角距比之和,即对等半径确定圆域an所划分的弯制点整体弯制复杂程度的量化描述,当变半径划分圆域an内的弯制点分别为tq+1、tq+2、tq+3、…、
Figure BDA0002518648960000075
时,规定
Figure BDA0002518648960000076
圆域弯制点个数
Figure BDA0002518648960000077
为半径为
Figure BDA0002518648960000078
的等半径确定圆域an内的弯制点个数,q表示正畸弓丝曲线上在等半径确定圆域an之前生成的n-1个圆域内所有弯制点的个数,即
Figure BDA0002518648960000079
以上依次提及的弯制点角距比E、圆域弯制点角距比和
Figure BDA00025186489600000710
圆域弯制点个数
Figure BDA00025186489600000711
三个参数统称为圆域限制参数,跳转至步骤三;
步骤三、试划分等半径确定圆域:
在以第一个弯制点t1为起点、以最后一个弯制点ti为终点的正畸弓丝曲线段上选取n+1个点作为圆域形成点,首个圆域形成点为弯制点t1所在的点,末尾圆域形成点为弯制点ti所在的点,使得每个圆域形成点与相邻圆域形成点间连接所得的n条直线段长度都相等,规定水平向右的向量按顺时针方向扫过的直线段依次用
Figure BDA00025186489600000712
表示,且存在
Figure BDA00025186489600000713
试划分等半径确定圆域以弯制点t1所在的圆域形成点为起始点进行,依次以
Figure BDA00025186489600000714
的中点作为圆心,以
Figure BDA00025186489600000715
作为半径生成n个等半径确定圆域,每个等半径确定圆域的边界线通过两个圆域形成点,相邻的两个等半径确定圆域的边界线相交于一个共用的圆域形成点,即第n-1个等半径确定圆域an-1右端的圆域形成点恰好是第n个等半径确定圆域an左端的圆域形成点,规定圆域边界线所截的正畸弓丝曲线段所包含的弯制点被该圆域所划分,当两个等半径确定圆域的边界线共用的圆域形成点所在的点恰好为正畸弓丝曲线上的一个弯制点,则规定该交点所在的弯制点被前一个等半径确定圆域所划分,假如第n-1个等半径确定圆域an-1与第n个等半径确定圆域an共用的圆域形成点为弯制点tj,则弯制点tj被等半径确定圆域an-1所划分,跳转至步骤四;
步骤四、寻找最佳试划分个数:
分别计算步骤三中生成的n个等半径确定圆域的圆域弯制点个数
Figure BDA00025186489600000716
可得到圆域弯制点个数集
Figure BDA00025186489600000717
对圆域弯制点个数集Q中的n个圆域弯制点个数进行降序排列,取出其中最大的圆域弯制点个数,记为Qam,根据所要求的圆域弯制点个数上限值Qmax,Qmax表示在正畸弓丝曲线上将要划分的任意一个等半径确定圆域an内的圆域弯制点个数
Figure BDA0002518648960000081
所要求的上限值,特别地,Qmax=5,判断是否存在Qam≤5,
具体为:
如果Qam≤5不成立,说明所生成的n个等半径确定圆域中存在不符合圆域弯制点个数上限值Qmax要求的圆域,其中Qmax=5,可知此时n值不是最佳试划分个数,则需要通过改变圆域的个数,从而改变圆域半径的大小,进行重新试划分等半径确定圆域,令n=n+1,即下一次试划分等半径确定圆域时在此次划分个数的基础上增加一个,进而跳转至步骤三;
如果Qam≤5成立,说明所生成的n个等半径确定圆域都符合圆域弯制点个数上限值Qmax的要求,进而按照
Figure BDA0002518648960000082
计算步骤三中生成的n个等半径确定圆域的圆域弯制点角距比和
Figure BDA0002518648960000083
可得到圆域弯制点角距比和集
Figure BDA0002518648960000084
对圆域弯制点角距比和集M中的n个圆域弯制点角距比和进行降序排列,取出其中最大的圆域弯制点角距比和,记为(∑E)am,根据所要求的圆域弯制点角距比和上限值(∑E)max,在Qam≤5成立的情况下判断是否存在(∑E)am≤(∑E)max
具体为:
如果(∑E)am≤(∑E)max成立,说明步骤三中所生成的n个等半径确定圆域都符合圆域弯制点角距比和上限值(∑E)max的要求,即此时所有的等半径确定圆域均符合划分要求,可知此时n值恰好为最佳试划分个数,称n个圆域a1、a2、…、an都为合理角距比等半径圆域,跳转至步骤五;
如果(∑E)am≤(∑E)max不成立,则说明步骤三中所生成的n个等半径确定圆域存在不符合圆域弯制点角距比和上限值(∑E)max要求的圆域,可知此时n值不是最佳试划分个数,则需要通过改变圆域的个数,从而改变圆域半径的大小,进行重新试划分等半径确定圆域,令n=n+1,即下一次试划分等半径确定圆域时在此次划分个数的基础上增加一个,进而跳转至步骤三;
步骤五、输出角距比合理等半径圆域划分半径
得到步骤四中输出的n个合理角距比等半径圆域及每个合理角距比等半径圆域的划分半径,划分半径值依次为
Figure BDA0002518648960000085
Figure BDA0002518648960000086
则requal即为可以在正畸弓丝曲线上划分出n个合理角距比等半径圆域的通用划分半径,输出合理角距比等半径圆域划分半径requal,程序结束。
实施例2:如图2所示,存在某一条包含i=17个的弯制点的个性化正畸弓丝曲线,所有单位圆域弯制点密度中的最大值
Figure BDA0002518648960000087
都不大于所规定的单位圆域弯制点密度的上限值
Figure BDA0002518648960000088
对该条个性化正畸弓丝曲线进行基于弯制点角距比和的平面等半径圆域划分半径确定的过程中,如图3所示,假设经过步骤二的计算可知,首次划分n=5个等半径确定圆域,继续进行步骤三,定义首次试划分的5个等半径圆域的圆心及半径,并在正畸弓丝曲线上生成半径相同的5个圆域对弯制点进行划分,进行步骤四验证所得,n=5时等半径确定圆域不符合所设定的圆域限制参数的要求,则继续增加等半径圆域的试划分个数n,如图4所示,直至等半径圆域的试划分个数n=9时,此时进行步骤四验证所得等半径圆域符合所设定的圆域限制参数的要求,将此次划分的9个等半径圆域定义为合理角距比等半径圆域,进而跳转至步骤五,最终输出合理角距比等半径圆域划分半径requal,程序结束。

Claims (1)

1.一种基于正畸弓丝弯制点角距比和的平面等半径圆域划分半径确定方法,其特征在于:所述方法的具体实现过程为:
步骤一、等半径确定圆域划分数据导入及正畸弓丝曲线转换:
根据患者有i个弯制点的正畸弓丝曲线,计算并输入正畸弓丝曲线弯制点信息集T={t1,t2,t3,...,ti},ti=(xi,yi,zi)′为每个正畸弓丝曲线弯制点的坐标,在每个弯制点ti上机器人执行不同的弯制运动,每一个正畸弓丝曲线弯制点ti均对应一个弯制点机器人弯制信息单元ri,输入弯制点的机器人弯制信息集为R={r1,r2,r3,...,ri},ri=(xi,yi,zi,αi)′表示机器人在弯制该点时的弯制点坐标及弯制角度,αi为机器人作用在弯制点ti上的弯制角度;
Figure FDA0002518648950000011
Figure FDA0002518648950000012
将个性化正畸弓丝曲线成形控制点信息集中各弯制点的坐标ti=(xi,yi,zi)′中的zi赋值为0,即令zi=0,获得正畸弓丝曲线转换平面正畸弓丝曲线T′;
Figure FDA0002518648950000013
步骤二、计算等半径确定圆域初始试划分个数:
按照
Figure FDA0002518648950000014
预先计算正畸弓丝曲线上i个弯制点的弯制点角距比,其中第j个弯制点的弯制点角距比规定
Figure FDA0002518648950000015
Ej是对第j个弯制点弯制复杂程度的量化描述,αj为作用在弯制点tj处的弯制角度,
Figure FDA0002518648950000016
表示作用在弯制点tj处弯制距离,即弯制点tj-1与tj之间曲线段的长度,特别地,由于第一个弯制点t1无需弯制,规定弯制点t1的弯制点角距比E1=0;首次将试划分n个等半径确定圆域,n的初始值为n=[∑E/(∑E)max]+1,其中[∑E/(∑E)max]表示对式∑E/(∑E)max所计算结果的取整,(∑E)max表示在正畸弓丝曲线上将要划分的任意一个等半径确定圆域an内的圆域弯制点角距比和
Figure FDA0002518648950000017
所要求的上限值,
Figure FDA0002518648950000018
表示正畸弓丝曲线上第n个等半径确定圆域an内的所有弯制点的弯制点角距比之和,即对等半径确定圆域an所划分的弯制点整体弯制复杂程度的量化描述,当变半径划分圆域an内的弯制点分别为tq+1、tq+2、tq+3、…、
Figure FDA0002518648950000021
时,规定
Figure FDA0002518648950000022
圆域弯制点个数
Figure FDA0002518648950000023
为半径为
Figure FDA0002518648950000024
的等半径确定圆域an内的弯制点个数,q表示正畸弓丝曲线上在等半径确定圆域an之前生成的n-1个圆域内所有弯制点的个数,即
Figure FDA0002518648950000025
以上依次提及的弯制点角距比E、圆域弯制点角距比和
Figure FDA0002518648950000026
圆域弯制点个数
Figure FDA0002518648950000027
三个参数统称为圆域限制参数,跳转至步骤三;
步骤三、试划分等半径确定圆域:
在以第一个弯制点t1为起点、以最后一个弯制点ti为终点的正畸弓丝曲线段上选取n+1个点作为圆域形成点,首个圆域形成点为弯制点t1所在的点,末尾圆域形成点为弯制点ti所在的点,使得每个圆域形成点与相邻圆域形成点间连接所得的n条直线段长度都相等,规定水平向右的向量按顺时针方向扫过的直线段依次用
Figure FDA0002518648950000028
表示,且存在
Figure FDA0002518648950000029
试划分等半径确定圆域以弯制点t1所在的圆域形成点为起始点进行,依次以
Figure FDA00025186489500000210
的中点作为圆心,以
Figure FDA00025186489500000211
作为半径生成n个等半径确定圆域,每个等半径确定圆域的边界线通过两个圆域形成点,相邻的两个等半径确定圆域的边界线相交于一个共用的圆域形成点,即第n-1个等半径确定圆域an-1右端的圆域形成点恰好是第n个等半径确定圆域an左端的圆域形成点,规定圆域边界线所截的正畸弓丝曲线段所包含的弯制点被该圆域所划分,当两个等半径确定圆域的边界线共用的圆域形成点所在的点恰好为正畸弓丝曲线上的一个弯制点,则规定该交点所在的弯制点被前一个等半径确定圆域所划分,假如第n-1个等半径确定圆域an-1与第n个等半径确定圆域an共用的圆域形成点为弯制点tj,则弯制点tj被等半径确定圆域an-1所划分,跳转至步骤四;
步骤四、寻找最佳试划分个数:
分别计算步骤三中生成的n个等半径确定圆域的圆域弯制点个数
Figure FDA00025186489500000212
得到圆域弯制点个数集
Figure FDA00025186489500000213
对圆域弯制点个数集Q中的n个圆域弯制点个数进行降序排列,取出其中最大的圆域弯制点个数,记为Qam,根据所要求的圆域弯制点个数上限值Qmax,Qmax表示在正畸弓丝曲线上将要划分的任意一个等半径确定圆域an内的圆域弯制点个数
Figure FDA00025186489500000214
所要求的上限值,特别地,Qmax=5,判断是否存在Qam≤5,
具体为:
如果Qam≤5不成立,说明所生成的n个等半径确定圆域中存在不符合圆域弯制点个数上限值Qmax要求的圆域,其中Qmax=5,可知此时n值不是最佳试划分个数,则需要通过改变圆域的个数,从而改变圆域半径的大小,进行重新试划分等半径确定圆域,令n=n+1,即下一次试划分等半径确定圆域时在此次划分个数的基础上增加一个,进而跳转至步骤三;
如果Qam≤5成立,说明所生成的n个等半径确定圆域都符合圆域弯制点个数上限值Qmax的要求,进而按照
Figure FDA0002518648950000031
计算步骤三中生成的n个等半径确定圆域的圆域弯制点角距比和
Figure FDA0002518648950000032
得到圆域弯制点角距比和集
Figure FDA0002518648950000033
对圆域弯制点角距比和集M中的n个圆域弯制点角距比和进行降序排列,取出其中最大的圆域弯制点角距比和,记为(∑E)am,根据所要求的圆域弯制点角距比和上限值(∑E)max,在Qam≤5成立的情况下判断是否存在(∑E)am≤(∑E)max
具体为:
如果(∑E)am≤(∑E)max成立,说明步骤三中所生成的n个等半径确定圆域都符合圆域弯制点角距比和上限值(∑E)max的要求,即此时所有的等半径确定圆域均符合划分要求,可知此时n值恰好为最佳试划分个数,称n个圆域a1、a2、…、an都为合理角距比等半径圆域,跳转至步骤五;
如果(∑E)am≤(∑E)max不成立,则说明步骤三中所生成的n个等半径确定圆域存在不符合圆域弯制点角距比和上限值(∑E)max要求的圆域,可知此时n值不是最佳试划分个数,则需要通过改变圆域的个数,从而改变圆域半径的大小,进行重新试划分等半径确定圆域,令n=n+1,即下一次试划分等半径确定圆域时在此次划分个数的基础上增加一个,进而跳转至步骤三;
步骤五、输出合理角距比等半径圆域划分半径
得到步骤四中输出的n个合理角距比等半径圆域及每个合理角距比等半径圆域的划分半径,划分半径值依次为
Figure FDA0002518648950000034
Figure FDA0002518648950000035
则requal即为在正畸弓丝曲线上划分出n个合理角距比等半径圆域的通用划分半径,输出合理角距比等半径圆域划分半径requal,程序结束。
CN202010484782.6A 2020-06-01 2020-06-01 一种基于正畸弓丝弯制点角距比和的平面等半径圆域划分半径确定方法 Active CN111588497B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010484782.6A CN111588497B (zh) 2020-06-01 2020-06-01 一种基于正畸弓丝弯制点角距比和的平面等半径圆域划分半径确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010484782.6A CN111588497B (zh) 2020-06-01 2020-06-01 一种基于正畸弓丝弯制点角距比和的平面等半径圆域划分半径确定方法

Publications (2)

Publication Number Publication Date
CN111588497A true CN111588497A (zh) 2020-08-28
CN111588497B CN111588497B (zh) 2021-06-04

Family

ID=72186103

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010484782.6A Active CN111588497B (zh) 2020-06-01 2020-06-01 一种基于正畸弓丝弯制点角距比和的平面等半径圆域划分半径确定方法

Country Status (1)

Country Link
CN (1) CN111588497B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115040275A (zh) * 2022-04-18 2022-09-13 哈尔滨理工大学 一种基于空间平移子坐标系卦限判断的正畸弓丝评价方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3729824A (en) * 1970-12-07 1973-05-01 G A C International Inc Orthodontic arch wire
CN107647925A (zh) * 2017-11-14 2018-02-02 哈尔滨理工大学 一种用于正畸弓丝弯制规划的圆域划分方法
CN107714203A (zh) * 2017-11-14 2018-02-23 哈尔滨理工大学 一种等角度划分的正畸弓丝弯制顺序规划方法
CN107822722A (zh) * 2017-11-14 2018-03-23 哈尔滨理工大学 一种用于正畸弓丝弯制运动规划的有限点寻回展成方法
US20180147026A1 (en) * 2007-03-14 2018-05-31 Orthoaccel Technologies, Inc. Pulsatile orthodontic device
CN110200710A (zh) * 2019-04-17 2019-09-06 广东工业大学 一种基于三维成像与实时建模的口腔修复方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3729824A (en) * 1970-12-07 1973-05-01 G A C International Inc Orthodontic arch wire
US20180147026A1 (en) * 2007-03-14 2018-05-31 Orthoaccel Technologies, Inc. Pulsatile orthodontic device
CN107647925A (zh) * 2017-11-14 2018-02-02 哈尔滨理工大学 一种用于正畸弓丝弯制规划的圆域划分方法
CN107714203A (zh) * 2017-11-14 2018-02-23 哈尔滨理工大学 一种等角度划分的正畸弓丝弯制顺序规划方法
CN107822722A (zh) * 2017-11-14 2018-03-23 哈尔滨理工大学 一种用于正畸弓丝弯制运动规划的有限点寻回展成方法
CN110200710A (zh) * 2019-04-17 2019-09-06 广东工业大学 一种基于三维成像与实时建模的口腔修复方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
姜金刚,韩英帅,张永德,刘怡,温馥嘉,蒋济雄: "机器人弯制正畸弓丝成形控制点规划及实验研究", 《仪器仪表学报》 *
左思浩: "基于人手仿生运动原理的正畸弓丝机器人弯制规划", 《中国优秀硕士学位论文全文数据库 信息科技辑》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115040275A (zh) * 2022-04-18 2022-09-13 哈尔滨理工大学 一种基于空间平移子坐标系卦限判断的正畸弓丝评价方法
CN115040275B (zh) * 2022-04-18 2023-06-20 哈尔滨理工大学 一种基于空间平移子坐标系卦限判断的正畸弓丝评价方法

Also Published As

Publication number Publication date
CN111588497B (zh) 2021-06-04

Similar Documents

Publication Publication Date Title
CN111588502B (zh) 一种基于弯制点角距比和的正畸弓丝变半径圆域划分方法
CN107647925B (zh) 一种用于正畸弓丝弯制规划的圆域划分方法
CN111588499B (zh) 一种基于正畸弓丝弯制点密度的平面等半径圆域划分半径确定方法
CN107714203B (zh) 一种等角度划分的正畸弓丝弯制顺序规划方法
CN107822722B (zh) 一种用于正畸弓丝弯制运动规划的有限点寻回展成方法
CN111588497B (zh) 一种基于正畸弓丝弯制点角距比和的平面等半径圆域划分半径确定方法
CN111588504B (zh) 一种用于正畸弓丝弯制顺序规划的空间变半径球域划分方法
CN111588496B (zh) 一种用于正畸弓丝弯制规划的平面变角度划分方法
CN111588494B (zh) 一种基于弯制点密度的正畸弓丝变角度划分方法
CN111588501B (zh) 一种正畸弓丝弯制规划等半径圆域划分半径确定方法
CN114972184B (zh) 一种基于权值比例法的正畸弓丝误差评价方法
CN111588503B (zh) 一种基于弯制点密度的正畸弓丝变半径圆域划分方法
CN111588491B (zh) 一种基于正畸弓丝弯制点密度的空间等半径球域划分半径确定方法
CN111588505B (zh) 一种用于正畸弓丝弯制顺序规划的平面变半径圆域划分方法
CN115035196A (zh) 一种基于弯制点复杂度判断的正畸弓丝误差率评价方法
CN115024840A (zh) 一种基于归一化弯制点密度的正畸弓丝误差率评价方法
CN115024841A (zh) 一种基于复杂度区分法的正畸弓丝误差评价方法
CN115040275A (zh) 一种基于空间平移子坐标系卦限判断的正畸弓丝评价方法
CN111588493B (zh) 一种基于弯制点角距比和的正畸弓丝变角度划分方法
CN111588500B (zh) 一种正畸弓丝弯制顺序规划等角度划分角度确定方法
CN111588495B (zh) 一种基于正畸弓丝弯制点单位角距比和的等角度划分角度确定方法
CN111588498B (zh) 一种基于正畸弓丝弯制点密度的等角度划分角度确定方法
CN117093819B (zh) 一种用于正畸弓丝弯制规划的弯制参数权重因子及限制参数确定方法
CN114943058B (zh) 一种基于位置误差判断的正畸弓丝误差波动度评价方法
CN115690009A (zh) 一种基于共面等角向量的正畸弓丝评价方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant