CN113510705B - 一种机器人关节加速度约束规划方法和*** - Google Patents
一种机器人关节加速度约束规划方法和*** Download PDFInfo
- Publication number
- CN113510705B CN113510705B CN202110766074.6A CN202110766074A CN113510705B CN 113510705 B CN113510705 B CN 113510705B CN 202110766074 A CN202110766074 A CN 202110766074A CN 113510705 B CN113510705 B CN 113510705B
- Authority
- CN
- China
- Prior art keywords
- connecting rod
- joint
- coordinate system
- motor shaft
- acceleration
- 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
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B25—HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
- B25J—MANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
- B25J9/00—Programme-controlled manipulators
- B25J9/16—Programme controls
- B25J9/1656—Programme controls characterised by programming, planning systems for manipulators
- B25J9/1664—Programme controls characterised by programming, planning systems for manipulators characterised by motion, path, trajectory planning
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B25—HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
- B25J—MANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
- B25J17/00—Joints
- B25J17/02—Wrist joints
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B25—HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
- B25J—MANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
- B25J9/00—Programme-controlled manipulators
- B25J9/16—Programme controls
- B25J9/1679—Programme controls characterised by the tasks executed
Landscapes
- Engineering & Computer Science (AREA)
- Robotics (AREA)
- Mechanical Engineering (AREA)
- Manipulator (AREA)
- Numerical Control (AREA)
Abstract
本发明公开了一种机器人关节加速度约束规划方法和***,属于机器人运动轨迹运动学规划领域。包括:建立机器人各连杆的D‑H模型,得到D‑H连杆坐标系的齐次变换矩阵和旋转变换矩阵;分别将各个连杆的质心坐标从基坐标系变换到各自对应的D‑H连杆坐标系下,进而获取当前位姿下各个关节的电机轴上的等效转动惯量;考虑负载的影响,更新等效转动惯量;建立各个连杆的动力学方程,得到加速度与电机输出转矩的函数关系;根据等效转动惯量和所述函数关系,得到加速度的幅值约束。本发明能综合考虑机器人的连杆重力和末端负载对各个转动关节的电机轴的等效转动惯量影响,分析该姿态下关节的电机轴加速度的最值,从而提高了运动效率。
Description
技术领域
本发明属于机器人运动轨迹运动学规划领域,更具体地,涉及一种机器人关节加速度约束规划方法和***。
背景技术
随着工业机器人技术的发展,其在制造业中的应用越来越广泛。针对简单的物体搬运任务,通常要求机器人在任务中的轨迹达到加速度连续,同时满足对应的加速度幅值约束。根据机器人运动学,其末端加速度又受到各个关节的电机轴的加速度幅值约束。在进行机器人轨迹的运动学规划时,通常给定关节的电机轴的加速度幅值约束为定值,从而去规划整个轨迹的运动学。但是,在实际的任务中,各个关节的电机轴的加减速能力受到不同位姿下各个关节的电机轴上的等效转动惯量和末端负载的影响。若保守的将每个关节的电机轴的加速度幅值约束作为定值,则在不同的位姿下不能完全发挥对应的各个关节的电机轴的加减速能力,从而影响机器人的运动效率。因此,研究考虑不同位姿下机器人的各个关节的电机轴等效转动惯量对加速度约束的影响,研究不同姿态下各关节的电机轴的最大加速度,充分发挥机器人的加减速能力,对于提高机器人的运动效率很有必要。
发明内容
针对相关技术的缺陷,本发明的目的在于提供一种机器人关节加速度约束规划方法和***,旨在解决机器人在运动过程中各个关节的电机轴的加减速能力没有得到充分利用的问题。
为实现上述目的,本发明的一个方面提供了一种机器人关节加速度约束规划方法,包括以下步骤:
建立机器人各连杆的D-H模型,得到D-H连杆坐标系的齐次变换矩阵和旋转变换矩阵;
根据所述齐次变换矩阵和旋转变换矩阵,分别将各个连杆的质心坐标从基坐标系变换到各自对应的D-H连杆坐标系下,进而获取当前位姿下各个关节的电机轴上的等效转动惯量;
考虑负载的影响,更新各个关节的电机轴上的等效转动惯量;
建立各个连杆的动力学方程,得到加速度与电机输出转矩的函数关系;
根据各个关节的电机轴上的等效转动惯量和所述函数关系,得到加速度的幅值约束,从而规划机器人运动时各个关节的电机轴的加速度。
其中,为绕轴Zi旋转θi的旋转变换,为沿着Zi轴方向移动di的平移变换,为绕轴Xi旋转αi的旋转变换,为沿着Xi轴方向移动ai的平移变换,为第i个D-H连杆坐标系相对于第i-1个D-H连杆坐标系的旋转变换矩阵,为第i个D-H连杆坐标系的原点在第i-1个D-H连杆坐标系中的坐标。
进一步地,根据所述齐次变换矩阵和旋转变换矩阵,分别将各个连杆的质心坐标从基坐标系变换到各自对应的D-H连杆坐标系下包括:
其中θi为第i个连杆的D-H连杆坐标系Zi轴的关节转角,i=1,2,3,4,5,6,各个连杆在基坐标系下的质心坐标定义为0pi,各个连杆在对应D-H连杆坐标系下的质心坐标定义为ipi。
进一步地,所述当前位姿下各个关节的电机轴上的等效转动惯量为:
进一步地,考虑负载的影响,更新各个关节的电机轴上的等效转动惯量包括:
将负载的重量折算到第六连杆上,
m′6=m6+m
6x6′=(6x6·m6+6x·m)/(m6+m)
6y6′=(6y6·m6+6y·m)/(m6+m)
6z6′=(6z6·m6+6z·m)/(m6+m)
将得到的新的第六连杆质心的质量m6′和坐标(6x6′,6y6′,6z6′)代替原来第六连杆质心的质量m6和坐标(6x6,6y6,6z6),重新计算各个关节的电机轴上的等效转动惯量。
进一步地,所述加速度与电机输出转矩的函数关系为:
本发明的另一方面提供了一种机器人关节加速度约束规划***,包括:计算机可读存储介质和处理器;
所述计算机可读存储介质用于存储可执行指令;
所述处理器用于读取所述计算机可读存储介质中存储的可执行指令,执行上述的机器人关节加速度约束规划方法。
通过本发明所构思的以上技术方案,与现有技术相比,本发明针对机器人在运动过程中的加速度约束,使得机器人在不同位姿下,能综合考虑机器人的连杆重力和末端负载对各个转动关节的电机轴的等效转动惯量影响,分析该姿态下关节的电机轴加速度的最值,从而分析规划运动过程中加速度的约束,提高了机器人的运动效率。
附图说明
图1是本发明实施例加速度规划流程图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
为实现上述目的,本发明实施例提供了一种根据机器人不同位姿下各个关节电机轴的等效转动惯量变化的关节加速度约束规划方法,包括下列步骤:
(a)根据机器人各个连杆的参数,建立对应的D-H(Denavit-Hartenberg)模型。建立D-H连杆坐标系的原则为:
1.Zi轴沿关节i+1的轴向。
2.原点Oi为Zi-1和Zi轴的交点或其公垂线与关节的电机轴Zi的交点。
3.Xi轴沿Zi-1和Zi轴的公垂线方向。
4.Yi轴按照右手定则确定。
每个参数的含义:
ai:连杆长度ai定义为从Zi-1移动到Zi的距离,沿Xi轴指向为正。其实质为公垂线的长度。
αi:连杆扭角αi定义为Zi-1旋转到Zi的角度,绕Xi轴正向旋转为正。
di:关节偏移量di定义为Xi-1移动到Xi的距离,沿Zi轴指向为正。其实质为两条公垂线之间的距离。
θi:关节转角θi定义为从Xi-1旋转到Xi的角度,绕Zi轴正向旋转为正。
按照D-H模型,可以得到机器人相邻D-H连杆坐标系的齐次变换矩阵:
其中为绕轴Zi旋转θi的旋转变换,为沿着Zi轴方向移动di的平移变换,为绕轴Xi旋转αi的旋转变换,为沿着Xi轴方向移动ai的平移变换,为第i个D-H连杆坐标系相对于第i-1个D-H连杆坐标系的旋转变换矩阵,i-1pi为第i个D-H连杆坐标系的原点在第i-1个D-H连杆坐标系中的坐标。
(b)各个连杆质心坐标到各个D-H连杆坐标系下的变换。
在某一个给定位姿q=[θ1,θ2,θ3,θ4,θ5,θ6]T下(其中θ1-θ6为第一至第六连杆的D-H连杆坐标系Zi轴的关节转角),各个连杆在基坐标系下的质心坐标定义为0pi,各个连杆在对应D-H连杆坐标系下的质心坐标定义为ipi。则可以通过旋转变换矩阵和齐次变换矩阵,求出质心在对应D-H连杆坐标系下的坐标:
(c)计算在当前位姿下各个关节的电机轴上的等效转动惯量。
对于每个关节的电机轴,其绕该连杆坐标系的Zi轴旋转。在第i个D-H连杆坐标系下,第j个连杆的质心坐标为:
所以第i个关节的电机轴上的等效转动惯量Ji为:
以第六连杆为例,根据平行轴定理,由第六关节的电机轴驱动的连杆相对于第六关节的电机轴的转动惯量为:
对第五关节的电机轴,其驱动了第六连杆和第五连杆。首先计算第六连杆在第五连杆坐标系下的质心坐标:
则可以计算出连杆第五关节的电机轴Z5的等效转动惯量:
(d)考虑末端负载。
在实际的任务中,末端一般会加上负载。当末端加上负载时,为了简化,可以将负载的重量m折算到第六连杆上。
在第六个D-H连杆坐标系下考虑:
m′6=m6+m (8)
6x6′=(6x6·m6+6x·m)/(m6+m) (9)
6y6′=(6y6·m6+6y·m)/(m6+m) (10)
6z6′(6z6·m6 6z·m)/(m6+m) (11)
在将负载的重量折算到第六连杆上之后,将得到的新质心的质量和坐标(m6′,6x6′,6y6′和6z6′)代替原来第六连杆质心的质量和坐标(m6,6x6,6y6和6z6)。然后代入公式(4)中,则可以得到考虑末端负载时各个关节的电机轴上等效转动惯量。
(e)根据电机的输出转矩来计算加速度幅值约束。
对于每一个关节,其除了驱动自身的连杆外,还驱动其后一直到末端的所有连杆。第i个关节电机的输出转矩Mi可以由下式计算:
设连杆的角速度方矢量为ωi,其方向也沿着Zi轴的轴向。Mi=Mi·zi为关节电机的驱动力矩,其正负根据连杆的运动方向决定,当ωi和zi同向时取正,反之取负。Mf,i=Mf,i·zi为摩擦力矩,其正负同样和关节的运动方向相关,当ωi和zi同向时取负,反之取正。
在方程(13)两边取绝对值可得:
从上式可得加速度的值是该位姿下关节电机的输出转矩Mi的一次函数。当ωi和zi同向时,有0≤Mi≤Mmax。反之-Mmax≤Mi≤0。Mmax为关节电机输出转矩的最大值。可以根据当前关节的运动情况,将Mi的最值代入方程(14)得到对应的姿态下加速度值的最值,从而规划机器人运动时的各个关节的电机轴的加速度值的最大值约束来充分发挥其在不同姿态下的加减速能力。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种机器人关节加速度约束规划方法,其特征在于,包括以下步骤:
建立机器人各连杆的D-H模型,得到D-H连杆坐标系的齐次变换矩阵和旋转变换矩阵;
根据所述齐次变换矩阵和旋转变换矩阵,分别将各个连杆的质心坐标从基坐标系变换到各自对应的D-H连杆坐标系下,进而获取当前位姿下各个关节的电机轴上的等效转动惯量,其中,第i个关节的电机轴上的等效转动惯量为其中ixj,iyj为第j个连杆的质心在第i个D-H连杆坐标系下的横纵坐标,为第j个连杆绕过其质心且与Zi轴平行的轴旋转的转动惯量,Zi轴沿关节i+1的轴向,mj为第j个连杆的质量;
考虑负载的影响,将负载的重量折算到第六连杆上,将得到的新的第六连杆质心的质量和坐标代替原来第六连杆质心的质量和坐标,更新各个关节的电机轴上的等效转动惯量;
建立各个连杆的动力学方程,得到加速度与电机输出转矩的函数关系为其中,Mi为第i个关节电机的输出转矩,为第i个关节的电机轴的角加速度矢量且方向沿着Zi轴的轴向,Mf,i为第i个关节的摩擦转矩,为第i个关节的电机轴驱动的所有连杆的重力矩之和;
根据各个关节的电机轴上的等效转动惯量和所述函数关系,得到加速度的幅值约束,从而规划机器人运动时各个关节的电机轴的加速度。
4.如权利要求3所述的机器人关节加速度约束规划方法,其特征在于,考虑负载的影响,更新各个关节的电机轴上的等效转动惯量包括:
将负载的重量折算到第六连杆上,
m′6=m6+m
6x6′=(6x6·m6+6x·m)/(m6+m)
6y6′=(6y6·m6+6y·m)/(m6+m)
6z6′=(6z6·m6+6z·m)/(m6+m)
将得到的新的第六连杆质心的质量m6′和坐标(6x6′,6y6′,6z6′)代替原来第六连杆质心的质量m6和坐标(6x6,6y6,6z6),重新计算各个关节的电机轴上的等效转动惯量。
5.一种机器人关节加速度约束规划***,其特征在于,包括:计算机可读存储介质和处理器;
所述计算机可读存储介质用于存储可执行指令;
所述处理器用于读取所述计算机可读存储介质中存储的可执行指令,执行权利要求1至4任一项所述的机器人关节加速度约束规划方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110766074.6A CN113510705B (zh) | 2021-07-07 | 2021-07-07 | 一种机器人关节加速度约束规划方法和*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110766074.6A CN113510705B (zh) | 2021-07-07 | 2021-07-07 | 一种机器人关节加速度约束规划方法和*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113510705A CN113510705A (zh) | 2021-10-19 |
CN113510705B true CN113510705B (zh) | 2022-05-17 |
Family
ID=78066968
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110766074.6A Active CN113510705B (zh) | 2021-07-07 | 2021-07-07 | 一种机器人关节加速度约束规划方法和*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113510705B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103495977A (zh) * | 2013-09-29 | 2014-01-08 | 北京航空航天大学 | 一种6r型工业机器人负载识别方法 |
CN109514602A (zh) * | 2018-11-08 | 2019-03-26 | 东南大学 | 一种基于负载自适应辨识的工业机器人力矩补偿控制方法 |
CN110450154A (zh) * | 2019-07-18 | 2019-11-15 | 广州弘度信息科技有限公司 | 一种基于动力学的机械臂柔性控制方法 |
CN111002313A (zh) * | 2019-12-20 | 2020-04-14 | 华中科技大学 | 一种机器人模态参数辨识与动态特性分析的方法 |
CN111515956A (zh) * | 2020-05-13 | 2020-08-11 | 中科新松有限公司 | 杆件及关节柔性的机器人运动学标定方法 |
-
2021
- 2021-07-07 CN CN202110766074.6A patent/CN113510705B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103495977A (zh) * | 2013-09-29 | 2014-01-08 | 北京航空航天大学 | 一种6r型工业机器人负载识别方法 |
CN109514602A (zh) * | 2018-11-08 | 2019-03-26 | 东南大学 | 一种基于负载自适应辨识的工业机器人力矩补偿控制方法 |
CN110450154A (zh) * | 2019-07-18 | 2019-11-15 | 广州弘度信息科技有限公司 | 一种基于动力学的机械臂柔性控制方法 |
CN111002313A (zh) * | 2019-12-20 | 2020-04-14 | 华中科技大学 | 一种机器人模态参数辨识与动态特性分析的方法 |
CN111515956A (zh) * | 2020-05-13 | 2020-08-11 | 中科新松有限公司 | 杆件及关节柔性的机器人运动学标定方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113510705A (zh) | 2021-10-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110065070B (zh) | 一种基于动力学模型的机器人自适应阻抗控制*** | |
CN110666794B (zh) | 一种机器人末端工具负载力和力矩补偿的方法 | |
CN112415086B (zh) | 一种基于遥操作飞行机械臂的高空金属管道探伤*** | |
CN107169196B (zh) | 空间机器人由末端执行器向基座的动力学建模方法 | |
CN109746913B (zh) | 一种机器人姿势保持拖动示教的方法及*** | |
CN108638070A (zh) | 基于动态平衡的机器人负载重力参数辨识方法 | |
CN109015658A (zh) | 一种用于捕获翻滚目标的空间双臂机器人控制方法 | |
CN112817328B (zh) | 一种四旋翼变绳长吊挂***的路径规划与减摆控制方法 | |
CN111610721A (zh) | 模型参数完全未知的带负载四旋翼无人机速度控制方法 | |
CN108340352A (zh) | 基于示教关节臂的工业机器人远程实时操控方法 | |
CN111177941A (zh) | 一种机器人摩擦力辨识方法 | |
CN113031655B (zh) | 一种伴有负载升降的飞行吊运***及其控制方法 | |
CN113510705B (zh) | 一种机器人关节加速度约束规划方法和*** | |
Saramago et al. | An optimum robot path planning with payload constraints | |
CN110450154A (zh) | 一种基于动力学的机械臂柔性控制方法 | |
CN109366486A (zh) | 柔性机器人逆运动学求解方法、***、设备、存储介质 | |
Hu et al. | Control and simulation of cable-driven parallel robots in offshore cargo handling | |
Bulut et al. | Computed torque control of an aerial manipulation system with a quadrotor and a 2-dof robotic arm | |
Li et al. | Neural network method for robot arm of service robot based on DH model | |
CN110815230B (zh) | 一种基于虚拟采样的六轴码垛机器人位姿适配方法 | |
Lintao et al. | Kinematic analysis and speed control of 3SPS-1S parallel mechanism for end actuator of segment erector | |
Liu et al. | Rigid-Flexible Coupling Simulation and Vibration Analysis of Flexible Robot | |
bin Hasnan et al. | A hardware-in-the-loop simulation and test for unmanned ground vehicle on indoor environment | |
Gardecki et al. | Multirotor aerial platform with manipulation system-static disturbances | |
CN116079743B (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 |