CN102230783B - 一种用于工业机器人的空间立体网格精度补偿方法 - Google Patents

一种用于工业机器人的空间立体网格精度补偿方法 Download PDF

Info

Publication number
CN102230783B
CN102230783B CN201110113246A CN201110113246A CN102230783B CN 102230783 B CN102230783 B CN 102230783B CN 201110113246 A CN201110113246 A CN 201110113246A CN 201110113246 A CN201110113246 A CN 201110113246A CN 102230783 B CN102230783 B CN 102230783B
Authority
CN
China
Prior art keywords
robot
cube grid
industrial robot
cube
space
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
CN201110113246A
Other languages
English (en)
Other versions
CN102230783A (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 University of Aeronautics and Astronautics
Chengdu Aircraft Industrial Group Co Ltd
Original Assignee
Nanjing University of Aeronautics and Astronautics
Chengdu Aircraft Industrial Group 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 University of Aeronautics and Astronautics, Chengdu Aircraft Industrial Group Co Ltd filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201110113246A priority Critical patent/CN102230783B/zh
Publication of CN102230783A publication Critical patent/CN102230783A/zh
Application granted granted Critical
Publication of CN102230783B publication Critical patent/CN102230783B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Numerical Control (AREA)

Abstract

本发明公开了一种用于工业机器人的空间立体网格精度补偿方法,属于工业机器人逆标定技术领域。该方法利用工业机器人具有较高重复定位精度的特性,采用激光跟踪仪来建立理论坐标与实际定位坐标之间的关系,对于包络空间内划分的某个立方体网格内的任一点,采用空间插值法来对机器人的理论坐标进行修正,完成机器人在该点的绝对定位精度补偿。本发明方法计算过程简单迅速,能显著提高机器人的绝对定位精度,使机器人能适应更广泛的应用场合。

Description

一种用于工业机器人的空间立体网格精度补偿方法
技术领域
本发明涉及一种机器人定位补偿方法,尤其涉及一种用于工业机器人的空间立体网格精度补偿方法,属于工业机器人逆标定技术领域。
背景技术
机器人的精度是反映机器人性能的一个重要指标,它包括绝对定位精度和重复定位精度。绝对定位精度误差是机器人实际运动与期望运动之间的偏差,由确定性原始误差(如连杆参数误差、运动副间隙等)所产生;重复定位精度误差是机器人重复执行同一期望运动时,机器人的实际运动之间的相互离散程度,由随机性原始误差(如关节伺服定位误差等)所产生。
通常,一般的工业机器人具有较高的重复定位精度,然而它的绝对定位精度却比较差,重复定位精度达到0.1mm的机器人,其绝对定位精度误差却能有2-3mm。示教利用了机器人的重复定位精度可以达到较高的绝对定位精度,但在实际应用中,很多场合对工业机器人的绝对定位精度要求很高,通常又不能通过示教的方式来完成。所以,对机器人的绝对定位精度进行补偿就显得至关重要。
文献“夏凯,陈崇端,洪涛,等.补偿机器人定位误差的神经网络[J].机器人,1995,17(3):171-176.”为了提高机器人的定位精度,提出了补偿机器人定位误差的神经网络方法。文章针对RM-501五自由度机器人,通过两种方式用多层感知器神经网络补偿机器人运动学方程定位的误差,分别是基于关节坐标的神经网络补偿和基于直角坐标位置的神经网络补偿。该方法的核心是利用人工神经网络具有很强的自学习、自适应能力,通过训练得到机器人运动学几何参数误差、非几何参数误差等误差源的作用规律,对期望关节角或直角坐标进行补偿,从而来提高机器人的绝对定位精度。该方法避免了其他传统标定方法繁琐的建模及参数辨识过程,但在实际应用中存在以下不足:
1)为了使得训练后的网络能够达到一定的精度以及适应机器人包络范围内的所有点,训练神经网络需要有大量的学习样本,因此测量工作量大;
2)该方法需要将直角坐标转换成关节坐标,因此需要对机器人运动学方程求逆解,这个过程计算量大,且在奇异点附近效果不理想;
3)试验结果表明机器人在标定后的定位精度仍不够理想。
发明内容
本发明为提高工业机器人的绝对定位精度,针对现有技术存在的不足,而提出一种用于工业机器人的空间立体网格精度补偿方法。
该方法包括如下步骤:
步骤1:在工业机器人的包络线范围内,按给定的最大步长把整个包络空间划分成一系列的立方体网格;
步骤2:通过激光跟踪仪测量并建立机器人基坐标系,用步骤1中划分的每个立方体网格的八个顶点的理论坐标来控制机器人进行定位,并用激光跟踪仪测量并记录实际定位坐标;
步骤3:对于期望到达包络线范围内的任一点P的理论坐标(X,Y,Z),
1)查找该点P所在的立方体网格;
2)计算该点P与所在立方体网格的八个顶点Ki的距离di,用算得的距离di进行反距离加权求得八个顶点Ki相对于该点P的权值qi,其中i=1,2,...8,下同;
3)用求得的权值qi来对八个顶点Ki的X、Y、Z三个方向上的定位误差分别进行空间插值,预算出该点P三个方向上的误差;
4)用求得的误差对该点P的理论坐标(X,Y,Z)进行反向修正,完成机器人在该点P的定位精度补偿。
本发明具有如下技术效果:
i)针对不同型号的工业机器人,确定了划分网格的最大步长,减少了工作量,有利于工业机器人快速地投入应用。
ii)本方法是在笛卡尔坐标系中进行的,与通常的机器人标定方法相比不需要进行机器人运动学的正解和逆解,计算过程简单迅速,可以实现在线补偿。
iii)可以显著提高机器人的绝对定位精度,使得标定后的工业机器人能适应更广泛的应用场合。
附图说明
图1为本发明补偿方法的原理示意图。
图2为本发明补偿方法中的空间插值示意图。
具体实施方式
如图1所示,本发明用于工业机器人的空间立体网格精度补偿方法的步骤如下:
步骤1:在工业机器人的包络线范围内,在笛卡尔坐标系中按给定的最大步长把整个包络空间划分成一系列的立方体网格;
步骤2:通过激光跟踪仪测量并建立机器人基坐标系,用步骤1中划分的每个立方体网格的八个顶点的理论坐标来控制机器人进行定位,机器人在定位立方体网格的八个顶点处的姿态一致或偏差在±10°范围内,然后用激光跟踪仪测量并记录实际定位坐标;
建立激光跟踪仪与机器人基坐标系之间关联的步骤为:
①将球形固定反射器SMR固定在末端执行器的TCP上,且保持A2到A6轴的位置(角度)固定不变,通过旋转A1轴,用FARO激光跟踪仪测量一系列位于圆周上的点;
②利用FARO自带的CAM2Measure软件,将步骤①得到的测量点依次拟合出一个平面和一个圆,从而得到圆心的理论坐标;
③测量机器人底座平面上的一系列点并拟合出一个平面,再对该平面做偏移量为SMR半径的偏移;
④把步骤②中得到的圆心投影到步骤③偏移得到的平面中,得到机器人的原点位置;
⑤测量机器人处于零点位置时机器人法兰盘上两个水平对称的安装孔,然后算出这两个测量点的中点,再把所得中点向步骤③中得到的偏移平面做投影,得到处于X轴上的一点;
⑥从机器人机械零点开始按逆时针方向旋转A1轴(小于90度的任意值)并测量,再把该点向步骤③中得到的偏移平面做投影,得到处于+XY平面上的一点;
⑦利用原点以及步骤⑤和⑥得到的两点构造出坐标系,该坐标系即为机器人坐标系。
步骤3:对于期望到达包络线范围内的任一点P的理论坐标(X,Y,Z),
1)查找该点P所在的立方体网格;
2)计算该点P与所在立方体网格的八个顶点Ki的距离di,用算得的距离di进行反距离加权求得八个顶点Ki相对于该点P的权值qi,其中i=1,2,...8,下同;
3)用求得的权值qi来对八个顶点Ki的X、Y、Z三个方向上的定位误差分别进行空间插值,预算出该点P三个方向上的误差;
4)用求得的误差对该点P的理论坐标(X,Y,Z)进行反向修正,完成机器人在该点P的定位精度补偿。
所述步骤3中,机器人在立方体网格中任一点的姿态与它在相应立方体网格八个顶点处的姿态保持一致,或允许偏差在±10°范围内。
所述步骤3中的计算如下:
d i = ( X - X i ′ ) 2 + ( Y - Y i ′ ) 2 + ( Z - Z i ′ ) 2 , i = 1,2 , . . . 8 - - - ( 1 )
q i = 1 d i 1 d 1 + 1 d 2 + 1 d 3 + 1 d 4 + 1 d 5 + 1 d 6 + 1 d 7 + 1 d 8 , i = 1,2 , . . . 8 - - - ( 2 )
ΔX = Σ i = 1 8 Δ X i q i ; ΔY = Σ i = 1 8 Δ Y i q i ; ΔZ = Σ i = 1 8 Δ Z i q i - - - ( 3 )
X=X+ΔX;Y=Y+ΔY;Z=Z+ΔZ                   (4)
式中:(X,Y,Z)是点P的理论坐标;(Xi′,Yi′,Zi′)为包含点P的立方体网格八个顶点的实际定位坐标;(ΔXi,ΔYi,ΔZi)为包含点P的立方体网格八个顶点的实际坐标与理论坐标偏差;(ΔX,ΔY,ΔZ)为(X,Y,Z)的修正量;(X,Y,Z)为(X,Y,Z)修正后的理论坐标。
下面以KUKA150-2机器人为例来说明本发明的具体实施步骤。
步骤1:在工业机器人的包络线范围内,按给定的最大步长lmax把机器人的工作空间划分成一系列的立方体网格;
这里在沿着机器人坐标系X方向上选取其中5个立方体网格,它们的八个顶点的理论坐标如下表所示:
Figure BDA0000059012540000044
Figure BDA0000059012540000051
步骤2:通过激光跟踪仪测量并建立机器人基坐标系,用步骤1中划分的每个立方体网格的八个顶点的理论坐标来控制机器人进行定位,并用激光跟踪仪测量并记录实际定位坐标;
经过测量,步骤1中选取的5个立方体网格的理论坐标和实际定位坐标如下表所示:
Figure BDA0000059012540000061
步骤3:对于期望到达包络线范围内的任一点P的理论坐标(X,Y,Z),如图2所示,
1)查找该点P所在的立方体网格;
2)计算该点P与所在立方体网格的八个顶点Ki的距离di,用算得的距离di进行反距离加权求得八个顶点Ki相对于该点P的权值qi,其中i=1,2,...8,下同;
3)用求得的权值qi来对八个顶点Ki的X、Y、Z三个方向上的定位误差分别进行空间插值,预算出该点P三个方向上的误差;
4)用求得的误差对该点P的理论坐标(X,Y,Z)进行反向修正,完成机器人在该点P的定位精度补偿。
步骤3中涉及的计算见式(1)~(4)。
对于所选取的5个立方体网格,分别在其中任取一点,进行空间插值求得的权值、修正量以及最终的补偿效果如以下各表所示:
Figure BDA0000059012540000071
计算5个任选点相应于所处立方体网格八个顶点的权值,结果如下表所示:
  q1   0.1141637   0.08445   0.343196   0.099888   0.149655
  q2   0.0831362   0.109828   0.082428   0.148499   0.099179
  q3   0.3404843   0.109224   0.113861   0.099998   0.149682
  q4   0.114108   0.195453   0.11303   0.14922   0.099301
  q5   0.0829874   0.08461   0.113981   0.100482   0.151292
  q6   0.0684113   0.110074   0.082584   0.151229   0.099765
  q7   0.113865   0.109563   0.082835   0.100337   0.151268
  q8   0.0828442   0.1968   0.06809   0.150347   0.099858
用所得权值计算得到的5个任选点理论坐标修正量△X,△Y,△Z分别为:
  △X   △Y   △Z
  -0.012784443   -0.106889088   -0.237533551
  -0.96534   0.488261   0.757579
  -1.41879   0.145454   0.394456
  -1.832755509   0.668394075   1.815461812
  -1.4033791   0.16524051   1.91628559
修正后的5个任选点的坐标为:
  X   Y   Z
  999.9872156   -200.1068891   2799.762466
  1299.034659   -699.511739   2400.757579
  1598.58121   100.145454   2300.39446
  2098.167244   0.668394075   2101.815462
  2198.59662   300.165241   1501.91629
用修正后的坐标控制机器人进行定位得到实际定位坐标,用它与期望坐标相比较,结果如下表所示:
Figure BDA0000059012540000072
Figure BDA0000059012540000081
通过对比结果可知,本发明提出的用于工业机器人的空间立体网格精度补偿方法可显著提高工业机器人的绝对定位精度。

Claims (4)

1.一种用于工业机器人的空间立方体网格精度补偿方法,其特征在于:
该方法包括如下步骤:
步骤1:在工业机器人的包络线范围内,按给定的最大步长把整个包络空间划分成一系列的立方体网格;
步骤2:通过激光跟踪仪测量并建立机器人基坐标系,用步骤1中划分的每个立方体网格的八个顶点的理论坐标来控制机器人进行定位,并用激光跟踪仪测量并记录实际定位坐标;
步骤3:对于期望到达包络线范围内的任一点P的理论坐标(X,Y,Z),
1)查找该点P所在的立方体网格;
2)计算该点P与所在立方体网格的八个顶点Ki的距离di,用算得的距离di进行反距离加权求得八个顶点Ki相对于该点P的权值qi,其中i=1,2,...8,下同;
3)用求得的权值qi来对八个顶点Ki的X、Y、Z三个方向上的定位误差分别进行空间插值,预算出该点P三个方向上的误差;
4)用求得的误差对该点P的理论坐标(X,Y,Z)进行反向修正,完成机器人在该点P的定位精度补偿。
2.根据权利要求1所述的用于工业机器人的空间立方体网格精度补偿方法,其特征在于:所述步骤1中立方体网格的划分是在笛卡尔坐标系中进行的。
3.根据权利要求1所述的用于工业机器人的空间立方体网格精度补偿方法,其特征在于:所述步骤2中在机器人进行定位时,机器人在定位立方体网格的八个顶点处的姿态一致或偏差在±10°范围内。
4.根据权利要求1所述的用于工业机器人的空间立方体网格精度补偿方法,其特征在于:所述步骤3中,机器人在立方体网格中任一点的姿态与它在相应立方体网格八个顶点处的姿态一致或偏差在±10°范围内。
CN201110113246A 2011-05-04 2011-05-04 一种用于工业机器人的空间立体网格精度补偿方法 Active CN102230783B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201110113246A CN102230783B (zh) 2011-05-04 2011-05-04 一种用于工业机器人的空间立体网格精度补偿方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201110113246A CN102230783B (zh) 2011-05-04 2011-05-04 一种用于工业机器人的空间立体网格精度补偿方法

Publications (2)

Publication Number Publication Date
CN102230783A CN102230783A (zh) 2011-11-02
CN102230783B true CN102230783B (zh) 2012-09-26

Family

ID=44843374

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201110113246A Active CN102230783B (zh) 2011-05-04 2011-05-04 一种用于工业机器人的空间立体网格精度补偿方法

Country Status (1)

Country Link
CN (1) CN102230783B (zh)

Families Citing this family (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102607552B (zh) * 2012-01-11 2014-12-10 南京航空航天大学 基于神经网络的工业机器人空间网格精度补偿方法
CN102693296B (zh) * 2012-05-16 2014-07-30 南京信息工程大学 一种大批量二维点数据中坐标快速匹配的方法
CN103144109B (zh) * 2013-02-25 2015-03-11 南京航空航天大学 一种用于附加外部轴的机器人***分站式精度补偿方法
CN104535027B (zh) * 2014-12-18 2017-06-13 南京航空航天大学 一种变参数误差辨识的机器人精度补偿方法
CN105066831A (zh) * 2015-09-09 2015-11-18 大族激光科技产业集团股份有限公司 一种单一或多机器人***协同工作坐标系的标定方法
CN106001927B (zh) * 2016-07-05 2018-03-23 温州大学激光与光电智能制造研究院 一种测量加工一体化的激光平整化抛光方法
CN107796276B (zh) * 2017-09-12 2020-05-15 武汉工程大学 一种估算工业机器人绝对定位精度的装置及方法
CN109895082A (zh) * 2017-12-07 2019-06-18 中国科学院沈阳自动化研究所 一种应用于航天装配设备的控制***
CN108489392B (zh) * 2018-04-11 2019-10-08 上海交通大学 机械臂末端重复定位精度测量装置与方法
CN108943024A (zh) * 2018-08-16 2018-12-07 北京航空航天大学 一种六自由度工业机器人绝对位置误差估计方法
CN110161850B (zh) * 2019-04-24 2020-04-07 南京航空航天大学 一种工业机器人变参数刚度辨识与建模方法
CN111055273B (zh) * 2019-05-16 2022-10-11 广西大学 一种用于机器人的两步误差补偿方法
CN110146017B (zh) * 2019-06-08 2021-02-02 青岛科技大学 工业机器人重复定位精度测量方法
CN110345867A (zh) * 2019-07-15 2019-10-18 北京工业大学 一种基于反距离权重法修正cmm空间任意点坐标误差的方法
CN110421566B (zh) * 2019-08-08 2020-10-27 华东交通大学 一种基于近似度加权平均插值法的机器人精度补偿方法
CN111390968B (zh) * 2019-12-31 2022-05-10 宁波大正工业机器人技术有限公司 一种料盘及其点位补偿方法
CN111590569A (zh) * 2020-05-14 2020-08-28 河北工业大学 一种工业机器人定位误差补偿方法
CN111633645A (zh) * 2020-05-15 2020-09-08 成都飞机工业(集团)有限责任公司 一种用于移动式机器人***的精度补偿方法
CN112720480B (zh) * 2020-12-22 2022-08-16 华中科技大学 一种基于分级误差的机器人轨迹修正方法及***
CN113427487B (zh) * 2021-07-09 2022-03-25 华南理工大学 一种基于电磁波测距的dh参数校准方法及***
CN114740801B (zh) * 2022-03-21 2023-09-29 成都飞机工业(集团)有限责任公司 一种用于数控设备群协同生产线安装的基坐标系创建方法
CN114932542B (zh) * 2022-06-20 2024-02-23 昆明理工大学 一种工业机器人距离误差补偿方法、***
CN116038717A (zh) * 2023-03-31 2023-05-02 珞石(北京)科技有限公司 光伏插片中机械臂快速定位的方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4530155A (en) * 1983-03-26 1985-07-23 Dr. Johannes Heidenhain Gmbh Method and apparatus for reproducing a reference position in an incremental measuring device
CN1994689A (zh) * 2005-12-28 2007-07-11 松下电器产业株式会社 机器人及机器人检测自动化方法
CN101077578A (zh) * 2007-07-03 2007-11-28 北京控制工程研究所 一种基于二元环境信息的移动机器人局部路径规划方法
JP2009233776A (ja) * 2008-03-26 2009-10-15 Yaskawa Electric Corp ロボット及びロボットの接触検出方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4530155A (en) * 1983-03-26 1985-07-23 Dr. Johannes Heidenhain Gmbh Method and apparatus for reproducing a reference position in an incremental measuring device
CN1994689A (zh) * 2005-12-28 2007-07-11 松下电器产业株式会社 机器人及机器人检测自动化方法
CN101077578A (zh) * 2007-07-03 2007-11-28 北京控制工程研究所 一种基于二元环境信息的移动机器人局部路径规划方法
JP2009233776A (ja) * 2008-03-26 2009-10-15 Yaskawa Electric Corp ロボット及びロボットの接触検出方法

Also Published As

Publication number Publication date
CN102230783A (zh) 2011-11-02

Similar Documents

Publication Publication Date Title
CN102230783B (zh) 一种用于工业机器人的空间立体网格精度补偿方法
CN102607552B (zh) 基于神经网络的工业机器人空间网格精度补偿方法
US9658055B2 (en) Accuracy traceability method based on precision coordinate control network for workshop measurement positioning system
CN110193829B (zh) 一种耦合运动学与刚度参数辨识的机器人精度控制方法
CN106625573B (zh) 一种五自由度混联机器人直接误差补偿技术
CN103218475B (zh) 一种基于在机测评的复杂空间型面误差反馈补偿方法
CN104237849B (zh) 一种双五元十字阵被动声定位融合方法
CN108705531A (zh) 工业机器人的运动学参数标定方法、标定***、电子设备
CN103389038A (zh) 激光跟踪仪定目标多站测量的数控机床几何精度检测方法
CN108908335B (zh) 基于改进差分进化算法的机器人标定方法
CN106199519A (zh) 一种超短基线五基元立体空间基阵及其水声定位方法
CN114161425B (zh) 一种工业机器人的误差补偿方法
CN107609227B (zh) 一种基于最大熵理论的装配工艺优化方法
CN103144109B (zh) 一种用于附加外部轴的机器人***分站式精度补偿方法
CN104977816B (zh) 基于Compact Particle Swarm Optimization算法的光刻机掩模台微动台的机械参数软测量方法
CN109520466A (zh) 一种基于全盘位测量的自由设站方法
TWI506243B (zh) 探針校準路徑類比系統及方法
CN111891176A (zh) 基于轨道几何状态测量仪的轨道测量数据整体平差方法
CN112069612A (zh) 一种齿轮测量中心测量不确定度评定方法
CN112720480B (zh) 一种基于分级误差的机器人轨迹修正方法及***
CN110967778B (zh) 一种动态坐标系多面体剖分重力布格校正方法
CN104834790A (zh) 一种空间尺寸链公差分析方法
Milovzorov et al. Angular installation options errors correction for three-component vector-measuring transducers with accelerometer at calibration phase
CN115415853A (zh) 一种五轴数控机床摆头几何误差辨识方法及***
CN110087307A (zh) 基于测距修正的水下传感器网络定位方法

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