CN113686278B - 一种高精度工业机器人工具tcp标定方法 - Google Patents

一种高精度工业机器人工具tcp标定方法 Download PDF

Info

Publication number
CN113686278B
CN113686278B CN202110972379.2A CN202110972379A CN113686278B CN 113686278 B CN113686278 B CN 113686278B CN 202110972379 A CN202110972379 A CN 202110972379A CN 113686278 B CN113686278 B CN 113686278B
Authority
CN
China
Prior art keywords
robot
degree
freedom
shafts
coordinate data
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
CN202110972379.2A
Other languages
English (en)
Other versions
CN113686278A (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.)
Nanjing Iungo Technology Co ltd
Original Assignee
Nanjing Iungo Technology Co ltd
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 Nanjing Iungo Technology Co ltd filed Critical Nanjing Iungo Technology Co ltd
Priority to CN202110972379.2A priority Critical patent/CN113686278B/zh
Publication of CN113686278A publication Critical patent/CN113686278A/zh
Application granted granted Critical
Publication of CN113686278B publication Critical patent/CN113686278B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B21/00Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Manipulator (AREA)
  • Numerical Control (AREA)

Abstract

本发明公开了一种高精度工业机器人工具TCP标定方法,包括设有1轴、2轴、3轴、4轴、5轴及6轴的六自由度机器人,利用安装在六自由度机器人末端的靶座、与测量靶座坐标的测量设备,对机器人工具TCP进行标定,六自由度机器人5轴、6轴的杆长小于其他各轴的杆长,所述测量设备为可测量空间坐标的仪器。本发明中TCP标定方法,减小了对机器人绝对定位精度的依赖,在机器人绝对定位精度较低时,也可以标定出具有较高精度的机器人工具TCP。

Description

一种高精度工业机器人工具TCP标定方法
技术领域
本发明涉及工业机器人工具TCP标定技术技术领域,具体为一种高精度工业机器TCP的标定方法。
背景技术
机器人工具TCP是机器人工具坐标系相对于机器人末端法兰坐标系的位置关系。在机器人作业时,需要将机器人工具坐标系下的工作点位通过工具TCP及机器人各轴的姿态计算出工作点位相对于机器人基坐标系或世界坐标系下的位姿。传统的TCP标定方法依赖于机器人各个轴的绝对定位精度,因此,当机器人的绝对定位精度较低时,传统方法标定出的TCP将含有较大的误差,难以应用于有较高精度要求的场合。
发明内容
本发明的目的在于提供一种高精度工业机器人工具TCP标定方法,操作简单,便于达到高精度工具TCP标定。
为了解决上述技术问题,本发明提供如下技术方案:一种高精度工业机器人工具TCP标定方法,包括设有1轴、2轴、3轴、4轴、5轴及6轴的六自由度机器人,利用安装在六自由度机器人末端的靶座、与测量靶座坐标的测量设备,对机器人工具TCP进行标定,包括以下步骤:
步骤1,将六自由度机器人的5轴旋转至零位,其他各轴保持静止,只转动六自由度机器人的6轴,使用测量设备测出六自由度机器人的6轴转动过程中m个位置的坐标数据,并记录为P1,1、P1,2、…、P1,m
步骤2,只转动六自由度机器人的5轴,并保持其他各轴静止,测量设备测出六自由度机器人的5轴转动过程中m个位置的坐标数据,并记录为P2,1、P2,2、…、P2,m
步骤3,通过步骤1中的坐标数据拟合一个空间圆,并过该空间圆的圆心、方向为该空间圆所在平面的法线方向作一条直线;
步骤4,通过步骤2中的坐标数据再次拟合一个空间圆,并过该空间圆的圆心、方向为该空间圆所在平面的法线方向作另外一条直线;
步骤5,计算步骤3和步骤4中所得两条直线的交点,以该交点Pb(xb,yb,zb)为坐标原点、六自由度机器人的5轴方向(ax,ay,az)为z轴、六自由度机器人的6轴方向(ox,oy,oz)为x轴,并按照右手法则确定y轴方向(nx,ny,nz),建立坐标系,得出该坐标系相对于测量设备的位姿矩阵Tb为:
Figure BDA0003226179940000021
步骤6,建立六自由度机器人5轴及六自由度机器人6轴的位姿变换矩阵T5、T6,得出六自由度机器人的正运动学方程:
TbT5T6Ptool=P
式中,P为测量设备所测的坐标值,Ptool为工具TCP;
步骤7,将步骤1中的和步骤2中的测得的坐标数据利用最小二乘法计算工具TCP的值Ptool
Figure BDA0003226179940000022
式中,Pi为第i次测量设备测得的坐标,T5i为第i次测量时六自由度机器人5轴的位姿矩阵,T6i为第i次测量时六自由度机器人6轴的位姿矩阵,n为总共测量的坐标数。
进一步的,步骤2、步骤3中空间圆拟合方法步骤如下:
根据同一组的坐标数据点拟合出平面方程,
设平面方程为
Z=aX+bY+c
式中,X、Y、Z为平面中任一点的坐标,方程中的a、b、c为待定系数,使用最小二乘法拟合,得出目标函数:
Figure BDA0003226179940000031
式中,Xi、Yi、Zi为平面中第i点的坐标;
根据目标函数取得极小值的必要条件,进一步得出:
Figure BDA0003226179940000032
解关于a、b、c的线性方程组即可得到平面方程;
通过将同一组的坐标数据点投影至上述所得平面上,利用投影点拟合该平面上的圆,设圆方程为:
x2+y2+Ax+By+C=0
式中,x、y为圆上任一点的坐标,A、B、C为待定系数,使用最小二乘法拟合,得出目标函数:
Figure BDA0003226179940000033
式中,xi、yi为圆上第i点的坐标;
根据该目标函数取得极小值的必要条件,进一步可得:
Figure BDA0003226179940000034
解上述关于A、B、C的线性方程组即可得到空间圆的方程。
进一步的,步骤1中“零位”指六自由度机器人5轴对应的0刻度角度值。
进一步的,六自由度机器人5轴、6轴的杆长小于其他各轴的杆长。
进一步的,所述测量设备为可测量空间坐标的仪器。
与现有技术相比,本发明所达到的有益效果是:
1、通过测量设备分别测量5轴和6轴转动过程中的多个坐标数据,同时根据两组数据坐标点利用最小二乘法计算工具TCP的值;
2、六自由度机器人的5轴、6轴两轴的杆设置较短,因此5轴、6轴两轴两轴对TCP的误差影响很小,从而达到高精度工具TCP标定;
3、该方法操作过程及后续数据处理过程简单,便于快速完成TCP的标定。
附图说明
附图用来提供对本发明的进一步理解,并且构成说明书的一部分,与本发明的实施例一起用于解释本发明,并不构成对本发明的限制。在附图中:
图1是本发明六自由度机器人、测量设备的结构示意图;
图2是本发明步骤5中的坐标系建立示意图;
图中:1、六自由度机器人;101、1轴;102、2轴;103、3轴;104、4轴;105、5轴;106、6轴;2、靶座;3、测量设备。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
请参阅图1-2,本发明提供技术方案:一种高精度工业机器人工具TCP标定方法,包括设有1轴101、2轴102、3轴103、4轴104、5轴105及6轴106的六自由度机器人1,利用安装在六自由度机器人末端的靶座2、与测量靶座2坐标的测量设备3,对机器人工具TCP进行标定,包括以下步骤:
步骤1,将六自由度机器人1的5轴105旋转至零位,同时其他各轴保持静止,只转动六自由度机器人1的6轴106,测量设备3测出六自由度机器人1的6轴106转动过程中m个位置的坐标数据,并记录为P1,1、P1,2、…、P1,m。此步骤中的“将六自由度机器人1的5轴105旋转至零位”,是指将5轴105旋转至其所对应0刻度角度值。
步骤2,只转动六自由度机器人1的5轴105,并保持其他各轴静止,测量设备3测出六自由度机器人1的5轴105转动过程中m个位置的坐标数据,并记录为P2,1、P2,2、…、P2,m
步骤3,通过步骤1中的坐标数据拟合一个空间圆,并过该空间圆的圆心、方向为该空间圆所在平面的法线方向作一条直线;
步骤4,通过步骤2中的坐标数据再次拟合一个空间圆,并过该空间圆的圆心、方向为该空间圆所在平面的法线方向作另外一条直线;
步骤5,计算步骤3和步骤4所得两条直线的交点,以该交点Pb(xb,yb,zb)为坐标原点、六自由度机器人1的5轴方向(ax,ay,az)为z轴、6轴方向(ox,oy,oz)为x轴,并按照右手法则确定y轴方向(nx,ny,nz),建立坐标系base5(如图2中所示坐标系),可得base5相对于测量设备3的位姿矩阵Tb为:
Figure BDA0003226179940000051
步骤6,建立六自由度机器人1的5轴105及六自由度机器人1的6轴106的位姿变换矩阵T5、T6,得出六自由度机器人1的正运动学方程:
TbT5T6Ptool=P
式中,P为测量设备3所测的坐标值,Ptool为待求得工具TCP;
步骤7,将步骤1中的和步骤2中的测得的坐标数据利用最小二乘法计算工具TCP的值Ptool:
Figure BDA0003226179940000052
式中,Pi为第i次测量设备3测得的坐标,T5i为第i次测量时六自由度机器人1的5轴的位姿矩阵,T6i为第i次测量时六自由度机器人1的6轴的位姿矩阵,n为总共测量的坐标数。
同时,步骤2、步骤3中空间圆拟合方法步骤如下:
根据同一组的坐标数据点拟合出平面方程,
设平面方程为:
Z=aX+bY+c
方程中的X、Y、Z为平面中任一点的坐标,方程中的a、b、c为待定系数,使用最小二乘法拟合,得出目标函数:
Figure BDA0003226179940000061
式中,Xi、Yi、Zi为平面中第i点的坐标;
根据目标函数取得极小值的必要条件,进一步得出:
Figure BDA0003226179940000062
Figure BDA0003226179940000063
解关于a、b、c的线性方程组即可得到平面方程;
通过将同一组的坐标数据点投影至上述所得平面上,利用投影点拟合该平面上的圆,设圆方程为:
x2+y2+Ax+By+C=0
式中,x、y为圆上任一点的坐标,A、B、C为待定系数,使用最小二乘法拟合,得出目标函数:
Figure BDA0003226179940000071
式中,xi、yi为圆上第i点的坐标;
根据该目标函数取得极小值的必要条件,进一步可得:
Figure BDA0003226179940000072
Figure BDA0003226179940000073
解上述关于A、B、C的线性方程组即可得到圆的方程,从而可得到步骤2和步骤3中的拟合空间圆。
所述测量设备3为可测量空间坐标的仪器,用于测量靶座2在各位置的空间坐标数据,架设该测量设备3后,需要将六自由度机器人1的六个轴转动到合适位置,以保证步骤2和步骤3中外部测量设备3能够方便地测量出点位坐标数据。
本发明TCP标定方法不依赖于六自由度机器人1前四个轴的精度,而六自由度机器人5轴105、6轴106的杆长远小于其他各轴的杆长,将5轴105、6轴106两轴的杆长设置为较短,因此5轴105、6轴106对TCP的误差影响很小,从而达到高精度工具TCP标定。此外,该方法操作及后续数据处理都较为简单,TCP精度可达±0.2mm以内,本发明TCP标定方法减小了对机器人绝对定位精度的依赖,在机器人绝对定位精度较低时,也可以标定出具有较高精度的机器人工具TCP。
需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。
最后应说明的是:以上所述仅为本发明的优选实施例而已,并不用于限制本发明,尽管参照前述实施例对本发明进行了详细的说明,对于本领域的技术人员来说,其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (4)

1.一种高精度工业机器人工具TCP标定方法:包括设有1轴、2轴、3轴、4轴、5轴及6轴的六自由度机器人,利用安装在六自由度机器人末端的靶座、与测量靶座坐标的测量设备,对机器人工具TCP进行标定,其特征在于,包括以下步骤:
步骤1,将六自由度机器人的5轴旋转至零位,其他各轴保持静止,只转动六自由度机器人的6轴,使用测量设备测出六自由度机器人的6轴转动过程中m个位置的坐标数据,并记录为P1,1、P1,2、…、P1,m
步骤2,只转动六自由度机器人的5轴,并保持其他各轴静止,测量设备测出六自由度机器人的5轴转动过程中m个位置的坐标数据,并记录为P2,1、P2,2、…、P2,m
步骤3,通过步骤1中的坐标数据拟合一个空间圆,并过该空间圆的圆心、方向为该空间圆所在平面的法线方向作一条直线;
步骤4,通过步骤2中的坐标数据再次拟合一个空间圆,并过该空间圆的圆心、方向为该空间圆所在平面的法线方向作另外一条直线;
步骤5,计算步骤3和步骤4中所得两条直线的交点,以该交点Pb(xb,yb,zb)为坐标原点、六自由度机器人的5轴方向(ax,ay,az)为z轴、六自由度机器人的6轴方向(ox,oy,oz)为x轴,并按照右手法则确定y轴方向(nx,ny,nz),建立坐标系,得出该坐标系相对于测量设备的位姿矩阵Tb为:
Figure FDA0003651904470000011
步骤6,建立六自由度机器人5轴及六自由度机器人6轴的位姿变换矩阵T5、T6,得出六自由度机器人的正运动学方程:
TbT5T6Ptool=P
式中,P为测量设备所测的坐标值,Ptool为工具TCP;
步骤7,将步骤1中的和步骤2中的测得的坐标数据利用最小二乘法计算工具TCP的值Ptool
Figure FDA0003651904470000021
式中,Pi为第i次测量设备测得的坐标,T5i为第i次测量时六自由度机器人5轴的位姿矩阵,T6i为第i次测量时六自由度机器人6轴的位姿矩阵,n为总共测量的坐标数;
六自由度机器人5轴、6轴的杆长小于其他各轴的杆长。
2.根据权利要求1所述的一种高精度工业机器人工具TCP标定方法其特征在于:步骤2、步骤3中空间圆拟合方法步骤如下,
根据将同一组的坐标数据点拟合平面方程,
设平面方程为:
Z=aX+bY+c
式中,X、Y、Z为平面中任一点的坐标,方程中的a、b、c为待定系数,使用最小二乘法拟合,得出目标函数:
Figure FDA0003651904470000022
式中,Xi、Yi、Zi为平面中第i点的坐标;
根据目标函数取得极小值的必要条件,进一步得出:
Figure FDA0003651904470000023
解关于a、b、c的线性方程组即可得到平面方程;
通过将同一组的坐标数据点投影至上述所得平面上,利用投影点拟合该平面上的圆,设圆方程为:
x2+y2+Ax+By+C=0
式中,x、y为圆上任一点的坐标,A、B、C为待定系数,使用最小二乘法拟合,得出目标函数:
Figure FDA0003651904470000031
式中,xi、yi为圆上第i点的坐标;
根据该目标函数取得极小值的必要条件,进一步可得:
Figure FDA0003651904470000032
解上述关于A、B、C的线性方程组即可得到空间圆的方程。
3.根据权利要求1所述的一种高精度工业机器人工具TCP标定方法其特征在于:步骤1中“零位”指六自由度机器人5轴对应的0刻度角度值。
4.根据权利要求1所述的一种高精度工业机器人工具TCP标定方法其特征在于:所述测量设备为可测量空间坐标的仪器。
CN202110972379.2A 2021-08-24 2021-08-24 一种高精度工业机器人工具tcp标定方法 Active CN113686278B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110972379.2A CN113686278B (zh) 2021-08-24 2021-08-24 一种高精度工业机器人工具tcp标定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110972379.2A CN113686278B (zh) 2021-08-24 2021-08-24 一种高精度工业机器人工具tcp标定方法

Publications (2)

Publication Number Publication Date
CN113686278A CN113686278A (zh) 2021-11-23
CN113686278B true CN113686278B (zh) 2022-12-30

Family

ID=78581662

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110972379.2A Active CN113686278B (zh) 2021-08-24 2021-08-24 一种高精度工业机器人工具tcp标定方法

Country Status (1)

Country Link
CN (1) CN113686278B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10859997B1 (en) * 2017-12-04 2020-12-08 Omax Corporation Numerically controlled machining

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10859997B1 (en) * 2017-12-04 2020-12-08 Omax Corporation Numerically controlled machining

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104457645B (zh) * 2014-11-27 2017-08-04 中南大学 一种利用二维测量功能平板的机器人工具中心点标定方法
CN105716525B (zh) * 2016-03-30 2018-03-30 西北工业大学 基于激光跟踪仪的机器人末端执行器坐标系标定方法
CN107121108B (zh) * 2017-06-06 2019-05-10 湖北工业大学 一种快捷自判误机器人工具坐标系标定方法
CN107738254B (zh) * 2017-08-25 2019-12-24 中国科学院光电研究院 一种机械臂坐标系的转换标定方法与***
CN107866823B (zh) * 2017-10-24 2019-10-11 南京工程学院 一种基于位置矢量法的工业机器人几何参数标定方法
CN110640747B (zh) * 2019-11-07 2023-03-24 上海电气集团股份有限公司 机器人的手眼标定方法、***、电子设备和存储介质
CN111203861B (zh) * 2020-01-06 2021-10-19 珠海格力电器股份有限公司 一种机器人工具坐标系的标定方法及标定***
CN112833792B (zh) * 2021-02-24 2022-07-05 南京邮电大学 一种用于六自由度机械臂的精度标定和验证方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10859997B1 (en) * 2017-12-04 2020-12-08 Omax Corporation Numerically controlled machining

Also Published As

Publication number Publication date
CN113686278A (zh) 2021-11-23

Similar Documents

Publication Publication Date Title
US11543234B2 (en) Measurement method for geometric errors of numerical control turntable based on four-station laser tracer system
CN113739717B (zh) 在机测量***中的线激光传感器位姿标定方法
CN106338990B (zh) 基于激光跟踪仪的工业机器人dh参数标定与零位标定方法
CN107042528A (zh) 一种工业机器人的运动学标定***及方法
CN110757504B (zh) 高精度可移动机器人的定位误差补偿方法
CN111862221B (zh) Uvw平台标定方法、设备、纠偏方法、装置及对位***
CN111941425B (zh) 基于激光跟踪仪和双目相机的机器人铣削***工件快速定位方法
CN107560538A (zh) 基于激光跟踪仪的六自由度机器人工具坐标系的标定方法
CN113146613B (zh) 一种工业机器人d-h参数三维自标定校准装置及方法
CN103659806A (zh) 一种工业机器人零位标定方法
CN114012585B (zh) 一种双摆轴式五轴磁流变机床抛光点位置标定方法
CN102151866A (zh) 一种基于三球的加工中心多工位坐标统一方法
CN111678471B (zh) 圆柱坐标测量机转台误差辨识与补偿方法
CN113686278B (zh) 一种高精度工业机器人工具tcp标定方法
CN107471257A (zh) 基于单拉线编码器的机器人几何标定方法
CN112304218A (zh) 工业机器人的工具中心点位置标定方法及***
CN115179323A (zh) 基于远心视觉约束的机器末位姿测量装置及精度提升方法
CN108582047A (zh) 一种六自由度混联抛光机器人位姿精度校准装置及方法
CN113814870B (zh) 一种磁流变抛光工件位姿测算方法及抛光方法
CN111360585A (zh) 一种机器人铣削***中刀具端实时位置误差的获取方法
JPH06249641A (ja) 座標測定機用ロータリテーブルの校正方法
CN112277002B (zh) 基于非完整位姿信息的机器人运动学标定装置及标定方法
CN115200475B (zh) 一种臂载式多视觉传感器快速校正方法
CN113733155B (zh) 六轴工业机器人标定装置和标定方法
CN111006626A (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