CN109521731A - 一种基于公差带的G2连续Bézier刀具轨迹平滑算法 - Google Patents

一种基于公差带的G2连续Bézier刀具轨迹平滑算法 Download PDF

Info

Publication number
CN109521731A
CN109521731A CN201710845700.4A CN201710845700A CN109521731A CN 109521731 A CN109521731 A CN 109521731A CN 201710845700 A CN201710845700 A CN 201710845700A CN 109521731 A CN109521731 A CN 109521731A
Authority
CN
China
Prior art keywords
point
curve
continuous
weight
tolerance range
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
CN201710845700.4A
Other languages
English (en)
Other versions
CN109521731B (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.)
Shenyang Zhongke Cnc Technology Co ltd
Original Assignee
Shenyang Gaojing Numerical Control Intelligent Technology Co Ltd
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 Shenyang Gaojing Numerical Control Intelligent Technology Co Ltd filed Critical Shenyang Gaojing Numerical Control Intelligent Technology Co Ltd
Priority to CN201710845700.4A priority Critical patent/CN109521731B/zh
Publication of CN109521731A publication Critical patent/CN109521731A/zh
Application granted granted Critical
Publication of CN109521731B publication Critical patent/CN109521731B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/41Numerical 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 interpolation, e.g. the computation of intermediate points between programmed end points to define the path to be followed and the rate of travel along that path
    • G05B19/4103Digital interpolation
    • 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/4097Numerical 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 using design data to control NC machines, e.g. CAD/CAM
    • G05B19/4099Surface or curve machining, making 3D objects, e.g. desktop manufacturing
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00Program-control systems
    • G05B2219/30Nc systems
    • G05B2219/33Director till display
    • G05B2219/33003Algorithm, hashing algorithm
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00Program-control systems
    • G05B2219/30Nc systems
    • G05B2219/34Director, elements to supervisory
    • G05B2219/34085Software interpolator

Landscapes

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

Abstract

本发明涉及一种基于公差带的G2连续Bézier刀具轨迹平滑算法,包括以下步骤:识别连续加工区域;在连续加工区域中,对离散的指令点进行局部插值,将折线加工路径转化为G1连续的二次有理Bézier曲线;调整相邻两条二次有理Bézier曲线的权值和连接点处的切线方向,使得加工路径达到G2连续性;通过建立公差带,将不满足精度要求的二次有理Bézier曲线进行重构。本发明方法实时性好,生成的加工轨迹满足G2连续性和精度要求,能够有效地提高加工效率。

Description

一种基于公差带的G2连续Bézier刀具轨迹平滑算法
技术领域
本发明涉及数控系统中参数样条曲线的拟合及拟合精度的控制,属于数控加工技术领域。
背景技术
在复杂工件的数控加工中,CAM系统通常将CAD系统设计生成的复杂图形转化为由大量微小线段构成的数控加工路径。如果数控系统按照生成的加工路径直接进行插补,由于频繁的加/减速控制,会引起机床振动,降低加工效率和加工质量。一些高档数控系统如FANUC和SIEMENS提出了纳米平滑与压缩器技术,其目的是对离散的小线段进行平滑处理以形成光滑的样条曲线,从而实现高速高精加工,但这些技术都作为商业机密并未公开。因此,如何在满足轮廓误差和机床动态性能等的前提下,将折线加工路径转化为光滑的样条曲线,是数控系统运动轨迹平滑技术的重要研究内容。
近年来,国内外学者针对轨迹平滑技术做了大量研究,如利用C样条拟合方法,将折线加工路径转化为C1连续的样条曲线,其具有连续的切线方向,能够保证机床各驱动轴速度的连续,但由于加速度和加加速度存在频繁突变,会对机床产生冲击,引起振动。因此,提出了一种使用三次Bézier曲线,以曲线张力能量最小为优化函数,曲线间一阶、二阶连续为限制条件的全局插值方法,生成了G2连续加工路径。算法所构造的G2/C2样条曲线,不仅具有连续的切线方向,而且具有连续的曲率,能够保证机床各驱动轴速度和加速度的连续,有效地抑制了机床的振动。为了进一步提高样条曲线的光滑性,提出了一种利用五次Bézier样条、九次和七次校正多项式,生成具有C3连续性加工轨迹的算法。对于上述所有算法,无论是插值方式,还是拟合方式,都只能保证离散指令点到样条曲线的最短距离满足轮廓误差,却不能保证由相邻离散指令点构成的小线段到样条曲线的最远距离不超过设定值。同时,因为拟合与全局插值方式需要进行迭代处理或求解复杂非线性方程组,计算量会随离散指令点的增多而增加,实时性和局部性较差。
发明内容
为了提高加工路径的光滑性和精度,在分析现有样条拟合方法不足的基础上,本发明的目的是提供一种刀具轨迹平滑算法,通过对连续微小线段加工区域中离散指令点进行局部插值、样条曲线调整和样条曲线重构,生成满足G2连续性和精度要求的加工轨迹,有效地提高加工效率。
本发明解决其技术问题所采用的技术方案是:一种基于公差带的G2连续Bézier刀具轨迹平滑算法,包括以下步骤:
识别连续加工区域;
在连续加工区域中,对离散的指令点进行局部插值,将折线加工路径转化为G1连续的二次有理Bézier曲线;
调整相邻两条二次有理Bézier曲线的权值和连接点处的切线方向,使得加工路径达到G2连续性;
通过建立公差带,将不满足精度要求的二次有理Bézier曲线进行重构。
所述二次有理Bézier曲线的调整包括以下步骤:
假设Ci(u)和Ci+1(u)是相邻的两条二次有理Bézier曲线,Pi-1、Qi、Pi、Qi+1、Pi+1为控制点,wi、wi+1为权值,且在连接点Pi处,具有相同的切矢,向量之间的夹角为α,向量之间的夹角为β。
2-1)为了构造出G2连续的加工轨迹,需要将直线QiQi+1绕点Pi进行旋转,使得到的直线与原直线QiQi+1的夹角为θ,与直线Pi-1Qi、Pi+1Qi+1的交点分别为点新的控制点与原曲线控制点之间存在以下关系,
将上式化简后,参数λ、μ可以由下式确定,
相邻样条曲线连接点处的切线方向经过改变后,由新的控制点和旧的权值wi、wi+1构成了两条新的二次有理Bézier曲线其起点和终点处新的曲率值由下式确定,
其中,0表示起点,1表示终点;ki(0)、ki(1)与ki+1(0)、ki+1(1)分别表示样条曲线Ci(u)和Ci+1(u)两端点处的曲率;分别表示样条曲线两端点处的曲率;
2-2)为了保证加工轨迹的G2连续性,需要调整曲线的权值wi,令其起点处的曲率值等于前一Ci-1(u)曲线终点处的曲率值,可得下式,
将上式化简后,可以求得调整后曲线的权值如下式所示,
调整后曲线终点处的曲率值kl由下式得到,
2-3)为了保证加工轨迹的G2连续性,需要调整曲线的权值wi+1,令其终点处的曲率值等于ki+1(1),可得下式,
将上式化简后,可以求得调整后曲线的权值如下式所示,
调整后曲线起点处的曲率值kr由下式得到,
2-4)为了使相邻两样条曲线在连接点Pi处具有相同的曲率值,令kl=kr,可得下式,
将上式化简后,旋转的角度θ可以由下式得到,
θ=arctan(-B/A)
其中,
A=cosα+s2cosβ+s0+s1s2
B=sinα-s2sinβ
然后,由2-1),得到参数λ和μ。最后,由2-1)、2-2)、2-3),求得调整后二次有理Bézier曲线的控制点和权值从而构造出G2连续的加工轨迹。
所述通过建立公差带,将不满足精度要求的二次有理Bézier曲线进行重构包括以下步骤:
3-1)识别不满足精度要求的曲线。对于标准型二次有理Bézier曲线,肩点S处的切矢平行于直线P0P2,因此,曲线C(u)上到线段P0P2距离最大的点为S=C(1/2)。设二次有理Bézier曲线Ci(u)的肩点S到线段Pi-1Pi的距离为d,系统设置的最大轮廓误差为δ,如果d>δ,那么曲线Ci(u)不满足精度要求;否则,为满足精度要求;
3-2)基于公差带的曲线重构。
所述基于公差带的曲线重构包括以下步骤:
4-1)对于不满足精度要求的样条曲线Ci(u),建立直线Pi-1Pi的平行线令它们之间的距离为δ,直线与线段Pi-1Qi、PiQi相交于点此时,直线Pi-1Pi和它们之间的部分便构成了误差为δ的公差带;
4-2)为了使得重构后的曲线既满足精度要求,又满足G2连续性,需要在数据点Pi-1、Pi之间构造两条样条曲线,并使得它们在连接点处具有相同的曲率值。
假设第一条重构曲线Cil(u)的控制点为Pi-1、Qil权值为wl;第二条重构曲线Cir(u)的控制点为Pi、Qir权值为wr;点Qil、Qir分别在线段Pi-1Qi、PiQi、QilQir上,且它们满足如下关系,
4-3)样条曲线Cil(u)在点Pi-1处G2连续。令ki(0)=kil(0),可得下式,
根据上式和4-2)中的关系式,可以得到系数k0,如下式所示,
4-4)样条曲线Cir(u)在点Pi处G2连续。令ki(1)=kir(1),可得下式,
根据上式和4-2)中的关系式,可以得到系数k2,如下式所示,
4-5)Cil(u)和Cir(u)在点处G2连续。令kil(1)=kir(0),可得下式,
根据上式和4-2)中的关系式,可以得到系数k1,如下式所示,
4-6)令点Qil、Qir分别与点重合,点位于线段上,由此,可以得到下式,
由k0、k2和4-5)可以求得k1,再由4-2),便能够确定控制点Qil、Qir的位置。最后,根据4-3)、4-4),可以得到两条重构曲线的权值wl、wr。此时,便完成了样条曲线的重构,重构后的曲线被严格限制在公差带内,且满足G2连续性。
本发明具有以下有益效果及优点:
1.本发明方法能够构造出G2连续的二次有理Bézier曲线,其具有连续的切线方向和曲率,能够有效地提高数控机床的加工效率。
2.本发明方法采用建立公差带的方式,对不满足精度要求的二次有理Bézier曲线进行重构,将其严格限制在公差带内,不仅能够保证离散指令点到样条曲线的最短距离满足轮廓误差,也能够保证由相邻离散指令点构成的小线段到样条曲线的最远距离满足轮廓误差。
3.本发明方法采用局部插值和调整的方式进行轨迹平滑,具有较好的局部性,即某段样条曲线的调整不会影响加工轨迹的整体形状。又因为其不涉及迭代处理和复杂数值计算,算法效率高。因此,本发明方法具有较好的实时性,能够实现在线轨迹平滑处理。
附图说明
图1为本发明方法的流程图;
图2为连续微小线段加工区域识别的示意图;
图3为二次有理Bézier曲线局部插值的示意图;
图4为二次有理Bézier曲线调整的示意图;
图5为二次有理Bézier曲线重构的示意图;
具体实施方式
下面结合实施案例对本发明做进一步的详细说明。
本发明是一种基于公差带的G2连续Bézier刀具轨迹平滑算法,根据双弓高误差限制,从由离散小线段构成的加工路径中识别出连续微小线段加工区域。对于连续微小线段加工区域,首先,对离散的指令点进行局部插值,将折线加工路径转化为G1连续的二次有理Bézier曲线;然后,调整相邻两条二次有理Bézier曲线的权值和连接点处的切线方向,使得加工路径达到G2连续性;最后,通过建立公差带,将不满足精度要求的二次有理Bézier曲线进行重构。
本发明提供了一种基于公差带的G2连续Bézier刀具轨迹平滑算法,包括以下步骤:
1.加工区域的识别,根据双弓高误差判断条件,从由离散小线段构成的加工路径中识别出连续微小线段加工区域。
2.二次有理Bézier曲线的局部插值,对连续微小线段加工区域中的离散指令点进行局部插值,将折线加工路径转化为G1连续的二次有理Bézier曲线。
3.二次有理Bézier曲线的调整,调整相邻两条二次有理Bézier曲线的权值和连接点处的切线方向,使得加工路径达到G2连续性。
4.二次有理Bézier曲线的重构,根据离散数据点建立公差带,在满足G2连续的条件下,对不满足轮廓误差要求的二次有理Bézier曲线进行重构,将其严格限制在公差带内。
如图1所示,本发明提出了一种基于公差带的G2连续Bézier刀具轨迹平滑算法,解决了加工路径光滑性和精度的问题,方法由加工区域的识别、二次有理Bézier曲线的局部插值、二次有理Bézier曲线的调整和二次有理Bézier曲线的重构4个部分组成,提高了加工的效率。
根据双弓高误差判断条件,从由离散小线段构成的加工路径中识别出连续微小线段加工区域。
对连续区内的离散数据点进行局部插值,将由指令点构成的折线加工路径转化为G1连续的二次有理Bézier曲线加工路径。
通过调整相邻两条二次有理Bézier曲线的权值及其连接点处的切线方向,使得加工路径达到G2连续。
根据离散数据点建立公差带,在满足G2连续的条件下,对不满足轮廓误差要求的二次有理Bézier曲线进行重构,将其严格限制在公差带内。
本发明的具体步骤如下:
1.连续微小线段加工区域的识别
如图2所示,Pi-1、Pi和Pi+1为相邻的三个数据点,l1和l2为小线段的段长,θ为小线段间的拐角,双弓高误差由下式确定,
δ1=R(1-cosφ1)
δ2=R(1-cosφ2)=R(1-cos(π-θ-φ1))
其中,δ1、δ2分别为小线段Pi-1Pi和PiPi+1的弓高误差,φ1为OPi-1和OPi夹角的一半。φ2为OPi和OPi+1夹角的一半。R表示曲率半径。
如果δ1或δ2大于设定的最大弓高误差值δmax,则Pi为断点。那么,两个断点连同它们之间的数据点就构成了一个连续微小线段加工区域。
2.二次有理Bézier曲线的局部插值
标准型二次有理Bézier曲线的参数形式如下式所示,
其中,R0、R1、R2为控制点,w1为权值,u为参变量,且u∈[0,1]。
由图3可知,当已知控制点R0、R2,以及两点处的切矢T0、T2时,可以得到控制点R1,再给定一点S就能求得权值w1,从而唯一确定该曲线。
令点M为线段R0R2的中点,点S0为线段R1M和∠R1R0R2的角平分线的交点,点S2为R1M和∠R1R2R0的角平分线的交点,则点S=(S0+S2)/2。
当u=1/2时,点S=C(1/2)被称为肩点,将u=1/2代入曲线参数公式,可得下式,
设s是点R1和点M之间线性插值的参数,那么,存在某个s使得下式成立,
S=(1-s)M+sR1
根据上述两个S的关系式,权值w1可以由下式得到,
由此,便确定了该二次有理Bézier曲线。
3.二次有理Bézier曲线的调整
如图4所示,假设Ci(u)和Ci+1(u)是相邻的两条二次有理Bézier曲线,Pi-1、Qi、Pi、Qi+1、Pi+1为控制点,wi、wi+1为权值,且在连接点Pi处,具有相同的切矢,向量之间的夹角为α,向量之间的夹角为β。
(1)为了构造出G2连续的加工轨迹,需要将直线QiQi+1绕点Pi进行旋转,使得到的直线与原直线QiQi+1的夹角为θ,与直线Pi-1Qi、Pi+1Qi+1的交点分别为点新的控制点与原曲线控制点之间存在以下关系,
将上式化简后,参数λ、μ可以由下式确定,
相邻样条曲线连接点处的切线方向经过改变后,由新的控制点和旧的权值wi、wi+1构成了两条新的二次有理Bézier曲线其起点和终点处新的曲率值由下式确定,
其中,0表示起点,1表示终点;ki(0)、ki(1)与ki+1(0)、ki+1(1)分别表示样条曲线Ci(u)和Ci+1(u)两端点处的曲率;分别表示样条曲线两端点处的曲率;
(2)为了保证加工轨迹的G2连续性,需要调整曲线的权值wi,令其起点处的曲率值等于前一Ci-1(u)曲线终点处的曲率值,可得下式,
将上式化简后,可以求得调整后曲线的权值如下式所示,
调整后曲线终点处的曲率值kl由下式得到,
(3)为了保证加工轨迹的G2连续性,需要调整曲线的权值wi+1,令其终点处的曲率值等于ki+1(1),可得下式,
将上式化简后,可以求得调整后曲线的权值如下式所示,
调整后曲线起点处的曲率值kr由下式得到,
(4)为了使相邻两样条曲线在连接点Pi处具有相同的曲率值,令kl=kr,可得下式,
其中,ki-1(1)表示样条曲线Ci-1(u)的终点曲率值;
将上式化简后,旋转的角度θ可以由下式得到,
θ=arctan(-B/A)
其中,
A=cosα+s2cosβ+s0+s1s2
B=sinα-s2sinβ
然后,由(1),得到参数λ和μ。最后,由(1)、(2)、(3),求得调整后二次有理Bézier曲线的控制点和权值,从而构造出G2连续的加工轨迹。
4.二次有理Bézier曲线的重构
(1)识别不满足精度要求的曲线。如图3所示,对于标准型二次有理Bézier曲线,肩点S处的切矢平行于直线P0P2,因此,曲线C(u)上到线段P0P2距离最大的点为S=C(1/2)。设二次有理Bézier曲线Ci(u)的肩点S到线段Pi-1Pi的距离为d,系统设置的最大轮廓误差为δ,如果d>δ,那么曲线Ci(u)不满足精度要求;
(2)基于公差带的曲线重构。
1)如图5所示,对于不满足精度要求的样条曲线Ci(u),建立直线Pi-1Pi的平行线令它们之间的距离为δ,直线与线段Pi-1Qi、PiQi相交于点此时,直线Pi- 1Pi和它们之间的部分便构成了误差为δ的公差带;
2)为了使得重构后的曲线既满足精度要求,又满足G2连续性,需要在数据点Pi-1、Pi之间构造两条样条曲线,并使得它们在连接点处具有相同的曲率值。
假设第一条重构曲线Cil(u)的控制点为Pi-1、Qil权值为wl;第二条重构曲线Cir(u)的控制点为Pi、Qir权值为wr;点Qil、Qir分别在线段Pi-1Qi、PiQi、QilQir上,且它们满足如下关系,
3)样条曲线Cil(u)在点Pi-1处G2连续。令ki(0)=kil(0),可得下式,
根据上式和2)中的关系式,可以得到系数k0,如下式所示,
4)样条曲线Cir(u)在点Pi处G2连续。令ki(1)=kir(1),可得下式,
根据上式和2)中的关系式,可以得到系数k2,如下式所示,
5)Cil(u)和Cir(u)在点处G2连续。令kil(1)=kir(0),可得下式,
根据上式和2)中的关系式,可以得到系数k1,如下式所示,
6)令点Qil、Qir分别与点重合,点位于线段上,由此,可以得到下式,
由k0、k2和5)可以求得k1,再由2),便能够确定控制点Qil、Qir的位置。最后,根据3)、4),可以得到两条重构曲线的权值wl、wr。此时,便完成了样条曲线的重构,重构后的曲线被严格限制在公差带内,且满足G2连续性。

Claims (4)

1.一种基于公差带的G2连续Bézier刀具轨迹平滑算法,其特征在于包括以下步骤:
识别连续加工区域;
在连续加工区域中,对离散的指令点进行局部插值,将折线加工路径转化为G1连续的二次有理Bézier曲线;
调整相邻两条二次有理Bézier曲线的权值和连接点处的切线方向,使加工路径达到G2连续性;
通过建立公差带,将不满足精度要求的二次有理Bézier曲线进行重构。
2.根据权利要求1所述的一种基于公差带的G2连续Bézier刀具轨迹平滑算法,其特征在于所述调整相邻两条二次有理Bézier曲线的权值和连接点处的切线方向包括以下步骤:
假设Ci(u)和Ci+1(u)是相邻的两条二次有理Bézier曲线,Pi-1、Qi、Pi、Qi+1、Pi+1为控制点,wi、wi+1为权值,且在连接点Pi处,具有相同的切矢,向量之间的夹角为α,向量之间的夹角为β;
2-1)将直线QiQi+1绕点Pi进行旋转,使得到的直线与原直线QiQi+1的夹角为θ,与直线Pi- 1Qi、Pi+1Qi+1的交点分别为点新的控制点与原曲线控制点之间存在以下关系,
将上式化简后,参数λ、μ由下式确定,
相邻二次有理Bézier曲线连接点处的切线方向经过改变后,由新的控制点和旧的权值wi、wi+1构成了两条新的二次有理Bézier曲线其起点和终点处新的曲率值由下式确定,
其中,0表示起点,1表示终点;ki(0)、ki(1)与ki+1(0)、ki+1(1)分别表示样条曲线Ci(u)和Ci+1(u)两端点处的曲率;分别表示样条曲线两端点处的曲率;
2-2)调整曲线的权值wi,令其起点处的曲率值等于前一Ci-1(u)曲线终点处的曲率值,得下式,
将上式化简后,求得调整后曲线的权值如下式所示,
调整后曲线终点处的曲率值kl由下式得到,
2-3)调整曲线的权值wi+1,令其终点处的曲率值等于ki+1(1),可得下式,
将上式化简后,可以求得调整后曲线的权值如下式所示,
调整后曲线起点处的曲率值kr由下式得到,
2-4)为了使相邻两样条曲线在连接点Pi处具有相同的曲率值,令kl=kr,得下式,
将上式化简后,旋转的角度θ由下式得到,
θ=arctan(-B/A)
其中,
A=cosα+s2cosβ+s0+s1s2
B=sinα-s2sinβ
然后,由2-1),得到参数λ和μ;最后,由2-1)、2-2)、2-3),求得调整后二次有理Bézier曲线的控制点和权值,从而构造出G2连续的加工轨迹。
3.根据权利要求1所述的一种基于公差带的G2连续Bézier刀具轨迹平滑算法,其特征在于所述通过建立公差带,将不满足精度要求的二次有理Bézier曲线进行重构包括以下步骤:
3-1)识别不满足精度要求的曲线:对于标准型二次有理Bézier曲线,肩点S处的切矢平行于直线P0P2,因此,曲线C(u)上到线段P0P2距离最大的点为S=C(1/2);设二次有理Bézier曲线即样条曲线Ci(u)的肩点S到线段Pi-1Pi的距离为d,数控系统设置的最大轮廓误差为δ,如果d>δ,那么曲线Ci(u)不满足精度要求;
3-2)基于公差带的曲线重构。
4.根据权利要求3所述的一种基于公差带的G2连续Bézier刀具轨迹平滑算法,其特征在于所述基于公差带的曲线重构包括以下步骤:
4-1)对于不满足精度要求的样条曲线Ci(u),建立直线Pi-1Pi的平行线令它们之间的距离为δ,直线与线段Pi-1Qi、PiQi相交于点此时,直线Pi-1Pi和它们之间的部分便构成了误差为δ的公差带;
4-2)为了使得重构后的曲线既满足精度要求,又满足G2连续性,在数据点Pi-1、Pi之间构造两条样条曲线,并使得它们在连接点处具有相同的曲率值;
假设第一条重构曲线Cil(u)的控制点为Pi-1、Qil权值为wl;第二条重构曲线Cir(u)的控制点为Pi、Qir权值为wr;点Qil、Qir分别在线段Pi-1Qi、PiQi、QilQir上,且它们满足如下关系,
4-3)样条曲线Cil(u)在点Pi-1处G2连续;令ki(0)=kil(0),可得下式,
根据上式和4-2)中的关系式,可以得到系数k0,如下式所示,
4-4)样条曲线Cir(u)在点Pi处G2连续;令ki(1)=kir(1),可得下式,
根据上式和4-2)中的关系式,可以得到系数k2,如下式所示,
4-5)Cil(u)和Cir(u)在点P处G2连续;令kil(1)=kir(0),可得下式,
根据上式和4-2)中的关系式,可以得到系数k1,如下式所示,
4-6)令点Qil、Qir分别与点重合,点位于线段上,得到下式,
由k0、k2和4-5)求得k1,再由4-2),确定控制点Qil、Qir的位置;最后,根据4-3)、4-4),得到两条重构曲线的权值wl、wr;此时,便完成了样条曲线的重构,重构后的曲线被限制在公差带内,且满足G2连续性。
CN201710845700.4A 2017-09-19 2017-09-19 一种基于公差带的G2连续Bézier刀具轨迹平滑算法 Active CN109521731B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710845700.4A CN109521731B (zh) 2017-09-19 2017-09-19 一种基于公差带的G2连续Bézier刀具轨迹平滑算法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710845700.4A CN109521731B (zh) 2017-09-19 2017-09-19 一种基于公差带的G2连续Bézier刀具轨迹平滑算法

Publications (2)

Publication Number Publication Date
CN109521731A true CN109521731A (zh) 2019-03-26
CN109521731B CN109521731B (zh) 2021-07-30

Family

ID=65767946

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710845700.4A Active CN109521731B (zh) 2017-09-19 2017-09-19 一种基于公差带的G2连续Bézier刀具轨迹平滑算法

Country Status (1)

Country Link
CN (1) CN109521731B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110865610A (zh) * 2019-10-18 2020-03-06 中国工程物理研究院机械制造工艺研究所 一种基于机床振动抑制的刀具轨迹插值和速度规划方法
CN112506139A (zh) * 2020-11-23 2021-03-16 西北工业大学 一种五轴短直线段轨迹的局部拐角光顺方法
CN113721551A (zh) * 2020-05-26 2021-11-30 沈机(上海)智能系统研发设计有限公司 数控加工方法和数控加工设备

Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1116186A1 (en) * 1998-09-24 2001-07-18 Brigit Ananya Computer curve construction system and method
EP1220069A1 (de) * 2000-12-27 2002-07-03 Siemens Aktiengesellschaft Verfahren zur Erstellung flexibler Kurvenscheibenfunktionen über das Steuerungs- oder Anwenderprogramm
DE69727200T2 (de) * 1996-12-05 2004-10-14 Adobe Systems, Inc., San Jose Verfahren und Program zur Erzeugung und Modifizierung von Kurven auf einen Rechnerbildschirm
CN101539769A (zh) * 2009-04-28 2009-09-23 中国科学院数学与系统科学研究院 基于二次b样条曲线对g01代码的拟合及插补方法
CN101907876A (zh) * 2010-05-28 2010-12-08 沈阳高精数控技术有限公司 适用于数控装置的指令点整形压缩插补方法
CN102147600A (zh) * 2011-04-30 2011-08-10 上海交通大学 实时生成曲率连续路径的数控插补系统
CN103631198A (zh) * 2013-11-04 2014-03-12 上海交通大学 基于G2连续Bézier曲线的刀具轨迹压缩方法
CN103646150A (zh) * 2013-12-23 2014-03-19 北京航空航天大学 一种基于5次Bezier曲线的G3连续过渡曲线构造方法
CN105955194A (zh) * 2016-05-10 2016-09-21 大连理工大学 一种离散加工路径的拐点平滑方法
CN106663432A (zh) * 2014-07-02 2017-05-10 杜比国际公司 对压缩的hoa表示解码的方法和装置以及对压缩的hoa表示编码的方法和装置
CN106651987A (zh) * 2016-10-10 2017-05-10 腾讯科技(深圳)有限公司 路径规划方法及装置
CN106774145A (zh) * 2016-12-23 2017-05-31 华中科技大学 一种用于生成无干涉的五轴加工轨迹的投影算法
CN106802631A (zh) * 2017-03-13 2017-06-06 北京航空航天大学 曲面加工的方法与装置
CN106826829A (zh) * 2017-02-22 2017-06-13 武汉工程大学 一种可控误差的工业机器人光顺运动轨迹生成方法

Patent Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE69727200T2 (de) * 1996-12-05 2004-10-14 Adobe Systems, Inc., San Jose Verfahren und Program zur Erzeugung und Modifizierung von Kurven auf einen Rechnerbildschirm
EP1116186A1 (en) * 1998-09-24 2001-07-18 Brigit Ananya Computer curve construction system and method
EP1220069A1 (de) * 2000-12-27 2002-07-03 Siemens Aktiengesellschaft Verfahren zur Erstellung flexibler Kurvenscheibenfunktionen über das Steuerungs- oder Anwenderprogramm
CN101539769A (zh) * 2009-04-28 2009-09-23 中国科学院数学与系统科学研究院 基于二次b样条曲线对g01代码的拟合及插补方法
CN101907876A (zh) * 2010-05-28 2010-12-08 沈阳高精数控技术有限公司 适用于数控装置的指令点整形压缩插补方法
CN102147600A (zh) * 2011-04-30 2011-08-10 上海交通大学 实时生成曲率连续路径的数控插补系统
CN103631198A (zh) * 2013-11-04 2014-03-12 上海交通大学 基于G2连续Bézier曲线的刀具轨迹压缩方法
CN103646150A (zh) * 2013-12-23 2014-03-19 北京航空航天大学 一种基于5次Bezier曲线的G3连续过渡曲线构造方法
CN106663432A (zh) * 2014-07-02 2017-05-10 杜比国际公司 对压缩的hoa表示解码的方法和装置以及对压缩的hoa表示编码的方法和装置
CN105955194A (zh) * 2016-05-10 2016-09-21 大连理工大学 一种离散加工路径的拐点平滑方法
CN106651987A (zh) * 2016-10-10 2017-05-10 腾讯科技(深圳)有限公司 路径规划方法及装置
CN106774145A (zh) * 2016-12-23 2017-05-31 华中科技大学 一种用于生成无干涉的五轴加工轨迹的投影算法
CN106826829A (zh) * 2017-02-22 2017-06-13 武汉工程大学 一种可控误差的工业机器人光顺运动轨迹生成方法
CN106802631A (zh) * 2017-03-13 2017-06-06 北京航空航天大学 曲面加工的方法与装置

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110865610A (zh) * 2019-10-18 2020-03-06 中国工程物理研究院机械制造工艺研究所 一种基于机床振动抑制的刀具轨迹插值和速度规划方法
CN113721551A (zh) * 2020-05-26 2021-11-30 沈机(上海)智能系统研发设计有限公司 数控加工方法和数控加工设备
CN112506139A (zh) * 2020-11-23 2021-03-16 西北工业大学 一种五轴短直线段轨迹的局部拐角光顺方法
CN112506139B (zh) * 2020-11-23 2022-02-25 西北工业大学 一种五轴短直线段轨迹的局部拐角光顺方法

Also Published As

Publication number Publication date
CN109521731B (zh) 2021-07-30

Similar Documents

Publication Publication Date Title
CN109521731A (zh) 一种基于公差带的G2连续Bézier刀具轨迹平滑算法
WO2022027836A1 (zh) 针对五轴加工刀尖点与刀轴方向非线性误差的控制方法
CN104102171A (zh) 一种球头刀多轴加工刀轴矢量优化方法
CN105955194B (zh) 一种离散加工路径的拐点平滑方法
CN103631198B (zh) 基于G2连续Bézier曲线的刀具轨迹压缩方法
CN103433804B (zh) 基于五轴无干涉刀轴控制线的叶轮加工刀轴矢量控制方法
CN104635619A (zh) 基于刀具矢量插补的五轴数控加工方法
CN108062073A (zh) 一种用于高质量加工的圆弧平滑压缩插补算法
CN106736846B (zh) 一种车铣复合机床车刀刀具半径补偿方法
CN110488747A (zh) 基于机床旋转轴角加速度最小的刀矢全局优化方法
CN103537743A (zh) 多轴曲面数控加工复杂曲面零件的方法
CN108073138B (zh) 适用于高速高精加工的椭圆弧平滑压缩插补算法
CN112162527A (zh) 适用于五轴数控装置的刀具路径拐角平滑过渡方法
CN112525133B (zh) 法向公差向极半径公差的转换方法
CN113146631B (zh) 一种异型预制体机器人针刺成形路径规划方法
CN112847364A (zh) 机器人铣削加工刀具路径c3连续的在线实时解析光顺方法
CN104238456A (zh) 一种非球头刀铣削加工自由曲面的方法
CN110399681A (zh) 一种圆弧头立铣刀刀刃曲线的参数化建模方法
CN106094737B (zh) 一种指定加工误差条件下的数控加工速度优化控制方法
CN110716493A (zh) 一种五轴微线段加工路径光顺方法
CN116400646A (zh) 一种基于奇异补偿的五轴数控机床通用后置处理方法
CN107085412B (zh) 一种控制曲线车削加工切削力突变的方法
CN102354151B (zh) 一种应用于多层鞋革数控裁割机床的切向跟随插补方法
CN107798730A (zh) 一种结构网格附面层网格自动生成方法
WO2024124794A1 (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
CP01 Change in the name or title of a patent holder
CP01 Change in the name or title of a patent holder

Address after: No. 16-2, Nanping East Road, Dongling District, Shenyang City, Liaoning Province, 110168

Patentee after: Shenyang Zhongke CNC Technology Co.,Ltd.

Address before: No. 16-2, Nanping East Road, Dongling District, Shenyang City, Liaoning Province, 110168

Patentee before: Shenyang Golding Nc Intelligence Tech.co.,ltd.