CN101556472B - 数控系统基于多周期最优拐角的小直线段插补方法 - Google Patents

数控系统基于多周期最优拐角的小直线段插补方法 Download PDF

Info

Publication number
CN101556472B
CN101556472B CN2009100839504A CN200910083950A CN101556472B CN 101556472 B CN101556472 B CN 101556472B CN 2009100839504 A CN2009100839504 A CN 2009100839504A CN 200910083950 A CN200910083950 A CN 200910083950A CN 101556472 B CN101556472 B CN 101556472B
Authority
CN
China
Prior art keywords
turning
interpolation
line segment
straight
speed
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.)
Expired - Fee Related
Application number
CN2009100839504A
Other languages
English (en)
Other versions
CN101556472A (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.)
Academy of Mathematics and Systems Science of CAS
Original Assignee
Academy of Mathematics and Systems Science of CAS
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 Academy of Mathematics and Systems Science of CAS filed Critical Academy of Mathematics and Systems Science of CAS
Priority to CN2009100839504A priority Critical patent/CN101556472B/zh
Publication of CN101556472A publication Critical patent/CN101556472A/zh
Application granted granted Critical
Publication of CN101556472B publication Critical patent/CN101556472B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Numerical Control (AREA)

Abstract

一种数控系统基于多周期最优拐角的小直线段插补方法,步骤如下:先基于加工精度和最大加工速度的限制,根据加工路径上拐角处的几何参数、机床各驱动轴的最大加速度和优化目标,确定加工路径上的每个拐角多周期过渡的最优拐角插补参数。再对每个拐角的最优拐角插补参数进行调整,使每个小直线段的两端加工速度满足可达性要求。然后对加工路径的各小直线段进行直线段和拐角的插补,并对其直线段插补点序列和拐角处的插补点序列进行衔接,实时顺序输出插补点序列,驱动数控机床进行实际加工。本发明方法能有效提高整体加工速度,计算速度快,满足实时加工需求,工作性能稳定可靠,实用性强,能适用于多种不同型号的三轴数控机床。

Description

数控系统基于多周期最优拐角的小直线段插补方法
技术领域
本发明涉及一种数控系统基于多周期最优拐角的小直线段插补方法,属于数控系统的数字控制加工技术领域。
背景技术
数控机床在加工复杂曲面时,通常是在设定的加工精度范围内,把被加工曲面离散为大量的空间小直线段。在对小直线段进行加工过程中,数控机床的加工速度与加工精度互相制约。因为实际加工时,数控机床的加工方向要频繁改变,造成机床也频繁起停,限制了加工速度,同时容易产生振动,从而影响到被加工曲面的加工精度和表面光洁度。尤其是在拐角处,由于速度方向的突然变化,如果高速通过拐角,可能超出机床各驱动轴的最大加速度,或者超出机床驱动电机所能承受的最大扭矩而产生振动,缩短机床寿命和降低被加工产品的质量。而降低拐角处的加工速度,因加工轨迹被离散成大量的小直线段,使得整体加工速度受到极大限制而延长加工时间,降低生产效率。
因此,如何在满足产品的加工精度和机床各驱动轴最大加速度的约束规定下,尽可能地提高加工速度是数控机床实现高速和高精度加工的奋斗目标。
目前,数控机床对小直线段进行插补加工的方法主要有以下几种:
参见图1,介绍传统的小直线段插补方法的速度-时间示意图,其中Vm是最大加工速度,其特点是每个待插补小直线段的起点和终点的速度均是零。如果采用该插补控制方法会造成机床的频繁起停,严重限制加工速度的提高。
参见图2,介绍二维小直线段拐角处等速率过渡插补方法的速度-时间示意图,其特点是待插补小线段的起点和终点的速度不为零,而是根据机床各个驱动轴的加速度限制确定其拐角的过渡速度,且拐角前速度和拐角后速度大小相等,在拐角处采用一个插补周期进行过渡。该方法虽然避免了数控机床的频繁起停,但是,由该方法得到的拐角速度一般很小,且因限制拐角前速度和拐角后速度的大小相等,没有充分利用机床各驱动轴的加速能力。因此其整体加工速度还有提高的余地。
参见图3,介绍采用拐角圆弧过渡方法示意图,图中的虚线是原始加工路径,实线是在拐角处插入圆弧的过渡方法的加工路径,相当于把一个拐角处速度方向的改变分散到多个插补周期中进行过渡,以提高拐角加工速度。然而,该方法直接确定了刀具在拐角处的加工轨迹,同样没有充分利用数控机床各个驱动轴的最大加速能力。
发明内容
有鉴于此,本发明的目的是提供一种数控系统基于多周期最优拐角的小直线段插补方法,作为一种新的针对三维空间小直线段的插补方法,本发明充分利用数控机床各驱动轴的最大加速能力,在拐角处速度方向的改变采用多个插补周期进行过渡,同时采用前瞻处理算法,在保证加工精度的前提下,实现了拐角处的最优插补,提高了整体的加工速度,且计算速度快,满足实时加工要求,工作性能稳定、可靠。再者,本发明方法实用性很强,能够适用于多种不同型号的三轴数控机床。
为了达到上述目的,本发明提供了一种数控系统基于多周期最优拐角的小直线段插补方法,其特征在于,所述方法包括下列操作步骤:
(1)根据数控机床各驱动轴的最大加速度、加工路径拐角处的几何参数、“拐角前速度和拐角后速度之和为最大”的优化目标和每个拐角的过渡周期数的初值n,确定加工路径上的每个拐角在采用多个插补周期来改变加工速度方向的拐角总加速度;再根据该拐角的拐角总加速度分别计算其拐角速度和拐角距离,其中拐角速度包括拐角前速度和拐角后速度,拐角距离包括拐角前距离和拐角后距离,并将该拐角的拐角过渡周期数、拐角速度和拐角距离以及拐角总加速度,作为该拐角在其加工路径上采用多个插补周期来改变加工速度方向的过渡方式的拐角处的最优插补参数;
(2)进行前瞻处理:对每个拐角在所述拐角处的的最优插补参数进行调整,使每个小直线段的两端加工速度满足可达性要求,即在该小直线段实际插补长度范围内,其起点速度能以该直线段插补时的最大加速度通过加速或减速方式,达到终点速度;也就是完成下述小直线段的前瞻处理操作内容:
(21)根据系统配置,顺序将加工路径上每m个连续小直线段作为一组进行前瞻处理,该m的数值大小应满足数控机床的实时加工要求,且规定该m个待插补小直线段中的第一个小直线段的起始速度和第m个小直线段的终止速度都为零,且其对应拐角的其他最优插补参数也都为零;
(22)顺次对一组待前瞻处理的小直线段的每个小直线段的两端速度进行可达性判断:即在该小直线段实际插补长度范围内,其起点速度能够以该直线段插补时的最大加速度通过加速或减速方式,达到终点速度;直到完成该m个小直线段的前瞻处理;
(23)继续顺序对后面的每组需要进行前瞻处理的直线段执行步骤(22)的操作;直到加工路径上的所有小直线段都完成前瞻处理后,保存前瞻处理后每个拐角的最优插补参数,以供后续步骤使用。
(3)根据前瞻处理后的最优插补参数、被加工小直线段的长度和机床各驱动轴的最大加速度限制,对加工路径的各个小直线段分别进行直线段的插补和拐角的插补,并保证每个小直线段的直线段插补点序列和该直线段拐角的插补点序列准确衔接,实时顺序输出插补点序列,驱动数控机床进行实际加工;即对所有的待插补的小直线段按照小直线段的输入顺序依次执行下述操作内容:
(31)顺序读取一个待插补小直线段及其拐角的前瞻处理后的最优插补参数,先进行直线段上插补,确定该直线段上插补的插补点序列:直线段上插补的插补点序列的起点是该直线段起点拐角处插补的最后一个插补点;根据直线段两端的速度,即前述步骤已经调整好的该小直线段的起点拐角处的拐角后速度和终点拐角处的拐角前速度,先计算在该直线段上插补时以直线加减速方式运动时所能达到的最大速度Vmm,该Vmm数值应不大于最大加工速度;再计算该直线段上的包括加速时间、匀速时间和减速时间的各个插补时间,最后,根据该直线段的起始速度、加速方式和相应加速方式下的时间进行该直线段的插补,得到该直线段的插补点序列,并实时输出插补点驱动数控机床进行实际加工,直线段插补完之后,接着执行步骤(32)的操作;
(32)在拐角处插补,确定该拐角处插补的插补点序列:先根据该拐角的拐角插补的开始位置、拐角前速度和该拐角每个插补周期的加速度:该拐角的拐角总加速度除以该拐角的过渡周期数的商,确定该拐角处的插补点序列;实际插补时,在该小直线段的直线段上插补的最后一个插补点和其拐角插补的第一个插补点不能在同一点重合时,则将该直线段上插补的最后一个插补点设置为该拐角插补点序列的第一个插补点,并在上述步骤所确定的拐角处插补点序列的基础上,根据残留比重新确定新的拐角处插补点序列,并实时输出插补点驱动数控机床进行实际加工。
所述拐角是加工路径中不在同一直线上的两个相邻小直线段的连接点;所述拐角的多周期过渡是一种把拐角处加工速度方向在一个插补周期的改变分散到多个插补周期进行逐渐改变的加工方式;所述最优插补参数是拐角在多周期过渡时的下述六个加工参数:拐角的过渡周期数、拐角前速度V1、拐角后速度V2、拐角开始距离、拐角结束距离、和拐角总加速度An,其中,拐角的过渡周期数是改变拐角速度方向所用的插补周期个数,其初值n是根据经验选取的一个不小于2的自然数,且后续操作步骤中,该数值会进行相应调整;拐角前速度V1是在拐角处插补开始时的加工速度;拐角后速度V2是在拐角处插补结束时的加工速度;拐角开始距离是拐角插补的开始位置与拐角的距离,拐角结束距离是拐角插补的结束位置与拐角的距离,拐角总加速度An是在一个插补周期内,由拐角前速度过渡到拐角后速度所需的数控机床各驱动轴加速度的矢量和。
所述方法在执行步骤(1)以后,还要继续执行下述操作内容:
(1A)分别判断前述步骤确定的拐角的拐角前速度和拐角后速度是否大于最大加工速度,如果是,则采用降低该拐角的过渡周期数来降低其拐角速度,直到该拐角的拐角前速度和拐角后速度都不大于最大加工速度为止;同时,对该拐角的其它三个最优插补参数:拐角前距离、拐角后距离和拐角总加速度也根据新的过渡周期数做相应调整,接着执行步骤(1B)操作;如果否,则不对该拐角的最优拐角插补参数作调整,直接执行步骤(1B)操作;
(1B)根据步骤(1A)得到的最优插补参数,确定该拐角处插补的插补轨迹,并判断该拐角处的插补因偏离原始加工路径所产生的尺寸误差是否超出加工精度要求;如果是,则降低该拐角的过渡周期数,并对该拐角的其它最优插补参数也根据新的过渡周期数做相应调整,直到满足加工精度要求为止;并将最后得到的最优插补参数进行存储,以备后用,接着执行步骤(1C)操作;如果否,则不对该拐角的最优插补参数作调整,直接执行步骤(1C)操作;
(1C)对每个拐角顺序执行步骤(1)、(1A)、(1B)、(1C)的操作,直到对所有的拐角均确定其最优插补参数。
所述步骤(1)进一步包括下列操作内容:
(11)读取直线插补方式的数控加工代码,先初步确定拐角多周期过渡时该拐角总加速度的取值范围,即由该拐角的拐角前速度方向矢量e1和拐角后速度方向矢量e2所在的平面,与以该拐角为中心、并由数控机床的三个驱动轴最大加速度值的n倍:ax∈(-nAxm,nAxm),ay∈(-nAym,nAym),az∈(-nAzm,nAzm)为三维尺寸而确定的立方体相交而得到的总加速度六边形或四边形,其中Axm,Aym,Azm分别表示机床三个驱动轴的最大加速度;
(12)根据拐角速度大于零的条件,进一步确定该拐角总加速度的实际取值范围是:由该拐角的拐角前速度方向矢量e1和拐角后速度方向矢量e2所形成的小于180度角的扇形与该拐角总加速度六边形或四边形的相交区域;
(13)确定该拐角总加速度An:将该拐角总加速度实际取值范围中的除拐角总加速度六边形或四边形的中心点以外的其他各个顶点所对应的加速度值分别代入优化函数
Figure GSB00000125661600051
,式中,T是插补周期,即每次插补的单位时间间隔;点代表数乘,叉号表示叉乘,再选取其中f数值最大的点所对应的加速度为该拐角的总加速度An
(14)将该拐角总加速度代入下述公式,分别确定该拐角的拐角前速度V1和拐角后速度V2,以及该拐角的拐角前距离SP和拐角后距离EP:
V 1 = A n × e 2 e 2 × e 1 · T , V 2 = A n × e 1 e 2 × e 1 · T , SP = n - 1 2 · V 1 · T , EP = n + 1 2 · V 2 · T .
所述步骤(2)中的(22)进一步包括下列操作内容:
(22A)顺次判断每个小直线段两端速度是否满足可达性要求,如果是,则继续对后续的小直线段进行上述判断;否则,直接执行步骤(22B)操作;
(22B)对不满足上述要求的待插补小直线段的一端或两端的加工速度进行调整:降低该直线段中速度较快一端的拐角的过渡周期数来降低其相应的拐角速度,即每次比较两端速度大小并将速度较快一端的拐角的过渡周期数减1,并根据新的过渡周期数调整其拐角速度和其他最优插补参数,反复进行调整和判断后,直到该小直线段的两端速度满足上述要求后,执行步骤(22C)操作;
(22C)判断是否需要回溯处理:如果在步骤(22B)中,该小直线段的起始速度未被调整,则不需进行回溯处理,继续对下一个小直线段执行步骤(22A)的操作;否则,从该小直线段开始进行回溯处理,回溯处理完成后,返回到回溯开始时的小直线段的下一个小直线段,执行步骤(22A)的判断;所述回溯处理是从当前插补的小直线段开始,按照从后向前的顺序依次对其之前已经判断过可达性要求的小直线段再次进行可达性判断,对不满足可达性要求的小直线段执行步骤(22B)的调整;直到遇到第一次可达性判断就能满足要求的小直线段时,终止回溯处理。
所述步骤(22A)中,判断是否满足可达性要求的操作包括下列内容:
(22A1)计算每个小直线段的实际插补长度:该小直线段的长度Li分别减去该小直线段起点和终点处因拐角插补所占用的相应的拐角后距离EPi-1和拐角前距离SPi的差,即第i个小直线段的实际插补长度Lia=Li-EPi-1-WPi
(22A2)计算受限于数控机床各驱动轴的最大加速度的每个小直线段进行插补时的最大加速度为
Figure GSB00000125661600071
,式中,x,y,z分别是数控机床三个驱动轴的运动方向单位向量,u是该小直线段的单位方向矢量。
与现有技术相比,本发明数控系统基于多周期最优拐角对被加工曲面在精度范围内离散成的众多小直线段进行插补的方法具有以下创新的技术优点:
在满足加工精度的要求下,本发明能够有效提高在拐角处的切削加工速度。在相同的加工条件下,由于在小直线段加工中众多拐角处的速度方向的频繁变化,使得拐角速度的允许值较低,成为制约数控机床提高加工速度的瓶颈,而本发明所解决的提高拐角速度的直接后果,就是明显提高了数控机床的加工速度。
由于本发明采用前瞻处理方法,使得整个加工过程具有预见性,不会因为突然的过大的加、减速而产生机床振动,提高了被加工产品的精度并可延长机床的使用寿命;结合前述拐角多周期过渡的方法,能够整体上提高机床的加工速度,从而有效提高生产效率。
本发明方法中的操作步骤比较简单,且所有计算都是线性运算,复杂度低,运算速度快,能满足数控系统的实时加工要求,故适应性强。
再者,本发明具有很好的通用性能,只需改变数控系统的外部控制参数(例如最大加速性能指标、最大允许加工速度、插补周期、最大前瞻段数等),就可实现本发明方法在不同型号三轴数控机床间的移植应用。因此,本发明具有很好的推广应用前景。
附图说明
图1是传统小直线段插补算法的速度-时间示意图。
图2是二维小直线段拐角处等速率过渡插补算法的速度-时间示意图。
图3是现有技术中的拐角圆弧过渡方法示意图。
图4是对拐角采用多周期过渡的示意图。
图5是本发明数控系统基于多周期最优拐角的小直线段插补方法操作步骤流程图。
图6是本发明方法的步骤(1)确定最优拐角插补参数的操作步骤流程图。
图7(A)、(B)分别是拐角总加速度的两种取值范围的示意图。
图8是拐角总加速度实际取值范围的示例图。
图9是本发明方法的步骤(2)前瞻处理的操作步骤流程图。
图10是本发明方法的步骤(3)插补处理操作步骤流程图。
图11是本发明直线段上插补点序列和拐角处插补点序列衔接处理的示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面结合附图对本发明作进一步的详细描述。
首先,参见图4,具体介绍本发明提出的多周期拐角过渡的几个名词和概念:
本发明提出的多周期最优拐角过渡的思路是把原来在一个插补周期实现的拐角处加工速度方向的改变分散到多个插补周期中逐渐改变的一种加工方式,用于提高拐角处受限于机床各驱动轴的最大加速度和被加工路径拐角处几何参数的拐角速度。这里的拐角是加工路径中不在同一直线上的两个相邻小直线段的连接点(如图4中P点)。ε是拐角多周期过渡的实际插补路径偏离原始路径产生的误差。插补是机床数控系统依照设定方法确定刀具运动轨迹的过程,即按照某种算法计算已知点之间的中间点的方法,也被称为“数据点的密化”。插补周期是把加工一个小直线段的全部时间分成的多个相等的单位时间间隔T,每个插补周期T完成一次插补计算和加工。直线上插补是仍然在原始加工路径上的插补过程。拐角处插补则是根据该拐角的最优拐角插补参数确定在该拐角附近插补时的插补点序列的插补过程,图4中的折线SP-PE表示原始加工路径,折线SA-AB-BC-CE表示本发明采用拐角多周期过渡方法时拐角处的插补路径。点S是拐角插补的开始位置,点E是拐角插补的结束位置,P是拐角。
最优拐角插补参数包括拐角多周期过渡时的下述六个加工参数:
拐角的过渡周期数是改变拐角速度方向所用的插补周期个数,其初值n是根据经验选取的一个≥2的自然数,且后续操作步骤中,该数值会进行相应调整。
两个拐角速度,其中拐角前速度(用V1表示)是在拐角处插补开始时的加工速度,拐角后速度(用V2表示)是在拐角处插补结束时的加工速度。
两个拐角距离,其中拐角开始距离(图4中的SP)是拐角插补的开始位置与拐角的距离,拐角结束距离(图4中的EP)是拐角插补的结束位置与拐角的距离。
拐角总加速度(用An表示)是在一个插补周期内,由拐角前速度过渡到拐角后速度所需的数控机床各驱动轴加速度的矢量和。
需要说明的是,上述各个字符的涵义在后续使用时,意义相同。
参见图5,介绍本发明方法中的三个操作步骤的具体工作流程:
步骤1、基于加工精度和最大加工速度的限制,根据加工路径上拐角处的几何参数、数控机床各个驱动轴的最大加速度和实现“拐角前速度和拐角后速度之和为最大”的优化目标,确定加工路径上的每个拐角多周期过渡的最优拐角插补参数。
该步骤进一步包括下列操作内容(参见图6所示):
(11)确定加工路径上一个拐角在n个插补周期过渡时的拐角总加速度An(n的初值是根据经验预设的一个自然数)。
该步骤要先确定该拐角的拐角总加速度的取值范围。拐角是局限在其相邻两直线段所确定的一个平面内,称为拐角平面。图7(A)和(B)是以拐角为中心,机床三个驱动轴最大加速度的n倍为三维尺寸的长方体与拐角平面的示意图。拐角平面与长方体相交得到的六边形或四边形称为该拐角处速度方向的改变在n个插补周期过渡时拐角总加速度取值范围的六边形或者四边形。
然后确定该拐角总加速度的实际取值区域。拐角总加速度An还要满足下述约束:拐角前速度和拐角后速度都不能小于零,设拐角前速度方向矢量为e1,拐角后速度方向矢量为e2,则拐角总加速度An的取值范围是由方向矢量e1、e2所形成的小于180度角的扇形与拐角总加速度六边形或四边形的相交区域,被称为拐角总加速度实际取值区域。图8给出了一种拐角总加速度取值范围的示例图,图中带有斜线的多边形OABCD所围区域是拐角总加速度实际取值区域。
接着,根据优化目标确定该拐角总加速度。由于优化函数:
Figure GSB00000125661600101
是关于拐角总加速度An的线性函数,式中,T是插补周期;因此函数f的极大值在总加速度实际取值区域除中心点O以外的其他顶点处达到。分别计算函数f在各个顶点的取值,选取其中使得f最大的拐角总加速度为An
最后,根据该拐角的拐角总加速度An,分别确定该拐角的其他最优拐角插补参数:拐角速度V1,V2和拐角距离SP,EP,计算公式如下:
V 1 = A n × e 2 e 2 × e 1 · T , V 2 = A n × e 1 e 2 × e 1 · T , SP = n - 1 2 · V 1 · T , EP = n + 1 2 · V 2 · T .
(12)分别判断前述步骤确定的拐角的拐角前速度和拐角后速度是否大于最大加工速度,如果是,则降低该拐角的过渡周期数的初值n,直到其两个拐角速度都满足不大于最大加工速度要求为止;同时,对该该拐角的其它三个最优拐角插补参数:拐角前距离、拐角后距离和拐角总加速度也根据新的过渡周期数做相应调整;如果否,则不对该拐角的最优拐角插补参数作调整,直接执行下述步骤操作。
(13)判断拐角处多周期过渡插补的加工误差是否满足加工精度要求。根据拐角处插补时拐角插补的开始位置、拐角前速度和每个插补周期的加速度:该拐角的拐角总加速度除以该拐角的过渡周期数的商,确定该拐角插补时各个插补周期的插补点坐标,再根据插补点坐标判断因多周期过渡所引起的加工误差ε是否在产品加工允许的误差E范围内。如果否,则降低该拐角的拐角过渡周期数,并根据新的过渡周期数调整该拐角的其他最优拐角插补参数,直到其加工误差满足要求为止;如果是,则不对该拐角的最优拐角插补参数作调整,直接执行下述步骤操作。
其中加工误差ε的计算方法是:当拐角处的不在原来小直线段上的插补点个数是偶数或奇数时,该加工误差ε分别为该拐角到中间两个插补点连线的距离或该拐角到中间插补点的距离;
(14)对每个拐角根据加工顺序执行上述步骤(11)~(13)的操作,直到所有待加工小直线段均完成其拐角处最优拐角插补参数的确定,并保存最终得到的各个拐角的最优拐角插补参数。
步骤2、进行前瞻处理,对每个拐角的最优拐角插补参数进行调整,使每个小直线段的两端加工速度满足可达性要求,即在该小直线段实际插补长度范围内,其起点速度能够以该直线段插补时的最大加速度通过加速或减速方式,达到终点速度。
参见图9,介绍该步骤执行的前瞻处理方法的具体操作过程如下:
(21)根据系统配置,顺序选取加工路径上的每m个小直线段作为一组,分别进行前瞻处理,该m数值大小应满足数控机床的实时加工要求,并规定该m个待插补小直线段中的第一个小直线段的起始速度和第m个小直线段的终止速度都为零,且其对应拐角的其他最优拐角插补参数也都设置为零。
(22)顺次对一组待前瞻处理的小直线段中的每个小直线段的两端速度进行可达性判断。对于两端速度不满足可达性要求的小直线段,降低该小直线段一端或两端的拐角速度,直到两端速度满足可达性要求为止;同时对被调整拐角速度的该拐角的其他最优拐角插补参数也根据新的过渡周期数做相应调整。如果小直线段的起始速度被调整,则需要回溯处理。
该步骤是前瞻处理的关键操作,包括下列三个操作内容(参见图9):
(A)判断小直线段两端速度是否满足可达性要求,如果是,则执行步骤(B)的操作;否则,执行步骤(F)的操作。
(B)存储此时的最优拐角插补参数,并判断直线段当前的起始速度是否与读取该直线段时的起始速度值相同,如果是,则执行步骤(C)的操作,否则,执行步骤(E)的操作。
(C)判断当前直线段是否是m个前瞻处理小直线段的最后一个,如果是,执行步骤(23)的操作;否则,执行操作步骤(D)的操作。
(D)读取下一个小直线段,返回执行步骤(A)的操作。
(E)进行回溯处理,回溯完成后返回到回溯开始时的小直线段的下一个小直线段,并执行步骤(A)的操作;。
(F)将直线段中速度较快一端的拐角的过渡周期数减1,并据新的过渡周期数重新确定该拐角其他最优拐角参数,返回执行步骤(A)的操作。
所谓回溯处理即从当前插补的小直线段开始依次对其之前已经判断过可达性要求的小直线段再次进行可达性判断,对不满足可达性要求的小直线段进行与前瞻处理相同的调整方法;如此反复,直到遇到第一次可达性判断就能满足要求的小直线段时,终止回溯处理。
(23)从后续加工路径中,继续读取下一组m个小直线段,并对该m个小直线段执行步骤(22)的操作;直到加工路径上的所有小直线段都完成前瞻处理。
步骤3、根据前瞻处理后的最优拐角插补参数、被加工小直线段的长度和机床各驱动轴的最大加速度限制,对加工路径的各个小直线段分别进行直线段的插补和拐角的插补,并保证每个小直线段的直线段插补点序列和该直线段拐角的插补点序列准确衔接,实时顺序输出插补点序列,驱动数控机床进行实际加工。
参见图10,介绍该步骤3根据最优拐角插补参数确定直线段的插补点序列和拐角的插补点序列的相应操作内容。步骤如下:
(31)顺序读取一个待插补小直线段,确定该直线段上插补的插补点序列:
直线段上插补的插补点序列的起点是该直线段起点拐角处插补的最后一个插补点;根据直线段起始速度和终止速度:前述步骤已经调整好的起点拐角处的拐角后速度和终点拐角处的拐角前速度,先计算在该直线段上插补时以直线加减速方式运动时所能达到的最大速度Vmm,该Vmm数值不能大于最大加工速度Vm;再计算该直线段上的插补时间(包括加速时间、匀速时间和减速时间),最后,根据该直线段的起始速度、加速方式和相应加速方式下的时间进行该直线段的插补,并输出该直线段的插补点序列。
(32)在拐角处插补,确定该拐角处插补的插补点序列,并处理直线段上插补与拐角处插补的衔接问题(参见图11):
根据调整好的最优拐角插补参数,先由该拐角的拐角插补开始的位置、拐角前速度和该拐角每个插补周期的加速度确定该拐角处的插补点序列(如图11中的A2...A6)。实际插补时,在该小直线段的直线段上插补的最后一个插补点(如图11中B2点)和其拐角处的第一个插补点(如图11中A2点)不能在同一点重合时,即B2到A2的距离与以已规划好的速度在一个插补周期所能够加工的距离不相等时,则为了实现两者的准确衔接,以该直线段插补的最后一个插补点作为拐角插补的第一个插补点,再在上述规划好的插补点序列的基础上,根据残留比重新确定实际拐角插补点序列(如图11中B2...B6),实际拐角插补点序列的计算公式为:Bi=Ai-1+(1-δ)(Ai-Ai-1),其中,Bi是拐角处实际插补点,Ai是前述根据最优拐角插补参数确定拐角处的插补点,
Figure GSB00000125661600131
是残留比。
(33)对所有的待插补的小直线段重复执行步骤(31)和(32)的操作,直到所有的待插补微小直线段全部插补完毕。
(34)根据前述步骤得到的插补点序列,驱动数控机床进行实际加工。
本发明已经进行了多次实施例试验,实施例的多次仿真实验得到的实际数据证明:在满足精度误差的前提下,与图1所示的传统加工方式相比较,本发明方法能够提高数控机床整体加工速度达到100%,与图2所示的匀速率过渡方式相比较,也能提高数控机床整体加工速度的50%,因此,实施例的试验结果是成功的,实现了发明目的。

Claims (6)

1.一种数控系统基于多周期最优拐角的小直线段插补方法,其特征在于,所述方法包括下列操作步骤:
(1)根据数控机床各驱动轴的最大加速度、加工路径拐角处的几何参数、
“拐角前速度和拐角后速度之和为最大”的优化目标和每个拐角的过渡周期数的初值n,确定加工路径上的每个拐角在采用多个插补周期来改变加工速度方向的拐角总加速度;再根据该拐角的拐角总加速度分别计算其拐角速度和拐角距离,其中拐角速度包括拐角前速度和拐角后速度,拐角距离包括拐角前距离和拐角后距离,并将该拐角的拐角过渡周期数、拐角速度和拐角距离以及拐角总加速度,作为该拐角在其加工路径上采用多个插补周期来改变加工速度方向的过渡方式的拐角处的最优插补参数;
(2)进行前瞻处理:对每个拐角在所述拐角处的最优插补参数进行调整,使每个小直线段的两端加工速度满足可达性要求,即在该小直线段实际插补长度范围内,其起点速度能以该直线段插补时的最大加速度通过加速或减速方式,达到终点速度;也就是完成下述小直线段的前瞻处理操作内容:
(21)根据系统配置,顺序将加工路径上每m个连续小直线段作为一组进行前瞻处理,该m的数值大小应满足数控机床的实时加工要求,且规定该m个待插补小直线段中的第一个小直线段的起始速度和第m个小直线段的终止速度都为零,且其对应拐角的其他最优插补参数也都为零;
(22)顺次对一组待前瞻处理的小直线段的每个小直线段的两端速度进行可达性判断:即在该小直线段实际插补长度范围内,其起点速度能够以该直线段插补时的最大加速度通过加速或减速方式,达到终点速度;直到完成该m个小直线段的前瞻处理;
(23)继续顺序对后面的每组需要进行前瞻处理的直线段执行步骤(22)的操作;直到加工路径上的所有小直线段都完成前瞻处理后,保存前瞻处理后每个拐角的最优插补参数,以供后续步骤使用;
(3)根据前瞻处理后的最优插补参数、被加工小直线段的长度和机床各驱动轴的最大加速度限制,对加工路径的各个小直线段分别进行直线段的插补和拐角的插补,并保证每个小直线段的直线段插补点序列和该直线段拐角的插补点序列准确衔接,实时顺序输出插补点序列,驱动数控机床进行实际加工;即对所有的待插补的小直线段按照小直线段的输入顺序依次执行下述操作内容:
(31)顺序读取一个待插补小直线段及其拐角的前瞻处理后的最优插补参数,先进行直线段上插补,确定该直线段上插补的插补点序列:直线段上插补的插补点序列的起点是该直线段起点拐角处插补的最后一个插补点;根据直线段两端的速度,即前述步骤已经调整好的该小直线段的起点拐角处的拐角后速度和终点拐角处的拐角前速度,先计算在该直线段上插补时以直线加减速方式运动时所能达到的最大速度Vmm,该Vmm数值应不大于最大加工速度;再计算该直线段上的包括加速时间、匀速时间和减速时间的各个插补时间,最后,根据该直线段的起始速度、加速方式和相应加速方式下的时间进行该直线段的插补,得到该直线段的插补点序列,并实时输出插补点驱动数控机床进行实际加工,直线段插补完之后,接着执行步骤(32)的操作;
(32)在拐角处插补,确定该拐角处插补的插补点序列:先根据该拐角的拐角插补的开始位置、拐角前速度和该拐角每个插补周期的加速度:该拐角的拐角总加速度除以该拐角的过渡周期数的商,确定该拐角处的插补点序列;实际插补时,在该小直线段的直线段上插补的最后一个插补点和其拐角插补的第一个插补点不能在同一点重合时,则将该直线段上插补的最后一个插补点设置为该拐角插补点序列的第一个插补点,并在上述步骤所确定的拐角处插补点序列的基础上,根据残留比重新确定新的拐角处插补点序列,并实时输出插补点驱动数控机床进行实际加工。
2.根据权利要求1所述的方法,其特征在于:所述拐角是加工路径中不在同一直线上的两个相邻小直线段的连接点;所述拐角的多周期过渡是一种把拐角处加工速度方向在一个插补周期的改变分散到多个插补周期进行逐渐改变的加工方式;所述最优插补参数是拐角在多周期过渡时的下述六个加工参数:拐角的过渡周期数、拐角前速度V1、拐角后速度V2、拐角开始距离、拐角结束距离、和拐角总加速度An,其中,拐角的过渡周期数是改变拐角速度方向所用的插补周期个数,其初值n是根据经验选取的一个不小于2的自然数,且后续操作步骤中,该数值会进行相应调整;拐角前速度V1是在拐角处插补开始时的加工速度;拐角后速度V2是在拐角处插补结束时的加工速度;拐角开始距离是拐角插补的开始位置与拐角的距离,拐角结束距离是拐角插补的结束位置与拐角的距离,拐角总加速度An是在一个插补周期内,由拐角前速度过渡到拐角后速度所需的数控机床各驱动轴加速度的矢量和。
3.根据权利要求1所述的方法,其特征在于:所述方法在执行步骤(1)以后,还要继续执行下述操作内容:
(1A)分别判断前述步骤确定的拐角的拐角前速度和拐角后速度是否大于最大加工速度,如果是,则采用降低该拐角的过渡周期数来降低其拐角速度,直到该拐角的拐角前速度和拐角后速度都不大于最大加工速度为止;同时,对该拐角的其它三个最优插补参数:拐角前距离、拐角后距离和拐角总加速度也根据新的过渡周期数做相应调整,接着执行步骤(1B)操作;如果否,则不对该拐角的最优插补参数作调整,直接执行步骤(1B)操作;
(1B)根据步骤(1A)得到的最优插补参数,确定该拐角处插补的插补轨迹,并判断该拐角处的插补因偏离原始加工路径所产生的尺寸误差是否超出加工精度要求;如果是,则降低该拐角的过渡周期数,并对该拐角的其它最优插补参数也根据新的过渡周期数做相应调整,直到满足加工精度要求为止,并将最后得到的最优插补参数进行存储,以备后用,接着执行步骤(1C)操作;如果否,则不对该拐角的最优插补参数作调整,直接执行步骤(1C)操作;
(1C)对每个拐角顺序执行步骤(1)、(1A)、(1B)、(1C)的操作,直到对所有的拐角均确定其最优插补参数。
4.根据权利要求1所述的方法,其特征在于:所述步骤(1)进一步包括下列操作内容:
(11)读取直线插补方式的数控加工代码,先初步确定拐角多周期过渡时该拐角总加速度的取值范围,即由该拐角的拐角前速度方向矢量e1和拐角后速度方向矢量e2所在的平面,与以该拐角为中心、并由数控机床的三个驱动轴最大加速度值的n倍:ax∈(-nAxm,nAxm),ay∈(-nAym,nAym),az∈(-nAzm,nAzm)为三维尺寸而确定的立方体相交而得到的总加速度六边形或四边形,其中Axm,Aym,Azm分别表示机床三个驱动轴的最大加速度;
(12)根据拐角速度大于零的条件,进一步确定该拐角总加速度的实际取值范围是:由该拐角的拐角前速度方向矢量e1和拐角后速度方向矢量e2所形成的小于180度角的扇形与该拐角总加速度六边形或四边形的相交区域;
(13)确定该拐角总加速度An:将该拐角总加速度实际取值范围中的除拐角总加速度六边形或四边形的中心点以外的其他各个顶点所对应的加速度值分别代入优化函数
Figure FSB00000125661500041
式中,T是插补周期,即每次插补的单位时间间隔;点代表数乘,叉号表示叉乘;再选取其中f数值最大的点所对应的加速度为该拐角总加速度An
(14)将该拐角总加速度代入下述公式,分别确定该拐角的拐角前速度V1和拐角后速度V2,以及该拐角的拐角前距离SP和拐角后距离EP:
Figure FSB00000125661500042
Figure FSB00000125661500043
Figure FSB00000125661500044
5.根据权利要求1所述的方法,其特征在于:所述步骤(2)中的(22)进一步包括下列操作内容:
(22A)顺次判断每个小直线段两端速度是否满足可达性要求,如果是,则继续对后续各个小直线段分别进行上述判断;否则,直接执行步骤(22B)操作;
(22B)对不满足上述要求的待插补小直线段的一端或两端的加工速度进行调整:降低该直线段中速度较快一端的拐角的过渡周期数来降低其相应的拐角速度,即每次比较两端速度大小并将速度较快一端的拐角的过渡周期数减1,并跟据新的过渡周期数调整其拐角速度和其他最优插补参数,反复进行调整和判断后,直到该小直线段的两端速度满足上述要求后,执行步骤(22C)操作;
(22C)判断是否需要回溯处理:如果在步骤(22B)中,该小直线段的起始速度未被调整,则不需进行回溯处理,继续对下一个小直线段执行步骤(22A)的操作;否则,从该小直线段开始进行回溯处理,回溯处理完成后,返回到回溯开始时的小直线段的下一个小直线段,执行步骤(22A)的判断;所述回溯处理是从当前插补的小直线段开始,按照从后向前的顺序依次对其之前已经判断过可达性要求的小直线段再次进行可达性判断,对不满足可达性要求的小直线段执行步骤(22B)的调整;直到遇到第一次可达性判断就能满足要求的小直线段时,终止回溯处理。
6.根据权利要求5所述的方法,其特征在于:所述步骤(22A)中,判断是否满足可达性要求的操作包括下列内容:
(22A1)计算每个小直线段的实际插补长度:该小直线段的长度Li分别减去该小直线段起点和终点处因拐角插补所占用的相应的拐角后距离EPi-1和拐角前距离SPi的差,即第i个小直线段的实际插补长度Lia=Li-EPi-1-SPi
(22A2)计算受限于数控机床各驱动轴的最大加速度的每个小直线段进行插补时的最大加速度为
Figure FSB00000125661500051
式中,x,y,z分别是数控机床三个驱动轴的运动方向单位向量,u是该小直线段的单位方向矢量。
CN2009100839504A 2009-05-13 2009-05-13 数控系统基于多周期最优拐角的小直线段插补方法 Expired - Fee Related CN101556472B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2009100839504A CN101556472B (zh) 2009-05-13 2009-05-13 数控系统基于多周期最优拐角的小直线段插补方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2009100839504A CN101556472B (zh) 2009-05-13 2009-05-13 数控系统基于多周期最优拐角的小直线段插补方法

Publications (2)

Publication Number Publication Date
CN101556472A CN101556472A (zh) 2009-10-14
CN101556472B true CN101556472B (zh) 2010-08-25

Family

ID=41174610

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009100839504A Expired - Fee Related CN101556472B (zh) 2009-05-13 2009-05-13 数控系统基于多周期最优拐角的小直线段插补方法

Country Status (1)

Country Link
CN (1) CN101556472B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107272758A (zh) * 2017-08-01 2017-10-20 深圳市雷赛控制技术有限公司 绕线设备效率及平稳性的提升方法及装置

Families Citing this family (26)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102010001781A1 (de) * 2010-02-10 2011-08-11 Siemens Aktiengesellschaft, 80333 Verfahren zur Bewegung eines Maschinenelements einer Maschine aus der Automatisierungstechnik und Steuereinrichtung
JP5255108B2 (ja) 2011-11-04 2013-08-07 ファナック株式会社 指令経路速度条件による速度制御を行う数値制御装置
CN102722140B (zh) * 2012-06-21 2013-12-25 中国科学院数学与系统科学研究院 基于s曲线加减速控制的多周期拐角小直线段插补方法
CN102799145B (zh) * 2012-08-13 2014-08-06 中国科学院数学与系统科学研究院 拐角多周期恒加加速度过渡的s曲线加减速直线插补方法
CN103454968B (zh) * 2012-12-10 2016-01-20 深圳信息职业技术学院 拐角高速铣削的刀轨生成方法
CN103135501B (zh) * 2013-02-01 2015-10-21 北京配天技术有限公司 基于s型曲线的加减速控制方法和装置及数控机床
CN104252153A (zh) * 2013-06-28 2014-12-31 鸿富锦精密工业(深圳)有限公司 Cnc加工程序生成系统及方法
CN103777573B (zh) * 2013-12-10 2016-06-08 固高科技(深圳)有限公司 减小轮廓加工误差的方法和系统
CN105607577B (zh) * 2014-11-20 2018-07-13 沈阳高精数控智能技术股份有限公司 一种基于轴向精度的小线段过渡方法
CN106054817B (zh) * 2016-06-08 2019-01-11 上海交通大学 局部最优小线段路径解析光顺的实时前瞻插补方法及系统
CN106227152A (zh) * 2016-10-17 2016-12-14 合肥工业大学 一种cnc加工的平滑过渡方法及其平滑过渡装置
CN108375957B (zh) * 2017-01-31 2021-03-26 兄弟工业株式会社 数值控制装置和控制方法
CN107329458A (zh) * 2017-07-10 2017-11-07 深圳市旗众智能自动化有限公司 一种连续小线段轨迹的动态自适应速度前瞻控制方法
CN108227630B (zh) * 2018-01-18 2020-08-18 合肥工业大学 一种采用时间参数多项式插补的自由曲面数控加工方法
CN108661850B (zh) * 2018-04-13 2020-03-24 中国航天空气动力技术研究院 一种系绳式洋流发电机运动轨迹控制方法
CN109048091B (zh) * 2018-07-17 2020-12-25 大族激光科技产业集团股份有限公司 激光切割速度规划方法、装置、存储介质及计算机设备
CN110874082B (zh) * 2018-08-31 2021-08-20 大族激光科技产业集团股份有限公司 工件拐角的切割路径确定方法、设备和存储介质
JP7037457B2 (ja) * 2018-09-05 2022-03-16 オークマ株式会社 数値制御装置および数値制御方法
CN108829045B (zh) * 2018-09-06 2019-12-20 深圳市雷赛控制技术有限公司 连续微直线段的衔接速度的优化方法及系统
CN111390902B (zh) * 2020-01-06 2022-09-06 广州视源电子科技股份有限公司 轨迹规划方法、轨迹规划装置、设备和存储介质
CN111880484B (zh) * 2020-07-15 2021-09-21 中国科学院大学 一种数控机床加工前瞻处理方法及系统
CN111736532A (zh) * 2020-07-15 2020-10-02 中国科学院大学 一种基于直线段和圆弧路径的多周期最优拐角的插补方法
CN111736537A (zh) * 2020-07-21 2020-10-02 天津大学 一种自由曲面加工中双nurbs路径极限速度的计算方法
CN112157260B (zh) * 2020-09-30 2021-08-27 中南大学 一种用于激光熔覆3d打印转角优化的方法
CN115616982B (zh) * 2021-07-12 2024-06-04 大族激光科技产业集团股份有限公司 一种拐弯轨迹的过度方法及激光装置和计算机设备
CN114609977B (zh) * 2022-03-23 2023-08-04 大族激光科技产业集团股份有限公司 速度规划方法、装置、设备及存储介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4866631A (en) * 1986-10-21 1989-09-12 Sony Corporation Method for generating offset surface data
US5276383A (en) * 1990-08-30 1994-01-04 Kabushiki Kaisha Kobe Seiko Sho Route interpolation method for robot
EP1398680A2 (en) * 1999-11-19 2004-03-17 Siemens Energy & Automation, Inc. Apparatus and method for smooth cornering in a motion control system
CN1971457A (zh) * 2005-11-25 2007-05-30 中国科学院沈阳计算技术研究所有限公司 一种用于数控机床的加减速控制方法
CN101211177A (zh) * 2006-12-29 2008-07-02 中国科学院沈阳计算技术研究所有限公司 一种基于滤波技术的数控系统加减速控制方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4866631A (en) * 1986-10-21 1989-09-12 Sony Corporation Method for generating offset surface data
US5276383A (en) * 1990-08-30 1994-01-04 Kabushiki Kaisha Kobe Seiko Sho Route interpolation method for robot
EP1398680A2 (en) * 1999-11-19 2004-03-17 Siemens Energy & Automation, Inc. Apparatus and method for smooth cornering in a motion control system
CN1971457A (zh) * 2005-11-25 2007-05-30 中国科学院沈阳计算技术研究所有限公司 一种用于数控机床的加减速控制方法
CN101211177A (zh) * 2006-12-29 2008-07-02 中国科学院沈阳计算技术研究所有限公司 一种基于滤波技术的数控系统加减速控制方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
JP特开2003-295918A 2003.10.17
张园,陈友东,黄荣瑛,魏洪兴,邹勇.高速加工中连续微小线段的前瞻自适应插补算法.机床与液压.2008,36(6),1-4. *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107272758A (zh) * 2017-08-01 2017-10-20 深圳市雷赛控制技术有限公司 绕线设备效率及平稳性的提升方法及装置
CN107272758B (zh) * 2017-08-01 2020-08-07 深圳市雷赛控制技术有限公司 绕线设备效率及平稳性的提升方法及装置

Also Published As

Publication number Publication date
CN101556472A (zh) 2009-10-14

Similar Documents

Publication Publication Date Title
CN101556472B (zh) 数控系统基于多周期最优拐角的小直线段插补方法
CN102722140A (zh) 基于s曲线加减速控制的多周期拐角小直线段插补方法
CN101539769B (zh) 基于二次b样条曲线对g01代码的拟合及插补方法
EP0145967B1 (en) Curvilinear interpolation system and method
CN102419570B (zh) 数控机床高速加工的加减速前瞻控制方法
CN104615084B (zh) 加工进给速度优化的刀轨曲线轮廓误差补偿方法
CN102640067A (zh) 刀具路径的形成方法和装置
CN101653921A (zh) 凸轮轴数控磨削轮廓误差补偿方法
CN108062073A (zh) 一种用于高质量加工的圆弧平滑压缩插补算法
CN106094737B (zh) 一种指定加工误差条件下的数控加工速度优化控制方法
CN104155916B (zh) 一种高精度快速空间圆弧插补方法
CN115122329B (zh) 一种便捷的工业机器人关节空间轨迹规划方法
CN111736532A (zh) 一种基于直线段和圆弧路径的多周期最优拐角的插补方法
KR900007161B1 (ko) 영역 가공 방법
Min et al. An improved b-spline fitting method with arc-length parameterization, g 2-continuous blending, and quality refinement
CN102414633A (zh) 数值控制装置
Chu et al. Generation of reciprocating tool motion in 5-axis flank milling based on particle swarm optimization
JP2006227701A (ja) 円弧加工指令作成装置及びその方法並びにプログラム
CN112506143A (zh) 基于s型曲线的高质量加工五次多项式速度规划方法
KR100880734B1 (ko) 기계 동작특성에 근거한 5축 밀링가공의 속도향상을 위한 보간길이의 최적화 방법
CN111487928A (zh) 一种基于刀位点增删改指令的数控加工轨迹平滑方法
CN111708327B (zh) 一种具有g2连续性的ph样条过渡线性路径加工方法
Zhao et al. An analytical local tool path smoothing method with curvature control
JP2003223208A (ja) 数値制御システム
CN111880484A (zh) 一种数控机床加工前瞻处理方法及系统

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20100825

Termination date: 20170513