CN108723888A - 直线进给轴刀尖点相对于光栅检测点位置偏差估计方法 - Google Patents
直线进给轴刀尖点相对于光栅检测点位置偏差估计方法 Download PDFInfo
- Publication number
- CN108723888A CN108723888A CN201810574797.4A CN201810574797A CN108723888A CN 108723888 A CN108723888 A CN 108723888A CN 201810574797 A CN201810574797 A CN 201810574797A CN 108723888 A CN108723888 A CN 108723888A
- Authority
- CN
- China
- Prior art keywords
- point
- knife
- straight
- data
- line feed
- 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
- 238000000034 method Methods 0.000 title claims abstract description 53
- 238000006073 displacement reaction Methods 0.000 claims abstract description 61
- 230000033001 locomotion Effects 0.000 claims abstract description 14
- 230000008569 process Effects 0.000 claims abstract description 9
- 230000005284 excitation Effects 0.000 claims description 33
- 230000005540 biological transmission Effects 0.000 claims description 14
- 238000012546 transfer Methods 0.000 claims description 14
- 238000012850 discrimination method Methods 0.000 claims description 11
- 238000005070 sampling Methods 0.000 claims description 9
- 238000012360 testing method Methods 0.000 claims description 7
- 230000001133 acceleration Effects 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 230000009471 action Effects 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 108010076504 Protein Sorting Signals Proteins 0.000 claims description 2
- 238000005086 pumping Methods 0.000 claims description 2
- 239000004744 fabric Substances 0.000 claims 1
- 241001269238 Data Species 0.000 description 2
- 238000013480 data collection Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B23—MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
- B23Q—DETAILS, COMPONENTS, OR ACCESSORIES FOR MACHINE TOOLS, e.g. ARRANGEMENTS FOR COPYING OR CONTROLLING; MACHINE TOOLS IN GENERAL CHARACTERISED BY THE CONSTRUCTION OF PARTICULAR DETAILS OR COMPONENTS; COMBINATIONS OR ASSOCIATIONS OF METAL-WORKING MACHINES, NOT DIRECTED TO A PARTICULAR RESULT
- B23Q15/00—Automatic control or regulation of feed movement, cutting velocity or position of tool or work
- B23Q15/007—Automatic control or regulation of feed movement, cutting velocity or position of tool or work while the tool acts upon the workpiece
- B23Q15/013—Control or regulation of feed movement
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B23—MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
- B23Q—DETAILS, COMPONENTS, OR ACCESSORIES FOR MACHINE TOOLS, e.g. ARRANGEMENTS FOR COPYING OR CONTROLLING; MACHINE TOOLS IN GENERAL CHARACTERISED BY THE CONSTRUCTION OF PARTICULAR DETAILS OR COMPONENTS; COMBINATIONS OR ASSOCIATIONS OF METAL-WORKING MACHINES, NOT DIRECTED TO A PARTICULAR RESULT
- B23Q15/00—Automatic control or regulation of feed movement, cutting velocity or position of tool or work
- B23Q15/20—Automatic control or regulation of feed movement, cutting velocity or position of tool or work before or after the tool acts upon the workpiece
- B23Q15/22—Control or regulation of position of tool or workpiece
- B23Q15/225—Control or regulation of position of tool or workpiece in feed control, i.e. approaching of tool or work in successive decreasing velocity steps
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B23—MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
- B23Q—DETAILS, COMPONENTS, OR ACCESSORIES FOR MACHINE TOOLS, e.g. ARRANGEMENTS FOR COPYING OR CONTROLLING; MACHINE TOOLS IN GENERAL CHARACTERISED BY THE CONSTRUCTION OF PARTICULAR DETAILS OR COMPONENTS; COMBINATIONS OR ASSOCIATIONS OF METAL-WORKING MACHINES, NOT DIRECTED TO A PARTICULAR RESULT
- B23Q17/00—Arrangements for observing, indicating or measuring on machine tools
- B23Q17/24—Arrangements for observing, indicating or measuring on machine tools using optics or electromagnetic waves
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Physics & Mathematics (AREA)
- Optics & Photonics (AREA)
- Numerical Control (AREA)
Abstract
本发明公开了一种直线进给轴刀尖点相对于光栅检测点位置偏差估计方法,直线进给轴用以提供高精度的直线进给运动,是多轴数控机床重要组成部分,其控制精度直接影响机床加工精度。当机床直线进给轴上串联其它零部件时,直线进给轴光栅检测点至刀尖点间的柔性结合部如导轨滑块结合部、栓接结合部等,导致进给轴在运动过程中尖点在进给方向位移偏离光栅检测点位移,并造成机床加工精度降低的问题。为降低甚至消除上述偏差,本发明公开了一种直线进给轴刀尖点相对于光栅检测点位置偏差估计方法。可根据直线进给轴运动过程中的光栅数据估计刀尖点相对于光栅检测点位移偏差,为机床高精控制提供技术基础。
Description
技术领域
本发明属于数控机床进给***控制领域,涉及一种直线进给轴刀尖点相对于光栅检测点位置偏差估计方法。
背景技术
数控机床广泛应用于冲压模具、轮船螺旋桨、飞机结构件、汽轮机叶片等复杂型面零件的加工及制造,是制造业的核心装备。提高数控机床的精度是机床设计及制造者不懈努力的目标。
数控机床中,对于安装有其它进给轴、主轴或其他零部件的直线进给轴而言,进给轴光栅位置与刀尖位置在该轴进方向上存在柔性环节,如导轨滑块结合部、栓接结合部等。该直线进给轴在运动过程中,上述柔性环节在电机驱动力作用下产生弹性变形,引起进给轴运动部件产生倾斜、俯仰、偏摆等位移;导致刀尖点在被测实际位置偏离光栅检测点位置,形成位置偏差;并进一步造成机床加工精度降低。
当前,机床进给轴控制器控制的主要目标是降低进给轴跟随误差,即降低插补指令至光栅检测点位移之间的偏差。而控制环外的位置偏差,即光栅检测点至刀尖点的位置偏差,通常认为该偏差很小而很少考虑。这在机床低速运动或执行缓变指令时是成立的。然而,当机床联动高速大曲率轨迹时,进给轴驱动力大幅增加,控制环外位置偏差也大幅增加,导致其对机床加工精度的影响已不能忽略。为降低甚至消除上述偏差,本发明,通过辨识进给轴光栅位移至刀尖点在该轴进给方向位移的传递函数,估计刀尖点相对于光栅检测点位移偏差,可为进给***高精控制提供技术基础。
发明内容
本发明的目的是针对数控机床直线进给轴刀尖点在该轴进给方向位移与该轴光栅测量位移存在偏差的问题,提供一种直线进给轴刀尖点相对于光栅检测点位置偏差估计方法,用以实时估计刀尖点相对于光栅检测点在进给方向的位移偏差,为进给轴高精控制提供技术基础。
为达到上述目的,本发明采用以下技术方案予以实现:
直线进给轴刀尖点相对于光栅检测点位置偏差估计方法,包括以下步骤:
步骤1:激励进给***,采集被测轴光栅及刀尖点在该轴进给方向的位移数据;
步骤2,使用所采集的位移数据辨识获取光栅检测点至刀尖点位移传递函数;
步骤3,估计直线进给轴刀尖点相对于光栅检测点位置偏差。
本发明进一步的改进在于:
步骤1中激励进给***、采集被测轴光栅及刀尖点在该轴进给方向位移数据的具体方法如下:
1-1)生成G代码序列:使用M序列及随机数序列生成激励指令,并将激励指令生成G代码序列,用以激励进给***;
1-2)测前准备:将生成的G代码序列传入机床数控***(3);增大数控***(3)及驱动器(4)中最大加速度、位置环及速度环增益参数值,保证激励频宽;架设激光干涉仪,用以测量并采集刀尖点在被测轴进给方向位移数据;将数据采集卡与被测轴光栅连接,用以采集光栅位移数据;将数据采集卡(9)与激光干涉仪(12)采样频率设置相同,同时配置二者的触发方式均为外部触发,以保证二者数据采集的同步性;
1-3)实施测试:首先,使用触发信号使激光干涉以及数据采集卡同时开始数据采集工作;而后,运行所传入的G代码序列,使进给轴在激励指令的作用下运动;最后,在G代码序列执行完成且机床进给轴停止运动后,使用触发信号使激光干涉仪及数据采集卡同时停止数据采集工作,并存储采集得到的光栅位移数据yopt及刀尖位移数据ytcp。
用以生成G代码序列的激励指令包括位移及速度指令,且均由M序列及随机数序列生成。其中,M序列幅值为1,是由多级线性移位寄存器产生;随机数序列为(0,1)间均匀分布的随机数,是通过混合同余法产生。
由M序列及随机数序列生激励指令的方法如(3)式所示:
式中:xjl为指令坐标位置,vjl为指令进给率,xm(1),xm(2),…,xm(N)为幅值为1的M序列,xr(1),xr(2),…,xr(N)为(0,1)间均匀分布的随机数序列,N为序列长度,Am为幅值,v为匀速信号;激励指令生成过程中,令以防止机床频繁启停;其中,L为进给轴行程,amax为最大加速度。
由激励指令生成G代码序列的方法为:以生成的激励信号序列xjl(1),…,xjl(N)为被辨识轴坐标点,以vjl(1),…,vjl(N)为进给速度,以数控机床快速移动指令G00为进给轴移动G指令,生成为G代码序列,用以激励机床中被辨识的进给轴。
步骤2中光栅检测点位移至刀尖点进给方向位移幅频曲线,其输入输出关系由(4)式所示离散传递函数表示:
式中,ci和dj分别为离散传递函数分子及分母的系数,Nz和Mz分别为分子分母阶数;这些参数是未知的,需要通过辨识来获得。
获取光栅检测点至刀尖点位移传递函数的具体方法如下:
Step 1:设定传函模型分子分母阶数的变化范围Mz∈(Mz_s,Mz_e),Nz∈(Nz_s,Nz_e);其中,Mz_s为传函模型分母阶数下限;Mz_e为传函模型分母阶数上限;Nz_s为传函模型分子阶数下限;Nz_e为传函模型分子阶数上限;且Mz,Mz_s,Mz_e,Nz,Nz_s,Nz_e均为正整数,且满足Mz_s<Mz_e,Nz_s<Nz_e;而后,令i=0;
Step 2:令Mz=Mz_s+i,j=0,开启第一层循环;
Step 3:令Nz=Nz_s+j,开启第二层循环;
Step 4:根据Mz,Nz值生成参数向量、输出向量及测量矩阵;未知的分子及分母系数生成的参数向量如(5)式所示:
若采集的光栅数据为yopt(i),刀尖点在进给方向上的位移数据为ytcp(i),其中i=0,1,2,…,Ndata;将其微分后得到速度数据vopt(i)及vtcp(i),则生成的输出向量为:
生成的测量矩阵如(7)式所示,式中Ts为采样时间;
而后,根据最小二乘法辨识参数向量,最小二乘法辨识参数向量的计算式为:
θ=(WT·W)-1·WT·Y (8)
最后,计算模型误差量;将(4)式变换为差分方程,如下式所示:
式中,tk=k·Ts,k=0,1,2,…,Ndata根据上式,及将辨识得到的参数向量θ,估计***输出,如下式所示:
上式中,为使用辨识模型估计得到的输出;
测量得到的输出量与通过模型估计得到的输出量之差即为模型误差量,其计算式如下所示:
式中:ytcp为测量的刀尖位移数据;为估计的刀尖位移数据;Ndata为刀尖位移数据采样点数;Ts为采样时间;
Step 5:根据模型误差量,存储使模型误差量最小的辨识数据;具体实施方法为:当i=0,j=0时,令error_b=error,best=[Mz,Nz,θ];否则当error<error_b时更新辨识结果error_b=error,best=[Mz,Nz,θ];如此,循环完成后,best值即为使模型误差error最小的模型参数;
Step 6:令j=j+1,并判断Nz是否等于Nz_e,若相等则程序继续向下执行;若不相等则返回Step 3,并继续执行;
Step 7:令i=i+1,并判断Mz是否等于Mz_e,若相等则程序继续向下执行;若不相等则返回Step 2,并继续执行;
Step 8:程序结束,输出辨识结果best值。
Step 4中传递函数系数的辨识方法为最小二乘法;传递函数阶数的辨识方法为搜索法,具体为:设置阶数值定义域,在该定义域内辨识不同阶数的传递函数,选择估计误差最小的阶数值作为阶数参数辨识结果。
步骤3估计直线进给轴刀尖点相对于光栅检测点位置偏差的具体方法如下:
在直线进给轴运动过程中实时采集光栅检测数据,并将其存入寄存器;而后,根据辨识的传递函数模型完成偏差的估计,如(12)式所示:
式中:yopt为光栅测量的位移数据;为根据模型估计的刀尖点相对于光栅检测点的位置偏差。
与现有技术相比,本发明具有以下有益效果:
本发明针对光栅检测点至刀尖点存在柔性环节的数控机床直线进给轴,使用光栅数据估计刀尖点相对于光栅检测点位置偏差。本发明以G代码指令激励进给轴沿进给方向运动,同时采集刀尖点及光栅检测点在进给方向位移数据;使用搜索法及最小二乘法辨识光栅位移至刀尖进给方向位移的传递函数。使用辨识的传递函数估计刀尖点位置相对于光栅检测位置的位移偏差。可根据光栅测量数据实时估计刀尖点相对于光栅检测点位置偏差,为进给***高精控制提供技术基础。
进一步的,本发明所使用的辨识方法,与传统的最小二乘法相比不仅可以辨识里传递函数系数还可以辨识离散传递函数阶数,可以有效防止由于传统辨识方法中,由于阶数选择不合理导致的辨识结果不准确的问题。
附图说明
图1为本发明的指令进给率图;
图2为本发明的指令位移图;
图3为本发明刀尖位移数据及光栅数据测试方案的示意图;
图4为本发明机床进给轴光栅至刀尖点进给方向位移频响图;
图5为本发明传递函数辨识方法流程图。
其中,1-柔性环节;2-激励指令;3-数控***;4-驱动器;5-电机;6-进给轴运动方向;7-光栅检测位移数据;8-光栅检测点;9-数据采集卡;10-光栅;11-刀尖点;12-激光干涉仪;13-光栅检测位移数据;14-同步触发信号;15-计算机;16-进给方向刀尖点位移数据。
具体实施方式
下面结合附图对本发明做进一步详细描述:
参见图1-5,本发明直线进给轴刀尖点相对于光栅检测点位置偏差估计方法,包括以下步骤:
步骤1:激励进给***,采集被测轴光栅及刀尖点在该轴进给方向的位移数据
本步骤的具体实施过程为:
1)生成激励指令:使用M序列及随机数序列生成激励指令,并将激励指令生成G代码序列;
2)测前准备:将生成的G代码序列传入机床数控***3;增大数控***3及驱动器4中最大加速度、位置环及速度环增益参数值,保证激励频宽;架设激光干涉仪,用以测量并采集刀尖点在被测轴进给方向位移数据;将数据采集卡与被测轴光栅连接,用以采集光栅位移数据;将数据采集卡9与激光干涉仪12采样频率设置相同,同时配置二者的触发方式均为外部触发,以保证二者数据采集的同步性;
3)实施测试:首先,使用触发信号使激光干涉以及数据采集卡同时开始数据采集工作;而后,运行所传入的G代码序列,使进给轴在激励指令的作用下运动;最后,在G代码序列执行完成且机床进给轴停止运动后,使用触发信号使激光干涉仪及数据采集卡同时停止数据采集工作,并存储采集得到的光栅及刀尖位移数据。
上述过程中,激励进给***所用的激励指令由位移指令及速度指令组成,且均由M序列及随机数序列生成。所用的M序列幅值为1,是由多级线性移位寄存器产生,如(1)式所示。所用的随机数序列取值范围为(0,1),是通过混合同余法产生,如(2)式所示。
xm(1,…,N)=idinput(N) (1)
xr(1,…,N)=rand(N,1) (2)
其中,N为序列长度。
由M序列及随机数序列生激励指令,包括坐标位置及指令进给率的方法如(3)式所示,生成的指令进给率如图1所示,生成的指令坐标位置图2所示。
式中:xjl为激励指令位移序列,vjl为激励指令速度序列,xm(1),xm(2),…,xm(N)为幅值为1的M序列,xr(1),xr(2),…,xr(N)为(0,1)间均匀分布的随机数序列,Am为幅值,v为匀速信号。激励指令生成过程中,令以防止机床频繁启停;其中,L为进给轴行程,amax为最大加速度。
由激励指令生成G代码序列的方法为:以生成的激励指令位置序列xjl(1),…,xjl(N)为被辨识轴运动坐标点,以激励指令速度序列vjl(1),…,vjl(N)为进给速度,以快速移动指令G00为坐标移动G指令,生成G代码序列,用以激励机床中被辨识的进给轴。
获取光栅检测点位移数据的方法具体如图3所示。将G代码序列输入到机床数控***3中,并运行该G代码程序,使进给轴按指令坐标点及指令进给率运动。同时,以激光干涉仪12测量并采集刀尖点在被测直线进给轴进给方向位移数据ytcp,使用光栅及数据采集卡9分别测量及采集被测直线进给轴光栅检测点进给方向位移数据yopt,用以传递函数的辨识。测试时,为保证激光干涉仪数据与采集卡数据的同步性,将采集卡与激光干涉仪12采样频率设置相同;并使用相同的触发信号使激光干涉仪12及采集卡同时开始及停止数据采集。
步骤2,传递函数辨识
获取光栅检测点至刀尖点位移传递函数的方法为,使用激励指令激励进给轴运动,同时采集光栅检测点及刀尖点位移数据,再使用所采集数据辨识获得。如图4所示,光栅检测点位移至刀尖点进给方向位移幅频曲线,二者间的输入输出关系可由(4)式所示离散传递函数表示。
式中,ci和dj分别为离散传递函数分子及分母的系数,Nz和Mz分别为分子分母阶数;这些参数是未知的,需要通过辨识来获得,辨识流程如图5所示,具体实施步骤为:
Step 1:设定传函模型分子分母阶数的变化范围Mz∈(Mz_s,Mz_e),Nz∈(Nz_s,Nz_e)。其中,Mz_s为传函模型分母阶数下限;Mz_e为传函模型分母阶数上限;Nz_s为传函模型分子阶数下限;Nz_e为传函模型分子阶数上限;且Mz,Mz_s,Mz_e,Nz,Nz_s,Nz_e均为正整数,且满足Mz_s<Mz_e,Nz_s<Nz_e。而后,令i=0。
Step 2:令Mz=Mz_s+i,j=0,开启第一层循环。
Step 3:令Nz=Nz_s+j,开启第二层循环。
Step 4:根据Mz,Nz值生成参数向量、输出向量及测量矩阵。未知的分子及分母系数生成的参数向量如(5)式所示:
若采集的光栅数据为yopt(i)刀尖点在进给方向上的位移数据为ytcp(i),其中i=0,1,2,…,Ndata。将其微分后得到速度数据vopt(i)及vtcp(i),则生成的输出向量为:
生成的测量矩阵如(7)式所示,式中Ts为采样时间。
而后,根据最小二乘法辨识参数向量,最小二乘法辨识参数向量的计算式为:
θ=(WT·W)-1·WT·Y (8)
最后,计算模型误差量。将(4)式变换为差分方程,如下式所示:
式中,tk=k·Ts,k=0,1,2,…,Ndata根据上式,及将辨识得到的参数向量θ,可估计***输出,如下式所示:
上式中,为使用辨识模型估计得到的输出。
传递函数系数的辨识方法为最小二乘法;传递函数阶数的辨识方法为搜索法,具体为:设置阶数值定义域,在该定义域内辨识不同阶数的传递函数,选择估计误差最小的阶数值作为阶数参数辨识结果。测量得到的输出量与通过模型估计得到的输出量之差即为模型误差量,其计算式如下所示:
式中:ytcp为测量的刀尖位移数据;为估计的刀尖位移数据;Ndata为刀尖位移数据采样点数;Ts为采样时间。
Step 5:根据模型误差量,存储使模型误差量最小的辨识数据。具体实施方法为:当i=0,j=0时,令error_b=error,best=[Mz,Nz,θ];否则当error<error_b时更新辨识结果error_b=error,best=[Mz,Nz,θ];如此,循环完成后,best值即为使模型误差error最小的模型参数。
Step 6:令j=j+1,并判断Nz是否等于Nz_e,若相等则程序继续向下执行;若不相等则返回Step 3,并继续执行。
Step 7:令i=i+1,并判断Mz是否等于Mz_e,若相等则程序继续向下执行;若不相等则返回Step 2,并继续执行。
Step 8:程序结束,输出辨识结果best值。
所使用的辨识方法,与传统的最小二乘法相比不仅可以辨识里传递函数系数还可以辨识离散传递函数阶数,可以有效防止由于传统辨识方法中,由于阶数选择不合理导致的辨识结果不准确的问题。
步骤3,估计直线进给轴刀尖点相对于光栅检测点位置偏差
首先,在直线进给轴运动过程中实时采集光栅检测数据,并将其存入寄存器。而后,根据辨识的传递函数模型即可完成偏差的估计,如(12)式所示
式中:yopt为光栅测量的位移数据;为根据模型估计的刀尖点相对于光栅检测点的位置偏差。
以上内容仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明权利要求书的保护范围之内。
Claims (9)
1.直线进给轴刀尖点相对于光栅检测点位置偏差估计方法,其特征在于,包括以下步骤:
步骤1:激励进给***,采集被测轴光栅及刀尖点在该轴进给方向的位移数据;
步骤2,使用所采集的位移数据辨识获取光栅检测点至刀尖点位移传递函数;
步骤3,估计直线进给轴刀尖点相对于光栅检测点位置偏差。
2.根据权利要求1所述的直线进给轴刀尖点相对于光栅检测点位置偏差估计方法,其特征在于,步骤1中激励进给***、采集被测轴光栅及刀尖点在该轴进给方向位移数据的具体方法如下:
1-1)生成G代码序列:使用M序列及随机数序列生成激励指令,并将激励指令生成G代码序列,用以激励进给***;
1-2)测前准备:将生成的G代码序列传入机床数控***(3);增大数控***(3)及驱动器(4)中最大加速度、位置环及速度环增益参数值,保证激励频宽;架设激光干涉仪,用以测量并采集刀尖点在被测轴进给方向位移数据;将数据采集卡与被测轴光栅连接,用以采集光栅位移数据;将数据采集卡(9)与激光干涉仪(12)采样频率设置相同,同时配置二者的触发方式均为外部触发,以保证二者数据采集的同步性;
1-3)实施测试:首先,使用触发信号使激光干涉以及数据采集卡同时开始数据采集工作;而后,运行所传入的G代码序列,使进给轴在激励指令的作用下运动;最后,在G代码序列执行完成且机床进给轴停止运动后,使用触发信号使激光干涉仪及数据采集卡同时停止数据采集工作,并存储采集得到的光栅位移数据yopt及刀尖位移数据ytcp。
3.根据权利要求2所述的直线进给轴刀尖点相对于光栅检测点位置偏差估计方法,其特征在于,用以生成G代码序列的激励指令包括位移及速度指令,且均由M序列及随机数序列生成;其中,M序列幅值为1,是由多级线性移位寄存器产生;随机数序列为(0,1)间均匀分布的随机数,是通过混合同余法产生。
4.根据权利要求2所述的直线进给轴刀尖点相对于光栅检测点位置偏差估计方法,其特征在于,由M序列及随机数序列生激励指令的方法如(3)式所示:
式中:xjl为指令坐标位置,vjl为指令进给率,xm(1),xm(2),…,xm(N)为幅值为1的M序列,xr(1),xr(2),…,xr(N)为(0,1)间均匀分布的随机数序列,N为序列长度,Am为幅值,v为匀速信号;激励指令生成过程中,令以防止机床频繁启停;其中,L为进给轴行程,amax为最大加速度。
5.根据权利要求4所述的直线进给轴刀尖点相对于光栅检测点位置偏差估计方法,其特征在于,由激励指令生成G代码序列的方法为:以生成的激励信号序列xjl(1),…,xjl(N)为被辨识轴坐标点,以vjl(1),…,vjl(N)为进给速度,以数控机床快速移动指令G00为进给轴移动G指令,生成为G代码序列,用以激励机床中被辨识的进给轴。
6.根据权利要求1所述的直线进给轴刀尖点相对于光栅检测点位置偏差估计方法,其特征在于,步骤2中光栅检测点位移至刀尖点进给方向位移的输入输出关系由(4)式所示离散传递函数表示:
式中,ci和dj分别为离散传递函数分子及分母的系数,Nz和Mz分别为分子分母阶数;这些参数是未知的,需要通过辨识来获得。
7.根据权利要求6所述的直线进给轴刀尖点相对于光栅检测点位置偏差估计方法,其特征在于,获取光栅检测点至刀尖点位移传递函数的具体方法如下:
Step1:设定传函模型分子分母阶数的变化范围Mz∈(Mz_s,Mz_e),Nz∈(Nz_s,Nz_e);其中,Mz_s为传函模型分母阶数下限;Mz_e为传函模型分母阶数上限;Nz_s为传函模型分子阶数下限;Nz_e为传函模型分子阶数上限;且Mz,Mz_s,Mz_e,Nz,Nz_s,Nz_e均为正整数,且满足Mz_s<Mz_e,Nz_s<Nz_e;而后,令i=0;
Step2:令Mz=Mz_s+i,j=0,开启第一层循环;
Step3:令Nz=Nz_s+j,开启第二层循环;
Step4:根据Mz,Nz值生成参数向量、输出向量及测量矩阵;未知的分子及分母系数生成的参数向量如(5)式所示:
若采集的光栅数据为yopt(i),刀尖点在进给方向上的位移数据为ytcp(i),其中i=0,1,2,…,Ndata;将其微分后得到速度数据vopt(i)及vtcp(i),则生成的输出向量为:
生成的测量矩阵如(7)式所示,式中Ts为采样时间;
而后,根据最小二乘法辨识参数向量,最小二乘法辨识参数向量的计算式为:
θ=(WT·W)-1·WT·Y (8)
最后,计算模型误差量;将(4)式变换为差分方程,如下式所示:
式中,tk=k·Ts,k=0,1,2,…,Ndata根据上式,及将辨识得到的参数向量θ,估计***输出,如下式所示:
上式中,为使用辨识模型估计得到的输出;
测量得到的输出量与通过模型估计得到的输出量之差即为模型误差量,其计算式如下所示:
式中:ytcp为测量的刀尖位移数据;为估计的刀尖位移数据;Ndata为刀尖位移数据采样点数;Ts为采样时间;
Step5:根据模型误差量,存储使模型误差量最小的辨识数据;具体实施方法为:当i=0,j=0时,令error_b=error,best=[Mz,Nz,θ];否则当error<error_b时更新辨识结果error_b=error,best=[Mz,Nz,θ];如此,循环完成后,best值即为使模型误差error最小的模型参数;
Step6:令j=j+1,并判断Nz是否等于Nz_e,若相等则程序继续向下执行;若不相等则返回Step3,并继续执行;
Step7:令i=i+1,并判断Mz是否等于Mz_e,若相等则程序继续向下执行;若不相等则返回Step2,并继续执行;
Step8:程序结束,输出辨识结果best值。
8.根据权利要求7所述的直线进给轴刀尖点相对于光栅检测点位置偏差估计方法,其特征在于,Step4中传递函数系数的辨识方法为最小二乘法;传递函数阶数的辨识方法为搜索法,具体为:设置阶数值定义域,在该定义域内辨识不同阶数的传递函数,选择估计误差最小的阶数值作为阶数参数辨识结果。
9.根据权利要求1所述的直线进给轴刀尖点相对于光栅检测点位置偏差估计方法,其特征在于,步骤3估计直线进给轴刀尖点相对于光栅检测点位置偏差的具体方法如下:
在直线进给轴运动过程中实时采集光栅检测数据,并将其存入寄存器;而后,根据辨识的传递函数模型完成偏差的估计,如(12)式所示:
式中:yopt为光栅测量的位移数据;为根据模型估计的刀尖点相对于光栅检测点的位置偏差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810574797.4A CN108723888B (zh) | 2018-06-05 | 2018-06-05 | 直线进给轴刀尖点相对于光栅检测点位置偏差估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810574797.4A CN108723888B (zh) | 2018-06-05 | 2018-06-05 | 直线进给轴刀尖点相对于光栅检测点位置偏差估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108723888A true CN108723888A (zh) | 2018-11-02 |
CN108723888B CN108723888B (zh) | 2020-07-28 |
Family
ID=63932259
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810574797.4A Expired - Fee Related CN108723888B (zh) | 2018-06-05 | 2018-06-05 | 直线进给轴刀尖点相对于光栅检测点位置偏差估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108723888B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112327757A (zh) * | 2020-11-19 | 2021-02-05 | 清华大学 | 使用在线测量的数控机床加减速智能优化方法及装置 |
CN113310673A (zh) * | 2021-04-02 | 2021-08-27 | 深圳市世宗自动化设备有限公司 | 重复精度检测方法、装置、计算机设备及其存储介质 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH04331041A (ja) * | 1991-04-25 | 1992-11-18 | Okuma Mach Works Ltd | 工作機械における刃先位置補正装置 |
CN101968641A (zh) * | 2010-07-08 | 2011-02-09 | 西华大学 | 一种机床xy平面误差修正*** |
CN104317246A (zh) * | 2014-08-27 | 2015-01-28 | 华中科技大学无锡研究院 | 一种对弱刚性刀具多轴加工路径进行让刀补偿的方法 |
CN105334804A (zh) * | 2015-12-01 | 2016-02-17 | 华中科技大学 | 一种基于数控编程的数控机床结构参数的辨识方法 |
CN105881103A (zh) * | 2014-05-09 | 2016-08-24 | 科德数控股份有限公司 | 利用激光直接反馈的直线误差补偿方法 |
CN205817105U (zh) * | 2016-05-31 | 2016-12-21 | 于浩源 | 一种用于五轴激光切割机床的误差补偿装置 |
-
2018
- 2018-06-05 CN CN201810574797.4A patent/CN108723888B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH04331041A (ja) * | 1991-04-25 | 1992-11-18 | Okuma Mach Works Ltd | 工作機械における刃先位置補正装置 |
CN101968641A (zh) * | 2010-07-08 | 2011-02-09 | 西华大学 | 一种机床xy平面误差修正*** |
CN105881103A (zh) * | 2014-05-09 | 2016-08-24 | 科德数控股份有限公司 | 利用激光直接反馈的直线误差补偿方法 |
CN104317246A (zh) * | 2014-08-27 | 2015-01-28 | 华中科技大学无锡研究院 | 一种对弱刚性刀具多轴加工路径进行让刀补偿的方法 |
CN105334804A (zh) * | 2015-12-01 | 2016-02-17 | 华中科技大学 | 一种基于数控编程的数控机床结构参数的辨识方法 |
CN205817105U (zh) * | 2016-05-31 | 2016-12-21 | 于浩源 | 一种用于五轴激光切割机床的误差补偿装置 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112327757A (zh) * | 2020-11-19 | 2021-02-05 | 清华大学 | 使用在线测量的数控机床加减速智能优化方法及装置 |
CN113310673A (zh) * | 2021-04-02 | 2021-08-27 | 深圳市世宗自动化设备有限公司 | 重复精度检测方法、装置、计算机设备及其存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN108723888B (zh) | 2020-07-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107553220B (zh) | 一种数控机床综合误差实时补偿方法 | |
CN101797704B (zh) | 数控滚齿机热变形误差补偿方法 | |
EP2664891B1 (en) | Shape measuring apparatus | |
TWI540402B (zh) | 馬達控制裝置 | |
CN108459559B (zh) | 五轴联动加工机床轨迹误差实时监测方法 | |
CN102944197B (zh) | 一种双转台结构的五轴加工中心精度检测方法 | |
CN109500619B (zh) | 机床的数控装置和数控方法 | |
CN107044837B (zh) | 用于标定检测工具坐标系的方法、装置以及控制设备 | |
CN102785128A (zh) | 面向数控车床的零件加工精度在线检测***及检测方法 | |
CN109782815B (zh) | 基于多轴联动***的复杂型面自适应测量路径规划方法 | |
CN110837248B (zh) | 机床健康度测评方法、***、综合***及数控机床 | |
CN109682626B (zh) | 机器人速度测试方法和装置 | |
CN109968402B (zh) | 回程间隙测量方法及控制方法、装置、***和存储介质 | |
CN108723888A (zh) | 直线进给轴刀尖点相对于光栅检测点位置偏差估计方法 | |
CN104808694B (zh) | 马达控制装置及该装置的校正数据生成方法 | |
CN109116808A (zh) | 一种基于数控机床大数据的联动轨迹轮廓精度预测方法 | |
CN106181527B (zh) | 数控机床用滚珠丝杠运动*** | |
CN107957234A (zh) | 一种测量自由曲面任意点处法向矢量的方法及装置 | |
Miller et al. | Improved machine tool linear axis calibration through continuous motion data capture | |
CN101187538A (zh) | 螺纹量规的检测方法 | |
CN111504223B (zh) | 一种基于线激光传感器的叶片轮廓测量方法及装置和*** | |
CN112008490B (zh) | 一种全闭环伺服***的数控机床反向间隙检测方法 | |
CN102658503A (zh) | 一种基于内置传感器的数控机床进给***的模态测试方法 | |
CN110977612B (zh) | Cnc数控加工在线测量误差修正方法及*** | |
CN108896049A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200728 |
|
CF01 | Termination of patent right due to non-payment of annual fee |