CN106020125B - 一种基于速度连续的工业机器人运动学算法 - Google Patents

一种基于速度连续的工业机器人运动学算法 Download PDF

Info

Publication number
CN106020125B
CN106020125B CN201610460099.2A CN201610460099A CN106020125B CN 106020125 B CN106020125 B CN 106020125B CN 201610460099 A CN201610460099 A CN 201610460099A CN 106020125 B CN106020125 B CN 106020125B
Authority
CN
China
Prior art keywords
formula
industrial robot
speed
matrix
error
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
CN201610460099.2A
Other languages
English (en)
Other versions
CN106020125A (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.)
Eft Intelligent Equipment Ltd By Share Ltd
Original Assignee
Eft Intelligent Equipment Ltd By Share 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 Eft Intelligent Equipment Ltd By Share Ltd filed Critical Eft Intelligent Equipment Ltd By Share Ltd
Priority to CN201610460099.2A priority Critical patent/CN106020125B/zh
Publication of CN106020125A publication Critical patent/CN106020125A/zh
Application granted granted Critical
Publication of CN106020125B publication Critical patent/CN106020125B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B19/00Programme-control systems
    • G05B19/02Programme-control systems electric
    • G05B19/18Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form
    • G05B19/414Structure of the control system, e.g. common controller or multiprocessor systems, interface to servo, programmable interface controller

Landscapes

  • Engineering & Computer Science (AREA)
  • Human Computer Interaction (AREA)
  • Manufacturing & Machinery (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Manipulator (AREA)
  • Numerical Control (AREA)

Abstract

本发明涉及一种基于速度连续的工业机器人运动学算法,包括以下两个部分:(1)对工业机器人轴进行速度和位置的规划;(2)误差控制。本发明以6自由度工业机器人为例,依据D‑H参数,计算机器人的雅克比矩阵,在笛卡尔空间进行TCP速度Ve的规划,反雅克比运算,再反雅克比矩阵中增加反奇异位的变量α,求出插补时间的关节速度然后对各轴速度进行积分出各轴的位置量q,运算过程严密,保证了运算结果的准确度。

Description

一种基于速度连续的工业机器人运动学算法
技术领域
本发明涉及工业机器人控制方法技术领域,具体的说是一种基于速度连续的工业机器人运动学算法。
背景技术
随着工业自动化的发展,工业机器人的使用领域越来越大,现场工艺对机器人运动要求越来越严格,工业机器人不仅仅是关节运动,针对工业机器人的结构不同,笛卡尔坐标系下运动存在奇异位置,有些工艺要求必须速度连续地通过奇异位,因此人们需要一种速度连续的工业机器人运动学算法,帮助工业机器人迅速通过笛卡尔坐标系下的奇异位置。
发明内容
针对上述技术的缺陷,本发明提出一种基于速度连续的工业机器人运动学算法。
一种基于速度连续的工业机器人运动学算法,包括以下两个部分:
(1)对工业机器人轴进行速度和位置的规划,具体算法为:
①依据工业机器人D-H矩阵参数求出工业机器人的雅克比矩阵J(q),引入反奇异位变量α,得出反雅克比矩阵J*与反奇异位变量α的关系,见公式(1):
J*=JT(JJT2I)-1 (1)
式中:I为单位阵;JT为J的旋转矩阵。
②按照公式(2)和公式(3)计算出工业机器人的关节速度
式中:Ve为机器人笛卡尔坐标系下的速度量。
③按照公式(4)计算出关节位置q:
式中:Tk+1和Tk为时间周期。
(2)误差控制:反奇异位变量α的引入,会同时引入位姿的误差,即控制位姿与实际位姿的差值(ep,eo),其中:误差(ep,eo)由公式(5)(6)(7)定义,为了减小误差,需要对误差进行抑制;
ep=pd-pe(q) (5)
式中:ep为位置误差;pd为命令位置,pe为实际位置。
eo=γsinθ (6)
式中:eo为姿态误差,γ为旋转轴,θ为旋转角度。
R(θ,γ)=RdRe T(q) (7)
式中:R为姿态矩阵,Re为实际姿态矩阵,Rd为命令姿态矩阵。
具体抑制方法为:建立一个控制位姿与实际位姿的差值(ep,eo)、增益K、反雅克比矩阵J*之间的闭环控制模型,通过调整参数变量α与增益K的值,满足机器人奇异位的速度连续,同时保证轨迹精度。
本发明的有益效果是:
本发明以6自由度工业机器人为例,依据D-H参数,计算机器人的雅克比矩阵,在笛卡尔空间进行TCP速度Ve的规划,反雅克比运算,再反雅克比矩阵中增加反奇异位的变量α,求出插补时间的关节速度然后对各轴速度进行积分出各轴的位置量q,运算过程严密,通过调整增益,可以控制路径的轨迹误差精度在0.005mm以内。
附图说明
下面结合附图和实施例对本发明进一步说明。
图1是本发明的差值(ep,eo)、增益以及反雅克比矩阵J*之间的闭环控制模型图。
具体实施方式
为了使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解,下面对本发明进一步阐述。
如图1所示,一种基于速度连续的工业机器人运动学算法,一种基于速度连续的工业机器人运动学算法,包括以下两个部分:
(3)对工业机器人轴进行速度和位置的规划,具体算法为:
①依据工业机器人D-H矩阵参数求出工业机器人的雅克比矩阵J(q),引入反奇异位变量α,得出反雅克比矩阵J*与反奇异位变量α的关系,见公式(1):
J*=JT(JJT2I)-1 (1)
②结合公式(1),按照公式(2)和公式(3)计算出工业机器人的关节速度
③结合公式(3),按照公式(4)计算出关节位置q:
(4)误差控制:令控制位姿与实际位姿的差值为(ep,eo),其中:误差(ep,eo)由公式(5)(6)(7)定义;
ep=pd-pe(q)(5)
eo=γsinθ (6)
R(θ,γ)=RdRe T(q) (7)
对误差进行抑制的具体方法为:建立一个控制位姿与实际位姿的差值(ep,eo)、增益K、增益的比例系数K(·)、反雅克比矩阵J*之间的闭环控制模型,如图1,通过调整参数变量α与增益K的值,满足机器人奇异位的速度连续,同时保证轨迹精度。
以上显示和描述了本发明的基本原理、主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明内。本发明要求保护范围由所附的权利要求书及其等效物界定。

Claims (1)

1.一种基于速度连续的工业机器人运动学算法,其特征在于:包括以下两个部分:
(1)对工业机器人轴进行速度和位置的规划,具体算法为:
①依据工业机器人D-H矩阵参数求出工业机器人的雅克比矩阵J(q),引入反奇异位变量α,得出反雅克比矩阵J*与反奇异位变量α的关系,见公式(1):
J*=JT(JJT2I)-1 (1)
式中:I为单位阵;JT为J的旋转矩阵;
②按照公式(2)和公式(3)计算出工业机器人的关节速度
式中:Ve为机器人笛卡尔坐标系下的速度量;
③按照公式(4)计算出关节位置q:
式中:Tk+1和Tk为时间周期;
(2)误差控制:反奇异位变量α的引入,会同时引入位姿的误差,即控制位姿与实际位姿的差值(ep,eo),其中:误差(ep,eo)由公式(5)(6)(7)定义,为了减小误差,需要对误差进行抑制;
ep=pd-pe(q) (5)
式中:ep为位置误差;pd为命令位置,pe为实际位置;
eo=γsinθ (6)
式中:eo为姿态误差,γ为旋转轴,θ为旋转角度;
R(θ,γ)=RdRe T(q) (7)
式中:R为姿态矩阵,Re为实际姿态矩阵,Rd为命令姿态矩阵;
具体抑制方法为:建立一个控制位姿与实际位姿的差值(ep,eo)、增益K、反雅克比矩阵J*之间的闭环控制模型,通过调整参数变量α与增益K的值,满足机器人奇异位的速度连续,同时保证轨迹精度。
CN201610460099.2A 2016-06-23 2016-06-23 一种基于速度连续的工业机器人运动学算法 Active CN106020125B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610460099.2A CN106020125B (zh) 2016-06-23 2016-06-23 一种基于速度连续的工业机器人运动学算法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610460099.2A CN106020125B (zh) 2016-06-23 2016-06-23 一种基于速度连续的工业机器人运动学算法

Publications (2)

Publication Number Publication Date
CN106020125A CN106020125A (zh) 2016-10-12
CN106020125B true CN106020125B (zh) 2018-10-02

Family

ID=57086436

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610460099.2A Active CN106020125B (zh) 2016-06-23 2016-06-23 一种基于速度连续的工业机器人运动学算法

Country Status (1)

Country Link
CN (1) CN106020125B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109375584B (zh) * 2018-09-28 2020-08-07 湖南希法工程机械有限公司 臂架关节控制方法及系统、工程机械及机器可读存储介质
CN111590592A (zh) * 2020-06-05 2020-08-28 中科新松有限公司 机械臂的控制方法及装置、电子设备、存储介质
CN111736514B (zh) * 2020-06-10 2020-12-04 杭州凯尔达机器人科技股份有限公司 基于通用计算机的机器人控制系统

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103909522A (zh) * 2014-03-19 2014-07-09 华南理工大学 一种六自由度工业机器人通过奇异域的方法
CN104035083A (zh) * 2014-06-20 2014-09-10 电子科技大学 一种基于量测转换的雷达目标跟踪方法
CN104483967A (zh) * 2014-11-11 2015-04-01 浙江师范大学 一种基于节能考虑的轮式移动机器人轨迹跟踪控制方法
US9008363B1 (en) * 2013-01-02 2015-04-14 Google Inc. System and method for computing optical flow
CN104655043A (zh) * 2014-10-28 2015-05-27 电子科技大学 一种基于机器视觉火车导轨轮廓测量中轮廓配准方法
CN105082134A (zh) * 2015-08-06 2015-11-25 珞石(北京)科技有限公司 一种基于多判据的六自由度串联机器人奇异性处理方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9008363B1 (en) * 2013-01-02 2015-04-14 Google Inc. System and method for computing optical flow
CN103909522A (zh) * 2014-03-19 2014-07-09 华南理工大学 一种六自由度工业机器人通过奇异域的方法
CN104035083A (zh) * 2014-06-20 2014-09-10 电子科技大学 一种基于量测转换的雷达目标跟踪方法
CN104655043A (zh) * 2014-10-28 2015-05-27 电子科技大学 一种基于机器视觉火车导轨轮廓测量中轮廓配准方法
CN104483967A (zh) * 2014-11-11 2015-04-01 浙江师范大学 一种基于节能考虑的轮式移动机器人轨迹跟踪控制方法
CN105082134A (zh) * 2015-08-06 2015-11-25 珞石(北京)科技有限公司 一种基于多判据的六自由度串联机器人奇异性处理方法

Also Published As

Publication number Publication date
CN106020125A (zh) 2016-10-12

Similar Documents

Publication Publication Date Title
CN109159151B (zh) 一种机械臂空间轨迹跟踪动态补偿方法和系统
CN104932506B (zh) 基于快速终端滑模的轮式移动机器人轨迹跟踪方法
CN105353725B (zh) 用于工业机器人的过辅助点姿态空间圆弧插补方法
CN106020125B (zh) 一种基于速度连续的工业机器人运动学算法
CN109857100B (zh) 一种基于反演法和快速终端滑模的复合轨迹跟踪控制算法
CN104589304A (zh) 机器人控制装置以及机器人
CN107589720B (zh) 一种等效平面交叉耦合控制方法
CN108621162A (zh) 一种机械臂运动规划方法
CN109782759B (zh) 一种轮式移动机器人的近似解耦、快速轨迹跟踪控制方法
CN102528554B (zh) 一种五轴加工奇异区域的轨迹优化方法
CN106965171A (zh) 具备学习功能的机器人装置
CN107206588A (zh) 机械手的偏移自动调整装置及机械手的偏移自动调整方法
CN105929791B (zh) 平面直角坐标运动系统的直接轮廓控制方法
JP5790840B2 (ja) ロボットの制御装置及びロボットの姿勢補間方法
CN107186753B (zh) 工业机器人性能测试的工作空间确定方法
CN108406765A (zh) 一种开链式多臂机器人阻抗控制方法
CN106097390B (zh) 一种基于卡尔曼滤波的机器人运动学参数标定方法
CN105563482A (zh) 用于工业机器人的末端执行器的旋转运动规划方法
CN107145640B (zh) 中性浮力实验中漂浮基座与机械臂的动态标度规划方法
CN206224153U (zh) 一种适用于不平衡转矩伺服系统的控制装置
CN208384428U (zh) 基于反向位移补偿的稳定转台系统
CN109866224A (zh) 一种机器人雅可比矩阵计算方法、装置及存储介质
CN112650217A (zh) 一种基于评估函数的机器人轨迹跟踪策略动态优化方法
CN107807516A (zh) 一种基于模糊自适应pid控制的模拟转台控制系统
CN107398903A (zh) 工业机械手臂执行端的轨迹控制方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant