CN113770814B - 基于向量方向测量的数控机床平动轴几何误差辨识方法 - Google Patents

基于向量方向测量的数控机床平动轴几何误差辨识方法 Download PDF

Info

Publication number
CN113770814B
CN113770814B CN202110987050.3A CN202110987050A CN113770814B CN 113770814 B CN113770814 B CN 113770814B CN 202110987050 A CN202110987050 A CN 202110987050A CN 113770814 B CN113770814 B CN 113770814B
Authority
CN
China
Prior art keywords
machine tool
numerical control
measurement
laser tracker
translational
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.)
Active
Application number
CN202110987050.3A
Other languages
English (en)
Other versions
CN113770814A (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.)
Southwest Jiaotong University
Original Assignee
Southwest 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 Southwest Jiaotong University filed Critical Southwest Jiaotong University
Priority to CN202110987050.3A priority Critical patent/CN113770814B/zh
Publication of CN113770814A publication Critical patent/CN113770814A/zh
Application granted granted Critical
Publication of CN113770814B publication Critical patent/CN113770814B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B23MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
    • B23QDETAILS, 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/00Arrangements for observing, indicating or measuring on machine tools
    • B23Q17/24Arrangements for observing, indicating or measuring on machine tools using optics or electromagnetic waves
    • B23Q17/2428Arrangements for observing, indicating or measuring on machine tools using optics or electromagnetic waves for measuring existing positions of tools or workpieces
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/02Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]

Landscapes

  • Physics & Mathematics (AREA)
  • Optics & Photonics (AREA)
  • Engineering & Computer Science (AREA)
  • Mechanical Engineering (AREA)
  • Numerical Control (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明公开了基于向量方向测量原理的数控机床平动轴几何误差辨识方法,在机床运动部件上设置不同测量点,并将这些测量点依次相连构成空间向量,利用在机床运动过程中构成空间向量方向的改变仅与机床轴的角位移误差有关,而与线位移误差无关这一特性,利用激光跟踪仪采用多站分时测量原理对机床平动轴运动过程中测量点构成向量方向变化的测量,从而依次辨识出机床平动轴的角位移误差和线位移误差。建立基于空间向量方向测量原理的数控机床平动轴几何误差测量的数学模型,依次推导出了基于精密数控转台测量的激光跟踪仪基站位置自标定算法、测量点确定算法、机床平动轴单项误差分离算法。该方法具有精度高、效率快、成本较低、算法较简单的优点。

Description

基于向量方向测量的数控机床平动轴几何误差辨识方法
技术领域
本发明属于激光精密测量技术领域,具体涉及基于向量方向测量的数控机床平动轴几何误差辨识方法。
背景技术
随着现代制造技术的不断发展,对高精度数控机床的需求量越来越大。精度是数控机床重要的性能指标,直接影响着工件的加工质量。如何经济地提高机床的加工精度是一个关键问题,国内外学者在这方面开展了广泛而深入的研究。相对机床误差硬件补偿技术,误差测量与补偿作为一种提高机床精度的软技术,具有适用性广,成本低等优点,在机床精度补偿领域得到了广泛的应用。机床加工工件的过程中,存在着多项误差因素,如几何误差、热变形误差、切削力误差、装夹误差等,这些误差共同影响着工件最终的加工精度。在这些误差因素中,几何误差是一个主要误差因素,同时几何误差相对稳定,在机床较长工作时间内具有可重复性,易于测量与补偿。通过对机床几何误差的测量与补偿是目前提高机床加工精度的有效途径之一。
目前测量机床平动轴几何误差的方法有很多,主要有以下几种:实物基准法、球杆仪测量法、正交光栅测量法,激光干涉测量法等。近年来,激光跟踪仪作为一种新型三维测量仪器应用到机床和三坐标测量机的误差检测上。激光跟踪仪通常采用球坐标的测量方法,对于这种测量模式而言,角度测量误差是引起坐标测量误差的主要来源。假定激光跟踪仪测角误差为1弧度每秒,由此引起的误差约为5μm每米。目前,对机床的检测精度要求越来越高,现有的这种基于球坐标的测量方法已难以满足机床精度的检测要求。为了提高测量精度,同时考虑测量成本,采用了一种多站分时测量方法。该方法采用GPS测量原理,在测量过程中仅涉及距离测量,有效避免了测角误差对测量结果的影响。但该方法也存在一定的不足:1)测量时,要求机床在3D空间内进给,测量轨迹较为复杂。当机床运动轨迹为2D平面和1D直线时,无法实现对机床误差的辨识;2)在基站标定过程中,基站的标定精度会受到机床自身误差的影响,从而在一定程度上限值了该方法测量精度的进一步提高。如何简化测量过程和误差辨识算法和克服基站标定精度受机床自身误差影响的不足,这对于进一步提高多轴机床整体加工精度有着积极的作用。
发明内容
本发明的目的是解决上述问题,提供一种能够快速、准确辨识出数控机床几何误差的基于向量方向测量的数控机床平动轴几何误差辨识方法。
为解决上述技术问题,本发明的技术方案是:基于向量方向测量的数控机床平动轴几何误差辨识方法,包括以下步骤:
S1、将精密数控转台安装在机床主轴附近,并可跟随主轴或者工作台一起平动;猫眼安装在精密数控转台上,并可跟随精密数控转台一起转动;
S2、将激光跟踪仪分时安装在机床正前方附近,将激光跟踪仪的位置定义为基站位置;
S3、控制精密数控转台每转过一定的角度θ,记下当前激光跟踪仪的测距数据,利用激光跟踪仪测量得到的测量点的距离数据,构建基于测量点欧矢距离的最小二乘模型,确定出激光跟踪仪自身仪器坐标系与精密数控转台坐标系间的齐次变换矩阵,从而确定出激光跟踪仪在转台上的各基站初始位置坐标;
S4、控制机床平动部件按预定的路径在沿机床轴线进给,机床每运动一定距离进行一次测量,安装在机床主轴上的精密转台依次旋转120度,测量点分别命名为A、B、C,激光跟踪仪先后在四个基站位置进行测量;
S5、利用激光跟踪仪测量得到的测量点与各基站的距离数据,建立非线性冗余方程组,按最小二成原理进行求解,确定出测量点的空间坐标;
S6、利用步骤S5中得到的平动部件沿轴线运动时的一系列测量点的空间坐标,通过A、B、C三点两两构成的3个向量
Figure BDA0003231060850000021
向量的方向将会在机床运动过程中发生改变,向量方向的改变仅与机床运动中的角位移误差有关,根据角位移误差引起的齐次变换矩阵和测量点构成一系列向量建立方程组,进行求解,即可辨识出各项角位移误差;
S7、利用步骤S5中得到的平动部件沿轴线运动时的一系列测量点的空间坐标和步骤S6中的到的各项角位移误差,通过取A、B中点为T,B、C中点为R,C、A中点为S,确定一系列新的点,对T、R、S分别建立运动误差方程,构成冗余方程组,带入各项角位移误差求解,即可辨识出各项线位移误差,从而完成完整的数控机床几何误差辨识。
进一步地,所述步骤S2中,至少选定四个基站位置,且任意三者不共线。
进一步地,所述步骤S3还包括以下步骤:
S31、猫眼中心到转台回转中心距离初值的确定;
S32、激光跟踪仪所在基站位置初值的确定;
S33、基站标定。
进一步地,所述步骤S3中当精密数控转台每转过一定的角度θ时,控制转动部件停止运动,并记下当前位置激光跟踪仪的测量结果,然后将激光跟踪仪移动到下一个基站位置,重复上述测量过程,直至在所有基站位置都完成测量。
进一步地,所述步骤S4中,在每个基站位置至少测量四个不同的猫眼位置,以获取数据冗余,提高测量***的可靠性和测量精度。
进一步地,所述步骤S4中,当平动部件运动到各测量点位置时,控制平动部件停止运动;当精密数控转台的转动部件运动到测量位置时,控制转台停止运动。
进一步地,所述步骤S1中将精密数控转台安装在车铣加工中心主轴附近的合适位置,并跟随机床工作台沿数控机床轴线方向一起平动,控制数控机床平动部件沿预先设定的路径沿机床轴线方向进给,在运动区域内,激光跟踪仪实时对平动部件的运动进行跟踪测量。
进一步地,所述数控机床平动部件沿机床轴线方向进给时,每运动一定距离分别对三个测量点进行测量。
本发明的有益效果是:本发明所提供的基于向量方向测量的数控机床平动轴几何误差辨识方法,激光跟踪仪先后在多轴机床转台上至少四个基站位置对多轴机床平动部件相同的点跟踪测量,利用测量得到的测量点距离数据,分别确定出精密数控转台上一系列测量点的空间坐标。通过得到的一系列测量点的空间坐标,先由机床平动件相同位置得到的测量点两两构成空间向量,利用这些向量方向的改变仅与角位移误差有关先辨识出角位移误差;再分别拟合出机床平轴线方向向量,带入角位移误差辨识出线位移误差,从而实现完整的数控机床几何误差辨识。该方法解决了目前利用激光跟踪仪测量数控机床几何误差的难点问题。同时该方法基于GPS原理,在测量过程中,只对距离量进行测量,有效避免了激光跟踪仪测角误差对整体测量精度的影响,从而极大提高了激光跟踪仪现场整体测量精度。利用激光跟踪仪快速测量和精密数控转台高精度的优点,实现数控机床几何误差的高效高精度测量。该方法具有快速、精度高、误差辨识简单等优点,能够数控机床几何误差的快速、准确辨识。为进一步提高数控机床整体加工精度奠定基础。
附图说明
图1是本发明中基站位置、精密数控转台和猫眼相对于数控机床位置的示意图;
图2是本发明数控机床角位移误差辨识方法的原理图;
图3是本发明转台坐标系示意图;
图4是本发明坐标系的变化示意图;
图5是本发明基站标定坐标系示意图。
具体实施方式
下面结合附图和具体实施例对本发明做进一步的说明:
如图1到图5所示,本发明提供的基于向量方向测量的数控机床平动轴几何误差辨识方法,包括以下步骤:
S1、将精密数控转台安装在机床主轴附近,并可跟随主轴或者工作台一起平动;猫眼安装在精密数控转台上,并可跟随精密数控转台一起转动。
在本步骤中将精密数控转台安装在车铣加工中心主轴附近的合适位置,并跟随机床工作台沿数控机床轴线方向一起平动,控制数控机床平动部件沿预先设定的路径沿机床轴线进给,在运动区域内,激光跟踪仪实时对平动部件的运动进行跟踪测量。
S2、将激光跟踪仪分时安装在机床正前方附近,即不同位置处,将激光跟踪仪的位置定义为基站位置。
在本步骤中,至少选定四个基站位置,且任意三者不共线。
S3、控制精密数控转台每转过一定的角度θ,记下当前激光跟踪仪的测距数据,利用激光跟踪仪测量得到的测量点的距离数据,构建基于测量点欧矢距离的最小二乘模型,确定出激光跟踪仪自身仪器坐标系与精密数控转台坐标系间的齐次变换矩阵,从而确定出激光跟踪仪在转台上的各基站初始位置坐标。
在本步骤中当精密数控转台每转过一定的角度θ时,控制转动部件停止运动,并记下当前位置激光跟踪仪的测量结果,然后将激光跟踪仪移动到下一个基站位置,重复上述测量过程,直至在所有基站位置都完成测量。
在基站标定过程中,利用测量得到的大量数据,可建立基站标定的非线性冗余方程组。在求解非线性冗余方程组时,首先应将其线性化,将其转化成线性冗余方程求解。目前线性化的主要途径是将函数在某一点位置处进行泰勒展开,并略去了一阶偏导数以后的各项,这里涉及到选定哪个位置处进行泰勒展开,即确定出泰勒展开位置的初值。初值位置选定的准确与否,直接影响到计算精度和效率。当选定的初值与其真值相差较远时,有可能导致迭代计算不收敛,从而无法得到求解结果,因此如何确定出基站位置的初值是一个关键问题。利用激光跟踪仪测量得到的大量测量点空间坐标,构建基于测量点欧矢距离的激光跟踪仪仪器坐标系与四轴机床自身机床坐标系变换最小二乘模型,从而确定出激光跟踪仪在转台上的各基站初始位置。
步骤S3具体包括以下步骤:
S31、猫眼中心到转台回转中心距离初值的确定。
测量时,利用激光跟踪仪测量得到的转台上一系列测量点在仪器坐标系下的坐标为A′i(x′i,y′i,z′i)如图5所示,一系列测量点分别为K、G、J、I、M……,A′i指代这些点,x′i、y′i、z′i就指这些点的空间坐标,通过空间平面拟合出这些测量点在仪器坐标系下的空间平面方程,假定些测量点所在的空间平面方程为z′=Ax′+By′+C,A、B、C分别为常数。标定的时候转台每转过一定角度测量一次,一共有N个测量点。根据最小二乘拟合原理,目标函数为:
Figure BDA0003231060850000051
根据极值原理,预使F(A,B,C)为最小值,则必有:
Figure BDA0003231060850000052
其中
Figure BDA0003231060850000053
就是根据极值原理求偏导,
Figure BDA0003231060850000054
为偏导数的表示符号,各点坐标是测量出来的已知量,A、B、C是未知量,F是关于A、B、C的函数,
把式1带入式2,由式2可得:
Figure BDA0003231060850000055
通过式(3)可以确定出测量点所在的空间平面方程。
设转台上相邻3点Mi(x′i,y′i,z′i),Mi+1(x′i+1,y′i+1,z′i+1),Mi+2(x′i+2,y′i+2,z′i+2)指三个相邻测量点Mi、Mi+1、Mi+2的空间坐标,如图5中的K、G、J三点确定的空间圆的圆心即转台回转中心在激光跟踪仪仪器测量坐标系下的坐标为O′(u,v,w),易知圆心O′在式(3)确定的平面内,则可以建立如下方程组:
Figure BDA0003231060850000056
通过求解式(4)即可确定出转台回转中心在激光跟踪仪仪器测量坐标系下的位置坐标O′(u,v,w),则转台上一系列测量点所对应的回转半径为:R0是指猫眼距离转台旋转中心的回转半径,猫眼固定在转台上。
Figure BDA0003231060850000057
通过式(5)即可确定出猫眼中心到转台回转中心的初始距离。步骤S31中点的坐标都是在仪器坐标系下的。
S32、激光跟踪仪所在基站位置初值的确定。
激光跟踪仪给出的测量点空间坐标是在仪器坐标系下的坐标,并将该坐标系记为旧坐标系X′Y′Z′。为了确定基站位置的初值,需要将测量点在仪器坐标系下的坐标转化到转台坐标系下,将转台坐标系记为新坐标系O1X1Y1Z1,如图3所示。坐标系的变化如图4所示。
通过坐标系的旋转和平移变换可将仪器坐标系和转台坐标系重合。假定仪器坐标系沿X轴移动Δx,沿Y轴移动Δy,沿Z轴移动Δz,绕X轴转动角度α,绕Y轴转动角度β,绕Z轴转动γ时,能够使仪器坐标系与转台坐标系重合,此时对应的齐次变换矩阵T与式27同理。假定各测量点在仪器坐标系下的坐标为A′i(x′i,y′i,z′i)和步骤S31里面的一样,在转台坐标下各测量点的理论参考坐标为Ai(xi,yi,zi)这是由仪器坐标系切换到转台坐标系下的空间坐标,其中xi=Rcosθi,yi=Rsinθi,zi=0。令U′=[x′i y′i z′i 1]T,U=[xi yi zi 1]T,构建基于测量点欧矢距离的最小二乘模型:
Figure BDA0003231060850000061
在式(6)中,涉及到6个未知参数(Δx,Δy,Δz,α,β,γ),这六个参数是空间左边变换的未知量,分别指沿x、y、z方向的位移量和绕x、y、z轴线的转动量。这里采用非线性最小二乘算法进行求解即可确定出这些未知参数,从而确定出仪器坐标系变换到转台坐标系时所对应的齐次变换矩阵,则激光跟踪仪在转台坐标系下的初始位置也能够很容易确定出来。
S33、基站标定。
测量时,把猫眼安装在精密转台上,猫眼中心距转台回转中心距离为R。
转台每转过一定角度θ,控制转台停止运动,并记下当前激光跟踪仪的测距数据。猫眼中心距转台回转中心距离为R,转台每转θ角度进行一次测量,便可以获得G,H,I等一系列测量点。以精密转台回转中心O1为坐标原点,以O1G所在的直线为Y1轴,向上方向为Y1轴正方向,以过O1点垂直于O1G的直线为X1轴,向后方向为X1轴正方向。根据右手螺旋法则确定出Z1轴,从而在转台上建立基站标定坐标系,如图5所示,其中OXYZ为机床坐标系。
假定激光跟踪仪在转台坐标系下的位置坐标为P1(xp1,yp1,zp1)激光跟踪仪在转台坐标系下的空间坐标,P1下标1指在第一个基站位置,转台依次转过θ角度时,激光跟踪仪的测距读数为l1i,对于测量点可建立如下方程组:l1i下标1指基站位置1,i对应不同的测量点,指在基站位置1时测量第i个点时激光跟踪仪的距离读数
Figure BDA0003231060850000062
式(7)为一非线性冗余方程组,可按最小二乘原理进行求解。
记残差为:fi为残差,i指的是第i个测量点,θi就是第i个测量点队应的转过的角度值
Figure BDA0003231060850000071
取x0,y0,z0,R0为xp1,yp1,zp1,R的近似值,即:
xp1=x0+Δx,yp1=y0+Δy,zp1=z0+Δz,R=R0+ΔR (9)
将式(8)按照Taylor级数在(x0,y0,z0,R0)处展开,为了消除非线性项,略去了一阶偏导数以后的项,可得到
Figure BDA0003231060850000072
Figure BDA0003231060850000073
经过整理将式8、9、11带入10可得:
Figure BDA0003231060850000074
令:
Figure BDA0003231060850000075
则:
Figure BDA0003231060850000076
式(14)是一个线性冗余方程方程组,根据最小二乘原理,目标函数为:
Figure BDA0003231060850000077
由极值原理,预使H为最小,则必有:
Figure BDA0003231060850000078
同时,
Figure BDA0003231060850000079
由此可知,式(16)中各方程求得的极值是最小值,满足最小二乘条件,通过求解式(16),可以得到Δx,Δy,Δz如式(18)所示。
Figure BDA0003231060850000081
当得到Δx,Δy,Δz,ΔR后,可依次根据式(19)(20)确定出激光跟踪仪所在的第一个基站位置以及猫眼中心到转台回转中心的距离。
Figure BDA0003231060850000082
R=R0+ΔR (20)
通过上述过程就可标定出激光跟踪仪在转台坐标系下的第一个基站位置P1,重复上述过程,即可标定出激光跟踪仪所在的其它基站P2、P3、P4的位置坐标。
S4、控制机床平动部件按预定的路径在沿机床轴线进给,机床每运动一定距离进行一次测量,安装在机床主轴上的精密转台依次旋转120度,测量点分别命名为A、B、C,激光跟踪仪先后在四个基站位置进行测量。
在本步骤中,在每个基站位置至少测量四个不同的猫眼位置,以获取数据冗余,提高测量***的可靠性和测量精度。当平动部件运动到各测量点位置时,控制平动部件停止运动;当精密数控转台的转动部件运动到测量位置时,控制转台停止运动。一定距离为50mm。
激光跟踪仪位于第一个基站位置P1处,将目标靶镜猫眼移动到A、B、C,并依次测量A、B、C到基站P1的距离,测量时控制精密数控转台停止转动、机床停止轴线运动,然后控制托板沿x轴移动到下一位置,并重复上述测量过程直至在第一个基站位置P1处完成对机床运动的测量。接下来将激光跟踪仪分别移动到P2、P3、P4位置处,并重复在基站P1位置处的测量过程,直至激光跟踪仪在所有基站位置处完成对机床运动的测量。
S5、利用激光跟踪仪测量得到的测量点与各基站的距离数据,建立非线性冗余方程组,按最小二成原理进行求解,确定出测量点的空间坐标。
假定标定得到的基站位置分别为P1(xp1,yp1,zp1)、P2(xp2,yp2,zp2)、P3(xp3,yp3,zp3)、P4(xp4,yp4,zp4),机床运动过程中测量点Mi到各基站的距离分别为l′1i下标1指基站位置1,i对应不同的测量点,指在基站位置1时测量第i个点时激光跟踪仪的距离读数,l′2i,l′3i,l′4i。根据GPS测量原理,对测量过程中的Mi(xmi,ymi,zmi)可建立如下方程:
Figure BDA0003231060850000091
采用类似方程组(7)的求解方法,即可确定出机床运动过程中测量点的实际坐标Mi(xmi,ymi,zmi)。
S6、利用步骤S5中得到的平动部件沿轴线运动时的一系列测量点的空间坐标,通过A、B、C三点两两构成的3个向量
Figure BDA0003231060850000092
向量的方向将会在机床运动过程中发生改变,向量方向的改变仅与机床运动中的角位移误差有关,根据角位移误差引起的齐次变换矩阵和测量点构成一系列向量建立方程组,进行求解,即可辨识出各项角位移误差。其中的字母A、B、C为点的坐标、角度或直线度误差。
在初始位置处,假定测量点A(xa0,ya0,za0),B(xb0,yb0,zb0),C(xc0,yc0,zc0),此时对应的方向向量依次为:
Figure BDA0003231060850000093
当托板从初始位置沿着x轴运动距离L时,点A,B,C依次运动到A′,B′,C′位置处。
该运动过程中,其理论齐次变换矩阵为:
Figure BDA0003231060850000094
托板运动过程中运动过程中,角度误差εx(x)、εy(x)和εz(x)引起的齐次变换矩阵分别为:
Figure BDA0003231060850000095
定位误差δx(x)、直线度误差δy(x)、δz(x)引起的齐次变换矩阵为:
Figure BDA0003231060850000101
总的误差变换矩阵为:
P=T1T2T3T4T5T6 (26)
经过整理P可化简为:
Figure BDA0003231060850000102
该托板沿x轴运动距离L时,则总的齐次变换矩阵为:
Figure BDA0003231060850000103
假定A′(xa1,ya1,za1),B′(xb1,yb1,zb1),C′(xc1,yc1,zc1),则点A′与点A间应存在如下关系式:
Figure BDA0003231060850000104
通过式(29),我们可以计算出测量点A′的空间坐标。同样的,我们也可以依次计算出测量点B′,C′的空间坐标。
经过整理可得
Figure BDA0003231060850000105
同理可依次计算出B′,C′的空间坐标。
此时
Figure BDA0003231060850000106
Figure BDA0003231060850000111
则在托板运动过程中向量
Figure BDA0003231060850000112
方向的偏差为
Figure BDA0003231060850000113
通过式(32)可以发现,在托板运动过程中向量
Figure BDA0003231060850000114
方向的改变仅与角位移误差εx(x),εy(x),εz(x)有关,而与线位移误差δx(x),δy(x),δz(x)无关。
在图2中,假定在初始位置处
Figure BDA0003231060850000115
Figure BDA0003231060850000116
当机床运动部件沿x轴运动L时,依次得到的新向量
Figure BDA0003231060850000117
Figure BDA0003231060850000118
对于向量
Figure BDA0003231060850000119
与向量
Figure BDA00032310608500001110
应存在如下关系式
Figure BDA00032310608500001111
经过整理可得:
Figure BDA00032310608500001112
同理可建立向量
Figure BDA00032310608500001113
Figure BDA00032310608500001114
以及向量
Figure BDA00032310608500001115
与向量
Figure BDA00032310608500001116
间的相互关系,并整理成矩阵形式可得:
Figure BDA0003231060850000121
在(35)式中的各向量的方向值可由激光跟踪仪测量出的测量点坐标计算得到。通过采用最小二乘方法对该方程组进行求解,即可辨识出各项角运动误差。该方法增加了冗余数据量,从而能够提高各项误差的辨识精度。
S7、利用步骤S5中得到的平动部件沿轴线运动时的一系列测量点的空间坐标和步骤S6中的到的各项角位移误差,通过取A、B中点为T,B、C中点为R,C、A中点为S,确定一系列新的点,对T、R、S分别建立运动误差方程,构成冗余方程组,带入各项角位移误差求解,即可辨识出各项线位移误差,从而完成完整的数控机床几何误差辨识。
在初始位置处,假定测量点A,B中点坐标为T(xt,yt,zt)。当机床运动部件沿着x轴运动距离L,点T运动到T′(xt′,yt′,zt′)位置处,则相应的运动误差为:
Δxt=x′t-xt,Δyt=y′t-yt,Δzt=z′t-zt (36)
这里的Δxt,Δyt和Δzt应由两部分组成,一部分是点T在运动过程中由角运动误差引起的在x轴、y轴、z轴的位移误差tx1,ty1,tz1,另一部分是在运动过程中的线位移误差如δx(x),δy(x),δz(x)。通过前面辨识出运动过程中的角运动误差εx(x),εy(x),εz(x),可计算出当点T运动到T′时,由角运动误差引起的在x轴、y轴、z轴的位移误差tx1,ty1,tz1。接下来可以建立如下方程组:
Figure BDA0003231060850000122
同理可建立B,C中点R和C,A中点S和运动到R′和S′处的运动误差方程,并与等式(37)构成冗余方程组,通过求解便可辨识出各项线位移误差。
在上述误差辨识过程中,激光跟踪仪只需对机床沿单轴运动测量即可辨识出相应的几何误差,机床运动轨迹无需在3D空间,从而简化了测量轨迹。
测量时,将精密数控转台安装在车铣加工中心主轴附近的合适位置,并跟随机床工作台沿数控机床轴线方向一起平动。控制数控机床平动部件沿预先设定的路径沿机床轴线方向进给,运动区域为450mm×550mm×350mm,激光跟踪仪实时对平动部件的运动进行跟踪测量。
数控机床平动部件沿机床轴线方向进给时,每运动一定距离分别对三个测量点进行测量。
平动部件沿x轴方向进给时,每运动50mm设置一个测量点,总测量点数30个。当平动部件运动到各测量点位置时,机床停止运动,安装在机床主轴上的精密转台依次旋转120度,每个位置控制其停下5秒,并记录猫眼在转台上不同位置处的激光跟踪仪测距数据。
如图2所示,当激光跟踪仪在第一个基站位置处完成对各测量点位置测量后,将激光跟踪仪移动到下一个基站位置并重复上述测量过程,直至在所有基站位置完成对测量点的测量。为了减小测量过程中随机误差对测量结果的影响,在每个基站位置对测量点进行三次测量,总测量时间约为3小时左右,测量效率较高。利用上面推导出的测量算法和机床单项误差分离算法,便可辨识出数控机床x轴的各单项几何误差。
本发明实现了机床平动轴几何误差的快速、准确标定,在机床运动部件上设置不同测量点,并将这些测量点依次相连构成空间向量,利用在机床运动过程中构成空间向量方向的改变仅与机床轴的角位移误差有关,而与线位移误差无关这一特性,利用激光跟踪仪采用多站分时测量原理对机床平动轴运动过程中测量点构成向量方向变化的测量,从而依次辨识出机床平动轴的角位移误差和线位移误差。建立基于空间向量方向测量原理的数控机床平动轴几何误差测量的数学模型,依次推导出了基于精密数控转台测量的激光跟踪仪基站位置自标定算法、测量点确定算法、机床平动轴单项误差分离算法。该方法简化了机床平动轴的测量轨迹,降低了误差辨识的复杂性,克服了基站标定精度受机床自身误差影响的不足。具有精度高、效率快、成本较低、算法较简单的优点。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。

Claims (8)

1.基于向量方向测量的数控机床平动轴几何误差辨识方法,其特征在于,包括以下步骤:
S1、将精密数控转台安装在机床主轴附近,并可跟随主轴或者工作台一起平动;猫眼安装在精密数控转台上,并可跟随精密数控转台一起转动;
S2、将激光跟踪仪分时安装在机床正前方附近,将激光跟踪仪的位置定义为基站位置;
S3、控制精密数控转台每转过一定的角度θ,记下当前激光跟踪仪的测距数据,利用激光跟踪仪测量得到的测量点的距离数据,构建基于测量点欧矢距离的最小二乘模型,确定出激光跟踪仪自身仪器坐标系与精密数控转台坐标系间的齐次变换矩阵,从而确定出激光跟踪仪在转台上的各基站初始位置坐标;
S4、控制机床平动部件按预定的路径在沿机床轴线进给,机床每运动一定距离进行一次测量,安装在机床主轴上的精密转台依次旋转120度,测量点分别命名为A、B、C,激光跟踪仪先后在四个基站位置进行测量;
S5、利用激光跟踪仪测量得到的测量点与各基站的距离数据,建立非线性冗余方程组,按最小二乘原理进行求解,确定出测量点的空间坐标;
S6、利用步骤S5中得到的平动部件沿轴线运动时的一系列测量点的空间坐标,通过A、B、C三点两两构成的3个向量
Figure FDA0003231060840000011
向量的方向将会在机床运动过程中发生改变,向量方向的改变仅与机床运动中的角位移误差有关,根据角位移误差引起的齐次变换矩阵和测量点构成一系列向量建立方程组,进行求解,即可辨识出各项角位移误差;
S7、利用步骤S5中得到的平动部件沿轴线运动时的一系列测量点的空间坐标和步骤S6中的到的各项角位移误差,通过取A、B中点为T,B、C中点为R,C、A中点为S,确定一系列新的点,对T、R、S分别建立运动误差方程,构成冗余方程组,带入各项角位移误差求解,即可辨识出各项线位移误差,从而完成完整的数控机床几何误差辨识。
2.根据权利要求1所述的基于向量方向测量的数控机床平动轴几何误差辨识方法,其特征在于,所述步骤S2中,至少选定四个基站位置,且任意三者不共线。
3.根据权利要求1所述的基于向量方向测量的数控机床平动轴几何误差辨识方法,其特征在于,所述步骤S3还包括以下步骤:
S31、猫眼中心到转台回转中心距离初值的确定;
S32、激光跟踪仪所在基站位置初值的确定;
S33、基站标定。
4.根据权利要求1所述的基于向量方向测量的数控机床平动轴几何误差辨识方法,其特征在于,所述步骤S3中当精密数控转台每转过一定的角度θ时,控制转动部件停止运动,并记下当前位置激光跟踪仪的测量结果,然后将激光跟踪仪移动到下一个基站位置,重复上述测量过程,直至在所有基站位置都完成测量。
5.根据权利要求1所述的基于向量方向测量的数控机床平动轴几何误差辨识方法,其特征在于,所述步骤S4中,在每个基站位置至少测量四个不同的猫眼位置,以获取数据冗余,提高测量***的可靠性和测量精度。
6.根据权利要求1所述的基于向量方向测量的数控机床平动轴几何误差辨识方法,其特征在于,所述步骤S4中,当平动部件运动到各测量点位置时,控制平动部件停止运动;当精密数控转台的转动部件运动到测量位置时,控制转台停止运动。
7.根据权利要求1所述的基于向量方向测量的数控机床平动轴几何误差辨识方法,其特征在于,所述步骤S1中将精密数控转台安装在车铣加工中心主轴附近的合适位置,并跟随机床工作台沿数控机床轴线方向一起平动,控制数控机床平动部件沿预先设定的路径沿机床轴线方向进给,在运动区域内,激光跟踪仪实时对平动部件的运动进行跟踪测量。
8.根据权利要求7所述的基于向量方向测量的数控机床平动轴几何误差辨识方法,其特征在于,所述数控机床平动部件沿机床轴线方向进给时,每运动一定距离分别对三个测量点进行测量。
CN202110987050.3A 2021-08-26 2021-08-26 基于向量方向测量的数控机床平动轴几何误差辨识方法 Active CN113770814B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110987050.3A CN113770814B (zh) 2021-08-26 2021-08-26 基于向量方向测量的数控机床平动轴几何误差辨识方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110987050.3A CN113770814B (zh) 2021-08-26 2021-08-26 基于向量方向测量的数控机床平动轴几何误差辨识方法

Publications (2)

Publication Number Publication Date
CN113770814A CN113770814A (zh) 2021-12-10
CN113770814B true CN113770814B (zh) 2022-05-31

Family

ID=78839323

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110987050.3A Active CN113770814B (zh) 2021-08-26 2021-08-26 基于向量方向测量的数控机床平动轴几何误差辨识方法

Country Status (1)

Country Link
CN (1) CN113770814B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115383518B (zh) * 2022-06-17 2023-07-21 中航西安飞机工业集团股份有限公司 数控铣床转台误差快速测量方法
CN116659549B (zh) * 2023-05-29 2024-03-29 浙江省计量科学研究院 一种基于激光跟踪干涉仪的三轴转台垂直度和相交度检测方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102430959A (zh) * 2011-10-10 2012-05-02 西安交通大学 数控机床转台运动误差的快速检测方法
CN103447884B (zh) * 2013-08-02 2016-01-20 西安交通大学 数控机床平动轴几何误差的测量装置及测量与辨识方法
CN103737426B (zh) * 2013-12-24 2015-10-21 西安交通大学 一种数控机床旋转轴几何误差三线测量法
CN106141814B (zh) * 2016-08-04 2017-12-29 天津大学 基于LaserTRACER的数控机床平动轴几何误差检测与辨识方法
CN106363465B (zh) * 2016-11-02 2018-11-13 西南交通大学 多轴数控机床平动轴与旋转轴相互位置关系辨识方法
ES2912525T3 (es) * 2017-12-22 2022-05-26 Ldi Finances Diana tridimensional con doble estructura, dispositivo y procedimiento de medición óptica con dicha diana

Also Published As

Publication number Publication date
CN113770814A (zh) 2021-12-10

Similar Documents

Publication Publication Date Title
CN110524309B (zh) 基于四基站激光追踪***的数控转台几何误差测量方法
CN106363465B (zh) 多轴数控机床平动轴与旋转轴相互位置关系辨识方法
CN113770814B (zh) 基于向量方向测量的数控机床平动轴几何误差辨识方法
Tsutsumi et al. Identification of angular and positional deviations inherent to 5-axis machining centers with a tilting-rotary table by simultaneous four-axis control movements
CN106441117B (zh) 基于多站etalon激光跟踪***的转台误差检测方法
CN102200429B (zh) 基于激光跟踪复合式测量的数控机床精度检测方法
JP2589034B2 (ja) コンピュータ数値制御機械の回転軸に対する加工品の直線変位及び傾斜オフセットを補償する方法、並びにコンピュータ数値制御機械のテーブルの回転軸に対する加工品の直線変位及び傾斜オフセットを補償するシステム
US20160116275A1 (en) Geometric-error identification system and geometric-error identification method
CN111678472B (zh) 四轴坐标测量机回转台误差辨识方法
EP2470968A1 (en) Machine tool calibration method
CN102062575A (zh) 基于激光多路分时测量的数控机床几何精度检测方法
CN110161965B (zh) 一种大型航天机匣斜孔的在机测量方法
CN109884658B (zh) 基于激光追踪仪多站位测量***的激光追踪仪站位方法
CN109884659B (zh) 基于激光追踪仪多站位测量***的大型精密转台标定方法
CN112229321B (zh) 基于lasso算法求解三坐标测量机21项几何误差的方法
Wang et al. The method of geometric error measurement of NC machine tool based on the principle of space vector’s direction measurement
Wang et al. The detection of rotary axis of NC machine tool based on multi-station and time-sharing measurement
CN110686655B (zh) 三维曲面工件的定位方法、装置及***
Zhang et al. Geometric error measuring, modeling, and compensation for CNC machine tools: a review
Lee et al. Optimal on-machine measurement of position-independent geometric errors for rotary axes in five-axis machines with a universal head
Wang et al. A novel approach to separate geometric error of the rotary axis of multi-axis machine tool using laser tracker
Peng et al. Position-dependent geometric errors measurement and identification for rotary axis of multi-axis machine tools based on optimization method using double ball bar
Huang et al. Identification of integrated geometric errors of rotary axis and setup position errors for 5-axis machine tools based on machining test
Guo et al. Continuous measurements with single setup for position-dependent geometric errors of rotary axes on five-axis machine tools by a laser displacement sensor
CN113601269B (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