CN106113034B - 一种六自由度机械臂考虑力约束的轨迹规划方法 - Google Patents

一种六自由度机械臂考虑力约束的轨迹规划方法 Download PDF

Info

Publication number
CN106113034B
CN106113034B CN201610421404.7A CN201610421404A CN106113034B CN 106113034 B CN106113034 B CN 106113034B CN 201610421404 A CN201610421404 A CN 201610421404A CN 106113034 B CN106113034 B CN 106113034B
Authority
CN
China
Prior art keywords
contact
mechanical arm
force
transition
constraint
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
CN201610421404.7A
Other languages
English (en)
Other versions
CN106113034A (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.)
Harbin Engineering University
Original Assignee
Harbin Engineering University
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 Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN201610421404.7A priority Critical patent/CN106113034B/zh
Publication of CN106113034A publication Critical patent/CN106113034A/zh
Application granted granted Critical
Publication of CN106113034B publication Critical patent/CN106113034B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B25HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
    • B25JMANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
    • B25J9/00Programme-controlled manipulators
    • B25J9/16Programme controls
    • B25J9/1656Programme controls characterised by programming, planning systems for manipulators
    • B25J9/1664Programme controls characterised by programming, planning systems for manipulators characterised by motion, path, trajectory planning
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00Program-control systems
    • G05B2219/30Nc systems
    • G05B2219/40Robotics, robotics mapping to robotics vision
    • G05B2219/40519Motion, trajectory planning

Landscapes

  • Engineering & Computer Science (AREA)
  • Robotics (AREA)
  • Mechanical Engineering (AREA)
  • Manipulator (AREA)

Abstract

本发明属于机械臂考虑力约束的轨迹规划领域,具体涉及一种六自由度机械臂考虑力约束的轨迹规划方法。本发明包括:(1)将力约束转化为接触形变和接触运动速度的约束;(2)优化动态接触冲击的机械臂运动速度;(3)基于力约束进行轨迹规划。本发明将机械臂力约束转化为接触形变和接触运动速度的约束,再根据接触运动学,研究考虑动态接触冲击的机械臂运动速度优化方法和考虑接触形变的基于不同接触边缘的连续轨迹规划方法,此方法为考虑力约束的机械臂轨迹规划提供了新方法。此方法避免了机械臂末端与外界环境的接触碰撞过程在极短时间内完成,产生极大瞬间冲击,导致接触力超出安全范围。

Description

一种六自由度机械臂考虑力约束的轨迹规划方法
技术领域
本发明属于机械臂迹规划领域,具体涉及一种六自由度机械臂考虑力约束的轨迹规划方法。
背景技术
机械臂轨迹规划是指在规定时间内,机械臂按照一定的速度和加速度完成从初始位置姿态运动到期望位置姿态的运动。机械臂考虑力约束的轨迹规划是使得机械臂在期望力的约束下,让六自由度机械臂实现期望轨迹的规划。基于力约束的机械臂轨迹规划有很多方法,如4-3-4轨迹,3-5-3轨迹以及5段3次轨迹等,这些方法均可进行处理机械臂动力约束下的轨迹规划,但面向关节空间的笛卡尔路径会损失部分采样点。
本发明提出的方法针对一类典型的六自由度机械臂,工业机械臂大部分属于此类机械臂,因此对于这类机械臂考虑力约束的轨迹规划的研究具有很大的意义。
发明内容
本发明的目的在于提供一种考虑接触形变基于不同接触边缘的六自由度机械臂考虑力约束的轨迹规划方法。
本发明的目的是这样实现的:
1.一种六自由度机械臂考虑力约束的轨迹规划方法,其特征在于,包括如下步骤:
(1)将力约束转化为接触形变和接触运动速度的约束;
根据接触模型将机械臂末端与环境的接触分为静态接触和动态接触;
(1.1)静态接触过程中,将力的约束转化为沿法线方向的接触位移的约束,即
其中,Ks是接触刚度;XE是机械臂末端沿接触面法向的形变量;F表示机械末端力传感器测量的接触面法线方向上的接触力;
(1.2)动态接触也就是冲击,机械臂末端与外界环境的接触碰撞过程在极短的时间内完成,将力约束转化为接触速度的约束,来实现满足机械臂轨迹规划过程中冲击最优;
(2)优化动态接触冲击的机械臂运动速度;
(2.1)采用S型加减速曲线使加速度连续;
采用S型速度变化使得运动过程中加速度连续,减小各关节的冲击,使机械臂运动和作业满足力的约束;
优化后的S型速度曲线共分为七个过程:
(2.1.1)加加速度一定、加速度匀速增加过程,加速度从初始值加速到最大加速度;
(2.1.2)匀加加速度加速过程,保持最大加速度不变,速度匀速增加;
(2.1.3)减加速度加速过程,速度继续增加,加速度以恒定的加加速度减小直到速度加到最大速度值;
(2.1.4)匀速运动过程,保持最大速度运动,加速度为零;
(2.1.5)减减速度减速过程,速度减小,减速度以恒定的减减速度减小;
(2.1.6)匀减减速度减速过程,保持减速度不变,速度匀速减小;
(2.1.7)减减速度一定减速度匀速减小过程;
(2.2)归一化时间算子;
笛卡尔空间位置插补的表达式为:
P(t)=Ps+l(t)(Pe-Ps),
其中,P(t)表示插补点的坐标;l(t)为归一化时间算子;Ps为期望位置;Pe为初始位置;
要对时间算子进行平滑处理,使得速度曲线平滑,定义归一化时间算子的表达式为:
其中,i=0,1,2,3,4,5,6;l(0)=0;l(T)=1;l(t)在区间内单调递增;S(t)为各时间段内对应的总位移;
(3)基于力约束进行轨迹规划;
(3.1)基于规则接触边缘具有力约束的连续轨迹规划:
规则接触边缘具有力约束的连续轨迹采用圆弧过渡的方法代替直线-直线连续轨迹规划方法进行平滑处理,达到运动在接触临界接触面时不具有沿接触面法线方向的位移、速度和加速度分量,过渡圆弧的联结点分别选取在前一段直线的减速段和后一段直线的加速段;
(3.1.1)在接触模型有效的情况下,得到机械臂末端与外界环境的接触刚度为K,在已知约束力为F的情况下,可求得接触形变Xe
Xe=F/K;
(3.1.2)通过接触形变,得到临界接触点B点的空间位置坐标p2
p2=pq+Xe
轨迹终点M的空间位置坐标p3为:
p3=p2+L;
(3.1.3)过渡精度的确定,过渡精度决定过渡圆弧的大小,决定实际轨迹与期望轨迹逼近程度;设过渡力大小为F',则过渡精度re为:
re=F-F'
当力传感器示数时,进行空间直线轨迹运动;当力传感器示数时,为过渡圆弧的起点;当力传感器示数时,进行圆弧过渡运运动;当力传感器示数时,为过渡圆弧的终点,之后机械臂末端与外界环境保持力F进行沿接触面切平面方向的运动;
(3.1.4)确定过渡圆弧起始联结点pe1和过渡圆弧终止联结点pe2
当力传感器示数到达过渡力F'时开始圆弧过渡运动,当力传感器示数到达约束力F时结束圆弧过渡运动;即过渡圆弧联结点为接触力为F'和F时机械臂末端的位置点;
当力传感器示数时,通过位置编码器可知当前机械臂六个关节关节角θe11~θe16,通过机械臂正运动学式得到此刻对应的机械臂末端实际位置即为过渡圆弧起始联结点pe1;当力传感器示数时,通过位置编码器可知当前机械臂六个关节关节角θe21~θe26,通过机械臂正运动学式得到此刻对应的机械臂末端实际位置即为过渡圆弧终止联结点pe2
(3.1.5)轨迹生成器通过路径点和过渡圆弧起始联结点pe1和终止联结点pe2的空间位置坐标,采用具有S型速度曲线拟合的空间直线和圆弧规划方法生成运动轨迹;
(3.1.6)运动过程中用力传感器进行监测,看整个运动过程是否满足期望力F的约束,如果时,重新规划。
(3.2)基于曲率接触边缘具有力约束的连续轨迹规划:
(3.2.1)机械臂实现沿接触面法线方向直线运动和沿临界安全接触面的接触运动和沿接触面法线方向远离临界安全接触面的直线运动,是一个直线-圆弧-直线连续的轨迹规划过程,确定运动过程中在安全接触圆柱面上的具有力约束运动中的特殊点M、N点;
(3.2.2)采用圆弧过渡来解决直线-圆弧连续运动时产生的尖角问题;使得末端在沿安全接触面的运动过程中只有沿接触面切线方向位移和速度,法线方向速度为零;
(3.2.3)通过调节过渡精度来调整实际轨迹与期望轨迹的精度,满足接触模型的基于圆柱接触面的连续轨迹规划方法与基于规则接触面连续轨迹规划方法中过渡圆弧联结点求取方法一致,过渡完后机械臂末端是在一定接触力下沿着接触面的圆弧轨迹。
本发明的有益效果在于:
本发明将机械臂力约束转化为接触形变和接触运动速度的约束,再根据接触运动学,提出了动态接触冲击的机械臂运动速度优化方法和考虑接触形变的基于不同接触边缘的连续轨迹规划方法,此方法为考虑力约束的机械臂轨迹规划提供了新方法。此方法避免了机械臂末端与外界环境的接触碰撞过程在极短时间内完成,产生极大瞬间冲击,导致接触力超出安全范围。
附图说明
图1是机械臂具有力约束轨迹规划系统结构图;
图2是机械臂末端与规则接触面示意图;
图3是基于规则接触面具有力约束轨迹规划示意图;
图4是机械臂末端与具有曲率接触面示意图。
具体实施方式
下面结合附图对本发明做进一步描述。
本发明提供的是一类六自由度机械臂考虑力约束的轨迹规划方法,属于机械臂轨迹规划领域。传统六自由度机械臂轨迹规划时,只考虑机械臂末端的位置而对力约束的轨迹规划研究很少,但此类问题尤为重要,针对这一现状,本发明提出一类六自由度机械臂考虑力约束的轨迹规划方法。本发明提出的机械臂考虑力约束的轨迹方法如下:首先利用接触模型和接触运动学将力约束转化为接触形变和接触运动速度的约束,然后在优化动态接触冲击的轨迹规划中,对运动速度进行S型曲线处理,最后提出对基于规则接触边缘和基于曲率接触边缘的具有力约束的连续轨迹规划方法,一种在力的过渡阈值范围内,采用法线方向接触运动向切线方向接触运动过渡的且具有S型速度曲线的连续轨迹规划方法。此方法避免机械臂与外界环境的接触过程中产生过大冲击造成超出约束要求,且适用于具有力约束的六自由度机械臂的轨迹规划,在不同接触边缘提出不同轨迹规划,更大程度的满足现代工程项目对工业机械臂轨迹规划的要求。
一类六自由度机械臂考虑力约束的轨迹规划方法,包括如下步骤:
1、将力约束转化为接触形变和接触运动速度的约束
接触力不仅与接触形变有关,还与接触过程中的冲击有关,于是接触模型可将机械臂末端与环境的接触分为静态接触和动态接触。
⑴静态接触作用时间较长,接触状态变化缓慢匀速或保持不变,接触过程中,可将力的约束转化为沿法线方向的接触位移的约束,即
其中,Ks是接触刚度,与材料有关;XE是机械臂末端沿接触面法向的形变量;F表示机械末端力传感器测量的接触面法线方向上的接触力;
⑵动态接触也就是所说的冲击,机械臂末端与外界环境的接触碰撞过程在极短的时间内完成,在此过程中可将力约束转化为接触速度的约束,来实现满足机械臂轨迹规划过程中冲击最优。
2、优化动态接触冲击的机械臂运动速度
⑴采用S型加减速曲线
采用S型速度变化使得运动过程中加速度连续,进而减小各关节的冲击,使机械臂运动和作业满足力的约束。
优化后的S型速度曲线共分为七个过程:
①加加速度一定加速度匀速增加过程,加速度从初始值加速到最大加速度;
②匀加加速度加速过程,此过程中保持最大加速度不变,速度匀速增加;
③减加速度加速过程,此过程中速度继续增加,加速度以恒定的加加速度减小直到速度加到最大速度值;
④匀速运动过程,此过程中保持最大速度运动,加速度为零;
⑤减减速度减速过程,此过程中速度减小,减速度以恒定的减减速度减小;
⑥匀减减速度减速过程,此过程中保持减速度不变,速度匀速减小;
⑦减减速度一定减速度匀速减小过程;
⑵归一化时间算子
笛卡尔空间位置插补的表达式为:
P(t)=Ps+l(t)(Pe-Ps)
其中,P(t)表示插补点的坐标;l(t)为归一化时间算子;Ps为期望位置;Pe为初始位置。
要对时间算子进行平滑处理,使得速度曲线平滑。定义归一化时间算子的表达式为:
其中,i=0,1,2,3,4,5,6;l(0)=0;l(T)=1;l(t)在区间内单调递增;S(t)为各时间段内对应的总位移;
3、基于力约束进行轨迹规划
⑴基于规则接触边缘具有力约束的连续轨迹规划:
规则接触边缘具有力约束的连续轨迹采用圆弧过渡的方法代替直线-直线连续轨迹规划方法进行平滑处理。达到运动在接触临界接触面时不具有沿接触面法线方向的位移、速度和加速度分量,过渡圆弧的联结点分别选取在前一段直线的减速段和后一段直线的加速段。
①在接触模型的有效的情况下,可得机械臂末端与外界环境的接触刚度为K,在已知约束力为F的情况下,可求得接触形变Xe
Xe=F/K
②通过接触形变,可得临界接触点B点的空间位置坐标p2
p2=pq+Xe
轨迹终点M的空间位置坐标p3为:
p3=p2+L
③过渡精度的确定,过渡精度决定过渡圆弧的大小,决定实际轨迹与期望轨迹逼近程度。设过渡力大小为F',则过渡精度re为:
re=F-F'
当力传感器示数时,进行空间直线轨迹运动;当力传感器示数时,为过渡圆弧的起点;当力传感器示数时,进行圆弧过渡运运动;当力传感器示数时,为过渡圆弧的终点,之后机械臂末端与外界环境保持力F进行沿接触面切平面方向的运动。
④确定过渡圆弧联结点pe1和pe2
当力传感器示数到达过渡力F'时开始圆弧过渡运动,当力传感器示数到达约束力F时结束圆弧过渡运动。即过渡圆弧联结点为接触力为F'和F时机械臂末端的位置点。
当力传感器示数时,通过位置编码器可知当前机械臂六个关节关节角θe11~θe16,通过机械臂正运动学式可得此刻对应的机械臂末端实际位置即为过渡圆弧起始联结点pe1。同理,当力传感器示数时,通过位置编码器可知当前机械臂六个关节关节角θe21~θe26,通过机械臂正运动学式可得此刻对应的机械臂末端实际位置即为过渡圆弧终止联结点pe2
⑤轨迹生成器通过路径点和过渡圆弧起始联结点pe1和终止联结点pe2的空间位置坐标,采用具有S型速度曲线拟合的空间直线和圆弧规划方法生成运动轨迹。
⑥运动过程中用力传感器进行监测,看整个运动过程是否满足期望力F的约束,如果时,重新规划。
⑵基于曲率接触边缘具有力约束的连续轨迹规划:
①机械臂实现沿接触面法线方向直线运动和沿临界安全接触面的接触运动和沿接触面法线方向远离临界安全接触面的直线运动,是一个直线-圆弧-直线连续的轨迹规划过程,运动过程中M、N点在安全接触圆柱面上,是具有力约束运动中的特殊点。
②采用圆弧过渡来解决直线-圆弧连续运动时产生的尖角问题。使得末端在沿安全接触面的运动过程中只有沿接触面切线方向位移和速度,法线方向速度为零,保证运动满足安全约束。
③通过调节过渡精度来调整实际轨迹与期望轨迹的精度。满足接触模型的基于圆柱接触面的连续轨迹规划方法与基于规则接触面连续轨迹规划方法中过渡圆弧联结点求取方法一致,不同的是过渡完后机械臂末端是在一定接触力下沿着接触面的圆弧轨迹。
1、具有接触形变约束的轨迹规划方法
接触模型的接触刚度稳定不变,将力约束转化为接触位移的约束来实现机械臂末端与外界环境之间的安全接触运动。具有力约束轨迹规划系统结构图如附1图所示,整个系统结构由接触模型、轨迹生成器、力位混合控制器及机械臂执行器组成。其中,Fd为期望力约束;Xe为通过接触模型得到的接触形变;Y为实际位置;F为机械臂腕部六维力传感器测量的机械臂末端施加的力。
1)、接触模型
将力约束转化为接触形变的约束进行机械臂与外界环境之间接触运动。由接触模型式得出,接触力不仅与接触形变有关还与接触刚度有关。机械臂末端与外界环境之间的接触刚度与机械臂自身末端刚度和环境刚度有关。对于机械臂末端的刚度与其姿态有关,在保持末端位置不变的情况下不同姿态所对应的末端刚度不同;环境刚度与接触环境的材料有关。对于刚度较大的接触,较小的形变就会产生较大的接触力。因此,要想实现力约束的转化就先要对接触模型进行验证。
本发明将采用实验测量的方法对接触模型进行验证,具体步骤如下:
(1)选运动方向沿接触面法线方向,设定机械臂末端起始位置X0,通过角位移编码器反馈读取当前实际六个关节角角位移值θ1~θ6,根据机械臂正运动学可知机械臂末端实际位置
(2)保持其他方向位移不变,让机械臂沿接触面法线方向运动期望距离ΔXi,与外界环境产生接触,读取当前力传感器示数Fi和当前实际六个关节角角位移值θi1~θi6,根据机械臂正运动学可知机械臂末端实际位置进行多组数据测量,每次的期望运动距离可不同。
(3)对采集到的Fi进行数据处理,验证两者是否满足线性关系,若满足则证明接触模型正确,进而可以得到接触刚度值,实现力约束向接触形变约束的转化。
2)、轨迹生成器
轨迹生成器可以根据所给的路径点生成期望轨迹。对于不同的环境不同的作业任务,机械臂末端与外界物体之间接触运动的轨迹不同,所选用的轨迹生成方法也不同。若要实现接触运动直线轨迹,则选用轨迹生成器中直线插补方法;若要求在接触运动过程中实现圆弧轨迹,则选择圆弧插补方法;要若要在接触过程中实现自由曲线轨迹,则选用样条函数插值方法。
接触运动学,当接触运动过程是沿切线方向运动和法线方向运动的组合运动时,不仅需要考虑轨迹的连续性还要使得整个运动过程在力的安全约束范围内。对于此要求,本发明将采用运动过渡的方法实现具有力约束的轨迹规划。在接触模型成立的基础上,可以将接触力的约束转化为接触形变约束,进而可以将力安全阈值范围转化为过渡轨迹。在有力约束的情况下,对于沿接触面法线方向运动与沿接触面切平面方向运动轨迹的过渡过程中,解决轨迹尖角问题是实现满足要求运动的关键。对于连续轨迹中的尖角问题处理有圆弧过渡和样条函数二次插补过渡两种方法。根据过渡轨迹要求不同可采取不同的过渡方法,样条函数二次插补适合对过渡轨迹精度要求较高或过渡轨迹曲线不确定的运动,相比圆弧过渡,样条函数二次插补过渡方法计算量大,效率低。因此,本文将采用圆弧过渡的方法实现具有力约束的轨迹规划。
由于机械臂末端刚度与机械臂末端姿态有关,姿态的变化会影响接触刚度的变化,导致末端接触力变化,末端姿态的变化还会产生与接触面之间滚动接触运动,使得机械末端与外界环境之间不仅有力的作用还有力矩的作用,是一个十分复杂的问题。本文这里假设运动过程中机械臂末端相对于接触环境的姿态保持不变即保持接触刚度不变,只考虑具有沿接触面法线方向力约束情况和与接触面之间滑动接触运动。
3)、力位混合控制器
轨迹生成器根据给定的路径点生成笛卡尔空间期望轨迹,并通过机械臂逆运动学转化为关节空间轨迹。将对应的关节变量输入给力位混合控制器,通过控制器实现机械臂运动。
2、基于规则接触边缘具有力约束的连续轨迹规划方法
本发明选择空间平面为规则接触面进行研究,图3机械臂末端与接触面的空间示意图。运动情况描述为:机械臂末端从空间起始点沿接触面法线方向运动到规则接触面,在具有一定接触力情况下沿着接触面运动。为了方面研究,将机械臂末端与环境的接触看作为点接触,接触过程中假设接触刚度不发生变化,研究滑动接触运动过程中连续轨迹规划方法。
在接触模型成立的情况下,由赫兹接触理论接触力的计算公式可知,在接触刚度一定,力指数为1的情况下,接触形变与接触力成线性关系。对于具有规则接触面的平面,已知最大接触力约束,可得规则接触平面对应的临界安全接触面是一个与规则接触面平行的平面如图2所示。图2中,A为起始点,Q为接触面上的点,B、M为临界安全接触面上的点,要求机械臂末端实现沿接触面法线和接触面切平面连续直线轨迹规划任务,即AB-BM连续直线轨迹。在AB-BM连续轨迹中B点在临界安全接触面上,是满足力约束运动中的一个特殊点。
对于空间连续直线运动,理论上可以采取单独规划每条独立直线的方法,但是这样会使得电机在各段直线间不停的启动和停止,实际中这种操作会产生较大的振动和摩擦[51]。若进行连续的规划,由图2可以看出,在直线与直线的联结处即B点处会有尖角产生,尖角的存在会使得位移、速度方向发生突变,不但轨迹不平滑而且会造成运动超出安全接触形变的约束。对于此种情况,本文将采用圆弧过渡的方法对直线-直线连续轨迹进行平滑处理,过渡圆弧的联结点分别在前一段直线的减速段和后一段直线的加速段,圆弧过渡示意图如图2中轨迹所示,这样不仅可以保证B点的特殊要求,并使得运动在接触临界接触面时不具有沿接触面法线方向的位移、速度和加速度分量。
对于图3中所描述的运动,满足接触模型的基于规则接触平面的连续轨迹规划方法如下:
已知起始点A的位置坐标为p1,外界接触物体上一点Q坐标pq,期望约束力为F,沿接触面切平面运动距离为L。
(1)在接触模型的有效的情况下,可得机械臂末端与外界环境的接触刚度为K,在已知约束力为F的情况下,可求得接触形变Xe
Xe=F/K
(2)通过接触形变,可得临界接触点B点的空间位置坐标p2
p2=pq+Xe
轨迹终点M的空间位置坐标p3为:
p3=p2+L
(3)过渡精度的确定,过渡精度决定过渡圆弧的大小,决定实际轨迹与期望轨迹逼近程度。设过渡力大小为F',则过渡精度re为:
re=F-F'
当力传感器示数时,进行空间直线轨迹运动;当力传感器示数时,为过渡圆弧的起点;当力传感器示数时,进行圆弧过渡运运动;当力传感器示数时,为过渡圆弧的终点,之后机械臂末端与外界环境保持力F进行沿接触面切平面方向的运动。
(4)确定过渡圆弧联结点pe1和pe2
当力传感器示数到达过渡力F'时开始圆弧过渡运动,当力传感器示数到达约束力F时结束圆弧过渡运动。即过渡圆弧联结点为接触力为F'和F时机械臂末端的位置点。
当力传感器示数时,通过位置编码器可知当前机械臂六个关节关节角θe11~θe16,通过机械臂正运动学式可得此刻对应的机械臂末端实际位置即为过渡圆弧起始联结点pe1。同理,当力传感器示数时,通过位置编码器可知当前机械臂六个关节关节角θe21~θe26,通过机械臂正运动学式可得此刻对应的机械臂末端实际位置即为过渡圆弧终止联结点pe2
(5)轨迹生成器通过路径点A、Q、M和过渡圆弧起始联结点pe1和终止联结点pe2的空间位置坐标,采用具有S型速度曲线拟合的空间直线和圆弧规划方法生成运动轨迹。
(6)运动过程中用力传感器进行监测,看整个运动过程是否满足期望力F的约束,如果时,重新规划。
需要说明的是,在机械臂末端接触运动过程中会产生摩擦力,本文暂不考虑摩擦力的影响;对于不同接触刚度的物体,力约束下的接触形变大小不同,运动中可以通过调整过渡精度来调整跟期望轨迹的逼近程度。
3.基于曲率接触边缘具有力约束的连续轨迹规划方法
本发明选择圆柱体为具有一定曲率的接触面进行研究,图4为机械臂末端与接触面的空间示意图。运动情况描述为:机械臂末端从空间中起始点沿接触法线方向运动到规则接触面,在具有一定接触力情况下沿着接触面运动。为了方面研究,将机械臂末端与环境的接触看作为点接触,接触过程中假设接触刚度不发生变化,研究滑动接触运动过程中连续轨迹规划方法。
由赫兹接触理论接触力的计算公式可知,在接触刚度一定,力指数为1的情况下,接触形变与接触力成线性关系。对于圆柱接触面,已知最大安全接触力,可知临界安全接触面也是一个与圆柱面通圆心的圆柱面投影图4所示。只要接触运动过程中产生的位移不超过临界安全接触面,即可使运动满足力的约束。图中A为起始点,Q为接触面上的点,M、U、N为临界安全接触面上的点,要求机械臂实现沿接触面法线方向直线运动和沿临界安全接触面的接触运动和沿接触面法线方向远离临界安全接触面的直线运动,即AM-MN-NB运动,是一个直线-圆弧-直线连续的轨迹规划过程,运动过程中M、N点在安全接触圆柱面上,是具有力约束运动中的特殊点。规划的关键是解决直线与圆弧连接时在接触点出现的尖角问题,进而解决特殊点位置、速度、加速度约束问题,使得运动满足安全要求。本发明采用圆弧过渡来解决直线-圆弧连续运动时产生的尖角问题,如图4中蓝色轨迹所示。这样使得末端在沿安全接触面的运动过程中只有沿接触面切线方向位移和速度,法线方向速度为零,保证运动满足安全约束。可通过调节过渡精度来调整实际轨迹与期望轨迹的精度。满足接触模型的基于圆柱接触面的连续轨迹规划方法与基于规则接触面连续轨迹规划方法中过渡圆弧联结点求取方法一致,不同的是过渡完后机械臂末端是在一定接触力下沿着接触面的圆弧轨迹。

Claims (1)

1.一种六自由度机械臂考虑力约束的轨迹规划方法,其特征在于,包括如下步骤:
(1)将力约束转化为接触形变和接触运动速度的约束;
根据接触模型将机械臂末端与环境的接触分为静态接触和动态接触;
(1.1)静态接触过程中,将力的约束转化为沿法线方向的接触位移的约束,即
其中,Ks是接触刚度;XE是机械臂末端沿接触面法向的形变量;F表示机械末端力传感器测量的接触面法线方向上的接触力;
(1.2)动态接触也就是冲击,机械臂末端与外界环境的接触碰撞过程在极短的时间内完成,将力约束转化为接触速度的约束,来实现满足机械臂轨迹规划过程中冲击最优;
(2)优化动态接触冲击的机械臂运动速度;
(2.1)采用S型加减速曲线使加速度连续;
采用S型速度变化使得运动过程中加速度连续,减小各关节的冲击,使机械臂运动和作业满足力的约束;
优化后的S型速度曲线共分为七个过程:
(2.1.1)加加速度一定、加速度匀速增加过程,加速度从初始值加速到最大加速度;
(2.1.2)匀加加速度加速过程,保持最大加速度不变,速度匀速增加;
(2.1.3)减加速度加速过程,速度继续增加,加速度以恒定的加加速度减小直到速度加到最大速度值;
(2.1.4)匀速运动过程,保持最大速度运动,加速度为零;
(2.1.5)减减速度减速过程,速度减小,减速度以恒定的减减速度减小;
(2.1.6)匀减减速度减速过程,保持减速度不变,速度匀速减小;
(2.1.7)减减速度一定减速度匀速减小过程;
(2.2)归一化时间算子;
笛卡尔空间位置插补的表达式为:
P(t)=Ps+l(t)(Pe-Ps),
其中,P(t)表示插补点的坐标;l(t)为归一化时间算子;Ps为期望位置;Pe为初始位置;
要对时间算子进行平滑处理,使得速度曲线平滑,定义归一化时间算子的表达式为:
其中,i=0,1,2,3,4,5,6;l(0)=0;l(T)=1;l(t)在区间内单调递增;S(t)为各时间段内对应的总位移;
(3)基于力约束进行轨迹规划;
(3.1)基于规则接触边缘具有力约束的连续轨迹规划:
规则接触边缘具有力约束的连续轨迹采用圆弧过渡的方法代替直线-直线连续轨迹规划方法进行平滑处理,达到运动在接触临界接触面时不具有沿接触面法线方向的位移、速度和加速度分量,过渡圆弧的联结点分别选取在前一段直线的减速段和后一段直线的加速段;
(3.1.1)在接触模型的有效的情况下,得到机械臂末端与外界环境的接触刚度为K,在已知约束力为F的情况下,可求得接触形变Xe
Xe=F/K;
(3.1.2)通过接触形变,得到临界接触点B点的空间位置坐标p2
p2=pq+Xe
轨迹终点M的空间位置坐标p3为:
p3=p2+L;
(3.1.3)过渡精度的确定,过渡精度决定过渡圆弧的大小,决定实际轨迹与期望轨迹逼近程度;设过渡力大小为F',则过渡精度re为:
re=F-F'
当力传感器示数时,进行空间直线轨迹运动;当力传感器示数时,为过渡圆弧的起点;当力传感器示数时,进行圆弧过渡运动;当力传感器示数时,为过渡圆弧的终点,之后机械臂末端与外界环境保持力F进行沿接触面切平面方向的运动;
(3.1.4)确定过渡圆弧起始联结点pe1和过渡圆弧终止联结点pe2
当力传感器示数到达过渡力F'时开始圆弧过渡运动,当力传感器示数到达约束力F时结束圆弧过渡运动;即过渡圆弧联结点为接触力为F'和F时机械臂末端的位置点;
当力传感器示数时,通过位置编码器可知当前机械臂六个关节关节角θe11~θe16,通过机械臂正运动学式得到此刻对应的机械臂末端实际位置即为过渡圆弧起始联结点pe1;当力传感器示数时,通过位置编码器可知当前机械臂六个关节关节角θe21~θe26,通过机械臂正运动学式得到此刻对应的机械臂末端实际位置即为过渡圆弧终止联结点pe2
(3.1.5)轨迹生成器通过路径点和过渡圆弧起始联结点pe1和终止联结点pe2的空间位置坐标,采用具有S型速度曲线拟合的空间直线和圆弧规划方法生成运动轨迹;
(3.1.6)运动过程中用力传感器进行监测,看整个运动过程是否满足期望力F的约束,如果时,重新规划;
(3.2)基于曲率接触边缘具有力约束的连续轨迹规划:
(3.2.1)机械臂实现沿接触面法线方向直线运动和沿临界安全接触面的接触运动和沿接触面法线方向远离临界安全接触面的直线运动,是一个直线-圆弧-直线连续的轨迹规划过程,确定运动过程中在安全接触圆柱面上的具有力约束运动中的特殊点M、N点;
(3.2.2)采用圆弧过渡来解决直线-圆弧连续运动时产生的尖角问题;使得末端在沿安全接触面的运动过程中只有沿接触面切线方向位移和速度,法线方向速度为零;
(3.2.3)通过调节过渡精度来调整实际轨迹与期望轨迹的精度,满足接触模型的基于圆柱接触面的连续轨迹规划方法与基于规则接触面连续轨迹规划方法中过渡圆弧联结点求取方法一致,过渡完后机械臂末端是在一定接触力下沿着接触面的圆弧轨迹。
CN201610421404.7A 2016-06-12 2016-06-12 一种六自由度机械臂考虑力约束的轨迹规划方法 Active CN106113034B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610421404.7A CN106113034B (zh) 2016-06-12 2016-06-12 一种六自由度机械臂考虑力约束的轨迹规划方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610421404.7A CN106113034B (zh) 2016-06-12 2016-06-12 一种六自由度机械臂考虑力约束的轨迹规划方法

Publications (2)

Publication Number Publication Date
CN106113034A CN106113034A (zh) 2016-11-16
CN106113034B true CN106113034B (zh) 2018-06-12

Family

ID=57270136

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610421404.7A Active CN106113034B (zh) 2016-06-12 2016-06-12 一种六自由度机械臂考虑力约束的轨迹规划方法

Country Status (1)

Country Link
CN (1) CN106113034B (zh)

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10456103B2 (en) 2016-06-14 2019-10-29 Shanghai United Imaging Healthcare Co., Ltd. Unified trajectory generation process and system
CN107139171B (zh) * 2017-05-09 2019-10-22 浙江工业大学 一种基于力矩控制的工业机器人避障轨迹规划方法
CN107290959B (zh) * 2017-06-08 2020-06-02 浙江理工大学 一种基于位移等效的机器人笛卡尔空间速度优化方法
CN107436555B (zh) * 2017-07-28 2019-12-13 珞石(山东)智能科技有限公司 非零初始状态s型停止曲线速度规则及其在线求解方法
WO2019127443A1 (en) * 2017-12-29 2019-07-04 Shenzhen United Imaging Healthcare Co., Ltd. System and method for synchronous motion optimization of device with moving components
CN108489652B (zh) * 2018-04-11 2020-04-21 中铁第四勘察设计院集团有限公司 第三轨与受流器之间接触力的确定方法及装置
CN110524371B (zh) * 2019-08-26 2021-05-04 南京工业大学 基于实时力控的复杂曲面恒定切除率的机器人打磨方法
CN112157661B (zh) * 2020-12-02 2021-03-05 成都卡诺普自动化控制技术有限公司 一种基于机器人动力学的关节运动轨迹优化方法
CN112506140B (zh) * 2020-12-02 2021-11-12 宁夏大学 一种五轴联动水切割机床的空间圆弧插补方法及系统
CN112775976B (zh) * 2021-02-05 2022-05-10 深圳市优必选科技股份有限公司 任务执行控制方法、装置、控制设备及可读存储介质
CN113040915A (zh) * 2021-03-12 2021-06-29 杭州柳叶刀机器人有限公司 机器人安全边界控制装置及方法、电子设备和存储介质
CN113017941A (zh) * 2021-03-12 2021-06-25 杭州柳叶刀机器人有限公司 机械臂作用力交互控制方法及装置、电子设备和存储介质
CN116512254B (zh) * 2023-04-11 2024-01-23 中国人民解放军军事科学院国防科技创新研究院 基于方向的机械臂智能控制方法及系统、设备、存储介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101733750A (zh) * 2009-11-05 2010-06-16 浙江大学 基于动态响应谱的关节含间隙机器人轨迹规划方法
CN104020770A (zh) * 2014-06-13 2014-09-03 哈尔滨工程大学 基于多项式的uuv空间轨迹规划方法
CN105234930A (zh) * 2015-10-15 2016-01-13 哈尔滨工程大学 基于构形平面的冗余机械臂运动控制方法
CN105583824A (zh) * 2016-01-26 2016-05-18 清华大学 力控牵引和摆位的多自由度机械臂控制装置及其控制方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3790816B2 (ja) * 2002-02-12 2006-06-28 国立大学法人 東京大学 人型リンク系の運動生成方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101733750A (zh) * 2009-11-05 2010-06-16 浙江大学 基于动态响应谱的关节含间隙机器人轨迹规划方法
CN104020770A (zh) * 2014-06-13 2014-09-03 哈尔滨工程大学 基于多项式的uuv空间轨迹规划方法
CN105234930A (zh) * 2015-10-15 2016-01-13 哈尔滨工程大学 基于构形平面的冗余机械臂运动控制方法
CN105583824A (zh) * 2016-01-26 2016-05-18 清华大学 力控牵引和摆位的多自由度机械臂控制装置及其控制方法

Also Published As

Publication number Publication date
CN106113034A (zh) 2016-11-16

Similar Documents

Publication Publication Date Title
CN106113034B (zh) 一种六自由度机械臂考虑力约束的轨迹规划方法
CN106041926B (zh) 一种基于卡尔曼滤波器的工业机械臂力/位置混合控制方法
CN106625671B (zh) 一种空间机器人抓捕翻滚目标的最优轨迹规划方法
CN107139171B (zh) 一种基于力矩控制的工业机器人避障轨迹规划方法
CN104932506B (zh) 基于快速终端滑模的轮式移动机器人轨迹跟踪方法
CN110653826B (zh) 一种面向传送带目标跟随的机器人实时轨迹规划方法
JP5327722B2 (ja) ロボットの負荷推定装置及び負荷推定方法
JP2019010713A (ja) ハンド制御装置、ハンド制御方法、およびハンドのシミュレーション装置
CN104570735A (zh) 一种码垛机器人弧形运动轨迹算法
CN106227154A (zh) 一种两轴联动机械臂运动控制的同步误差补偿方法
CN101791800A (zh) 一种双轮差动式机器人运动控制方法
CN111469129A (zh) 一种基于双机械臂的船体外板曲面成形无碰撞运动规划方法
CN104647379A (zh) 一种驱动器非线性条件下的双臂机器人运动控制方法
CN104090492B (zh) 基于指数函数的scara机器人ptp轨迹规划
Ahmad et al. Shape recovery from robot contour-tracking with force feedback
TW201426229A (zh) 控制裝置、控制方法及位置命令補償方法
CN111515954B (zh) 一种机械臂高质量运动路径生成方法
CN105904458B (zh) 一种基于复杂操作任务的非完整遥操作约束控制方法
CN104564913A (zh) 并联六自由度转台液压控制系统
CN109483551B (zh) 一种控制机器人多轴运动的方法、装置及系统
CN116300445A (zh) 一种汽车电泳涂装输送混联机构的任务空间新型综合误差快速连续滑模控制方法
Hu et al. Oscillation suppressing for an energy efficient bridge crane using input shaping
CN115030706A (zh) 基于锚固工艺的钻锚机器人双钻臂协同控制系统及方法
CN116749150B (zh) 一种多轴机器人系统运动规划方法、电子设备及介质
CN113946136A (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