CN113319855B - 一种多关节诊疗机器人柔顺力控制模式下的重力补偿方法 - Google Patents

一种多关节诊疗机器人柔顺力控制模式下的重力补偿方法 Download PDF

Info

Publication number
CN113319855B
CN113319855B CN202110719915.8A CN202110719915A CN113319855B CN 113319855 B CN113319855 B CN 113319855B CN 202110719915 A CN202110719915 A CN 202110719915A CN 113319855 B CN113319855 B CN 113319855B
Authority
CN
China
Prior art keywords
diagnosis
gravity
component
treatment
moment
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
CN202110719915.8A
Other languages
English (en)
Other versions
CN113319855A (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.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
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 Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN202110719915.8A priority Critical patent/CN113319855B/zh
Publication of CN113319855A publication Critical patent/CN113319855A/zh
Application granted granted Critical
Publication of CN113319855B publication Critical patent/CN113319855B/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
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B25HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
    • B25JMANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
    • B25J9/00Programme-controlled manipulators
    • B25J9/16Programme controls
    • B25J9/1694Programme controls characterised by use of sensors other than normal servo-feedback from position, speed or acceleration sensors, perception control, multi-sensor controlled systems, sensor fusion

Abstract

一种多关节诊疗机器人柔顺力控制模式下的重力补偿方法,属于诊疗机器人力控领域,为了解决现有的重力补偿方法对诊疗机器人重力补偿精度不足的问题。本发明实时记录机器人多个姿态下的姿态矩阵以及传感器的测量值,构造线性方程组,通过最小二乘法求解诊疗工具的参数向量;在诊疗工具处于初始位置时以及按照大地坐标系调整诊疗工具的姿态,使诊疗工具的重力只剩沿六维力与力矩传感器坐标的y轴方向时,分别记录传感器的测量值;在传感器坐标系下,利用降维解析法计算诊疗工具重心的位置;计算出诊疗工具的重力分量与力矩分量补偿值,实现对诊疗工具的重力补偿。有益效果为实现对多种诊疗工具在多姿态工作模式下的高精度重力补偿。

Description

一种多关节诊疗机器人柔顺力控制模式下的重力补偿方法
技术领域
本发明属于诊疗机器人力控领域。
背景技术
诊疗机器人在重力环境下结合六维力与力矩传感器通过人手接触控制柔顺运动时,进入控制闭环的力为诊疗机器人末端与人手之间的接触力,换言之,要实现人手力对机器人的控制,首先就要精确检测出人手力的大小和方向;实际上,安装于诊疗机器人末端的六维力传感器测量到的力,不仅包含机器人末端与人手间的接触力,还包含安装于力传感器上诊疗工具的重力;但是传感器无法自主区分人手力和诊疗工具的重力,诊疗工具的重力属于干扰力,因此有必要进行精确的重力补偿;如果不进行重力补偿,系统的控制精度必定大大受其影响;其次,根据诊疗机器人的实际需求,机器人末端需要在超声阵列探头、诊疗光纤和治疗光纤这三种诊疗工具中来回切换,而这三种诊疗工具的质量特性各异,重心也随自身姿态的变化而变化;现有的重力补偿方法无法满足诊疗机器人多模式、多姿态和高精度的工作需求,每次对末端工具进行更换或者重新安装后,都需要对其进行精密测量来测定其质量特性,操作繁琐,不适用于诊疗机器人的工作环境。
发明内容
本发明的目的是为了解决现有的重力补偿方法对诊疗机器人重力补偿精度不足的问题,提出了一种多关节诊疗机器人柔顺力控制模式下的重力补偿方法。
本发明所述的一种多关节诊疗机器人柔顺力控制模式下的重力补偿方法通过以下步骤实现的:
步骤一、将六维力与力矩传感器安装在多关节诊疗机器人末端,选取多关节诊疗机器人的多个姿态,并实时记录关节诊疗机器人多个姿态下对应的姿态矩阵以及记录关节诊疗机器人多个姿态下六维力与力矩传感器的测量值,构造线性方程组,通过最小二乘法求解诊疗工具的参数向量;
步骤二、在诊疗工具处于初始位置时,记录六维力与力矩传感器的测量值;
步骤三、按照大地坐标系调整诊疗工具的姿态,使诊疗工具的重力只剩沿六维力与力矩传感器坐标的y轴方向,并记录此时六维力与力矩传感器的测量值;
步骤四、根据步骤二记录的测量值与步骤三记录的测量值,在传感器坐标系下,利用降维解析法计算诊疗工具重心的位置;
步骤五、将步骤四计算出的诊疗工具重心的位置投影到六维力与力矩传感器坐标系上,并结合步骤一求解出的诊疗工具参数向量,计算出诊疗工具的重力分量与力矩分量补偿值,实现对诊疗工具的重力补偿。
本发明的有益效果是:该重力补偿方法能够实时提供重力补偿,来抵消多关节诊疗机器人末端诊疗工具自身重力的影响,使得补偿后的多关节诊疗机器人在任何末端位姿下不受外力作用时,多关节诊疗机器人在基坐标系下沿各方向的力和力矩都趋于零,因此该重力补偿方法重力补偿算法精度较高。
附图说明
图1为为具体实施方式所述的一种多关节诊疗机器人柔顺力控制模式下的重力补偿方法流程图;
图2为具体实施方式一中多关节诊疗机器人的实物图;
图3为具体实施方式二中基于坐标系变换的重力计算示意图;
图4为具体实施方式七中对计算出的诊疗工具重力分量与力矩分量补偿值进行验证的流程图。
具体实施方式
具体实施方式一:结合图1至图2说明本实施方式,本实施方式所述的一种多关节诊疗机器人柔顺力控制模式下的重力补偿方法通过以下步骤实现的:
步骤一、将六维力与力矩传感器安装在多关节诊疗机器人末端,选取多关节诊疗机器人的多个姿态,并实时记录关节诊疗机器人多个姿态下对应的姿态矩阵以及记录关节诊疗机器人多个姿态下六维力与力矩传感器的测量值,构造线性方程组,通过最小二乘法求解诊疗工具的参数向量;
步骤二、在诊疗工具处于初始位置时,记录六维力与力矩传感器的测量值;
步骤三、按照大地坐标系调整诊疗工具的姿态,使诊疗工具的重力只剩沿六维力与力矩传感器坐标的y轴方向,并记录此时六维力与力矩传感器的测量值;
步骤四、根据步骤二记录的测量值与步骤三记录的测量值,在传感器坐标系下,利用降维解析法计算诊疗工具重心的位置;
步骤五、将步骤四计算出的诊疗工具重心的位置投影到六维力与力矩传感器坐标系上,并结合步骤一求解出的诊疗工具参数向量,计算出诊疗工具的重力分量与力矩分量补偿值,实现对诊疗工具的重力补偿。
在本实施方式中,采用降维解析法通过记录几个特殊位置的参数,结合多关节诊疗机器人编程和多关节诊疗机器人运动学计算,在短时间内计算出诊疗工具的重力和重心相对于机器人末端中心的位置;该方法在降低计算量的同时提高了计算精度,解决了不同诊疗工具的重力大小及其在不同工况下重心位置计算的问题。
本实施方式中,重力补偿方法的基本思想是,根据多关节诊疗机器人末端诊疗工具的工具参数和机器人姿态矩阵,求解六维力与力矩传感器上的重力作用分量和力矩分量,对实际测量的力觉信息进行补偿,从而获得传感器坐标系下的实际接触力状态;因此,在补偿过程中应首先确定末端被安装工具重力的大小和重心位置;重力补偿主要由两个步骤实现:首先用最小二乘法定位诊疗机器人基坐标系、大地坐标系和工具坐标系的关系,求解诊疗工具的理论参数向量(G,OGF),其次用特殊位置法求出机器人末端诊疗工具的重力和重心位置。诊疗机器人如图2所示。
具体实施方式二:结合图3说明本实施方式,本实施方式是对具体实施方式一所述的一种多关节诊疗机器人柔顺力控制模式下的重力补偿方法进一步限定,在本实施方式中,步骤一中通过最小二乘法求解诊疗工具的参数向量的具体方法为:
设诊疗工具的重力为G,则重力向量gG为(0,0,G)T,则诊疗工具重力作用在六维力与力矩传感器上的三维力向量gF和三维力矩向量gM分别为(gFX,gFY,gFZ)T和(gMX,gMY,gMZ)T
根据如图3所示的基于坐标系变换的工具重力计算方法,三维力向量gF可表示为:
Figure GDA0003502003720000031
式中,
Figure GDA0003502003720000032
为力觉传感器坐标系OF到重力坐标系OG的姿态矩阵;多关节诊疗机器人基坐标系OB到机器人末端关节坐标系OE的姿态矩阵为
Figure GDA0003502003720000033
并且
Figure GDA0003502003720000034
由多关节诊疗机器人正运动学方程确定;
Figure GDA0003502003720000035
为力觉传感器参数姿态矩阵;OB到OG的姿态矩阵
Figure GDA0003502003720000036
用旋转角(α,β,γ)表示;
由于重力坐标系ZG轴与重力方向平行,XGOGYG平面可绕ZG轴进行任意旋转,因此设γ=0,旋转角为Q(α,β),
Figure GDA0003502003720000037
可表示为:
Figure GDA0003502003720000038
求解出诊疗工具的参数向量为(G,Q,OGF)。
具体实施方式三:本实施方式是对具体实施方式二所述的一种多关节诊疗机器人柔顺力控制模式下的重力补偿方法进一步限定,在本实施方式中,步骤二中六维力与力矩传感器的测量值分别为:fFX0、fFY0、fFZ0、fMX0、fMY0、fMZ0
其中,fFX0为:诊疗工具处于初始位置时,六维力与力矩传感器测得的沿X轴方向力的分量;fFY0为:诊疗工具处于初始位置时,六维力与力矩传感器测得的沿Y方向力的分量;fFZ0为:诊疗工具处于初始位置时,六维力与力矩传感器测得的沿Z轴方向力的分量;fMX0为:诊疗工具处于初始位置时,六维力与力矩传感器测得的沿X轴方向力矩的分量;fMY0为:诊疗工具处于初始位置时,六维力与力矩传感器测得的沿Y方向力矩的分量;fMZ0为:诊疗工具处于初始位置时,六维力与力矩传感器测得的沿Z轴方向力矩的分量。
具体实施方式四:本实施方式是对具体实施方式三所述的一种多关节诊疗机器人柔顺力控制模式下的重力补偿方法进一步限定,在本实施方式中,步骤三中六维力与力矩传感器的测量值分别为:fFX00、fFY00、fFZ00、fMX00、fMY00、fMZ00
其中,fFX00为:诊疗工具按照大地坐标系调整后,六维力与力矩传感器测得的沿X轴方向力的分量;fFY00为:诊疗工具按照大地坐标系调整后,六维力与力矩传感器测得的沿Y方向力的分量;fFZ00为:诊疗工具按照大地坐标系调整后,六维力与力矩传感器测得的沿Z轴方向力的分量;fMX00为:诊疗工具按照大地坐标系调整后,六维力与力矩传感器测得的沿X轴方向力矩的分量;fMX00为:诊疗工具按照大地坐标系调整后,六维力与力矩传感器测得的沿Y方向力矩的分量;fMZ00为:诊疗工具按照大地坐标系调整后,六维力与力矩传感器测得的沿Z轴方向力矩的分量。
具体实施方式五:本实施方式是对具体实施方式四所述的一种多关节诊疗机器人柔顺力控制模式下的重力补偿方法进一步限定,在本实施方式中,步骤四中利用降维解析法计算诊疗工具重心位置的具体方法为:
步骤四一、利用公式(3)计算出六维力与力矩传感器坐标系下的诊疗工具重心坐标OGF
所述公式(3)为:
Figure GDA0003502003720000041
其中,XGF为诊疗工具重心坐标OGF的X轴坐标,YGF为诊疗工具重心坐标OGF的Y轴坐标,ZGF为诊疗工具重心坐标OGF的Z轴坐标;
步骤四二、求解出诊疗工具重心位置为:
Figure GDA0003502003720000051
具体实施方式六:本实施方式是对具体实施方式五所述的一种多关节诊疗机器人柔顺力控制模式下的重力补偿方法进一步限定,在本实施方式中,步骤五中计算出诊疗工具的重力分量与力矩分量补偿值的具体方法为:
由于多关节诊疗机器人的姿态决定了重力在六维力与力矩传感器上的分量,因此将旋转角Q换算成坐标系旋转矩阵,将重力投影到六维力与力矩传感器坐标系上就得到重力的三个分量和力矩的三个分量:
Figure GDA0003502003720000052
其中,gFX重力在X轴上的分量,gFY为重力在Y轴上的分量,gFZ为重力在Z轴上的分量,gMX为重力力矩在X轴上的分量,gMY为重力力矩在Y轴上的分量,gMZ为重力力矩在Z轴上的分量;
从而得出多关节诊疗机器人末端诊疗工具重力分量与力矩分量的补偿值:
Figure GDA0003502003720000053
其中,fCF为诊疗工具的重力分量补偿值,fCM为诊疗工具的重力力矩分量补偿值,fFX为诊疗工具所受的力在X轴上的分量,fFY为诊疗工具所受的力在Y轴上的分量,fFZ为诊疗工具所受的力在Z轴上的分量,fMX为诊疗工具所受的力矩在X轴上的分量,fMY为诊疗工具所受的力矩在Y轴上的分量,fMZ为诊疗工具所受的力矩在Z轴上的分量。
具体实施方式七:结合图4说明本实施方式,本实施方式是对具体实施方式六所述的一种多关节诊疗机器人柔顺力控制模式下的重力补偿方法进一步限定,在本实施方式中,该重力补偿方法还包括步骤六、对步骤五计算出的诊疗工具的重力分量与力矩分量补偿值进行验证;
对步骤五计算出的诊疗工具的重力分量与力矩分量补偿值进行验证的具体方法为:
步骤六一、对六维力与力矩传感器进行校零;
步骤六二、安装诊疗工具,记录误差补偿值;
步骤六三、计算诊疗工具的重心位置,得出相对误差;
步骤六四、计算重力补偿值;
步骤六五、判断诊疗工具是否测试结束,如果测试未结束,则执行步骤六六,否则执行步骤六七;
步骤六六、更换诊疗工具的种类;
步骤六七、验证结束。
具体实施方式八:本实施方式是对具体实施方式七所述的一种多关节诊疗机器人柔顺力控制模式下的重力补偿方法进一步限定,在本实施方式中,步骤六一中对六维力与力矩传感器进行校零的具体方法为:
将六维力与力矩传感器安装在多关节诊疗机器人末端的法兰上,且六维力与力矩传感器的三个坐标轴与多关节诊疗机器人末端的三个坐标轴重合;
用示教器控制关节诊疗机器人,使六维力与力矩传感器坐标系下的Z轴平行于基座坐标系的Z轴,并且六维力与力矩传感器坐标系下的Z轴竖直向下,第一次点击六维力与力矩传感器的校零按钮;
沿六维力与力矩传感器坐标系下的Z轴负方向施加一个大小恒定的力,再次点击六维力与力矩传感器的校零按钮,完成对六维力与力矩传感器的校零。
具体实施方式九:本实施方式是对具体实施方式七所述的一种多关节诊疗机器人柔顺力控制模式下的重力补偿方法进一步限定,在本实施方式中,步骤六二中记录误差补偿值的具体方法为:
将诊疗工具固定在六维力与力矩传感器上;
用示教器控制多关节诊疗机器人,使多关节诊疗机器人各轴旋转角度都为0°,即初始位置,记录六维力与力矩传感器测得的力在Y、Z方向和力矩在X、Z方向上的分量:Fy、Fz、Tx、Tz
用示教器控制多关节诊疗机器人,使多关节诊疗机器人出各个轴转动90°,记录六维力与力矩传感器测得的力在X方向上的分量Fx和力矩在Y方向上的分量Ty
具体实施方式十:本实施方式是对具体实施方式七所述的一种多关节诊疗机器人柔顺力控制模式下的重力补偿方法进一步限定,在本实施方式中,步骤六三中计算诊疗工具的重心位置,得出相对误差的具体步骤为:
步骤六三一、用示教器移动多关节诊疗机器人到一个任意位置,用六维力与力矩传感器示值减去误差值得到的力与力矩:Fx0、Fy0、Fz0、Tx0、Ty0、Tz0
步骤六三二、用示教器按照多关节诊疗机器人基坐标系调整多关节诊疗机器人末端姿态,使诊疗工具重力只剩传感器坐标的Y轴分量,用六维力与力矩传感器示值减去误差值得到的力与力矩:Fx00、Fy00、Fz00、Tx00、Ty00、Tz00
步骤六三三、得出诊疗工具重心对于六维力与力矩传感器坐标系的位置,即测量值:
Figure GDA0003502003720000071
步骤六三四、比较理论值OGF与测量值,计算相对误差。
在本实施方式中,计算重力补偿值具体方法为:
1)在步骤2中诊疗工具固定位置的基础上,用示教器控制多关节诊疗机器人到一个任意位置,读取多关节诊疗机器人末端的诊疗工具的Q向量,在已知诊疗工具重力的情况下,计算该位置时诊疗工具重力在六维力与力矩传感器坐标系上沿X、Y、Z方向的投影Fgx、Fgy、Fgz
2)读取此时六维力与力矩传感器测得的力在X、Y、Z方向上的分量Fx、Fy、Fz,减去步骤2中的误差值,得到重力补偿的实际值,与步骤1)中的理论值进行对比,计算相对误差;
3)重复步骤1)、2)操作,得到多个不同位置时,重力补偿的理论值和实际值的相对误差。

Claims (6)

1.一种多关节诊疗机器人柔顺力控制模式下的重力补偿方法,其特征在于,该重力补偿方法通过以下步骤实现的:
步骤一、将六维力与力矩传感器安装在多关节诊疗机器人末端,选取多关节诊疗机器人的多个姿态,并实时记录关节诊疗机器人多个姿态下对应的姿态矩阵以及记录关节诊疗机器人多个姿态下六维力与力矩传感器的测量值,构造线性方程组,通过最小二乘法求解诊疗工具的参数向量;
步骤二、在诊疗工具处于初始位置时,记录六维力与力矩传感器的测量值;
步骤三、按照大地坐标系调整诊疗工具的姿态,使诊疗工具的重力只剩沿六维力与力矩传感器坐标的y轴方向,并记录此时六维力与力矩传感器的测量值;
步骤四、根据步骤二记录的测量值与步骤三记录的测量值,在传感器坐标系下,利用降维解析法计算诊疗工具重心的位置;
步骤五、将步骤四计算出的诊疗工具重心的位置投影到六维力与力矩传感器坐标系上,并结合步骤一求解出的诊疗工具参数向量,计算出诊疗工具的重力分量与力矩分量补偿值,实现对诊疗工具的重力补偿;
步骤六、对步骤五计算出的诊疗工具的重力分量与力矩分量补偿值进行验证;
对步骤五计算出的诊疗工具的重力分量与力矩分量补偿值进行验证的具体方法为:
步骤六一、对六维力与力矩传感器进行校零;
步骤六二、安装诊疗工具,记录误差补偿值;
步骤六三、计算诊疗工具的重心位置,得出相对误差;
步骤六四、计算重力补偿值;
步骤六五、判断诊疗工具是否测试结束,如果测试未结束,则执行步骤六六,否则执行步骤六七;
步骤六六、更换诊疗工具的种类;
步骤六七、验证结束;
步骤六一中对六维力与力矩传感器进行校零的具体方法为:
将六维力与力矩传感器安装在多关节诊疗机器人末端的法兰上,且六维力与力矩传感器的三个坐标轴与多关节诊疗机器人末端的三个坐标轴重合;
用示教器控制关节诊疗机器人,使六维力与力矩传感器坐标系下的Z轴平行于基座坐标系的Z轴,并且六维力与力矩传感器坐标系下的Z轴竖直向下,第一次点击六维力与力矩传感器的校零按钮;
沿六维力与力矩传感器坐标系下的Z轴负方向施加一个大小恒定的力,再次点击六维力与力矩传感器的校零按钮,完成对六维力与力矩传感器的校零;
步骤六二中记录误差补偿值的具体方法为:
将诊疗工具固定在六维力与力矩传感器上;
用示教器控制多关节诊疗机器人,使多关节诊疗机器人各轴旋转角度都为0°,即初始位置,记录六维力与力矩传感器测得的力在Y、Z方向和力矩在X、Z方向上的分量:Fy、Fz、Tx、Tz
用示教器控制多关节诊疗机器人,使多关节诊疗机器人出各个轴转动90°,记录六维力与力矩传感器测得的力在X方向上的分量Fx和力矩在Y方向上的分量Ty
步骤六三中计算诊疗工具的重心位置,得出相对误差的具体步骤为:
步骤六三一、用示教器移动多关节诊疗机器人到一个任意位置,用六维力与力矩传感器示值减去误差值得到的力与力矩:Fx0、Fy0、Fz0、Tx0、Ty0、Tz0
步骤六三二、用示教器按照多关节诊疗机器人基坐标系调整多关节诊疗机器人末端姿态,使诊疗工具重力只剩传感器坐标的Y轴分量,用六维力与力矩传感器示值减去误差值得到的力与力矩:Fx00、Fy00、Fz00、Tx00、Ty00、Tz00
步骤六三三、得出诊疗工具重心对于六维力与力矩传感器坐标系的位置,即测量值:
Figure FDA0003489914230000021
步骤六三四、比较理论值OGF与测量值,计算相对误差。
2.根据权利要求1所述的一种多关节诊疗机器人柔顺力控制模式下的重力补偿方法,其特征在于,步骤一中通过最小二乘法求解诊疗工具的参数向量的具体方法为:
设诊疗工具的重力为G,则重力向量gG为(0,0,G)T,则诊疗工具重力作用在六维力与力矩传感器上的三维力向量gF和三维力矩向量gM分别为(gFX,gFY,gFZ)T和(gMX,gMY,gMZ)T
基于坐标系变换的工具重力计算方法,三维力向量gF可表示为:
Figure FDA0003489914230000022
式中,
Figure FDA0003489914230000023
为力觉传感器坐标系OF到重力坐标系OG的姿态矩阵;多关节诊疗机器人基坐标系OB到机器人末端关节坐标系OE的姿态矩阵为
Figure FDA0003489914230000031
并且
Figure FDA0003489914230000032
由多关节诊疗机器人正运动学方程确定;
Figure FDA0003489914230000033
为力觉传感器参数姿态矩阵;OB到OG的姿态矩阵
Figure FDA0003489914230000034
用旋转角(α,β,γ)表示;
由于重力坐标系ZG轴与重力方向平行,XGOGYG平面可绕ZG轴进行任意旋转,因此设γ=0,旋转角为Q(α,β),
Figure FDA0003489914230000035
可表示为:
Figure FDA0003489914230000036
求解出诊疗工具的参数向量为(G,Q,OGF)。
3.根据权利要求2所述的一种多关节诊疗机器人柔顺力控制模式下的重力补偿方法,其特征在于,步骤二中六维力与力矩传感器的测量值分别为:fFX0、fFY0、fFZ0、fMX0、fMY0、fMZ0
其中,fFX0为:诊疗工具处于初始位置时,六维力与力矩传感器测得的沿X轴方向力的分量;fFY0为:诊疗工具处于初始位置时,六维力与力矩传感器测得的沿Y方向力的分量;fFZ0为:诊疗工具处于初始位置时,六维力与力矩传感器测得的沿Z轴方向力的分量;fMX0为:诊疗工具处于初始位置时,六维力与力矩传感器测得的沿X轴方向力矩的分量;fMY0为:诊疗工具处于初始位置时,六维力与力矩传感器测得的沿Y方向力矩的分量;fMZ0为:诊疗工具处于初始位置时,六维力与力矩传感器测得的沿Z轴方向力矩的分量。
4.根据权利要求3所述的一种多关节诊疗机器人柔顺力控制模式下的重力补偿方法,其特征在于,步骤三中六维力与力矩传感器的测量值分别为:fFX00、fFY00、fFZ00、fMX00、fMY00、fMZ00
其中,fFX00为:诊疗工具按照大地坐标系调整后,六维力与力矩传感器测得的沿X轴方向力的分量;fFY00为:诊疗工具按照大地坐标系调整后,六维力与力矩传感器测得的沿Y方向力的分量;fFZ00为:诊疗工具按照大地坐标系调整后,六维力与力矩传感器测得的沿Z轴方向力的分量;fMX00为:诊疗工具按照大地坐标系调整后,六维力与力矩传感器测得的沿X轴方向力矩的分量;fMX00为:诊疗工具按照大地坐标系调整后,六维力与力矩传感器测得的沿Y方向力矩的分量;fMZ00为:诊疗工具按照大地坐标系调整后,六维力与力矩传感器测得的沿Z轴方向力矩的分量。
5.根据权利要求4所述的一种多关节诊疗机器人柔顺力控制模式下的重力补偿方法,其特征在于,步骤四中利用降维解析法计算诊疗工具重心位置的具体方法为:
步骤四一、利用公式(3)计算出六维力与力矩传感器坐标系下的诊疗工具重心坐标OGF
所述公式(3)为:
Figure FDA0003489914230000041
其中,XGF为诊疗工具重心坐标OGF的X轴坐标,YGF为诊疗工具重心坐标OGF的Y轴坐标,ZGF为诊疗工具重心坐标OGF的Z轴坐标;
步骤四二、求解出诊疗工具重心位置为:
Figure FDA0003489914230000042
6.根据权利要求4所述的一种多关节诊疗机器人柔顺力控制模式下的重力补偿方法,其特征在于,步骤五中计算出诊疗工具的重力分量与力矩分量补偿值的具体方法为:
由于多关节诊疗机器人的姿态决定了重力在六维力与力矩传感器上的分量,因此将旋转角Q换算成坐标系旋转矩阵,将重力投影到六维力与力矩传感器坐标系上就得到重力的三个分量和力矩的三个分量:
Figure FDA0003489914230000043
其中,gFX重力在X轴上的分量,gFY为重力在Y轴上的分量,gFZ为重力在Z轴上的分量,gMX为重力力矩在X轴上的分量,gMY为重力力矩在Y轴上的分量,gMZ为重力力矩在Z轴上的分量;
从而得出多关节诊疗机器人末端诊疗工具重力分量与力矩分量的补偿值:
Figure FDA0003489914230000051
其中,fCF为诊疗工具的重力分量补偿值,fCM为诊疗工具的重力力矩分量补偿值,fFX为诊疗工具所受的力在X轴上的分量,fFY为诊疗工具所受的力在Y轴上的分量,fFZ为诊疗工具所受的力在Z轴上的分量,fMX为诊疗工具所受的力矩在X轴上的分量,fMY为诊疗工具所受的力矩在Y轴上的分量,fMZ为诊疗工具所受的力矩在Z轴上的分量。
CN202110719915.8A 2021-06-28 2021-06-28 一种多关节诊疗机器人柔顺力控制模式下的重力补偿方法 Active CN113319855B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110719915.8A CN113319855B (zh) 2021-06-28 2021-06-28 一种多关节诊疗机器人柔顺力控制模式下的重力补偿方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110719915.8A CN113319855B (zh) 2021-06-28 2021-06-28 一种多关节诊疗机器人柔顺力控制模式下的重力补偿方法

Publications (2)

Publication Number Publication Date
CN113319855A CN113319855A (zh) 2021-08-31
CN113319855B true CN113319855B (zh) 2022-04-12

Family

ID=77424973

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110719915.8A Active CN113319855B (zh) 2021-06-28 2021-06-28 一种多关节诊疗机器人柔顺力控制模式下的重力补偿方法

Country Status (1)

Country Link
CN (1) CN113319855B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115847392B (zh) * 2021-09-26 2023-10-27 苏州艾利特机器人有限公司 一种力传感器空载读数补偿方法、装置及存储介质
CN114310878B (zh) * 2021-12-22 2024-03-15 上海智能制造功能平台有限公司 基于力控的大质量零件人工引导装配方法、系统及介质

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4267027B2 (ja) * 2006-12-07 2009-05-27 ファナック株式会社 ロボット制御装置
CN103600354B (zh) * 2013-11-08 2016-10-05 北京卫星环境工程研究所 航天器机械臂柔性随动控制重力补偿方法
CN105643641B (zh) * 2014-11-11 2017-08-11 沈阳新松机器人自动化股份有限公司 力传感器标定装置、标定方法及力控制机器人
CN109129475B (zh) * 2018-08-15 2021-02-02 珠海格力电器股份有限公司 机械臂的重力补偿方法、装置、系统及存储介质
CN110103229B (zh) * 2019-06-06 2022-09-13 上海电气集团股份有限公司 设于机器人末端的工具的重力补偿方法和系统
CN110802585B (zh) * 2019-11-20 2021-04-16 山东大学 机械臂末端传感器补偿方法及接触力/力矩测量方法
CN111230879B (zh) * 2020-02-20 2023-07-11 佛山科学技术学院 一种基于力传感器的机器人末端接触力补偿方法及系统

Also Published As

Publication number Publication date
CN113319855A (zh) 2021-08-31

Similar Documents

Publication Publication Date Title
CN109822574B (zh) 一种工业机器人末端六维力传感器标定的方法
CN107042528B (zh) 一种工业机器人的运动学标定系统及方法
Huang et al. A novel modeling of volumetric errors of three-axis machine tools based on Abbe and Bryan principles
JP5618066B2 (ja) 力制御ロボットのキャリブレーション装置と方法
Zargarbashi et al. Single setup estimation of a five-axis machine tool eight link errors by programmed end point constraint and on the fly measurement with Capball sensor
CN113319855B (zh) 一种多关节诊疗机器人柔顺力控制模式下的重力补偿方法
CN109238199B (zh) 一种机器人旋转轴运动学参数标定方法
WO2018196232A1 (zh) 机器人和末端执行器的自动标定方法及系统
Gaudreault et al. Local and closed-loop calibration of an industrial serial robot using a new low-cost 3D measuring device
JP6960893B2 (ja) 工作機械の計測誤差評価方法及びプログラム
Li et al. A dynamic pre-travel error prediction model for the kinematic touch trigger probe
Botero-Valencia et al. A low-cost platform based on a robotic arm for parameters estimation of Inertial Measurement Units
EP3807058A1 (en) Estimation of payload attached to a robot arm
CN109062139B (zh) 一种基于数据驱动的机器人直线轴定位误差补偿方法
Cho et al. Screw based kinematic calibration method for robot manipulators with joint compliance using circular point analysis
Xiao et al. A new fixed axis-invariant based calibration approach to improve absolute positioning accuracy of manipulators
Majarena et al. Modelling and calibration of parallel mechanisms using linear optical sensors and a coordinate measuring machine
Santolaria et al. Self-alignment of on-board measurement sensors for robot kinematic calibration
Traslosheros et al. One camera in hand for kinematic calibration of a parallel robot
Li et al. End Position Detection of Industrial Robots Based on Laser Tracker.
Traslosheros et al. A method for kinematic calibration of a parallel robot by using one camera in hand and a spherical object
CN113878586B (zh) 机器人运动学标定装置、方法及系统
Barreto et al. Inertia parameter identification using a Stewart platform
TWI747079B (zh) 機械手臂的定位精度量測系統與方法
CN115674171A (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
GR01 Patent grant
GR01 Patent grant