CN109807899B - 针对拖动示教的协作机器人摩擦力矩补偿方法 - Google Patents

针对拖动示教的协作机器人摩擦力矩补偿方法 Download PDF

Info

Publication number
CN109807899B
CN109807899B CN201910176132.2A CN201910176132A CN109807899B CN 109807899 B CN109807899 B CN 109807899B CN 201910176132 A CN201910176132 A CN 201910176132A CN 109807899 B CN109807899 B CN 109807899B
Authority
CN
China
Prior art keywords
joint
friction torque
speed
robot
test
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
CN201910176132.2A
Other languages
English (en)
Other versions
CN109807899A (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.)
Hefei Hagong Renji Intelligent Technology Co., Ltd
Original Assignee
HRG International Institute for Research and Innovation
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 HRG International Institute for Research and Innovation filed Critical HRG International Institute for Research and Innovation
Priority to CN201910176132.2A priority Critical patent/CN109807899B/zh
Publication of CN109807899A publication Critical patent/CN109807899A/zh
Application granted granted Critical
Publication of CN109807899B publication Critical patent/CN109807899B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Manipulator (AREA)

Abstract

本发明公开了一种针对拖动示教的协作机器人摩擦力矩补偿方法,涉及协作机器人技术领域,本发明建立起机器人连杆重力矩和机器人关节摩擦力矩之间的数学联系,通过机器人连杆重力矩对应的关节电枢电流估算机器人关节摩擦力矩对应的关节电枢电流。本发明的优点在于:提出了一种简单可靠的摩擦力矩估计方法,跟系统辨识方法相比,此方法原理简单,成本低廉,使用这种方法,机器人关节摩擦力矩的补偿变得不再艰难,拖动外力比无摩擦力矩补偿时有了显著减小。

Description

针对拖动示教的协作机器人摩擦力矩补偿方法
技术领域
本发明涉及协作机器人技术领域,尤其涉及一种针对拖动示教的协作机器人摩擦力矩补偿方法。
背景技术
随着工业机器人技术的日益成熟和应用的不断普及,对人机协作要求越来越高。传统的示教方式主要依赖于示教器,工作效率较低,对操作人员知识水平要求较高。而采用直接力矩控制的拖动示教技术可以帮助操作人员将机器人末端操作器直接拖动到任意期望位置,大大降低了对操作人员的要求,提高了工作效率。
拖动示教技术要求机器人各连杆处于或接近零外力静止悬浮状态。此时,机械臂各关节位置几乎保持恒定,故而机器人各连杆的惯性力矩、科氏和离心力矩都接近于零,机械臂各关节的输出转矩接近完美地抵消了机械臂各连杆的自身重力矩和关节摩擦力矩。当机器人处于零外力静止悬浮状态时,很小的外作用力就可以将机器人本体拖动到工作空间内的任意位姿。
机器人各连杆的重力矩跟各关节位置相关,它们可以经由简单的计算求得。然而,摩擦力矩由于受到机器人关节位姿,关节速度、负载大小和温度等很多因素的影响,很难通过实验直接测量。如何测量或估计机器人关节摩擦力矩已经成为了协作机器人研发过程中的瓶颈问题,不解决这个问题,协作机器人性能就不可能有质的提高。
理论上,协作机器人关节摩擦力矩可以用系统辨识的方法间接求得,但这种方法理论艰深,步骤繁琐,成本很高,需要海量的实验数据,并且需要实验人员具有极强的对实验数据进行综合分析的能力,不便于工程应用。
发明内容
本发明所要解决的技术问题是:提供一种简单实用的协作机器人摩擦力矩补偿方法。
本发明是通过以下技术方案解决上述技术问题的:针对拖动示教的协作机器人摩擦力矩补偿方法,包括以下步骤:
步骤一,针对多轴协作机器人,依次对需要进行摩擦力矩补偿的关节进行单独测试;
在关节位置空间内等间隔选取多个测试角度,在关节速度空间内等间隔选取多个测试速度;
步骤二,从最小的测试速度开始测试,每个测试速度关节在位置空间内往复转动一定次数,每次往复转动间隔一定时间;
依照上述方法遍历关节速度空间内的所有测试速度,记录下关节位置、速度和电流数据;
步骤三,选取每个测试速度的中间若干次往复运动的数据做处理,在每个测试角度,分别将顺时针转动和逆时针转动时的关节电枢电流数据取平均,得到每个测试角度顺时针转动和逆时针转动时的关节电枢电流值I-和I+
步骤四,关节在匀速转动时所需要的电枢电流为
Figure GDA0002608317030000021
其中Ig为补偿重力矩所需要的电枢电流,
Figure GDA0002608317030000022
为补偿摩擦力矩所需要的电枢电流,If为正数;
以逆时针方向为正,逆时针转动时所需要的电枢电流为I+=Ig+If,顺时针转动时所需要的电枢电流为I-=Ig-If
假设顺时针转动和逆时针转动时的摩擦力矩相同,则在所测试的关节角度下,If=(I+-I-)/2,Ig=(I++I-)/2;
步骤五,在所测试的关节角度下,将机器人各连杆相对于各关节的重力矩G(q)建模并计算出结果;
假设某一个关节电机转矩常数为Kt,该关节的减速比为n,那么重力矩对应的关节电枢电流Ig=G(q)/(Kt·n);
步骤六,在所选取的测试速度下,以Ig为自变量,If为因变量,用二次曲线进行拟合,得到两者之间的关系
Figure GDA0002608317030000031
由于在不同测试速度下拟合得到的If和Ig二次曲线系数不同,故对不同测试速度下的系数进行拟合,得到
Figure GDA0002608317030000032
作为优化的技术方案,步骤二中,每个测试速度关节在位置空间内往复转动10次,每次往复转动时间间隔为5秒。
作为优化的技术方案,步骤二中,转动关节时,关节先加速运动,达到测试速度,然后以测试速度匀速运动,最后减速到0。
作为优化的技术方案,步骤三中,选取每个测试速度的中间6次往复运动的数据做处理。
作为优化的技术方案,还包括以下步骤:步骤七,在步骤四中,假设了顺时针转动和逆时针转动时的摩擦力矩相同,为此,将摩擦力矩对应的电枢电流写为λIf,其中λ为调节系数,通过对实际机器人的零力控制进行测试,调节系数λ。
作为优化的技术方案,还包括以下步骤:步骤八,将符号函数取值如下:当
Figure GDA0002608317030000033
时,
Figure GDA0002608317030000034
Figure GDA0002608317030000035
时,
Figure GDA0002608317030000036
否则,
Figure GDA0002608317030000037
这里σ为正实数,具体取值根据速度信号噪声大小决定。
本发明的优点在于:本发明针对协作机器人的拖动示教功能,提出了一种简单可靠的摩擦力矩估计方法,跟系统辨识方法相比,此方法原理简单,操作步骤简单直接,不需要大量的实验数据,不易出错;不需要昂贵的实验设备,成本低廉;使用这种方法,机器人关节摩擦力矩的补偿变得不再艰难,拖动外力比无摩擦力矩补偿时有了显著减小;大量的实验结果表明,此方法效果显著,非常适合工程应用。
附图说明
图1是本发明针对拖动示教的协作机器人摩擦力矩补偿方法应用于实施例中的摩擦力矩补偿电流If和重力矩补偿电流Ig之间的关系函数图。
具体实施方式
针对拖动示教的协作机器人摩擦力矩补偿方法,包括以下步骤:
步骤一,针对多轴协作机器人,依次对需要进行摩擦力矩补偿的关节进行单独测试。
在关节位置空间内等间隔选取多个测试角度,在关节速度空间内等间隔选取多个测试速度,关节速度空间的范围根据实际拖动示教速度范围选取。
步骤二,从最小的测试速度开始测试,每个测试速度关节在位置空间内往复转动10次,每次往复转动时间间隔为5秒。
转动关节时,关节先加速运动,达到测试速度,然后以测试速度匀速运动,最后减速到0。
关节运动经过所有测试角度,且在每个测试角度都处于匀速运动段。
依照上述方法遍历关节速度空间内的所有测试速度,记录下关节位置、速度和电流数据。
步骤三,选取每个测试速度的中间6次往复运动的数据做处理,在每个测试角度,分别将顺时针转动和逆时针转动时的关节电枢电流数据取平均,得到每个测试角度顺时针转动和逆时针转动时的关节电枢电流值I-和I+
步骤四,关节在匀速转动时所需要的电枢电流为
Figure GDA0002608317030000041
其中Ig为补偿重力矩所需要的电枢电流,
Figure GDA0002608317030000042
为补偿摩擦力矩所需要的电枢电流,If为正数。
以逆时针方向为正,逆时针转动时所需要的电枢电流为I+=Ig+If,顺时针转动时所需要的电枢电流为I-=Ig-If
假设顺时针转动和逆时针转动时的摩擦力矩相同,则在所测试的关节角度下,If=(I+-I-)/2,Ig=(I++I-)/2。
在实际应用中,我们只利用I+和I-求得If,而不用它们来估算Ig,因为Ig可以通过更为系统化的方法求得。
步骤五,在所测试的关节角度下,将机器人各连杆相对于各关节的重力矩G(q)建模并计算出结果。
假设某一个关节电机转矩常数为Kt,该关节的减速比为n,那么重力矩对应的关节电枢电流Ig=G(q)/(Kt·n)。
步骤六,在所选取的测试速度下,以Ig为自变量,If为因变量,用二次曲线进行拟合,得到两者之间的关系
Figure GDA0002608317030000051
由于在不同测试速度下拟合得到的If和Ig二次曲线系数不同,故对不同测试速度下的系数进行拟合,得到
Figure GDA0002608317030000052
注意重力矩对应的电枢电流Ig只与关节位置q有关,而与关节速度
Figure GDA0002608317030000053
无关,关节摩擦力矩对应的电枢电流与关节位置q和关节速度
Figure GDA0002608317030000054
都有关系。
步骤七,在步骤四中,假设了顺时针转动和逆时针转动时的摩擦力矩相同,为此,将摩擦力矩对应的电枢电流写为λIf,其中λ为调节系数,通过对实际机器人的零力控制进行测试,小幅度调节系数λ,使得零力控制达到较为理想的状态。
步骤八,由于实际的速度信号有一定的噪声,如果用传统符号函数的定义,可能会出现在速度较小时,符号函数的值在正负之间反复跳变,导致摩擦力矩补偿不准确,故这里将符号函数取值如下:当
Figure GDA0002608317030000055
时,
Figure GDA0002608317030000056
Figure GDA0002608317030000057
时,
Figure GDA0002608317030000058
否则,
Figure GDA0002608317030000059
这里σ为正实数,具体取值根据速度信号噪声大小决定。
按上述步骤,在自行研制的六轴协作机器人上进行了多次实验验证。
这里给出了特定转速下,在自行研制的六轴协作机器人的第二个关节上的测试结果,如图1所示,摩擦力矩补偿电流If和重力矩补偿电流Ig之间的关系为
Figure GDA00026083170300000510
本发明的关键点是建立起机器人连杆重力矩和机器人关节摩擦力矩之间的数学联系,通过机器人连杆重力矩对应的关节电枢电流估算机器人关节摩擦力矩对应的关节电枢电流。为了建立起这个联系,分析了大量的机器人连杆重力矩数据和机器人关节摩擦力矩数据,考虑一定范围内速度变化对摩擦力矩的影响,并测试了不同的拟合曲线,最终得到了十分满意的结果。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (6)

1.一种针对拖动示教的协作机器人摩擦力矩补偿方法,其特征在于,包括以下步骤:
步骤一,针对多轴协作机器人,依次对需要进行摩擦力矩补偿的关节进行单独测试;
在关节位置空间内等间隔选取多个测试角度,在关节速度空间内等间隔选取多个测试速度;
步骤二,从最小的测试速度开始测试,每个测试速度关节在位置空间内往复转动一定次数,每次往复转动间隔一定时间;
依照上述方法遍历关节速度空间内的所有测试速度,记录下关节位置、速度和电流数据;
步骤三,选取每个测试速度的中间若干次往复运动的数据做处理,在每个测试角度,分别将顺时针转动和逆时针转动时的关节电枢电流数据取平均,得到每个测试角度顺时针转动和逆时针转动时的关节电枢电流值I_和I+
步骤四,关节在匀速转动时所需要的电枢电流为
Figure FDA0002608317020000011
其中Ig为补偿重力矩所需要的电枢电流,
Figure FDA0002608317020000012
为补偿摩擦力矩所需要的电枢电流,If为正数;
以逆时针方向为正,逆时针转动时所需要的电枢电流为I+=Ig+If,顺时针转动时所需要的电枢电流为I-=Ig-If
假设顺时针转动和逆时针转动时的摩擦力矩相同,则在所测试的关节角度下,If=(I+-I-)/2,Ig=(I++I-)/2;
步骤五,在所测试的关节角度下,将机器人各连杆相对于各关节的重力矩G(q)建模并计算出结果;
假设某一个关节电机转矩常数为Kt,该关节的减速比为n,那么重力矩对应的关节电枢电流Ig=G(q)/(Kt·n);
步骤六,在所选取的测试速度下,以Ig为自变量,If为因变量,用二次曲线进行拟合,得到两者之间的关系
Figure FDA0002608317020000021
由于在不同测试速度下拟合得到的If和Ig二次曲线系数不同,故对不同测试速度下的系数进行拟合,得到
Figure FDA0002608317020000022
2.如权利要求1所述的针对拖动示教的协作机器人摩擦力矩补偿方法,其特征在于:步骤二中,每个测试速度关节在位置空间内往复转动10次,每次往复转动时间间隔为5秒。
3.如权利要求1所述的针对拖动示教的协作机器人摩擦力矩补偿方法,其特征在于:步骤二中,转动关节时,关节先加速运动,达到测试速度,然后以测试速度匀速运动,最后减速到0。
4.如权利要求1所述的针对拖动示教的协作机器人摩擦力矩补偿方法,其特征在于:步骤三中,选取每个测试速度的中间6次往复运动的数据做处理。
5.如权利要求1所述的针对拖动示教的协作机器人摩擦力矩补偿方法,其特征在于,还包括以下步骤:步骤七,在步骤四中,假设了顺时针转动和逆时针转动时的摩擦力矩相同,为此,将摩擦力矩对应的电枢电流写为λIf,其中λ为调节系数,通过对实际机器人的零力控制进行测试,调节系数λ。
6.如权利要求1所述的针对拖动示教的协作机器人摩擦力矩补偿方法,其特征在于,还包括以下步骤:步骤八,将符号函数取值如下:当
Figure FDA0002608317020000023
时,
Figure FDA0002608317020000024
Figure FDA0002608317020000025
时,
Figure FDA0002608317020000026
否则,
Figure FDA0002608317020000027
这里σ为正实数,具体取值根据速度信号噪声大小决定。
CN201910176132.2A 2019-03-08 2019-03-08 针对拖动示教的协作机器人摩擦力矩补偿方法 Active CN109807899B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910176132.2A CN109807899B (zh) 2019-03-08 2019-03-08 针对拖动示教的协作机器人摩擦力矩补偿方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910176132.2A CN109807899B (zh) 2019-03-08 2019-03-08 针对拖动示教的协作机器人摩擦力矩补偿方法

Publications (2)

Publication Number Publication Date
CN109807899A CN109807899A (zh) 2019-05-28
CN109807899B true CN109807899B (zh) 2020-09-29

Family

ID=66608527

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910176132.2A Active CN109807899B (zh) 2019-03-08 2019-03-08 针对拖动示教的协作机器人摩擦力矩补偿方法

Country Status (1)

Country Link
CN (1) CN109807899B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2021037582A (ja) * 2019-09-03 2021-03-11 ファナック株式会社 アームを有するロボットの制御装置
CN113442111B (zh) * 2020-03-27 2022-07-12 北京配天技术有限公司 机器人及其力矩补偿方法
CN115674190B (zh) * 2022-09-30 2024-05-07 深圳市越疆科技股份有限公司 协作机械臂及其运动控制方法、碰撞检测方法、控制系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106994686A (zh) * 2016-12-01 2017-08-01 遨博(北京)智能科技有限公司 关节外力力矩的计算方法及装置、机器人
CN107097233A (zh) * 2017-04-17 2017-08-29 珞石(山东)智能科技有限公司 一种无力矩传感器的工业机器人拖动示教方法
CN108189036A (zh) * 2018-01-17 2018-06-22 广州视源电子科技股份有限公司 力矩控制方法、装置、机器人及存储介质
CN108582069A (zh) * 2018-04-17 2018-09-28 上海达野智能科技有限公司 机器人拖动示教系统和方法、存储介质、操作系统
CN108789363A (zh) * 2018-05-25 2018-11-13 雅客智慧(北京)科技有限公司 一种基于力传感器的直接拖动示教系统及方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6582483B2 (ja) * 2015-03-26 2019-10-02 セイコーエプソン株式会社 ロボット制御装置およびロボットシステム

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106994686A (zh) * 2016-12-01 2017-08-01 遨博(北京)智能科技有限公司 关节外力力矩的计算方法及装置、机器人
CN107097233A (zh) * 2017-04-17 2017-08-29 珞石(山东)智能科技有限公司 一种无力矩传感器的工业机器人拖动示教方法
CN108189036A (zh) * 2018-01-17 2018-06-22 广州视源电子科技股份有限公司 力矩控制方法、装置、机器人及存储介质
CN108582069A (zh) * 2018-04-17 2018-09-28 上海达野智能科技有限公司 机器人拖动示教系统和方法、存储介质、操作系统
CN108789363A (zh) * 2018-05-25 2018-11-13 雅客智慧(北京)科技有限公司 一种基于力传感器的直接拖动示教系统及方法

Also Published As

Publication number Publication date
CN109807899A (zh) 2019-05-28

Similar Documents

Publication Publication Date Title
CN109807899B (zh) 针对拖动示教的协作机器人摩擦力矩补偿方法
CN109927032B (zh) 一种基于高阶滑模观测器的机械臂轨迹跟踪控制方法
CN106994686B (zh) 关节外力力矩的计算方法及装置、机器人
Wang et al. Uncalibrated visual tracking control without visual velocity
CN112241124B (zh) 一种自适应反演积分非奇异快速终端滑模控制器设计方法
Wu et al. Parametric identification and structure searching for underwater vehicle model using symbolic regression
Sauvée et al. Image based visual servoing through nonlinear model predictive control
CN110394801B (zh) 一种机器人的关节控制系统
CN107702709B (zh) 一种非合作目标运动与惯性参数的时频域混合辨识方法
CN110941184A (zh) 一种电磁轴承柔性转子不同位系统滑模振动主动控制方法
Leite et al. Passivity‐based adaptive 3D visual servoing without depth and image velocity measurements for uncertain robot manipulators
CN111624878A (zh) 自主式水面机器人轨迹跟踪的积分滑模获取方法及系统
CN110908389B (zh) 一种针对不确定水下机器人的自适应快速速度跟踪控制方法
Sun Kinematics model identification and motion control of robot based on fast learning neural network
Jia et al. Improved dynamic parameter identification method relying on proprioception for manipulators
Mustafa et al. Adaptive backstepping sliding mode control based on nonlinear disturbance observer for trajectory tracking of robotic manipulator
CN108445778B (zh) 一种针对空间非合作目标非完全约束组合体的动力学建模方法
Liu Parameter identification for LuGre friction model using genetic algorithms
CN109783931B (zh) 一种基于高斯过程回归的永磁球形电机电磁转矩建模方法
CN109048995B (zh) 一种三自由度Delta并联机器人的非线性关节摩擦力补偿方法
CN108646563B (zh) 一种基于均值耦合的多机械臂系统固定时间参数辨识与位置同步控制方法
CN112925206B (zh) 一种非线性多倒立摆互联系统分布式鲁棒故障诊断方法
CN113954077A (zh) 带有能量优化的水下游动机械臂轨迹跟踪控制方法及装置
WO2021040706A1 (en) Qualitative mechanics based system modeling
CN114454161B (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
TR01 Transfer of patent right

Effective date of registration: 20200930

Address after: 230000 Room 501, building 3, intelligent equipment science and Technology Park, 3963 Susong Road, Hefei Economic and Technological Development Zone, Anhui Province

Patentee after: Hefei Hagong Renji Intelligent Technology Co., Ltd

Address before: 236000 Anhui city of Hefei Province Economic and Technological Development Zone Cuiwei Road No. 6 Haiheng building room 6012

Patentee before: HRG INTERNATIONAL INSTITUTE FOR RESEARCH & INNOVATION

TR01 Transfer of patent right
CP02 Change in the address of a patent holder

Address after: 230000 room 510, building 3, intelligent equipment science and Technology Park, 3963 Susong Road, Hefei Economic and Technological Development Zone, Anhui Province

Patentee after: Hefei Hagong Renji Intelligent Technology Co., Ltd

Address before: 230000 Room 501, building 3, intelligent equipment science and Technology Park, 3963 Susong Road, Hefei Economic and Technological Development Zone, Anhui Province

Patentee before: Hefei Hagong Renji Intelligent Technology Co., Ltd

CP02 Change in the address of a patent holder