第五章机器人动力学ppt课件

1、 前面我们所研讨的机器人运动学都是在稳态下进展的，没前面我们所研讨的机器人运动学都是在稳态下进展的，没有思索机器人运动的动态过程。实践上，机器人的动态性能不仅有思索机器人运动的动态过程。实践上，机器人的动态性能不仅与运动学相对位置有关，还与机器人的构造方式、质量分布、执与运动学相对位置有关，还与机器人的构造方式、质量分布、执行机构的位置、传动安装等因案有关。机器人动态性能由动力学行机构的位置、传动安装等因案有关。机器人动态性能由动力学方程描画，动力学是思索上述要素，研讨机器人运动与关节力方程描画，动力学是思索上述要素，研讨机器人运动与关节力(力矩力矩)间的动态关系。描画这种动态关系的微分方程称

2、为机器人间的动态关系。描画这种动态关系的微分方程称为机器人动力学方程。机器人动力学要处理两类问题：动力学方程。机器人动力学要处理两类问题： 动力学正问题和逆问题。动力学正问题和逆问题。 动力学正问题是动力学正问题是根据关节驱动力矩或力，计算机器人根据关节驱动力矩或力，计算机器人的运动的运动(关节位移、速度和加速度关节位移、速度和加速度)； 动力学逆问题是动力学逆问题是知轨迹对应的关节位移、速度和加速知轨迹对应的关节位移、速度和加速度，求出所需求的关节力矩或力。度，求出所需求的关节力矩或力。 不思索机电控制安装的惯性、摩擦、间隙、饱和等要素时不思索机电控制安装的惯性、摩擦、间隙、饱和等要素时，n

3、 自在度机器人动力方程为自在度机器人动力方程为n个二阶耦合非线性微分方程。个二阶耦合非线性微分方程。方程中包括惯性力方程中包括惯性力/力矩、哥氏力力矩、哥氏力/力矩、离心力力矩、离心力/力矩及重力力矩及重力/力力矩，是一个耦合的非线性多输入多输出系统。对机器人动力学矩，是一个耦合的非线性多输入多输出系统。对机器人动力学的研讨，所采用的方法很多，有拉格朗日的研讨，所采用的方法很多，有拉格朗日(Lagrange)方法、牛方法、牛顿一欧拉顿一欧拉(NewtonEuler)、高斯、高斯(Gauss)、凯恩、凯恩(Kane)、旋量、旋量对偶数、罗伯逊一魏登堡对偶数、罗伯逊一魏登堡(RobersonWit

4、tenburg)等方法。等方法。 研讨机器人动力学的目的是多方面的。研讨机器人动力学的目的是多方面的。 动力学正问题与机器人的仿真有关；动力学正问题与机器人的仿真有关； 逆问题是为了实时控制的需求，利用动力学模型，实现最逆问题是为了实时控制的需求，利用动力学模型，实现最优控制，以期到达良好的动态性能和最优目的。在设计中需优控制，以期到达良好的动态性能和最优目的。在设计中需根据连杆质量、运动学和动力学参数、传动机构特征和负载根据连杆质量、运动学和动力学参数、传动机构特征和负载大小进展动态仿真，从而决议机器人的构造参数和传动方案大小进展动态仿真，从而决议机器人的构造参数和传动方案，验算设计方案的合

5、理性和可行性，以及构造优化程度。，验算设计方案的合理性和可行性，以及构造优化程度。 在离线编程时，为了估计机器人高速运动引起的动载荷和在离线编程时，为了估计机器人高速运动引起的动载荷和途径偏向，要进展途径控制仿真和动态模型仿真。这些都需途径偏向，要进展途径控制仿真和动态模型仿真。这些都需求以机器人动力学模型为根底。求以机器人动力学模型为根底。研讨机器人动力学的目的研讨机器人动力学的目的5.1 5.1 机器人静力学机器人静力学 机器人静力学研讨机器人静止或者缓慢运动时作机器人静力学研讨机器人静止或者缓慢运动时作用在手臂上的力和力矩问题，特别是当手端与外界环用在手臂上的力和力矩问题，特别是当手端与

6、外界环境有接触力时，各关节力矩与接触力的关系。境有接触力时，各关节力矩与接触力的关系。 以下图表示作用在机器人手臂杆件以下图表示作用在机器人手臂杆件i i上的力和力矩上的力和力矩。其。其i-1fii-1fi为杆件为杆件i-1i-1对杆对杆i i的作用力，的作用力，-ifi+1-ifi+1为杆为杆i+1i+1对杆对杆i i的作用力，的作用力，i-1Nii-1Ni为杆件为杆件i-1i-1对杆对杆i i的作用力矩，的作用力矩，- -iNi+1iNi+1为杆为杆i+1i+1对杆对杆i i的作用力矩，的作用力矩，cici为杆为杆i i质心。质心。作用在杆作用在杆i i的的力和力矩力和力矩根据力、力矩平衡

7、原理有根据力、力矩平衡原理有5.2 5.2 机器人动力学正问题机器人动力学正问题 机器人动力学正问题研讨机器人手臂在关节力机器人动力学正问题研讨机器人手臂在关节力矩作用下的动态呼应。其主要内容是如何建立机器矩作用下的动态呼应。其主要内容是如何建立机器人手臂的动力学方程。建立机器人动力学方程的方人手臂的动力学方程。建立机器人动力学方程的方法有牛顿法有牛顿欧拉法和拉格朗日法等。欧拉法和拉格朗日法等。1、牛顿、牛顿欧拉法方程欧拉法方程 在思索速度与加在思索速度与加速度影响的情况下，速度影响的情况下，作用在机器人手臂杆作用在机器人手臂杆i上的力和力矩如右图上的力和力矩如右图所示。其中所示。其中vci和

8、和i分分别为杆别为杆i质心的平移速质心的平移速度向量和此杆的角速度向量和此杆的角速度向量。度向量。 根据力、力矩平衡根据力、力矩平衡原理有原理有:5-15-2称称5-1为牛顿方程，为牛顿方程，5-2为欧拉方程。为欧拉方程。其中其中IiIi为杆为杆i i绕其质心的惯性张量绕其质心的惯性张量2 2、 拉格朗日方程拉格朗日方程 牛顿一欧拉运动学方程是基于牛顿第二定律牛顿一欧拉运动学方程是基于牛顿第二定律和欧拉方程，利用达朗伯原理，将动力学问题变和欧拉方程，利用达朗伯原理，将动力学问题变成静力学问题求解。该方法计算快。拉格朗日动成静力学问题求解。该方法计算快。拉格朗日动力学那么是基于系统能量的概念，以

9、简单的方式力学那么是基于系统能量的概念，以简单的方式求得非常复杂的系统动力学方程，并具有显式构求得非常复杂的系统动力学方程，并具有显式构造，物理意义比较明确。造，物理意义比较明确。 (1) (1) 拉格朗日函数拉格朗日函数 对于任何机械系统，拉格朗日函数对于任何机械系统，拉格朗日函数L L定义为系统总的定义为系统总的动能动能EkEk与总的势能与总的势能EpEp之差，即之差，即: :( , )( , )( )kpL q qEq qEq12nqqqq12nqqqq表示动能与势能的广义坐标表示动能与势能的广义坐标相应的广义速度相应的广义速度 (2) (2) 机器人系统动能机器人系统动能 在机器人中，

10、连杆是运动部件，连杆在机器人中，连杆是运动部件，连杆i i的动能的动能EkiEki为连杆质心线速度引起的动能和连杆角速度产生的动为连杆质心线速度引起的动能和连杆角速度产生的动能之和，即能之和，即: :1122TiTiikiiciciiiiEmI 1nkkiiEE1( , )( )2TkE q qq D q q系统的动能为系统的动能为n n个连杆的动能之和，即：个连杆的动能之和，即： 由于由于 和和 是关节变量是关节变量 和关节速度和关节速度 的函数，因此，从上式可知，机器人的动能是关节变的函数，因此，从上式可知，机器人的动能是关节变量和关节速度的标量函数，记为量和关节速度的标量函数，记为 ，可

11、表示成，可表示成：1( , )( )2TkE q qq D q qciiiqq ( , )kE q q 式中，式中， 是是nxnnxn阶的机器人惯性矩阵阶的机器人惯性矩阵( )D q 3机器人系统势能机器人系统势能 设连杆设连杆i的势能为的势能为 ，连杆，连杆i的质心在的质心在O坐标系中的位置矢坐标系中的位置矢量为量为 ，重力加速度矢量在坐标系中为，重力加速度矢量在坐标系中为g，那么：，那么： 机器人系统的势能为各连杆的势能之和，即：机器人系统的势能为各连杆的势能之和，即： 它是它是q的标量函数。的标量函数。TpiiciEm g p 1npipiiEEpiEcip 4拉格朗日方程拉格朗日方程

12、系统的拉格朗日方程为：系统的拉格朗日方程为： 上式又称为拉格朗日上式又称为拉格朗日欧拉方程，简称欧拉方程，简称LE方程。式方程。式中，中， 是是n个关节的驱动力或力矩矢量，上式可写成：个关节的驱动力或力矩矢量，上式可写成：pkkEEEddtqqqdLLdtqq 例例平面平面RP机器人如下图，连杆机器人如下图，连杆l和连杆和连杆2的质量分别为的质量分别为m1和和m2，质心的位置由，质心的位置由l1和和d2所规定，惯量矩阵为：所规定，惯量矩阵为：1111100000 xxyyzzIIIiI2222200000 xxyyzzIIIiI (1) 取坐标，确定关节变量和驱动力或力矩取坐标，确定关节变量和

13、驱动力或力矩 建立连杆建立连杆D-H坐标系如上图所示，关节变量为坐标系如上图所示，关节变量为1+/2为求为求解方便，此处取关节变量为解方便，此处取关节变量为1和和d2，关节驱动力矩，关节驱动力矩l和力和力f2。(2)(2)系统动能系统动能 由式由式(1)(1)，分别得，分别得22211 111 11122kyyEmlI2222222122111()22kyyEm ddI22221 1122212211()22kyyyyEmlIIm dm d1122TiTiikiiciciiiiEmI 1总动能为：总动能为：(3)(3)系统势能系统势能 由于：由于：00Tgg11 11 10Tcpl cl s1

14、1111 1TpcEm g pm gl s 22222 1TpcEm g pm gd s 1 1221()pEg mlm ds那么：那么：总势能为：总势能为：221 11222122()yyyykmlIIm dEqm d222 10kEm dq1 12212 1()pEg mlm d cgm sq(4)(4)偏导数偏导数2211 11222122 121 1221222222 121()2()yyyymlIIm dm ddg mlm d cm dm dm gs(5)(5)拉格朗日动力学方程拉格朗日动力学方程 将偏导数代入拉格朗日将偏导数代入拉格朗日方程，得到平面方程，得到平面RPRP机器人的机

15、器人的动力学方程的封锁方式：动力学方程的封锁方式：pkkEEEddtqqq拉格朗日方程拉格朗日方程2( )( , )( )D q qH q qG q221 1122220( )0yyyymlIIm dD qm22 12222 12( , )m ddH q qm d1 122121()( )mlm dgcG qm gs (1) (1)关节空间动力学方程关节空间动力学方程 将式将式2 2写成矩阵方式：写成矩阵方式：3、关节空间与操作空间动力学、关节空间与操作空间动力学式中式中3 式式(3)为机器人在关节空间中的动力学方程封锁方式的普通构造式。为机器人在关节空间中的动力学方程封锁方式的普通构造式。它

16、反映了关节力或力矩与关节变量、速度和加速度之间的函数关系。对于它反映了关节力或力矩与关节变量、速度和加速度之间的函数关系。对于n个关节的机器人，个关节的机器人， 是是nn正定对称矩阵，是正定对称矩阵，是q的函数，称为机器人惯性的函数，称为机器人惯性矩阵；矩阵； 是是n1的离心力和哥氏力向量；的离心力和哥氏力向量； 是是n1重力矢量，与机器人重力矢量，与机器人的形位的形位q有关。有关。( )D q( , )H q q ( )G q( )( , )( )xxxFMq xUq qG q 2 2操作空间动力学方程操作空间动力学方程 与关节空间动力学方程相对应，在笛卡尔操作空间中，与关节空间动力学方程相

17、对应，在笛卡尔操作空间中，操作力操作力F F与末端加速度与末端加速度 之间的关系可表示为：之间的关系可表示为：x ( )xMq( , )xUq q ( )xG q操作空间中的惯性矩阵操作空间中的惯性矩阵离心力和哥氏力矢量离心力和哥氏力矢量重力矢量重力矢量xF广义操作力矢量广义操作力矢量机器人末机器人末端位姿向量端位姿向量( )TJq F( )xJ q q( )( )( )( , )rxJ q qJ q qJ q qa q q( )( )( ) ( )TxD qJq Mq J q( , )( )( , )( ) 9 )( , )TTxrH q qJq Uq qJqq a q q( )( )( )TxG qJq G q由上一章可知，广义操作力和关节力之间的关系为：由上一章可知，广义操作力和关节力之间的关系为：操作空间与关节空间之间的速度与加速度的关系：操作空间与关节空间之间的速度与加速度的关系：比较关节空间与操作空间动力学方程，可以得到：比较关节空间与操作空间动力学方程，可以得到：( )( )( , )( )TxxxJq Mq