考虑负载时变的线路巡检机器人动态性能分析

李小彭, 樊 星, 李 凯, 张凌越

(东北大学 机械工程与自动化学院, 辽宁 沈阳 110819)

传统对输电线路定期巡检的方法大多为人工检测,包括:爬行检测、地面监测和使用望远镜检测等.但随着电力系统的发展,这些检测方法表现出较大的限制性[1].相比人工检测,使用机器人检测的成本更低、速度更快、安全性更高,已成为国内外许多学者和机构的研究热点[2].目前对巡检机器人的研究多数集中于构型创新设计[3-5]与电力检测识别方面(如:电力系统组件[6]、线路[7]、电力塔[8]).考虑到机器人在工作中需要通过变换关节(如:旋转关节(R)、移动关节(P))姿态改变位姿来跨越诸如防振锤、悬垂线夹、绝缘子等线路障碍,而这个过程作为一个动态过程,受参数时变影响较大,所以需要对越障过程的动态性能进行分析.

对机器人的动态性能分析的传统方法通常为提出适当的指标评价惯量用以衡量系统的动态性能,如用于评价机械臂动力学性能的GIE方法[9]和DME方法[10].虽然惯量对机器人来说是影响动态性能的重要因素[11],但机器人作为一个机电系统,工作时的动态性能除受到动力学特性如惯性特性的影响,更与控制特性(如:控制结构与控制参数)密切相关[12],所以上述方法无法很好地使用在巡检机器人上.伺服电机作为机器人的驱动器,与控制特性密切相关.对伺服电机的很多研究,虽然传统控制模型中考虑了负载惯性,但只是将其作为扰动,很大程度上不会影响机电系统的稳定性和性能[13-14],所以使用这一类传统的模型的动态性能分析准确性不够.相比传统的控制模型,一些研究使用负载惯性控制为半闭环控制的双惯量模型更符合实际情况[15-16].以巡检机器人为例,负载惯性在关节空间与机器人的姿态有关,具有时变性,但基于双惯量模型同时考虑负载参数时变性方面研究却很少.

本文利用双惯量控制模型,在考虑负载惯量时变特性的情况下,研究了一种完整的动态性能分析方法.首先对机器人的关节臂建立了基于双惯量系统的控制模型与动力学方程.使用改进D-H参数法与拉格朗日法进行了机器人的运动学、动力学建模,从而使机器人越障时负载的时变性得到考虑.分析时主要通过研究负载惯量的变化对定常系统的影响来对动态性能进行分析,并给出了分析的步骤流程.最后,进行了包括速度端阶跃响应、位移端轨迹跟踪、系统关键参数变化等一系列仿真.

1 机构描述及控制系统模型建立

1.1 机构描述

图1为双臂输电线巡检机器人的三维模型.整机由两个关节臂与连接部件组成,每个关节臂包含了行走、旋转、升降、回转、平移机构;连接部件则由机架与重心调节机构组成.各机构由伺服电机驱动,通过传感器测量各机构实际的运行,控制器再通过误差对电机发出指令.整机共包含7个关节,可实现7个自由度的活动.

在巡检机器人执行输电线巡检任务时,需要跨越诸如防震锤、悬垂线夹、空中警示球等障碍物.越障作为一个动态过程,存在很强的时变性.

图1 越障中的双臂输电线巡检机器人

1.2 基于双惯量系统的关节控制模型的建立

关节空间控制为机器人主要的控制方法之一[13],为分析越障时负载端时变性对机器人造成的影响,首先需要对回转关节建立基于双惯量系统的控制模型,如图2所示.

图2 回转机构的双惯量模型

基于双惯量模型,得到如式(1)所示的动力学方程.

(1)

进一步对动力学方程拉氏变换,可推导出从TM到ωM与TM到ωL的传递函数,如式(2)、式(3)所示.

(2)

(3)

2 动力学建模

2.1 改进D-H参数法建模

双臂式巡检机器人多数为开环链式结构.建立坐标系时,将越障时抓住导线的行走机构视为基座,此时另一端脱离线路越过障碍的行走机构视为末端执行器.使用改进D-H参数法由基座起始依次建立共7个连杆坐标系,如图3所示.

图3 双臂巡检机器人D-H连杆坐标系

改进D-H参数法是通过4个参数即ai-1,αi-1,θi,di来描述从第i-1个连杆到第i个连杆的齐次变换矩阵,其中矩阵的形式如式(4)所示.

(4)

齐次变换矩阵连续乘积可将坐标系从基座变换到末端执行器,得到的总变换表达式如式(5)所示.

(5)

其中:p=[pxpypz]代表末端执行器的位置向量;R=[noa]表示末端方向的3个向量.

2.2 动力学方程

通过拉格朗日法对机器人的动力学进行分析.首先定义拉格朗日函数如式(6)所示.

L=Ek-Ep.

(6)

其中,Ek和Ep分别表示机构的总动能和总势能,具体形式如式(7)、式(8)所示.

(7)

(8)

由广义坐标下的拉格朗日方程可建立机器人的动力学方程,如式(9)所示.

(9)

根据式(9),得到巡检机器人的动力学方程如式(10)所示.

(10)

式中:

3 考虑负载时变性的动态性能分析

机器人回转机构的双惯量模型下的控制结构框图如图4所示,具体包含位置环、速度环和电流环.图中速度环控制器和电流环控制器为PI控制器,位置环控制器为P控制器.由于速度外环与电流内环的带宽差距,可认为在速度外环调节时,电流内环的调节已经完成,所以在下文中将Kt视为1进行分析[16].

根据图4得到从ωL到ω*的传递函数,如式(11)所示.

(11)

图4 双惯量模型的控制结构框图

其中:Ksp,Tsi分别代表速度环控制器的增益与时间积分常数;ω*为期望角速度.

同样位置控制下从θL到θ*的传递函数,如式(12)所示.

(12)

式中:e=Kpp/(s+Kpp);Kpp表示位置环控制器增益;θ*为期望位移.

由式(11)、式(12)可知转速ωL、转角θL与反谐振频率ωa、谐振频率ωn相关.所以具有时变性的负载参数会影响控制效果而影响机器人稳定.

对动态性能变化分析前需要计算控制器参数.将闭环传递函数式(11)化简为形式为式(13)所示传递函数与形式为分母包含2个二阶环节组合的如式(14)所示传递函数.

(13)

式中:a=JMTsi;b=KspKtTsi;c=JMTsiωn+Kt;

(14)

其中:ω1和ω2代表自然角频率;ξ1和ξ2为自然角频率对应的阻尼系数.

令传递函数(13)、(14)的特征方程逐项相等得到式(15)～式(18).其中,式(15)、式(16)为控制器参数计算表达式;式(17)、式(18)为极点的限制方程.因为方程存在4个极点,而控制器参数只有2个,所以极点不能自由配置,需要其他限制条件.

(15)

(16)

(17)

(18)

考虑二阶系统受阻尼系数ξ影响较大,使用极点配置法中使两对极点阻尼系数相同的方法作为限制条件来进一步确定控制器参数.

将ξ1=ξ2=ξ代入限制方程(17)、方程(18)可解出ω1,ω2表达式如式(19)所示.

(19)

式中:A=4ξ2Kt;B=(Kt-1)2.因为自然角频率ω1,ω2均不小于0,由此可得阻尼系数ξ表达式如式(20)所示.

(20)

位置环开环传递函数如式(21)所示,其中参数Kpp的值可通过根轨迹法求得.

Gp(s)=

(21)

假设此时在某个确定位姿,通过上述步骤计算可得控制器参数,进一步计算传递函数式(13)特征方程式的各项系数,求得特征方程的根,结果如式(22)所示.

(22)

因为特征方程含义等效,当用式(14)计算某控制器参数下特征方程的根,得到的含阻尼系数和自然角频率参数的解的形式如式(23)所示.

(23)

(24)

4 数值仿真

巡检机器人在实际执行巡检任务时需要跨越诸如悬垂线夹、警示球、防震锤等障碍.其跨越悬垂线夹具体越障过程如图5所示,即机器人要先移动重心调节机构,抬升关节臂使行走机构脱离输电线,再通过两侧关节臂的回转机构同时旋转翻越障碍.在单臂对悬垂线夹的跨越完成后,另一只臂循环上述步骤,完成整机越障.

越障作为一个动态过程受负载时变性影响较大.由图5可知在越障第②步姿态变化时,线路上仅有一个行走机构,此时受负载时变影响较大,容易造成工作时动态性能的不稳定.

图5 机器人越障工况示意图

4.1 负载惯量的参数计算

将动力学方程中惯量矩阵M(q)最大奇异值作为不同姿态下伺服电机的负载转动惯量值[13],计算得到以关节3与关节5角度值为自变量参数下的负载惯量值与惯量比的变化结果,如图6所示.

由图可知,在越障时发生的机器人姿态变化会导致负载转动惯量的变化,具体变化范围为从小于0.01 kg·m2到大于0.03 kg·m2;对应惯量比的变化范围为从小于1.5到大于3.5.即越障时负载惯量会因机器人姿态的变化而表现出强烈的参数时变性.

图6 惯量及惯量比计算结果

4.2 负载端速度的仿真分析

仿真前先计算控制器参数,由图6按照最大、中等、最小的次序分别以负载转动惯量JL1=0.032 2 kg·m2(R=3.788),JL2=0.016 7 kg·m2(R=1.965)及JL3=0.009 8 kg·m2(R=1.153)为参考,对不同惯量下负载端速度变化进行仿真分析.

仿真所需各项参数如表1所示.以最大负载转动惯量和最小负载转动惯量为参考数值计算控制器参数分别如表中参数1、参数2所示.

表1 仿真参数

仿真负载端速度阶跃响应如图7所示.图7为在参数1、参数2下的控制器参数进行仿真的阶跃响应图.

由图可知,无论采用哪一组控制参数,越障时负载参数的变化会导致机器人动态性能都会发生变化,尤其是振动幅值的变化;对比图7a与图7b,可得较大惯量参数计算得到的仿真参数下的阶跃响应相比另一种,总体振动幅值更小.

图7 速度阶跃响应

4.3 负载端位移的仿真分析

根据位置环开环传递函数绘制系统根轨迹图,如图8所示.图中箭头方向为增益参数变大方向,依据根轨迹图选择能满足系统稳定性的比例增益参数,结果为Kpp=40.由图可知,虽然使用的控制结构相同,但受负载惯性参数变化的影响,根轨迹是不同的,动态性能也是不同的.

输入正弦函数作为参考轨迹,位置环的轨迹跟踪误差曲线如图9所示.受负载惯量时变性的影响,两图在最大转动惯量即R=3.788的初始阶段存在较为明显的振动;此外,由较大的负载惯量确定的控制参数下系统的跟踪性能更好.

图8 根轨迹图

4.4 考虑负载时变下的关键参数变化分析

机器人在越障时发生姿态变化会导致动态性能的变化.观察系统特征方程可知,其根本原因在于系统的阻尼比与固有角频率的变化.

由式(22)～式(24)计算机器人越障时整个工作空间中在表1中参数1、参数2下负载时变引起的系统阻尼系数ξ和固有角频率ω的变化规律,结果分别如图10、图11所示.

由图10可知,参数1下阻尼系数ξ和自然角频率ω会随关节角的变化而变化.其中阻尼系数变化范围分别为从0.8到大于2和从小于0.2到0.8;自然角频率的变化范围分别为从小于60到大于70和小于120到大于180.图11中参数2阻尼系数ξ和自然角频率ω同样随关节角的变化而变化.其中阻尼系数变化范围分别为从小于0.2到0.4和0.4到大于0.5;自然角频率的变化范围分别为从小于65到大于125和小于225到大于240.

所以固定控制参数不能保证系统的动态性能,由姿态变化引起的阻尼系数和自然角频率的变化会导致特征多项式根的变化,从而进一步影响机器人的动态性能.

图9 位移轨迹误差曲线

图10 参数1下关键参数变化图

图11 参数2下关键参数变化图

5 结 论

1) 研究了一种用于双臂式输电线巡检机器人在越障时,分析负载惯量参数变化如何对机器人动态性能影响的方法.具体包括基于双惯量控制系统的建模、机器人整机动力学建模、控制器参数设计计算以及系统参数变化的分析.

2) 负载惯量的时变性对机器人动态性能有较大影响,具体表现为越障时位姿变化导致的负载惯量变化会加大机身振动程度,所以固定控制参数不能很好保证系统的动态性能;此外选择控制器参数时,计算范围内较大的阻尼系数下的控制器参数总体控制效果更好.

3) 研究成果可用于评价越障时机器人的动力学特性,并可用于机器人电机伺服参数调试过程以提高机器人的动态性能.