CN105573143A - 用于六自由度的工业机器人的逆运动学求解方法 - Google Patents
用于六自由度的工业机器人的逆运动学求解方法 Download PDFInfo
- Publication number
- CN105573143A CN105573143A CN201510848596.5A CN201510848596A CN105573143A CN 105573143 A CN105573143 A CN 105573143A CN 201510848596 A CN201510848596 A CN 201510848596A CN 105573143 A CN105573143 A CN 105573143A
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- joint
- industrial robot
- transformation matrix
- joints
- 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
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
- G05B17/00—Systems involving the use of models or simulators of said systems
- G05B17/02—Systems involving the use of models or simulators of said systems electric
-
- 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
- G05B2219/00—Program-control systems
- G05B2219/30—Nc systems
- G05B2219/39—Robotics, robotics to robotics hand
- G05B2219/39417—6-DOF
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Manipulator (AREA)
- Numerical Control (AREA)
Abstract
本发明提出了一种用于六自由度的工业机器人的逆运动学求解方法,包括:建立连杆坐标系,其中,所述连杆坐标系包括:X-Z轴平面,所述工业机器人的六个关节对应坐标系和所述工业机器人末端执行器的坐标系;在所述连杆坐标系基础上,根据所述工业机器人的六个关节之间的结构几何参数计算欧拉角变换矩阵;利用所述欧拉叫变换矩阵,按照预设顺序依次求解所述工业机器人的六个关节的旋转轴转动角度;根据所述六个关节的旋转轴转动角度对应的八个求解结果,解离上一关节空间所处位置对应的各个关节旋转轴角度差值范数总和,依据范数综合最小原则选择最优解。本发明具有求解精度高、求解速度快、求解过程更加简单易懂,计算量小的特点。
Description
技术领域
本发明涉及无人机技术领域,特别涉及一种用于六自由度的工业机器人的逆运动学求解方法。
背景技术
工业机器人的逆向运动学问题是在给定执行器末端相对于基座的位置和姿态,以及所有连杆几何参数的情况下,求取所有关节转动角度值,是正向运动学的逆过程。正向运动学根据几何参数可以得到前后相邻关节之间的欧拉角变换矩阵T,当知道每个关节转动角度θ,通过变换矩阵T简单运算即可得到末端执行器的位置和姿态,且解唯一;而逆向运动学的求解则相对复杂,并且可能出现多解或无解的情况。
现在市场上大部分六自由度工业机器人都是属于Pieper提出的一类腕部相邻的三个关节旋转轴相交于同一点(俗称腕点)的工业机器人,现有的通用坐标系建模都是采用DH参数法,同时现有逆运动学求解方法大多数都是利用代数法或几何分析法单独实现的。
现有工业机器人的逆运动学求解主要存在缺陷及不足:针对特殊构型机器人,通用求解方法坐标系建模方法单一、求解过程复杂难懂、计算量大且求解速度慢。
发明内容
本发明的目的旨在至少解决所述技术缺陷之一。
为此,本发明的目的在于提出一种用于六自由度的工业机器人的逆运动学求解方法,具有求解精度高、求解速度快、求解过程更加简单易懂,计算量小的特点。
为了实现上述目的,本发明的实施例提供一种用于六自由度的工业机器人的逆运动学求解方法,包括如下步骤:
步骤S1,建立连杆坐标系,其中,所述连杆坐标系包括:X-Z轴平面,所述工业机器人的六个关节对应坐标系和所述工业机器人末端执行器的坐标系;
步骤S2,在所述连杆坐标系基础上,根据所述工业机器人的六个关节之间的结构几何参数计算欧拉角变换矩阵;
步骤S3,利用所述欧拉叫变换矩阵,按照预设顺序依次求解所述工业机器人的六个关节的旋转轴转动角度;
步骤S4,根据所述六个关节的旋转轴转动角度对应的八个求解结果,解离上一关节空间所处位置对应的各个关节旋转轴角度差值范数总和,依据范数综合最小原则选择最优解。
进一步,在所述步骤S2中,
根据所述工业机器人的六个关节之间的结构几何参数d1计算姿态变换矩阵和位置变换矩阵其中i=0,1,2,…,9;
计算坐标系{j}相对于坐标系{i}的变换矩阵
进一步,所述工业机器人的六个关节间的变换矩阵满足如下:
其中,为末端执行器的坐标系相对于基础坐标系的变换矩阵,为第一关节的坐标系相对于基础坐标系的变换矩阵,为第二关节的坐标系相对于第一关节坐标系的变换矩阵;为第三关节的坐标系相对于第二关节坐标系的变换矩阵;为第四关节的坐标系相对于第三关节坐标系的变换矩阵;为第五关节的坐标系相对于第四关节坐标系的变换矩阵;为第六关节的坐标系相对于第五关节坐标系的变换矩阵。
进一步,在所述步骤S3中,依次计算所述工业机器人的第一关节、第五关节、第六关节、第三关节、第二关节和第四关节的旋转轴转动角度。
根据本发明实施例的用于六自由度的工业机器人的逆运动学求解方法,针对该类特殊构型工业机器人,坐标系建模方法简单易懂。并且利用几何分析法和欧拉角变换法联合求解机器人逆运动学解,能够保证该类工业机器人逆运动学解的求解精度。本发明提出的逆运动学求解过程相比于单一代数法和几何分析法要简单许多,能够保证该类工业机器人逆运动学求解速度,求解速度快。本发明中提出的逆运动学求解过程更加简单易懂,计算量小。
本发明附加的方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解,其中:
图1为根据本发明实施例的用于六自由度的工业机器人的逆运动学求解方法的流程图;
图2为根据本发明实施例的六自由度特殊构型工业机器人的二维平面图;
图3为根据本发明实施例的坐标系{5}的原点在{0}坐标系的X-Y平面的投影示意图;
图4为根据本发明实施例的都坐标系{1}的X-Y平面示意图;
图5为根据本发明实施例的坐标系{1}的Y-Z平面示意图。
具体实施方式
下面详细描述本发明的实施例,所述实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。
本发明提出一种用于六自由度的工业机器人的逆运动学求解方法,利用几何分析法和欧拉角变换法结合求解该类特殊构型工业机器人逆运动学解的方法。该类工业机器人特殊构型在于坐标系{3},{4},{5}的旋转轴不相交于同一点,这种机器人腕部不存在奇异性。
如图1所示,本发明实施例的用于六自由度的工业机器人的逆运动学求解方法,包括如下步骤:
步骤S1,建立连杆坐标系
如图2所示,连杆坐标系包括:X-Z轴平面,Y轴根据右手法则确定,工业机器人的六个关节对应坐标系{0},{1},{2},{3},{4},{5},其中{0}为固定基础坐标系,工业机器人末端执行器的坐标系{6}。图2中带箭头射线围绕的坐标轴为相应关节旋转轴,依次设该类特殊构型工业机器人的六个关节旋转角度分别为变量θ1,θ2,θ3,θ4,θ5,θ6。
步骤S2,在连杆坐标系基础上,根据工业机器人的六个关节之间的结构几何参数计算欧拉角变换矩阵。
根据工业机器人的六个关节之间的结构几何参数d1计算姿态变换矩阵和位置变换矩阵其中i=0,1,2,…,9;
计算坐标系{j}相对于坐标系{i}的变换矩阵
工业机器人的六个关节间的变换矩阵满足如下:
其中,为末端执行器的坐标系相对于基础坐标系的变换矩阵,为第一关节的坐标系相对于基础坐标系的变换矩阵,为第二关节的坐标系相对于第一关节坐标系的变换矩阵;为第三关节的坐标系相对于第二关节坐标系的变换矩阵;为第四关节的坐标系相对于第三关节坐标系的变换矩阵;为第五关节的坐标系相对于第四关节坐标系的变换矩阵;为第六关节的坐标系相对于第五关节坐标系的变换矩阵。
步骤S3,利用欧拉叫变换矩阵,按照预设顺序依次求解工业机器人的六个关节的旋转轴转动角度。
在本发明的一个实施例中,依次计算工业机器人的第一关节、第五关节、第六关节、第三关节、第二关节和第四关节的旋转轴转动角度。
第一步,分析该特殊构型机器人几何结构及投影关系,求解逆运动学解的θ1。
设坐标系{0}的原点到坐标系{5}原点矢量分离出 其中为3x3的旋转矩阵,得到等式
根据工业机器人逆运动学解法,及几何尺寸均已知。
参考图3,计算得到
根据几何对称特性,θ1有两个解,分别对应工业机器人领域中的“左肩”,“右肩”。
第二步,求解逆运动学解的θ5。
参考图4,以坐标系{1}为参考基础坐标系。坐标系{4}旋转θ5,坐标系{5}的中d7长度在坐标系{4}的X轴方向的投影长度为d7cos(θ5),然后得到关系式,
联立求解公式(3),得到
θ5有两个解,分别对应工业机器人领域中的“腕上”,“腕下”。
第三步,求解逆运动学解的θ6。
求解出关节转动角θ1,θ5后,可以求出再根据坐标系{1},{2},{3}的X轴方向始终平行的特点,可以得到下面两个关系式:
cos(θ1)r12+sin(θ1)r22=-sin(θ5)cos(θ6),(5)
cos(θ1)r13+sin(θ1)r23=sin(θ5)sin(θ6),(6)
联立上述两个关系式(5)和(6),推导得到即得到
但当sin(θ5)=0或者(1)=0,(2)=0时,θ6是没有明确定义的。
第四步,求解逆运动学解的θ3及θ2。
结合以上已经求解出的关节转动角度θ1,θ5和θ6,根据欧拉角变换法,可以得到
参考图5,OA=d2,AB=d4,
得到等式
然后计算得到
θ3有两个解,分别对应工业机器人领域中的“上肘”,“下肘”。、
根据正弦定理,可以得到以下关系式:参考图5得到根据关节二转动角度θ2=η-β关系式,求解得到
第五步,求解逆运动学解的θ4。
根据前几步已经求解出各关节转动角度,求出变换矩阵
步骤S4,根据六个关节的旋转轴转动角度对应的八个求解结果,解离上一关节空间所处位置对应的各个关节旋转轴角度差值范数总和,依据范数综合最小原则选择最优解。
从步骤S3中求出的八组可能的逆运动学解中选择最优解,根据八组可能的解离起始(上一关节空间所处位置,即六个转动轴旋转角度θ1,θ2,θ3,θ4,θ5,θ6)对应各关节旋转轴角度差值的范数总和最小值原则进行选择。
根据本发明实施例的用于六自由度的工业机器人的逆运动学求解方法,针对该类特殊构型工业机器人,坐标系建模方法简单易懂。并且利用几何分析法和欧拉角变换法联合求解机器人逆运动学解,能够保证该类工业机器人逆运动学解的求解精度。本发明提出的逆运动学求解过程相比于单一代数法和几何分析法要简单许多,能够保证该类工业机器人逆运动学求解速度,求解速度快。本发明中提出的逆运动学求解过程更加简单易懂,计算量小。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。本发明的范围由所附权利要求极其等同限定。
Claims (4)
1.一种用于六自由度的工业机器人的逆运动学求解方法,其特征在于,包括如下步骤:
步骤S1,建立连杆坐标系,其中,所述连杆坐标系包括:X-Z轴平面,所述工业机器人的六个关节对应坐标系和所述工业机器人末端执行器的坐标系;
步骤S2,在所述连杆坐标系基础上,根据所述工业机器人的六个关节之间的结构几何参数计算欧拉角变换矩阵;
步骤S3,利用所述欧拉叫变换矩阵,按照预设顺序依次求解所述工业机器人的六个关节的旋转轴转动角度;
步骤S4,根据所述六个关节的旋转轴转动角度对应的八个求解结果,解离上一关节空间所处位置对应的各个关节旋转轴角度差值范数总和,依据范数综合最小原则选择最优解。
2.如权利要求1所述的用于六自由度的工业机器人的逆运动学求解方法,其特征在于,在所述步骤S2中,
根据所述工业机器人的六个关节之间的结构几何参数di计算姿态变换矩阵和位置变换矩阵其中i=0,1,2,...,9;
计算坐标系{j}相对于坐标系{i}的变换矩阵
3.如权利要求2所述的用于六自由度的工业机器人的逆运动学求解方法,其特征在于,所述工业机器人的六个关节间的变换矩阵满足如下:
其中,为末端执行器的坐标系相对于基础坐标系的变换矩阵,为第一关节的坐标系相对于基础坐标系的变换矩阵,为第二关节的坐标系相对于第一关节坐标系的变换矩阵;为第三关节的坐标系相对于第二关节坐标系的变换矩阵;为第四关节的坐标系相对于第三关节坐标系的变换矩阵;为第五关节的坐标系相对于第四关节坐标系的变换矩阵;为第六关节的坐标系相对于第五关节坐标系的变换矩阵。
4.如权利要求1所述的用于六自由度的工业机器人的逆运动学求解方法,其特征在于,在所述步骤S3中,依次计算所述工业机器人的第一关节、第五关节、第六关节、第三关节、第二关节和第四关节的旋转轴转动角度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510848596.5A CN105573143B (zh) | 2015-11-30 | 2015-11-30 | 用于六自由度的工业机器人的逆运动学求解方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510848596.5A CN105573143B (zh) | 2015-11-30 | 2015-11-30 | 用于六自由度的工业机器人的逆运动学求解方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105573143A true CN105573143A (zh) | 2016-05-11 |
CN105573143B CN105573143B (zh) | 2019-05-31 |
Family
ID=55883407
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510848596.5A Active CN105573143B (zh) | 2015-11-30 | 2015-11-30 | 用于六自由度的工业机器人的逆运动学求解方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105573143B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106844951A (zh) * | 2017-01-18 | 2017-06-13 | 哈尔滨工业大学深圳研究生院 | 基于分段几何法求解超冗余机器人逆运动学的方法及系统 |
CN106845037A (zh) * | 2017-03-21 | 2017-06-13 | 山东科技大学 | 一种五自由度串联机器人的逆运动学通用求解方法 |
CN107016207A (zh) * | 2017-04-17 | 2017-08-04 | 珞石(北京)科技有限公司 | 基于特定关节运动的工业机器人负载动力学参数辨识方法 |
CN107685330A (zh) * | 2017-10-18 | 2018-02-13 | 佛山华数机器人有限公司 | 一种六自由度手腕偏置串联机器人的运动学逆解求解方法 |
CN107791248A (zh) * | 2017-09-28 | 2018-03-13 | 浙江理工大学 | 基于不满足pipper准则的六自由度串联机器人的控制方法 |
CN109693235A (zh) * | 2017-10-23 | 2019-04-30 | 中国科学院沈阳自动化研究所 | 一种仿人眼视觉跟踪装置及其控制方法 |
CN113334393A (zh) * | 2021-08-06 | 2021-09-03 | 成都博恩思医学机器人有限公司 | 一种机械臂控制方法、系统、机器人及存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102509025A (zh) * | 2011-11-25 | 2012-06-20 | 苏州大学 | 一种六自由度仿人灵巧臂逆运动学的快速求解方法 |
CN102637158A (zh) * | 2012-04-28 | 2012-08-15 | 谷菲 | 一种六自由度串联机器人运动学逆解的求解方法 |
JP2014014920A (ja) * | 2012-06-15 | 2014-01-30 | Kumamoto Univ | 機構制御装置、プログラム及び機構制御方法 |
CN103968761A (zh) * | 2014-05-28 | 2014-08-06 | 中科华赫(北京)科技有限责任公司 | 串联关节式机器人绝对定位误差校准方法及标定系统 |
-
2015
- 2015-11-30 CN CN201510848596.5A patent/CN105573143B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102509025A (zh) * | 2011-11-25 | 2012-06-20 | 苏州大学 | 一种六自由度仿人灵巧臂逆运动学的快速求解方法 |
CN102637158A (zh) * | 2012-04-28 | 2012-08-15 | 谷菲 | 一种六自由度串联机器人运动学逆解的求解方法 |
JP2014014920A (ja) * | 2012-06-15 | 2014-01-30 | Kumamoto Univ | 機構制御装置、プログラム及び機構制御方法 |
CN103968761A (zh) * | 2014-05-28 | 2014-08-06 | 中科华赫(北京)科技有限责任公司 | 串联关节式机器人绝对定位误差校准方法及标定系统 |
Non-Patent Citations (2)
Title |
---|
姜宏超等: "六自由度模块化机械臂的逆运动学分析", 《浙江大学学报(工学版)》 * |
朱齐丹 等: "六自由度机械臂逆运动学算法", 《机器人技术与应用》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106844951A (zh) * | 2017-01-18 | 2017-06-13 | 哈尔滨工业大学深圳研究生院 | 基于分段几何法求解超冗余机器人逆运动学的方法及系统 |
CN106844951B (zh) * | 2017-01-18 | 2020-08-18 | 哈尔滨工业大学深圳研究生院 | 基于分段几何法求解超冗余机器人逆运动学的方法及系统 |
CN106845037A (zh) * | 2017-03-21 | 2017-06-13 | 山东科技大学 | 一种五自由度串联机器人的逆运动学通用求解方法 |
CN106845037B (zh) * | 2017-03-21 | 2018-07-24 | 山东科技大学 | 一种五自由度串联机器人的逆运动学通用求解方法 |
CN107016207A (zh) * | 2017-04-17 | 2017-08-04 | 珞石(北京)科技有限公司 | 基于特定关节运动的工业机器人负载动力学参数辨识方法 |
CN107791248A (zh) * | 2017-09-28 | 2018-03-13 | 浙江理工大学 | 基于不满足pipper准则的六自由度串联机器人的控制方法 |
CN107685330A (zh) * | 2017-10-18 | 2018-02-13 | 佛山华数机器人有限公司 | 一种六自由度手腕偏置串联机器人的运动学逆解求解方法 |
CN109693235A (zh) * | 2017-10-23 | 2019-04-30 | 中国科学院沈阳自动化研究所 | 一种仿人眼视觉跟踪装置及其控制方法 |
CN113334393A (zh) * | 2021-08-06 | 2021-09-03 | 成都博恩思医学机器人有限公司 | 一种机械臂控制方法、系统、机器人及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN105573143B (zh) | 2019-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105573143A (zh) | 用于六自由度的工业机器人的逆运动学求解方法 | |
Zein et al. | Non-singular assembly-mode changing motions for 3-RPR parallel manipulators | |
CN111267073B (zh) | 一种基于增强现实技术的工业机器人示教系统及方法 | |
Xu et al. | Dual arm-angle parameterisation and its applications for analytical inverse kinematics of redundant manipulators | |
CN104965517A (zh) | 一种机器人笛卡尔空间轨迹的规划方法 | |
CN107791248B (zh) | 基于不满足Pieper准则的六自由度串联机器人的控制方法 | |
Hyatt et al. | Developable mechanisms on right conical surfaces | |
CN111496783B (zh) | 一种6r工业机器人逆运动学求解方法 | |
Bouzgou et al. | Workspace analysis and geometric modeling of 6 dof fanuc 200ic robot | |
Ren et al. | Determination of the base position and working area for mobile manipulators | |
Sreenath et al. | Monocular tracking of human hand on a smart phone camera using mediapipe and its application in robotics | |
CN105404174B (zh) | 一种六自由度串联机器人运动学反解的求解方法 | |
Zhang et al. | Kuka youBot arm shortest path planning based on geodesics | |
JP6891773B2 (ja) | ロボットの制御装置及びロボットの逆変換処理方法 | |
JP2017131990A (ja) | 干渉回避方法 | |
Simas et al. | Extended Jacobian for Redundant Robots Obtained From the Kinematics Constraints | |
CN115481489A (zh) | 基于增强现实的白车身与生产线适配性验证系统及方法 | |
JP6897520B2 (ja) | ロボットの制御装置及びロボットの逆変換処理方法 | |
JP6891774B2 (ja) | ロボットの制御装置及びロボットの逆変換処理方法 | |
Wang et al. | Workspace Calculation of Free-Floating Planar Dual-Arm Space Robot | |
Qin et al. | Where to install the manipulator: optimal installation pose planning based on whale algorithm | |
Vogt et al. | Automatic end tool alignment through plane detection with a RANSAC-algorithm for robotic grasping | |
Xia et al. | ABB-IRB120 Robot Modeling and Simulation Based on MATLAB | |
Han et al. | Path planning for robotic manipulator in narrow space with search algorithm | |
Wang et al. | Efficient Collision Detection for Industrial Robot Simulation System |
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 | ||
TR01 | Transfer of patent right |
Effective date of registration: 20200527 Address after: No.321, Xinghu street, Suzhou Industrial Park, Suzhou, Jiangsu Province Patentee after: SUZHOU KONUO ROBOT Co.,Ltd. Address before: 100085 Beijing city Haidian District xi'erqi Lingxiu new Silicon Valley D District No. 57 building four unit 502 Patentee before: ROKAE, Inc. |
|
TR01 | Transfer of patent right |