CN111061213B - 一种基于Bezier曲线转角平滑过渡算法的加工方法 - Google Patents

一种基于Bezier曲线转角平滑过渡算法的加工方法 Download PDF

Info

Publication number
CN111061213B
CN111061213B CN201911240783.XA CN201911240783A CN111061213B CN 111061213 B CN111061213 B CN 111061213B CN 201911240783 A CN201911240783 A CN 201911240783A CN 111061213 B CN111061213 B CN 111061213B
Authority
CN
China
Prior art keywords
transition
curve
transition curve
corner
data points
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.)
Active
Application number
CN201911240783.XA
Other languages
English (en)
Other versions
CN111061213A (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.)
Tianjin University
Original Assignee
Tianjin University
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 Tianjin University filed Critical Tianjin University
Priority to CN201911240783.XA priority Critical patent/CN111061213B/zh
Publication of CN111061213A publication Critical patent/CN111061213A/zh
Application granted granted Critical
Publication of CN111061213B publication Critical patent/CN111061213B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B19/00Programme-control systems
    • G05B19/02Programme-control systems electric
    • G05B19/18Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form
    • G05B19/19Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form characterised by positioning or contouring control systems, e.g. to control position from one programmed point to another or to control movement along a programmed continuous path
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/30Computing systems specially adapted for manufacturing

Landscapes

  • Engineering & Computer Science (AREA)
  • Human Computer Interaction (AREA)
  • Manufacturing & Machinery (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Image Generation (AREA)

Abstract

本发明公开了一种基于Bezier曲线转角平滑过渡算法的加工方法,步骤包括:先用CAD/CAM软件获取小线段轨迹离散数据点;然后读入三个数据点,并标记中间点为拐点,计算拐点两侧微线段长度、转角以及单位方向矢量;之后针对圆心角的取值范围构建不同的具有G2连续特性的转角圆弧过渡矢量模型;再之后计算过渡误差以及曲率极值,并依据预设的近似误差约束和最大曲率限制确定转角两侧转接点的位置和过渡段长度,代入模型求解过渡曲线;从而完成一个转角的光顺过渡;重复读入数据点,依次完成所有转角的光顺过渡,并输出采用过渡曲线实现小线段间光滑衔接的加工路径。

Description

一种基于Bezier曲线转角平滑过渡算法的加工方法
技术领域
本发明属于机械工程领域,主要应用在智能制造行业的CNC数控加工,涉及一种用于小线段高速加工中的路径轨迹优化方法,特别涉及一种基于Bezier曲线转角平滑过渡算法的加工方法。
背景技术
目前,在航空航天、汽车轮船与高端模具制造等领域,对于具有复杂曲线曲面造型的零件加工,通常的做法是利用CAD/CAM软件将其表征为一系列巨量小线段路径轨迹后再进行数控加工,这种做法简单高效适用范围广、刀具路径生成速度较快,而且容易实现偏置处理和插补,因此仍然被绝大多数的数控系统广泛使用。而小线段形式的数控程序存在数据量大、数据传输易出错、加减速频繁易引起冲击、转折点多易造成过切与欠切,严重制约了复杂曲线曲面零件的高性能制造发展。所以,针对小线段的数控加工程序进行轨迹优化,即用平滑的曲线来替代转角,以消除转角处的切向不连续进而实现轨迹局部光顺过渡,对于减小速度和加速度波动、改善零件表面质量的光洁度以及提高机床的加工效率和运动平稳性有着十分重要的意义。
发明内容
本发明的目的旨在克服现有小线段路径直接插补技术的不足,提供一种基于Bezier曲线转角平滑过渡算法的加工方法。
本发明为解决公知技术中存在的技术问题所采取的技术方案是:一种基于Bezier曲线转角平滑过渡算法的加工方法,将原始的加工轮廓轨迹转化为一系列由连续小线段路径轨迹表征的离散数据点;依次读入三个相邻的数据点,并标记中间点为拐点,计算拐点两侧小线段长度、转角以及单位方向矢量;采用Bezier曲线构建具有G2连续特性的转角圆弧过渡曲线模型;并依据预设的过渡误差约束和最大曲率约束来确定转角两侧转接点的坐标位置和过渡曲线长度,将确定的坐标位置和过渡曲线长度代入转角圆弧过渡曲线模型来求解过渡曲线,从而完成一个转角的光顺过渡;重复读入数据点,依次完成所有转角的光顺过渡,输出采用过渡曲线实现小线段间光滑衔接的加工路径。
进一步地,将原始的加工轮廓轨迹转化为一系列由连续小线段路径轨迹表征的离散数据点的方法为:先利用CAD软件设计零件模型,之后利用CAM软件采用等误差法将该零件模型细化为顺序直线连接成多个连续小线段的离散数据点,并存储在txt文件中。
进一步地,求解过渡曲线的算法用C语言编程实现,嵌入到数控系统当中,并将最终输出的加工路径轨迹转化为数控系统所能识别的NC代码,进行数控加工。
进一步地,采用三次有理Bezier曲线构建具有G2连续特性的转角圆弧曲线矢量模型。
进一步地,采用Bezier曲线构建具有G2连续特性的转角圆弧过渡曲线模型的方法为:
设三个相邻的离散数据点分别为Q0、Q1、Q2,θ为线段Q0Q1与Q1Q2的转角,预设过渡误差目标值ε和最大曲率目标值K;
当θ∈(0,π/2]时,构造一条使线段Q0Q1与Q1Q2的转角平滑过渡的过渡曲线模型A,设P0、P12、P3为过渡曲线模型A的三个控制点,P0在线段Q0Q1上,P3在线段Q1Q2上,P12和Q1重合,T1、T2分别为经过P0和P3的过渡曲线模型A在P0和P3处的单位切矢,O1为过渡曲线模型A的圆心,设O1、Q1两点连接成的直线与过渡曲线模型A交于E,使EO1的长度≤ε,且过渡曲线模型A的曲率≤K;
当θ∈(π/2,π)时,构造一条使线段Q0Q1与Q1Q2的转角平滑过渡的过渡曲线模型B,过渡曲线模型B包括前半段过渡曲线B1和后半段过渡曲线B2,设
Figure BDA0002306161490000021
P3 1为前半段过渡曲线B1的三个控制点;设
Figure BDA0002306161490000022
P3 2为后半段过渡曲线B2的三个控制点,且P3 1
Figure BDA0002306161490000023
重合,
Figure BDA0002306161490000024
在线段Q0Q1上,
Figure BDA0002306161490000025
P3 2在线段Q1Q2上,T1 1
Figure BDA0002306161490000026
为经过
Figure BDA0002306161490000027
和P3 1的前半段过渡曲线B1
Figure BDA0002306161490000028
和P3 1处的单位切矢;T1 2
Figure BDA0002306161490000029
为经过
Figure BDA00023061614900000210
和P3 2的后半段过渡曲线B2
Figure BDA00023061614900000211
和P3 2处的单位切矢;
Figure BDA00023061614900000212
与T1 2相等;O2为前半段过渡曲线B1和后半段过渡曲线B2的共同圆心,P3 1在O2、Q1两点连接成的直线上,使P3 1、Q1两点的长度≤ε;且前半段过渡曲线B1和后半段过渡曲线B2的曲率≤K。
进一步地,该方法包括如下具体步骤:
步骤一、获得一系列形成连续小线段加工轨迹的离散数据点,并存储在txt文件中;
步骤二、从起始点开始依次连续读取三个数据点信息,将第一个数据点作为首点;
步骤三、读取的三个数据点中,中间数据点标记为拐点,并计算拐点两侧相邻小线段的长度l1、l2、单位矢量T1、T2以及两矢量之间的夹角θ;
步骤四、根据θ的不同,构建不同的具有G2连续特性的转角圆弧过渡曲线模型,当θ∈(0,π/2]时,构建过渡曲线模型A,当θ∈(π/2,π)时,构建过渡曲线模型B;
步骤五、针对不同的过渡曲线矢量模型,分别计算原始轮廓与过渡曲线之间在拐点处的过渡误差ε以及过渡曲线的最大曲率κ,并依据预设的过渡误差约束和最大曲率约束来确定转角两侧转接点的坐标位置和过渡段长度,然后代入模型,即可唯一确定过渡曲线;
步骤六、判断未读取数据点个数是否≥1;
若是,舍弃三个数据点的第一个数据点,保留后两个数据点,并将之前标记为拐点的数据点作为新的首点,继续读取下一个数据点的位置矢量信息,返回执行步骤三;
若否,则结束算法,并输出经过光顺过渡处理后的加工路径。
步骤七、将上述算法流程用C语言编程实现,嵌入到数控系统当中,并将算法最终输出的曲率连续光顺路径转化为数控系统所能识别的NC代码,进行数控加工。
本发明具有的优点和积极效果是:本发明可应用在中高档数控系统中,其具备误差可控,且能实现路径局部光顺和转折点平滑过渡,并为下一步的在线速度规划和实时插补提供便利。现有的技术方法相比,本发明的优势和积极效果是:此方法能够有效地解决因小线段路径中相邻线段衔接处不具G1及以上连续而造成进给速度和加速度突变的问题,为国产的中高档数控机床,提供了一种可以直接集成到数控系统当中的具备误差可控且能在线路径光顺的参数曲线转角平滑过渡算法。
附图说明
图1是本发明所采用的标准型的三次有理Bezier曲线简图;
图2是本发明所采用的三次有理Bezier曲线简图;
图3是本发明相邻小线段的转角求解示意图;
图4是本发明有理Bezier曲线圆弧示意图;
图5是本发明当θ∈(0,π/2]时的圆弧过渡模型构建示意图;
图6是本发明当θ∈(π/2,π)时的圆弧过渡模型构建示意图;
图7是本发明当θ∈(0,π/2]时G2连续Bezier过渡曲线的曲率梳状线;
图8是本发明当θ∈(π/2,π)时G2连续Bezier过渡曲线的曲率梳状线;
图9是本发明近似误差约束下过渡长度的求解示意图;
图10是本发明相邻小线段间转角平滑过渡算法的具体实现流程框图。
图中的代号、符号说明如下:
图1、2中:P0、P12、P3为三次有理Bezier曲线的控制点,w12为中间控制点P12所对应的权因子,向量V0、V3为曲线的首末两点出的切向矢量,Ps为曲线上的肩点,Pm为投影点;
图4中:θ为有理Bezier曲线描述一段圆弧时所对应的圆心角,R为半径,A、B、C、D分别为三个控制点以及肩点。
图5、6中:Q0、Q1、Q2为连续的三个离散数据点,中间点记为拐点,P0、P12、P3为过渡曲线模型A的三个控制点,T1、T2分为过渡曲线在P0和P3处的单位切矢,θ为线段Q0Q1与Q1Q2的转角,其值同圆心角大小一样,O1、O2为圆心,ε为过渡误差。
Figure BDA0002306161490000041
P3 1
Figure BDA0002306161490000042
P3 2分别为模型B中前半段和后半段Bezier圆弧过渡曲线的三个控制点,T1 1、T1 2
Figure BDA0002306161490000043
分别为模型B中过渡曲线在
Figure BDA0002306161490000044
P3 1
Figure BDA0002306161490000045
P3 2处的单位切矢。
图7、8、9中:Qi为离散数据点,Pi为控制点,κ(u)为曲率,ε为过渡误差,
Figure BDA0002306161490000046
为拐点到弦Pi-1Pi+1的距离。
具体实施方式
为能进一步了解本发明的发明内容、特点及功效,兹列举以下实施例,并配合附图详细说明如下:
请参见图1至图10,一种基于Bezier曲线转角平滑过渡算法的加工方法,将原始的加工轮廓轨迹转化为一系列由连续小线段路径轨迹表征的离散数据点;依次读入三个相邻的数据点,并标记中间点为拐点,计算拐点两侧小线段长度、转角以及单位方向矢量;采用Bezier曲线构建具有G2连续特性的转角圆弧过渡曲线模型;并依据预设的过渡误差约束和最大曲率约束来确定转角两侧转接点的坐标位置和过渡曲线长度,将确定的坐标位置和过渡曲线长度代入转角圆弧过渡曲线模型来求解过渡曲线,从而完成一个转角的光顺过渡;重复读入数据点,依次完成所有转角的光顺过渡,输出采用过渡曲线实现小线段间光滑衔接的加工路径。
优选地,将原始的加工轮廓轨迹转化为一系列由连续小线段路径轨迹表征的离散数据点的方法可为:可先利用CAD软件设计零件模型,之后可利用CAM软件采用等误差法将该零件模型细化为顺序直线连接成多个连续小线段的离散数据点,并可存储在txt文件中。
优选地,求解过渡曲线的算法可用C语言编程实现,可嵌入到数控系统当中,并可将最终输出的加工路径轨迹转化为数控系统所能识别的NC代码,进行数控加工。
优选地,可采用三次有理Bezier曲线构建具有G2连续特性的转角圆弧曲线矢量模型。
优选地,采用Bezier曲线构建具有G2连续特性的转角圆弧过渡曲线模型的方法可为:
可设三个相邻的离散数据点分别为Q0、Q1、Q2,θ为线段Q0Q1与Q1Q2的转角,预设过渡误差目标值ε和最大曲率目标值K;
当θ∈(0,π/2]时,可构造一条使线段Q0Q1与Q1Q2的转角平滑过渡的过渡曲线模型A,可设P0、P12、P3为过渡曲线模型A的三个控制点,P0在线段Q0Q1上,P3在线段Q1Q2上,P12和Q1重合,T1、T2分别为经过P0和P3的过渡曲线模型A在P0和P3处的单位切矢,O1为过渡曲线模型A的圆心,可设O1、Q1两点连接成的直线与过渡曲线模型A交于E,可使EO1的长度≤ε,且过渡曲线模型A的曲率≤K;
当θ∈(π/2,π)时,可构造一条使线段Q0Q1与Q1Q2的转角平滑过渡的过渡曲线模型B,过渡曲线模型B包括前半段过渡曲线B1和后半段过渡曲线B2,可设
Figure BDA0002306161490000051
P3 1为前半段过渡曲线B1的三个控制点;可设
Figure BDA0002306161490000052
P3 2为后半段过渡曲线B2的三个控制点,且P3 1和P0 2重合,
Figure BDA0002306161490000053
在线段Q0Q1上,
Figure BDA0002306161490000054
P3 2在线段Q1Q2上,T1 1
Figure BDA0002306161490000055
为经过
Figure BDA0002306161490000056
和P3 1的前半段过渡曲线B1
Figure BDA0002306161490000057
和P3 1处的单位切矢;T1 2
Figure BDA0002306161490000058
为经过
Figure BDA0002306161490000059
和P3 2的后半段过渡曲线B2
Figure BDA00023061614900000510
和P3 2处的单位切矢;
Figure BDA00023061614900000511
与T1 2相等;O2为前半段过渡曲线B1和后半段过渡曲线B2的共同圆心,P3 1在O2、Q1两点连接成的直线上,可使P3 1、Q1两点的长度≤ε;且前半段过渡曲线B1和后半段过渡曲线B2的曲率≤K。
优选地,该方法可包括如下具体步骤:
步骤一、可获得一系列形成连续小线段加工轨迹的离散数据点,并存储在txt文件中;
步骤二、可从起始点开始依次连续读取三个数据点信息,将第一个数据点作为首点;
步骤三、读取的三个数据点中,中间数据点可标记为拐点,并计算拐点两侧相邻小线段的长度l1、l2、单位矢量T1、T2以及两矢量之间的夹角θ;
步骤四、可根据θ的不同,构建不同的具有G2连续特性的转角圆弧过渡曲线模型,当θ∈(0,π/2]时,构建过渡曲线模型A,当θ∈(π/2,π)时,构建过渡曲线模型B;
步骤五、针对不同的过渡曲线矢量模型,可分别计算原始轮廓与过渡曲线之间在拐点处的过渡误差ε以及过渡曲线的最大曲率κ,并可依据预设的过渡误差约束和最大曲率约束来确定转角两侧转接点的坐标位置和过渡段长度,然后代入模型,即可唯一确定过渡曲线;
步骤六、判断未读取数据点个数是否≥1;
若是,可舍弃三个数据点的第一个数据点,保留后两个数据点,并可将之前标记为拐点的数据点作为新的首点,继续读取下一个数据点的位置矢量信息,返回执行步骤三;
若否,则可结束算法,并可输出经过光顺过渡处理后的加工路径。
步骤七、将上述算法流程用C语言编程实现,嵌入到数控系统当中,并将算法最终输出的曲率连续光顺路径转化为数控系统所能识别的NC代码,进行数控加工。
以下结合本发明的一个优选实施例来说明本发明的工作原理:
一种基于Bezier曲线转角平滑过渡算法的小线段轨迹优化方法,该方法用于小线段高速加工,该方法基于三次有理Bezier曲线转角平滑过渡算法,该方法的工作原理说明如下:
请见图1至图4,该方法基于三次有理Bezier曲线,标准型的三次有理Bezier曲线方程一般表达式为:
Figure BDA0002306161490000061
式中:
P0、P1、P2、P3为:三次标准型的有理Bezier曲线的四个控制点;
w0、w1、w2、w3为:控制点P0、P1、P2、P3所对应的权因子;
u为Bezier曲线的参数。
调节中间两个控制点,使得P1=P2=P12,如图2所示,变成了三次有理Bezier曲线的特例,取w0=w3=1,w1=w2=w12>0,代入式上式,则本文所采用的三次有理Bezier曲线参数表达式为:
Figure BDA0002306161490000071
式中:
P12为:中间控制点;w12为:中间控制点P12所对应的权因子。
如图4所示,三次有理Bezier曲线C(u)恰好表示一段圆弧时,中间的未知权因子w12以及过渡段长的约束关系应同时满足以下条件:
Figure BDA0002306161490000072
||AC||=||BC||、w0=w3=1。
该方法的具体步骤可如下:
步骤一、利用CAD/CAM软件,并采用等误差法将原始零件模型细化为顺序直线连接成多个连续小线段的离散数据点,误差大小可参考加工精度;
步骤二、可从起始点开始依次连续读取三个数据点信息,将第一个数据点作为首点;
步骤三、读取的三个数据点中,中间数据点可标记为拐点,并计算拐点两侧相邻小线段的长度l1、l2、单位矢量T1、T2以及两矢量之间的夹角θ;
如图4以及图5、6所示,l1、l2、T1、T2、θ求解公式如下:
Figure BDA0002306161490000073
其中,l1=Q0Q1、l2=Q1Q2
Figure BDA0002306161490000074
步骤四、可根据θ的不同,构建不同的具有G2连续特性的转角圆弧曲线矢量模型,当θ∈(0,π/2]时,构建过渡曲线模型A,当θ∈(π/2,π)时,构建过渡曲线模型B;
可设三个相邻的离散数据点分别为Q0、Q1、Q2,θ为线段Q0Q1与Q1Q2的转角,预设过渡误差目标值ε和最大曲率目标值K;
当θ∈(0,π/2]时,可构造一条使线段Q0Q1与Q1Q2的转角平滑过渡的过渡曲线模型A,可设P0、P12、P3为过渡曲线模型A的三个控制点,P0在线段Q0Q1上,P3在线段Q1Q2上,P12和Q1重合,T1、T2分别为经过P0和P3的过渡曲线模型A在P0和P3处的单位切矢,O1为过渡曲线模型A的圆心,可设O1、Q1两点连接成的直线与过渡曲线模型A交于E,可使EO1的长度≤ε,且过渡曲线模型A的曲率≤K;
请参见图5,图中曲线对应过渡曲线模型A,过渡曲线模型A的矢量表达式如下:
Figure BDA0002306161490000081
Figure BDA0002306161490000082
式中:
θ为:线段Q0Q1与Q1Q2的转角,其值刚好等于圆弧过渡曲线所对应的圆心角大小;
u为:Bezier曲线的参数;
h为:过渡曲线的过渡段长度;
P0、P1、P2、P3为:三次有理Bezier过渡曲线的四个控制点,其中P1、P2重合并重新记为P12,称作中间控制点;
Q1为:拐点,即:线段Q0Q1与Q1Q2的交点;
T1、T2为:过渡曲线在P0和P3处的单位切矢。
当θ∈(π/2,π)时,可构造一条使线段Q0Q1与Q1Q2的转角平滑过渡的过渡曲线模型B,过渡曲线模型B包括前半段过渡曲线B1和后半段过渡曲线B2,可设
Figure BDA0002306161490000083
P3 1为前半段过渡曲线B1的三个控制点;可设
Figure BDA0002306161490000084
P3 2为后半段过渡曲线B2的三个控制点,且P3 1
Figure BDA0002306161490000085
重合,
Figure BDA0002306161490000086
在线段Q0Q1上,
Figure BDA0002306161490000087
P3 2在线段Q1Q2上,T1 1
Figure BDA0002306161490000088
为经过
Figure BDA0002306161490000089
和P3 1的前半段过渡曲线B1
Figure BDA00023061614900000810
和P3 1处的单位切矢;T1 2
Figure BDA00023061614900000811
为经过
Figure BDA00023061614900000812
和P3 2的后半段过渡曲线B2
Figure BDA00023061614900000813
和P3 2处的单位切矢;
Figure BDA00023061614900000814
与T1 2相等;O2为前半段过渡曲线B1和后半段过渡曲线B2的共同圆心,P3 1在O2、Q1两点连接成的直线上,可使P3 1、Q1两点的长度≤ε;且前半段过渡曲线B1和后半段过渡曲线B2的曲率≤K。
请参见图6,图中曲线对应过渡曲线模型B,过渡曲线模型B的矢量表达式如下:
Figure BDA00023061614900000815
Figure BDA0002306161490000091
式中:
θ为:线段Q0Q1与Q1Q2的转角,其值刚好等于圆弧过渡曲线所对应的圆心角大小;
u为:Bezier曲线的参数;
h为:过渡曲线的过渡段长度;
Figure BDA0002306161490000092
P3 1为:前半段过渡曲线的三个控制点;
Figure BDA0002306161490000093
P3 2为:后半段过渡曲线的三个控制点;
T1 1
Figure BDA0002306161490000094
为:前半段过渡曲线在
Figure BDA0002306161490000095
和P3 1处的单位切矢;
T1 2
Figure BDA0002306161490000096
为:后半段过渡曲线在
Figure BDA0002306161490000097
和P3 2处的单位切矢。
步骤五、针对不同的过渡模型,分别计算原始轮廓与过渡曲线之间在拐点处的过渡误差ε以及过渡曲线上的最大曲率κ,并依据预设的过渡误差约束和最大曲率约束来确定转角两侧转接点的位置和过渡段长度,然后代入模型,即可唯一确定过渡曲线,执行下一步。
过渡误差约束以及最大曲率约束条件下过渡长度的求解原理请参见图7至图9。
本发明算法中所采用的转角圆弧过渡曲线,其过渡误差以及曲率极值均可解析表达,表达式如下,
Figure BDA0002306161490000098
Figure BDA0002306161490000099
κ2(0)=κ2(0.5)=κ2(1)=0
Figure BDA00023061614900000910
ε2=tan(θ/4)*h,θ∈(π/2,π)
式中:
κ1(u)|u=0:0.5:1为:过渡曲线模型A上,曲线参数u=0,0.5,1时的曲率;
κ2(u)|u=0:0.25:1为:过渡曲线模型B上,曲线参数u=0,0.25,0.5,0.75,1时的曲率;
θ为:过渡曲线所对应的圆心角;h为:过渡段长度;
ε1、ε2为:不同的圆心角取值范围求解出来的过渡误差。
所以当θ∈(0,π/2]以及θ∈(π/2,π)时,过渡长度h的取值范围分别为:
Figure BDA0002306161490000101
Figure BDA0002306161490000102
断点两侧曲线转接点的位置由过渡长度h、中间拐点处的位置矢量坐标以及两侧小线段的单位方向矢量来确定。
步骤六、判断待读数据点个数是否≥1;
若是,则继续读取下一个数据点的位置矢量信息,舍弃第一个数据点保留后两个数据点,并将之前标记的中间拐点作为新的首点,返回执行步骤三。
若否,则结束算法,并输出经过光顺过渡处理后的加工路径。
步骤七、将上述算法流程用C语言编程实现,嵌入到数控系统当中,并将算法最终输出的曲率连续光顺路径转化为数控系统所能识别的NC代码,进行数控加工。
以上所述的实施例仅用于说明本发明的技术思想及特点,其目的在于使本领域内的技术人员能够理解本发明的内容并据以实施,不能仅以本实施例来限定本发明的专利范围,即凡本发明所揭示的精神所作的同等变化或修饰,仍落在本发明的专利范围内。

Claims (5)

1.一种基于Bezier曲线转角平滑过渡算法的加工方法,其特征在于,将原始的加工轮廓轨迹转化为一系列由连续小线段路径轨迹表征的离散数据点;依次读入三个相邻的数据点,并标记中间点为拐点,计算拐点两侧小线段长度、转角以及单位矢量;采用Bezier曲线构建具有G2连续特性的转角圆弧过渡曲线模型;并依据预设的过渡误差约束和最大曲率约束来确定转角两侧转接点的坐标位置和过渡曲线长度,将确定的坐标位置和过渡曲线长度代入转角圆弧过渡曲线模型来求解过渡曲线,从而完成一个转角的光顺过渡;重复读入数据点,依次完成所有转角的光顺过渡,输出采用过渡曲线实现小线段间光滑衔接的加工路径;该方法包括如下具体步骤:
步骤一、获得一系列形成连续小线段加工轨迹的离散数据点,并存储在txt文件中;
步骤二、从起始点开始依次连续读取三个数据点信息,将第一个数据点作为首点;
步骤三、读取的三个数据点中,中间数据点标记为拐点,并计算拐点两侧相邻小线段的长度l1、l2、单位矢量T1、T2以及两矢量之间的转角θ;
步骤四、根据θ的不同,构建不同的具有G2连续特性的转角圆弧过渡曲线模型,当θ∈(0,π/2]时,构建过渡曲线模型A,当θ∈(π/2,π)时,构建过渡曲线模型B;
步骤五、针对不同的过渡曲线矢量模型,分别计算原始轮廓与过渡曲线之间在拐点处的过渡误差ε以及过渡曲线的最大曲率κ,并依据预设的过渡误差约束和最大曲率约束来确定转角两侧转接点的坐标位置和过渡段长度,然后代入模型,即可唯一确定过渡曲线;
步骤六、判断未读取数据点个数是否≥1;
若是,舍弃三个数据点的第一个数据点,保留后两个数据点,并将之前标记为拐点的数据点作为新的首点,继续读取下一个数据点的位置矢量信息,返回执行步骤三;
若否,则结束算法,并输出经过光顺过渡处理后的加工路径;
步骤七、将上述算法流程用C语言编程实现,嵌入到数控系统当中,并将算法最终输出的曲率连续光顺路径转化为数控系统所能识别的NC代码,进行数控加工。
2.根据权利要求1所述的基于Bezier曲线转角平滑过渡算法的加工方法,其特征在于,将原始的加工轮廓轨迹转化为一系列由连续小线段路径轨迹表征的离散数据点的方法为:先利用CAD软件设计零件模型,之后利用CAM软件采用等误差法将该零件模型细化为顺序直线连接成多个连续小线段的离散数据点,并存储在txt文件中。
3.根据权利要求1所述的基于Bezier曲线转角平滑过渡算法的加工方法,其特征在于,求解过渡曲线的算法用C语言编程实现,嵌入到数控系统当中,并将最终输出的加工路径轨迹转化为数控系统所能识别的NC代码,进行数控加工。
4.根据权利要求1所述的基于Bezier曲线转角平滑过渡算法的加工方法,其特征在于,采用三次有理Bezier曲线构建具有G2连续特性的转角圆弧曲线矢量模型。
5.根据权利要求1所述的基于Bezier曲线转角平滑过渡算法的加工方法,其特征在于,采用Bezier曲线构建具有G2连续特性的转角圆弧过渡曲线模型的方法为:
设三个相邻的离散数据点分别为Q0、Q1、Q2,θ'为线段Q0Q1与Q1Q2的转角,预设过渡误差ε和最大曲率目标值K;
当θ'∈(0,π/2]时,构造一条使线段Q0Q1与Q1Q2的转角平滑过渡的过渡曲线模型A,设P0、P12、P3为过渡曲线模型A的三个控制点,P0在线段Q0Q1上,P3在线段Q1Q2上,P12和Q1重合,T1'、T2'分别为经过P0和P3的过渡曲线模型A在P0和P3处的单位切矢,O1为过渡曲线模型A的圆心,设O1、Q1两点连接成的直线与过渡曲线模型A交于E,使EO1的长度≤ε,且过渡曲线模型A的曲率≤K;
当θ'∈(π/2,π)时,构造一条使线段Q0Q1与Q1Q2的转角平滑过渡的过渡曲线模型B,过渡曲线模型B包括前半段过渡曲线B1和后半段过渡曲线B2,设
Figure FDA0003644498060000021
为前半段过渡曲线B1的三个控制点;设
Figure FDA0003644498060000022
为后半段过渡曲线B2的三个控制点,且P3 1
Figure FDA00036444980600000217
重合,
Figure FDA0003644498060000024
在线段Q0Q1上,
Figure FDA0003644498060000025
在线段Q1Q2上,T1 1
Figure FDA0003644498060000026
为经过
Figure FDA0003644498060000027
和P3 1的前半段过渡曲线B1
Figure FDA0003644498060000028
Figure FDA00036444980600000216
处的单位切矢;T1 2
Figure FDA0003644498060000029
为经过
Figure FDA00036444980600000210
和P3 2的后半段过渡曲线B2
Figure FDA00036444980600000211
Figure FDA00036444980600000212
处的单位切矢;
Figure FDA00036444980600000213
与T1 2相等;O2为前半段过渡曲线B1和后半段过渡曲线B2的共同圆心,
Figure FDA00036444980600000214
在O2、Q1两点连接成的直线上,使
Figure FDA00036444980600000215
Q1两点的长度≤ε;且前半段过渡曲线B1和后半段过渡曲线B2的曲率≤K。
CN201911240783.XA 2019-12-04 2019-12-04 一种基于Bezier曲线转角平滑过渡算法的加工方法 Active CN111061213B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911240783.XA CN111061213B (zh) 2019-12-04 2019-12-04 一种基于Bezier曲线转角平滑过渡算法的加工方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911240783.XA CN111061213B (zh) 2019-12-04 2019-12-04 一种基于Bezier曲线转角平滑过渡算法的加工方法

Publications (2)

Publication Number Publication Date
CN111061213A CN111061213A (zh) 2020-04-24
CN111061213B true CN111061213B (zh) 2022-08-09

Family

ID=70300068

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911240783.XA Active CN111061213B (zh) 2019-12-04 2019-12-04 一种基于Bezier曲线转角平滑过渡算法的加工方法

Country Status (1)

Country Link
CN (1) CN111061213B (zh)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114115112B (zh) * 2020-09-01 2024-02-09 大族激光科技产业集团股份有限公司 数控加工方法及装置
CN112276949B (zh) * 2020-10-21 2022-03-11 哈工大机器人(合肥)国际创新研究院 一种相邻关节空间-笛卡尔空间轨迹过渡方法及装置
CN112356025B (zh) * 2020-10-29 2022-05-03 中国科学院软件研究所 一种用双圆弧近似nurbs曲线的机器人运动方法和装置
CN112346406B (zh) * 2020-11-03 2022-03-01 西北工业大学 五轴数控机床刀具轨迹的光顺方法
CN112276957B (zh) * 2020-11-05 2022-03-11 哈工大机器人(合肥)国际创新研究院 一种直线段和圆弧的平滑过渡方法及系统
CN112965443B (zh) * 2021-02-08 2022-03-01 浙江工业大学 一种裁床拐角轨迹跟踪高精度插值控制方法
CN113341886B (zh) * 2021-05-25 2022-12-27 深圳市汇川技术股份有限公司 平滑进给速度规划方法、设备及计算机可读存储介质
CN113359607B (zh) * 2021-07-06 2022-08-09 深圳数马电子技术有限公司 一种应用于五轴数控机的拐角过渡的轨迹确定方法
CN113589761B (zh) * 2021-07-16 2024-06-07 天津大学 一种基于三角函数加减速控制的直接速度过渡光顺方法
CN113467384B (zh) * 2021-07-30 2022-09-06 深圳数马电子技术有限公司 一种应用于五轴数控机床的拐角过渡方法
CN113791576B (zh) * 2021-08-19 2023-10-17 五邑大学 一种轨迹间的局部光顺过渡方法、装置、设备及存储介质
CN113759830B (zh) * 2021-09-02 2023-04-18 浙江大学 基于等效加速度的线性路径数控加工进给速度控制方法
CN113900413B (zh) * 2021-11-24 2024-06-18 南京达风数控技术有限公司 一种数控系统平顺速度控制方法
CN116501071B (zh) * 2023-06-30 2023-09-05 苏州铼钠克信息技术有限公司 近似回折路径的规划方法和装置、电子设备和存储介质
CN117666475A (zh) * 2023-12-07 2024-03-08 齐齐哈尔大学 一种连续短线段拐角加工路径平滑方法

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4620287A (en) * 1983-01-20 1986-10-28 Dicomed Corporation Method and apparatus for representation of a curve of uniform width
JP2003015714A (ja) * 2001-06-29 2003-01-17 Toshiba Tungaloy Co Ltd 数値制御工作機械
US9551885B2 (en) * 2012-12-21 2017-01-24 Eyeprint Prosthetics Llc Prosthetic lenses and methods of making the same
CN107292938B (zh) * 2016-03-31 2021-06-18 北京神州泰岳软件股份有限公司 一种生成虚圆滑曲线的方法和装置
CN108062073B (zh) * 2016-11-08 2020-06-23 沈阳高精数控智能技术股份有限公司 一种用于高质量加工的圆弧平滑压缩插补方法
CN108132645B (zh) * 2016-12-01 2020-05-19 华中科技大学 一种保证刀具轨迹整体g2连续的曲线拟合方法
CN107255998B (zh) * 2017-07-05 2019-06-25 大连理工大学 误差可控且过给定刀位点的离散加工路径的局部光顺方法
CN108568817B (zh) * 2018-03-13 2021-05-11 同济大学 一种基于贝塞尔曲线的Delta机器人轨迹连接控制方法
CN109799792B (zh) * 2019-01-08 2021-06-01 华侨大学 一种面向高速数控加工微小线段平滑过渡的方法

Also Published As

Publication number Publication date
CN111061213A (zh) 2020-04-24

Similar Documents

Publication Publication Date Title
CN111061213B (zh) 一种基于Bezier曲线转角平滑过渡算法的加工方法
CN109571473B (zh) 一种误差可控的小线段轨迹光顺方法
CN101539769B (zh) 基于二次b样条曲线对g01代码的拟合及插补方法
CN108062073B (zh) 一种用于高质量加工的圆弧平滑压缩插补方法
CN103064344B (zh) 一种基于nurbs曲线插补的速度平滑控制方法
CN106054817B (zh) 局部最优小线段路径解析光顺的实时前瞻插补方法及系统
CN108073138B (zh) 适用于高速高精加工的椭圆弧平滑压缩插补算法
CN100517143C (zh) 数控机床系统中实现对刀路进行插补控制的方法
CN106774144B (zh) 一种基于工业机器人的智能cnc处理方法
CN106647623B (zh) 一种几何精度及衔接速度最优化的五轴联动平滑插补方法
Bi et al. A general, fast and robust B-spline fitting scheme for micro-line tool path under chord error constraint
CN101907876B (zh) 适用于数控装置的指令点整形压缩插补方法
CN108319224B (zh) 一种基于径向曲线插值的多轴数控加工螺旋路径生成方法
CN108568817A (zh) 一种基于贝塞尔曲线的Delta机器人轨迹连接控制方法
CN111966047B (zh) 基于三角函数加减速控制的三轴微线段直接速度过渡方法
CN109976262A (zh) 一种针对微线段加工的全局曲率连续光顺方法
CN108227630A (zh) 一种采用时间参数多项式插补的自由曲面数控加工方法
CN111506988B (zh) 微小线段刀具路径c2连续压缩平滑方法
CN109270892B (zh) 一种数控加工中非圆曲线的最小二乘螺旋线逼近算法
CN115179306A (zh) 一种复杂木模工业机器人铣削及控制方法
CN113359607A (zh) 一种应用于五轴数控机的拐角过渡的轨迹确定方法
CN111610751B (zh) 过点集nurbs插值曲线的插值误差多次细分迭代计算方法
Wu et al. Fast NURBS interpolation based on the biarc guide curve
CN106444636A (zh) 一种适用于连续混合曲线的速度规划参数选择方法
CN114019911B (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