二自由度机械臂动力学分析

平面二自由度机械臂动力学分析姓名：黄辉龙专业年级：13级机电单位：汕头大学摘要：机器臂是一个非线性的复杂动力学系统。动力学问题的求解比较困难，而且需要较长的运算时间，因此，这里主要对平面二自由度机械臂进行动力学研究。拉格朗日方程在多刚体系统动力学的应用方法分析平面二自由度机械臂的正向动力学。经过分析，得出平面二自由度机械臂的动力学方程，为后续更深入研究做铺垫。关键字：平面二自由度动力学方程拉格朗日方程相关介绍机器人动力学的研究有牛顿-欧拉（Newton-Euler）法、拉格朗日(Langrange)法、高斯（Gauss）法等，但一般在构建机器人动力学方程中，多采用牛顿-欧拉法及拉格朗日法。欧拉方程又称牛顿-欧拉方程，应用欧拉方程建立机器人机构的动力学方程是指研究构件质心的运动使用牛顿方程，研究相对于构件质心的转动使用欧拉方程，欧拉方程表征了力、力矩、惯性张量和加速度之间的关系。在机器人的动力学研究中，主要应用拉格朗日方程建立机器人的动力学方程，这类方程可直接表示为系统控制输入的函数，若采用齐次坐标，递推的拉格朗日方程也可以建立比较方便且有效的动力学方程。在求解机器人动力学方程过程中，其问题有两类：1)给出已知轨迹点上、及、，即机器人关节位置、速度和加速度，求相应的关节力矩矢量。这对实现机器人动态控制是相当有用的。2)已知关节驱动力矩，求机器人系统相应各瞬时的运动。也就是说，给出关节力矩矢量，求机器人所产生的运动、及、。这对模拟机器人的运动是非常有用的。平面二自由度机械臂动力学方程分析及推导过程1、机器人是结构复杂的连杆系统，一般采用齐次变换的方法，用拉格朗日方程建立其系统动力学方程，对其位姿和运动状态进行描述。机器人动力学方程的具体推导过程如下：1)选取坐标系，选定完全而且独立的广义关节变量nr,,2,1,r。2)选定相应关节上的广义力rF：当r是位移变量时，rF为力；当r是角度变量时，rF为力矩。3)求出机器人各构件的动能和势能，构造拉格朗日函数。4)代入拉格朗日方程求得机器人系统的动力学方程。2、下面以图1所示说明机器人二自由度机械臂动力学方程的推导过程。1）如图1，设21,是广义坐标，21,QQ是广义力。2）分别求出两杆的动能和势能11112111111sin,2121:1cccTcglmUIvvmE杆（1-1）]sin[,2121:22112222122222）（）（杆lgmUIvvmEccTc（1-2）式中，1cv是杆1质心),(111ccyxC的速度向量，2cv是杆2质心),(222ccyxC的速度向量。它们可以根据质心21,CC的位置方程导出。3）分别求出两杆的速度)sin()cos(1111111ldtdldtddtdydtdxvccc（1-3）)sin(sin)cos(cos2121121211222ccccclldtdlldtddtdydtdxv（1-4）4）代入拉格朗日方程求得机械臂动力学方程根据具有完整理想约束的有N个广义坐标系统的拉格朗日方程nrQqUqEqEdtdrrrr,2,1,（1-5）式中，rq是第r个广义坐标，E是系统动能，U是系统势能，rQ是对第r个广义坐标的广义力。该问题为二自由度的动力学研究，所以n=2，由于势能函数U与广义速度无关，即rqU=0。由（1-5）式可写成：rrrQqLqLdtd（1-6）其中，是拉格朗日算子LUEL,可知在这里拉格朗日算子为：2121UUEEL代入式（1-6）可导出相应的式子，经过整理得：QgCM)(),()(（1-7）式中21212122211211,)(,),(,)(QQQgggCCCMMMMM)cos()cos(coscossin),2(sin,)cos()cos2(21222212112111122122212221212222221221222122221222122212121111cccccccccccccglmgllgmglmgllmCllmCIlmMMMIlllmMIllllmIlmM式（1-7）为机械臂在关节空间的动力学方程的一般结构形式，它反映了关节力矩和关节变量、速度、加速度之间的函数关系。对于n个关节的机械臂，)(M是nn的正定对称矩阵，是θ的函数，称为机械臂的惯性矩阵；),(C是nx1的离心力和科氏力矢量；)(g是nx1的重力矢量，与机械臂的形位θ有关。