CN108073138A - 适用于高速高精加工的椭圆弧平滑压缩插补算法 - Google Patents
适用于高速高精加工的椭圆弧平滑压缩插补算法 Download PDFInfo
- Publication number
- CN108073138A CN108073138A CN201610978793.3A CN201610978793A CN108073138A CN 108073138 A CN108073138 A CN 108073138A CN 201610978793 A CN201610978793 A CN 201610978793A CN 108073138 A CN108073138 A CN 108073138A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- elliptic arc
- arc
- point
- 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
Links
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B19/00—Programme-control systems
- G05B19/02—Programme-control systems electric
- G05B19/18—Numerical 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/41—Numerical 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/4103—Digital interpolation
Landscapes
- Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Theoretical Computer Science (AREA)
- Human Computer Interaction (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Numerical Control (AREA)
Abstract
本发明涉及一种适用于高速高精加工的椭圆弧平滑压缩插补算法,包括以下步骤:识别连续微小线段加工区域;在连续微小线段加工区域中选取型值点并对其进行拟合得到二次有理Bézier曲线;根据曲线特征识别出椭圆弧,并对其进行几何形式转换;将属于同一椭圆的相邻椭圆弧合并,得到插补曲线;在插补曲线上进行椭圆弧插补。本发明方法利用样条曲线的特征,识别出加工路径中的椭圆弧,用其几何形式进行插补,能够精确地计算出弧长对应的插补参数,降低计算量和加工中速度的频繁波动,实现高速高精加工。
Description
技术领域
本发明涉及高质量加工中参数样条曲线的拟合及椭圆弧段的识别与合并,属于数控加工技术领域。
背景技术
椭圆是圆锥曲线中的一种,在工件加工程序中较为常见。由于椭圆弧长的计算涉及到椭圆积分问题,且其原函数不能用有限形式表示,所以无法直接对椭圆进行精确插补。又因其具有较强的工程实用性,国内外的学者对此进行了广泛的研究。
如映射法,将椭圆映射到某一平面,使其投影为圆,通过对圆进行插补,然后利用坐标变换得到椭圆上的插补点。离心角增量法,在弧长增量相对椭圆周长非常小的情况下,通过下一插补点与当前插补点和离心角增量的关系进行插补。上述方法虽简单可行,但都存在轮廓误差和速度波动无法控制的问题。于是提出了等误差方法,虽然提高了椭圆的加工精度,但仍无法有效地控制速度波动。因此大量学者对数值方法进行了研究,通过较为精确地推导出插补点与弧长的关系,来有效地降低速度波动,但这些方法因涉及到多次迭代与复杂数值计算,实时性较差,所以均无法应用于数控***。
随着设计和制造技术的发展,越来越多的人使用计算机辅助设计(CAD)***来进行复杂零部件的设计。但由于绝大多数的数控***不支持参数样条数据的传输,通常利用计算机辅助制造(CAM)***将CAD设计的自由曲线或曲面,在特定的容差范围内用一系列折线去覆盖,从而生成由大量指令点构成的数控加工程序。
通过对微小直线段进行数控加工,可以不用考虑复杂的椭圆弧长计算问题。但由于加速度和加加速度的频繁变化,会引起机床的震动,降低加工的效率和质量。而现在对于微小线段高速加工的研究主要分为两种方法。一种是在相邻微小线段的拐角处***过渡样条曲线。例如,通过在拐角处***三次B样条曲线或四次Bézier曲线的方式,提高拐角处的速度,从而提高加工效率,但由于插补点在直线段和过渡样条曲线段上的循环出现,使得插补步长不一致,从而导致加工速度的波动,如果指令点越密集,速度的波动就会越频繁。另一种则是通过插值或者逼近的方式将离散的指令点拟合成平滑的加工路径。例如,将由连续微小线段指定的加工路径转化为由二次Bézier曲线表示的加工路径,通过对Bézier曲线插补来实现对自由曲线的高速高精加工。虽然这种方式能较好地逼近原设计曲线,但由于拟合曲线较复杂,无法精确计算插补步长对应的插补参数,导致加工速度波动,降低了加工精度。
发明内容
为了克服直接插补算法和现有小线段插补算法不能高精加工和存在速度波动的不足,本发明的目的是提供一种参数样条拟合方法,通过参数样条曲线特征,识别加工路径中的椭圆弧段,通过对椭圆弧段进行插补,来保证加工的精度,同时降低加工速度的波动。
本发明解决其技术问题所采用的技术方案是:一种适用于高速高精加工的椭圆弧平滑压缩插补算法,包括以下步骤:
识别连续微小线段加工区域;
在连续微小线段加工区域中选取型值点并对其进行拟合得到二次有理Bézier曲线;
根据曲线特征识别椭圆弧及其几何形式转换;
将属于同一椭圆的相邻椭圆弧合并,得到插补曲线;
在插补曲线上进行椭圆弧插补。
所述进行拟合包括以下步骤:
2-1)求取二次有理Bézier曲线的权重;
2-2)根据权重通过求取平均值得到二次有理Bézier的拟合曲线。
所述权重通过下式得到:
其中,w1为权重,P0为首型值点、P2为末型值点、P1为控制点、P为指令点,、u为P0Q与QP2的比值,Q为以P1为投影中心,直线段[P0 P2]投影到二次有理Bézier拟合曲线的投影点。
所述根据权重通过求取平均值得到二次有理Bézier的拟合曲线包括以下步骤:
4-1)根据权重的求取公式得到求取二次有理Bézier曲线对应的权值wk及相应的肩点坐标sk=wk/(1+wk),k=i+1,...,j-1;
4-2)通过对sk求平均值得到平均肩点坐标s和平均权值w=s/(1-s);
4-3)根据P0,P1,P2和w确定型值点Qi与Qj之间二次有理Bézier的拟合曲线。
所述根据曲线特征识别椭圆弧及其几何形式转换包括以下步骤:
5-1)当|P0P1|!=|P1P2|且0<w1<1时,则该曲线为椭圆弧;
5-2)通过以下步骤得到椭圆弧的几何信息:
Ci(u)为加工路径中某段的标准型二次有理Bézier曲线,且其对应于一个椭圆弧段,根据Ci(u)可得到该椭圆弧段的起点坐标p_start和终点坐标p_end;该椭圆弧段的中心坐标为,
P1+ε(S+T)
该椭圆弧段的长、短半径分别为,
假设λ2>λ1>0,且是如下二次方程的根,
2δλ2-(kη+4β)λ+2(k-1)=0
δ=|S×T|2,η=|S-T|2,β=S·T
该椭圆弧段长轴上的两个点为,
Q1=P1+(ε+r1x0)S+(ε+r1y0)T
Q2=P1+(ε-r1x0)S+(ε-r1y0)T
根据Q1、Q2,可以得到长轴的斜率kl,和长轴与x轴正半轴的夹角d_kl;以长轴为x′轴,短轴为y′轴建立局部坐标系,利用Q1、Q2、中心坐标、椭圆弧段起点坐标和椭圆弧段终点坐标,求得椭圆弧段起点与中心连线、椭圆弧段终点与中心连线分别相对于局部坐标系x′正半轴的起始角度d_start、终点角度d_end,以及椭圆弧的方向turn;
当根据曲线特征识别椭圆弧及其几何形式转换结束后,得到一个椭圆弧段数组Ellipse_Arcs[],数组中每个数据的结构如下:
本发明具有以下有益效果及优点:
1.本发明方法控制简单,可以有效地降低椭圆弧加工速度的频繁波动,实现椭圆弧的高质量加工。
2.压缩量大,平滑度高。本发明方法能够识别出数控加工程序中的椭圆弧,并能用椭圆弧的几何形式进行表示,极大地减少程序段的数量,同时椭圆弧的几何表示形式提高了加工路径的光滑度。
3.加工精度和加工效率高。本发明方法用几何形式的椭圆弧进行插补,能够较为精确地计算出弧长对应的插补参数,降低了计算的复杂度和加工中速度的频繁波动,提高了加工质量和加工效率。
附图说明
图1为本发明的方法流程图
图2为连续微小线段加工区域的识别示意图;
图3为局部曲率最大值的识别示意图;
图4为二次有理Bézier曲线示意图;
图5为型值点切矢的计算原理图;
图6为椭圆弧几何信息的计算原理图;
图7为区域及变量的划分示意图。
具体实施方式
下面结合实施例对本发明做进一步的详细说明。
本发明是一种适用于高速高精加工的椭圆弧平滑压缩插补算法,根据双弓高误差限制将加工路径划分为非连续微小线段加工区域和连续微小线段加工区域。对于非连续微小线段加工区域,直接在相邻指令点构成的直线段上进行插补计算,以保证加工精度。对于连续微小线段加工区域,根据离散指令点的曲率值,对曲率极值点和拐点进行拟合,将折线加工路径转化为平滑的二次有理Bézier曲线加工路径;然后,利用二次有理Bézier曲线特征识别出椭圆弧,并转换为几何形式;最后,将相邻椭圆弧段合并后,进行插补计算。
本发明提供了一种适用于高速高精加工的椭圆弧平滑压缩插补算法,包括以下步骤:
1.加工区域的识别,根据双弓高误差判断条件,将相邻两个不满足条件的点及它们之间的指令点称为连续微小线段加工区域。
2.型值点的选取,对于连续微小线段加工区域中的指令点,通过离散点曲率计算公式,计算其曲率值,并根据相邻指令点的曲率值和判断条件,找出加工路径中局部曲率最大值点和拐点,将连续微小线段加工区域的两端点、局部曲率最大值点和拐点标记为型值点。
3.型值点的拟合,在保证加工精度的条件下,根据型值点的坐标值和单位切矢,将指令点指定的折线加工路径转化为平滑的二次有理Bézier曲线加工路径。
4.椭圆弧的识别,根据二次有理Bézier曲线的特征识别出椭圆弧段,并将椭圆弧的二次有理Bézier形式转换为几何形式。
5.椭圆弧段的合并,对相邻的椭圆弧段,根据其几何信息,识别他们是否属于同一椭圆,对于属于同一椭圆的相邻椭圆弧段进行合并。
6.椭圆弧的插补,通过在几何形式的椭圆弧上进行插补计算,实现椭圆弧的高速高精加工。
如图1所示,本发明提出了一种适用于高速高精加工的椭圆弧平滑压缩插补算法,解决了小线段插补的问题,方法由加工区域的识别、型值点的选取、型值点的拟合、椭圆弧的识别与几何形式的转换、椭圆弧段的合并和椭圆弧的插补6个部分组成,提高了加工的质量和效率。
根据双弓高误差判断条件,将相邻两个不满足条件的点及它们之间的指令点称为连续微小线段加工区域。
对于连续微小线段加工区域中的指令点,通过离散点曲率计算公式,计算其曲率值,并根据相邻指令点的曲率值和判断条件,找出加工路径中局部曲率最大值点和拐点,将连续微小线段加工区域的两端点、局部曲率最大值点和拐点标记为型值点。
保证加工精度的条件下,根据型值点的坐标值和单位切矢,将指令点指定的折线加工路径转化为平滑的二次有理Bézier曲线加工路径。
根据二次有理Bézier曲线的特征识别出椭圆弧段,并将椭圆弧的二次有理Bézier形式转换为几何形式。
对相邻的椭圆弧段,根据其几何信息,识别他们是否属于同一椭圆,对于属于同一椭圆的相邻椭圆弧段进行合并。
椭圆弧的插补通过在几何形式的椭圆弧上进行插补计算,实现椭圆弧的高速高精加工。
本发明的具体步骤如下:
1.连续微小线段加工区域的识别
如图2所示,Pi-1、Pi和Pi+1为顺序相邻的三个指令点,l1、l2为小线段的段长,θ为小线段间的拐角,双弓高误差判断条件如下所示,
其中,δ1、δ2分别为小线段Pi-1Pi和PiPi+1的弓高误差,φ1为OPi-1和OPi夹角的一半。φ2为OPi和OPi+1夹角的一半。
如果δ1或δ2大于设定的最大弓高误差值δmax,则Pi为断点;若两个相邻的断点间存在指令点,那么两个断点连同它们之间的指令点就构成了一个连续微小线段加工区域。
2.型值点的选取
为了减少对连续微小线段加工区域的拟合次数、增加程序段的压缩量,通过如下三个步骤来选取型值点。
(1)将连续加工区域的开始点和结束点,即断点,标记为型值点。
(2)将局部曲率最大值点标记为型值点。
如图2所示,三个相邻指令点的坐标为Pi-1(xi-1,yi-1)、Pi(xi,yi)和Pi+1(xi+1,yi+1),离散指令点Pi的曲率值由下式确定,
其中ΔPi-1PiPi+1为带符号的三角型面积由下式确定,
假设kl为Pi左边的局部曲率最小值,kr为Pi右边的局部曲率最小值,如果满足下面两个条件,则Pi为局部曲率最大值点,如图3所示。
1)|ki|>|kl|并且|ki|>|kr|
2)|ki|-|kl|≥δf或者|ki|-|kr|≥δf,δf为设定的最大曲率差值。
(3)将加工路径弯曲方向改变的点,即拐点,标记为型值点。
利用第(2)步计算出的指令点曲率值进行判断,如果ki-1ki>0并且kiki+1<0,ki-1、ki+1分别为指令点Pi-1和Pi+1的曲率值,则Pi为拐点。
3.型值点的拟合
对于连续微小线段加工区域中的n个型值点,可用一条n-1段的二次有理Bézier曲线进行拟合,以达到平滑压缩加工路径的目的。
标准型的二次有理Bézier曲线如下所示,
其中,P0、P1和P2为控制点,w1为权值,u为变量,范围为[0,1]。
由图4可知,当给定首末端点P0和P2,以及这两点处的切线方向T0和T2时,可以很容易地通过直线[P0 T0]和[P2 T2]的交点得到P1,再给定一点P就能唯一确定该曲线,因而确定了w1。
将所求的曲线看成是由点P0,P1和P2确定的抛物线的投影,P1为投影中心。如图4所示,将直线段[P0 P2]投影到要求的曲线上,则点P和Q为相对应的投影点。令w1=0,得到直线段L(u)=[P0 P2],即
L(u)是点P0和P2的凸组合,因此|P0Q|和|QP2|的比值为u2:(1-u)2,从而推出
将u和P带入二次有理Bézier曲线,得到w1,从而得到所求曲线。
假定型值点Qi与Qj之间由指令点Qi+1,Qi+2,…,Qj-2与Qj-1指定的折线加工路径组成,且型值点切矢已知,则型值点Qi与Qj之间的拟合详细步骤如下所述。
(1)构造以P0=Qi,P1和P2=Qj为控制点的拟合曲线,使其分别插值于Qk(k=i+1,...,j-1)。这将产生中间的权值wk及相应的肩点坐标sk=wk/(1+wk),k=i+1,...,j-1;
(2)通过对sk求平均值得到逼近曲线的肩点坐标,即
则中间的权值为w=s/(1-s);
(3)根据P0,P1,P2和w即可确定型值点Qi与Qj之间的拟合曲线。
4.型值点切矢的计算
由于数控加工程序中并不提供指令点处的切矢量,但可以通过型值点以及它周围的四个指令点来计算型值点处的切矢量。如图5所示,Qi为型值点,Qi-2、Qi-1、Qi+1、Qi+2为它周围的四个指令点,Qi的单位切矢表达式如下
由于计算公式中未使用参数u,所求的单位切矢只能被看作切矢的方向。利用下式,即可得到T0,T1和Tn-1,Tn。
q0=2q1-q2,q-1=2q0-q1
qn+2=2qn+1-qn,qn+1=2qn-qn-1
5.拟合精度的控制
虽然曲线段经过型值点Qi与Qj,但不能保证它到Qi与Qj之间的所有指令点的距离都满足最大拟合误差。因此,将Qi与Qj之间的所有指令点投影到拟合的曲线上,检查它们的偏差是否都在允许的误差范围内,即小于等于***设定的最大轮廓误差δc,如果满足误差要求,进行下一段的曲线拟合;否则,将偏差最大的指令点设置为新的型值点,用新的型值点重新进行曲线拟合,再次检验拟合精度,重复这个过程,直至满足误差要求。
6.椭圆弧的识别与几何形式的转换
由于标准型的二次有理Bézier曲线只有一个权因子,其表达能力比非有理形式的Bézier曲线强,可以表达多种曲线,当|P0P1|!=|P1P2|且0<w1<1时,则所表达的曲线为椭圆弧。
如果Ci(u)为加工路径中某段的标准型二次有理Bézier曲线,且其对应于一个椭圆弧段,如图6所示,根据Ci(u)可得到该椭圆弧段的起点坐标p_start和终点坐标p_end。由椭圆的标准型二次有理Bézier形式与其几何形式的关系可知,该椭圆弧段的中心坐标为,
P1+ε(S+T)
该椭圆弧段的长、短半径分别为,
假设λ2>λ1>0,且是如下二次方程的根,
2δλ2-(kη+4β)λ+2(k-1)=0
δ=|S×T|2,η=|S-T|2,β=S·T
该椭圆弧段长轴上的两个点为,
Q1=P1+(ε+r1x0)S+(ε+r1y0)T
Q2=P1+(ε-r1x0)S+(ε-r1y0)T
根据Q1、Q2,可以得到长轴的斜率kl,和长轴与x轴正半轴的夹角d_kl;以长轴为x′轴,短轴为y′轴建立局部坐标系,利用Q1、Q2、中心坐标、椭圆弧段起点坐标和椭圆弧段终点坐标,求得椭圆弧段起点与中心连线、椭圆弧段终点与中心连线分别相对于局部坐标系x′正半轴的起始角度d_start、终点角度d_end,以及椭圆弧的方向turn;x0、y0为椭圆弧段起点坐标。
当连续微小线段加工区域完成椭圆弧的识别和几何形式转换后,会得到一个椭圆弧段数组Ellipse_Arcs[],其中椭圆弧段的数据结构如下:
7.椭圆弧的合并
为了提高椭圆弧加工的效率和精度,将连续微小线段加工区域内的椭圆弧段合并,如果相邻的椭圆弧段具有相同的中心center、长半径r_l、短半径r_s、方向turn、长轴与x正半轴的夹角d_kl,并且上一椭圆弧段的终点坐标p_endi-1与下一椭圆弧段的起点坐标p_starti相等,则它们属于同一椭圆弧上连续的椭圆弧段,通过修改上一椭圆弧段的终点角度d_endi-1和终点坐标p_endi-1将它们合并为一个椭圆弧段。重复上述步骤,直至椭圆弧段不能再合并。
8.椭圆弧的插补
8.1椭圆弧长的计算
假设椭圆方程为其中a为长半轴半径,b为短半轴半径,如图7所示,将椭圆划分为四个区域,每个区域对应不同的参数作为变量。
以第一象限y为变量为例,则该区域内任意两点间的弧长为,
其中yi、yi+1分别表示要求该区域内任意两点间弧长的起点横坐标和终点横坐标。
根据高斯勒让德求积公式,可将上式转换为,
其它象限x或y为变量的情况,其弧长的推导方法相同。
8.2弧长与变量的关系
以第一象限y为变量,从点E0(a,0)到点的椭圆弧为例,其弧长与变量的关系求解步骤如下:
(1)在变量y的区间内,取两个三等分点y1和y2。
(2)根据8.1的弧长公式分别计算yi+1=0、y1、y2和到椭圆弧的起点yi=0的弧长,s0、s1、s2、s3。
(3)将(s0,0)、(s1,y1)、(s2,y2)和代入如下三次多项式中进行插值,求得多项式的系数,即得到弧长与变量的关系。
y=a0+a1s+a2s2+a3s3,s∈[s0,s3]
其中,a0~a3为系数。其它象限x或y为变量的情况,其弧长与变量关系的推导方法相同。
8.3椭圆弧的插补
通过上述过程便可将数控加工路径中的椭圆弧识别出来,并用几何的形式表示。对于识别出来的椭圆弧段,按照四个区域对其进行划分,并分别计算出各个区域内弧长与变量的关系,然后利用这些弧长与变量的关系来计算第i个插补周期的插补点坐标(xi′,yi′),具体步骤如下所示。
假定插补周期为T,P0′(x′i-1,y′i-1)和Si-1分别为第i-1个插补周期的插补点坐标与弧长,第i个插补周期的进给速度为vi-1,则弧长增量为ΔS=vi-1T,第i个插补周期的弧长为Si=Si-1+ΔS,将Si代入弧长与变量的关系中,得到对应变量xi′或yi′的值,再将其代入椭圆方程中,从而求得第i个插补周期的插补点坐标P′(xi′,yi′)。
Claims (6)
1.一种适用于高速高精加工的椭圆弧平滑压缩插补算法,其特征在于包含以下步骤:
识别连续微小线段加工区域;
在连续微小线段加工区域中选取型值点并对其进行拟合得到二次有理Bézier曲线;
根据曲线特征识别椭圆弧及其几何形式转换;
将属于同一椭圆的相邻椭圆弧合并,得到插补曲线;
在插补曲线上进行椭圆弧插补。
2.根据权利要求1所述的一种适用于高速高精加工的椭圆弧平滑压缩插补算法,其特征在于所述进行拟合包括以下步骤:
2-1)求取二次有理Bézier曲线的权重;
2-2)根据权重通过求取平均值得到二次有理Bézier的拟合曲线。
3.根据权利要求2所述的一种适用于高速高精加工的椭圆弧平滑压缩插补算法,其特征在于所述权重通过下式得到:
<mrow>
<msub>
<mi>w</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mi>P</mi>
<mo>-</mo>
<msub>
<mi>P</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<mi>P</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msup>
<mi>u</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mi>P</mi>
<mo>-</mo>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<mi>P</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>2</mn>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>|</mo>
<mrow>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<mi>P</mi>
</mrow>
<mo>|</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mrow>
其中,w1为权重,P0为首型值点、P2为末型值点、P1为控制点、P为指令点,、u为P0Q与QP2的比值,Q为以P1为投影中心,直线段[P0 P2]投影到二次有理Bézier拟合曲线的投影点。
4.根据权利要求2所述的一种适用于高速高精加工的椭圆弧平滑压缩插补算法,其特征在于所述根据权重通过求取平均值得到二次有理Bézier的拟合曲线包括以下步骤:
4-1)根据权重的求取公式得到求取二次有理Bézier曲线对应的权值wk及相应的肩点坐标sk=wk/(1+wk),k=i+1,...,j-1;
4-2)通过对sk求平均值得到平均肩点坐标s和平均权值w=s/(1-s);
4-3)根据P0,P1,P2和w确定型值点Qi与Qj之间二次有理Bézier的拟合曲线。
5.根据权利要求1所述的一种适用于高速高精加工的椭圆弧平滑压缩插补算法,其特征在于所述根据曲线特征识别椭圆弧及其几何形式转换包括以下步骤:
5-1)当|P0P1|!=|P1P2|且0<w1<1时,则该曲线为椭圆弧;
5-2)通过以下步骤得到椭圆弧的几何信息:
Ci(u)为加工路径中某段的标准型二次有理Bézier曲线,且其对应于一个椭圆弧段,根据Ci(u)可得到该椭圆弧段的起点坐标p_start和终点坐标p_end;该椭圆弧段的中心坐标为,
P1+ε(S+T)
S=P0-P1,T=P2-P1,
该椭圆弧段的长、短半径分别为,
<mrow>
<msub>
<mi>r</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<msqrt>
<mfrac>
<mi>&epsiv;</mi>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
</mfrac>
</msqrt>
<mo>,</mo>
<msub>
<mi>r</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<msqrt>
<mfrac>
<mi>&epsiv;</mi>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
</mfrac>
</msqrt>
</mrow>
假设λ2>λ1>0,且是如下二次方程的根,
2δλ2-(kη+4β)λ+2(k-1)=0
δ=|S×T|2,η=|S-T|2,β=S·T
该椭圆弧段长轴上的两个点为,
Q1=P1+(ε+r1x0)S+(ε+r1y0)T
Q2=P1+(ε-r1x0)S+(ε-r1y0)T
根据Q1、Q2,可以得到长轴的斜率kl,和长轴与x轴正半轴的夹角d_kl;以长轴为x′轴,短轴为y′轴建立局部坐标系,利用Q1、Q2、中心坐标、椭圆弧段起点坐标和椭圆弧段终点坐标,求得椭圆弧段起点与中心连线、椭圆弧段终点与中心连线分别相对于局部坐标系x′正半轴的起始角度d_start、终点角度d_end,以及椭圆弧的方向turn。
6.根据权利要求5所述的一种适用于高速高精加工的椭圆弧平滑压缩插补算法,其特征在于当根据曲线特征识别椭圆弧及其几何形式转换结束后,得到一个椭圆弧段数组Ellipse_Arcs[],数组中每个数据的结构如下:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610978793.3A CN108073138B (zh) | 2016-11-08 | 2016-11-08 | 适用于高速高精加工的椭圆弧平滑压缩插补算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610978793.3A CN108073138B (zh) | 2016-11-08 | 2016-11-08 | 适用于高速高精加工的椭圆弧平滑压缩插补算法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108073138A true CN108073138A (zh) | 2018-05-25 |
CN108073138B CN108073138B (zh) | 2020-08-11 |
Family
ID=62153978
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610978793.3A Active CN108073138B (zh) | 2016-11-08 | 2016-11-08 | 适用于高速高精加工的椭圆弧平滑压缩插补算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108073138B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108908350A (zh) * | 2018-07-16 | 2018-11-30 | 苏州大学张家港工业技术研究院 | 一种交互式机器人的控制***及其控制方法 |
CN109799792A (zh) * | 2019-01-08 | 2019-05-24 | 华侨大学 | 一种面向高速数控加工微小线段平滑过渡的方法 |
WO2020133049A1 (zh) * | 2018-12-27 | 2020-07-02 | 深圳配天智能技术研究院有限公司 | 数控加工方法、数控机床及计算机存储介质 |
CN111552288A (zh) * | 2020-04-28 | 2020-08-18 | 西南交通大学 | 一种移动机器人路径平滑方法 |
CN111610751A (zh) * | 2020-05-21 | 2020-09-01 | 天津工业大学 | 过点集nurbs插值曲线的插值误差多次细分迭代计算方法 |
CN113644917A (zh) * | 2021-06-25 | 2021-11-12 | 河北省科学院应用数学研究所 | 用于脑电、心电数据压缩的方法、终端及存储介质 |
CN113835397A (zh) * | 2021-06-30 | 2021-12-24 | 浙江大学 | 基于b样条曲线和路径积分的线性数控加工路径平滑方法 |
CN113970305A (zh) * | 2021-06-16 | 2022-01-25 | 广西大学 | 一种通过挠度测量压杆轴向位移的方法 |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4243924A (en) * | 1977-06-28 | 1981-01-06 | Oki Electric Industry Co., Ltd. | System for interpolating an arc for a numerical control system |
JPS59229618A (ja) * | 1983-06-10 | 1984-12-24 | Nec Corp | 運動制御装置 |
JPH06150013A (ja) * | 1992-11-11 | 1994-05-31 | Tokyo Electric Co Ltd | 楕円弧発生装置 |
JP2003331296A (ja) * | 2002-05-10 | 2003-11-21 | Ricoh Co Ltd | 円弧描画方法 |
CN101140460A (zh) * | 2006-09-06 | 2008-03-12 | 北京数码大方科技有限公司 | 用于驱动数控机床进行空间圆弧插补加工的方法和装置 |
CN101539769A (zh) * | 2009-04-28 | 2009-09-23 | 中国科学院数学与***科学研究院 | 基于二次b样条曲线对g01代码的拟合及插补方法 |
CN101751006A (zh) * | 2008-12-18 | 2010-06-23 | 陈学恭 | 一种插补方法 |
CN101907876A (zh) * | 2010-05-28 | 2010-12-08 | 沈阳高精数控技术有限公司 | 适用于数控装置的指令点整形压缩插补方法 |
CN102109834A (zh) * | 2011-01-10 | 2011-06-29 | 浙江工业大学 | 基于nurbs曲线插补方法的运动控制*** |
CN102147599A (zh) * | 2011-03-30 | 2011-08-10 | 华南理工大学 | 椭圆弧插补方法 |
CN103699056A (zh) * | 2013-12-02 | 2014-04-02 | 嘉兴学院 | 高速高精度数控加工的小线段实时平滑过渡插补方法 |
CN103809521A (zh) * | 2012-11-14 | 2014-05-21 | 中国科学院沈阳计算技术研究所有限公司 | 基于弦截法的样条曲线插补方法 |
CN104268330A (zh) * | 2014-09-19 | 2015-01-07 | 华南理工大学 | 一种复杂曲线的Bezier拟合方法 |
CN104597847A (zh) * | 2013-10-31 | 2015-05-06 | 中国科学院沈阳计算技术研究所有限公司 | 基于Akima样条曲线拟合的前瞻插补方法 |
-
2016
- 2016-11-08 CN CN201610978793.3A patent/CN108073138B/zh active Active
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4243924A (en) * | 1977-06-28 | 1981-01-06 | Oki Electric Industry Co., Ltd. | System for interpolating an arc for a numerical control system |
JPS59229618A (ja) * | 1983-06-10 | 1984-12-24 | Nec Corp | 運動制御装置 |
JPH06150013A (ja) * | 1992-11-11 | 1994-05-31 | Tokyo Electric Co Ltd | 楕円弧発生装置 |
JP2003331296A (ja) * | 2002-05-10 | 2003-11-21 | Ricoh Co Ltd | 円弧描画方法 |
CN101140460A (zh) * | 2006-09-06 | 2008-03-12 | 北京数码大方科技有限公司 | 用于驱动数控机床进行空间圆弧插补加工的方法和装置 |
CN101751006A (zh) * | 2008-12-18 | 2010-06-23 | 陈学恭 | 一种插补方法 |
CN101539769A (zh) * | 2009-04-28 | 2009-09-23 | 中国科学院数学与***科学研究院 | 基于二次b样条曲线对g01代码的拟合及插补方法 |
CN101907876A (zh) * | 2010-05-28 | 2010-12-08 | 沈阳高精数控技术有限公司 | 适用于数控装置的指令点整形压缩插补方法 |
CN102109834A (zh) * | 2011-01-10 | 2011-06-29 | 浙江工业大学 | 基于nurbs曲线插补方法的运动控制*** |
CN102147599A (zh) * | 2011-03-30 | 2011-08-10 | 华南理工大学 | 椭圆弧插补方法 |
CN103809521A (zh) * | 2012-11-14 | 2014-05-21 | 中国科学院沈阳计算技术研究所有限公司 | 基于弦截法的样条曲线插补方法 |
CN104597847A (zh) * | 2013-10-31 | 2015-05-06 | 中国科学院沈阳计算技术研究所有限公司 | 基于Akima样条曲线拟合的前瞻插补方法 |
CN103699056A (zh) * | 2013-12-02 | 2014-04-02 | 嘉兴学院 | 高速高精度数控加工的小线段实时平滑过渡插补方法 |
CN104268330A (zh) * | 2014-09-19 | 2015-01-07 | 华南理工大学 | 一种复杂曲线的Bezier拟合方法 |
Non-Patent Citations (2)
Title |
---|
单东日: "《基于二次代数样条差分插补的平面参数曲线数控加工》", 《制造技术与机床》 * |
张晓辉等: "《数控加工中的平滑压缩插补算法研究》", 《机械工程学报》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108908350A (zh) * | 2018-07-16 | 2018-11-30 | 苏州大学张家港工业技术研究院 | 一种交互式机器人的控制***及其控制方法 |
WO2020133049A1 (zh) * | 2018-12-27 | 2020-07-02 | 深圳配天智能技术研究院有限公司 | 数控加工方法、数控机床及计算机存储介质 |
CN111684374A (zh) * | 2018-12-27 | 2020-09-18 | 深圳配天智能技术研究院有限公司 | 数控加工方法、数控机床及计算机存储介质 |
CN111684374B (zh) * | 2018-12-27 | 2022-05-31 | 深圳配天智能技术研究院有限公司 | 数控加工方法、数控机床及计算机存储介质 |
CN109799792A (zh) * | 2019-01-08 | 2019-05-24 | 华侨大学 | 一种面向高速数控加工微小线段平滑过渡的方法 |
CN111552288A (zh) * | 2020-04-28 | 2020-08-18 | 西南交通大学 | 一种移动机器人路径平滑方法 |
CN111610751A (zh) * | 2020-05-21 | 2020-09-01 | 天津工业大学 | 过点集nurbs插值曲线的插值误差多次细分迭代计算方法 |
CN113970305A (zh) * | 2021-06-16 | 2022-01-25 | 广西大学 | 一种通过挠度测量压杆轴向位移的方法 |
CN113644917A (zh) * | 2021-06-25 | 2021-11-12 | 河北省科学院应用数学研究所 | 用于脑电、心电数据压缩的方法、终端及存储介质 |
CN113835397A (zh) * | 2021-06-30 | 2021-12-24 | 浙江大学 | 基于b样条曲线和路径积分的线性数控加工路径平滑方法 |
CN113835397B (zh) * | 2021-06-30 | 2024-02-09 | 浙江大学 | 基于b样条曲线和路径积分的线性数控加工路径平滑方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108073138B (zh) | 2020-08-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108073138A (zh) | 适用于高速高精加工的椭圆弧平滑压缩插补算法 | |
CN111061213B (zh) | 一种基于Bezier曲线转角平滑过渡算法的加工方法 | |
CN108062073B (zh) | 一种用于高质量加工的圆弧平滑压缩插补方法 | |
Lai et al. | On the development of a parametric interpolator with confined chord error, feedrate, acceleration and jerk | |
CN106774153B (zh) | 一种s形球头铣刀的轮廓度补偿方法 | |
CN113467384B (zh) | 一种应用于五轴数控机床的拐角过渡方法 | |
CN109214032B (zh) | 一种空心叶片的自适应加工方法 | |
CN114545863B (zh) | 一种基于b样条曲线拟合的数控加工的轨迹平滑方法 | |
CN103149879B (zh) | 一种基于弧长的数控***椭圆插补方法 | |
CN111310106B (zh) | 一种基于逐步逼近原始轮廓的切割轮廓拟合优化方法 | |
CN112284324B (zh) | 一种涡轮导向叶片几何喉部面积的求解方法 | |
CN115202291A (zh) | 一种基于椭圆弧拟合的nurbs曲线插补方法 | |
CN112486101A (zh) | Nurbs曲线自适应前瞻插补方法 | |
CN106844966B (zh) | 一种螺旋桨叶面叶背精确建模方法 | |
CN108196511B (zh) | 一种螺旋刀具车削轨迹的误差预测方法 | |
CN113504764B (zh) | 基于位置矢量加权积分的连续线段数控加工路径平滑方法 | |
CN112051803B (zh) | 一种基于空间平面法向量的小线段光顺方法 | |
CN111610751B (zh) | 过点集nurbs插值曲线的插值误差多次细分迭代计算方法 | |
TW200821784A (en) | Processing path method of converting straight line to NURBS curve of computer numerical control | |
CN109697272B (zh) | 一种简易的二次b样条曲线拟合方法 | |
CN116484635A (zh) | 基于改进圆弧拟合算法的正弦弧线拟合局部光顺方法 | |
Wu et al. | Adaptive sampling method for thin-walled parts based on on-machine measurement | |
CN108628255B (zh) | 一种指令点纠正平滑处理方法 | |
CN113467376B (zh) | 一种面向多加工场景的多轴轨迹压缩方法 | |
CN113835397B (zh) | 基于b样条曲线和路径积分的线性数控加工路径平滑方法 |
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 |
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. |
|
CP01 | Change in the name or title of a patent holder |