CN113741341A - 一种严格定义下切削加工的零件轮廓误差预测方法和系统 - Google Patents

一种严格定义下切削加工的零件轮廓误差预测方法和系统 Download PDF

Info

Publication number
CN113741341A
CN113741341A CN202110902889.2A CN202110902889A CN113741341A CN 113741341 A CN113741341 A CN 113741341A CN 202110902889 A CN202110902889 A CN 202110902889A CN 113741341 A CN113741341 A CN 113741341A
Authority
CN
China
Prior art keywords
interpolation
theoretical
actual
contour error
error
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
CN202110902889.2A
Other languages
English (en)
Other versions
CN113741341B (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.)
Xian Jiaotong University
Original Assignee
Xian Jiaotong 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 Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN202110902889.2A priority Critical patent/CN113741341B/zh
Publication of CN113741341A publication Critical patent/CN113741341A/zh
Application granted granted Critical
Publication of CN113741341B publication Critical patent/CN113741341B/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/404Numerical 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 control arrangements for compensation, e.g. for backlash, overshoot, tool offset, tool wear, temperature, machine construction errors, load, inertia
    • 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/35Nc in input of data, input till input file format
    • G05B2219/35408Calculate new position data from actual data to compensate for contour error

Landscapes

  • Engineering & 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

本发明公开了一种严格定义下切削加工的零件轮廓误差预测方法和系统,属于数控加工技术领域。首先给出轮廓误差的严格定义,避免了在大曲率零件计算过程中定义失效的问题;其次采用系统辨识的方法代替复杂的坐标变换,建立进给伺服系统模型;接下来采用局部搜索的方法确定目标插补点,同时采用局部直线插补的方式拟合加工轨迹,并且根据理论插补点与实际插补点的位置关系分别求解轮廓误差,提高了计算效率,适用于复杂零件轮廓误差的预测。本发明方法同时考虑了轮廓误差的预测效率与计算准确性,适用于大曲率零件轮廓误差的预测。

Description

一种严格定义下切削加工的零件轮廓误差预测方法和系统
技术领域
本发明属于数控加工技术领域,涉及一种严格定义下切削加工的零件轮廓误差预测方法和系统。
背景技术
近年在航空航天、汽车、船舶等领域,逐渐使用大量高速高精数控加工中心进行精密、复杂零部件的加工制造。特别是一些关键的零件,逐渐向型面复杂化、结构轻量化和制造精密化发展,其几何构型复杂且难加工,并且对外形协调性要求较高,对多轴联动数控机床的工作性能提出新的更高要求。而在数控加工过程中,由于机床进给轴伺服系统滞后、多轴动态特性不匹配以及外部扰动因素的存在,导致实际加工轨迹与理论轨迹产生偏差,即轮廓误差。因此,研究轮廓误差预测方法,对提高数控系统跟踪精度,实现高精密数控加工具有重要意义。
在轮廓误差研究方面,专利号CN108803487A公布了一种零件侧铣表面的点位轮廓误差预测方法,重构加工表面,计算实际刀位点到理想加工表面的距离,得到点位轮廓误差;专利号CN109240214A公布了一种面向多轴数控加工的轮廓误差估计与可视化方法,将加工轨迹划分为直线与曲线,曲线轨迹采用线性样条插值为直线,计算各刀位点距离理想轨迹的距离;专利号CN106843146A公布了一种自适应变增益轮廓误差补偿方法,根据数控参数曲线插补器所生成的自由曲线与实际刀位点之间的几何位置关系,利用切向误差逆推策略快速寻找曲线轮廓上距离实际刀位点最近的垂足点,进而实现轮廓误差矢量的高精度估计;专利号CN110262394A公布了一种数控加工中轮廓误差的补偿方法,通过实时采集指令轮廓和实时响应轮廓以估计当前时刻的轮廓误差;李学伟等(李学伟,赵万华,卢秉恒.轨迹误差建模的多轴联动机床轮廓误差补偿技术[J].西安交通大学学报,2012,46(03):47-52.)计算实际点到理论曲线上最近的三个点组成的两条直线的最短距离,选取其中的最小值作为轮廓误差。
综上所述,目前有关轮廓误差预测的研究基本思路为:根据理想插补点数据生成加工轨迹,计算与理想插补点对应的实际位置点坐标,求解实际位置点到理想加工轨迹的距离作为轮廓误差。这种做法的缺点是轮廓误差定义不严格,当加工轨迹曲率变化较大时,会产生较大偏差。同时在轮廓误差预测算法方面,大部分研究聚焦于提高预测精度而忽略了算法计算复杂、计算时间长的问题,很难具体应用到实时轮廓误差计算的任务中,尤其是在零件比较复杂的情况下。因此,亟需发展一种新的关于轮廓误差计算的方法以解决上述问题。
发明内容
本发明的目的在于克服上述现有技术中,轮廓误差测量方法的误差定义不严格且算法计算过程复杂、耗时较长的缺点,提供一种严格定义下切削加工的零件轮廓误差预测方法和系统。
为了达到上述目的,本发明采用以下技术方案予以实现:
一种严格定义下切削加工的零件轮廓误差预测方法,包括如下步骤:
步骤1)获取机床传递函数,并基于零件加工的理论插补指令位置,得到零件加工的实际插补指令位置;
步骤2)基于零件加工特征,为轮廓误差进行严格定义,基于严格定义后的轮廓误差,获取距离实际插补指令位置最近的理论插补指令位置;
严格定义后的轮廓误差为时间跨度最小的前提下,实际插补指令位置到理想加工轨迹的法向距离;
步骤3)基于步骤2)的理论插补指令位置计算得到零件的理想加工轨迹,基于理想加工轨迹得到零件轮廓误差计算需要的理论插补指令位置;
步骤4)基于实际插补指令位置与零件轮廓误差计算需要的理论插补指令位置,获取实际插补指令位置处对应的轮廓误差;
步骤5)重复步骤1)~步骤4),得到每个实际插补指令位置处的轮廓误差。
优选地,步骤1)中机床传递函数是通过系统辨识的方法对机床伺服系统进行辨识得到的。
优选地,机床传递函数的建立过程具体为:
首先获取能够激励机床各进给轴运动的激励信号;
之后将激励信号转化为激励代码输入机床数控系统中,使机床各进给轴进行激励运动;
在机床各进给轴激励运动的过程中,采集理论插补指令位置、光栅检测位置数据;
基于理论插补指令位置和光栅检测位置数据,经过系统辨识,建立机床传递函数。
优选地,步骤2)中,距离实际插补指令位置最近的理论插补指令位置是通过中心窗口法或系统遍历法得到的。
优选地,步骤2)中,通过中心窗口法获得距离实际插补指令位置最近的理论插补指令位置时,具体操作过程为:
在每一个插补点指定一个移动窗口,分别计算实际插补指令位置与窗口内各插补点的距离,从而确定距离实际插补指令位置最近的理论插补指令位置。
优选地,步骤3)中,理想加工轨迹是通过直线插补的方法得到的。
优选地,步骤4)中,轮廓误差的具体计算过程为:
首先获取距离实际插补指令位置最近的理论插补指令位置,计算理论插补指令位置的方向向量、实际插补指令位置与理论插补指令位置连线的方向向量;
根据两个方向向量数量积的符号,确定轮廓误差的具体计算情况;
作实际指令位置到理论插补指令点的方向向量的垂线,得到垂足的位置坐标;
基于轮廓误差定义,得到轮廓误差矢量,进一步计算得到轮廓误差。
进一步优选地,步骤4)中,轮廓误差的具体计算情况包括以下三种:
情况1:当
Figure BDA0003200408280000041
计算Pa
Figure BDA0003200408280000042
的距离作为轮廓误差值;
情况2:当
Figure BDA0003200408280000043
计算Pa
Figure BDA0003200408280000044
的距离作为轮廓误差值;
情况3:当
Figure BDA0003200408280000045
计算Pa到Pr的距离作为轮廓误差值;
其中,Pa为实际插补指令位置;Pr为距离实际插补指令位置最近的理论插补指令位置;Pr+1为Pr后一个理论插补指令位置;Pr-1为Pr前一个理论插补指令位置。
一种严格定义下切削加工的零件轮廓误差预测系统,包括:
传递函数模块,用于获取机床传递函数,基于机床传递函数和加工的理论插补指令位置,得到零件加工的实际插补指令位置;
理想加工轨迹模块,与传递函数模块相交互,用于基于严格定义后的轮廓误差,获取距离实际插补指令位置最近的理论插补指令位置,并进一步计算得到零件的理想加工轨迹;
轮廓误差计算模块,与理想加工轨迹模块相交互,用于计算实际插补指令位置到理想加工轨迹的法向距离。
与现有技术相比,本发明具有以下有益效果:
本发明公开了一种严格定义下切削加工的零件轮廓误差预测方法,首先给出轮廓误差的严格定义,避免了在大曲率零件计算过程中定义失效的问题;其次采用系统辨识的方法代替复杂的坐标变换,建立进给伺服系统模型;接下来采用局部搜索的方法确定目标插补点,同时采用局部直线插补的方式拟合加工轨迹,并且根据理论插补点与实际插补点的位置关系分别求解轮廓误差,提高了计算效率,适用于复杂零件轮廓误差的预测。本发明方法同时考虑了轮廓误差的预测效率与计算准确性,适用于大曲率零件轮廓误差的预测。
本发明还公开了一种严格定义下切削加工的零件轮廓误差预测系统,是基于上述方法建立的,包括三个模块,分别为:传递函数模块,用于获取机床传递函数,基于机床传递函数和加工的理论插补指令位置,得到零件加工的实际插补指令位置;理想加工轨迹模块,与传递函数模块相交互,用于基于严格定义后的轮廓误差,获取距离实际插补指令位置最近的理论插补指令位置,并进一步计算得到零件的理想加工轨迹;轮廓误差计算模块,与理想加工轨迹模块相交互,用于计算实际插补指令位置到理想加工轨迹的法向距离。
附图说明
图1为本发明零件轮廓误差预测方法的流程图;
图2为本发明机床伺服进给系统模型;
图3为本发明轮廓误差严格定义示意图;
图4为本发明距离实际指令位置最近的理论插补点求解原理;
图5为本发明确定轮廓误差所需理论插补点的示意图;
图6为本发明轮廓误差计算示意图。
具体实施方式
下面结合附图对本发明做进一步详细描述:
实施例1
一种严格定义下切削加工的零件轮廓误差预测方法,包括如下步骤:
步骤1)获取机床传递函数,并基于零件加工的理论插补指令位置,得到零件加工的实际插补指令位置;
步骤2)基于零件加工特征,为轮廓误差进行严格定义;严格定义后的轮廓误差为时间跨度最小的前提下,实际插补指令位置到理想加工轨迹的法向距离;
基于严格定义后的轮廓误差,获取距离实际插补指令位置最近的理论插补指令位置;
步骤3)基于步骤2)的理论插补指令位置计算得到零件的理想加工轨迹,基于理想加工轨迹得到零件轮廓误差计算需要的理论插补指令位置;
步骤4)基于实际插补指令位置与理论插补指令位置,获取实际插补指令位置处对应的轮廓误差;
步骤5)重复步骤1)~步骤4),得到每个实际插补指令位置处的轮廓误差实施例2
一种严格定义下切削加工的零件轮廓误差预测方法,包括如下步骤:
步骤1:辨识机床各进给轴传递函数,将零件加工的理论插补指令输入至传递函数,得到实际插补指令位置。采用系统辨识的方法对机床伺服系统的传递函数进行辨识,辨识步骤为:首先获取激励信号,以充分激励起机床各进给轴运动;其次将激励信号转化为激励G代码输入机床数控系统中,使机床各进给轴激励运动,并采集理论插补指令位置、光栅检测位置数据;最后经过系统辨识,建立机床进给系统的传递函数。
接下来获取零件加工的实际插补指令位置。首先获取零件加工的理论插补指令位置,将理论插补指令位置输入至机床传递函数,即可得到对应的实际插补指令位置。
步骤2:根据零件加工特征,给出轮廓误差的严格定义。轮廓误差一般定义为实际插补指令位置距离理想加工轨迹的偏差,常见的定义有两种,即实际插补指令位置到理想轨迹的法向距离(定义1)以及实际插补指令位置到理想轨迹的最短距离(定义2)。当零件曲率较大时,按照定义1计算,可能会产生多个轮廓误差结果,按照定义2计算,可能会在加工轨迹曲率较大处产生过切,因此,上述定义均不能准确计算轮廓误差。基于此,本文给出轮廓误差严格定义:轮廓误差为实际插补指令位置到理想加工轨迹的法向距离,且时间跨度最小,在严格定义下,计算轮廓误差时可以避免之前定义的问题。
步骤3:求解距离实际插补指令位置Pa最近的理论插补指令位置Pr。为了后续轮廓误差计算,需要获得理想加工轨迹上距离实际插补指令位置Pa最近的理论插补指令位置Pr。假定(Pi,Oi)表示指令轨迹上任一点刀具位姿,Pi表示参考刀尖点位置,Oi表示刀轴方向;(Pa,Oa)表示实际轨迹上与Pi点对应的插补点位姿;Pc、Pc1分别为Pa在理想轨迹上的垂足;Pr为在指定插补周期内,距离Pa最近的插补点;μp描述实际刀尖点位置到刀尖点位置参考轨迹的位置偏差,而μo描述实际刀轴方向到刀轴方向参考轨迹的角度偏差。
采用局部逐点比较的方法,即为每一个插补点Pi指定一个移动窗口k,分别计算实际插补指令位置Pa与窗口内各插补点的距离,从而确定距离实际插补指令位置Pa最近的理论插补指令位置Pr
步骤4:生成理想加工轨迹,确定轮廓误差计算所需的插补点。计算轮廓误差的关键是计算实际插补指令位置Pa到理想轨迹的法向距离,但是考虑计算效率,本方法采用直线插补的方式生成理想加工轨迹,用线段Pr-1Pr作为Pr-1与Pr之间的理想轨迹,线段PrPr+1作为Pr与Pr+1之间的理想轨迹。同时根据实际插补指令位置Pa与Pr-1、Pr、Pr+1的相对距离情况,确定后续轮廓误差计算所需的理论插补指令位置以及理想轨迹。
根据实际插补指令位置Pa与Pr-1、Pr、Pr+1的相对距离情况,可以分为以下三种情况:
情况1:此时
Figure BDA0003200408280000081
计算Pa到
Figure BDA0003200408280000082
的距离,作为轮廓误差值。
情况2:此时
Figure BDA0003200408280000083
计算Pa到
Figure BDA0003200408280000084
的距离,作为轮廓误差值。
情况3:此时
Figure BDA0003200408280000085
计算Pa到Pr的距离,作为轮廓误差值。
步骤5:根据实际插补指令位置Pa与理论插补指令位置信息,求解Pa距离理想加工轨迹的距离。
以步骤4中的情况1为例,轮廓误差计算步骤如下:
设各插补指令点坐标为:Pr-1(xr-1,yr-1,zr-1),Pr(xr,yr,zr),Pr+1(xr+1,yr+1,zr+1);实际插补指令位置点坐标:Pa(xa,ya,za),理论插补指令位置的方向向量为:
Figure BDA0003200408280000086
其中:X=xr-xr-1,Y=yr-yr-1,Z=zr-zr-1
作点Pa到直线
Figure BDA0003200408280000087
的垂线,垂足为Pr'(xr',yr',zr'),则直线
Figure BDA0003200408280000088
的方向矢量:
Figure BDA0003200408280000091
其中:X'=xa-xr',Y=ya-yr',Z=za-zr'
由向量
Figure BDA0003200408280000092
和向量
Figure BDA0003200408280000093
垂直可得:
X·X′+Y·Y′+Z·Z=0
且点Pr’在直线Pr-1Pr上,进而求得Pi'(xi',yi',zi')坐标如下:
Figure BDA0003200408280000094
Figure BDA0003200408280000095
Figure BDA0003200408280000096
其中,Ex=xa-xr,Ey=ya-yr,Ez=za-zr
根据轮廓误差定义,可以得到刀尖位置轮廓误差矢量为:
Figure BDA0003200408280000097
进而可以得到:
Figure BDA0003200408280000098
Figure BDA0003200408280000099
Figure BDA00032004082800000910
此时,刀尖轮廓误差的模长为:
Figure BDA00032004082800000911
步骤6:重复步骤1~5,得到每个实际插补指令位置处的轮廓误差。
实施例3
一种严格定义下切削加工的零件轮廓误差预测方法,如图1所示,包括以下步骤:
步骤1:辨识机床各进给轴传递函数,根据理论插补指令位置预测实际插补指令位置。
机床在加工过程中,由于存在机床结构误差、切削过程以及驱动系统动态特性的影响,导致各进给系统实际位置与指令位置之间存在跟随误差,在多轴联动时便会产生轮廓误差,这里的位置指的是刀位点信息,因此为了对轮廓误差进行预测,首先需要获得零件加工指令的理想刀位点信息与实际刀位点信息。参照图2,加工指令经过机床的伺服控制系统,会输出实际的刀位点信息。基于此,采用系统辨识的方法对机床伺服系统的传递函数进行辨识,辨识步骤为:首先获取激励信号,由于M序列自相关性较好,具有伪随机性,容易产生和复制,因此选择变幅值M序列信号,以充分激励起机床各进给轴运动;其次将激励信号转化为激励G代码输入机床数控系统中,使机床各进给轴激励运动,并采集理论插补指令位置、光栅检测位置数据;最后经过系统辨识,建立各进给轴传递函数,其离散表达形式为:
Figure BDA0003200408280000101
其中:B(z-1)为光栅检测位置;A(z-1)为插补指令位置;bi和a j分别为离散传递函数分子及分母的系数,na和nb分别为离散传递函数分子及分母的阶数。根据机床伺服系统的特点,将传递函数简化为二阶模型,利用Matlab系统辨识工具箱即可得到各进给轴的传递函数。辨识实验所用机床为五轴数控加工中心DMU50,刀柄型号为BT50;刀具为整体式硬质合金刀具,直径20mm,4齿,螺旋角为30°,刀具悬长78mm。辨识得到x轴以及y轴传递函数如下:
Figure BDA0003200408280000111
Figure BDA0003200408280000112
接下来获取实际插补指令位置。首先获取零件加工的理论插补指令位置,将理论插补指令位置输入至机床传递函数中,即可得到实际插补指令位置。
步骤2:根据加工零件特征,给出轮廓误差的严格定义。轮廓误差一般定义为实际插补指令位置距离理想加工轨迹的偏差,常见的定义有两种,即实际插补指令位置到理想轨迹的法向距离(定义1)以及实际插补指令位置到理想轨迹的最短距离(定义2)。参照图3,(Pi,Oi)表示指令轨迹上任一点刀具位姿,Pi表示参考刀尖点位置,Oi表示刀轴方向;(Pa,Oa)表示实际轨迹上与Pi点对应的插补点位姿;Pc、Pc1分别为Pa在理想轨迹上的垂足;Pr为在指定插补周期内,距离Pa最近的插补点;μp描述实际刀尖点位置到刀尖点位置参考轨迹的位置偏差,而μo描述实际刀轴方向到刀轴方向参考轨迹的角度偏差。
此时加工轨迹曲率较大,如果按照定义1计算,存在两个法向距离||PaPc||与||PaPc1||,无法判断轮廓误差的取值;若按照定义2计算,此时的轮廓误差为||Pa Pc1||,则插补点Pi与Pc1之间的轨迹都会被切除,导致工件失效。因此上述定义均不能准确计算轮廓误差。基于此,本文给出轮廓误差严格定义:轮廓误差为实际插补指令位置到理想加工轨迹的法向距离,且时间跨度最小。在严格定义下,Pa点的轮廓误差为||PaPc||,此时避免了上述定义不严格的问题。
步骤3:求解距离实际插补指令位置Pa最近的理论插补指令位置Pr。根据步骤2轮廓误差严格定义,Pa点的轮廓误差为该点到最近插补点轨迹的法向距离,因此需要首先获取距离实际插补指令位置Pa最近的理论插补指令位置Pr
参照图4,采用局部逐点比较的方法,即为每一个插补点Pi指定一个移动窗口k,分别计算实际插补指令位置Pa与窗口内各插补点的距离,从而确定距离实际插补指令位置最近的理论插补指令位置Pr
步骤4:生成理想加工轨迹,确定轮廓误差计算所需的插补点。计算轮廓误差的关键是计算实际插补指令位置Pa到理想轨迹的法向距离,但是考虑计算效率,采用直线插补的方式生成理想加工轨迹,参照图6,用线段Pr-1Pr作为Pr-1与Pr之间的理想轨迹,线段PrPr+1作为Pr与Pr+1之间的理想轨迹。
同时根据实际插补指令位置Pa与Pr-1、Pr、Pr+1的相对距离情况,确定后续轮廓误差计算所需的理想轨迹。
参照图5,根据实际插补指令位置Pa与Pr-1、Pr、Pr+1的相对距离情况,可以分为以下三种情况:
情况1:此时
Figure BDA0003200408280000121
计算Pa到
Figure BDA0003200408280000122
的距离,作为轮廓误差值。
情况2:此时
Figure BDA0003200408280000123
计算Pa到
Figure BDA0003200408280000124
的距离,作为轮廓误差值。
情况3:此时
Figure BDA0003200408280000125
计算Pa到Pr的距离,作为轮廓误差值。
步骤5:根据实际指令位置Pa与插补指令的位置信息,求解Pa距离理想加工轨迹的距离。
以步骤4中的情况1为例,轮廓误差计算步骤如下:
设各插补指令点坐标为:Pr-1(xr-1,yr-1,zr-1),Pr(xr,yr,zr),Pr+1(xr+1,yr+1,zr+1);实际指令位置点坐标:Pa(xa,ya,za),插补指令点的方向向量为:
Figure BDA0003200408280000126
其中:X=xr-xr-1,Y=yr-yr-1,Z=zr-zr-1
作点Pa到直线
Figure BDA0003200408280000131
的垂线,垂足为Pr'(xr',yr',zr'),则直线
Figure BDA0003200408280000132
的方向矢量:
Figure BDA0003200408280000133
其中:X'=xa-xr',Y=ya-yr',Z=za-zr'
由向量
Figure BDA0003200408280000134
和向量
Figure BDA0003200408280000135
垂直可得:
X·X′+Y·Y′+Z·Z=0
且点Pr’在直线Pr-1Pr上,进而求得Pi'(xi',yi',zi')坐标如下:
Figure BDA0003200408280000136
Figure BDA0003200408280000137
Figure BDA0003200408280000138
其中,Ex=xa-xr,Ey=ya-yr,Ez=za-zr
根据轮廓误差定义,可以得到刀尖位置轮廓误差矢量为:
Figure BDA0003200408280000139
进而可以得到:
Figure BDA00032004082800001310
Figure BDA00032004082800001311
Figure BDA00032004082800001312
此时,刀尖轮廓误差的模长为:
Figure BDA00032004082800001313
步骤6:重复步骤1~5,得到每个实际插补指令位置处的轮廓误差。
实施例4
一种严格定义下切削加工的零件轮廓误差预测系统,包括:
传递函数模块,用于获取机床传递函数,基于机床传递函数和加工的理论插补指令位置,得到零件加工的实际插补指令位置;
理想加工轨迹模块,与传递函数模块相交互,用于基于严格定义后的轮廓误差,获取距离实际插补指令位置最近的理论插补指令位置,并进一步计算得到零件的理想加工轨迹;
轮廓误差计算模块,与理想加工轨迹模块相交互,用于计算实际插补指令位置到理想加工轨迹的法向距离。
综上所述,本发明的零件轮廓误差预测方法,首先辨识机床各进给轴的传递函数,将零件加工的理论插补指令位置输入至传递函数,得到实际插补指令位置;接下来给出轮廓误差的严格定义,并利用中心窗口法求解距离实际插补指令位置最近的理论插补指令位置;根据理论插补指令位置拟合得到理想加工轨迹,并计算实际插补指令位置距离理想加工轨迹的偏差,即可得到每个实际插补指令位置处的轮廓误差。本发明同时考虑了轮廓误差的计算准确性与预测效率,可以实现零件轮廓误差的准确快速预估,能够为优化加工参数、保证零件轮廓精度提供理论依据。
以上内容仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明权利要求书的保护范围之内。

Claims (9)

1.一种严格定义下切削加工的零件轮廓误差预测方法,其特征在于,包括如下步骤:
步骤1)获取机床传递函数,并基于零件加工的理论插补指令位置,得到零件加工的实际插补指令位置;
步骤2)基于零件加工特征,为轮廓误差进行严格定义,基于严格定义后的轮廓误差,获取距离实际插补指令位置最近的理论插补指令位置;
严格定义后的轮廓误差为时间跨度最小的前提下,实际插补指令位置到理想加工轨迹的法向距离;
步骤3)基于步骤2)的理论插补指令位置计算得到零件的理想加工轨迹,基于理想加工轨迹得到零件轮廓误差计算需要的理论插补指令位置;
步骤4)基于实际插补指令位置与零件轮廓误差计算需要的理论插补指令位置,获取实际插补指令位置处对应的轮廓误差;
步骤5)重复步骤1)~步骤4),得到每个实际插补指令位置处的轮廓误差。
2.根据权利要求1所述的严格定义下切削加工的零件轮廓误差预测方法,其特征在于,步骤1)中机床传递函数是通过系统辨识的方法对机床伺服系统进行辨识得到的。
3.根据权利要求1所述的严格定义下切削加工的零件轮廓误差预测方法,其特征在于,机床传递函数的建立过程具体为:
首先获取能够激励机床各进给轴运动的激励信号;
之后将激励信号转化为激励代码输入机床数控系统中,使机床各进给轴进行激励运动;
在机床各进给轴激励运动的过程中,采集理论插补指令位置、光栅检测位置数据;
基于理论插补指令位置和光栅检测位置数据,经过系统辨识,建立机床传递函数。
4.根据权利要求1所述的严格定义下切削加工的零件轮廓误差预测方法,其特征在于,步骤2)中,距离实际插补指令位置最近的理论插补指令位置是通过中心窗口法或系统遍历法得到的。
5.根据权利要求4所述的严格定义下切削加工的零件轮廓误差预测方法,其特征在于,步骤2)中,通过中心窗口法获得距离实际插补指令位置最近的理论插补指令位置时,具体操作过程为:
在每一个插补点指定一个移动窗口,分别计算实际插补指令位置与窗口内各插补点的距离,从而确定距离实际插补指令位置最近的理论插补指令位置。
6.根据权利要求1所述的严格定义下切削加工的零件轮廓误差预测方法,其特征在于,步骤3)中,理想加工轨迹是通过直线插补的方法得到的。
7.根据权利要求1所述的严格定义下切削加工的零件轮廓误差预测方法,其特征在于,步骤4)中,轮廓误差的具体计算过程为:
首先获取距离实际插补指令位置最近的理论插补指令位置,计算理论插补指令位置的方向向量、实际插补指令位置与理论插补指令位置连线的方向向量;
根据两个方向向量数量积的符号,确定轮廓误差的具体计算情况;
作实际指令位置到理论插补指令点的方向向量的垂线,得到垂足的位置坐标;
基于轮廓误差定义,得到轮廓误差矢量,进一步计算得到轮廓误差。
8.根据权利要求7所述的严格定义下切削加工的零件轮廓误差预测方法,其特征在于,步骤4)中,轮廓误差的具体计算情况包括以下三种:
情况1:当
Figure FDA0003200408270000031
计算Pa
Figure FDA0003200408270000032
的距离作为轮廓误差值;
情况2:当
Figure FDA0003200408270000033
计算Pa
Figure FDA0003200408270000034
的距离作为轮廓误差值;
情况3:当
Figure FDA0003200408270000035
计算Pa到Pr的距离作为轮廓误差值;
其中,Pa为实际插补指令位置;Pr为距离实际插补指令位置最近的理论插补指令位置;Pr+1为Pr后一个理论插补指令位置;Pr-1为Pr前一个理论插补指令位置。
9.一种严格定义下切削加工的零件轮廓误差预测系统,其特征在于,包括:
传递函数模块,用于获取机床传递函数,基于机床传递函数和加工的理论插补指令位置,得到零件加工的实际插补指令位置;
理想加工轨迹模块,与传递函数模块相交互,用于基于严格定义后的轮廓误差,获取距离实际插补指令位置最近的理论插补指令位置,并进一步计算得到零件的理想加工轨迹;
轮廓误差计算模块,与理想加工轨迹模块相交互,用于计算实际插补指令位置到理想加工轨迹的法向距离。
CN202110902889.2A 2021-08-06 2021-08-06 一种严格定义下切削加工的零件轮廓误差预测方法和系统 Active CN113741341B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110902889.2A CN113741341B (zh) 2021-08-06 2021-08-06 一种严格定义下切削加工的零件轮廓误差预测方法和系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110902889.2A CN113741341B (zh) 2021-08-06 2021-08-06 一种严格定义下切削加工的零件轮廓误差预测方法和系统

Publications (2)

Publication Number Publication Date
CN113741341A true CN113741341A (zh) 2021-12-03
CN113741341B CN113741341B (zh) 2023-03-31

Family

ID=78730379

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110902889.2A Active CN113741341B (zh) 2021-08-06 2021-08-06 一种严格定义下切削加工的零件轮廓误差预测方法和系统

Country Status (1)

Country Link
CN (1) CN113741341B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117908464A (zh) * 2024-03-20 2024-04-19 山东大学 一种机床轮廓误差预测方法、系统、装置及可读存储介质

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030081722A1 (en) * 2001-08-27 2003-05-01 Nikon Corporation Multilayer-film mirrors for use in extreme UV optical systems, and methods for manufacturing such mirrors exhibiting improved wave aberrations
KR20050103627A (ko) * 2004-04-27 2005-11-01 김두진 수치제어밀링의 윤곽오차 계산 및 제어 방법
CN101342664A (zh) * 2008-08-20 2009-01-14 天津大学 一种复杂空间型面在机质量检测系统
CN103197598A (zh) * 2013-03-13 2013-07-10 西北工业大学 基于非参数多元统计过程控制的叶身型面加工误差控制方法
CN104483897A (zh) * 2014-08-27 2015-04-01 渤海大学 一种直驱龙门式运动平台轮廓控制装置及方法
CN109116808A (zh) * 2018-08-23 2019-01-01 西安交通大学 一种基于数控机床大数据的联动轨迹轮廓精度预测方法
CN110161966A (zh) * 2019-07-04 2019-08-23 合肥工业大学 一种多轴联动轮廓误差的实时检测方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030081722A1 (en) * 2001-08-27 2003-05-01 Nikon Corporation Multilayer-film mirrors for use in extreme UV optical systems, and methods for manufacturing such mirrors exhibiting improved wave aberrations
KR20050103627A (ko) * 2004-04-27 2005-11-01 김두진 수치제어밀링의 윤곽오차 계산 및 제어 방법
CN101342664A (zh) * 2008-08-20 2009-01-14 天津大学 一种复杂空间型面在机质量检测系统
CN103197598A (zh) * 2013-03-13 2013-07-10 西北工业大学 基于非参数多元统计过程控制的叶身型面加工误差控制方法
CN104483897A (zh) * 2014-08-27 2015-04-01 渤海大学 一种直驱龙门式运动平台轮廓控制装置及方法
CN109116808A (zh) * 2018-08-23 2019-01-01 西安交通大学 一种基于数控机床大数据的联动轨迹轮廓精度预测方法
CN110161966A (zh) * 2019-07-04 2019-08-23 合肥工业大学 一种多轴联动轮廓误差的实时检测方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
吕盾等: "五轴机床加工零件轮廓误差预测方法", 《西安交通大学学报》 *
李祥飞等: "面向伺服动态特性匹配的轮廓误差补偿控制研究", 《机械工程学报》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117908464A (zh) * 2024-03-20 2024-04-19 山东大学 一种机床轮廓误差预测方法、系统、装置及可读存储介质

Also Published As

Publication number Publication date
CN113741341B (zh) 2023-03-31

Similar Documents

Publication Publication Date Title
Yuwen et al. Path, feedrate and trajectory planning for free-form surface machining: A state-of-the-art review
CN112698625B (zh) 一种面向联动轨迹误差预测的五轴机床数字孪生建模方法
CN101493687B (zh) 实时前瞻全程加减速控制的nurbs曲线自适应分段插补方法
CN102591257B (zh) 面向参数曲线刀具轨迹的数控系统轮廓误差控制方法
CN103064344B (zh) 一种基于nurbs曲线插补的速度平滑控制方法
CN101539769B (zh) 基于二次b样条曲线对g01代码的拟合及插补方法
CN110262394B (zh) 一种数控加工中轮廓误差的补偿方法
Li et al. Dual sliding mode contouring control with high accuracy contour error estimation for five-axis CNC machine tools
CN105242637A (zh) 一种航空薄壁叶片补偿加工方法
CN108803487B (zh) 一种零件侧铣表面的点位轮廓误差预测方法
CN111427308B (zh) 一种用于数控平台轨迹规划的误差补偿综合控制方法
Feng et al. A novel feature-guided trajectory generation method based on point cloud for robotic grinding of freeform welds
CN110531709B (zh) 一种分析曲面零件加工误差与进给率关系的方法
CN113741341B (zh) 一种严格定义下切削加工的零件轮廓误差预测方法和系统
Liping et al. Prediction of cutting forces in flank milling of parts with non-developable ruled surfaces
Ye et al. Geometric parameter optimization in multi-axis machining
CN114839921A (zh) 一种基于数据驱动的五轴轮廓控制方法
CN114019910A (zh) 一种小线段刀具轨迹实时全局光顺方法
CN104865898A (zh) 数控机床中处理微小路径段的方法及设备
CN110161966B (zh) 一种多轴联动轮廓误差的实时检测方法
WO2024124794A1 (zh) 一种五轴联动同步刀具路径插补方法和系统
CN112731865A (zh) 一种基于插补数据的多轴加工轮廓误差预补偿方法
Lu et al. Collision-free and smooth path generation method of five-axis ball-end milling considering effective tool diameter change
KR100871456B1 (ko) 5축 밀링가공 시간 예측방법
Zhang et al. Intelligent machining of complex aviation components

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