CN105426627B - 机器人工装安装偏差计算***及其应用 - Google Patents
机器人工装安装偏差计算***及其应用 Download PDFInfo
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/02—Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]
Landscapes
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Numerical Control (AREA)
- Manipulator (AREA)
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表示:
为激光坐标系(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即为机器人机械臂自由末端的位置;
为产品坐标系(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表示:
为入射点局部坐标系(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个机械臂的角度;
机器人运动方程为:
其中,j≥1;
其中,θj,αj,dj和aj为机器人D-H参数,θj为机械臂之间的夹角,包含qi,j;αj为机械臂扭转角;dj为机械臂之间的距离;aj为机械臂长度;上述参数代入式Ⅱ中即可计算每个激光入射点对应的机器人机械臂角度值qi,j,i≥2;j≥1;
5)工装安装偏差的计算
工装坐标系与机器人末端TCP坐标系并不重合,假设工装坐标系在机器人末端TCP坐标系下的位置及姿态用一个4×4的旋转变换矩阵Qerr表示:
假设“O”型标记的位置为(pxi,pyi),然后发现当机器人多轴角度设置为计算值时,激光入射到“×”型标记,测量得到“×”型标记的位置为(pxi’,pyi’),经计算简化得下式:
已知(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表示:
3)切割变换矩阵的建立
获得待加工产品的三维数据,并从中提取出产品坐标系(3)下的激光入射点P1、P2、P3、……、Pn,以及对应的激光入射点法线方向N1、N2、N3、……、Nn;根据激光入射点和激光入射点法线方向建立局部坐标系(4),以该局部坐标系的原点为激光入射点,坐标系相对于产品坐标系(3)的位置及姿态可以用一个4×4的旋转变换矩阵Qmark2plate表示:
4)计算理想状态下多轴机器人角度
利用机器人运动学方程计算当激光垂直入射到“O”型标记时多轴机器人的角度,根据机器人的D-H参数,计算对应的机械臂可到达角度{qi,j},其中,qi,j表示产品的第i个激光入射点位置对应的机器人第j个机械臂的角度;
机器人运动方程为:
其中,θ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参数
5)工装安装偏差的计算
工装坐标系与机器人末端TCP坐标系并不重合,假设工装坐标系在机器人末端TCP坐标系下的位置及姿态用一个4×4的旋转变换矩阵Qerr表示:
假设“O”型标记的位置为(pxi,pyi),然后发现当机器人多轴角度设置为计算值时,激光入射到“×”型标记,测量得到“×”型标记的位置为(pxi’,pyi’),经计算简化得下式:
已知(pxi,pyi)和(pxi’,pyi’),最终拟合得到x0,y0以及θ。
上述机器人工件安装偏差计算***应用于机器人工件偏差的矫正,在进行切割或激光打标时能够更精准的定位工件的位置。
以上所述,仅是本发明的较佳实施例而已,并非是对本发明作其它形式的限定,任何熟悉本专业的技术人员可能利用上述技术内容作为启示加以变更或改型为等同变化的等效实施例。但是凡是未脱离本发明技术方案内容,依据本发明的技术实质对以上实施例所作出的简单修改,等同变化与改型,仍属于本发明权利要求的保护范围。
Claims (4)
1.一种机器人工装安装偏差计算***,其特征在于它按照以下步骤顺序进行:
1)建立激光切割硬件设施
机器人激光切割机包括底座(11),串联在其上连续的N个机械臂(12),机械臂自由末端的抓取工装(13)以及激光器(14);
2)基础变换矩阵的建立
建立激光坐标系(1),该坐标系原点为激光在产品上的入射点,激光坐标系在机器人坐标系(2)下的位置及姿态可以用一个4×4的旋转变换矩阵Qlaser2base表示:
其中:
为激光坐标系(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即为机器人机械臂自由末端的位置;
其中:
为产品坐标系(3)在TCP坐标系(5)下的位置,其中,px,py,pz为该位置在x,y,z轴方向的分量;
为产品坐标系(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表示:
其中:
为每个入射点局部坐标系(4)在产品坐标系(3)下的位置,其中,px,py,pz为该位置在x,y,z轴方向的分量;
为入射点局部坐标系(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个机械臂的角度;
机器人运动方程为:
其中,j≥1;
其中,θj,αj,dj和aj为机器人D-H参数,θj为机械臂之间的夹角,包含qi,j;αj为机械臂扭转角;dj为机械臂之间的距离;aj为机械臂长度;上述参数代入中即可计算每个激光入射点对应的机器人机械臂角度值qi,j,i≥2;j≥1;5)工装安装偏差的计算
工装坐标系与机器人末端TCP坐标系并不重合,假设工装坐标系在机器人末端TCP坐标系下的位置及姿态用一个4×4的旋转变换矩阵Qerr表示:
当θ为小角度时,经计算,
其中x0,y0为中心平移的位置,θ为中心旋转的角度;
假设“O”型标记的位置为(pxi,pyi),然后发现当机器人多轴角度设置为计算值时,激光入射到“×”型标记,测量得到“×”型标记的位置为(pxi , ,pyi , ),经计算简化得下式:
已知(pxi,pyi)和(pxi , ,pyi , ),最终拟合得到x0,y0以及θ。
2.根据权利要求1所述的机器人工装安装偏差计算***,其特征在于:所述的机械臂个数N≥1。
3.根据权利要求2所述的机器人工装安装偏差计算***,其特征在于:所述的N=6。
4.一种如权利要求1-3中任一项所述的机器人工装安装偏差计算***的应用,
其特征在于:所述的计算***应用于机器人工件偏差的矫正。
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)
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 | 北京配天大富精密机械有限公司 | 工件坐标系的标定方法、装置及工件加工处理方法、装置 |
-
2015
- 2015-12-09 CN CN201510907608.7A patent/CN105426627B/zh active Active
Patent Citations (3)
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 |
---|---|---|
CN106166749B (zh) | 一种空间移动多臂机器人的移动轨迹规划方法 | |
US20190099887A1 (en) | System and method for calibrating tool center point of robot | |
CN110757454B (zh) | 一种双机器人协同旋转的路径规划方法和装置 | |
CN103942427B (zh) | 一类六自由度机械臂运动学逆解的快速简便求法 | |
Tian et al. | Auto-normalization algorithm for robotic precision drilling system in aircraft component assembly | |
CN113001535B (zh) | 机器人工件坐标系自动校正***与方法 | |
CN104408299B (zh) | 基于距离识别冗余运动学参数的机器人位置误差补偿方法 | |
CN106864776B (zh) | 一种基于对接环的捕获目标卫星的方法与*** | |
CN104385283B (zh) | 一种六自由度机械臂奇异位形的快速判断方法 | |
CN105856231B (zh) | 一种特定构型六轴工业机器人的运动控制方法 | |
CN105938512B (zh) | 机器人激光切割路径的计算方法及其应用 | |
CN106844951B (zh) | 基于分段几何法求解超冗余机器人逆运动学的方法及*** | |
CN107791248B (zh) | 基于不满足Pieper准则的六自由度串联机器人的控制方法 | |
CN112123341A (zh) | 机器人双臂协调运动控制方法、装置和电子设备 | |
CN110480638A (zh) | 一种多关节机器人位姿误差自补偿的码垛方法及其码垛*** | |
CN111844005A (zh) | 一种应用于隧道湿喷的2r-p-2r-p-2r机械臂运动规划方法 | |
Liu et al. | Laser tracker-based control for peg-in-hole assembly robot | |
JP5505155B2 (ja) | ロボットシステムおよびロボット制御方法 | |
CN110026986B (zh) | 多臂协同作业柔顺控制方法 | |
CN105426627B (zh) | 机器人工装安装偏差计算***及其应用 | |
CN112476435A (zh) | 重力加速度方向的标定方法、标定装置及存储介质 | |
Huang et al. | Robotic hand-eye calibration based on active vision | |
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 |
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 |