CN100398274C - 基于平板测量的机器人工具中心点三分量校准法 - Google Patents

基于平板测量的机器人工具中心点三分量校准法 Download PDF

Info

Publication number
CN100398274C
CN100398274C CNB200510087041XA CN200510087041A CN100398274C CN 100398274 C CN100398274 C CN 100398274C CN B200510087041X A CNB200510087041X A CN B200510087041XA CN 200510087041 A CN200510087041 A CN 200510087041A CN 100398274 C CN100398274 C CN 100398274C
Authority
CN
China
Prior art keywords
robot
coordinate
plane
point
normal vector
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.)
Expired - Fee Related
Application number
CNB200510087041XA
Other languages
English (en)
Other versions
CN1903525A (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.)
Tianjin Zhitong Information System Integration Co ltd
Original Assignee
ZHITONG ROBOT SYSTEM CO Ltd LANGFANG CITY
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 ZHITONG ROBOT SYSTEM CO Ltd LANGFANG CITY filed Critical ZHITONG ROBOT SYSTEM CO Ltd LANGFANG CITY
Priority to CNB200510087041XA priority Critical patent/CN100398274C/zh
Publication of CN1903525A publication Critical patent/CN1903525A/zh
Application granted granted Critical
Publication of CN100398274C publication Critical patent/CN100398274C/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Manipulator (AREA)
  • Length Measuring Devices With Unspecified Measuring Means (AREA)
  • Numerical Control (AREA)

Abstract

本发明涉及一种基于平板测量的机器人工具中心点三分量校准法,机器人末端固定一个测量用探头,一个表面平整的平板位于机器人的活动范围之内,包括以下步骤:步骤1:机器人末端保持一定的姿态平动,用探头接触平板上不在一条直线上的三个或三个以上的点,得到这些点的坐标和位置,利用这些坐标和位置计算表面所在平面的法向量;步骤2:机器人末端坐标系的X、Y和Z轴方向分别与表面所在平面的法向量方向相同时,分别用探头接触平板的一点,测量点的坐标;然后与法向量方向相反,分别用探头接触平板的一点,测量点的坐标;步骤3:利用步骤2测量出的点的坐标计算出机器人工具中心点在机器人末端坐标系下的位置。

Description

基于平板测量的机器人工具中心点三分量校准法
技术领域
本发明涉及机器人工具中心点三分量校准法,尤其涉及通过测量平板上的点的坐标的,基于几何平面的机器人工具中心点三分量单独校准的机器人工具中心点三分量校准法。
背景技术
机器人TCP(Tool Center Point)校准指的是校准工具中心点在机器人末端坐标系(对工具被机器人夹持而言)或机器人基坐标系(对工具固定的情况而言)下的位置,通常情况下,机器人TCP的各个参数可以通过两种方法得到。一种就是把TCP所有的变量都单独测量,然后利用约束条件将TCP的所有变量同时计算出来。另一种方法就是把TCP所有的变量都单独分开来测量,然后单独计算,每次测量只计算其中一个变量。
上述TCP校准法常用的是“四点法”,在该方法的校准过程中,操作者控制机器人以多个姿态(通常为四种)使机器人工具中心点达到空间的某一固定点,从而利用约束关系得到TCP校准结果。这种方法,所有待测变量在测量中都是相关的,每一个变量校准的精确程度处决于这个变量在测量中的相关程度。在操作上,该方法要求操作者以差异较大的几种姿态来控制机器人达到固定空间点,以减小求解的变量的相关性,但这种方法在操作上很难量化。并且由于校准过程中机器人各关节变化大,标定过程中引入了较多的机器人系统误差,这会影响标定精度。另外,该方法通常是由人眼判断工具中心点是否与空间固定点重合,这也使得标定过程中不可避免地引入了人眼判断误差,降低了标定精度。
发明内容
本发明的目的是把机器人TCP所有的变量都单独分开来测量和计算,从而容易控制测量过程中校准的精确程度。
为实现上述目的,本发明提供了一种基于平板测量的机器人工具中心点在机器人末端坐标系下的位置校准法,机器人末端固定一个测量用探头,一个表面平整的平板位于机器人的活动范围之内,本方法包括以下步骤:
步骤1:所述机器人末端保持一定的姿态平动,用所述探头接触所述平板上不在一条直线上的三个或三个以上的点,得到这些点的坐标和位置,利用这些坐标和位置计算所述表面所在平面的法向量;
步骤2:机器人末端坐标系的X、Y和Z轴方向分别与所述表面所在平面的法向量方向相同时,分别用所述探头接触所述平板的一点,测量所述点的坐标;然后与法向量方向相反,分别用所述探头接触所述平板的一点,测量所述点的坐标;
步骤3:利用步骤1得出的平面法向量和步骤2测量出的所有点的坐标计算出机器人工具中心点在机器人末端坐标系下的位置。
所述步骤2具体为,首先机器人末端坐标系的X轴方向与所述表面所在平面的法向量方向相同时,用所述探头接触所述平板的一点,测量所述点的坐标,然后机器人末端坐标系的Y轴方向与所述表面所在平面的法向量方向相同时,用所述探头接触所述平板的一点,测量所述点的坐标,接着机器人末端坐标系的Z轴方向与所述表面所在平面的法向量方向相同时,用所述探头接触所述平板的一点,测量所述点的坐标;再后机器人末端坐标系的X轴方向与所述表面所在平面的法向量方向相反时,用所述探头接触所述平板的一点,测量所述点的坐标,然后机器人末端坐标系的Y轴方向与所述表面所在平面的法向量方向相反时,用所述探头接触所述平板的一点,测量所述点的坐标,最后机器人末端坐标系的Z轴方向与所述表面所在平面的法向量方向相反时,用所述探头接触所述平板的一点,测量所述点的坐标。
所述步骤中的机器人为六自由度的机器人。所述步骤2中,机器人改变姿态,用所述探头接触所述平板,测量所述点的坐标。所述步骤中通过读取机器人控制器而得到所述点的坐标和位置。
本发明还提供了一种通过测平板上点的坐标的机器人的固定工具中心点在机器人基座坐标系下的位置校准法,机器人末端固定一个平板,该平板的两个表面均为光滑且平行,为第一表面和第二表面,一个测量用的探头固定于机器人的活动范围之内,本方法包括以下步骤:
步骤1:机器人末端保持一定的姿态平动,使得探头接触平板的第一表面上不在一条直线上的三个或三个以上的点,得到这些点的坐标和位置,利用这些坐标和位置计算该第一表面的法向量;
步骤2:机器人基座坐标系的X、Y和Z轴方向分别与所述第一表面所在平面的法向量方向相反时,分别用所述探头接触所述第一表面的一点,测量所述点的坐标;然后与法向量方向相同,分别用所述探头接触所述第二表面的一点,分别测量所述点的坐标;
步骤3:利用平板厚度值和步骤2测量出的所有点的坐标计算出机器人工具中心点在机器人基座坐标系下的位置。
所述步骤2具体为,首先所述第一表面所在平面的法向量方向与机器人基座坐标系的X轴方向相同时,用所述探头接触所述第一表面上的一点,测量所述点的坐标,接着所述第一表面所在平面的法向量方向与机器人基座坐标系的X轴方向相反时,用所述探头接触所述第二表面上的一点,测量所述点的坐标;然后所述第一表面所在平面的法向量方向与机器人基座坐标系的Y轴方向相同时,用所述探头接触所述第一表面的一点,测量所述点的坐标,接着将所述第一表面所在平面的法向量方向与机器人基座坐标系的Y轴方向相反时,用所述探头接触所述第二表面上的一点,测量所述点的坐标;最后所述第一表面所在平面的法向量方向与机器人基座坐标系的Z轴方向相同时,用所述探头接触所述第一表面上的一点,测量所述点的坐标,接着所述第一表面所在平面的法向量方向与机器人基座坐标系的Z轴方向相反时,用所述探头接触所述第二表面上的一点,测量所述点的坐标。
所述步骤中的机器人为六自由度的机器人。所述步骤2中,机器人改变姿态,用所述探头接触所述平板,测量所述点的坐标。所述步骤中通过读取机器人控制器而得到所述点的坐标和位置。
因此,本发明将TCP所有的变量都单独分开来测量和计算,每次测量只测量和计算其中一个变量。由于所有的变量都被单独分离开,所以很容易控制测量过程中校准的精确程度。
附图说明
图1为本发明基于平板测量的机器人工具中心点在机器人末端坐标系下的位置校准法实施例1的结构示意图。
图2为本发明基于平板测量的机器人工具中心点在机器人末端坐标系下的位置校准法实施例1的方法流程图。
图3为本发明基于平板测量的机器人工具中心点在机器人基座坐标系下的位置校准法实施例2的结构示意图。
图4为本发明基于平板测量的机器人工具中心点在机器人基座坐标系下的位置校准法实施例2的方法流程图。
具体实施方式
本发明的思想是把机器人TCP的所有的变量都单独分开来计算,每次测量只计算其中一个变量,基于所有的变量都被单独分离开,从而容易控制测量过程中校准的精确程度。
下面用实施例,对本发明的技术方案做进一步的详细描述。
实施例1:如图1所示,为本发明实施例1的结构示意图,机器人1为六自由度工业机器人,底座10为机器人基座,tool0为机器人末端11的坐标系,测量用的探头2固定在机器人末端上。还包括一个厚度已知的表面平整的平板3。如图2所示,为本发明机器人TCP在机器人末端坐标系下的位置校准法的流程图,包括如下步骤:
步骤101:将一个表面平整的平板置于机器人的活动范围内;
步骤102:机器人末端保持一定的姿态平动,测量所述表面上不在一条直线上的三个或三个以上的点的坐标,利用这些坐标计算所述表面所在平面的法向量;
可以证明:在对平面法向量的测量中,一旦机器人末端保持一定的姿态测量,平面法向量的精确值与TCP的初始值是没有关系的。因此既然平面法向量的测量与未知的TCP值没有关系,那么法向量方向的测量它就是准确的。
步骤103:首先机器人末端坐标系的X轴方向与所述表面所在平面的法向量方向相同时,用所述探头接触所述平板的一点,测量所述点的坐标,然后机器人末端坐标系的Y轴方向与所述表面所在平面的法向量方向相同时,用所述探头接触所述平板的一点,测量所述点的坐标,接着机器人末端坐标系的Z轴方向与所述表面所在平面的法向量方向相同时,用所述探头接触所述平板的一点,测量所述点的坐标;再后机器人末端坐标系的X轴方向与所述表面所在平面的法向量方向相反时,用所述探头接触所述平板的一点,测量所述点的坐标,然后机器人末端坐标系的Y轴方向与所述表面所在平面的法向量方向相反时,用所述探头接触所述平板的一点,测量所述点的坐标,最后机器人末端坐标系的Z轴方向与所述表面所在平面的法向量方向相反时,用所述探头接触所述平板的一点,测量所述点的坐标;
步骤104:利用步骤103测量出的六个点的坐标计算出机器人TCP的TX、TY和TZ三个分量。
下面证明步骤102中机器人末端保持同样的姿态进行测量,平面法向量的精确值与TCP的初始值是没有关系的。
假设机器人的TCP的值是(Tx、Ty、Tz)。用接触式触发探头测量的测量点相对于机器人的底座的坐标取决于以下公式:
x=r11*Tx+r12*Ty+r13*Tz+Px
y=r21*Tx+r22*Ty+r23*Tz+Py    (1)
z=r31*Tx+r32*Ty+r33*Tz+Pz
其中(rij,i=1,2,3;j=1,2,3)和(Px,Py,Pz)表示末端tool0相对于机器人的底座的方向和位置,这些值可以从机器人控制器上实时读取。
确定一个平面,至少需要测量不在同一直线上的三个点,此处用三个点分析。三个以上的点所决定的平面是基于最小二乘法得到。
因此平板上的三个测量点的坐标取决于以下公式:
xi=r11*Tx+r12*Ty+r13*Tz+Pxi
yi=r21*Tx+r22*Ty+r23*Tz+Pyi    (2)
zi=r31*Tx+r32*Ty+r33*Tz+Pzi
i=1,2,3;
在这三个测量点的测量过程中,机器人一直保持着相同的运动方向,因此rij是不变的。
平面的法向量计算如下:
假设平面的方程是
A*x+B*y+C*z=1;    (3)
平面的参数A、B和C由测量得到的三个测量点的坐标(x1,y1,z1),(x2,y2,z2),(x3,y3,z3)决定。
令:
A = Δ 1 Δ
B = Δ 2 Δ - - - ( 4 )
C = Δ 3 Δ
Δ 1 = 1 , y 1 , z 1 1 , y 2 , z 2 1 , y 3 , z 3 , Δ 2 = x 1 , 1 , z 1 x 2 , 1 , z 2 x 3 , 1 , z 3 , Δ 3 = x 1 , y 1 , 1 x 2 , y 2 , 1 x 3 , y 3 , 1
则平面法向量(nx,ny,nz)和偏移量是:
n x = Δ 1 ( Δ 1 ) 2 + ( Δ 2 ) 2 + ( Δ 3 ) 2
n y = Δ 2 ( Δ 1 ) 2 + ( Δ 2 ) 2 + ( Δ 3 ) 2 - - - ( 5 )
n z = Δ 3 ( Δ 1 ) 2 + ( Δ 2 ) 2 + ( Δ 3 ) 2
d = Δ ( Δ 1 ) 2 + ( Δ 2 ) 2 + ( Δ 3 ) 2
又, Δ 1 = 1 , y 1 , z 1 1 , y 2 , z 2 1 , y 3 , z 3 = 1 , y 1 , z 1 0 , y 2 - y 1 , z 2 - z 1 0 , y 3 - y 1 , z 3 - z 1 = ( y 2 - y 1 ) * ( z 3 - z 1 ) - ( z 2 - z 1 ) * ( y 3 - y 1 ) (6)p表示的
= ( p y 2 - p y 1 ) * ( p z 3 - p z 1 ) - ( p z 2 - p z 1 ) * ( p y 3 - p y 1 ) ;
是机器人末端点(Tool0)的坐标。
上面的数学公式推导过程说明了Δ1和(Tx、Ty、Tz)无关,同理可以证明Δ2和Δ3与(Tx、Ty、Tz)也无关,所以(5)式中的平面法向量(nx,ny,nz)也与(Tx、Ty、Tz)无关。
步骤103的测量和步骤104的计算方法如下:
利用(1)~(6)式可以证明平面的法向量与TCP的准确值无关,然后将法向量计算出来,现设为(nx,ny,nz),
使用机器人控制程序(如ABB机器人的Rapid程序)使机器人接触平板上的一点,它的TCP位置由下式决定:
x=r11*Tx+r12*Ty+r13*Tz+Px
y=r21*Tx+r22*Ty+r23*Tz+Py    (7)
z=r31*Tx+r32*Ty+r33*Tz+Pz
其中(rrij,i=1.2.3,j=1,2,3)和(Px,Py,Pz)表示点末端tool0相对于机器人底部的方向和位置,这些值从机器人控制器上实时读取。偏移量就可以用以下的公式来计算:
d=nx*x+ny*y+nz*z
=nx*(r11*Tx+r12*Ty+r13*Tz+px)+ny*(r21*Tx+r22*Ty+r23*Tz+Py)+nz*(r31*Tx+r32*Ty+r33*Tz+Pz);    (8)
1、Tx的校准:
第一步:使机器人末端坐标系的X轴与平板平整表面所在平面的法向量方向相同,即:
(r11,r21,r31)=(nx,ny,nz)
然后用探头接触该表面上的一点,测量该点,
把上式代入式(8),由于旋转矩阵rij是正交矩阵,Ty和Tz将被消去,于是方程(8)变成:
d1=Tx+nx*px1+ny*py1+nz*pz1    (9)
第二步:使机器人末端坐标系的X轴与平板平整表面所在平面的法向量方向相反,即:(r11,r21,r31)=-(nx,ny,nz)
然后用探头接触该表面上的一点,测量该点。
于是方程(8)将变为
d2=-Tx+nx*px2+ny*py2+nz*pz2
d2=-Tx+nx*px2+ny*py2+nz*pz2    (10)
联立方程(9)和(10),得到:
T x = ( d 1 - d 2 ) - n x * ( p x 1 - p x 2 ) - n y * ( p y 1 - p y 2 ) - n z * ( p z 1 - p z 2 ) 2 - - - ( 11 )
因为两次测量用的是同一平面,则d1=d2,于是:
T x = - n x * ( p x 1 - p x 2 ) - n y * ( p y 1 - p y 2 ) - n z * ( p z 1 - p z 2 ) 2 - - - ( 12 )
(px1,py1,pz1)是第一个点相对于机器人底部的坐标,从机器人控制器上实时读取,(px2,py2,pz2)是第二个点相对于机器人底部的坐标,从机器人控制器上实时读取,(nx,ny,nz)是平面的法向量,均已知。
2、Ty的校准:
同理,Ty的校准和Tx的校准相同,
第一步:使机器人末端坐标系的Y轴与平板平整表面所在平面的法向量方向相同,即:
(r11,r21,r31)=(nx,ny,nz)
然后用探头接触该表面上的一点,测量该点,
第二步:使机器人末端坐标系的Y轴与平板平整表面所在平面的法向量方向相反,即:
(r11,r21,r31)=-(nx,ny,nz)
然后用探头接触该表面上的一点,测量该点,
则结果: T y = - n x * ( p x 1 - p x 2 ) - n y * ( p y 1 - p y 2 ) - n z * ( p z 1 - p z 2 ) 2
3、Tz的校准:
同理,Tz的校准也和Tx的校准是一样的,在测量中,分别使:
第一步:使机器人末端坐标系的Y轴与平板平整表面所在平面的法向量方向相同,即:
(r11,r21,r31)=(nx,ny,nz)
用探头接触该表面上的一点,然后测量该点,
第二步:使机器人末端坐标系的Y轴与平板平整表面所在平面的法向量方向相反,即:
(r11,r21,r31)=-(nx,ny,nz)
然后用探头接触该表面上的一点,测量该点,
则结果: T z = - n x * ( p x 1 - p x 2 ) - n y * ( p y 1 - p y 2 ) - n z * ( p z 1 - p z 2 ) 2
在机器人标定其中最后一个量(如本实施例TCP中的Tz)时,为避免机器人关节运动过大,引入过大机器人系统误差,可以通过仅控制工具末端点接触到其中一个面来计算。
计算方法如下:
以最后计算Tz为例,在标定完Tx和Ty后,由式(8)可知,平板表面方程的偏移量是Tz的函数,改变机器人的姿态,使工具末端再次接触同一平板表面点,得到平板表面方程的偏移量的另外一个表达式。消去两式中的d,可建立方程求解Tz.
实施例2:
如图3所示,为本发明实施例2的结构示意图,机器人1为六自由度工业机器人,底座10为机器人基座,tool0为机器人末端11的坐标系,一个表面平整的平板3固定在机器人末端上,该平板的两个表面为第一表面和第二表面,而且相互平行,地面放置探头2。如图4所示,为本发明机器人的固定TCP在机器人基座坐标系下的位置校准法的流程示意图,包括如下步骤:
步骤201:将测量用的探头安置在机器人的活动范围内;
步骤202:机器人的末端保持一定的姿态平动,测量所述第一表面上不在一条直线上的三个或三个以上的点的坐标,利用这些坐标计算所述表面所在平面的法向量;
可以证明:在对平面法向量的测量中,一旦机器人的末端保持一定的姿态进行测量,平面法向量的精确值与TCP的初始值是没有关系的。因此既然平面法向量的测量与未知的TCP值没有关系,那么法向量方向的测量它就是准确的。
步骤203:首先所述第一表面所在平面的法向量方向与机器人基座坐标系的X轴方向相同时,用所述探头接触所述第一表面上的一点,测量所述点的坐标,接着所述第一表面所在平面的法向量方向与机器人基座坐标系的X轴方向相反时,用所述探头接触所述第二表面上的一点,测量所述点的坐标;然后所述第一表面所在平面的法向量方向与机器人基座坐标系的Y轴方向相同时,用所述探头接触所述第一表面的一点,测量所述点的坐标,接着将所述第一表面所在平面的法向量方向与机器人基座坐标系的Y轴方向相反时,用所述探头接触所述第二表面上的一点,测量所述点的坐标;最后所述第一表面所在平面的法向量方向与机器人基座坐标系的Z轴方向相同时,用所述探头接触所述第一表面上的一点,测量所述点的坐标,接着所述第一表面所在平面的法向量方向与机器人基座坐标系的Z轴方向相反时,用所述探头接触所述第二表面上的一点,测量所述点的坐标;
步骤204:利用步骤203测量出的六个点的坐标计算出机器人TCP的TX、TY和TZ三个分量。
下面证明步骤202中机器人末端保持同样的姿态进行测量,平面法向量的精确值与TCP的初始值是没有关系的。
假设机器人的TCP的值是T=(Tx,Ty,Tz)T。当平板和探头接触时,有下列关系成立:
R0·Xt+T0=T
其中,R0和T0分别是Tool0对于基座的旋转和平移变换,Xt是探头探测的点在Tool0坐标系下的坐标.
让机器人基座保持同一姿态平行探测三个不在同一条直线上的点,则有如下关系成立:
R 0 · X t 1 + T 01 = T R 0 · X t 2 + T 02 = T R 0 · X t 3 + T 03 = T
这里R0,T01,T02,T03可以从机器人的控制器上读出,Xt1,Xt2,Xt3为探头在平板面上探测到的三点在Tool0坐标系下的位置。
将上三等式两两相减,可以得到:
R 0 · ( X t 1 - X t 2 ) = T 02 - T 01 R 0 · ( X t 1 - X t 3 ) = T 03 - T 01
如此可以得到关系平板的法向量N:
N = ( X t 1 - X t 2 ) × ( X t 1 - X t 3 ) | | ( X t 1 - X t 2 ) × ( X t 1 - X t 3 ) | | = R 0 - 1 ( ( T 02 - T 01 ) × ( T 03 - T 01 ) ) | | R 0 - 1 ( ( T 02 - T 01 ) × ( T 03 - T 01 ) ) | |
此关系式表明,平板的法向量与TCP的位置没有关系.其中Xt1,Xt2,Xt3确定一个平面,至少需要测量不在同一直线上的三个点,此处用三个点分析。三个以上的点所决定的平面是基于最小二乘法得到。
固定TCP在机器人基座坐标系下的位置的校准原理如下描述.
假设平板的表面所在的平面可表示为:
N·X=d
其中N为平面的单位法向量,X为tool0下的坐标.
Tool 0 = R 0 T 0 0 1 为机器人末端到基座的变换, Too l 0 - 1 = R 0 - 1 - R 0 - 1 · T 0 0 1
同时令R0=(r1,r2,r3)-1,则有关系:
N · ( Tool 0 - 1 · T 1 ) = N · ( R 0 - 1 · T - R 0 - 1 · T 0 ) = d
即有关系:
N · ( ( r 1 , r 2 , r 3 ) · T x T y T z - ( r 1 , r 2 , r 3 ) · T 0 x T 0 y T 0 z ) = d
移动机器人,让r1=N,可以得到:
Tx-T01x=d1
再移动机器人,让r1=-N,可以得到:
-Tx+T02x=d2
如此,可以计算出:
T x = T 01 x + T 02 x + ( d 1 - d 2 ) 2
同样的方法,我们可以得到:
T y = T 01 y + T 02 y + ( d 1 - d 2 ) 2
T z = T 01 z + T 02 z + ( d 1 - d 2 ) 2
由于d1-d2为平板的厚度,为已知量,并且T01x,T02x,T01y,T02y,Y01z,T02z可从机器人控制器上读出,故我们可以得到分离地得到固定TCP的三分量.
具体的实施时候,可以用机器人控制程序(如ABB机器人的Rapid程序)使机器人接触平板上的一点,先校准平板的法向量然后再校准TCP位置.
步骤203的测量和步骤204的计算方法如下:
1、平板法向量的校准
使用机器人控制程序(如ABB机器人的Rapid程序)使机器人末端持的平板的第一表面接触到探头,并让探头在平面上移动(注意机器人末端的姿态保持不变),采集到三个不在同一条直线上的三点的坐标,然后由以上计算平板平面法向量的方法计算出平板的单位法向量.
2、固定TCP三分量TX,TY,Tz的校准:
2.1Tx的校准
第一步:使机器人基座坐标系的X轴与平板第一表面所在平面的法向量方向相同,即:r1=N
然后用探头接触该表面上的一点,测量该点,则该点满足关系式:
Tx-T01x=d1
第二步:使机器人基座坐标系的X轴与平板第一表面所在平面的法向量方向相反,即r1=-N
然后用探头接触第二表面上的一点,测量该点,同样的该点也满足关系式:
-Tx+T02x=d2
由上面的两个等式,可以得到:
T x = T 01 x + T 02 x + ( d 1 - d 2 ) 2
由于d1-d2为平板厚度,为已知量,且可以从控制器上读出机器人末端的T01x,T02x,则可以完成Tx的校准
2.2Ty的校准:
Ty的校准和Tx的校准原理相同,具体的操作如下.
第一步:使机器人基座坐标系的Y轴与平板平整表面所在平面的法向量方向相同,即:r2=N
然后用探头接触该表面上的一点,测量该点,则该点满足关系式:
Ty-T01y=d1
第二步:使机器人基座坐标系的Y轴与平板平整表面所在平面的法向量方向相反,即r2=-N
然后用探头接触第二表面上的一点,测量该点,同样的该点也满足关系式:
-Ty+T02y=d2
由上面的两个等式,可以得到:
T y = T 01 y + T 02 y + ( d 1 - d 2 ) 2
由于d1-d2为平板厚度,为已知量,且可以从控制器上读出机器人末端的T01y,T02y,则可以完成Ty的校准.
2.3Tz的校准:
Tz的校准也和Tx的校准是一样的,在测量中,分别使:
第一步:使机器人基座坐标系的Y轴与平板平整表面所在平面的法向量方向相同,即:r3=N
然后用探头接触该表面上的一点,测量该点,则该点满足关系式:
Tz-T01z=d1
第二步:使机器人基座坐标系的Y轴与平板平整表面所在平面的法向量方向相反,即r3=-N
然后用探头接触第二表面上的一点,测量该点,同样的该点也满足关系式:
-Tz+T02z=d2
由上面的两个等式,可以得到:
T z = T 01 z + T 02 z + ( d 1 - d 2 ) 2
由于d1-d2为平板厚度,为已知量,且可以从控制器上读出机器人末端的T01z,T02z,则可以完成Tz的校准.
在机器人标定其中最后一个量(如本实施例TCP中的Tz)时,为避免机器人关节运动过大,引入过大机器人系统误差,可以通过仅控制工具末端点接触到其中一个面来计算。
计算方法如下:
以最后计算Tz为例,在标定完Tx和Ty后,可知平板表面方程的偏移量d是Tz的函数:
d = N · ( R 0 - 1 · T x T y T z - R 0 - 1 · T 0 )
其中,N为平面的法向量,R0,T0为机器人末端相对于基座的旋转和平移关系。
改变机器人的姿态,使工具末端再次接触同一平板表面点,得到平板表面方程的偏移量的另外一个表达式。消去两式中的d,可建立方程求解Tz.
因此,本发明将TCP所有的变量都单独分开来计算,每次测量只计算其中一个变量。由于所有的变量都被单独分离开,所以很容易控制测量过程中校准的精确程度。
最后所应说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围。

Claims (10)

1.一种基于平板测量的机器人工具中心点在机器人末端坐标系下的位置校准法,机器人末端固定一个测量用探头,一个表面平整且厚度已知的平板置于机器人的活动范围之内,本方法包括以下步骤:
步骤1:所述机器人末端保持一定的姿态平动,用所述探头接触所述平板上不在一条直线上的三个或三个以上的点,得到这些点的坐标和位置,利用这些坐标和位置计算所述表面所在平面的法向量;
步骤2:机器人末端坐标系的X、Y和Z轴方向分别与所述表面所在平面的法向量方向相同时,分别用所述探头接触所述平板的一点,测量所述点的坐标;然后与法向量方向相反,分别用所述探头接触所述平板的一点,测量所述点的坐标;
步骤3:利用步骤1得出的平面法向量和步骤2测量出的所有点的坐标计算出机器人工具中心点在机器人末端坐标系下的位置。
2.根据权利要求1所述的基于平板测量的机器人工具中心点在机器人末端坐标系下的位置校准法,其特征在于:所述步骤2具体为,首先机器人末端坐标系的X轴方向与所述表面所在平面的法向量方向相同时,用所述探头接触所述平板的一点,测量所述点的坐标,然后机器人末端坐标系的Y轴方向与所述表面所在平面的法向量方向相同时,用所述探头接触所述平板的一点,测量所述点的坐标,接着机器人末端坐标系的Z轴方向与所述表面所在平面的法向量方向相同时,用所述探头接触所述平板的一点,测量所述点的坐标;再后机器人末端坐标系的X轴方向与所述表面所在平面的法向量方向相反时,用所述探头接触所述平板的一点,测量所述点的坐标,然后机器人末端坐标系的Y轴方向与所述表面所在平面的法向量方向相反时,用所述探头接触所述平板的一点,测量所述点的坐标,最后机器人末端坐标系的Z轴方向与所述表面所在平面的法向量方向相反时,用所述探头接触所述平板的一点,测量所述点的坐标。
3.根据权利要求1或2所述的基于平板测量的机器人工具中心点在机器人末端坐标系下的位置校准法,其特征在于:所述步骤中的机器人为六自由度的机器人。
4.根据权利要求1或2所述的基于平板测量的机器人工具中心点在机器人末端坐标系下的位置校准法,其特征在于:所述步骤2中,机器人改变姿态,用所述探头接触所述平板,测量所述点的坐标。
5.根据权利要求1或2所述的基于平板测量的机器人工具中心点在机器人末端坐标系下的位置校准法,其特征在于:所述步骤中通过读取机器人控制器而得到所述点的坐标和位置。
6.一种基于平板测量的机器人工具中心点在机器人基座坐标系下的位置校准法,机器人末端固定一个平板,该平板的两个表面均为光滑且平行,为第一表面和第二表面,一个测量用的探头固定于机器人的活动范围之内,本方法包括以下步骤:
步骤1:机器人末端保持一定的姿态平动,使得探头接触平板的第一表面上不在一条直线上的三个或三个以上的点,得到这些点的坐标和位置,利用这些坐标和位置计算该第一表面的法向量;
步骤2:机器人基座坐标系的X、Y和Z轴方向分别与所述第一表面所在平面的法向量方向相同时,分别用所述探头接触所述第一表面的一点,测量所述点的坐标;然后机器人基座坐标系的X、Y和Z轴方向分别与所述第一表面所在平面的法向量方向相反,分别用所述探头接触所述第二表面的一点,分别测量所述点的坐标;
步骤3:利用平板厚度值和步骤2测量出的所有点的坐标计算出机器人工具中心点在机器人基座坐标系下的位置。
7.根据权利要求6所述的基于平板测量的机器人工具中心点在机器人基座坐标系下的位置校准法,其特征在于:所述步骤2具体为,首先所述第一表面所在平面的法向量方向与机器人基座坐标系的X轴方向相同时,用所述探头接触所述第一表面上的一点,测量所述点的坐标,接着所述第一表面所在平面的法向量方向与机器人基座坐标系的X轴方向相反时,用所述探头接触所述第二表面上的一点,测量所述点的坐标;然后所述第一表面所在平面的法向量方向与机器人基座坐标系的Y轴方向相同时,用所述探头接触所述第一表面的一点,测量所述点的坐标,接着将所述第一表面所在平面的法向量方向与机器人基座坐标系的Y轴方向相反时,用所述探头接触所述第二表面上的一点,测量所述点的坐标;最后所述第一表面所在平面的法向量方向与机器人基座坐标系的Z轴方向相同时,用所述探头接触所述第一表面上的一点,测量所述点的坐标,接着所述第一表面所在平面的法向量方向与机器人基座坐标系的Z轴方向相反时,用所述探头接触所述第二表面上的一点,测量所述点的坐标。
8.根据权利要求6或7所述的基于平板测量的机器人工具中心点在机器人基座坐标系下的位置校准法,其特征在于:所述步骤中的机器人为六自由度的机器人。
9.根据权利要求6或7所述的基于平板测量的机器人工具中心点在机器人基座坐标系下的位置校准法,其特征在于:所述步骤2中,机器人改变姿态,用所述探头接触所述平板,测量所述点的坐标。
10.根据权利要求6或7所述的基于平板测量的机器人工具中心点在机器人基座坐标系下的位置校准法,其特征在于:所述步骤中通过读取机器人控制器而得到所述点的坐标和位置。
CNB200510087041XA 2005-07-25 2005-07-25 基于平板测量的机器人工具中心点三分量校准法 Expired - Fee Related CN100398274C (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB200510087041XA CN100398274C (zh) 2005-07-25 2005-07-25 基于平板测量的机器人工具中心点三分量校准法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB200510087041XA CN100398274C (zh) 2005-07-25 2005-07-25 基于平板测量的机器人工具中心点三分量校准法

Publications (2)

Publication Number Publication Date
CN1903525A CN1903525A (zh) 2007-01-31
CN100398274C true CN100398274C (zh) 2008-07-02

Family

ID=37672958

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB200510087041XA Expired - Fee Related CN100398274C (zh) 2005-07-25 2005-07-25 基于平板测量的机器人工具中心点三分量校准法

Country Status (1)

Country Link
CN (1) CN100398274C (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10926414B2 (en) 2017-09-29 2021-02-23 Industrial Technology Research Institute System and method for calibrating tool center point of robot
US11247340B2 (en) 2018-12-19 2022-02-15 Industrial Technology Research Institute Method and apparatus of non-contact tool center point calibration for a mechanical arm, and a mechanical arm system with said calibration function

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102909728B (zh) * 2011-08-05 2015-11-25 鸿富锦精密工业(深圳)有限公司 机器人工具中心点的视觉校正方法
CN103376080B (zh) * 2012-04-27 2016-09-07 上海西门子燃气轮机部件有限公司 一种用于校准机器人工作单元的方法
CN103021898B (zh) * 2012-12-17 2016-03-02 华中科技大学 测量芯片与基板相对倾角测量方法及系统
CN105157567B (zh) * 2015-05-15 2017-10-10 天津智通机器人有限公司 一种测量机器人的工具坐标系标定方法及系统
DE102016223841A1 (de) * 2016-11-30 2018-05-30 Siemens Healthcare Gmbh Berechnen eines Kalibrierungsparameters eines Roboterwerkzeugs
CN109579805B (zh) * 2018-11-26 2020-06-23 成都经纬时空科技有限公司 一种基线自校准测量方法
CN109900207A (zh) * 2019-03-12 2019-06-18 精诚工科汽车系统有限公司 机器人视觉工具的工具中心点标定方法及系统
CN110375688A (zh) * 2019-06-18 2019-10-25 宁波敏实汽车零部件技术研发有限公司 一种工业机器人工具坐标系姿态标定系统与方法
CN110794223A (zh) * 2019-11-13 2020-02-14 中国船舶重工集团公司第七二四研究所 基于工业机器人的数字多通道雷达天线标校系统及方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1033590A (zh) * 1987-07-28 1989-07-05 布莱阿姆青年大学 机器人误差矫正装置和方法
US5162713A (en) * 1989-10-20 1992-11-10 Hitachi, Ltd. Structural error correction method for SCARA robot
US5910719A (en) * 1996-09-17 1999-06-08 Cycle Time Corporation Tool center point calibration for spot welding guns
US5929584A (en) * 1998-04-27 1999-07-27 Pht, Inc. Tool center point calibration apparatus and method
US6070109A (en) * 1998-03-10 2000-05-30 Fanuc Robotics North America, Inc. Robot calibration system

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1033590A (zh) * 1987-07-28 1989-07-05 布莱阿姆青年大学 机器人误差矫正装置和方法
US5162713A (en) * 1989-10-20 1992-11-10 Hitachi, Ltd. Structural error correction method for SCARA robot
US5910719A (en) * 1996-09-17 1999-06-08 Cycle Time Corporation Tool center point calibration for spot welding guns
US6070109A (en) * 1998-03-10 2000-05-30 Fanuc Robotics North America, Inc. Robot calibration system
US5929584A (en) * 1998-04-27 1999-07-27 Pht, Inc. Tool center point calibration apparatus and method

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10926414B2 (en) 2017-09-29 2021-02-23 Industrial Technology Research Institute System and method for calibrating tool center point of robot
US11247340B2 (en) 2018-12-19 2022-02-15 Industrial Technology Research Institute Method and apparatus of non-contact tool center point calibration for a mechanical arm, and a mechanical arm system with said calibration function

Also Published As

Publication number Publication date
CN1903525A (zh) 2007-01-31

Similar Documents

Publication Publication Date Title
CN100398274C (zh) 基于平板测量的机器人工具中心点三分量校准法
CN107042528B (zh) 一种工业机器人的运动学标定系统及方法
Kunzmann et al. A uniform concept for calibration, acceptance test, and periodic inspection of coordinate measuring machines using reference objects
US9542355B2 (en) Method for recalibrating coordinate positioning apparatus
US7079969B2 (en) Dynamic artefact comparison
JP3005681B1 (ja) Cmm校正ゲージ及びcmmの校正方法
US7131207B2 (en) Workpiece inspection method
Li et al. An accurate probe pre-travel error compensation model for five-axis on-machine inspection system
JP2020046301A (ja) 工作機械の計測誤差評価方法及びプログラム
CN103439085B (zh) 一种接触式测量曲面棱镜参数的方法及装置
Forbes et al. Uncertainty associated with coordinate measurement in comparator mode
KR102035334B1 (ko) 4축 공작기계의 기하학적 오차 측정 방법
Furutani et al. Parameter calibration for non-cartesian CMM
NL2005591C2 (en) Apparatus and method for calibrating a coordinate measuring apparatus.
CN102445279B (zh) 一种测量干涉仪波长的装置及方法
Gąska et al. Identification and correction of coordinate measuring machine geometrical errors using lasertracer systems
CN108375337B (zh) 一种机器人及其工艺设备相对位姿的测定方法及测定装置
CN109062138A (zh) 一种基于立体标定块的五轴平台系统标定方案
Clarke et al. The case for a consistent method of verifying the performance of large volume metrology systems
RU2753150C1 (ru) Способ скалярной калибровки блока акселерометров
Laaouina et al. Application of CMM separation method for identifying absolute values of probe errors and machine errors
EP3745223B1 (en) Method for self-verification of mechatronic systems
Dueso et al. Thermal Error Compensation Strategy in Simultaneous Laser Multilateration for Volumetric Verification
De Vicente et al. Simplified statistical method for uncertainty estimation in coordinate metrology
Brosed Dueso et al. Thermal Error Compensation Strategy in Simultaneous Laser Multilateration for Volumetric Verification

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20180824

Address after: 300000 Tianjin Dongli District Tianjin free trade pilot area (Airport Economic Zone) Aviation Road 278 office building, block C.

Patentee after: TIANJIN ZHITONG INFORMATION SYSTEM INTEGRATION Co.,Ltd.

Address before: 065001 30 hongro Road, Langfang Development Zone, Hebei

Patentee before: INTERSMART ROBOTIC SYSTEMS Co.,Ltd.

TR01 Transfer of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20080702

CF01 Termination of patent right due to non-payment of annual fee