CN113878586A - 机器人运动学标定装置、方法及系统 - Google Patents

机器人运动学标定装置、方法及系统 Download PDF

Info

Publication number
CN113878586A
CN113878586A CN202111300417.6A CN202111300417A CN113878586A CN 113878586 A CN113878586 A CN 113878586A CN 202111300417 A CN202111300417 A CN 202111300417A CN 113878586 A CN113878586 A CN 113878586A
Authority
CN
China
Prior art keywords
robot
calibration
coordinate system
flange
relative
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.)
Granted
Application number
CN202111300417.6A
Other languages
English (en)
Other versions
CN113878586B (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.)
Shanghai Jingwu Intelligent Technology Co Ltd
Original Assignee
Shanghai Jingwu Intelligent Technology 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 Shanghai Jingwu Intelligent Technology Co Ltd filed Critical Shanghai Jingwu Intelligent Technology Co Ltd
Priority to CN202111300417.6A priority Critical patent/CN113878586B/zh
Publication of CN113878586A publication Critical patent/CN113878586A/zh
Application granted granted Critical
Publication of CN113878586B publication Critical patent/CN113878586B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B25HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
    • B25JMANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
    • B25J9/00Programme-controlled manipulators
    • B25J9/16Programme controls
    • B25J9/1628Programme controls characterised by the control loop
    • B25J9/1653Programme controls characterised by the control loop parameters identification, estimation, stiffness, accuracy, error analysis
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B25HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
    • B25JMANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
    • B25J9/00Programme-controlled manipulators
    • B25J9/16Programme controls
    • B25J9/1656Programme controls characterised by programming, planning systems for manipulators
    • B25J9/1664Programme controls characterised by programming, planning systems for manipulators characterised by motion, path, trajectory planning
    • 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

本发明提供了一种机器人运动学标定装置、方法及系统,包括机器人本体,所述机器人本体的末端设置有法兰,所述法兰上设置有标定传感装置,所述标定传感装置上安装有三个激光位移传感器;所述标定传感装置的一端安装在法兰上,另一端沿标定传感装置的圆周均匀分布三个激光位移传感器。本发明的用于机器人的运动学标定装置和方法,可以对机器人运动学参数进行校正,提高其工作的绝度定位精度。非接触式测量,精度较高;使用三个激光位移传感器和标定球,成本很低,操作简便,可以广泛应用于中小企业;采用的组成元件或设备均可从现有成熟商业产品中选择,在实际中具有良好的可行性。

Description

机器人运动学标定装置、方法及系统
技术领域
本发明涉及机器人标定的技术领域,具体地,涉及机器人运动学标定装置、方法及系统,尤其涉及一种机器人运动学标定装置和方法。
背景技术
随着机器人技术在各行各业中的广泛应用,对机器人的重复定位精度和绝对定位精度要求也越来越高。目前机器人重复定位精度比较高,绝对定位精度却很低,并且每台机器人的差别较大,严重制约了机器人的应用范围。机器人运动学标定一般包括四个步骤,建模,测量,误差辨识,补偿。目前数据测量一般要借助激光跟踪仪、三坐标测量机、球杆仪等昂贵的精密测量仪器,而且需要专业人员操作。
在公告号为CN107042528B的专利文献中公开了一种工业机器人的运动学标定系统及方法,包括放置在工作区域内的测量目标、机器人和设于机器人上的末端执行器,以及计算机。计算机设有读取机器人关节角数据的第一接收模块、读取三个位移传感器检测数据的第二接收模块、根据第一接收模块数据和第二接收模块数据来计算测量目标的名义坐标位置的第一计算模块、以及根据测量目标的名义距离与实际距离的误差来计算机器人补偿数据的第二计算模块。
在公开号为CN107042527A的专利文献中公开了一种基于三维力传感器的工业机器人标定装置及标定方法,装置有标定测量组件安装在机器人末端法兰上,标定球组件固定在机器人工作空间内,标定测量组件有连接底板,固定安装在连接底板一端的三维力传感器,通过螺栓固定安装在三维力传感器上端的测量球结构,连接底板的另一端形成有法兰安装结构;标定球组件有固定底座,固定安装在固定底座上端的连接件,通过螺栓固定安装在连接件上端且与测量球结构相对应的标定球结构。方法:安装标定测量组件和标定球组件;多次操作工业机器人使标定测量组件和标定球组件接触,记录每一次的接触力,根据记录建立标定约束方程,求得工业机器人运动学参数误差向量。
在公开号为CN110978059A的专利文献中公开了一种便携式六轴机械手标定装置,机械手端部设置标定传感装置;标定装置设有固定安装的标定球装置,标定球固定安装在标定球装置上;标定传感装置上设置四个或四个以上的激光位移传感器,通过激光位移传感器获取标定球球心在机械手处于不同位置和姿态下的坐标参数。
针对上述中的相关技术,发明人认为现有运动学标定存在设备昂贵、操作复杂的技术问题,因此,需要提出一种技术方案以改善上述技术问题。
发明内容
针对现有技术中的缺陷,本发明的目的是提供一种机器人运动学标定装置、方法及系统。
根据本发明提供的一种机器人运动学标定装置,包括机器人本体,所述机器人本体的末端设置有法兰,所述法兰上设置有标定传感装置,所述标定传感装置上安装有三个激光位移传感器。
优选地,所述标定传感装置的一端安装在法兰上,另一端固定三个激光位移传感器,激光位移传感器沿标定传感装置的圆周均匀分布。
本发明还提供一种机器人运动学标定方法,所述方法应用上述中的机器人运动学标定装置,所述方法包括如下步骤:
步骤S1:将标定传感装置安装在机器人末端法兰上;
步骤S2:将标定球固定在机器人工作空间内;
步骤S3:操作机器人移动,使三个激光位移传感器的激光光束均打在标定球上,记录三个激光位移传感器的读数l=(l1,l2,l3),同时记录机器人关节角度θi,根据标定球半径r计算标定球球心相对于机器人法兰坐标系的位置,其中l1、l2、l3分别为三个激光位移传感器得读数,表示传感器到标定球表面的距离;l表示记录三个激光位移传感器读数;
步骤S4:保持标定球位置不动,再次操作机器人移动,机器人处于不同的位姿,且三个激光位移传感器的激光光束均打在标定球上,记录三个激光位移传感器的读数(l1,l2,l3),同时记录机器人关节角度θi,根据标定球半径r计算标定球球心相对于机器人法兰坐标系的位置;
步骤S5:多次改变标定球在机器人工作空间的位置,重复步骤S3和S4,得到多组测量数据;
步骤S6:采用D-H法建立机器人的运动学模型,将连杆坐标系i-1到连杆坐标系i的齐次变换矩阵表示为
Figure BDA0003338163240000021
表示机器人连杆坐标系i-1到连杆坐标系i的齐次变换矩阵,则机器人末端法兰坐标系n相对于机器人基坐标系0的齐次变换矩阵为
Figure BDA0003338163240000031
Figure BDA0003338163240000032
表示机器人法兰坐标系n相对于基坐标系0的齐次变换矩阵根据标定球球心相对于机器人法兰的位置
Figure BDA0003338163240000033
表示标定球球心s相对与机器人法兰坐标系n的偏置,x表示标定球球心相对与机器人法兰坐标系x方向的偏置,y表示标定球球心相对与机器人法兰坐标系y方向的偏置,z表示标定球球心相对与机器人法兰坐标系z方向的偏置;得到标定球球心坐标系s相对于机器人基坐标系0的齐次变换矩阵为
Figure BDA0003338163240000034
表示标定球球心坐标系s相对于机器人基坐标系0的齐次变换矩阵,I为3×3的单位矩阵。
优选地,所述步骤S3中根据标定球半径r计算标定球球心相对于机器人法兰坐标系的位置的计算方法如下:
标定球半径r,激光位移传感器相对于机器人法兰坐标系的初始坐标分别为(xi,yi,zi),i=1,2,3,其中xi表示激光位移传感器在机器人法兰坐标系下x方向的偏置;yi表示激光位移传感器在机器人法兰坐标系下y方向的偏置;zi表示激光位移传感器在机器人法兰坐标系下z方向的偏置;i表示激光位移传感器编号;则此时标定球表面激光点的坐标为(xi,yi,zi-li),i=1,2,3,其中li表示第i个激光位移传感器的读数,对应l1、l2、l3;标定球球心相对于法兰坐标系的坐标(x,y,z)通过求解下面的方程组得到,其中x表示标定球球心相对于法兰坐标x方向的偏置,y表示标定球球心相对于法兰坐标y方向的偏置,z表示标定球球心相对于法兰坐标z方向的偏置:
(xi-x)2+(yi-y)2+(zi-z-li)2=r2,i=1,2,3
通过求解得到两组解,根据标定球表面激光点的坐标求得平面方程Ax+By+Cz+D=0,其中A、B、C、D是平面方程的表示参数,带入两组解,选择Ax+By+Cz+D>0的一组解。
优选地,所述步骤S6提取标定球球心在机器人基坐标系的位置分量为pn=f(θi,l),pn表示标定球球心在机器人基坐标系0的计算位置,θi表示记录的机器人关节角,l表示激光位移传感器读数,考虑机器人运动学误差,机器人标定球球心在机器人基坐标系的实际位置为p=pn+J*σ,p表示机器人标定球球心在机器人基坐标系的实际位置,J为标定雅可比矩阵,σ为运动学参数误差向量。
优选地,当操作机器人得到一组测量数据,得到两个位置方程:
Figure BDA0003338163240000035
Figure BDA0003338163240000036
表示标定球固定在一个地方的测量第一次所列方程,
Figure BDA0003338163240000037
表示标定球固定在一个地方的测量第二次所列方程,根据两次测量标定球球心在基坐标系的位置相同,得到
Figure BDA0003338163240000041
多次改变标定球得位置,得到m组数据,建立标定方程:A=B*σ,A代表m组计算位置差的向量,B代表m组雅可比矩阵差的向量,其中
Figure BDA0003338163240000042
B=[(J2-J1)1,(J2-J1)2,…,(J2-J1)m]。
优选地,根据最小二乘法计算运动学误差向量σ=(BTB)-1BTA。
本发明还提供一种机器人运动学标定系统,所述系统应用上述中的机器人运动学标定装置,所述系统包括如下模块:
模块M1:将标定传感装置安装在机器人末端法兰上;
模块M2:将标定球固定在机器人工作空间内;
模块M3:操作机器人移动,使三个激光位移传感器的激光光束均打在标定球上,记录三个激光位移传感器的读数l=(l1,l2,l3),同时记录机器人关节角度θi,根据标定球半径r计算标定球球心相对于机器人法兰坐标系的位置,其中l1、l2、l3分别为三个激光位移传感器得读数,表示传感器到标定球表面的距离;l表示记录三个激光位移传感器读数;
模块M4:保持标定球位置不动,再次操作机器人移动,机器人处于不同的位姿,且三个激光位移传感器的激光光束均打在标定球上,记录三个激光位移传感器的读数(l1,l2,l3),同时记录机器人关节角度θi,根据标定球半径r计算标定球球心相对于机器人法兰坐标系的位置;
模块M5:多次改变标定球在机器人工作空间的位置,重复模块M3和S4,得到多组测量数据;
模块M6:采用D-H法建立机器人的运动学模型,将连杆坐标系i-1到连杆坐标系i的齐次变换矩阵表示为
Figure BDA0003338163240000043
表示机器人连杆坐标系i-1到连杆坐标系i的齐次变换矩阵,则机器人末端法兰坐标系n相对于机器人基坐标系0的齐次变换矩阵为
Figure BDA0003338163240000044
Figure BDA0003338163240000045
表示机器人法兰坐标系n相对于基坐标系0的齐次变换矩阵根据标定球球心相对于机器人法兰的位置
Figure BDA0003338163240000046
表示标定球球心s相对与机器人法兰坐标系n的偏置,x表示标定球球心相对与机器人法兰坐标系x方向的偏置,y表示标定球球心相对与机器人法兰坐标系y方向的偏置,z表示标定球球心相对与机器人法兰坐标系z方向的偏置;得到标定球球心坐标系s相对于机器人基坐标系0的齐次变换矩阵为
Figure BDA0003338163240000047
表示标定球球心坐标系s相对于机器人基坐标系0的齐次变换矩阵,I为3×3的单位矩阵。
优选地,所述模块M3中根据标定球半径r计算标定球球心相对于机器人法兰坐标系的位置的计算系统如下:
标定球半径r,激光位移传感器相对于机器人法兰坐标系的初始坐标分别为(xi,yi,zi),i=1,2,3,其中xi表示激光位移传感器在机器人法兰坐标系下x方向的偏置;yi表示激光位移传感器在机器人法兰坐标系下y方向的偏置;zi表示激光位移传感器在机器人法兰坐标系下z方向的偏置;i表示激光位移传感器编号;则此时标定球表面激光点的坐标为(xi,yi,zi-li),i=1,2,3,其中li表示第i个激光位移传感器的读数,对应l1、l2、l3;标定球球心相对于法兰坐标系的坐标(x,y,z)通过求解下面的方程组得到,其中x表示标定球球心相对于法兰坐标x方向的偏置,y表示标定球球心相对于法兰坐标y方向的偏置,z表示标定球球心相对于法兰坐标z方向的偏置:
(xi-x)2+(yi-y)2+(zi-z-li)2=r2,i=1,2,3
通过求解得到两组解,根据标定球表面激光点的坐标求得平面方程Ax+By+Cz+D=0,其中A、B、C、D是平面方程的表示参数,带入两组解,选择Ax+By+Cz+D>0的一组解;
所述模块M6提取标定球球心在机器人基坐标系的位置分量为pn=f(θi,l),pn表示标定球球心在机器人基坐标系0的计算位置,θi表示记录的机器人关节角,l表示激光位移传感器读数,考虑机器人运动学误差,机器人标定球球心在机器人基坐标系的实际位置为p=pn+J*σ,p表示机器人标定球球心在机器人基坐标系的实际位置,J为标定雅可比矩阵,σ为运动学参数误差向量。
优选地,当操作机器人得到一组测量数据,得到两个位置方程:
Figure BDA0003338163240000051
Figure BDA0003338163240000052
表示标定球固定在一个地方的测量第一次所列方程,
Figure BDA0003338163240000053
表示标定球固定在一个地方的测量第二次所列方程,根据两次测量标定球球心在基坐标系的位置相同,得到
Figure BDA0003338163240000054
多次改变标定球得位置,得到m组数据,建立标定方程:A=B*σ,A代表m组计算位置差的向量,B代表m组雅可比矩阵差的向量,其中
Figure BDA0003338163240000055
B=[(J2-J1)1,(J2-J1)2,…,(J2-J1)m];根据最小二乘法计算运动学误差向量σ=(BTB)-1BTA。
与现有技术相比,本发明具有如下的有益效果:
1、本发明的用于机器人的运动学标定装置和方法,可以对机器人运动学参数进行校正,提高其工作的绝度定位精度。非接触式测量,精度较高;
2、本发明使用三个激光位移传感器和标定球,成本很低,操作简便,可以广泛应用于中小企业;
3、本发明采用的组成元件或设备均可从现有成熟商业产品中选择,在实际中具有良好的可行性。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:
图1为本发明的标定装置示意图。
其中:
机器人1 标定球装置3
标定传感装置2 激光位移传感器4
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变化和改进。这些都属于本发明的保护范围。
本发明提供的一种机器人运动学标定装置、方法及系统,包括安装在机器人末端法兰的标定传感装置,机器人工作区域内标定球,其中标定传感装置安装三个激光位移传感器,通过激光位移传感器获取标定球球心在机器人不同位姿下的坐标。
参照图1,标定传感装置,一端安装在机器人法兰上,另一端沿装置圆周均匀分布固定三个激光位移传感器。标定球放置在机器人工作区域,确保移动机器人,激光位移传感器能够测到球表面。
本发明提供的一种机器人运动学标定方法,包括以下步骤:
步骤S1:将标定传感装置安装在机器人末端法兰上;步骤S2:将标定球固定在机器人工作空间内。
步骤S3:操作机器人移动,使三个激光位移传感器的激光光束均能够打在标定球上,记录三个激光位移传感器的读数l=(l1,l2,l3),同时记录机器人关节角度θi,根据标定球半径r可以计算标定球球心相对于机器人法兰坐标系的位置,其中l1、l2、l3分别为三个激光位移传感器得读数,表示传感器到标定球表面的距离;l表示记录三个激光位移传感器读数,计算方法如下:
标定球半径r,激光位移传感器相对于机器人法兰坐标系的初始坐标分别为(xi,yi,zi),i=1,2,3,其中xi表示激光位移传感器在机器人法兰坐标系下x方向的偏置;yi表示激光位移传感器在机器人法兰坐标系下y方向的偏置;zi表示激光位移传感器在机器人法兰坐标系下z方向的偏置;i表示激光位移传感器编号;则此时标定球表面激光点的坐标为(xi,yi,zi-li),i=1,2,3,其中li表示第i个激光位移传感器的读数,对应l1、l2、l3;标定球球心相对于法兰坐标系的坐标(x,y,z)可以通过求解下面的方程组得到,其中x表示标定球球心相对于法兰坐标x方向的偏置,y表示标定球球心相对于法兰坐标y方向的偏置,z表示标定球球心相对于法兰坐标z方向的偏置:
(xi-x)2+(yi-y)2+(zi-z-li)2=r2,i=1,2,3
通过求解可以得到两组解,根据标定球表面激光点的坐标求得平面方程Ax+By+Cz+D=0,其中A、B、C、D是平面方程的表示参数,带入两组解,选择Ax+By+Cz+D>0的一组解。
步骤S4:保持标定球位置不动,再次操作机器人移动,机器人处于不同的位姿,且三个激光位移传感器的激光光束均能够打在标定球上,记录三个激光位移传感器的读数(l1,l2,l3),同时记录机器人关节角度θi,根据标定球半径r计算标定球球心相对于机器人法兰坐标系的位置。
步骤S5:多次改变标定球在机器人工作空间的位置,重复步骤S3和S4,得到多组测量数据。
步骤S6:采用D-H法建立机器人的运动学模型,将连杆坐标系i-1到连杆坐标系i的齐次变换矩阵表示为
Figure BDA0003338163240000071
表示机器人连杆坐标系i-1到连杆坐标系i的齐次变换矩阵,则机器人末端法兰坐标系n相对于机器人基坐标系0的齐次变换矩阵为
Figure BDA0003338163240000072
Figure BDA0003338163240000073
表示机器人法兰坐标系n相对于基坐标系0的齐次变换矩阵根据标定球球心相对于机器人法兰的位置
Figure BDA0003338163240000074
表示标定球球心s相对与机器人法兰坐标系n的偏置,x表示标定球球心相对与机器人法兰坐标系x方向的偏置,y表示标定球球心相对与机器人法兰坐标系y方向的偏置,z表示标定球球心相对与机器人法兰坐标系z方向的偏置;得到标定球球心坐标系s相对于机器人基坐标系0的齐次变换矩阵为
Figure BDA0003338163240000075
表示标定球球心坐标系s相对于机器人基坐标系0的齐次变换矩阵,I为3×3的单位矩阵。提取标定球球心在机器人基坐标系的位置分量为pn=f(θi,l),pn表示标定球球心在机器人基坐标系0的计算位置,θi表示记录的机器人关节角,l表示激光位移传感器读数,考虑机器人运动学误差,机器人标定球球心在机器人基坐标系的实际位置为p=pn+J*σ,p表示机器人标定球球心在机器人基坐标系的实际位置,J为标定雅可比矩阵,σ为运动学参数误差向量。当操作机器人得到一组测量数据,得到两个位置方程:
Figure BDA0003338163240000087
表示标定球固定在一个地方的测量第一次所列方程,
Figure BDA0003338163240000082
表示标定球固定在一个地方的测量第二次所列方程,根据两次测量标定球球心在基坐标系的位置相同,得到
Figure BDA0003338163240000083
Figure BDA0003338163240000084
多次改变标定球得位置,得到m组数据,建立标定方程:A=B*σ,A代表m组计算位置差的向量,B代表m组雅可比矩阵差的向量,其中
Figure BDA0003338163240000085
Figure BDA0003338163240000086
B=[(J2-J1)1,(J2-J1)2,…,(J2-J1)m]。根据最小二乘法计算运动学误差向量σ=(BTB)-1BTA。
本发明还提供一种机器人运动学标定系统,所述系统应用上述中的机器人运动学标定装置,所述系统包括如下模块:
模块M1:将标定传感装置安装在机器人末端法兰上;模块M2:将标定球固定在机器人工作空间内。
模块M3:操作机器人移动,使三个激光位移传感器的激光光束均打在标定球上,记录三个激光位移传感器的读数l=(l1,l2,l3),同时记录机器人关节角度θi,根据标定球半径r计算标定球球心相对于机器人法兰坐标系的位置,其中l1、l2、l3分别为三个激光位移传感器得读数,表示传感器到标定球表面的距离;l表示记录三个激光位移传感器读数;标定球半径r,激光位移传感器相对于机器人法兰坐标系的初始坐标分别为(xi,yi,zi),i=1,2,3,其中xi表示激光位移传感器在机器人法兰坐标系下x方向的偏置;yi表示激光位移传感器在机器人法兰坐标系下y方向的偏置;zi表示激光位移传感器在机器人法兰坐标系下z方向的偏置;i表示激光位移传感器编号;则此时标定球表面激光点的坐标为(xi,yi,zi-li),i=1,2,3,其中li表示第i个激光位移传感器的读数,对应l1、l2、l3;标定球球心相对于法兰坐标系的坐标(x,y,z)通过求解下面的方程组得到,其中x表示标定球球心相对于法兰坐标x方向的偏置,y表示标定球球心相对于法兰坐标y方向的偏置,z表示标定球球心相对于法兰坐标z方向的偏置:
(xi-x)2+(yi-y)2+(zi-z-li)2=r2,i=1,2,3
通过求解得到两组解,根据标定球表面激光点的坐标求得平面方程Ax+By+Cz+D=0,其中A、B、C、D是平面方程的表示参数,带入两组解,选择Ax+By+Cz+D>0的一组解。
模块M4:保持标定球位置不动,再次操作机器人移动,机器人处于不同的位姿,且三个激光位移传感器的激光光束均打在标定球上,记录三个激光位移传感器的读数(l1,l2,l3),同时记录机器人关节角度θi,根据标定球半径r计算标定球球心相对于机器人法兰坐标系的位置;模块M5:多次改变标定球在机器人工作空间的位置,重复模块M3和S4,得到多组测量数据。
模块M6:采用D-H法建立机器人的运动学模型,将连杆坐标系i-1到连杆坐标系i的齐次变换矩阵表示为
Figure BDA0003338163240000091
表示机器人连杆坐标系i-1到连杆坐标系i的齐次变换矩阵,则机器人末端法兰坐标系n相对于机器人基坐标系0的齐次变换矩阵为
Figure BDA0003338163240000092
Figure BDA0003338163240000093
表示机器人法兰坐标系n相对于基坐标系0的齐次变换矩阵根据标定球球心相对于机器人法兰的位置
Figure BDA0003338163240000094
表示标定球球心s相对与机器人法兰坐标系n的偏置,x表示标定球球心相对与机器人法兰坐标系x方向的偏置,y表示标定球球心相对与机器人法兰坐标系y方向的偏置,z表示标定球球心相对与机器人法兰坐标系z方向的偏置;得到标定球球心坐标系s相对于机器人基坐标系0的齐次变换矩阵为
Figure BDA0003338163240000095
表示标定球球心坐标系s相对于机器人基坐标系0的齐次变换矩阵,I为3×3的单位矩阵。提取标定球球心在机器人基坐标系的位置分量为pn=f(θi,l),pn表示标定球球心在机器人基坐标系0的计算位置,θi表示记录的机器人关节角,l表示激光位移传感器读数,考虑机器人运动学误差,机器人标定球球心在机器人基坐标系的实际位置为p=pn+J*σ,p表示机器人标定球球心在机器人基坐标系的实际位置,J为标定雅可比矩阵,σ为运动学参数误差向量。
当操作机器人得到一组测量数据,得到两个位置方程:
Figure BDA0003338163240000096
Figure BDA0003338163240000097
表示标定球固定在一个地方的测量第一次所列方程,
Figure BDA0003338163240000098
表示标定球固定在一个地方的测量第二次所列方程,根据两次测量标定球球心在基坐标系的位置相同,得到
Figure BDA0003338163240000099
多次改变标定球得位置,得到m组数据,建立标定方程:A=B*σ,A代表m组计算位置差的向量,B代表m组雅可比矩阵差的向量,其中
Figure BDA00033381632400000910
Figure BDA00033381632400000911
B=[(J2-J1)1,(J2-J1)2,…,(J2-J1)m];根据最小二乘法计算运动学误差向量σ=(BTB)-1BTA。
本发明的用于机器人的运动学标定装置和方法,可以对机器人运动学参数进行校正,提高其工作的绝度定位精度。非接触式测量,精度较高;本发明使用三个激光位移传感器和标定球,成本很低,操作简便,可以广泛应用于中小企业;本发明采用的组成元件或设备均可从现有成熟商业产品中选择,在实际中具有良好的可行性。
本领域技术人员知道,除了以纯计算机可读程序代码方式实现本发明提供的系统及其各个装置、模块、单元以外,完全可以通过将方法步骤进行逻辑编程来使得本发明提供的系统及其各个装置、模块、单元以逻辑门、开关、专用集成电路、可编程逻辑控制器以及嵌入式微控制器等的形式来实现相同功能。所以,本发明提供的系统及其各项装置、模块、单元可以被认为是一种硬件部件,而对其内包括的用于实现各种功能的装置、模块、单元也可以视为硬件部件内的结构;也可以将用于实现各种功能的装置、模块、单元视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
在本申请的描述中,需要理解的是,术语“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本申请和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本申请的限制。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变化或修改,这并不影响本发明的实质内容。在不冲突的情况下,本申请的实施例和实施例中的特征可以任意相互组合。

Claims (10)

1.一种机器人运动学标定装置,其特征在于,包括机器人本体,所述机器人本体的末端设置有法兰,所述法兰上设置有标定传感装置,所述标定传感装置上安装有三个激光位移传感器。
2.根据权利要求1所述的机器人运动学标定装置,其特征在于,所述标定传感装置的一端安装在法兰上,另一端固定三个激光位移传感器,激光位移传感器沿标定传感装置的圆周均匀分布。
3.一种机器人运动学标定方法,其特征在于,所述方法应用如权利要求1-2任一项所述的机器人运动学标定装置,所述方法包括如下步骤:
步骤S1:将标定传感装置安装在机器人末端法兰上;
步骤S2:将标定球固定在机器人工作空间内;
步骤S3:操作机器人移动,使三个激光位移传感器的激光光束均打在标定球上,记录三个激光位移传感器的读数l=(l1,l2,l3),同时记录机器人关节角度θi,根据标定球半径r计算标定球球心相对于机器人法兰坐标系的位置,其中l1、l2、l3分别为三个激光位移传感器得读数,表示传感器到标定球表面的距离;l表示记录三个激光位移传感器读数;
步骤S4:保持标定球位置不动,再次操作机器人移动,机器人处于不同的位姿,且三个激光位移传感器的激光光束均打在标定球上,记录三个激光位移传感器的读数(l1,l2,l3),同时记录机器人关节角度θi,根据标定球半径r计算标定球球心相对于机器人法兰坐标系的位置;
步骤S5:多次改变标定球在机器人工作空间的位置,重复步骤S3和S4,得到多组测量数据;
步骤S6:采用D-H法建立机器人的运动学模型,将连杆坐标系i-1到连杆坐标系i的齐次变换矩阵表示为
Figure FDA0003338163230000011
表示机器人连杆坐标系i-1到连杆坐标系i的齐次变换矩阵,则机器人末端法兰坐标系n相对于机器人基坐标系0的齐次变换矩阵为
Figure FDA0003338163230000012
Figure FDA0003338163230000013
Figure FDA0003338163230000014
表示机器人法兰坐标系n相对于基坐标系0的齐次变换矩阵根据标定球球心相对于机器人法兰的位置
Figure FDA0003338163230000015
Figure FDA0003338163230000016
表示标定球球心s相对与机器人法兰坐标系n的偏置,x表示标定球球心相对与机器人法兰坐标系x方向的偏置,y表示标定球球心相对与机器人法兰坐标系y方向的偏置,z表示标定球球心相对与机器人法兰坐标系z方向的偏置;得到标定球球心坐标系s相对于机器人基坐标系0的齐次变换矩阵为
Figure FDA0003338163230000021
Figure FDA0003338163230000022
表示标定球球心坐标系s相对于机器人基坐标系0的齐次变换矩阵,I为3×3的单位矩阵。
4.根据权利要求3所述的机器人运动学标定方法,其特征在于,所述步骤S3中根据标定球半径r计算标定球球心相对于机器人法兰坐标系的位置的计算方法如下:
标定球半径r,激光位移传感器相对于机器人法兰坐标系的初始坐标分别为(xi,yi,zi),i=1,2,3,其中xi表示激光位移传感器在机器人法兰坐标系下x方向的偏置;yi表示激光位移传感器在机器人法兰坐标系下y方向的偏置;zi表示激光位移传感器在机器人法兰坐标系下z方向的偏置;i表示激光位移传感器编号;则此时标定球表面激光点的坐标为(xi,yi,zi-li),i=1,2,3,其中li表示第i个激光位移传感器的读数,对应l1、l2、l3;标定球球心相对于法兰坐标系的坐标(x,y,z)通过求解下面的方程组得到,其中x表示标定球球心相对于法兰坐标x方向的偏置,y表示标定球球心相对于法兰坐标y方向的偏置,z表示标定球球心相对于法兰坐标z方向的偏置:
(xi-x)2+(yi-y)2+(zi-z-li)2=r2,i=1,2,3
通过求解得到两组解,根据标定球表面激光点的坐标求得平面方程Ax+By+Cz+D=0,其中A、B、C、D是平面方程的表示参数,带入两组解,选择Ax+By+Cz+D>0的一组解。
5.根据权利要求3所述的机器人运动学标定方法,其特征在于,所述步骤S6提取标定球球心在机器人基坐标系的位置分量为pn=f(θi,l),pn表示标定球球心在机器人基坐标系0的计算位置,θi表示记录的机器人关节角,l表示激光位移传感器读数,考虑机器人运动学误差,机器人标定球球心在机器人基坐标系的实际位置为p=pn+J*σ,p表示机器人标定球球心在机器人基坐标系的实际位置,J为标定雅可比矩阵,σ为运动学参数误差向量。
6.根据权利要求5所述的机器人运动学标定方法,其特征在于,当操作机器人得到一组测量数据,得到两个位置方程:
Figure FDA0003338163230000023
表示标定球固定在一个地方的测量第一次所列方程,
Figure FDA0003338163230000024
表示标定球固定在一个地方的测量第二次所列方程,根据两次测量标定球球心在基坐标系的位置相同,得到
Figure FDA0003338163230000025
Figure FDA0003338163230000026
多次改变标定球得位置,得到m组数据,建立标定方程:A=B*σ,A代表m组计算位置差的向量,B代表m组雅可比矩阵差的向量,其中
Figure FDA0003338163230000027
Figure FDA0003338163230000028
B=[(J2-J1)1,(J2-J1)2,...,(J2-J1)m]。
7.根据权利要求6所述的机器人运动学标定方法,其特征在于,根据最小二乘法计算运动学误差向量σ=(BTB)-1BTA。
8.一种机器人运动学标定系统,其特征在于,所述系统应用如权利要求1-2任一项所述的机器人运动学标定装置,所述系统包括如下模块:
模块M1:将标定传感装置安装在机器人末端法兰上;
模块M2:将标定球固定在机器人工作空间内;
模块M3:操作机器人移动,使三个激光位移传感器的激光光束均打在标定球上,记录三个激光位移传感器的读数l=(l1,l2,l3),同时记录机器人关节角度θi,根据标定球半径r计算标定球球心相对于机器人法兰坐标系的位置,其中l1、l2、l3分别为三个激光位移传感器得读数,表示传感器到标定球表面的距离;l表示记录三个激光位移传感器读数;
模块M4:保持标定球位置不动,再次操作机器人移动,机器人处于不同的位姿,且三个激光位移传感器的激光光束均打在标定球上,记录三个激光位移传感器的读数(l1,l2,l3),同时记录机器人关节角度θi,根据标定球半径r计算标定球球心相对于机器人法兰坐标系的位置;
模块M5:多次改变标定球在机器人工作空间的位置,重复模块M3和S4,得到多组测量数据;
模块M6:采用D-H法建立机器人的运动学模型,将连杆坐标系i-1到连杆坐标系i的齐次变换矩阵表示为
Figure FDA0003338163230000031
表示机器人连杆坐标系i-1到连杆坐标系i的齐次变换矩阵,则机器人末端法兰坐标系n相对于机器人基坐标系0的齐次变换矩阵为
Figure FDA0003338163230000032
Figure FDA0003338163230000033
Figure FDA0003338163230000034
表示机器人法兰坐标系n相对于基坐标系0的齐次变换矩阵根据标定球球心相对于机器人法兰的位置
Figure FDA0003338163230000035
Figure FDA0003338163230000036
表示标定球球心s相对与机器人法兰坐标系n的偏置,x表示标定球球心相对与机器人法兰坐标系x方向的偏置,y表示标定球球心相对与机器人法兰坐标系y方向的偏置,z表示标定球球心相对与机器人法兰坐标系z方向的偏置;得到标定球球心坐标系s相对于机器人基坐标系0的齐次变换矩阵为
Figure FDA0003338163230000037
Figure FDA0003338163230000038
表示标定球球心坐标系s相对于机器人基坐标系0的齐次变换矩阵,I为3×3的单位矩阵。
9.根据权利要求8所述的机器人运动学标定系统,其特征在于,所述模块M3中根据标定球半径r计算标定球球心相对于机器人法兰坐标系的位置的计算系统如下:
标定球半径r,激光位移传感器相对于机器人法兰坐标系的初始坐标分别为(xi,yi,zi),i=1,2,3,其中xi表示激光位移传感器在机器人法兰坐标系下x方向的偏置;yi表示激光位移传感器在机器人法兰坐标系下y方向的偏置;zi表示激光位移传感器在机器人法兰坐标系下z方向的偏置;i表示激光位移传感器编号;则此时标定球表面激光点的坐标为(xi,yi,zi-li),i=1,2,3,其中li表示第i个激光位移传感器的读数,对应l1、l2、l3;标定球球心相对于法兰坐标系的坐标(x,y,z)通过求解下面的方程组得到,其中x表示标定球球心相对于法兰坐标x方向的偏置,y表示标定球球心相对于法兰坐标y方向的偏置,z表示标定球球心相对于法兰坐标z方向的偏置:
(xi-x)2+(yi-y)2+(zi-z-li)2=r2,i=1,2,3
通过求解得到两组解,根据标定球表面激光点的坐标求得平面方程Ax+By+Cz+D=0,其中A、B、C、D是平面方程的表示参数,带入两组解,选择Ax+By+Cz+D>0的一组解;
所述模块M6提取标定球球心在机器人基坐标系的位置分量为pn=f(θi,l),pn表示标定球球心在机器人基坐标系0的计算位置,θi表示记录的机器人关节角,l表示激光位移传感器读数,考虑机器人运动学误差,机器人标定球球心在机器人基坐标系的实际位置为p=pn+J*σ,p表示机器人标定球球心在机器人基坐标系的实际位置,J为标定雅可比矩阵,σ为运动学参数误差向量。
10.根据权利要求9所述的机器人运动学标定系统,其特征在于,当操作机器人得到一组测量数据,得到两个位置方程:
Figure FDA0003338163230000041
表示标定球固定在一个地方的测量第一次所列方程,
Figure FDA0003338163230000042
表示标定球固定在一个地方的测量第二次所列方程,根据两次测量标定球球心在基坐标系的位置相同,得到
Figure FDA0003338163230000043
Figure FDA0003338163230000044
多次改变标定球得位置,得到m组数据,建立标定方程:A=B*σ,A代表m组计算位置差的向量,B代表m组雅可比矩阵差的向量,其中
Figure FDA0003338163230000045
Figure FDA0003338163230000046
B=[(J2-J1)1,(J2-J1)2,...,(J2-J1)m];根据最小二乘法计算运动学误差向量σ=(BTB)-1BTA。
CN202111300417.6A 2021-11-04 2021-11-04 机器人运动学标定装置、方法及系统 Active CN113878586B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111300417.6A CN113878586B (zh) 2021-11-04 2021-11-04 机器人运动学标定装置、方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111300417.6A CN113878586B (zh) 2021-11-04 2021-11-04 机器人运动学标定装置、方法及系统

Publications (2)

Publication Number Publication Date
CN113878586A true CN113878586A (zh) 2022-01-04
CN113878586B CN113878586B (zh) 2023-07-18

Family

ID=79016973

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111300417.6A Active CN113878586B (zh) 2021-11-04 2021-11-04 机器人运动学标定装置、方法及系统

Country Status (1)

Country Link
CN (1) CN113878586B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117140535A (zh) * 2023-10-27 2023-12-01 南湖实验室 一种基于单笔测量的机器人运动学参数标定方法及系统

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105136031A (zh) * 2015-08-25 2015-12-09 华中科技大学 一种五轴联动机床旋转轴的几何误差连续测量方法
US20160239013A1 (en) * 2015-02-12 2016-08-18 The Boeing Company Location Calibration for Automated Production Manufacturing
CN107042528A (zh) * 2017-06-01 2017-08-15 中国科学院宁波材料技术与工程研究所 一种工业机器人的运动学标定系统及方法
DE102016212695A1 (de) * 2016-05-31 2017-11-30 Siemens Aktiengesellschaft Industrieroboter
KR20170137353A (ko) * 2016-06-03 2017-12-13 전남대학교산학협력단 케이블 로봇 캘리브레이션 장치 및 방법
CN109676636A (zh) * 2019-03-06 2019-04-26 南京航空航天大学 一种工业机器人运动学标定系统及标定方法
CN110978059A (zh) * 2019-12-23 2020-04-10 芜湖哈特机器人产业技术研究院有限公司 一种便携式六轴机械手标定装置及其标定方法
CN111390914A (zh) * 2020-04-17 2020-07-10 上海智殷自动化科技有限公司 一种机器人零位和工具坐标标定方法
CN112847341A (zh) * 2020-12-25 2021-05-28 中国科学院宁波材料技术与工程研究所 工业机器人分步式标定系统及方法
CN113029200A (zh) * 2021-03-29 2021-06-25 上海景吾智能科技有限公司 基于机器人传感器测试航向角、精度的方法、系统及介质

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20160239013A1 (en) * 2015-02-12 2016-08-18 The Boeing Company Location Calibration for Automated Production Manufacturing
CN105136031A (zh) * 2015-08-25 2015-12-09 华中科技大学 一种五轴联动机床旋转轴的几何误差连续测量方法
DE102016212695A1 (de) * 2016-05-31 2017-11-30 Siemens Aktiengesellschaft Industrieroboter
KR20170137353A (ko) * 2016-06-03 2017-12-13 전남대학교산학협력단 케이블 로봇 캘리브레이션 장치 및 방법
CN107042528A (zh) * 2017-06-01 2017-08-15 中国科学院宁波材料技术与工程研究所 一种工业机器人的运动学标定系统及方法
CN109676636A (zh) * 2019-03-06 2019-04-26 南京航空航天大学 一种工业机器人运动学标定系统及标定方法
CN110978059A (zh) * 2019-12-23 2020-04-10 芜湖哈特机器人产业技术研究院有限公司 一种便携式六轴机械手标定装置及其标定方法
CN111390914A (zh) * 2020-04-17 2020-07-10 上海智殷自动化科技有限公司 一种机器人零位和工具坐标标定方法
CN112847341A (zh) * 2020-12-25 2021-05-28 中国科学院宁波材料技术与工程研究所 工业机器人分步式标定系统及方法
CN113029200A (zh) * 2021-03-29 2021-06-25 上海景吾智能科技有限公司 基于机器人传感器测试航向角、精度的方法、系统及介质

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117140535A (zh) * 2023-10-27 2023-12-01 南湖实验室 一种基于单笔测量的机器人运动学参数标定方法及系统
CN117140535B (zh) * 2023-10-27 2024-02-02 南湖实验室 一种基于单笔测量的机器人运动学参数标定方法及系统

Also Published As

Publication number Publication date
CN113878586B (zh) 2023-07-18

Similar Documents

Publication Publication Date Title
CN111660295B (zh) 工业机器人绝对精度标定系统及标定方法
CN109822574B (zh) 一种工业机器人末端六维力传感器标定的方法
CN107042528B (zh) 一种工业机器人的运动学标定系统及方法
CN109304730B (zh) 一种基于激光测距仪的机器人运动学参数标定方法
CN110695993B (zh) 一种柔性机械臂同步测量方法、系统及装置
CN110193829B (zh) 一种耦合运动学与刚度参数辨识的机器人精度控制方法
Joubair et al. Non-kinematic calibration of a six-axis serial robot using planar constraints
CN112873199B (zh) 基于运动学与空间插值的机器人绝对定位精度标定方法
Zhuang et al. Robot calibration with planar constraints
Santolaria et al. A self-centering active probing technique for kinematic parameter identification and verification of articulated arm coordinate measuring machines
CN107042527B (zh) 一种基于三维力传感器的工业机器人标定装置及标定方法
CN111203861B (zh) 一种机器人工具坐标系的标定方法及标定系统
CN113386136B (zh) 一种基于标准球阵目标估计的机器人位姿矫正方法及系统
KR101797122B1 (ko) 이동형 3차원 좌표 측정기(cmm)의 에러 측정 및 보상 방법
CN112847341A (zh) 工业机器人分步式标定系统及方法
CN110893619A (zh) 一种基于激光跟踪仪的工业机器人位姿校准装置
CN113146613B (zh) 一种工业机器人d-h参数三维自标定校准装置及方法
Yin et al. Real-time thermal error compensation method for robotic visual inspection system
CN111267143A (zh) 六自由度工业串联机器人关节刚度辨识方法及系统
Guo et al. A measurement method for calibrating kinematic parameters of industrial robots with point constraint by a laser displacement sensor
Santolaria et al. Self-alignment of on-board measurement sensors for robot kinematic calibration
CN113878586B (zh) 机器人运动学标定装置、方法及系统
JPH0445841B2 (zh)
Chen et al. Elasto-geometrical calibration of a hybrid mobile robot considering gravity deformation and stiffness parameter errors
TW202112514A (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
CB02 Change of applicant information

Address after: 311231 building 3, No. 477, Hongxing Road, Qiaonan block, economic and Technological Development Zone, Xiaoshan District, Hangzhou City, Zhejiang Province

Applicant after: Hangzhou Jingwu Intelligent Technology Co.,Ltd.

Address before: Room 12, 3rd floor, No.2 Lane 1446, Yunguan Road, Lingang New District, Pudong New Area pilot Free Trade Zone, Shanghai, 201306

Applicant before: Shanghai Jingwu Intelligent Technology Co.,Ltd.

CB02 Change of applicant information
GR01 Patent grant
GR01 Patent grant