第七讲机器人动力学

第七讲机器人动力学第1页，课件共78页，创作于2023年2月第四章机器人动力学9/7/2023动力学研究的问题：机器人各个关节的运动与关节需要的驱动力（矩）之间的关系。正问题：已知关节运动，求关节驱动力（矩）。逆问题：已知关节驱动力（矩），求关节运动。第2页，课件共78页，创作于2023年2月数学模型：关节运动→位移、速度、加速度变化→关节驱动力（矩）→驱动力或驱动力矩→τi动力学方程：

，i=1，…，n正问题：已知

，求τi。逆问题：已知τi，求。第四章机器人动力学9/7/2023第3页，课件共78页，创作于2023年2月4.1动力学模型1、力学分析2、拉格朗日方程法3、动力学模型9/7/2023第4页，课件共78页，创作于2023年2月4.1动力学模型1、力学分析9/7/2023（1）静力学分析机器人各个关节处于静止状态。当负载为一重物时：关节承受的力和力矩：关节需要的驱动力（矩）：mgf3=mgf2=mgf1=mgm2=mgl2l2l1m1=mg(l1+l2)τ1=0τ2=mgl2τ3=mg第5页，课件共78页，创作于2023年2月4.1动力学模型1、力学分析9/7/2023（1）静力学分析机器人各个关节处于静止状态。考虑杆件自重时：关节承受的力和力矩：关节需要的驱动力（矩）：mgf3=mgf2=mgf1=mgm2=mgl2l2l1m1=mg(l1+l2)τ1=0τ2=mgl2τ3=mgm3gm2gm1g第6页，课件共78页，创作于2023年2月4.1动力学模型1、力学分析9/7/2023（2）动力学分析机器人各个关节处于运动状态。当负载为一重物时：关节承受的力和力矩：关节需要的驱动力（矩）：f3f2f1m2l2l1m1τ1τ2τ3m3第7页，课件共78页，创作于2023年2月4.1动力学模型2、拉格朗日方程法9/7/2023拉格朗日方程的一般形式为：

式中，——广义力，它可以是力，也可以是力矩； ——系统选定的广义坐标； ——广义坐标对时间的一阶导数，即速度；——拉格朗日函数，又称为拉格朗日算子，它被定义为系统的动能与势能之差L=T-U。第8页，课件共78页，创作于2023年2月4.1动力学模型2、拉格朗日方程法9/7/2023

对给定的机器人，可以按以下几个步骤建立拉格朗日动力学方程：（1）选取完全并独立的广义坐标；（2）选定广义力；（3）求出系统的动能T和势能U，并用其构造拉格朗日函数L=T-U；（4）将以上结果代入拉格朗日方程式中，即可求得机器人的动力学方程。第9页，课件共78页，创作于2023年2月4.1动力学模型2、拉格朗日方程法9/7/2023例：已知二关节机器人如图所示，机器人的两个连杆长度分别为l1和l2，质量分别为m1和m2，且集中在各连杆的端部。若将机器人直接悬挂在加速度为g的重力场中，试用拉格朗日方程建立该机器人的动力学方程。解：①选取连杆绕关节的转角为变量θ1和θ2

，则系统的广义坐标就可以选为，即②转动关节对应的是力矩，所以广义力就选为，即。θ1θ2m1m2xy关节1关节2第10页，课件共78页，创作于2023年2月4.1动力学模型2、拉格朗日方程法9/7/2023求出各连杆的动能和势能：连杆l1的动能为：连杆l1的势能为：对连杆l2求动能和势能时，要先写出其质心在直角坐标系中的位置表达式：然后求微分，则其速度就为：由此可得连杆的速度平方值为：第11页，课件共78页，创作于2023年2月4.1动力学模型2、拉格朗日方程法9/7/2023求出各连杆的动能和势能：从而连杆l2的动能为：势能为：则可构造出拉格朗日函数为：

第12页，课件共78页，创作于2023年2月4.1动力学模型2、拉格朗日方程法9/7/2023求出机器人动力学方程：先将拉格朗日函数对和进行微分，即：第13页，课件共78页，创作于2023年2月4.1动力学模型2、拉格朗日方程法9/7/2023求出机器人动力学方程：再将拉格朗日函数对和进行微分，即：第14页，课件共78页，创作于2023年2月4.1动力学模型2、拉格朗日方程法9/7/2023求出机器人动力学方程：将以上结果代入方程即可得关节上的力矩分别为：

第15页，课件共78页，创作于2023年2月4.1动力学模型3、动力学模型9/7/2023将得到的机器人动力学方程简写为如下形式：当机器人有n个关节时，上式可推广为普遍形式：

第16页，课件共78页，创作于2023年2月4.1动力学模型3、动力学模型9/7/2023将上式进一步简化为如下所示的矩阵形式：

上式也称为机器人的动力学模型。式中：是机器人动力学模型中的惯性力项，表示机器人操作机的质量矩阵，它是n×n阶的对称矩阵，。

是n×1阶矩阵，表示机器人动力学模型中非线性的耦合力项，包括离心力（自耦力）和哥氏力（互耦力）。也是n×1阶矩阵，表示机器人动力学模型中的重力项。

第17页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法9/7/20231、牛顿——欧拉方程2、递推计算公式3、递推算法应用第18页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法9/7/2023牛顿——欧拉方程原理：将机器人的每个杆件看成刚体，并确定每个杆件质心的位置和表征其质量分布的惯性张量矩阵。当确定机器人坐标系后，根据机器人关节速度和加速度，则可先由机器人机座开始向手部杆件正向递推出每个杆件在自身坐标系中的速度和加速度，再用牛顿——欧拉方程得到机器人每个杆件上的惯性力和惯性力矩，然后再由机器人末端关节开始向第一个关节反向递推出机器人每个关节上承受的力和力矩，最终得到机器人每个关节所需要的驱动力（矩），这样就确定了机器人关节的驱动力（矩）与关节位移、速度和加速度之间的函数关系，即建立了机器人的动力学方程。第19页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法9/7/2023牛顿——欧拉方程递推过程：正向递推：已知机器人各个关节的速度和加速度→从1～n递推出机器人每个杆件在自身坐标系中的速度和加速度→机器人每个杆件质心上的速度和加速度→再用牛顿——欧拉方程得到机器人每个杆件质心上的惯性力和惯性力矩。反向递推：根据正向递推的结果→从n～1递推出机器人每个关节上承受的力和力矩→得到机器人每个关节所需要的驱动力（矩）。第20页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法9/7/20231、牛顿——欧拉方程（1）牛顿方程→惯性力①矢量。②质心上的线加速度。第21页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法9/7/20231、牛顿——欧拉方程（2）欧拉方程→惯性力矩①矢量。②质心上的惯性张量矩阵。第22页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法9/7/20231、牛顿——欧拉方程（2）欧拉方程→惯性力矩惯性张量矩阵：a.坐标系：与杆件坐标系同向，位于杆件质心上。b.元素名称：Icxx，Icyy，Iczz——惯性矩；Icxy=Icyx,Icyz=Iczy,Iczx=Icxz——惯性积。第23页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法9/7/20231、牛顿——欧拉方程（2）欧拉方程→惯性力矩惯性张量矩阵：理论计算方法：实验测试法：

惯量摆仪器。第24页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法9/7/20232、递推计算公式（1）正向递推：已知机器人各个关节的速度和加速度：①从1～n递推出机器人每个杆件在自身坐标系中的速度和加速度；②机器人每个杆件质心上的速度和加速度；③机器人每个杆件质心上的惯性力和惯性力矩。第25页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法9/7/20232、递推计算公式（1）正向递推：以第一种杆件坐标系为例①杆件速度和加速度递推计算公式建立相邻两个杆件的坐标系：{i-1}、{i}。Oiii-1关节iXi-1Zi-1Oi-1XiZi第26页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法2、递推计算公式（1）正向递推：①杆件速度和加速度递推计算公式已知：i关节速度和加速度i-1杆件速度和加速度计算：i杆件速度和加速度9/7/2023ii-1关节iXi-1Zi-1Oi-1XiZiOi第27页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法2、递推计算公式（1）正向递推：①杆件速度和加速度递推计算公式分析：I、坐标系:相邻杆件位姿矩阵II、关节速度和加速度的矢量化：9/7/2023ii-1关节iXi-1Zi-1Oi-1XiZiOi第28页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法2、递推计算公式（1）正向递推：①杆件速度和加速度递推计算公式9/7/2023ii-1关节iXi-1Zi-1Oi-1XiZiOi第29页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法2、递推计算公式（1）正向递推：①杆件速度和加速度递推计算公式9/7/2023ii-1关节iXi-1Zi-1Oi-1XiZiOi第30页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法2、递推计算公式（1）正向递推：①杆件速度和加速度递推计算公式9/7/2023关节iii-1Xi-1Zi-1Oi-1XiZiOi第31页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法2、递推计算公式（1）正向递推：②杆件质心上的速度和加速度9/7/2023iXiZiOi第32页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法2、递推计算公式（1）正向递推：③杆件质心上的惯性力和惯性力矩惯性力：惯性力矩：9/7/2023iXiZiOi第33页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法9/7/20232、递推计算公式（1）正向递推：以第二种杆件坐标系为例①杆件速度和加速度递推计算公式建立相邻两个杆件的坐标系：{i-1}、{i}。ii-1关节iXi-1Zi-1Oi-1XiZiOi第34页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法2、递推计算公式（1）正向递推：①杆件速度和加速度递推计算公式已知：i关节速度和加速度i-1杆件速度和加速度计算：i杆件速度和加速度9/7/2023ii-1关节iXi-1Zi-1Oi-1XiZiOi第35页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法2、递推计算公式（1）正向递推：①杆件速度和加速度递推计算公式分析：I、坐标系:相邻杆件位姿矩阵II、关节速度和加速度的矢量化：9/7/2023ii-1关节iXi-1Zi-1Oi-1XiZiOi第36页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法2、递推计算公式（1）正向递推：①杆件速度和加速度递推计算公式9/7/2023ii-1关节iXi-1Zi-1Oi-1XiZiOi第37页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法2、递推计算公式（1）正向递推：①杆件速度和加速度递推计算公式9/7/2023ii-1关节iXi-1Zi-1Oi-1XiZiOi第38页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法2、递推计算公式（1）正向递推：②杆件质心上的速度和加速度9/7/2023iXiZiOi第39页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法2、递推计算公式（1）正向递推：③杆件质心上的惯性力和惯性力矩惯性力：惯性力矩：9/7/2023iXiZiOi第40页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法9/7/20232、递推计算公式（2）反向递推：已知机器人各个杆件的惯性力和惯性力矩：①从n～1递推出机器人每个关节承受的力和力矩；②机器人每个关节的驱动力或驱动力矩。第41页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法9/7/20232、递推计算公式（2）反向递推：以第二种杆件坐标系为例①关节承受的力和力矩递推计算公式建立相邻两个杆件的坐标系：{i-1}、{i}ii-1关节iXi-1Zi-1Oi-1XiZiOi第42页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法2、递推计算公式（2）反向递推：①关节承受的力和力矩递推计算公式已知：i-1杆件的惯性力和惯性力矩i关节承受的力和力矩计算：i-1关节承受的力和力矩9/7/2023ii-1关节iXi-1Zi-1Oi-1XiZiOi第43页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法2、递推计算公式（2）反向递推：①关节承受的力和力矩递推计算公式分析：I、坐标系:相邻杆件位姿矩阵II、i-1杆件受力分析：9/7/2023ii-1关节iXi-1Zi-1Oi-1XiZiOi第44页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法2、递推计算公式（2）反向递推：①关节承受的力和力矩递推计算公式以i-1杆件为研究对象，由达朗贝尔原理可得：9/7/2023i-1关节iXi-1Zi-1Oi-1Oi第45页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法2、递推计算公式（2）反向递推：①关节承受的力和力矩递推计算公式以i-1杆件为研究对象，由达朗贝尔原理可得：9/7/2023i-1关节iXi-1Zi-1Oi-1Oi第46页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法2、递推计算公式（2）反向递推：②关节驱动力（矩）平移关节：回转关节：则关节驱动力（矩）为：9/7/2023i-1关节iXi-1Zi-1Oi-1第47页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法3、递推算法应用（1）递推初始条件①正向递推机座0的速度和加速度：*考虑杆件自重或手部负载为重物时：为描述在机座坐标系{0}中的标准重力加速度。9/7/2023gx0z0o0第48页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法3、递推算法应用（1）递推初始条件②反向递推机器人手部负载：9/7/2023第49页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法3、递推算法应用（2）递推应用条件①已知机器人的关节变量及其速度和加速度；②已知任一杆件i相对于与自身坐标系{i}方向相同的坐标系{Ci}所描述的惯性张量及其质心在自身坐标系{i}中的位置矢量（可用实验等方法确定）；③已知相邻杆件的位姿矩阵及必要的初始数据。9/7/2023第50页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法3、递推算法应用（2）递推应用条件第二种坐标系下递推算法正向递推:9/7/2023第51页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法3、递推算法应用（2）递推应用条件第二种坐标系下递推算法反向递推:9/7/2023第52页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法3、递推算法应用例：已知二自由度机器人如图所示，机器人两个杆件的长度分别为和，且其质量和都集中在杆件的端头。若用第二种方法建立机器人的坐标系，当机器人各个关节的位移、速度和加速度已知时，试用牛顿——欧拉递推算法计算各关节的驱动力矩。9/7/2023θ1θ2m1m2xy关节1关节2第53页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法3、递推算法应用解：建立坐标系如图所示。相邻杆件的位姿矩阵为：9/7/2023θ1θ2m1m2x0y0关节1关节2x1x2y1y2第54页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法3、递推算法应用解：(1)正向递推已知关节速度和加速度分别为，由于考虑杆件的重量，所以机座的运动参数（初始条件）设为：9/7/2023θ1θ2m1m2x0y0关节1关节2x1x2y1y2g第55页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法3、递推算法应用解：(1)正向递推9/7/2023θ1θ2m1m2x0y0关节1关节2x1x2y1y2g第56页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法3、递推算法应用解：(1)正向递推9/7/2023θ1θ2m1m2x0y0关节1关节2x1x2y1y2g第57页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法3、递推算法应用解：(1)正向递推9/7/2023θ1θ2m1m2x0y0关节1关节2x1x2y1y2g第58页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法3、递推算法应用解：(1)正向递推9/7/2023θ1θ2m1m2x0y0关节1关节2x1x2y1y2g第59页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法3、递推算法应用解：(1)正向递推9/7/2023θ1θ2m1m2x0y0关节1关节2x1x2y1y2g第60页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法3、递推算法应用解：(1)正向递推9/7/2023θ1θ2m1m2x0y0关节1关节2x1x2y1y2g第61页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法3、递推算法应用解：(1)正向递推9/7/2023θ1θ2m1m2x0y0关节1关节2x1x2y1y2g第62页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法3、递推算法应用解：(1)正向递推9/7/2023θ1θ2m1m2x0y0关节1关节2x1x2y1y2g第63页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法3、递推算法应用解：(1)正向递推9/7/2023θ1θ2m1m2x0y0关节1关节2x1x2y1y2g第64页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法3、递推算法应用解：(1)正向递推9/7/2023θ1θ2m1m2x0y0关节1关节2x1x2y1y2g第65页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法3、递推算法应用解：(1)正向递推9/7/2023θ1θ2m1m2x0y0关节1关节2x1x2y1y2g第66页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法3、递推算法应用解：(1)正向递推9/7/2023θ1θ2m1m2x0y0关节1关节2x1x2y1y2g第67页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法3、递推算法应用解：(2)反向递推由于机器人手部无负载，所以初始条件为：9/7/2023θ1θ2m1m2x0y0关节1关节2x1x2y1y2g第68页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法3、递推算法应用解：(2)反向递推关节2受的力和力矩：9/7/2023θ1θ2m1m2x0y0关节1关节2x1x2y1y2g第69页，课件共78页，创作于2023年2月4.2牛顿——欧拉方程法3、递推算法应用解：(2)反向递推关节2受的力和力矩：9/7/2