CN106020125B - 一种基于速度连续的工业机器人运动学算法 - Google Patents
一种基于速度连续的工业机器人运动学算法 Download PDFInfo
- 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
Links
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B19/00—Programme-control systems
- G05B19/02—Programme-control systems electric
- G05B19/18—Numerical 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/414—Structure 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(JJT+α2I)-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(JJT+α2I)-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(JJT+α2I)-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的值,满足机器人奇异位的速度连续,同时保证轨迹精度。
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)
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)
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 | 珞石(北京)科技有限公司 | 一种基于多判据的六自由度串联机器人奇异性处理方法 |
-
2016
- 2016-06-23 CN CN201610460099.2A patent/CN106020125B/zh active Active
Patent Citations (6)
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 |