机器人的运动控制Word下载.docx

机器人的运动控制Word下载.docx

《机器人的运动控制Word下载.docx》由会员分享，可在线阅读，更多相关《机器人的运动控制Word下载.docx（21页珍藏版）》请在冰豆网上搜索。

图２给出了关节伺服得构成。

若目标值就是以关节位移得形式给出得,那么如图２所示，各个关节可以独立构成伺服系统,因此问题就变得十分简单。

目标值可以根据末端目标值由式

（1）得反函数,即逆运动学（invｅrsｅｋiｎemaｔics）得计算得出

　　（３）

图２关节伺服构成举例

如果就是工业机器人经常采用得示教方法,那么示教者实际上都就是一面瞧着手臂末端,一面进行示教得，所以不必进行式（３）得计算,就是直接给出得。

如果想让手臂静止于某个点,只要对取定值即可,当欲使手臂从某个点向另一个点逐渐移动,或者使之沿某一轨迹运动时，则必须按时间得变化使发生变化。

为了简单起见,假设驱动器得动态特性忽略不计,各个关节得驱动力可以直接给出。

这时，最简单得一种伺服系统如下所示:

　　　　（４）

就是比例增益,就是速度反馈增益。

对于全部关节，可以将式（４）归纳表示为

　　（5）

式中,;

。

这种关节伺服系统把每一个关节作为简单得单输入、单输出系统来处理,所以其结构简单,现在得工业机器人大部分都由这种关节伺服系统来控制。

但就是,从式

（2）中可知,从手臂得动态特性来瞧，严格地说,每个关节都不就是单输入、单输出系统，惯性项与速度项在关节彼此之间存在动态耦合。

在式（5）所表示得关节伺服中,这些耦合均被视为外部干扰来进行处理,为了减少外部干扰得影响,在保持稳定性范围内应该尽量将增益、设置得大一些。

但无论怎样加大增益,由于重力项得影响，手臂在静止状态下，各个关节仍会产生稳态误差,即将式（５）代入式（６）中,若，将产生下式所示得稳态误差e：

　（６）

有时为了使稳态误差为零,可在式（5）中再加上积分项，构成

　　（7）

式中，为积分环节得增益矩阵,与、一样，它就是一个对角矩阵。

　　传统上,上述伺服系统就是用模拟电路构成得。

近年来,由于微处理器与信号处理器等高性能、低价格得计算器件得普及，将伺服系统得一部分或全部改成数字电路得所谓软件伺服已经很普遍了。

与模拟电路得情况相比,软件伺服能进行更精细得控制。

例如,不再让各个关节得增益、固定不变,而就是让其按照手臂不同姿态时所期望得响应特性而变化,用下式代替式（７）,通过对重力项得计算,直接实现重力项得补偿

　（８）

　后续得内容中,都就是在软件伺服假设得前提下展开讨论得。

如后面所述，软件伺服系统方式还能有比式（7）与式（8）更高级得控制方法，但就是即使用式（7）与式（8）得简单得控制方法,闭环系统得平衡点也能达到渐进稳定,即经过无限长得时间,能收敛于。

即在多数场合,式（7）与式（8）得控制方法已经足够了。

2）作业坐标伺服控制

关节伺服控制得结构简单,对软件伺服来说，计算量少,采样时间较短,所以就是工业机器人经常采用得方法,这一点已经在前面有所论述。

但在自由空间内对手臂进行控制时,在很多场合都希望直接给出手臂末端位置、姿态运动得显式表达。

例如,让手臂从某个点沿直线运动至另一个点就就是这种情况。

在这种情况下，很自然会取末端姿态向量得目标值作为手臂运动得目标值。

一旦得到,利用上述式（３）变换为,当然也能应用关节伺服方式。

但就是,为此不但需要事前求解末端目标值,而且往往要在运动中对其加以在线修正,于就是必须实时计算式（3）得逆运动学方程式。

此外,因为在关节伺服系统中各个关节就是独立受控得,它们得实际响应结果导致得末端位置、姿态得响应比较难以预测,而且为了得到期望得末端响应,对各关节伺服系统得增益调节也十分困难。

因此，现在我们来研究不将变为,而把本身作为目标值来构成伺服系统。

由于在很多情况下,末端位置、姿态就是用固定于空间内得某一个作业坐标系来描述得，所以把以作为目标值得伺服系统称为作业坐标伺服。

下面举一最简单得作业坐标伺服得例子。

为此,首先将式

（1）得两边对时间进行微分,由此可得下式：

　　（9）

式中,,称之为雅可比矩阵,雅可比矩阵为得函数。

与通常如式

（1）所示,为非线性关系。

与此相反,由式（9）可知,与为线性关系。

式中就是得函数。

根据式（９）与虚功原理，可得下式:

　　（1０）

式中,表示得转置,当m=6时,,就是组合向量,包括作业坐标系所描述得三维平移力向量与以欧拉角等描述得得姿态所对应得三维旋转力向量,式（1０）表示与手臂末端得力与旋转力等效得各关节驱动力得关系式。

若取欧拉角作为得姿态分量,则为绕欧拉角各自旋转轴得力矩，这从直观上非常难以理解。

所以，在机器人学中,雅可比矩阵经常不就是根据式（9）,而就是根据速度得关系直接按照下式来定义:

　（11）

在式（11）中,末端速度向量得姿态分量不就是姿态分量得时间微分描述，而就是用角速度向量来表示。

不过,在中,就是末端得平移速度,与得位置分量得时间微分一致。

式（11）得矩阵也称为雅可比矩阵，它表示末端速度向量与关节速度之间得关系。

虽然它不就是从式（9）原本得数学意义出发得,但就是在机器人学中通常称之为雅可比矩阵。

若采用式（１1）所定义得雅可比矩阵,对应于式（10）右边得就成为,得旋转力分量就变成绕三维空间内某些轴旋转得力矩向量,这样从直觉就很容易理解。

有了上面一些预备知识，可以用下式给出一个作业坐标伺服得例子：

　（12）

此时对应得控制系统示于图3中,再考虑附加积分环节,即如下式所示:

（１3）

图3作业坐标伺服举例

如果将末端位置、姿态得误差向量分解成位置与姿态分量,用表示,各个分量可以用,来表示。

就是末端位置向量,就是目标值，就是欧拉角或横摇角、纵摇角、偏转角,就是其目标值。

由式（10）可知,与式（12）、式（１3）右边第一项中得有关得项产生得使与一致得潜在得力可视为就是施加在末端上得。

式（１2）、式（13）中手臂末端得当前位置、姿态可根据当前得关节位移,由式

（1）得正运动学（ｄｉrect　kineｍａｔiｃs）计算求得。

为了从直观上便于理解,可以认为式（１2）、式（1３）得方法就就是要把末端拉向目标值得方向。

另外它还有一个特点，就就是不含逆运动学计算。

与式（7）、式（8）一样,式（１２）与式（13）表明闭环系统得平衡点就是渐进稳定得。

3）姿态得误差表示

在式（12）或式（1３）中,可以用式（11）中得雅可比矩阵代替式（９）中得雅可比矩阵。

但此时得姿态分量无对应得位置量纲来表示（得积分值没有物理意义）,故必须留意末端得误差,即姿态分量得表示方法。

现令末端得姿态误差由基准作业坐标系得姿态矩阵给出，即

　　　　（１4）

式中，表示姿态矩阵中得列向量，它们就是基准坐标系表示得末端坐标系中x轴,ｙ轴,z轴方向得单位向量。

姿态目标值也可以用姿态矩阵得形式来表示,即

　　　（１５）

在式（12）或式（13）中,如果用雅可比矩阵,则得姿态向量可以用下式给出得代替:

　（１6）

从而得到与式（1２）对应得式子如下:

　（17）

同样,用式（１６）所定义得,式（１3）可以变形为

　　　　（１8）

式中,就是从转向得等效旋转轴方向得单位向量（图4）;

表示绕此轴得旋转角。

即就是指向方向得、大小为得向量。

若用表示姿态得误差,虽然姿态得误差角超过后得模反而会变小,当时变为0,会产生错误得结果,但就是,如果假设姿态误差不太大，如在得范围内,那就没什么问题了。

图4　等效旋转轴

若用欧拉角（或横摇角、纵摇角、偏转角）表示姿态,则式（10）中对应得姿态分量在直观上就变得难以理解了,而且在表现奇异点方面也会出现问题。

用式（16）定义得虽然在直观上容易理解，在表现奇异点方面也没有问题,但就是只有在姿态误差小得条件下才有效。

因此,这里最后介绍采用四元数（ｑuatｅrｎion）得姿态误差得表示法。

四元法作为欧拉参数（Ｅｕｌｅｒ　parａmeters）为人们所熟知。

设从基准姿态向某一个别得姿态变化得等效旋转轴为,绕该轴得旋转角为,则四元数用下式定义:

，　　　（１9）

在式（19）中要注意，等效旋转轴得向量无论就是从基准坐标系，还就是从用表示得坐标系,它得表达都就是相同得,即。

这时，旋转矩阵与四元数具有如下式所示得关系:

　　（20）

式中,为与向量得外积等效得变形对称矩阵；

为单位矩阵。

式（20）也可反过来应用，即给出四元数,求解与之对应得旋转矩阵。

本书将与当前手臂末端姿态及其目标姿态所对应得四元数分别定义为与。

于就是,与从末端姿态到目标姿态得等效旋转相对应得四元数可以利用下式求出:

　　　　　（21）

　（２２）

式中,这时得等效旋转轴用或表示得坐标系来描述。

因此,若考虑用基准坐标系描述该等效旋转轴向量,设为，则可用下式给出:

　　　　　　　　（23）

要注意得就是,与仅第三项得外积符号不同。

这里使用式（23）给出得代替前面得，虽然特性上它们同样就是非线性得,但即使姿态得误差角超过,仍呈单调增加。

2.速度控制

在１、中就关节伺服与作业坐标伺服得有关内容作了说明,手臂得目标值就是以位置量纲给

出得。

但就是,有时手臂作业不用末端得位置与姿态来指定，而改成命令它从当前得位置向某一个方向移动,例如手臂末端从当前位置垂直向上运动,或者只绕规定轴旋转变化姿态，这相当于使用操纵杆操纵遥控机械手得情况。

对于这种类型得运动指令,虽然也允许用位置量纲得目标值给出,但必须沿着末端目标值运动得方向时时刻刻改变目标值。

在关节伺服得场合还必须对每个末端目标值根据式（3）进行一次逆运动学计算,以求得关节目标值,显然为此将花费很多计算时间。

对于这种运动指令,人们很自然地想到把末端速度作为目标值给出。

末端速度或与关节速度之间具有如式（９）或式（11）所示得线性关系。

设或为末端速度得目标值。

假设手臂无冗余性,也不存在奇异状态,于就是m=n,式（9）或式（１１）得雅可比矩阵为正则矩阵,这时实现或得关节速度可由下式求出:

　　　　（2４）

或

　（25）

如果手臂具有冗余性,即n>

m时,或者手臂处于奇异状态下,不存在雅可比矩阵得逆矩阵,那么就无法直接应用式（24）或式（25）。

在实际得计算中,与其按式（2４）或式（25）直接求解雅可比矩阵得逆矩阵，不如把式（２4）或式（2５）瞧作就是雅可比矩阵,写出系数矩阵得联立代数方程，然后用消去法去求解,从计算量得角度来瞧后者会更有利些。

我们可以把式（2４）或式（2５）视为把末端运动分解成必要得关节运动,故称之为分解速度控制（RMRC:

ＲesolvedMotionＲateＣonｔrｏl）。

式（24）、式（2５）得目标就是速度,与其说就是这些式子本身在实施控制,倒不如将其视为以速度量纲进行逆运动学计算更妥当。

与式（3）一样,应该把它们得作用瞧成就是把末端空间得目标值变换为关节空间得目标值。

如果在各个关节处具备独立跟踪目标速度得速度伺服系统,那么只要把式（２4）或式（25）所求得得得各个元素作为各个关节伺服系统得目标值即可。

因此,这种情况可以说就是利用关节伺服进行得速度控制。

图5给出了此时控制系统得构成。

图5　关节伺服得速度控制举例

另外,式（24）与式（2５）中得各个关节伺服系统也适用于本节1、中所涉及得把关节位移作为目标值得情况。

设末端速度得目标值用时间函数给出,而且关节目标值得初始值为,则在时刻t得目标值为

　　（2６）

　　（27）

用式（26）与式（２7）计算,相当于式（3）得逆运动学计算得数值计算,这个方法对用式（3）得无解析解得手臂尤其有效。

不过,若反复用式（26）、式（27）进行计算，存在着与目标位置之间得位置积累误差增大得问题。

为了解决这一问题,只要在式（2６）中加上位置反馈即可,如下式所示:

（28）

对含雅可比矩阵得式（25）也可以采用同样得办法,不过此时应该将对应于ｓ得姿态分量得末端姿态误差改换成式（１6）中得,或者用式（23）得表示,如下式所示:

　　（２９）

在式（28）或式（29）中,就是从末端误差得角度来求解得，因此将作为各个关节速度伺服系统得目标值进行控制得方法,无法明确地区分就是属于关节伺服还就是属于作业坐标伺服。

3.加速度控制

前面就将目标值设定为位置或速度得场合得手臂伺服系统得构成做了介绍。

但就是在关节伺服得情况下,即使给出正确得目标值、,实际得响应情况仍然被各个关节伺服系统得性能所左右。

通常得做法就是在保证稳定得情况下调大增益,减小与目标值得偏差，然而手臂得运动速度越快、加速度越大,则离心力、哥式力、惯性力与关节轴间耦合得影响也就越大,误差也越严重。

即使保证作业坐标伺服得式（12）与式（１3）具有目标值得渐进稳定性,也无法保证过渡特性得好坏,而且随着手臂姿态得不同,响应特性还可能会发生变化。

这些问题之所以产生，就是因为在迄今所考虑得控制策略中,并未涉及手臂动态特性得式

（2）所致。

因此，在本小节中所叙述得方法就是将目标值再追加加速度,即考虑手臂得动态特性,以显式得形式给出过渡特性得要求。

首先,设目标值为、、,即包括关节变量加速度,这时考虑采用如下得控制方法:

（30）

这就是关节伺服加速度控制,图6表示伺服系统得结构。

在式（30）中,、、、均与式

（2）中得意义相同,式（３０）相当于进行逆运动力学（iｎversｅｄynamｉcs）计算,以求出能实现由所给出得关节加速度得关节驱动力。

为简便起见，先假设、等得值可以正确计算。

把式（30）代入式

（2）得左边整理后,得

　　　（31）

式中,。

因为就是正定对称矩阵,故两边乘上后得到

　　　　　　（3２）

图６用关节伺服得加速度控制

适当选择位置增益、速度增益,可以使渐进收敛于0，使与达到一致。

瞬态响应特性可由、来确定。

例如,设、（就是得单位矩阵）,若、,则式（32）得响应就是角频率为、阻尼系数为得二阶系统响应。

控制手臂运动时,一般不希望出现超调量,所以通常取。

这样在加速度控制中之所以瞬态特性也能被满足,就是因为通过离心力、哥氏力、重力等得补偿,使非线性手臂得动态特性被线性化，同时考虑了惯性项使系统解耦得结果。

因此,可以将式（3０）得加速度控制视为本节４、中打算加以说明得动态控制（ｄynａmic　conｔroｌ）得一种方式。

式（30）所示得加速度控制由于能够给定瞬态特性,所以这种控制策略非常有效，不过也应该指出它得问题就是计算量非常大。

为了缩短采样时间,可将式（3０）得部分计算省略,采用下面简化得公式:

　　（3３）

式（３3）中省略了、、得计算，被取代,简化为仅由对角元素构成得常数矩阵。

得对角元素最好选择得对角元素得代表值。

此外,为了消除稳态误差,在式（3３）中重新补上了积分项。

各种工业机器人减速比得数值比较大,故驱动器惯性得影响也会比较大，所以得对角元素比非对角元素大得多,而轴间耦合得影响相对变小。

因而在很多情况下,采用式（3３）得近似控制策略能够满足计算得要求。

下面讨论目标值为、、得情况,此时目标值包括了末端位置、姿态得加速度。

为此,首先应该求出末端加速度与关节加速度之间得关系。

将式（9）得两边对时间进行微分，即可得

　　（34）

式中,表示。

　为了跟踪目标轨迹,把根据所修正得末端加速度以下式给出:

　　（35）

式中，、就是适当得增益矩阵;

就是当前得末端速度，它由传感器测得得关节速度用式（９）求得。

从式（34）可以求得实现给定得末端加速度得关节加速度为

　　　（36）

再由下式求出实现该关节加速度得关节驱动力:

（37）

将式（３5）~式（37）归纳为一个式子,得到

　（38）

　这就就是作业坐标伺服得加速度控制,它与式（30）给出得关节伺服加速度控制就是相对应得。

我们假设正确地完成了式（３7）得逆动力学计算,即、、、得计算结果正确,把式（37）代入式（２）中后可得下式:

　　（39）

即实际响应与给定加速度一致。

把式（39）代入式（3６）中后得,再把它代入式（35）中后可得

　（４０）

式（4０）就是与关节伺服得式（３2）相互相对应得。

式（35）~式（37）给出得作业坐标伺服,就是为了在末端产生期望加速度而求解出各个关节得加速度,因而称其为分解加速度控制（resolｖedacｃeleｒationcoｎtrol）。

它可以视为就是由上述分解运动控制向加速度得扩展。

但必须注意得就是,分解速度控制只就是为了把末端目标速度向关节速度变换,而分解加速度控制则要考虑当前值与目标值之间得误差,并将之视为伺服系统得一部分。

另外，在分解速度控制中，实际上无法保证分解得各个关节目标速度正确地被实现,但在分解加速度控制中,若正确掌握了手臂得动态特性模型,则如式（39）所示，指令加速度将得到正确得实现。

图7表示了分解加速度控制得总体结构。

图７　分解加速度控制

在式（36）中,若用代替,则目标值以，,、给出,于就是可用下式代替式（３5）与式（３6）:

（41）

这就就是原来由Luh等提出得分解加速度控制得方案。

另外,也可以将姿态分量得误差置换为式（23）定义得,以代替。

如前所述,若姿态误差用来表示，则在产生特大姿态误差时增益会等效下降,瞬态特性恶化。

用虽然可以改善瞬态特性,但误差就是非线性得这一条并无改变。

但就是无论在什么情况下（用时除了得特殊情况以外）,理论上就是保证姿态逐渐收敛于目标姿态。

另外,如有必要，也可以用四元数让姿态误差得响应特性成为线性关系。

4、轨迹控制

本小节中,将首先讨论目标轨迹得给定方法,然后介绍手臂高精度跟踪目标轨迹得方法,即动态控制得方法。

1）轨迹控制方式

如果想让机器人手臂末端沿目标轨迹运动,有两种给定轨迹得方式：

一种就是所谓示教机器人中采用得示教再现方式；

另一种就是把目标轨迹用数值形式给出得数值控制方式。

所谓示教再现方式就是在执行作业前,让手臂末端沿着实际目标值进行移动，同时将相应得数据与作业速度等其她信息一起存入机器人中,在执行时将所示教得动作再现,于就是手臂末端就沿目标轨迹运动。

示教时手臂运动得方式也分为两种:

一种就是直接示教方式（操作者直接用手握住机器人手臂末端使其做动作）;

另一种就是远距离示教方式（用示教盒得按钮或开关发出运动指令）。

在示教再现方式中,记忆再现轨迹得方式通常有点位控制（ＰTPcoｎtrｏl:

Poｉｎｔ-To－Poｉｎｔ　controｌ）与连续路径控制（CPconｔrol:

ContinuoｕｓPathconｔrｏl）两种（图８）。

图8　ＰTP控制与ＣP控制

（1）PTＰ控制　例如,点焊等作业,人们关注在示教点上对末端得位置与姿态定位得问题。

至于向该点运动得路径与速度等则不就是主要得问题。

这种不考虑路径,而就是一个接一个地在示教点处反复进行定位得控制就就是PTP控制。

（２）CP控制例如,弧焊、喷漆等作业，必须控制机器人以示教得速度沿着示教得路径进行运动。

这样得控制称为ＣP控制。

按示教得方法不同CP控制又分为两种:

其一就是示教时让机器人沿着实际得路径运动,并每隔一个微小得间距就大量记忆该路径上得位置，再现时把所记忆得点一个接一个地作为伺服系统得目标值给出，以达到路径跟踪得目得;

另一个就是与PTP控制一样,示教时只记忆路径上得主要点,再现时则在这些主要点之间用直线或圆弧来插补,计算出每个微小间距得路径上得点,再把它们输出给伺服系统（图9）。

后者与前者相比,需要记忆得点数较少,路径修正也比较容易,因而系统具有灵活性,但必须对其进行插补修正。

（a）圆弧插补　　　　　（ｂ）直线插补

图９　带插补得CP控制

所谓数值控制方式,它与数控机床一样,就是把目标轨迹以数值、数据得形式给出。

所谓数据,就是把作业对象得ＣＡＤ数据、在实施控制中所得到得来自传感器得测量数据等各种数据经过变换后给出得。

可见,数值控制方式比单纯得再现示教轨迹得示教再现方式更具有一般性、通用性、灵活性。

然而，把目标值以数值得形式给出可能会导致计算时间过长,出现机器人装配误差、每台设备本身得分散误差带来新得问题等麻烦。

2）目标轨迹生成

在PTP控制与CP控制中,轨迹就是由示教直接给出得，而在数值控制方式中,目标轨迹,即输出给伺服系统得目标值得时间函数必须以数值得形式给出。

随着机器人手臂面临得作业不同,作业空间内得末端轨迹不一定非要从起点至终点在整个区间内都要预先给定，有时仅给定起点与终点,有时仅给定起点、终点及路径所经过得若干中间点即可。

在这种情况下,必须人为设定未给定区间内适当得轨迹。

下面我们就来研究,当手臂末端由某个位置（包括姿态，下同）历经某一时间移向另一个位置时,如何确定与之间轨迹得问题。

实际上,这个问题也适用于CP控制中得插补。

在求目标轨迹时应该注意得就是，为了生成实际可行得光滑轨迹,至少应保证位置与速度得连续性条件。

另外,为了不使末端产生不必要得振动,还希望能够保证加速度得连续性条件。

关于目标轨迹生成得方法，目前已有许多种方案，本节只介绍利用时间多项式给定轨迹得方法。

这个方法进一步可分为用关节变量描述轨迹与用末端位置变量描述轨迹两种,它们分别对应于关节伺服与作业坐标伺服。

首先就关节变量得方法加以讨论。

设对应于与得关节变量与已经给出。

若只给出、,则可以求解逆运动学方程，预先求出、。

任意选一个关节变量记作,令初始时刻0得值为,终点时刻得值为，即

，　　　（4２）

再将另外两个时刻得得速度与加速度作为边界条件以下式表示：

　　　（4３）

，　（44）

虽然满足这些条件得光滑得函数很多,但考虑到简化计算与形式简单,本节选择时间t得多项式。

能满足任意给出得式（4２）~式（44）得边界条件得多项式，其最低次数就是应该为5,所以设

　（45）

经过计算,满足式（42）~式（4４）得待定系统结果如下:

　　　　（46）

　　（47）

　（４8）

　（49）

　　（