工程科学学报,第38卷,第11期:1643-1651,2016年11月 Chinese Journal of Engineering,Vol.38,No.11:1643-1651,November 2016 D0l:10.13374/j.issn2095-9389.2016.11.019:http://journals.ustb.edu.cn 高压输电线路带电检修机器人作业臂运动优化 江 维⑧,吴功平,樊飞,张颉 武汉大学动力与机械学院,武汉430072 ☒通信作者,E-mail:jiangwei2013@whu.cdu.cn 摘要针对完全依靠人工带电拧紧高压输电线路引流板螺栓作业效率低、劳动强度大、高空和高电压危险的问题,研制了 一种双臂和双机械手的引流板螺栓拧紧带电检修机器人·通过各关节的轨迹规划双作业臂及其末端由初始位姿运动到螺栓 螺母对准状态是机器人顺利完成检修作业的关键.针对现有多项式插值关节轨迹规划依赖于轨迹端点时刻,导致该方法实 用性差及关节轨迹运动过程中忽略驱动机构对关节状态制约的问题,提出基于M一Max时间标准化的改进多项式插值关节 轨迹规划方法.基于该方法,以关节运动时间为优化对象,提出满足关节轨迹运动全程约束条件的关节运动时间区间范围的 求取方法.通过仿真实验验证了改进算法关节轨迹仅与轨迹端点状态及运动时间有关,而与端点时刻无关,进一步淡化轨迹 端点时刻对关节轨迹的影响,提高算法的实用性.通过选取最优关节轨迹运动时间,满足了关节状态全程状态约束要求,从而避 免超调的发生,优化了各关节运动轨迹,提高了关节运动的效率.最后通过现场试验进一步验证改进算法的工程实用性 关键词高压输电线路:机器人;作业臂控制:运动规划:优化 分类号TP242.3 Working manipulator motion optimization of live maintenance robot for high voltage transmission line JIANG Wei,WU Gong-ping,FAN Fei,ZHANG Jie School of Power and Mechanical Engineering,Wuhan University,Wuhan 430072,China Corresponding author,E-mail:jiangwei2013@whu.edu.cn ABSTRACT In response to the problems of artificial low-operation efficiency,high labor intensity,high altitude and high risk in drainage board bolt tightening for high voltage transmission line,a bolt tightening live maintenance robot with double arms and double manipulators was designed in this paper.The key to the successful completion of maintenance work is that double working arms and their ends move from the initial position to the bolt and nut alignment state through the trajectory planning of each joint of the robot. The existed polynomial interpolation of joint trajectory planning depends on trajectory endpoint time,leading to some issues such as the low practicability of the method and the neglect of the constraints of the drive mechanism on the joint state in joint trajectory motion,so an improved polynomial interpolation joint trajectory planning method based on Min-Max time standardization was presented.Based on the method,a solving method of joint motion time ranges which meet the requirements of joint trajectory motion constraints was proposed with articulation motion time as the optimization target.Simulation results show that the joint trajectory of the improved algo- rithm is only with regard to trajectory endpoint state and motion time,but not to trajectory endpoint time,which further dilutes the influence of trajectory endpoint time on the joint trajectory,and the practicality of the algorithm is also improved.By selecting the opti- mal trajectory articulation motion time which meets the requirements of whole joint state constraints,the improved algorithm not only avoids the occurrence of overshooting and optimizes the trajectory of each joint,but also improves the efficiency of joint motion.Final- ly,the engineering practicality of the improved algorithm was verified by field operation test. 收稿日期:2015-12-31 基金项目:国家自然科学基金资助项目(51105281):中央高校基本科研业务费专项资金资助项目(2104005):湖南省电力公司科技项目 (5216A01400B1)
工程科学学报,第 38 卷,第 11 期: 1643--1651,2016 年 11 月 Chinese Journal of Engineering,Vol. 38,No. 11: 1643--1651,November 2016 DOI: 10. 13374 /j. issn2095--9389. 2016. 11. 019; http: / /journals. ustb. edu. cn 高压输电线路带电检修机器人作业臂运动优化 江 维,吴功平,樊 飞,张 颉 武汉大学动力与机械学院,武汉 430072 通信作者,E-mail: jiangwei2013@ whu. edu. cn 摘 要 针对完全依靠人工带电拧紧高压输电线路引流板螺栓作业效率低、劳动强度大、高空和高电压危险的问题,研制了 一种双臂和双机械手的引流板螺栓拧紧带电检修机器人. 通过各关节的轨迹规划双作业臂及其末端由初始位姿运动到螺栓 螺母对准状态是机器人顺利完成检修作业的关键. 针对现有多项式插值关节轨迹规划依赖于轨迹端点时刻,导致该方法实 用性差及关节轨迹运动过程中忽略驱动机构对关节状态制约的问题,提出基于 Min--Max 时间标准化的改进多项式插值关节 轨迹规划方法. 基于该方法,以关节运动时间为优化对象,提出满足关节轨迹运动全程约束条件的关节运动时间区间范围的 求取方法. 通过仿真实验验证了改进算法关节轨迹仅与轨迹端点状态及运动时间有关,而与端点时刻无关,进一步淡化轨迹 端点时刻对关节轨迹的影响,提高算法的实用性. 通过选取最优关节轨迹运动时间,满足了关节状态全程状态约束要求,从而避 免超调的发生,优化了各关节运动轨迹,提高了关节运动的效率. 最后通过现场试验进一步验证改进算法的工程实用性. 关键词 高压输电线路; 机器人; 作业臂控制; 运动规划; 优化 分类号 TP242. 3 收稿日期: 2015--12--31 基金项目: 国家自然科学基金资助项目( 51105281) ; 中央高校基本科研业务费专项资金资助项目( 2104005) ; 湖南省电力公司科技项目 ( 5216A01400B1) Working manipulator motion optimization of live maintenance robot for high voltage transmission line JIANG Wei ,WU Gong-ping,FAN Fei,ZHANG Jie School of Power and Mechanical Engineering,Wuhan University,Wuhan 430072,China Corresponding author,E-mail: jiangwei2013@ whu. edu. cn ABSTRACT In response to the problems of artificial low-operation efficiency,high labor intensity,high altitude and high risk in drainage board bolt tightening for high voltage transmission line,a bolt tightening live maintenance robot with double arms and double manipulators was designed in this paper. The key to the successful completion of maintenance work is that double working arms and their ends move from the initial position to the bolt and nut alignment state through the trajectory planning of each joint of the robot. The existed polynomial interpolation of joint trajectory planning depends on trajectory endpoint time,leading to some issues such as the low practicability of the method and the neglect of the constraints of the drive mechanism on the joint state in joint trajectory motion,so an improved polynomial interpolation joint trajectory planning method based on Min--Max time standardization was presented. Based on the method,a solving method of joint motion time ranges which meet the requirements of joint trajectory motion constraints was proposed with articulation motion time as the optimization target. Simulation results show that the joint trajectory of the improved algorithm is only with regard to trajectory endpoint state and motion time,but not to trajectory endpoint time,which further dilutes the influence of trajectory endpoint time on the joint trajectory,and the practicality of the algorithm is also improved. By selecting the optimal trajectory articulation motion time which meets the requirements of whole joint state constraints,the improved algorithm not only avoids the occurrence of overshooting and optimizes the trajectory of each joint,but also improves the efficiency of joint motion. Finally,the engineering practicality of the improved algorithm was verified by field operation test.
·1644 工程科学学报,第38卷,第11期 KEY WORDS high voltage transmission line;robot;manipulator control:motion planning:optimization 带电检修·是指在高压电气设备上不停电进行 状态及关节运动时间有关,而与作业臂关节轨迹端点 检修、测试的一种作业方法,其对电网稳定运行,确保 时刻没有直接联系,当端点时刻发生改变后,只需将关 稳定供电具有极其重要的意义.引流板螺栓紧固是电 节轨迹与标准化时间的映射关系转换为与物理时间的 力部门带电检修的重要内容之一,传统的螺栓拧紧方 映射就可以得到新的关节轨迹函数,因而减小计算量, 法主要是通过携带工器具登塔完全依靠人工独立 提高算法的实用性.为使关节轨迹运动规划更趋于关 完成,然而高空、高电压环境使得完全依靠人工的方式 节实际运动,解决关节轨迹全局位置约束问题,提出以 不仅劳动强度大和效率低,而且危险性高。基于此,本 轨迹运动时间为优化对象,轨迹最优运动时间所在区 文研制了一种高压线路引流板螺栓紧固带电检修机器 间的求取方法,通过关节运动时间的优化获得满足全 人,通过搭载的双作业臂及其作业末端套筒完成关键 局约束条件的最优关节轨迹,避免超调的发生,提高关 位置的定位,实现引流板螺栓的紧固作业.带电检修 节的运动效率和机器人作业的安全性.通过本文轨迹 机器人在作业过程中,作业臂由初始位姿到工作位姿 规划方法的推广可进一步解决关节速度和加速度状态 各关节位置、速度和加速度状态的瞬间突变对作业臂 约束问题 造成一定冲击,导致机器人振动的发生,给作业臂及其 1 带电检修机器人的实体结构 末端的运动定位精度带来一定影响,因此快速、准确、 稳定的作业臂多关节轨迹规划是保证作业任务顺利完 带电检修机器人由双行走臂、双作业臂及其末端 成的关键 及控制箱等组成,其实体结构如图1所示.双行走臂 机器人关节运动规划方法有A”算法仞、C空间 能够在线路上行走,作业臂1由旋转关节、伸缩关节和 法圆、人工势场法回、MARKOV法@等.在众多方法 纵移关节组成,作业臂2由横移关节、旋转关节、伸缩 中多项式插值算法应用最为普遍,常用的插值多项式 关节和纵移关节组成.双作业臂在线路上作业时都有 有样条函数、三次多项式、五次多项式和高阶多项式. 足够裕度,使其能满足检修作业空间要求.拧螺母作 三次多项式插值算法如简单但不能满足关节加速度 业末端套筒通过燕尾槽固定在作业臂纵移关节上, 的约束要求:高阶多项式插值算法网能够满足关节加 安装于作业臂2的作业末端套筒有大功率电机,可 速度甚至关节脉动的约束要求但算法过于复杂且容易 以实现螺栓紧固,末端采用六角套筒扳手形式以确 产生额外运动:结合机器人作业臂对位置、速度和加速 保套简与螺栓接触面足够大,此外末端内设复位弹 度的控制要求,五次多项式插值算法四及与其相结合 簧的虎克铰,能够适应耐张线夹螺栓一定角度的位 的智能优化算法能够满足作业臂关节运动规划要 置,减小紧固螺栓时对末端产生的内应力.输电导线 求.当前关节轨迹多项式插值规划研究成果在一定程 和耐张引流线通过引流板相联接,其中联接部分通 度上优化了机械臂的运动轨迹,但是存在两个问题:其 过螺栓固定起来.在进行螺栓紧固作业时,先利用机 一,关节轨迹与轨迹端点时刻有关,一旦轨迹端点时刻 械手1末端将螺栓头固定,再由机械手2末端套筒套 发生改变,轨迹函数需要重新计算,从而在一定程度上 住螺母,通过末端工具电机的转向控制实现螺栓的 降低算法的实用性;其二,轨迹规划过于理想化,只考 松开和拧紧操作 虑了轨迹端点的关节状态约束条件,而忽略了中间过 10.作业臂1纵移关节 11.拧螺栓机构 程的关节状态约束,即驱动机构会约束关节在任意时 4.偶臂行走轮 4.奇臂行走轮 刻的运动状态而绝非只是轨迹端点处.通常情况下, 5.夹爪 8.作业臂2 纵移关节 当关节低速运动时只需要保证关节位置不超限,关节 2.作业臂1 及其伸缩关节“ 2.作业营2及其 速度和加速度不易越过约束范围,但当关节高速运动 伸缩关节 1.偶臂 9.作业臂1旋转关节 时,关节速度和加速度很容易越过其约束范围,其结果 机器人机体1.奇臂、 是可能会导致关节电机驱动电流过大损坏机器人硬件 7.作业臂2 系统.因此,在关节轨迹规划时必须考虑各种约束条 旋转关节 6.作业臂2 件,否则可能会造成大的作业事故 横移关节 针对上述问题,为消除关节轨迹端点时刻对关节 轨迹的影响,本文将关节运动时间进行Min一Max标准 3机器人控制箱 化处理,引入标准化时间的概念,以标准化时间为插值 图1带电检修机器人实体结构图 多项式变量,这样作业臂关节运动轨迹仅与轨迹端点 Fig.1 Entity structure of the live maintenance robot
工程科学学报,第 38 卷,第 11 期 KEY WORDS high voltage transmission line; robot; manipulator control; motion planning; optimization 带电检修[1--3]是指在高压电气设备上不停电进行 检修、测试的一种作业方法,其对电网稳定运行,确保 稳定供电具有极其重要的意义. 引流板螺栓紧固是电 力部门带电检修的重要内容之一,传统的螺栓拧紧方 法[4--6]主要是通过携带工器具登塔完全依靠人工独立 完成,然而高空、高电压环境使得完全依靠人工的方式 不仅劳动强度大和效率低,而且危险性高. 基于此,本 文研制了一种高压线路引流板螺栓紧固带电检修机器 人,通过搭载的双作业臂及其作业末端套筒完成关键 位置的定位,实现引流板螺栓的紧固作业. 带电检修 机器人在作业过程中,作业臂由初始位姿到工作位姿 各关节位置、速度和加速度状态的瞬间突变对作业臂 造成一定冲击,导致机器人振动的发生,给作业臂及其 末端的运动定位精度带来一定影响,因此快速、准确、 稳定的作业臂多关节轨迹规划是保证作业任务顺利完 成的关键. 机器人关节运动规划方法有 A* 算法[7]、C 空间 法[8]、人工势场法[9]、MARKOV 法[10]等. 在众多方法 中多项式插值算法应用最为普遍,常用的插值多项式 有样条函数、三次多项式、五次多项式和高阶多项式. 三次多项式插值算法[11]简单但不能满足关节加速度 的约束要求; 高阶多项式插值算法[12]能够满足关节加 速度甚至关节脉动的约束要求但算法过于复杂且容易 产生额外运动; 结合机器人作业臂对位置、速度和加速 度的控制要求,五次多项式插值算法[13]及与其相结合 的智能优化算法[14]能够满足作业臂关节运动规划要 求. 当前关节轨迹多项式插值规划研究成果在一定程 度上优化了机械臂的运动轨迹,但是存在两个问题: 其 一,关节轨迹与轨迹端点时刻有关,一旦轨迹端点时刻 发生改变,轨迹函数需要重新计算,从而在一定程度上 降低算法的实用性; 其二,轨迹规划过于理想化,只考 虑了轨迹端点的关节状态约束条件,而忽略了中间过 程的关节状态约束,即驱动机构会约束关节在任意时 刻的运动状态而绝非只是轨迹端点处. 通常情况下, 当关节低速运动时只需要保证关节位置不超限,关节 速度和加速度不易越过约束范围,但当关节高速运动 时,关节速度和加速度很容易越过其约束范围,其结果 是可能会导致关节电机驱动电流过大损坏机器人硬件 系统. 因此,在关节轨迹规划时必须考虑各种约束条 件[15--16],否则可能会造成大的作业事故. 针对上述问题,为消除关节轨迹端点时刻对关节 轨迹的影响,本文将关节运动时间进行 Min--Max 标准 化处理,引入标准化时间的概念,以标准化时间为插值 多项式变量,这样作业臂关节运动轨迹仅与轨迹端点 状态及关节运动时间有关,而与作业臂关节轨迹端点 时刻没有直接联系,当端点时刻发生改变后,只需将关 节轨迹与标准化时间的映射关系转换为与物理时间的 映射就可以得到新的关节轨迹函数,因而减小计算量, 提高算法的实用性. 为使关节轨迹运动规划更趋于关 节实际运动,解决关节轨迹全局位置约束问题,提出以 轨迹运动时间为优化对象,轨迹最优运动时间所在区 间的求取方法,通过关节运动时间的优化获得满足全 局约束条件的最优关节轨迹,避免超调的发生,提高关 节的运动效率和机器人作业的安全性. 通过本文轨迹 规划方法的推广可进一步解决关节速度和加速度状态 约束问题. 1 带电检修机器人的实体结构 带电检修机器人由双行走臂、双作业臂及其末端 及控制箱等组成,其实体结构如图 1 所示. 双行走臂 能够在线路上行走,作业臂 1 由旋转关节、伸缩关节和 纵移关节组成,作业臂 2 由横移关节、旋转关节、伸缩 关节和纵移关节组成. 双作业臂在线路上作业时都有 足够裕度,使其能满足检修作业空间要求. 拧螺母作 业末端套筒通过燕尾槽固定在作业臂纵移关节上, 安装于作业臂 2 的作业末端套筒有大功率电机,可 以实现螺栓紧固,末端采用六角套筒扳手形式以确 保套筒与螺栓接触面足够大,此外末端内设复位弹 簧的虎克铰,能够适应耐张线夹螺栓一定角度的位 置,减小紧固螺栓时对末端产生的内应力. 输电导线 和耐张引流线通过引流板相联接,其中联接部分通 过螺栓固定起来. 在进行螺栓紧固作业时,先利用机 械手 1 末端将螺栓头固定,再由机械手 2 末端套筒套 住螺母,通过末端工具电机的转向控制实现螺栓的 松开和拧紧操作. 图 1 带电检修机器人实体结构图 Fig. 1 Entity structure of the live maintenance robot ·1644·
江维等:高压输电线路带电检修机器人作业臂运动优化 ·1645· 2改进N次多项式插值轨迹规划算法的基 i(.)=-Rnd.), ti-to 本原理 2.1常规N次多项式轨迹规划算法 x(t.)=7 -o-Rand(t), 4-o)2 常规多项式插值关节轨迹运动规划方法是通过轨 (5) 迹端点处关节的边界状态来确定唯一的一个插值多项 w)--合d. 式.设起始和终点时刻分别为和,起点的状态分别 为o,o,。,…,终点的状态分别为x,,元,…,插值 函数为x(),则x(t)满足式(1)的边界条件. -. x(la)=xo, n阶多项式插值轨迹需要n+1个初始条件,由于 x()=x 插值端点初始状态的对称性n一般为奇数,若定义关 (to)=o, 节轨迹的运行时间为lam=4-o,轨迹端点位置差 x(t)=, S=x:-x。,将关节变量两个端点的边界条件分别代入 式(4)和式(5)可得式(6). 花()=元o, (1) [x(to)=xo+S Rand(0)=x, (t)=, x(t)=xo+S Rand (1)=x x0()=xg0. i()=5-Raind(O) ()=,S-Rand()=, 由式(1)n+1个边界条件可确定唯一的一个关 节轨迹n(n为奇数)次多项式: x()=a.+a-1-+a-2-2+…+a,2+a,t+ao 上E(-之)=o: (6) (2) 将式(1)n+1个边界条件代入式(2)可得到n+ )-之Rmdm- 1元一次方程组,所以可解得插值多项式系数a,~an, 从而获取关节的运动轨迹多项式 2.2改进W次多项式轨迹规划算法 学)w=之a)(0=), 常规N次多项式插值轨迹是时间的函数.对于 x()(e)=sRad(学)()=x学) 不同的端点时刻有不同的边界状态,插值多项式系数 需要重新计算.为淡化端点时刻对插值多项式的影 由式(6)可得Rand(t.)函数的n+1个边界条件 响,减小插值系数的计算量,对时间变量t进行Min一 式(7). Mx标准化处理,将物理时间转化为一个无量纲的标 [Rand (0)=0, 准化时间变量t。如式(3),则关节轨迹函数可表示为 Rand(1)=1, 式(4),再根据轨迹端点处的关节状态值可求得关节 Raind (0)=- om 关于t。的多项式轨迹 t-o,0≤t.≤1. ta=I-to (3) Rand (1)= x(l.)=x。+(x-xo)·Rand(t.),0≤Rand(t.)≤l. Rand (0)= , (4) (7) 式(4)中Rand(t,)为D,l]之间的随机函数,可对 Band (1) 关节轨迹函数进行实时修正.当t=l。时Rand(t.)= S 0,此时表示轨迹起点;当t=t:时Rand(l.)=l,此时表 示轨迹终点:当0<Rand(t)<1时表示关节轨迹运动 Rm(学)O0-学g 过程中的任意一个中间点.对式(4)分别求一阶导数、 二阶导数、三阶导数至子(a-1)阶导数,可得关节轨 Band((1) 迹速度、加速度函数、加加速度等状态如式(5)所示
江 维等: 高压输电线路带电检修机器人作业臂运动优化 2 改进 N 次多项式插值轨迹规划算法的基 本原理 2. 1 常规 N 次多项式轨迹规划算法 常规多项式插值关节轨迹运动规划方法是通过轨 迹端点处关节的边界状态来确定唯一的一个插值多项 式. 设起始和终点时刻分别为 t0和 tf,起点的状态分别 为 x0,x · 0,x ·· 0,…,终点的状态分别为 xf,x · f,x ·· f,…,插值 函数为 x( t) ,则 x( t) 满足式( 1) 的边界条件. x( t0 ) = x0, x( tf ) = xf, x ·( t0 ) = x · 0, x ·( tf ) = x · f, x ··( t0 ) = x ·· 0, x ··( tf ) = x ·· f, … x( n) ( t0 ) = x( n) 0 . … ( 1) 由式( 1) n + 1 个边界条件可确定唯一的一个关 节轨迹 n ( n 为奇数) 次多项式: x( t) = an t n + an - 1 t n - 1 + an - 2 t n - 2 + … + a2 t 2 + a1 t + a0 . ( 2) 将式( 1) n + 1 个边界条件代入式( 2) 可得到 n + 1 元一次方程组,所以可解得插值多项式系数 a0 ~ an, 从而获取关节的运动轨迹多项式. 2. 2 改进 N 次多项式轨迹规划算法 常规 N 次多项式插值轨迹是时间 t 的函数. 对于 不同的端点时刻有不同的边界状态,插值多项式系数 需要重新计算. 为淡化端点时刻对插值多项式的影 响,减小插值系数的计算量,对时间变量 t 进行 Min-- Max 标准化处理,将物理时间转化为一个无量纲的标 准化时间变量 tα 如式( 3) ,则关节轨迹函数可表示为 式( 4) ,再根据轨迹端点处的关节状态值可求得关节 关于 tα 的多项式轨迹. tα = t - t0 tf - t0 ,0≤tα≤1. ( 3) x( tα ) = x0 + ( xf - x0 )·Rand( tα ) ,0≤Rand( tα ) ≤1. ( 4) 式( 4) 中 Rand( tα ) 为[0,1]之间的随机函数,可对 关节轨迹函数进行实时修正. 当 t = t0 时 Rand( tα ) = 0,此时表示轨迹起点; 当 t = tf 时 Rand( tα ) = 1,此时表 示轨迹终点; 当 0 < Rand( tα ) < 1 时表示关节轨迹运动 过程中的任意一个中间点. 对式( 4) 分别求一阶导数、 二阶导数、三阶导数至 1 2 ( n - 1) 阶导数,可得关节轨 迹速度、加速度函数、加加速度等状态如式( 5) 所示. x ·( tα ) = xf - x0 tf - t0 Rand · ( tα ) , x ··( tα ) = xf - x0 ( tf - t0 ) 2Rand ·· ( tα ) , x ···( tα ) = xf - x0 ( tf - t0 ) 3Rand ·· ( tα ) , … ( x n - 1 2 ) ( tα ) = xf - x0 ( tf - t0 ) n ( Rand n - 1 2 ) ( tα ) . ( 5) n 阶多项式插值轨迹需要 n + 1 个初始条件,由于 插值端点初始状态的对称性 n 一般为奇数,若定义关 节轨迹的运行时间为 tmotion = tf - t0,轨迹端点位置差 S = xf - x0,将关节变量两个端点的边界条件分别代入 式( 4) 和式( 5) 可得式( 6) . x( t0 ) = x0 + S Rand( 0) = x0, x( tf ) = x0 + S Rand( 1) = xf, x ·( t0 ) = S tmotion Rand · ( 0) = x · 0, x ·( tf ) = S tmotion Rand · ( 1) = x · f, x ··( t0 ) = S t 2 motion Rand ·· ( 0) = x ·· 0, x ··( tf ) = S t 2 motion Rand ·· ( 1) = x ·· f, … ( x n - 1 2 ) ( t0 ) = S t n motion ( Rand n - 1 2 ) ( 0) = x ( 0 n - 1 2 ) , ( x n - 1 2 ) ( tf ) = S t n motion ( Rand n - 1 2 ) ( 1) = x( f n - 1 2 ) . ( 6) 由式( 6) 可得 Rand( tα ) 函数的 n + 1 个边界条件 式( 7) . Rand( 0) = 0, Rand( 1) = 1, Rand · ( 0) = x · 0 tmotion S , Rand · ( 1) = x · f tmotion S , Rand ·· ( 0) = x ·· 0 t 2 motion S , Rand ·· ( 1) = x ·· f t 2 motion S , … ( Rand n - 1 2 ) ( 0) = x ( 0 n - 1 2 ) t n motion S , ( Rand n - 1 2 ) ( 1) = x( f n - 1 2 ) t n motion S . ( 7) ·1645·
·1646 工程科学学报,第38卷,第11期 由式(7)的n+1个边界条件可唯一确定Rand(.)的 将式(7)边界条件分别代入式(8)和式(9)可得式 一个n次多项式: (10) Rand(t.)=An(t.)"+An-i(t.)"-1+…+ Az (ta)2+At Ao. (8) 「Ao=Rand(o), 对式(8)求一阶导数、二阶导数至方(a-1)阶导 An+An-1+…+A1+A。=Rand(1), A.Raind (0), 数可得式(9). nA。+(n-1)An-1+…+2A2+A1=Rnd(1), [Rnd(l.)=nA.。+(n-1)An-ia-2+…+ 2A2。+A, Rand(t.)=n(n-1)A.&-2+ (2)1A(学=Ran()o, (n-1)(n-2)An--3+…+2A2, nm-)…(a-2)A+ Rand(学)(e.)= m-0a-2)…(n-2)4+…+ nm-)…(a-"2)A)(学+ (2)!A(学)=Rad()(). a-Da-2(a-” (10) A.)()+…+(受)!A(学) 为便于算法的推广及MATLAB计算将式(1O)改 写为矩阵形式: (9) 0 0 0 1 1 0 0 … … … 1 0 n-1 … … 2 1 … … … … … n(m-)…(n-"2)(a-a-2)…(n-2) (n+)×(n+) Rand (0) A.-1 Rand(1) A-2 Rand (0) Rand(1) … A(学) (11) Rand(学)(O) Rand(学)(1) A Rand(学)(O) LRand(学)(a).J 所以
工程科学学报,第 38 卷,第 11 期 由式( 7) 的 n + 1 个边界条件可唯一确定 Rand( tα ) 的 一个 n 次多项式: Rand( tα ) = An ( tα ) n + An - 1 ( tα ) n - 1 + … + A2 ( tα ) 2 + A1 tα + A0 . ( 8) 对式( 8) 求一阶导数、二阶导数至 1 2 ( n - 1) 阶导 数可得式( 9) . Rand · ( tα ) = nAn t n - 1 α + ( n - 1) An - 1 t n - 2 α + … + 2A2 tα + A1, Rand ·· ( tα ) = n( n - 1) An t n - 2 α + ( n - 1) ( n - 2) An - 1 t n - 3 α + … + 2A2, … ( Rand n - 1 2 ) ( tα ) = n( n - 1 ( ) … n - n - 3 ) 2 An ( tα ( ) n + 1 2 ) + ( n - 1) ( n - 2 ( ) … n - n - 1 ) 2 · An - 1 ( tα ( ) n - 1 2 ) + … ( + n - 1 ) 2 ! ( A n - 1 2 ) . ( 9) 将式( 7) 边界条件分别代入式( 8) 和式( 9) 可得式 ( 10) . A0 = Rand( 0) , An + An - 1 + … + A1 + A0 = Rand( 1) , A1 = Rand · ( 0) , nAn + ( n - 1) An - 1 + … + 2A2 + A1 = Rand · ( 1) , ( … n - 1 ) 2 ! ( A n - 1 2 ) ( = Rand n - 1 2 ) ( 0) , n( n - 1 ( ) … n - n - 3 ) 2 An + ( n - 1) ( n - 2 ( ) … n - n - 1 ) 2 An - 1 + … + ( n - 1 ) 2 ! ( A n - 1 2 ) ( = Rand n - 1 2 ) ( 1) . ( 10) 为便于算法的推广及 MATLAB 计算将式( 10) 改 写为矩阵形式: 0 0 … … … … … 0 1 1 1 … … … … … 1 1 0 0 … … … … … 1 0 n n - 1 … … … … … 2 1 … … … … … … … … … … … … … … … … … … … … … … … … … … … … … … … … … … … … n( n - 1 ( ) … n - n - 3 ) 2 ( n - 1) ( n - 2 ( ) … n - n - 1 ) 2 … … … ( n - 1 ) 2 ! … … … ( n + 1) × ( n + 1) · An An - 1 An - 2 … ( A n - 1 2 ) … A2 A1 A 0 = Rand( 0) Rand( 1) Rand · ( 0) Rand · ( 1) … ( Rand n - 3 2 ) ( 0) ( Rand n - 3 2 ) ( 1) ( Rand n - 1 2 ) ( 0) ( Rand n - 1 2 ) ( 1) . . ( 11) 所以, ·1646·
江维等:高压输电线路带电检修机器人作业臂运动优化 ·1647· A。 0 0 1 1 1 1 An-2 0 0 ……… 1 0 n n-1 A(学) … … … …… … A ……… …… A a-a-2)…(a-"号)》 《m+1》x(a+1) Rand (0) Rand (1) Rand (0) Rand (1) (12) Rand(学)(O) Rand(学)(a) Rand(学)(o) LRand()(1). 通过式(12)MATLAB矩阵运算即可解得 [4=6+-3i-3-05,2+0564, Rand(t.)的系数,从而可以获取关于t。的关节n次多 项式插值轨迹.由于对于一般的关节而言满足位置、 A4=-15+ 8e+7在m+l.5,2。- 速度和加速度的性能要求就可以达到较好的轨迹规划 效果,因此一般情况下采用五次多项式插值算法即可 4,=10+-6,-4-15元+05, 满足性能要求.本文以五次多项式插值为例可求得轨 迹多项式的系数为式(13). 0.5 A, TAs =-6Rand (0)+6Rand (1)-3 Rand (0)- 3Rand(1)-0.5Rand(0)+0.5Rand(1), A=S A =15Rand (0)-15Rand(1)+8 Raind (0)+ 1A=0. (14) 7 Rand (1)+1.5 Rand (0)-Rand (1), 将求得的系数式(14)代入式(8)可得Rand(t.)的 A3 =-10Rand (0)+10Rand (1)-6 Raind (0)- 五次多项式为式(15),将式(15)代入式(4)可得基于 4 Rand(1)-1.5 Rand (0)+0.5 Rand(1), 标准化时间的关节轨迹函数x(1)为式(16) A2=0.5Rand(0), Rand(t)= A =Raind (0), (6+-3-3-0+0)+ A。=Rand(0). (13) (-15+8+75-)+ 将式(7)边界条件代入式(13)得式(14). (0+-6-4-15+05)+ (15)
江 维等: 高压输电线路带电检修机器人作业臂运动优化 An An - 1 An - 2 … ( A n - 1 2 ) … A2 A1 A 0 = 0 0 … … … … … 0 1 1 1 … … … … … 1 1 0 0 … … … … … 1 0 n n - 1 … … … … … 2 1 … … … … … … … … … … … … … … … … … … … … … … … … … … … … … … … … … … … … n( n - 1 ( ) … n - n - 3 ) 2 ( n - 1) ( n - 2 ( ) … n - n - 1 ) 2 ( ) … … … n - 1 ) 2 ! … … … - 1 ( n + 1) × ( n + 1) · Rand( 0) Rand( 1) Rand · ( 0) Rand · ( 1) … ( Rand n - 3 2 ) ( 0) ( Rand n - 3 2 ) ( 1) ( Rand n - 1 2 ) ( 0) ( Rand n - 1 2 ) ( 1) . . ( 12) 通 过 式 ( 12 ) MATLAB 矩 阵 运 算 即 可 解 得 Rand( tα ) 的系数,从而可以获取关于 tα 的关节 n 次多 项式插值轨迹. 由于对于一般的关节而言满足位置、 速度和加速度的性能要求就可以达到较好的轨迹规划 效果,因此一般情况下采用五次多项式插值算法即可 满足性能要求. 本文以五次多项式插值为例可求得轨 迹多项式的系数为式( 13) . A5 = - 6Rand( 0) + 6Rand( 1) - 3 Rand · ( 0) - 3 Rand · ( 1) - 0. 5 Rand ·· ( 0) + 0. 5 Rand ·· ( 1) , A4 = 15Rand( 0) - 15Rand( 1) + 8 Rand · ( 0) + 7 Rand · ( 1) + 1. 5 Rand ·· ( 0) - Rand ·· ( 1) , A3 = - 10Rand( 0) + 10Rand( 1) - 6 Rand · ( 0) - 4 Rand · ( 1) - 1. 5 Rand ·· ( 0) + 0. 5 Rand ·· ( 1) , A2 = 0. 5 Rand ·· ( 0) , A1 = Rand · ( 0) , A0 = Rand( 0) . ( 13) 将式( 7) 边界条件代入式( 13) 得式( 14) . A5 =6 + -3x · 0 tmotion -3x · f tmotion -0. 5x ·· 0 t 2 motion +0. 5x ·· f t 2 motion S , A4 = - 15 + 8x · 0 tmotion + 7x · f tmotion + 1. 5x ·· 0 t 2 motion - x ·· f t 2 motion S , A3 =10 + -6x · 0 tmotion -4x · f tmotion -1. 5x ·· 0 t 2 motion +0. 5x ·· f t 2 motion S , A2 = 0. 5x ·· 0 t 2 motion S , A1 = x · 0 tmotion S , A0 = 0 . ( 14) 将求得的系数式( 14) 代入式( 8) 可得 Rand( tα ) 的 五次多项式为式( 15) ,将式( 15) 代入式( 4) 可得基于 标准化时间的关节轨迹函数 x( tα ) 为式( 16) . Rand( tα ) ( = 6 + -3x · 0 tmotion -3x · f tmotion -0. 5x ·· 0 t 2 motion +0. 5x ·· f t 2 motion ) S t 5 α ( + - 15 + 8x · 0 tmotion + 7x · f tmotion + 1. 5x ·· 0 t 2 motion - x ·· f t 2 motion ) S t 4 α ( + 10 + -6x · 0 tmotion -4x · f tmotion -1. 5x ·· 0 t 2 motion +0. 5x ·· f t 2 motion ) S t 3 α + 0. 5x ·· 0 t 2 motion S t 2 α + x · 0 tmotion S tα . ( 15) ·1647·
·1648· 工程科学学报,第38卷,第11期 x(l.)=(6S-3 o-3tmim-0.5元2ian+ 动到作业位姿,实现双作业手末端与引流板螺栓螺 0.5元r2ia)t2+(-15S+8元+7元liam+ 母的对准对接.为便于分析以定位到压接管的机器 L.5tot2im-花2ia)。+(10S+6 xoli- 人奇臂行走轮为空间坐标原点,沿导线方向为X轴、 4i1ain-l.5元of2in+0.5元f2ii)2+ 作业臂伸缩方向为Y轴、作业手纵移方向为Z轴建 0.3min+ (16) 立空间坐标系.设机器人作业手末端在空间坐标系 通过对式(16)求导可得关节轨迹速度与加速度 中的坐标为Po:(xo,ya,2)(i=1时表示作业手1, 的表达式.上述各式为轨迹的标准化时间函数,若要 i=2时表示作业手2),将P定义为机器人双作业手 的初始位置如图2(a)所示.设引流板螺栓螺母在空 转换为物理时间函数可将1。=仁(0≤1.≤1)代入 to-to 间坐标系中的坐标为Pa=(xya,z)(i=1时表示 各式即可. 螺栓,i=2时表示螺母),将P定义为机器人双作业 手的终止位姿如图2(b)所示.机器人作业手末端由 3带电检修机器人作业臂多关节轨迹优化 初始位姿到螺栓螺母的对准状态整个过程可抽象为 3.1双作业臂运动模型的建立 作业臂各关节轨迹规划控制问题,其实质是作业手 机器人双作业臂位姿的规划是依靠各关节的运动 在空间坐标系中由初始位姿P通过作业臂的旋转、 来实现的,因此必须将双作业臂轨迹转换到关节空间, 伸缩和横移,作业手的纵移关节运动调整双作业手 再在关节空间中进行轨迹规划.通过机器人双作业 在空间坐标系中的位姿,不断逼近理想位姿P:的 臂各关节的轨迹规划,使得双作业臂由初始位姿运 过程. (o Yor=n) 0.0.0 婚 螺母 (x- 螺栓头 导轨 耐张跳线 耐张跳线 带电检修机器人 带电检修机器人 图2机器人双作业臂运动规划.()初始位姿:(b)理想位姿 Fig.2 Motion planning of the robot's double manipulators:(a)initial pose:(b)ideal pose 3.2基于运动时间最优的关节轨迹优化 终点的状态及轨迹的运动时间有关,而与端点时刻无 理论上讲作业臂关节轨迹状态包括速度和加速度 关,因此可以考虑通过优化关节的运动时间,使得关节 都会由于机械设计及关节电机的因素而限定在一定范 速度轨迹在运动时间内呈单调递增或递减状态,以保 围之内.作业臂关节速度轨迹应该在关节电机输出的 证关节位置轨迹满足其约束条件.即只需要保证关节 最大转速与最小转速之间、关节加速度轨迹会受到关 速度在运动时间内恒为正或者恒为负,当速度大于零 节电机驱动力矩的限制,也应该在某一个区间范围内,但 时位置轨迹为单调递增状态,当速度小于零时位置轨 关节速度轨迹和加速度轨迹的制约均可以通过关节电机 迹为单调递减状态.因此保证关节位置轨迹满足其约 的选型解决,通过选择合适参数的电机来克服关节速度 束范围的条件为 和加速度的制约条件.因此关于作业臂轨迹优化,本文主 x(t)≥0或x(t)≤0 Rand(tn)≥0或Rnd(tn)≤0. 要研究的是关节位置约束优化问题,通过本文方法的推 (17) 广也可进一步研究关节速度和加速度约束优化问题 对式(15)求导可得 采用多项式插值法对作业臂进行运动规划,任意 关节轨迹的端点位置需要满足关节位置的约束要求, Rmd)=(0-15、玉.15ti_25三+ 在整个运动过程中如果关节位置从起点到终点呈单调 递增或递减运动状态,那么在满足这种条件状态下时 25+(-0与,, 可保证关节位置满足其运动约束条件.由改进多项式 30- 18_12主 插值轨迹规划算法可知,关节运动轨迹仅仅与起点和
工程科学学报,第 38 卷,第 11 期 x( tα ) = ( 6S - 3x · 0 tmotion - 3x · f tmotion - 0. 5x ·· 0 t 2 motion + 0. 5x ·· f t 2 motion ) t 5 α + ( - 15S + 8x · 0 tmotion + 7x · f tmotion + 1. 5x ·· 0 t 2 motion - x ·· f t 2 motion ) t 4 α + ( 10S + 6x · 0 tmotion - 4x · f tmotion - 1. 5x ·· 0 t 2 motion + 0. 5x ·· f t 2 motion ) t 3 α + 0. 5x ·· 0 t 2 motion t 2 α + 3x · 0 tmotion tα + x0 . ( 16) 通过对式( 16) 求导可得关节轨迹速度与加速度 的表达式. 上述各式为轨迹的标准化时间函数,若要 转换为物理时间函数可将 tα = t - t0 tf - t0 ( 0≤tα≤1) 代入 各式即可. 3 带电检修机器人作业臂多关节轨迹优化 3. 1 双作业臂运动模型的建立 机器人双作业臂位姿的规划是依靠各关节的运动 来实现的,因此必须将双作业臂轨迹转换到关节空间, 再在关节空间中进行轨迹规划. 通过机器人双作业 臂各关节的轨迹规划,使得双作业臂由初始位姿运 动到作业位姿,实现双作业手末端与引流板螺栓螺 母的对准对接. 为便于分析以定位到压接管的机器 人奇臂行走轮为空间坐标原点,沿导线方向为 X 轴、 作业臂伸缩方向为 Y 轴、作业手纵移方向为 Z 轴建 立空间坐标系. 设机器人作业手末端在空间坐标系 中的坐标为 P0i ( x0i,y0i,z0i ) ( i = 1 时表示作业手 1, i = 2 时表示作业手 2) ,将 P0i定义为机器人双作业手 的初始位置如图 2( a) 所示. 设引流板螺栓螺母在空 间坐标系中的坐标为 Pfi = ( xfi,yfi,zfi ) ( i = 1 时表示 螺栓,i = 2 时表示螺母) ,将 Pfi定义为机器人双作业 手的终止位姿如图 2( b) 所示. 机器人作业手末端由 初始位姿到螺栓螺母的对准状态整个过程可抽象为 作业臂各关节轨迹规划控制问题,其实质是作业手 在空间坐标系中由初始位姿 P0i通过作业臂的旋转、 伸缩和横移,作业手的纵移关节运动调整双作业手 在空间坐标系中的位姿,不断逼近理想位姿 Pfi 的 过程. 图 2 机器人双作业臂运动规划. ( a) 初始位姿; ( b) 理想位姿 Fig. 2 Motion planning of the robot's double manipulators: ( a) initial pose; ( b) ideal pose 3. 2 基于运动时间最优的关节轨迹优化 理论上讲作业臂关节轨迹状态包括速度和加速度 都会由于机械设计及关节电机的因素而限定在一定范 围之内. 作业臂关节速度轨迹应该在关节电机输出的 最大转速与最小转速之间、关节加速度轨迹会受到关 节电机驱动力矩的限制,也应该在某一个区间范围内,但 关节速度轨迹和加速度轨迹的制约均可以通过关节电机 的选型解决,通过选择合适参数的电机来克服关节速度 和加速度的制约条件. 因此关于作业臂轨迹优化,本文主 要研究的是关节位置约束优化问题,通过本文方法的推 广也可进一步研究关节速度和加速度约束优化问题. 采用多项式插值法对作业臂进行运动规划,任意 关节轨迹的端点位置需要满足关节位置的约束要求, 在整个运动过程中如果关节位置从起点到终点呈单调 递增或递减运动状态,那么在满足这种条件状态下时 可保证关节位置满足其运动约束条件. 由改进多项式 插值轨迹规划算法可知,关节运动轨迹仅仅与起点和 终点的状态及轨迹的运动时间有关,而与端点时刻无 关,因此可以考虑通过优化关节的运动时间,使得关节 速度轨迹在运动时间内呈单调递增或递减状态,以保 证关节位置轨迹满足其约束条件. 即只需要保证关节 速度在运动时间内恒为正或者恒为负,当速度大于零 时位置轨迹为单调递增状态,当速度小于零时位置轨 迹为单调递减状态. 因此保证关节位置轨迹满足其约 束范围的条件为 x ·( tα ) ≥0 或 x ·( tα ) ≤0Rand · ( tα ) ≥0 或Rand · ( tα ) ≤0. ( 17) 对式( 15) 求导可得 Rand · ( tα ) ( = 30 - 15tmotion x · 0 S - 15tmotion x · f S - 2. 5t 2 motion x ·· 0 S + 2. 5t 2 motion x ·· f ) S t 4 α ( + - 60 + 32tmotion x · 0 S + 28tmotion x · f S + 6t 2 motion x ·· 0 S - 4t 2 motion x ·· f ) S t 3 α ( + 30 - 18tmotion x · 0 S - 12tmotion x · f S - ·1648·
江维等:高压输电线路带电检修机器人作业臂运动优化 ·1649· i伯(amia≤tin(0≤i()(ma),i=0,l,…,7. (20) Max(()≤tot≤Min(mn(),i=l,2,…,7. (18) (21) 对于式(18)Rnd(t,)为一个关于t。的四次多项 4 式,由于0≤t。≤l,为保证Rnd(t.)在区间D,l]内恒 仿真实验与现场运行 为正或恒为负,则其零点必须在D,1]之外,由于、x 4.1仿真实验 和S均为已知的关节状态常量,因此当x,>x。位置函 设关节起点状态为xo=0.5rad、i。=l.5rads和 数为增函数时可令式(18)大于0求解出1i的范围, x。=0.52rads2,终点状态为x,=5.5rad、t,=2.3rad 当x,<x。位置函数为减函数时可令式(18)小于0求 s和=1.22rads2,端点时刻为。=0s,l,=5s],常 解出L的范围,在t的范围内选取其最优值以获 规算法与改进算法的轨迹运动仿真结果对比如图3所 得满足位置约束条件的关节最优轨迹.若将上述两种 示.其中常规五次多项式插值关节轨迹规划曲线如 情况归纳为一种统一的形式,可得关节轨迹运动最优 图3()所示:将端点时刻映射到D,1]区间范围后,改 时间所满足的条件为式(19),结合式(18)可将式(19) 进五次多项式插值关节轨迹规划曲线如图3(b)所示. 并代入边界条件参数可得关于t的不等式,同样易 保持原有的关节状态不变,在保证关节运动时间不变 求得最优ti的取值范围. 的前提下选取端点时刻为。=1.5s,=6.5s],利用 Rand(0)×Rand(1)≥0. (19) 改进算法获得的关节轨迹如图4(b)所示.由仿真结 3.3双作业臂多关节轨迹优化 果可知改进算法和常规算法轨迹规划曲线完全相同, 带电检修机器人的双作业臂由旋转、伸缩、横移和 仅仅只是坐标轴横轴的刻度不同,因此改进算法完全 纵移共七个关节构成,通过各关节的轨迹运动规划使 可以实现常规算法的功能。另外,改进算法在选取不 得双作业臂及其末端由初始位姿运动到螺栓螺母对准 同端点时刻时轨迹规划曲线也是完全相同的,验证了 状态,为保证各关节位置轨迹满足其约束条件,利用关 改进算法只与端点状态及执行时间有关,而与端点时 节运动时间最优轨迹规划原理对各关节的运动时间进 刻无关,当端点时刻发生改变后只需将关节轨迹与标 行优化处理以获得各关节最优运动轨迹.为提高作业 准化时间的映射关系转换为与物理时间的映射就可以 效率,各关节为联动模式,由于作业内容和对象已经确 得到新的关节轨迹函数,提高了算法的实用性. 定,要求各关节同时到达目标点,而各关节的起始时间 对于原有的关节状态,通过人工设置关节运动时 时刻是相同的,所以对于双作业臂各关节其运动时间 间和经过优化关节运动时间仿真结果对比如图5所 是相同的,即对七个关节设定一个相同的运动时间. 示,其中人工设定关节运动时间tm=1s.利用改进 然而由于机械特性及关节电机的驱动能力等因素,各 算法关节轨迹如图5(a)所示:利用优化后的执行时间 关节各有一个独立的最优执行时间范围,因此可对各 【m=0.9s改进算法获得的关节轨迹如图5(b)所示. 关节独立最优执行时间取交集以获得各关节整体最优 通过仿真结果可知:当人工设定执行时间时轨迹出现 执行时间,设各关节独立执行时间约束范围为式 超调,当轨迹运行时间为0.72s时,关节轨迹已经运动 (20),取各关节独立执行时间的交集可得各关节整体 到达指定位置,这样会造成在剩余时间里关节产生额外 最优执行时间la应满足式(21).根据各关节的初始 运动,甚至位置超出其约束范围:而经过算法优化得到轨 条件以及l可得双臂各关节满足位置约束的最优运 迹运行时间,关节轨迹在li=0.9s时刻恰好到达指定 动轨迹规划. 位置,避免了超调的发生,提高了关节的运动效率 6 a 6而 5 5 4 2 0.51.01.5 2.02.53.03.54.04.55.0 0 0.10.2 0.30.4050.60.70.80.91.0 时间/s 标准化时间 图3常规算法(a)和改进算法(b)关节轨迹曲线 Fig.3 Joint trajectory curves of the conventional algorithm (a)and the improved algorithm (b)
江 维等: 高压输电线路带电检修机器人作业臂运动优化 4. 5t 2 motion x ·· 0 S + 1. 5t 2 motion x ·· f ) S t 2 α + t 2 motion x ·· 0 S tα + tmotion x · 0 S . ( 18) 对于式( 18) Rand · ( tα ) 为一个关于 tα 的四次多项 式,由于 0≤tα≤1,为保证Rand · ( tα ) 在区间[0,1]内恒 为正或恒为负,则其零点必须在[0,1]之外,由于 x · 0、x · f 和 S 均为已知的关节状态常量,因此当 xf > x0 位置函 数为增函数时可令式( 18) 大于 0 求解出 tmotion的范围, 当 xf < x0 位置函数为减函数时可令式( 18) 小于 0 求 解出 tmotion的范围,在 tmotion的范围内选取其最优值以获 得满足位置约束条件的关节最优轨迹. 若将上述两种 情况归纳为一种统一的形式,可得关节轨迹运动最优 时间所满足的条件为式( 19) ,结合式( 18) 可将式( 19) 并代入边界条件参数可得关于 tmotion的不等式,同样易 求得最优 tmotion的取值范围. Rand · ( 0) × Rand · ( 1) ≥0. ( 19) 图 3 常规算法( a) 和改进算法( b) 关节轨迹曲线 Fig. 3 Joint trajectory curves of the conventional algorithm ( a) and the improved algorithm ( b) 3. 3 双作业臂多关节轨迹优化 带电检修机器人的双作业臂由旋转、伸缩、横移和 纵移共七个关节构成,通过各关节的轨迹运动规划使 得双作业臂及其末端由初始位姿运动到螺栓螺母对准 状态,为保证各关节位置轨迹满足其约束条件,利用关 节运动时间最优轨迹规划原理对各关节的运动时间进 行优化处理以获得各关节最优运动轨迹. 为提高作业 效率,各关节为联动模式,由于作业内容和对象已经确 定,要求各关节同时到达目标点,而各关节的起始时间 时刻是相同的,所以对于双作业臂各关节其运动时间 是相同的,即对七个关节设定一个相同的运动时间. 然而由于机械特性及关节电机的驱动能力等因素,各 关节各有一个独立的最优执行时间范围,因此可对各 关节独立最优执行时间取交集以获得各关节整体最优 执行时 间,设 各 关 节 独 立 执 行 时 间 约 束 范 围 为 式 ( 20) ,取各关节独立执行时间的交集可得各关节整体 最优执行时间 tbest应满足式( 21) . 根据各关节的初始 条件以及 tbest可得双臂各关节满足位置约束的最优运 动轨迹规划. tmotion( i) ( min) ≤tmotion( i) ≤tmotion( i) ( max) ,i = 0,1,…,7. ( 20) Max( tmotion( i) ( min) ) ≤tbest≤Min( tmotion( i) ( max) ) ,i = 1,2,…,7. ( 21) 4 仿真实验与现场运行 4. 1 仿真实验 设关节起点状态为 x0 = 0. 5 rad、x · 0 = 1. 5 rad·s - 1 和 x ·· 0 = 0. 52 rad·s - 2 ,终点状态为 xf = 5. 5 rad、x · f = 2. 3 rad· s 和 x ·· f = 1. 22 rad·s - 2 ,端点时刻为[t0 = 0 s,tf = 5 s],常 规算法与改进算法的轨迹运动仿真结果对比如图 3 所 示. 其中常规五次多项式插值关节轨迹规划曲线如 图 3( a) 所示; 将端点时刻映射到[0,1]区间范围后,改 进五次多项式插值关节轨迹规划曲线如图 3( b) 所示. 保持原有的关节状态不变,在保证关节运动时间不变 的前提下选取端点时刻为[t0 = 1. 5 s,tf = 6. 5 s],利用 改进算法获得的关节轨迹如图 4( b) 所示. 由仿真结 果可知改进算法和常规算法轨迹规划曲线完全相同, 仅仅只是坐标轴横轴的刻度不同,因此改进算法完全 可以实现常规算法的功能. 另外,改进算法在选取不 同端点时刻时轨迹规划曲线也是完全相同的,验证了 改进算法只与端点状态及执行时间有关,而与端点时 刻无关,当端点时刻发生改变后只需将关节轨迹与标 准化时间的映射关系转换为与物理时间的映射就可以 得到新的关节轨迹函数,提高了算法的实用性. 对于原有的关节状态,通过人工设置关节运动时 间和经过优化关节运动时间仿真结果对比如图 5 所 示,其中人工设定关节运动时间 tmotion = 1 s. 利用改进 算法关节轨迹如图 5( a) 所示; 利用优化后的执行时间 tmotion = 0. 9 s 改进算法获得的关节轨迹如图 5( b) 所示. 通过仿真结果可知: 当人工设定执行时间时轨迹出现 超调,当轨迹运行时间为 0. 72 s 时,关节轨迹已经运动 到达指定位置,这样会造成在剩余时间里关节产生额外 运动,甚至位置超出其约束范围; 而经过算法优化得到轨 迹运行时间,关节轨迹在 tmotion = 0. 9 s 时刻恰好到达指定 位置,避免了超调的发生,提高了关节的运动效率. ·1649·
·1650 工程科学学报,第38卷,第11期 6 6 (a) (b) 5 5 4 3 平 2 1 0.10.2 0.30.40.50.60.70.80.9 1.0 0.1 0.203 0.4050.6 0.70.80.91.0 标准化时间 标准化时间 图4不同端点时刻关节轨迹曲线.(a)0=0,4=5s:(b)0=1.5s,4=6.5s Fig.4 Joint trajectory curves at different endpoint times:(a)to =0,=5s:(b)to =1.5s,t =6.5s 12 b 10 以 8 3 2 2 00 0.1020.30.40.50.60.70.8 0.91.0 0.090.180.270.360.450.540.630.720.810.90 标准化时间 标准化时间 图5不同运动行时间关节轨迹曲线.(a)人工设定'aim=1s:(b)算法优化'aim=0.9s Fig.5 Joint trajectory curves at different motion times:(a)set value=1s:(b)optimized value=0.9s 4.2现场作业试验 螺栓试验,分别记录双臂的运行状态.通过各关节的运 为进一步验证改进算法的有效性及优越性,在现场 动规划,双作业臂由初始位姿到双臂实现引流板螺栓螺 线路上进行带电检修机器人的作业试验,用常规算法和 母对准对接位姿如图6所示,其中常规算法与改进算法 改进算法对机器人双作业臂进行轨迹优化,进行两次拧 作业臂1和作业臂2的轨迹运行状态如表1所示. 可 图6带电检修机器人现场运行.()初始状态:(b)对准状态 Fig.6 Field operation test of the robot:(a)initial state:(b)alignment state 表1作业臂运行状态对比表 通过现场作业试验机器人各关节协调运动,双作 Table 1 Comparison of robot operation states 业臂由初始状态进入到工作状态,再到螺栓螺母对准 作业手1 作业手2 对接状态,整个过程机器人双作业臂及其末端各关节 性能参数 常规改进常规改进 动作流畅、连续且平稳.由表1不同方法及不同被控 算法算法算法算法 作业臂在仿真与现场试验环境下的算法收敛速度可 仿真关节运动收敛时间/s 1.20.810.8 知:对于作业臂1仿真实验关节运动收敛时间,改进算 现场试验定位到作业对象时间/s 1028612498 法提升0.4s,实际运动作业臂1末端由初始位姿运动 平均对准逼近精度/% 90958898 到螺栓头对准位姿所耗费的时间缩短16s:对于作业
工程科学学报,第 38 卷,第 11 期 图 4 不同端点时刻关节轨迹曲线. ( a) t0 = 0,tf = 5 s; ( b) t0 = 1. 5 s,tf = 6. 5 s Fig. 4 Joint trajectory curves at different endpoint times: ( a) t0 = 0,tf = 5 s; ( b) t0 = 1. 5 s,tf = 6. 5 s 图 5 不同运动行时间关节轨迹曲线. ( a) 人工设定 tmotion = 1 s; ( b) 算法优化 tmotion = 0. 9 s Fig. 5 Joint trajectory curves at different motion times: ( a) set value tmotion = 1 s; ( b) optimized value tmotion = 0. 9 s 4. 2 现场作业试验 为进一步验证改进算法的有效性及优越性,在现场 线路上进行带电检修机器人的作业试验,用常规算法和 改进算法对机器人双作业臂进行轨迹优化,进行两次拧 螺栓试验,分别记录双臂的运行状态. 通过各关节的运 动规划,双作业臂由初始位姿到双臂实现引流板螺栓螺 母对准对接位姿如图 6 所示,其中常规算法与改进算法 作业臂 1 和作业臂 2 的轨迹运行状态如表 1 所示. 图 6 带电检修机器人现场运行. ( a) 初始状态; ( b) 对准状态 Fig. 6 Field operation test of the robot: ( a) initial state; ( b) alignment state 表 1 作业臂运行状态对比表 Table 1 Comparison of robot operation states 性能参数 作业手 1 作业手 2 常规 算法 改进 算法 常规 算法 改进 算法 仿真关节运动收敛时间/s 1. 2 0. 8 1 0. 8 现场试验定位到作业对象时间/s 102 86 124 98 平均对准逼近精度/% 90 95 88 98 通过现场作业试验机器人各关节协调运动,双作 业臂由初始状态进入到工作状态,再到螺栓螺母对准 对接状态,整个过程机器人双作业臂及其末端各关节 动作流畅、连续且平稳. 由表 1 不同方法及不同被控 作业臂在仿真与现场试验环境下的算法收敛速度可 知: 对于作业臂 1 仿真实验关节运动收敛时间,改进算 法提升 0. 4 s,实际运动作业臂 1 末端由初始位姿运动 到螺栓头对准位姿所耗费的时间缩短 16 s; 对于作业 ·1650·
江维等:高压输电线路带电检修机器人作业臂运动优化 ·1651· 臂2仿真实验关节运动收敛时间,改进算法提升0.2 Int J Sci Res Publ,2014,4(1)533 s,实际运动作业臂2末端由初始位姿运动到螺母对准 [5]Singh J,Gandhi K,Kapoor M,et al.New approaches for live 位姿所耗费的时间缩短26s.因此不论是何种工况下, wire maintenance of transmission lines.MIT Int J Electr Instrum 改进算法的收敛速度均快于常规算法.此外,为进一 Eng,2013,3(2):67 Song Y F,Wang H C,Jiang Y,et al.AApedD:a novel power 步量化末端与螺栓头及螺母的对准效果,在不同算法 transmission line maintenance robot for broken strand repair / 不同作业对象中引入平均对准逼近精度的性能指标 The 2nd International Conference on Applied Robotics for the Power 由表1可知:改进算法的平均对准逼近精度为98%和 Industry (CARPI 2012).Zurich,2012:117 95%,均高于常规算法,由于作业末端套筒具有柔性关 Zhan WW,Wang W,Chen N C,et al.Path planning strategies 节和一定自适应容错能力,允许误差在5%以内均可 for UAV based on improved A'algorithm.Geomat Inf Sci Wuhan Umi,2015,40(3):315 以实现末端与作业对象的对接,因此改进算法完全能 (占伟伟,王伟,陈能成,等.一种利用改进A*算法的无人 够满足对准对接要求;而常规算法由于平均对准逼近 机航迹规划.武汉大学学报(信息科学版),2015,40(3): 精度相对较低,有可能出现末端与作业对象无法对接 315) 的现象,这样需要再次人工干预调节双臂位姿,以提高 [8]Jin M J,Lou Y J,Liu G F,et al.A noval motion planning algo- 对准逼近精度,实现成功对接,由于常规算法需要人工 rithm with adaptive dynamic collision detection for industrial ro- 干预对准对接过程,在一定程度上降低作业效率.由 bots.J Unig Sci Technol China,2012.42(6):448 此可见,改进多项式插值关节轨迹规划方法具有较强 (金明杰,楼云江,刘冠峰,等.基于自适应动态碰撞检测的 工业机器人运动规划算法研究.中国科学技术大学学报, 的工程实用性,进一步提高作业效率,提高末端与作业 2012,42(6):448) 对象对准对接的可靠性,并且一定程度上体现了机器 9] Wang K M,Zhao Y F,Hou S P,et al.Dynamic obstacle avoid- 人作业的智能性 ance for unmanned underwater vehicle based on an improved artifi- 5结论 cial potential field.CAA/Trans Intell Syst,2014,9(1):47 (王奎民,赵玉飞,候恕萍,等.一种改进人工势场的UUV动 (1)研制出面向高压输电线路引流板螺栓紧固的 碍航物规避方法.智能系统学报,2014,9(1):47) 00] 带电检修机器人试验样机. Liu J F,Chen K.Su QQ,et al.Markov game based 3D path planning for palletizing robot.China Mech Eng,2012,23(7): (2)通过对轨迹时间变量的Min一Max标准化处 851 理,改进多项式插值算法轨迹仅仅只与轨迹端点状态 (刘久富,陈魁,苏青琴,等.基于Markov对策的码垛机器 和执行时间有关,而与端点时刻无关,提高了算法的实 人三维路径规划.中国机械工程,2012,23(7):851) 用性。 01] Feng A S,Wang Q H,Ge X Y.Trajectory planning of the live- (3)通过轨迹运动时间的优化,有效避免关节轨 line maintenance robot manipulator //2nd International Confer- 迹超调的发生,提高了关节运动的效率,同时满足轨迹 ence on Electronic and Mechanical Engineering and Information Technology (EMEIT-2012).Shengyang,2012:2089 状态全局约束要求,避免作业事故的发生 [12]Guan Y S,Yokoi K,Stasse 0,et al.On robotie trajectory plan- (4)现场作业试验进一步验证了改进算法的工程 ning using polynomial interpolations/IEEE International Con- 实用性. ference on Robotics and Biomimetics (ROBIO 2005).Shenzhen, 2005:111 参考文献 03] Wang X,Wang J W,Rao Z.An adaptive parametric interpolator [Hu Y.Research and development of live working technology on for trajectory planning.Adv Eng Sofucare,2010,41(2):180 transmission and distribution lines.High Voltage Eng,2006,32 [14]Zhang C.Li Q,Dong J Y.Switch strategy based on chaos parti- (11):1 cle swarm optimization and specialized genetic algorithm for path (胡毅.输配电线路带电作业技术的研究与发展.高电压技 planning of mobile robots.J Unie Sci Technol Beijing,2013,35 术,2006,32(11):1) (6):826 2]Lu S Y,Ma PS,Qi H,et al.Research on high voltage electric (张超,李擎,董冀媛,等.基于混沌粒子群专用遗传算法 power live line working robot.Autom Electr Power Syst,2003,27 切换策略的移动机器人路径规划.北京科技大学学报, (17):56 2013,35(6):826) (鲁守银,马培荪,威晖,等.高压带电作业机器人的研制 [15]Fahham H R,Farid M,Khooran M.Time optimal trajectory 电力系统自动化,2003,27(17):56) tracking of redundant planar cable-suspended robots considering 3]Lu S Y,Li Y P,Qi W.Robotie live-working for electric power both tension and velocity constraints.J Dyn Syst Meas Control, lines maintenances /ICIEA 2009,4th IEEE Conference on Indus- 2011,133(1):011004 trial Electronics and Applications.Xi'an,2009:1716 [16]Haddad M,Khalil W,Lehtihet H E.Traiectory planning of uni- 4]Chithamacharyulu P V,Sashanka D P,Kiran G U,et al.Design cycle mobile robots with a trapezoidal-velocity constraint.IEEE of automated hotline maintenance robot using haptic technology Trans Rob,2010,26(5):954
江 维等: 高压输电线路带电检修机器人作业臂运动优化 臂 2 仿真实验关节运动收敛时间,改进算法提升 0. 2 s,实际运动作业臂 2 末端由初始位姿运动到螺母对准 位姿所耗费的时间缩短 26 s. 因此不论是何种工况下, 改进算法的收敛速度均快于常规算法. 此外,为进一 步量化末端与螺栓头及螺母的对准效果,在不同算法 不同作业对象中引入平均对准逼近精度的性能指标. 由表 1 可知: 改进算法的平均对准逼近精度为 98% 和 95% ,均高于常规算法,由于作业末端套筒具有柔性关 节和一定自适应容错能力,允许误差在 5% 以内均可 以实现末端与作业对象的对接,因此改进算法完全能 够满足对准对接要求; 而常规算法由于平均对准逼近 精度相对较低,有可能出现末端与作业对象无法对接 的现象,这样需要再次人工干预调节双臂位姿,以提高 对准逼近精度,实现成功对接,由于常规算法需要人工 干预对准对接过程,在一定程度上降低作业效率. 由 此可见,改进多项式插值关节轨迹规划方法具有较强 的工程实用性,进一步提高作业效率,提高末端与作业 对象对准对接的可靠性,并且一定程度上体现了机器 人作业的智能性. 5 结论 ( 1) 研制出面向高压输电线路引流板螺栓紧固的 带电检修机器人试验样机. ( 2) 通过对轨迹时间变量的 Min--Max 标准化处 理,改进多项式插值算法轨迹仅仅只与轨迹端点状态 和执行时间有关,而与端点时刻无关,提高了算法的实 用性. ( 3) 通过轨迹运动时间的优化,有效避免关节轨 迹超调的发生,提高了关节运动的效率,同时满足轨迹 状态全局约束要求,避免作业事故的发生. ( 4) 现场作业试验进一步验证了改进算法的工程 实用性. 参 考 文 献 [1] Hu Y. Research and development of live working technology on transmission and distribution lines. High Voltage Eng,2006,32 ( 11) : 1 ( 胡毅. 输配电线路带电作业技术的研究与发展. 高电压技 术,2006,32( 11) : 1) [2] Lu S Y,Ma P S,Qi H,et al. Research on high voltage electric power live line working robot. Autom Electr Power Syst,2003,27 ( 17) : 56 ( 鲁守银,马培荪,戚晖,等. 高压带电作业机器人的研制. 电力系统自动化,2003,27( 17) : 56) [3] Lu S Y,Li Y P,Qi W. Robotic live-working for electric power lines maintenances / / ICIEA 2009,4th IEEE Conference on Industrial Electronics and Applications. Xi'an,2009: 1716 [4] Chithamacharyulu P V,Sashanka D P,Kiran G U,et al. Design of automated hotline maintenance robot using haptic technology. Int J Sci Res Publ,2014,4( 1) : 533 [5] Singh J,Gandhi K,Kapoor M,et al. New approaches for live wire maintenance of transmission lines. MIT Int J Electr Instrum Eng,2013,3( 2) : 67 [6] Song Y F,Wang H G,Jiang Y,et al. AApe-D: a novel power transmission line maintenance robot for broken strand repair / / The 2nd International Conference on Applied Robotics for the Power Industry ( CARPI 2012) . Zurich,2012: 117 [7] Zhan W W,Wang W,Chen N C,et al. Path planning strategies for UAV based on improved A* algorithm. Geomat Inf Sci Wuhan Univ,2015,40( 3) : 315 ( 占伟伟,王伟,陈能成,等. 一种利用改进 A* 算法的无人 机航迹规划. 武汉大学学报( 信息科学版) ,2015,40 ( 3) : 315) [8] Jin M J,Lou Y J,Liu G F,et al. A noval motion planning algorithm with adaptive dynamic collision detection for industrial robots. J Univ Sci Technol China,2012,42( 6) : 448 ( 金明杰,楼云江,刘冠峰,等. 基于自适应动态碰撞检测的 工业机器人运动规划算法研究. 中国科学技术大学学报, 2012,42( 6) : 448) [9] Wang K M,Zhao Y F,Hou S P,et al. Dynamic obstacle avoidance for unmanned underwater vehicle based on an improved artificial potential field. CAAI Trans Intell Syst,2014,9( 1) : 47 ( 王奎民,赵玉飞,候恕萍,等. 一种改进人工势场的 UUV 动 碍航物规避方法. 智能系统学报,2014,9( 1) : 47) [10] Liu J F,Chen K,Su Q Q,et al. Markov game based 3D path planning for palletizing robot. China Mech Eng,2012,23( 7) : 851 ( 刘久富,陈魁,苏青琴,等. 基于 Markov 对策的码垛机器 人三维路径规划. 中国机械工程,2012,23( 7) : 851) [11] Feng A S,Wang Q H,Ge X Y. Trajectory planning of the liveline maintenance robot manipulator / / 2nd International Conference on Electronic and Mechanical Engineering and Information Technology ( EMEIT-2012) . Shengyang,2012: 2089 [12] Guan Y S,Yokoi K,Stasse O,et al. On robotic trajectory planning using polynomial interpolations / / IEEE International Conference on Robotics and Biomimetics ( ROBIO 2005) . Shenzhen, 2005: 111 [13] Wang X,Wang J W,Rao Z. An adaptive parametric interpolator for trajectory planning. Adv Eng Software,2010,41( 2) : 180 [14] Zhang C,Li Q,Dong J Y. Switch strategy based on chaos particle swarm optimization and specialized genetic algorithm for path planning of mobile robots. J Univ Sci Technol Beijing,2013,35 ( 6) : 826 ( 张超,李擎,董冀媛,等. 基于混沌粒子群!专用遗传算法 切换策略 的 移 动 机 器 人 路 径 规 划. 北京科技大学学报, 2013,35( 6) : 826) [15] Fahham H R,Farid M,Khooran M. Time optimal trajectory tracking of redundant planar cable-suspended robots considering both tension and velocity constraints. J Dyn Syst Meas Control, 2011,133( 1) : 011004 [16] Haddad M,Khalil W,Lehtihet H E. Traiectory planning of unicycle mobile robots with a trapezoidal-velocity constraint. IEEE Trans Rob,2010,26( 5) : 954. ·1651·