一种基于泛克里金的机械臂绝对定位误差估计方法

申请号 CN201710942520.8 申请日 2017-10-11 公开(公告)号 CN107443388A 公开(公告)日 2017-12-08
申请人 北京航空航天大学; 发明人 袁培江; 蔡鹦; 陈冬冬; 曹双倩; 高豆豆; 吴旭磊;
摘要 本 发明 公开了一种基于泛克里金的 机械臂 绝对 定位 误差估计方法,属于 机器人 技术领域。本发明方法使用的装置包括机械臂、激光 跟踪 仪和靶球。靶球固定在机械臂末端,标定过程中靶球的中心点设为工具中心点。激光跟踪仪的测量范围应涵盖机械臂所规划的运动空间。本发明方法首先借助激光跟踪仪和靶球测出离线编程中规划出来的一系列点的实际坐标,从而求解理论坐标和实际坐标的 位置 误差。进一步,确定位置误差的漂移函数,通过求解位置误差和漂移函数的差值得到残余函数,进而求解得到残余函数的半变异函数。最后采用泛克里金插值方法求取权重系数,估计出机械臂运动空间某点处的绝对定位误差。本发明具有适应性广、 精度 高、可用于离线编程等特征。
权利要求

1.一种基于泛克里金的机械臂绝对定位误差估计方法,其特征在于:
发明方法所使用的装置包括机械臂、激光跟踪仪和靶球;
进一步地,靶球固定在机械臂末端,标定过程中靶球的中心点设为工具中心点(TCP);
进一步地,激光跟踪仪的测量范围应涵盖机械臂所规划的运动空间;
本发明方法首先借助激光跟踪仪和靶球测出离线编程中规划出来的一系列点的实际坐标,从而求解理论坐标和实际坐标的位置误差;进一步,确定位置误差的漂移函数,通过求解位置误差和漂移函数的差值得到残余函数,进而求解得到残余函数的半变异函数;最后采用泛克里金插值方法求取权重系数,估计出机械臂运动空间某点处的绝对定位误差;
本发明方法具体实施包括如下顺序步骤:
步骤1:统一激光跟踪仪的坐标系和机械臂的基坐标系;
步骤2:通过离线编程在机械臂的运动空间中随机规划出一系列点,作为TCP的理论位置;
步骤3:将机械臂依次驱动到这些点,借助于激光跟踪仪测量TCP在这些点的实际位置;
步骤4:根据这些点的理论位置和实际位置求解这些点的TCP的位置误差;
步骤5:求解位置误差中的确定性部分,即沿X、Y和Z轴的漂移函数;
步骤6:根据所求得的漂移函数,求解位置误差沿X、Y和Z轴的残余函数;
步骤7:根据这些点的理论位置和残余函数建立残余函数沿X、Y和Z轴的半变异函数;
步骤8:进行机械臂运动空间中某点的位置误差估计;
该基于泛克里金的机械臂绝对定位误差估计方法,其效益在于:
1)无需建立机械臂的运动学模型,适用于各种型号的机械臂;
2)考虑了待估计点与邻近观测数据点的空间位置以及各邻近点之间的位置关系,估计结果比传统插值方法更精确;
3)通过确定位置误差的漂移函数,从而对位置误差和漂移函数进行差值得到残余函数,进而求解残余函数的半变异函数,为漂移函数和残余函数的半变异函数同时未知的情况提供了解决办法;
4)基于机械臂运动空间相似性,无需在线测量,适用于离线编程。

说明书全文

一种基于泛克里金的机械臂绝对定位误差估计方法

技术领域

[0001] 本发明涉及机器人技术,尤其涉及一种基于泛克里金的机械臂绝对定位误差估计方法。

背景技术

[0002] 离线编程随着机器人的不断应用与发展已经成为主流,与此同时,机器人的绝对运动精度这一指标的要求在不断提高,特别是在航空领域,要求其在±0.5mm以内。由于机器人的绝对定位精度受到多种因素的影响,如机械加工误差、装配误差、温度变化和末端负载等,传统的基于运动学模型的标定方法只能辨识因几何误差造成的占总误差80%的误差量。
[0003] 基于此,本文所提出的一种基于泛克里金的机械臂绝对定位误差估计方法,通过对已测量点的分析,计算相应的权值,从而预测新的位置的误差。这种方法既不需要花费大量的时间去建立模型、改进模型,并且同时考虑了几何误差和非几何误差,标定精度高。

发明内容

[0004] 本发明为了估计机械臂的绝对定位误差,针对现有技术存在的不足,提出了一种基于泛克里金的机械臂绝对定位误差估计方法。
[0005] 本发明通过以下技术方案实现。
[0006] 本发明方法所使用的装置包括机械臂、激光跟踪仪和靶球。本发明方法首先借助激光跟踪仪和靶球测出离线编程中规划出来的一系列点的实际坐标,从而求解理论坐标和实际坐标的位置误差。进一步,确定位置误差的漂移函数,通过求解位置误差和漂移函数的差值得到残余函数,进而求解得到残余函数的半变异函数。最后采用泛克里金插值方法求取权重系数,估计出机械臂运动空间某点处的绝对定位误差。
[0007] 本发明所述的基于泛克里金的机械臂绝对定位误差估计方法的有益效果在于:
[0008] (1)无需建立机械臂的运动学模型,适用于各种型号的机械臂;
[0009] (2)考虑了待估计点与邻近观测数据点的空间位置以及各邻近点之间的位置关系,估计结果比传统插值方法更精确;
[0010] (3)通过确定位置误差的漂移函数,从而对位置误差和漂移函数进行差值得到残余函数,进而求解残余函数的半变异函数,为漂移函数和残余函数的半变异函数同时未知的情况提供了解决办法;
[0011] (4)基于机械臂运动空间相似性,无需在线测量,适用于离线编程。附图说明
[0012] 图1为本发明中测量机械臂末端位置误差所使用的装置图。
[0013] 图中:1.机械臂、2.激光跟踪仪、3.靶球。

具体实施方式

[0014] 如图1所示,本发明方法所使用的装置包括机械臂1、激光跟踪仪2和靶球3。
[0015] 进一步地,靶球3固定在机械臂1末端,标定过程中靶球3的中心点设为工具中心点(TCP)。
[0016] 进一步地,激光跟踪仪2的测量范围应涵盖机械臂1所规划的运动空间。
[0017] 本发明方法的具体实施步骤如下:
[0018] 步骤1:统一激光跟踪仪2的坐标系和机械臂1的基坐标系;
[0019] 步骤2:通过离线编程在机械臂1的运动空间中随机规划出n个点Pi(xi,yi,zi),i=1,2,…,n,作为TCP的理论位置;
[0020] 步骤3:将机械臂1依次驱动到这n个点,借助于激光跟踪仪2测量TCP在这n个点的实际位置,即Pi’(xi’,yi’,zi’),i=1,2,…,n;
[0021] 步骤4:根据这n个点理论位置和实际位置求解这n个点TCP的位置误差:
[0022] e(Pi)=(ex(Pi),ey(Pi),ez(Pi))=Pi'-Pi,i=1,2,…,n
[0023] 步骤5:求解位置误差中的确定性部分,即沿X、Y和Z轴的漂移函数mx(P)、my(P)和mz(P),以求解mx(P)为例:
[0024]
[0025] 一般地,fx(P)为X轴方向上的k次多项式函数,axl为第k次项fxl(P)的系数;
[0026] 步骤6:根据所求得的漂移函数,求解位置误差沿X、Y和Z轴的残余函数rx(P)、ry(P)和rz(P),以rx(P)为例:
[0027] rx(P)=ex(P)-mx(P)
[0028] 步骤7:根据这n个点的理论位置和残余函数建立残余函数沿X、Y和Z轴的半变异函数γrx(h)、γry(h)和γrz(h),以γrx(h)为例:
[0029]
[0030] 其中,N(h)为距离为h的点对总数。根据γrx(h)的计算结果从球面模型、指数模型、高斯模型和线性模型等模型中选取合适的半变异函数模型;
[0031] 步骤8:进行机械臂1运动空间中某点P0处的位置误差e’(P0)估计:
[0032]
[0033]
[0034]
[0035] 其中,ex’(P0)、ey’(P0)和ez’(P0)分别为P0处的位置误差e’(P0)沿X、Y和Z轴的分量,λxi、λyi和λzi分别为第i个测量点对ex’(P0)、ey’(P0)和ez’(P0)的权重,以λxi,i=1,2,…,n为例,其计算方法如下:
[0036]
[0037] 其中,μxi,i=1,2,…,k为未知,可通过上式算出。
QQ群二维码
意见反馈