CN105426627B - 机器人工装安装偏差计算系统及其应用 - Google Patents

机器人工装安装偏差计算系统及其应用 Download PDF

Info

Publication number
CN105426627B
CN105426627B CN201510907608.7A CN201510907608A CN105426627B CN 105426627 B CN105426627 B CN 105426627B CN 201510907608 A CN201510907608 A CN 201510907608A CN 105426627 B CN105426627 B CN 105426627B
Authority
CN
China
Prior art keywords
coordinate system
robot
laser
angle
product
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
CN201510907608.7A
Other languages
English (en)
Other versions
CN105426627A (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.)
Zhejiang Zhengya Dental Co ltd
Original Assignee
Zhejiang Zhengya Dental 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 Zhejiang Zhengya Dental Co ltd filed Critical Zhejiang Zhengya Dental Co ltd
Priority to CN201510907608.7A priority Critical patent/CN105426627B/zh
Publication of CN105426627A publication Critical patent/CN105426627A/zh
Application granted granted Critical
Publication of CN105426627B publication Critical patent/CN105426627B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • 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]

Abstract

本发明公开了一种机器人工装安装偏差计算系统,它按照以下步骤顺序进行:1)建立激光切割硬件设施;2)基础变换矩阵的建立;3)切割变换矩阵的建立;4)计算理想状态下多轴机器人角度;5)工装安装偏差的计算;采用对机器人工件的安装进行偏差计算,及时的将位置有偏差的工件进行矫正。该计算系统应用于激光切割或激光打标过程中对于工件位置偏差的矫正,能够提高最终产品的合格率。本发明适用于机器人工件安装时偏差的计算,尤其适用于机器人激光切割或激光打标时工件的矫正。

Description

机器人工装安装偏差计算系统及其应用
技术领域
本发明属于机械领域,涉及一种计算系统及其应用,具体涉及一种机器人工装安装偏差计算系统及其应用。
背景技术
目前,机器人运动学只涉及到物体的运动规律,不考虑产生运动的力和力矩。机器人正运动学所研究的内容是:给定机器人各关节的角度或唯一,求解计算机器人末端执行器相对于参考坐标系的位置和姿态问题。机械手是机器人上的重要部件,机械手是一系列由关节连接起来的机械臂构成的,每一个机械臂建立一个坐标系,并用齐次变换描述坐标系之间的相对位置和姿态。
当多个机械臂组成机械手时,就需要使用联体坐标,如联体坐标变换A、B、C,假设A是参考坐标系(基坐标系),则B相对于A的坐标变换以及C相对于B的坐标变换成为联体坐标变换。已知B在A中的表示为T1,C在B中的表示为T2,刚体在C中的表示为T3,则刚体在A中的表示为T=T1T2T3。但是当需要进行激光操作时,往往对于工件位置的精度具有较高的要求,激光切割部或激光打标部上的工件位置出现偏移时,就容易造成产品被切割坏或打标位置出现偏移等现象。
因此,研究一种能够及时矫正机器人工装安装偏差的计算系统,这对于切割或打标效率的提高具有重要的意义。
发明内容
本发明要解决的技术问题,是提供一种机器人工装安装偏差计算系统及其应用,采用对机器人工件的安装进行偏差计算,及时的将位置有偏差的工件进行矫正。该计算系统应用于激光切割或激光打标过程中对于工件位置偏差的矫正,能够提高最终产品的合格率。
为解决上述技术问题,本发明所采取的技术方案是:
一种机器人工装安装偏差计算系统,按照以下步骤顺序进行:
1)建立激光切割硬件设施
机器人激光切割机包括底座,串联在其上连续的N个机械臂,机械臂自由末端的抓取工装以及激光器;
2)基础变换矩阵的建立
建立激光坐标系(1),该坐标系原点为激光在产品上的入射点,激光坐标系在机器人坐标系(2)下的位置及姿态可以用一个4×4的旋转变换矩阵Qlaser2base表示:
Figure BDA0000871815530000021
其中:
Figure BDA0000871815530000022
为激光坐标系(1)在机器人坐标系(2)下的位置,其中,px,py,pz为该位 置在x,y,z轴方向的分量;
Figure BDA0000871815530000023
为激光坐标系(1)在机器人坐标系(2)下的姿态,其中,nx,ny,nz为该 姿态的X轴方向的x,y,z轴方向的分量;ox,oy,oz为该姿态的Y轴方向的x,y,z轴方向的分量; ax,ay,az为该姿态的Z轴方向的x,y,z轴方向的分量;
建立产品坐标系(3),其相对于机械臂自由末端的TCP坐标系(5)的位置及姿态可以用一个4×4的旋转变换矩阵Qplate2TCP表示,其中,TCP即为机器人机械臂自由末端的位置;
Figure BDA0000871815530000024
其中:
Figure BDA0000871815530000031
为产品坐标系(3)在TCP坐标系(5)下的位置,其中,px,py,pz为该位置在 x,y,z轴方向的分量;
Figure BDA0000871815530000032
为产品坐标系(3)在TCP坐标系(5)下的姿态;其中,nx,ny,nz为该姿 态的X轴方向的x,y,z轴方向的分量;ox,oy,oz为该姿态的x,y,z轴方向的分量;ax,ay,az为该 姿态的Z轴方向的x,y,z轴方向的分量;
3)切割变换矩阵的建立
获得待加工产品的三维数据,并从中提取出产品坐标系(3)下的激光入射点P1、P2、P3、……、Pn,以及对应的激光入射点法线方向N1、N2、N3、……、Nn;根据激光入射点和激光入射点法线方向建立局部坐标系(4),以该局部坐标系的原点为激光入射点,坐标系相对于产品坐标系(3)的位置及姿态可以用一个4×4的旋转变换矩阵Qmark2plate表示:
Figure BDA0000871815530000033
其中:
Figure BDA0000871815530000034
为每个入射点局部坐标系(4)在产品坐标系(3)下的位置,其中,px,py, pz为该位置在x,y,z轴方向的分量;
Figure BDA0000871815530000035
为入射点局部坐标系(4)在产品坐标系(3)下的位姿态,其中,nx,ny, nz为该姿态的X轴方向的x,y,z轴方向的分量;ox,oy,oz为该姿态的Y轴方向的x,y,z轴方向 的分量;ax,ay,az为该姿态的Z轴方向的x,y,z轴方向的分量,并且所述的Z轴方向即为入射 点的法线方向;
4)计算理想状态下多轴机器人角度
利用机器人运动学方程计算当激光垂直入射到“O”型标记时多轴机器人的角度,根据机器人的D-H参数,计算对应的机械臂可到达角度{qi,j},其中,qi,j表示产品的第i个激光入射点位置对应的机器人第j个机械臂的角度;
机器人运动方程为:
Figure BDA0000871815530000041
其中,j≥1;
其中
Figure BDA0000871815530000046
为机器第j各机械臂相对于第j-1个机械臂的旋转矩阵:
Figure BDA0000871815530000042
其中,θj,αj,dj和aj为机器人D-H参数,θj为机械臂之间的夹角,包含qi,j;αj为机械臂扭转角;dj为机械臂之间的距离;aj为机械臂长度;上述参数代入式Ⅱ中即可计算每个激光入射点对应的机器人机械臂角度值qi,j,i≥2;j≥1;
5)工装安装偏差的计算
工装坐标系与机器人末端TCP坐标系并不重合,假设工装坐标系在机器人末端TCP坐标系下的位置及姿态用一个4×4的旋转变换矩阵Qerr表示:
Figure BDA0000871815530000043
当θ为小角度时,经计算,
Figure BDA0000871815530000044
其中x0,y0为中心平移的位置,θ 为中心旋转的角度;
Figure BDA0000871815530000045
假设“O”型标记的位置为(pxi,pyi),然后发现当机器人多轴角度设置为计算值时,激光入射到“×”型标记,测量得到“×”型标记的位置为(pxi’,pyi’),经计算简化得下式:
Figure BDA0000871815530000051
已知(pxi,pyi)和(pxi’,pyi’),最终拟合得到x0,y0以及θ。
作为本发明的一种限定,所述的机械臂个数N≥1。
作为上述限定的进一步限定,所述的N=6。
本发明还提供了上述机器人工件安装偏差计算系统的应用,它可以应用于机器人工件偏差的矫正。
由于采用了上述的技术方案,本发明与现有技术相比,所取得的技术进步在于:
本发明采用对机器人工件的安装进行偏差计算,及时的将位置有偏差的工件进行矫正,如在进行激光打标过程或激光切割过程中对于工件的位置或者需要操作的工件位置精度要求非常高,如果有偏差,在实际的操作过程中就会出现不合格的产品,因此,能够及时的对工装的偏差进行矫正,具有重要的意义。
本发明适用于机器人工件安装时偏差的计算,尤其适用于机器人激光切割或激光打标时,工件的矫正。
本发明下面将结合说明书附图与具体实施例作进一步详细说明。
附图说明
图1为本发明实施例1的激光切割硬件设施的结构示意图;
图2为工件偏差示意图。
图中:11—底座,12—连杆,13—抓取工装,14—激光器;1—激光坐标系,2—机器人坐标系,3—产品坐标系,4—局部坐标系,5—TCP坐标系。
具体实施方式
实施例1一种机器人工装安装偏差计算系统及其应用
一种机器人工装安装偏差计算系统,按照以下步骤顺序进行:
1)建立激光切割硬件设施
机器人激光切割机包括底座11,串联在其上连续的6个机械臂12,机械臂自由末端的抓取工装13以及激光器14;
2)基础变换矩阵的建立
建立激光坐标系(1),该坐标系原点为激光在产品上的入射点,激光坐标系在机器人坐标系(2)下的位置及姿态可以用一个4×4的旋转变换矩阵Qlaser2base表示:
Figure BDA0000871815530000061
建立产品坐标系(3),其相对于机械臂自由末端的TCP 坐标系(5)的位置及姿态可以用一个4×4的旋转变换矩阵Qplate2TCP表示,其中,TCP即为机器 人机械臂自由末端的位置;
Figure BDA0000871815530000062
3)切割变换矩阵的建立
获得待加工产品的三维数据,并从中提取出产品坐标系(3)下的激光入射点P1、P2、P3、……、Pn,以及对应的激光入射点法线方向N1、N2、N3、……、Nn;根据激光入射点和激光入射点法线方向建立局部坐标系(4),以该局部坐标系的原点为激光入射点,坐标系相对于产品坐标系(3)的位置及姿态可以用一个4×4的旋转变换矩阵Qmark2plate表示:
对于切割点位置
Figure BDA0000871815530000063
而言:
Figure BDA0000871815530000071
4)计算理想状态下多轴机器人角度
利用机器人运动学方程计算当激光垂直入射到“O”型标记时多轴机器人的角度,根据机器人的D-H参数,计算对应的机械臂可到达角度{qi,j},其中,qi,j表示产品的第i个激光入射点位置对应的机器人第j个机械臂的角度;
机器人运动方程为:
Figure BDA0000871815530000072
其中,1≤j≤6;则
Figure BDA0000871815530000073
其中
Figure BDA0000871815530000074
为机器第j各机械臂相对于第j-1个机械臂的旋转矩阵:
Figure BDA0000871815530000075
其中,θj,αj,dj和aj为机器人D-H参数,θj为机械臂之间的夹角,包含qi,j;αj为机械臂扭转角;dj为机械臂之间的距离;aj为机械臂长度;上述参数代入式Ⅱ中即可计算每个激光切割点对应的机器人机械臂角度值qi,j,i≥2;j≥1;其中,当默认机器人每个轴运动范围在-180deg到180deg时,根据机器人实际轴运动范围从中选取合适解,本实施例选用的是ARB-120型机器人;
ARB-120型机器人D-H参数
Figure BDA0000871815530000076
Figure BDA0000871815530000086
将上述参数代入
Figure BDA0000871815530000081
的计算式中,求解出{qi,j};
5)工装安装偏差的计算
工装坐标系与机器人末端TCP坐标系并不重合,假设工装坐标系在机器人末端TCP坐标系下的位置及姿态用一个4×4的旋转变换矩阵Qerr表示:
Figure BDA0000871815530000082
当θ为小角度时,sinθ≈θ;cosθ≈1经计算得
Figure BDA0000871815530000083
其中x0,y0为 中心平移的位置,θ为中心旋转的角度;
Figure BDA0000871815530000084
假设“O”型标记的位置为(pxi,pyi),然后发现当机器人多轴角度设置为计算值时,激光入射到“×”型标记,测量得到“×”型标记的位置为(pxi’,pyi’),经计算简化得下式:
Figure BDA0000871815530000085
已知(pxi,pyi)和(pxi’,pyi’),最终拟合得到x0,y0以及θ。
上述机器人工件安装偏差计算系统应用于机器人工件偏差的矫正,在进行切割或激光打标时能够更精准的定位工件的位置。
以上所述,仅是本发明的较佳实施例而已,并非是对本发明作其它形式的限定,任何熟悉本专业的技术人员可能利用上述技术内容作为启示加以变更或改型为等同变化的等效实施例。但是凡是未脱离本发明技术方案内容,依据本发明的技术实质对以上实施例所作出的简单修改,等同变化与改型,仍属于本发明权利要求的保护范围。

Claims (4)

1.一种机器人工装安装偏差计算系统,其特征在于它按照以下步骤顺序进行:
1)建立激光切割硬件设施
机器人激光切割机包括底座(11),串联在其上连续的N个机械臂(12),机械臂自由末端的抓取工装(13)以及激光器(14);
2)基础变换矩阵的建立
建立激光坐标系(1),该坐标系原点为激光在产品上的入射点,激光坐标系在机器人坐标系(2)下的位置及姿态可以用一个4×4的旋转变换矩阵Qlaser2base表示:
Figure QLYQS_1
其中:
Figure QLYQS_2
为激光坐标系(1)在机器人坐标系(2)下的位置,其中,px,py,pz为该位置在x,y,z轴方向的分量;
Figure QLYQS_3
为激光坐标系(1)在机器人坐标系(2)下的姿态,其中,nx,ny,nz为该姿态的X轴方向的x,y,z轴方向的分量;ox,oy,oz为该姿态的Y轴方向的x,y,z轴方向的分量;ax,ay,az为该姿态的Z轴方向的x,y,z轴方向的分量;
建立产品坐标系(3),其相对于机械臂自由末端的TCP坐标系(5)的位置及姿态可以用一个4×4的旋转变换矩阵Qplate2TCP表示,其中,TCP即为机器人机械臂自由末端的位置;
Figure QLYQS_4
其中:
Figure QLYQS_5
为产品坐标系(3)在TCP坐标系(5)下的位置,其中,px,py,pz为该位置在x,y,z轴方向的分量;
Figure QLYQS_6
为产品坐标系(3)在TCP坐标系(5)下的姿态;其中,nx,ny,nz为该姿态的X轴方向的x,y,z轴方向的分量;ox,oy,oz为该姿态的x,y,z轴方向的分量;ax,ay,az为该姿态的Z轴方向的x,y,z轴方向的分量;
3)切割变换矩阵的建立
获得待加工产品的三维数据,并从中提取出产品坐标系(3)下的激光入射点P1、P2、P3、……、Pn,以及对应的激光入射点法线方向N1、N2、N3、……、Nn;根据激光入射点和激光入射点法线方向建立局部坐标系(4),以该局部坐标系的原点为激光入射点,坐标系相对于产品坐标系(3)的位置及姿态可以用一个4×4的旋转变换矩阵Qmark2plate表示:
Figure QLYQS_7
其中:
Figure QLYQS_8
为每个入射点局部坐标系(4)在产品坐标系(3)下的位置,其中,px,py,pz为该位置在x,y,z轴方向的分量;
Figure QLYQS_9
为入射点局部坐标系(4)在产品坐标系(3)下的位姿态,其中,nx,ny,nz为该姿态的X轴方向的x,y,z轴方向的分量;ox,oy,oz为该姿态的Y轴方向的x,y,z轴方向的分量;ax,ay,az为该姿态的Z轴方向的x,y,z轴方向的分量,并且所述的Z轴方向即为入射点的法线方向;
4)计算理想状态下多轴机器人角度
利用机器人运动学方程计算当激光垂直入射到“O”型标记时多轴机器人的角度,根据机器人的D-H参数,计算对应的机械臂可到达角度{qi,j},其中,qi,j表示产品的第i个激光入射点位置对应的机器人第j个机械臂的角度;
机器人运动方程为:
Figure QLYQS_10
其中,j≥1;
其中
Figure QLYQS_11
为机器第j个机械臂相对于第j-1个机械臂的旋转矩阵:
Figure QLYQS_12
其中,θj,αj,dj和aj为机器人D-H参数,θj为机械臂之间的夹角,包含qi,j;αj为机械臂扭转角;dj为机械臂之间的距离;aj为机械臂长度;上述参数代入
Figure QLYQS_13
中即可计算每个激光入射点对应的机器人机械臂角度值qi,j,i≥2;j≥1;5)工装安装偏差的计算
工装坐标系与机器人末端TCP坐标系并不重合,假设工装坐标系在机器人末端TCP坐标系下的位置及姿态用一个4×4的旋转变换矩阵Qerr表示:
Figure QLYQS_14
当θ为小角度时,经计算,
Figure QLYQS_15
其中x0,y0为中心平移的位置,θ为中心旋转的角度;
Figure QLYQS_16
假设“O”型标记的位置为(pxi,pyi),然后发现当机器人多轴角度设置为计算值时,激光入射到“×”型标记,测量得到“×”型标记的位置为(pxi , ,pyi , ),经计算简化得下式:
Figure QLYQS_17
已知(pxi,pyi)和(pxi , ,pyi , ),最终拟合得到x0,y0以及θ。
2.根据权利要求1所述的机器人工装安装偏差计算系统,其特征在于:所述的机械臂个数N≥1。
3.根据权利要求2所述的机器人工装安装偏差计算系统,其特征在于:所述的N=6。
4.一种如权利要求1-3中任一项所述的机器人工装安装偏差计算系统的应用,
其特征在于:所述的计算系统应用于机器人工件偏差的矫正。
CN201510907608.7A 2015-12-09 2015-12-09 机器人工装安装偏差计算系统及其应用 Active CN105426627B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510907608.7A CN105426627B (zh) 2015-12-09 2015-12-09 机器人工装安装偏差计算系统及其应用

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510907608.7A CN105426627B (zh) 2015-12-09 2015-12-09 机器人工装安装偏差计算系统及其应用

Publications (2)

Publication Number Publication Date
CN105426627A CN105426627A (zh) 2016-03-23
CN105426627B true CN105426627B (zh) 2023-06-20

Family

ID=55504836

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510907608.7A Active CN105426627B (zh) 2015-12-09 2015-12-09 机器人工装安装偏差计算系统及其应用

Country Status (1)

Country Link
CN (1) CN105426627B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101968341A (zh) * 2010-08-31 2011-02-09 南京理工大学 一种工业机器人零位自标定方法及装置
CN102087096A (zh) * 2010-11-12 2011-06-08 浙江大学 一种基于激光跟踪测量的机器人工具坐标系自动标定装置及方法
CN103322953A (zh) * 2013-05-22 2013-09-25 北京配天大富精密机械有限公司 工件坐标系的标定方法、装置及工件加工处理方法、装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101968341A (zh) * 2010-08-31 2011-02-09 南京理工大学 一种工业机器人零位自标定方法及装置
CN102087096A (zh) * 2010-11-12 2011-06-08 浙江大学 一种基于激光跟踪测量的机器人工具坐标系自动标定装置及方法
CN103322953A (zh) * 2013-05-22 2013-09-25 北京配天大富精密机械有限公司 工件坐标系的标定方法、装置及工件加工处理方法、装置

Also Published As

Publication number Publication date
CN105426627A (zh) 2016-03-23

Similar Documents

Publication Publication Date Title
US20190099887A1 (en) System and method for calibrating tool center point of robot
Tian et al. Auto-normalization algorithm for robotic precision drilling system in aircraft component assembly
CN110757454B (zh) 一种双机器人协同旋转的路径规划方法和装置
CN113001535B (zh) 机器人工件坐标系自动校正系统与方法
CN103942427A (zh) 一类六自由度机械臂运动学逆解的快速简便求法
CN112123341B (zh) 机器人双臂协调运动控制方法、装置和电子设备
CN104385283B (zh) 一种六自由度机械臂奇异位形的快速判断方法
Guo et al. Vision based navigation for Omni-directional mobile industrial robot
CN107791248B (zh) 基于不满足Pieper准则的六自由度串联机器人的控制方法
CN106844951B (zh) 基于分段几何法求解超冗余机器人逆运动学的方法及系统
Lee et al. High precision hand-eye self-calibration for industrial robots
Liu et al. Laser tracker-based control for peg-in-hole assembly robot
CN112476435B (zh) 重力加速度方向的标定方法、标定装置及存储介质
JP5505155B2 (ja) ロボットシステムおよびロボット制御方法
CN110026986B (zh) 多臂协同作业柔顺控制方法
Mustafa et al. Kinematic calibration of a 7-DOF self-calibrated modular cable-driven robotic arm
CN105426627B (zh) 机器人工装安装偏差计算系统及其应用
Huang et al. Robotic hand-eye calibration based on active vision
CN112847362B (zh) 一种基于图像矩的搬运机器人视觉伺服控制方法
Si et al. A complete solution to the inverse kinematics problem for 4-DOF manipulator robot
JP2021186929A (ja) 多軸ロボットの制御方法
Liu et al. Automatic calibration algorithm of robot TCP based on binocular vision
Ju et al. Inverse kinematics of 2-dof space manipulator based on dh parameters
Bian et al. A path planning method of robotic belt grinding system for grinding workpieces with complex shape surfaces
Deng et al. Calibration of a robot vision system coupled with structured light: method and experiments

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information

Address after: 314006 Building 2, 3556 Linggongtang Road, Daqiao Town, Nanhu District, Jiaxing City, Zhejiang Province

Applicant after: Zhejiang Zhengya Dental Co.,Ltd.

Address before: East of 5th Floor, No. 239 Asia Pacific Road, Nanhu District, Jiaxing City, Zhejiang Province, 314006

Applicant before: ZHEJIANG YINCHILI MEDICAL TECHNOLOGY Co.,Ltd.

CB02 Change of applicant information
GR01 Patent grant
GR01 Patent grant