CN111409106A - 工业机器人距离性能测试方法 - Google Patents

工业机器人距离性能测试方法 Download PDF

Info

Publication number
CN111409106A
CN111409106A CN202010236384.2A CN202010236384A CN111409106A CN 111409106 A CN111409106 A CN 111409106A CN 202010236384 A CN202010236384 A CN 202010236384A CN 111409106 A CN111409106 A CN 111409106A
Authority
CN
China
Prior art keywords
distance
robot
coordinate
coordinate system
points
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.)
Pending
Application number
CN202010236384.2A
Other languages
English (en)
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.)
Borunte Robot Co Ltd
Original Assignee
Borunte Robot 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 Borunte Robot Co Ltd filed Critical Borunte Robot Co Ltd
Priority to CN202010236384.2A priority Critical patent/CN111409106A/zh
Publication of CN111409106A publication Critical patent/CN111409106A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B25HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
    • B25JMANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
    • B25J19/00Accessories fitted to manipulators, e.g. for monitoring, for viewing; Safety devices combined with or specially adapted for use in connection with manipulators
    • B25J19/0095Means or methods for testing manipulators

Abstract

本发明涉及机器人的距离性能参数技术领域,尤其为一种工业机器人距离性能测试方法,包括:选取测试点位:根据机器人的工作空间,选取实验轨迹上P2与P4两个点,采用连续轨迹编程,使机器人控制器驱动机器人末端从P4开始,在P2、P4位姿间连续运动,循环30次,单方向进行测量,记录跟踪仪采集数据和指令位姿数据,用两点间的距离公式计算指令距离,用两点间的距离公式计算30次循环的实到距离之和后,除以循环次数,得到实际距离的平均值,通过位置距离准确度公式计算位置距离准确度,通过位置距离重复性公式计算位置距离重复性。本发明对计算原理进行了详细的推导说明,基于该方法可以开发相应的程序实现机器人的距离性能指标测试。

Description

工业机器人距离性能测试方法
技术领域
本发明涉及机器人的距离性能参数技术领域,具体为一种工业机器人距离性能测试方法。
背景技术
工业机器人因同时具备通用性、高柔性、高精度等诸多特点,在全球制造业智能化趋势的推动下蓬勃发展。目前,工业机器人大多采用开环控制,为保证机器人的末端精度,在出厂前或使用一段时间后需要对机器人性能起显著影响的特性进行测试与调试。
机器人的距离性能参数是其中之一特性。机器人按程序设定的轨迹从起始指令点运动到终止指时,便会出现距离偏差。机器人距离性能功能测试参数分为距离准确度(AD)测试和距离重复性测试(RD)。专门用于机器人的测量系统价使用和维护成本高昂,且测量的机理均不清楚。
发明内容
本发明的目的在于提供一种工业机器人距离性能测试方法,以解决上述背景技术中提出的问题。
为实现上述目的,本发明提供如下技术方案:
一种工业机器人距离性能测试方法,包括如下步骤:
步骤1、选取测试点位:根据机器人的工作空间,选取实验轨迹上P2与P4两个点;
步骤2、采用连续轨迹编程,使机器人控制器驱动机器人末端从P4开始,在P2、P4位姿间连续运动,记录跟踪仪采集数据和指令位姿数据;
步骤3、用两点间的距离公式:
Figure BDA0002431124770000011
计算指令距离Di
步骤4、用两点间的距离公式:
Figure BDA0002431124770000021
得到实际距离的平均值Da
步骤5、通过位置距离准确度公式:
Figure BDA0002431124770000022
式中:
Figure BDA0002431124770000023
Figure BDA0002431124770000024
Figure BDA0002431124770000025
其中:Xc1、Yc1、Zc1是在机器人控制中可用的Pc1坐标;
Xc2、Yc2、Zc2是在机器人控制中可用的Pc2坐标;
X1j、Y1j、Z1j是P1j的坐标;
X2j、Y2j、Z2j是P2j的坐标;
n是重复次数;
计算位置距离准确度:ADp=Da-Di
步骤6、通过位置距离重复性公式:
Figure BDA0002431124770000026
计算位置距离重复性。
进一步的,步骤2中,连续运动为,循环30次,单方向进行测量。
进一步的,步骤4中,计算30次循环的实到距离之和后,除以循环次数,得到实际距离的平均值Da
进一步的,选取测试点位需要先完成三维空间测量仪坐标系与机器人坐标系的坐标转换,包括:首先,三维空间测量仪采集机器人末端法兰盘的球极坐标点位数据方位角
Figure BDA0002431124770000031
仰角θ和距离γ数据,然后通过球极坐标系与直角坐标系转换,可以获取测量仪的直角坐标点位数据,最后,再完成三维空间测量仪坐标系与机器人坐标系的坐标转换。
进一步的,球极坐标系(γ、θ、
Figure BDA0002431124770000032
)与直角坐标系(x、y、z)转换公式如下:
Figure BDA0002431124770000033
Figure BDA0002431124770000034
z=γ·cosθ
三维空间测量仪坐标系与机器人坐标系坐标转换公式如下:
Pt=RPr+T
R——旋转矩阵,T——平移矩阵;
Pt——三维空间测量仪坐标系下的坐标点位;
Pr——机器人坐标系下的坐标点位;
任意点Pi坐标的矩阵表示:
三维空间测量仪坐标系下的坐标点位:
Figure BDA0002431124770000035
机器人坐标系下的坐标点位:
Figure BDA0002431124770000036
进一步的,采用SVD法算出R、T,包括:
(1)机器人坐标系下指令坐标点位的点集为:
Figure BDA0002431124770000037
指令坐标点位在测量仪下测得的实际坐标点位构成的点集:
Figure BDA0002431124770000038
(2)分别计算机器人和测量仪坐标系下的坐标点位点集Pr、Pt的重心,即坐标点位点集包含的全部点的坐标的平均值,分别为:
Figure BDA0002431124770000039
Figure BDA0002431124770000041
(3)将两个坐标点位点集的重心对齐重合,分别计算各点集相对重心的相对坐标构成新的点集:
Figure BDA0002431124770000042
i=1,2,3···n,记:
Figure BDA0002431124770000043
为:
Figure BDA0002431124770000044
Figure BDA0002431124770000045
i=1,2,3···n,记:
Figure BDA0002431124770000046
为:
Figure BDA0002431124770000047
(4)由坐标点位点集
Figure BDA0002431124770000048
构造协矩阵
Figure BDA0002431124770000049
对协矩阵进行SVD分解:
Figure BDA00024311247700000410
(5)旋转矩阵即为R3×3=VUT,平移矩阵为T3×1=μr-Rμt,当n≥3时既可求出R矩阵,R的各列为长度为3的单位向量,且两两相互垂直。
与现有技术相比,本发明的有益效果是:
本发明,所述的是测量机器人距离性能指标的方法,文中对计算原理进行了详细的推导说明,基于该方法可以开发相应的程序实现机器人的距离性能指标测试。机器人距离功能测试对机器人按照指令轨迹从P1运行到P2后的两组指令位姿与两组实到位姿均值之间的距离偏差和在这两个位姿间一系列重复移动的距离波动进行测试。姿态距离准确度计算相当于单轴距离准确度,姿态距离重复性相当于单轴距离重复性。
附图说明
图1为本发明测量点的选择示意图。
图2为本发明运动轨迹示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
在本发明的描述中,需要说明的是,术语“上/下端”、“内”、“外”“前端”、“后端”、“两端”、“一端”、“另一端”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性。
在本发明的描述中,需要说明的是,除非另有明确的规定和限定,术语“安装”、“设置/套设有”、“套接”、“连接”等,应做广义理解,例如“连接”,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本发明中的具体含义。
请参阅图1-2,本发明提供一种技术方案:
一、坐标转换
1.建立两个坐标系对应关系
工业机器人要进行位姿特性测量的话,需要先完成三维空间测量仪坐标系与机器人坐标系的坐标转换。首先,三维空间测量仪采集机器人末端法兰盘的球极坐标点位数据——方位角
Figure BDA0002431124770000051
仰角θ和距离γ数据。然后通过球极坐标系与直角坐标系转换,可以获取测量仪的直角坐标点位数据。最后,再完成三维空间测量仪坐标系与机器人坐标系的坐标转换。
球极坐标系(γ、θ、
Figure BDA0002431124770000052
)与直角坐标系(x、y、z)转换公式如下:
Figure BDA0002431124770000053
Figure BDA0002431124770000061
z=γ·cosθ
三维空间测量仪坐标系与机器人坐标系坐标转换公式如下:
Pt=RPr+T
R——旋转矩阵,T——平移矩阵;
Pt——三维空间测量仪坐标系下的坐标点位;
Pr——机器人坐标系下的坐标点位。
任意点Pi坐标的矩阵表示:
三维空间测量仪坐标系下的坐标点位:
Figure BDA0002431124770000062
机器人坐标系下的坐标点位:
Figure BDA0002431124770000063
2.SVD法(奇异值分解法)算出R、T
(1)机器人坐标系下指令坐标点位的点集为:
Figure BDA0002431124770000064
指令坐标点位在测量仪下测得的实际坐标点位构成的点集:
Figure BDA0002431124770000065
(2)分别计算机器人和测量仪坐标系下的坐标点位点集Pr、Pt的重心,即坐标点位点集包含的全部点的坐标的平均值,分别为:
Figure BDA0002431124770000066
Figure BDA0002431124770000067
(3)将两个坐标点位点集的重心对齐重合,分别计算各点集相对重心的相对坐标构成新的点集:
Figure BDA0002431124770000068
i=1,2,3···n,记:
Figure BDA0002431124770000069
为:
Figure BDA00024311247700000610
Figure BDA00024311247700000611
i=1,2,3···n,记:
Figure BDA00024311247700000612
为:
Figure BDA00024311247700000613
(4)由坐标点位点集
Figure BDA00024311247700000614
构造协矩阵
Figure BDA00024311247700000615
对协矩阵进行SVD分解:
Figure BDA0002431124770000071
(5)旋转矩阵即为R3×3=VUT,平移矩阵为T3×1=μr-Rμt,当n≥3时既可求出R矩阵,R的各列为长度为3的单位向量,且两两相互垂直。为保证拟合的优度,n取5。也可根据需要适当的增加点数。
二、测量方法
1、选取测试点位:根据机器人的工作空间,选取实验轨迹上P2与P4两个点,点位位置如附图1所示。
2、采用连续轨迹编程,使机器人控制器驱动机器人末端从P4开始,在P2、P4位姿间连续运动,循环30次,单方向进行测量,如附图2所示,记录跟踪仪采集数据和指令位姿数据。
3、用两点间的距离公式:
Figure BDA0002431124770000072
计算指令距离Di
4、用两点间的距离公式:
Figure BDA0002431124770000073
计算30次循环的实到距离之和后,除以循环次数,得到实际距离的平均值Da
5、通过位置距离准确度公式:
Figure BDA0002431124770000074
式中:
Figure BDA0002431124770000075
Figure BDA0002431124770000076
Figure BDA0002431124770000077
其中:Xc1、Yc1、Zc1是在机器人控制中可用的Pc1坐标;
Xc2、Yc2、Zc2是在机器人控制中可用的Pc2坐标;
X1j、Y1j、Z1j是P1j的坐标;
X2j、Y2j、Z2j是P2j的坐标;
n是重复次数。
计算位置距离准确度:ADp=Da-Di
6、通过位置距离重复性公式:
Figure BDA0002431124770000081
计算位置距离重复性。
本发明,图1、图2中C1-C8为机器人最大空间的立方体,P2、P4为测量点。
本发明,所述的是测量机器人距离性能指标的方法,文中对计算原理进行了详细的推导说明,基于该方法可以开发相应的程序实现机器人的距离性能指标测试。机器人距离功能测试对机器人按照指令轨迹从P1运行到P2后的两组指令位姿与两组实到位姿均值之间的距离偏差和在这两个位姿间一系列重复移动的距离波动进行测试。姿态距离准确度计算相当于单轴距离准确度,姿态距离重复性相当于单轴距离重复性。
尽管已经示出和描述了本发明的实施例,对于本领域的普通技术人员而言,可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由所附权利要求及其等同物限定。

Claims (6)

1.一种工业机器人距离性能测试方法,其特征在于,包括如下步骤:
步骤1、选取测试点位:根据机器人的工作空间,选取实验轨迹上P2与P4两个点;
步骤2、采用连续轨迹编程,使机器人控制器驱动机器人末端从P4开始,在P2、P4位姿间连续运动,记录跟踪仪采集数据和指令位姿数据;
步骤3、用两点间的距离公式:
Figure FDA0002431124760000011
计算指令距离Di
步骤4、用两点间的距离公式:
Figure FDA0002431124760000012
得到实际距离的平均值Da
步骤5、通过位置距离准确度公式:
Figure FDA0002431124760000013
式中:
Figure FDA0002431124760000014
Figure FDA0002431124760000015
Figure FDA0002431124760000016
其中:Xc1、Yc1、Zc1是在机器人控制中可用的Pc1坐标;
Xc2、Yc2、Zc2是在机器人控制中可用的Pc2坐标;
X1j、Y1j、Z1j是P1j的坐标;
X2j、Y2j、Z2j是P2j的坐标;
n是重复次数;
计算位置距离准确度:ADp=Da-Di
步骤6、通过位置距离重复性公式:
Figure FDA0002431124760000021
计算位置距离重复性。
2.根据权利要求1所述的一种工业机器人距离性能测试方法,其特征在于,步骤2中,连续运动为,循环30次,单方向进行测量。
3.根据权利要求1所述的一种工业机器人距离性能测试方法,其特征在于,步骤4中,计算30次循环的实到距离之和后,除以循环次数,得到实际距离的平均值Da
4.根据权利要求1所述的一种工业机器人距离性能测试方法,其特征在于,选取测试点位需要先完成三维空间测量仪坐标系与机器人坐标系的坐标转换,包括:首先,三维空间测量仪采集机器人末端法兰盘的球极坐标点位数据方位角
Figure FDA0002431124760000022
仰角θ和距离γ数据,然后通过球极坐标系与直角坐标系转换,可以获取测量仪的直角坐标点位数据,最后,再完成三维空间测量仪坐标系与机器人坐标系的坐标转换。
5.根据权利要求4所述的一种工业机器人距离性能测试方法,其特征在于,球极坐标系
Figure FDA0002431124760000023
与直角坐标系(x、y、z)转换公式如下:
Figure FDA0002431124760000024
Figure FDA0002431124760000025
z=γ·cosθ
三维空间测量仪坐标系与机器人坐标系坐标转换公式如下:
Pt=RPr+T
R——旋转矩阵,T——平移矩阵;
Pt——三维空间测量仪坐标系下的坐标点位;
Pr——机器人坐标系下的坐标点位;
任意点Pi坐标的矩阵表示:
三维空间测量仪坐标系下的坐标点位:
Figure FDA0002431124760000031
机器人坐标系下的坐标点位:
Figure FDA0002431124760000032
6.根据权利要求5所述的一种工业机器人距离性能测试方法,其特征在于,采用SVD法算出R、T,包括:
(1)机器人坐标系下指令坐标点位的点集为:
Figure FDA0002431124760000033
指令坐标点位在测量仪下测得的实际坐标点位构成的点集:
Figure FDA0002431124760000034
(2)分别计算机器人和测量仪坐标系下的坐标点位点集Pr、Pt的重心,即坐标点位点集包含的全部点的坐标的平均值,分别为:
Figure FDA0002431124760000035
Figure FDA0002431124760000036
(3)将两个坐标点位点集的重心对齐重合,分别计算各点集相对重心的相对坐标构成新的点集:
Figure FDA0002431124760000037
记:
Figure FDA0002431124760000038
为:
Figure FDA0002431124760000039
Figure FDA00024311247600000310
记:
Figure FDA00024311247600000311
为:
Figure FDA00024311247600000312
(4)由坐标点位点集
Figure FDA00024311247600000313
构造协矩阵
Figure FDA00024311247600000315
对协矩阵进行SVD分解:
Figure FDA00024311247600000316
(5)旋转矩阵即为R3×3=VUT,平移矩阵为T3×1=μr-Rμt,当n≥3时既可求出R矩阵,R的各列为长度为3的单位向量,且两两相互垂直。
CN202010236384.2A 2020-03-30 2020-03-30 工业机器人距离性能测试方法 Pending CN111409106A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010236384.2A CN111409106A (zh) 2020-03-30 2020-03-30 工业机器人距离性能测试方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010236384.2A CN111409106A (zh) 2020-03-30 2020-03-30 工业机器人距离性能测试方法

Publications (1)

Publication Number Publication Date
CN111409106A true CN111409106A (zh) 2020-07-14

Family

ID=71487884

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010236384.2A Pending CN111409106A (zh) 2020-03-30 2020-03-30 工业机器人距离性能测试方法

Country Status (1)

Country Link
CN (1) CN111409106A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112549019A (zh) * 2020-11-06 2021-03-26 北京工业大学 一种基于连续动态时间规整的工业机器人轨迹准确度分析方法
CN112797931A (zh) * 2020-12-28 2021-05-14 韩莉莉 工业机器人位姿准确度和位姿重复性检测装置及检测方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006053754A (ja) * 2004-08-11 2006-02-23 Honda Motor Co Ltd 平面検出装置及び検出方法
CN106289379A (zh) * 2016-08-12 2017-01-04 杭州亿恒科技有限公司 工业机器人性能测量方法
JP2018009917A (ja) * 2016-07-15 2018-01-18 株式会社リコー 測距装置、移動体、ロボット、3次元計測装置及び測距方法
CN107972070A (zh) * 2017-11-29 2018-05-01 上海新时达机器人有限公司 机器人性能的测试方法、测试系统及计算机可读存储介质
CN110161850A (zh) * 2019-04-24 2019-08-23 南京航空航天大学 一种工业机器人变参数刚度辨识与建模方法
CN110370314A (zh) * 2019-06-11 2019-10-25 杭州亿恒科技有限公司 基于多激光跟踪仪的双臂机器人性能测量方法
CN110497373A (zh) * 2019-08-07 2019-11-26 大连理工大学 一种移动作业机器人的三维激光雷达与机械臂间的联合标定方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006053754A (ja) * 2004-08-11 2006-02-23 Honda Motor Co Ltd 平面検出装置及び検出方法
JP2018009917A (ja) * 2016-07-15 2018-01-18 株式会社リコー 測距装置、移動体、ロボット、3次元計測装置及び測距方法
CN106289379A (zh) * 2016-08-12 2017-01-04 杭州亿恒科技有限公司 工业机器人性能测量方法
CN107972070A (zh) * 2017-11-29 2018-05-01 上海新时达机器人有限公司 机器人性能的测试方法、测试系统及计算机可读存储介质
CN110161850A (zh) * 2019-04-24 2019-08-23 南京航空航天大学 一种工业机器人变参数刚度辨识与建模方法
CN110370314A (zh) * 2019-06-11 2019-10-25 杭州亿恒科技有限公司 基于多激光跟踪仪的双臂机器人性能测量方法
CN110497373A (zh) * 2019-08-07 2019-11-26 大连理工大学 一种移动作业机器人的三维激光雷达与机械臂间的联合标定方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
杨书评: "《工业机器人性能规范及其试验方法(GB/T 12642-2013)》", 12 November 2013 *
谢锡麟: "《复旦先导讲义 -数理基础与程序设计》", 30 June 2018 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112549019A (zh) * 2020-11-06 2021-03-26 北京工业大学 一种基于连续动态时间规整的工业机器人轨迹准确度分析方法
CN112797931A (zh) * 2020-12-28 2021-05-14 韩莉莉 工业机器人位姿准确度和位姿重复性检测装置及检测方法

Similar Documents

Publication Publication Date Title
CN111660295B (zh) 工业机器人绝对精度标定系统及标定方法
CN110695993B (zh) 一种柔性机械臂同步测量方法、系统及装置
CN109813336B (zh) 惯性测量单元标定方法
CN110370314A (zh) 基于多激光跟踪仪的双臂机器人性能测量方法
RU2466858C1 (ru) Способ контроля точности контурных перемещений промышленных роботов
CN111409108B (zh) 工业机器人测量坐标系与机器人的指令坐标系转换方法
CN109514549A (zh) 一种可实现六自由度的tcp在线快速标定方法及装置
CN102636140A (zh) 测量空间六自由度运动的伸缩拉杆式并联装置
CN111409105A (zh) 工业机器人位姿特性中准确度、重复性的测量方法
CN111452048B (zh) 多机器人相对空间位置关系的标定方法和装置
CN106289378A (zh) 基于激光跟踪仪的工业机器人性能测量方法
CN111409106A (zh) 工业机器人距离性能测试方法
CN110253574A (zh) 一种多任务机械臂位姿检测和误差补偿方法
CN113843792A (zh) 一种手术机器人的手眼标定方法
CN113319855B (zh) 一种多关节诊疗机器人柔顺力控制模式下的重力补偿方法
CN111409103A (zh) 工业机器人位姿特性中互换性的测量方法
CN114029982A (zh) 一种相机在机器臂外的手眼标定装置及其标定方法
Lange et al. Learning accurate path control of industrial robots with joint elasticity
CN111409109B (zh) 工业机器人稳定性的测量和评估方法
JPH0445842B2 (zh)
CN113878586B (zh) 机器人运动学标定装置、方法及系统
CN111421573B (zh) 一种工业机器人拐角偏差性能测试方法
CN111336969A (zh) 工业机器人位姿特性中多方向位姿准确度变动的测量方法
CN111409104B (zh) 工业机器人位姿特性中位姿特性漂移的测量方法
CN112487615B (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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20200714