CN109240205A - 一种数控机床非线性误差补偿方法 - Google Patents

一种数控机床非线性误差补偿方法 Download PDF

Info

Publication number
CN109240205A
CN109240205A CN201811179059.6A CN201811179059A CN109240205A CN 109240205 A CN109240205 A CN 109240205A CN 201811179059 A CN201811179059 A CN 201811179059A CN 109240205 A CN109240205 A CN 109240205A
Authority
CN
China
Prior art keywords
error compensation
nonlinear error
control system
digital control
nonlinear
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
CN201811179059.6A
Other languages
English (en)
Other versions
CN109240205B (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.)
Hunan Institute of Technology
Original Assignee
Hunan Institute of Technology
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 Hunan Institute of Technology filed Critical Hunan Institute of Technology
Priority to CN201811179059.6A priority Critical patent/CN109240205B/zh
Publication of CN109240205A publication Critical patent/CN109240205A/zh
Application granted granted Critical
Publication of CN109240205B publication Critical patent/CN109240205B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/404—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 control arrangements for compensation, e.g. for backlash, overshoot, tool offset, tool wear, temperature, machine construction errors, load, inertia
    • 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
    • G05B2219/00—Program-control systems
    • G05B2219/30—Nc systems
    • G05B2219/33—Director till display
    • G05B2219/33133—For each action define function for compensation, enter parameters

Landscapes

  • Engineering & Computer Science (AREA)
  • Human Computer Interaction (AREA)
  • Manufacturing & Machinery (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

一种数控机床非线性误差补偿方法，是：在PC端生成非线性误差补偿激光测量数据报表，再通过管理软件下载至数控系统，数控系统控制端打开非线性误差补偿报表，自动设定测量轴非线性误差补偿参数，最终完成对光栅尺反馈坐标进行非线性误差补偿与坐标显示输出非线性误差逆向补偿。本发明使用自动变步长线性回归方法可以根据样本数据自动改变步长大小，快速得出满足相关系数的线性回归方程，实现数控机床各轴的非线性误差快速补偿。上述非线性误差补偿方法简单实用，可有效降低系统的设计难度，提高加工精度。

Description

一种数控机床非线性误差补偿方法
技术领域
本发明涉及数控技术与自动化控制技术领域，尤其指一种数控机床非线性误差补偿方法。
背景技术
在经济全球化和信息技术革命的影响下，国际制造业正在发生深刻变化，制造业的规模和水平已成为衡量一个国家综合实力的重要标志。数控(numerical control，NC)机床，特别是高档数控机床是国际装备制造业竞争的热点领域，加强数控系统技术领域的基础研究和共性关键问题的攻关，已成为装备制造业发展的当务之急。
数控机床加工过程中，产生误差是在所难免的，被加工零件与数控机床之间存在误差是必然的现象。这种一定存在的误差称为原始误差。所以，要想提高数控机床零件加工的精度，控制数控机床的原始误差是比较重要的对策，常见的测量方法是采用双频激光干涉仪检测。检测完成后一般需要对数控机床的定位精度进行误差补偿。而机床的定位误差一般是非线性的，造成非线性关系输出的原因很多，除了传感器自身的因素外，还有很多外界因素(例如漂移)。这些外界因素很复杂，不能直接计算得到。所以要进行非线性补偿，排除其他影响，才能得到合理的输出。
发明内容
本发明所要解决的技术问题是，提供一种数控机床非线性误差补偿方法，所述数控机床非线性误差补偿方法简单实用，可有效降低系统的设计难度，提高加工精度。
本发明解决其技术问题所采用的技术方案是：一种数控机床非线性误差补偿方法，是：在PC端生成非线性误差补偿激光测量数据报表，再通过管理软件下载至数控系统，数控系统控制端打开非线性误差补偿报表，自动设定测量轴非线性误差补偿参数，最终完成对光栅尺反馈坐标进行非线性误差补偿与坐标显示输出非线性误差逆向补偿。
所述的数控机床非线性误差补偿方法，具体包括以下步骤：
S101：正确设置数控系统非线性误差补偿参数(设置每个轴需要补偿的数据个数MAXL)，选择测量轴，开启非线性误差补偿功能；
S102：打开激光测量系统，激光测量准备就绪，进行非线性误差补偿激光数据处理；
S103：在激光测量系统中运行测量轴对应的通用数控程序；
S104：PC端的激光测量系统生成非线性误差补偿激光测量数据报表；
S105：下载非线性误差补偿数据报表至数控系统；
S106：数控系统控制端打开非线性误差补偿报表，自动设定测量轴非线性误差补偿参数；
S107：数控系统通过自动变步长线性回归的方法自动进行误差补偿，即对光栅尺反馈坐标进行非线性误差补偿与坐标显示输出非线性误差逆向补偿。
进一步，所述对光栅尺反馈坐标进行非线性误差补偿：数控系统读取各轴未补偿的机床坐标值(光栅尺反馈值)的数据报表，通过自动变步长线性回归方法计算输出补偿后的机床坐标值，对应激光测量值；
所述坐标显示输出非线性误差逆向补偿：数控系统读取当前各轴绝对坐标位置值(相当于理论激光测量正确值)，通过回归方程计算输出显示坐标的非线性误差逆向补偿值，相当于光栅尺反馈值。
进一步，所述自动变步长线性回归方法，包括以下步骤：
1)假设线性回归方程y＝ax+b；开始选择步长数据为n个(xC，yC)；
2)基于最小二乘法拟合直线，得出a,b的参数化方程如下：
其中：
为x0～xn相加所得总和的平均值；
为y0～yn相加所得总和的平均值；
为x0～xn的平方和的平均值；
为xi与yi的乘积总和的平均值(i＝[0,n])；
3)计算线性度量pearson系数；
其中：
Lxx为x的方差；
Lyy为y的方差；
Lxy为x与y的协方差；
4)自动变步长
当计算pearson(相关系数)系数小于系统设定值时，删除测量点第n个(即最后一个样本数据)点(x，y)；重新计算pearson(相关系数)系数，循环直至pearson(相关系数)系数达到要求；
5)同理采用自动变步长方法完成后面所有数据(数据为PC端生成非线性误差补偿激光测量数据报表)的线性回归方程。
本发明使用自动变步长线性回归方法可以根据样本数据自动改变步长大小，快速得出满足相关系数的线性回归方程，实现数控机床各轴的非线性误差快速补偿。
附图说明
图1为本发明中的数控机床非线性误差补偿方法的流程图；
图2为本发明中非线性误差补偿方法示意图。
具体实施方式
为了便于本领域技术人员的理解，下面结合实施例和附图对本发明作进一步的说明，实施方式提及的内容并非对本发明的限定。
参照图1，数控机床非线性误差补偿方法的流程图，包括以下步骤：
S101：正确设置数控系统非线性误差补偿参数(设置每个轴需要补偿的数据个数MAXL)，选择测量轴，开启非线性误差补偿功能；
S102：打开激光测量系统，激光测量准备就绪，进行非线性误差补偿激光数据处理；
S103：在激光测量系统中运行测量轴对应的通用数控程序；
S104：PC端的激光测量系统生成非线性误差补偿激光测量数据报表；
S105：下载非线性误差补偿数据报表至数控系统；
S106：数控系统控制端打开非线性误差补偿报表，自动设定测量轴非线性误差补偿参数；
S107：数控系统通过自动变步长线性回归的方法自动进行误差补偿，即对光栅尺反馈坐标进行非线性误差补偿与坐标显示输出非线性误差逆向补偿。
具体地，非线性误差补偿功能需要使用的软件及文件：PC端激光测量系统非线性误差补偿参数处理软件(包括正向和反向)(是现有通用技术，如雷尼绍(Renishaw)激光测量系统)；PC端发送至数控系统的管理软件(是现有通用技术)；测量轴专用的非线性误差补偿参数运行NC文件是现有通用技术)。
激光测量系统需要规范的测量轴非线性误差补偿参数运行数控程序文件是CAM(计算机辅助制造)软件产生的数控加工程序，由G代码组成，是现有通用技术。
基于上述的方法在PC端生成非线性误差补偿激光测量数据报表，如表1所示，再通过管理软件下载至数控系统，数控系统控制端打开非线性误差补偿报表，自动设定测量轴非线性误差补偿参数，最终完成对光栅尺反馈坐标进行非线性误差补偿与坐标显示输出非线性误差逆向补偿。
具体地，非线性误差补偿是基于自动变步长线性回归的方法实现的。自动变步长线性回归的方法如下，如表1所示数据。
1.1)假设线性回归方程y＝ax+b；开始选择步长数据为10个(xC，yC)；
式中，y是与自变量取值x相对应的因变量均值的估计；a和b分别为总体回归方程参数的估计量，b是样本回归方程的常数项，也就是样本回归直线在y轴上的截距，表示除自变量x以外的其他因素对因变量y的平均影响量；a是样本回归系数，也即样本回归直线的斜率，表示自变量x每增加一个单位时因变量y的平均增加量；
1.2)基于最小二乘法拟合直线，得出a，b的参数化方程如下：
即：
y＝1.00020363636363637x-2.36363636363665；
其中，表1为PC端生成非线性误差补偿激光测量数据报表；
表1
为XL[0]～XL[9]的平均值；
为XR[0]～XR[9]的平均值；
为XL[0]～XL[9]的平方和的平均值；
为XL[i]与XL[i]的乘积总和的平均值(i＝[0,9])；
1.3)计算线性度量pearson(相关系数)系数
其中：
Lxx为x的方差；
Lyy为y的方差；
Lxy为x与y的协方差；
1.4)自动变步长
取pearson(相关系数)系数系统设定值为0.98，计算pearson(相关系数)系数小于系统设定值，删除所取样本最后一个数据点(第10个点(90000,90020))；重新计算pearson(相关系数)系数。
重新计算pearson(相关系数)系数如下。
2.1)删除所取样本最后一个数据点(第10个点(90000,90020))后，步长数据为9个(xC，yC)；
2.2)基于最小二乘法拟合直线，得出a,b的参数化方程如下：
即：
y＝1.00017x-1.46667；
其中：
为XL[0]～XL[9]的平均值；
为XR[0]～XR[9]的平均值；
为XL[0]～XL[9]的平方和的平均值；
为XL[i]与XL[i]的乘积总和的平均值(i＝[0,9])；
2.3)计算线性度量pearson(相关系数)系数
其中：
2.4)自动变步长
pearson(相关系数)系数系统设定阈值为0.98，计算pearson(相关系数)系数为0.999999963，大于系统设定值，pearson(相关系数)系数计算值达到要求。即线性回归的步长自动变为9个点的值。
同理采用自动变步长方法完成后面所有样本数据的线性回归方程。
如图2所示，数控系统读取各轴未补偿的机床坐标值(光栅尺反馈值)的数据报表，通过上一步自动变步长计算出的线性回归方程，计算输出补偿后的机床坐标值，对应激光测量值。
具体地，以X轴为例，设定XL:光栅尺反馈值；XR:非线性补偿后的值；ECPT:光栅尺非线性补偿计算数组位置指针。
程序流程如下：
1)首先根据光栅尺原始数据判断当前位置处于哪个补偿区域。
2)ECPT＝0 THEN XR＝XL-XL[0]+XR[0]
注释：当光栅尺非线性补偿计算数组位置指针为0时，XR＝XL-XL[0]+XR[0]。XR[0]表示数组XR第一个数据，即PC端生成非线性误差补偿激光测量数据报表中当前绝对坐标位置值第一个数据；XL[0]表示数组XL第一个数据，即PC端生成非线性误差补偿激光测量数据报表中光栅尺反馈值第一个数据。
3)ECPT＝MAXL+1 THEN XR＝XL-XL[MAXL]+XR[MAXL]
注释：当光栅尺非线性补偿计算数组位置指针为最大值时，XR＝XL-XL[MAXL]+XR[MAXL]。XR[MAXL]表示数组XR最后一个数据，即PC端生成非线性误差补偿激光测量数据报表中当前绝对坐标位置值最后一个数据；XL[MAXL]表示数组XL最后一个数据，即PC端生成非线性误差补偿激光测量数据报表中光栅尺反馈值的最后一个数据。
4)ECPT∈(0,MAXL)THEN代入方程y＝ax+b
XR＝a*XL+b
注释：当光栅尺非线性补偿计算数组位置指针为最小值与最大值之间时，代入线性回归方程求得XR(a,b系数已求出，求解简单二元一次方程即可)。
以表1中的数据为例，当光栅尺反馈值为55000时，代入方程
XR＝1.00017*XL-1.46667
＝1.00017*55000-1.4667
＝55007.8833＝55008(取整)
计算输出补偿后的机床坐标值应为55008。即光栅尺反馈值为55000时，通过数控系统非线性误差补偿方法进行误差补偿后实际机床坐标值为55008。
如图2所示，数控系统读取当前各轴绝对坐标位置值(相当于理论激光测量正确值)，通过计算输出显示坐标的非线性误差逆向补偿值，相当于光栅尺反馈值。
具体地，以X轴为例，设定XR：当前绝对坐标位置值(相当于理论激光测量正确值)；XL：计算X轴显示坐标的非线性误差逆向补偿值，相当于光栅尺反馈值；POCPT：光栅尺非线性补偿计算数组位置指针。
程序流程如下：
1)首先根据当前插补点机床坐标位置值XR判断当前位置处于哪个补偿区域。
2)POCPT＝0THEN XL＝XR-XR[0]+XL[0]
注释：当光栅尺非线性补偿计算数组位置指针为0时，XL＝XR-XR[0]+XL[0]。
XR[0]表示数组XR第一个数据，即PC端生成非线性误差补偿激光测量数据报表中当前绝对坐标位置值第一个数据；XL[0]表示数组XL第一个数据，即PC端生成非线性误差补偿激光测量数据报表中光栅尺反馈值第一个数据。
3)POCPT＝MAXL+1THEN XL＝XR-XR[MAXL]+XL[MAXL]
注释：当光栅尺非线性补偿计算数组位置指针为最大值时，XL＝XR-XR[MAXL]+XL[MAXL]。XR[MAXL]表示数组XR最后一个数据，即PC端生成非线性误差补偿激光测量数据报表中当前绝对坐标位置值最后一个数据；XL[MAXL]表示数组XL最后一个数据，即PC端生成非线性误差补偿激光测量数据报表中光栅尺反馈值的最后一个数据。
4)
即：
注释：当光栅尺非线性补偿计算数组位置指针为最小值与最大值之间时，代入线性回归方程求得XL(a,b系数已求出，求解简单二元一次方程即可)。
以表1中的数据为例，当前绝对坐标位置值为55008时，代入线性回归方程：
当前光栅尺反馈值为55008，即理论激光测量值为55008时，通过数控系统进行非线性误差逆向补偿后实际光栅尺反馈值为55000。
上述实施例为本发明较佳的实现方案，除此之外，本发明还可以其它方式实现，在不脱离本技术方案构思的前提下任何显而易见的替换均在本发明的保护范围之内。
为了让本领域普通技术人员更方便地理解本发明相对于现有技术的改进之处，本发明的一些附图和描述已经被简化，并且为了清楚起见，本申请文件还省略了一些其它元素，本领域普通技术人员应该意识到这些省略的元素也可构成本发明的内容。

Claims (4)

1.一种数控机床非线性误差补偿方法，其特征在于，是：在PC端生成非线性误差补偿激光测量数据报表，再通过管理软件下载至数控系统，数控系统控制端打开非线性误差补偿报表，自动设定测量轴非线性误差补偿参数，最终完成对光栅尺反馈坐标进行非线性误差补偿与坐标显示输出非线性误差逆向补偿。
2.根据权利要求1所述的数控机床非线性误差补偿方法，其特征在于，具体包括以下步骤：
S101：正确设置数控系统非线性误差补偿参数，选择测量轴，开启非线性误差补偿功能；
S102：打开激光测量系统，激光测量准备就绪，进行非线性误差补偿激光数据处理；
S103：在激光测量系统中运行测量轴对应的通用数控程序；
S104：PC端的激光测量系统生成非线性误差补偿激光测量数据报表；
S105：下载非线性误差补偿数据报表至数控系统；
S106：数控系统控制端打开非线性误差补偿报表，自动设定测量轴非线性误差补偿参数；
S107：数控系统通过自动变步长线性回归的方法自动进行误差补偿，即对光栅尺反馈坐标进行非线性误差补偿与坐标显示输出非线性误差逆向补偿。
3.根据权利要求1所述的数控机床非线性误差补偿方法，其特征在于，所述对光栅尺反馈坐标进行非线性误差补偿：数控系统读取各轴未补偿的机床坐标值，即光栅尺反馈值的数据报表，通过自动变步长线性回归方法计算输出补偿后的机床坐标值，对应激光测量值；
所述坐标显示输出非线性误差逆向补偿：数控系统读取当前各轴绝对坐标位置值，即理论激光测量正确值，通过回归方程计算输出显示坐标的非线性误差逆向补偿值，相当于光栅尺反馈值。
4.根据权利要求3所述的数控机床非线性误差补偿方法，其特征在于，所述自动变步长线性回归方法，包括以下步骤：
1)假设线性回归方程y＝ax+b；开始选择步长数据为n个(xC，yC)；
2)基于最小二乘法拟合直线，得出a,b的参数化方程如下：
其中：
为x0～xn相加所得总和的平均值；
为y0～yn相加所得总和的平均值；
为x0～xn的平方和的平均值；
为xi与yi的乘积总和的平均值(i＝[0,n])；
3)计算线性度量pearson系数；
其中：
Lxx为x的方差；
Lyy为y的方差；
Lxy为x与y的协方差；
4)自动变步长
当计算pearson系数小于系统设定值时，删除测量点第n个点(x，y)；重新计算pearson系数，循环直至pearson系数达到要求；
5)同理采用自动变步长方法完成后面所有数据的线性回归方程。
CN201811179059.6A 2018-10-10 2018-10-10 一种数控机床非线性误差补偿方法 Active CN109240205B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811179059.6A CN109240205B (zh) 2018-10-10 2018-10-10 一种数控机床非线性误差补偿方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811179059.6A CN109240205B (zh) 2018-10-10 2018-10-10 一种数控机床非线性误差补偿方法

Publications (2)

Publication Number Publication Date
CN109240205A true CN109240205A (zh) 2019-01-18
CN109240205B CN109240205B (zh) 2020-02-14

Family

ID=65055645

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811179059.6A Active CN109240205B (zh) 2018-10-10 2018-10-10 一种数控机床非线性误差补偿方法

Country Status (1)

Country Link
CN (1) CN109240205B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3612840A (en) * 1969-09-08 1971-10-12 Kearney & Trecker Corp Linear error compensator for numerically controlled machine tools
JP3727005B2 (ja) * 1999-09-24 2005-12-14 義昭 垣野 パラレルメカニズムを用いた工作機械の誤差補正方法
CN1892175A (zh) * 2005-07-08 2007-01-10 鸿富锦精密工业（深圳）有限公司 三坐标测量机床误差补偿系统及方法
CN104181863A (zh) * 2014-07-28 2014-12-03 上海理工大学 基于位置反馈信号的数控机床误差补偿系统

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3612840A (en) * 1969-09-08 1971-10-12 Kearney & Trecker Corp Linear error compensator for numerically controlled machine tools
JP3727005B2 (ja) * 1999-09-24 2005-12-14 義昭 垣野 パラレルメカニズムを用いた工作機械の誤差補正方法
CN1892175A (zh) * 2005-07-08 2007-01-10 鸿富锦精密工业（深圳）有限公司 三坐标测量机床误差补偿系统及方法
CN104181863A (zh) * 2014-07-28 2014-12-03 上海理工大学 基于位置反馈信号的数控机床误差补偿系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
张虎: "多轴数控机床几何误差的软件补偿技术", 《机械工程学报》 *

Also Published As

Publication number Publication date
CN109240205B (zh) 2020-02-14

Similar Documents

Publication Publication Date Title
JP5400056B2 (ja) 座標測定機を校正するための方法
CN103707132B (zh) 一种基于对刀测量头的感知部位为球形触头的数控机床对刀方法
CN107870604B (zh) 数值控制装置
CN106767505B (zh) 大型复材高精度难点零件的mbd测量工艺
CN108287522B (zh) 一种基于多平台的自动化在机检测方法
JP6316513B1 (ja) 加工プログラム分析装置
CN104570940A (zh) Cnc加工调机系统及方法
JP2008114322A (ja) 工作機械における位置補正方法および位置補正装置
CN112147951B (zh) 机加工设备热误差补偿方法及其装置、系统、介质、终端
CN110531301B (zh) 射频电源输出功率的校准方法和校准装置
CN105823504A (zh) 一种编码器的越零点处理方法
CN107598676A (zh) 大型机床自动校正被测试件分中及加工余量确定方法
ATE319125T1 (de) Verfahren zur automatischen korrektur von systematischen fehlern von mess- und bearbeitungsmaschinen und vorrichtung dafür
CN109070299A (zh) 用于磨轮的基于反馈的修整的系统和方法
CN112101711A (zh) 一种间歇性生产中人工投料防错的流程控制方法
CN109240205A (zh) 一种数控机床非线性误差补偿方法
CN113172476B (zh) 数控机床直线运动重复定位精度快速检测装置及补偿方法
US20020131633A1 (en) System and method for machine vision inspection through platen
CN111552229B (zh) 单机免算点与自动找正加工方法
CN108303669B (zh) 一种提升智能电表校表效率的方法
KR101875391B1 (ko) Cnc 가공기계를 제어하기 위한 프로그램을 기록한 컴퓨터로 읽을 수 있는 매체
CN117709864A (zh) 一种镭雕打码顺序优化方法、系统、电子设备和存储介质
DE102008006927A1 (de) Anordnung zur Korrektur einer Positionssensoreinheit und korrespondierendes Positionierungssystem und Verfahren zur Korrektur eines Positionierungssystem
CN104678887A (zh) 混合料产品制造系统及方法
US20230418256A1 (en) Compensation method for a machining member

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