冗余度机器人最小力矩与回避奇异并行优化
- 格式:pdf
- 大小:219.17 KB
- 文档页数:4
第38卷第1期 2 0 0 6年1月 哈尔滨工业大学学报 JOURNAL OF HARBIN INSTITUTE OF rIECHNOLOGY V0I.38 No.1 Jan.2oo6
冗余度机器人最小力矩与回避奇异并行优化
刘 宇 ,姜艳妹 ,孙立宁
(1.哈尔滨工业大学机器人研究所,黑龙江哈尔滨150001,E—mail:lyu11@hit.edu.cn:
2.哈尔滨理工大学自动化系,黑龙江哈尔滨150080)
摘要:针对冗余度机器人梯度投影法最小力矩局部优化所发生的轨迹不稳定现象,经过讨论分析后,明确 指出这种不稳定是由关节构形奇异引起的,并提出了一种冗余度机器人最小力矩与回避奇异的’并行优化方 案.通过对一三自由度冗余度机器人的仿真,证实了所提优化方案是可行、有效的. 关键词:梯度投影;冗余度机器人;回避奇异;并行优化 中图分类号:TP242.2 文献标识码:A 文章编号:0367—6234(2006)01—0001—03
Research on parallel optimization scheme between singularity
avoidance and minimum Torque for a Redundant Manipulator
LIU Yu ,JIANG Yan-shu ,SUN Li.ning
1.Robotics Institute,Harbin Institute of Technology,Harbin 150001,China,E—email:lyul 1@hit.edu.cn; 2.Dept.of Automation,Harbin University of Science and Technology,Harbin 150080,China)
Abstract:Aiming at unstable trajectory of minimum torque optimization based on gradient projection,this pa—
per through analysis explicitly points out that instability is caused by joint configuration singularity.Therefore, it presents an optimization scheme between Singularity Avoidance and Minimum Torque for a Redundant Ma-
nipulator.At last,through simulation of a 3一DOF redundant manipulator,the scheme proves to be feasible
and effective.
Key words:gradient projection;redundant manipulators;singularity avoidance;Parallel Optimization
目前,冗余度机器人在回避障碍、回避奇异以
及回避关节极限等运动学优化方面取得了卓有成 效的进展¨ .与此同时,对于动力学优化也展开
了积极的研究.动力学优化主要包括关节力矩优
化和最小能量优化.关节力矩优化从合理分配冗
余度机器人载荷出发,利用冗余自由度使得机器
人各个关节不发生力矩超限或使关节力矩尽可能 介于上限和下限之间.最小能量优化从能耗降低
的角度来优化机器人的动力学性能,合理配置了
冗余度机器人各轴的轨迹输出,使惯量大的关节
获得较小的关节加速度.本文主要探讨冗余度机 器人最小关节力矩规划问题.
收稿日期:2004—09—12. 基金项目:国家高技术研究发展计划资助项目(2004AA421030) 作者简介:刘宇(1971一),男,博士,讲师; 孙立宁(1964一)。男。博士。教授。博士生导师. 1冗余度机器人局部最小关节力矩规划
对于一个冗余自由度机器人,其关节空间和
作业空间之间的速度描述可用下式表示:
=J . (1)
式中:戈为末端操作器的速度矢量;J为雅克比矩
阵; 为关节空间的速度矢量.由(1)式,容易得到 基于伪逆的关节速度通用局部最优解
0=., 戈+ (,一., .,)V H. (2) 式中: 为实标量放大系数;J 为雅克比伪逆阵; (,一., .,)为零空间梯度投影矩阵;7H为任意性
能准则梯度.
对式(1)两边求导得
=J +t=『0. (3) 根据式(3),可容易得到基于伪逆的关节加速度
通用局部最优解为 维普资讯 http://www.cqvip.com ・2・ 哈尔滨工业大学学报 第38卷
=J ( 一.7 )+(,一., J) . (4)
式中: 为任意矢量. 机器人动力学方程可用下式表述:
M(0)舀+C(0, )+N(0, )= . (5)
式中:M(0)为n×n维机器人惯性矩阵;C(0, )
为n维哥氏力及离心力矢量;Ⅳ(0, )为n维有势
力矢量; 为n维关节扭矩. Ballieul、Hollerbach等最先基于零空间梯度
投影研究了冗余度机器人局部最小关节力矩规划 问题,并提出了如下关节扭矩优化准则 :
rain G= T . (6) 将式(5)代入到式(6)中可得
G=ll M +c+Ⅳl (7) 将式(4)代入到式(7)中可得
G=ll』l ( 一.7 )+c+Ⅳ+M(I—J .,) ll .
对上式求最小2一范数解,可得零空间矢量 为
=一[ (,一., .,)] [』l ( 一.7扫)+c+Ⅳ], 代入到式(4)中,注意到(,一., .,)[ (,一., .,)] =[ (,一., .,)] 可得关节加速度
l;=., ( 一.7扫)一[ (,一., .,)] [』l ( 一
.7 )+C+Ⅳ]. 再代入到式(5)中,可得最小关节扭矩解为
=MJ ( 一.7扫)+C+Ⅳ一 [ (,一
., .,)] [MJ ( 一.7扫)+c+Ⅳ]. (8)
上面的求解并未考虑具体每个关节扭矩的差 异.因此,Hollerbach和Suh研究了一个具有n自由
度的冗余度机器人,其关节扭矩的上、下限分别为
r 和r一,期望的末端轨迹为 ( ),要求实现在满
足末端操作器作业轨迹的前提下,同时减少对关节
力矩的需求.为达到此目的,应尽可能使关节扭矩
接近于扭矩极限的中点1/2( + 一).因此提出
了rain F=ll 一( + 一)/2ll 优化准则 .
利用与上述类似的方法,可以得到相应于式
(8)的最小关节力矩解 = 。+ [ (,一.,
.,)] ( + i)/2.式中:r。=MJ ( 一.7扫)+
C(0,扫)+N(0).
2 冗余度机器人局部最小力矩解的
不稳定性
利用局部最小力矩解规划冗余度机器人的关 节加速度在实践中出现不稳定的问题,尤其对于
比较长的轨迹.经过进一步研究发现,最小力矩解 在长轨迹的情况下导致了很大的关节速度,而这 个大的速度又需要更大的扭矩来维持期望的末端
轨迹,结果超过关节扭矩极限,造成了运动抽搐现 象引起了不稳定.文献[4]提出了关节扭矩最小
化和关节加速度融合的一种平衡方法.文献[5]
提出了一种零空间阻尼的方法,并以所提出的方
法与未加阻尼的最小力矩规划方法进行了对比模
拟研究,结果发现,其所提出的零空间阻尼优化力 矩的方法在运动中表现出很强的稳定性,明显优
于未加阻尼的方式,而且对于周期循环的运动具
有关节保持性,但这种方法在一定程度上牺牲了 对关节力矩的优化.文献[6]认为不稳定性是由
齐次解部分引起的,因而提出了一种在齐次解和
最小二乘解之间的切换技术,其切换的阈值根据
关节加速度范数而定. 然而,这种力矩不稳定现象主要由关节速度
过大引起的,而过大的关节速度一般在关节奇异
构形附近才会发生.因此,这种最小力矩优化的不 稳定性源自关节的奇异构形.在最小力矩的优化
中,不妨加入回避关节奇异准则,实现运动学和动
力学的均衡优化,可获得稳定的运动轨迹.
3 回避奇异和最小力矩并行的优化方案
这里以条件数作为衡量奇异程度的可操作测 度.假定其性能函数为U,则U(0)= …/o" i .
式中: , 为.,的最大、最小奇值. 为了达到动力学运动规划的目的,应当将回
避奇异优化投影在加速度级上面.故此,由式(3)
容易得到相应于式(2)的回避奇异性能准则的加 速度级优化.
=J ( 一.7扫)+k(,一., .,)△ (0)扫. (9) 将式(3)代入到式(5)中可得
一 扫+.,。(C+N)=Jo . (10)
式中:., 为基于惯量的雅克比矩阵,J。=埘~.
由式(10)可直接推导得到基于伪逆的关节 扭矩通用局部最优解为
r= [ 一.7扫+ (c+Ⅳ)]+p(1一 )
(11)
式中:西为用来优化附加性能准则的任意矢量;P 为调节因子.将式(9)代入到式(5)中,可得
r=M[J ( 一.7 )+k(,一
., .,)△ (0)扫]+c+Ⅳ. (12)
维普资讯 http://www.cqvip.com 第1期 刘宇,等:冗余度机器人最小力矩与回避奇异并行优化
将式(12)代入到式(1 1)中可求得零空间矢量为
=(I—J:J MJ _jc +k(、l—J:Ju M(、I—J
J)AU +(,一.,:J )(C+N). (13)
这里要用到(,一 Ja) =(,一 Jo),(,一
J )(,一. J )=(,一.,:J ).
方程(13)是基于最小扭矩范数表述的优化 可操作测度所获得的零空间投影矢量.将式(13)
代入到式(11)中可得
.r={.,:+P(,一.,:J )MJ }( 一.:『 )+
P k(、I—J:J M(、l—J J AU +
{p,+(1一p).,:J }(C+Ⅳ)
根据J =JM 以及 =,,得到力矩优化 和可操作测度优化的平衡解
.r=P{MJ ( 一 )+kM(,一., .,)AU +C+
Ⅳ}+(1一P){.,:( 一.:『 )+. Jo(C+Ⅳ)}.
4 回避奇异和最小力矩并行的优化
方案仿真研究
为了证实上述所提方案的正确性,简化起见,
以一个三自由度冗余机器人为例进行动力学优化
仿真研究.如图l所示,容易得到此冗余度机器人
籁 世 的雅克比矩阵为
.,: +/2Cl2+/3 ̄123 12c ̄2+13C123 13C1 ̄1. t一2l l—l2s12—23 l 一l2 l2一l3sl丑 一l3 l J
中:c1,2.…,。=cos(0l+ +…+0 ), 1.2.…. = sin(0l+02+…+0 ).
O x/m 图l三自由度冗余机器人 假定机器人的臂长分别为z =z =z =l m, 质量均为6.4 kg,初始角0l= 6,0:= ,03=
2/5'rr,初始角速度均为零,加速度矢量为[A-/2,
/2]rad/s ,采用bang—bang类型,即在仿真中间
时刻加速度反向.仿真时间2s,采样时间2 ms.
根据上述初始条件,当P=0.25时,对回避奇