1、s analytical solution fast and precisely,anopti malm ethod of setting the coordinate system has been put for ward according to the theory of a robot?s4-DOF m e-chanical ar m.The kine m atics and inverse kinematics equations of the 4-DOF mechanical ar m were establishedbased on this coordinate system
2、.The positions of itsmoving elements and end-effectorsw ith their relationshipwerepresented;thesemade the robotmove in accordance w ith a scheduled program.Finally,the Jacobian matrix wassolved,providing theoretical basis for realizing themechanical ar m?s program design of velocity control in theCa
3、rte-sian space.Key words:mechanical ar m;D-H parameters;kinematics;inverse kinematics;Jacobian matrix?机器人运动学着重研究了机器人各个坐标系之间的运动关系。为了控制机器人的运动,首先需要在机器人中建立相应的坐标系 1。传统建立机械手臂坐标系的方法是将固定在大地上的笛卡尔坐标系作为参考坐标系 2,3。这种方法容易引入不必要的参数,增加计算量,引起计算误差增大,使机器人不能够精确控制。特别是在解机器人逆运动学方程中,因为引入了多余的参数,使方程复杂度增加,甚至会引起无解或者无穷多解的情况 4 7,
4、需要通过观察和经验判断才能得出正确的解,既给运动学分析带来了不便,又增加了机器人控制程序编写难度 8 11。笔者采用优化的坐标系建立方法,减少了参数的数量。通过这种方法,在求解正、逆运动学方程过程中,可以减少运算公式的参数计算,降低计算难度,求出逆运动学方程封闭的解析解,从而能够更好地建立机器人各运动构件与末端执行器在空间的位置、姿态之间的关系。并方便的求解出雅可比矩阵,为实现机械手臂末端在笛卡尔空间的速度控制程序第 6期吴?磊等:四自由度机械手臂运动学分析及雅可比矩阵求解设计提供理论依据。1?机械手臂运动学分析1?正向运动学分析该机械手臂由 4个转动关节和 3个连杆组成,终端有一个夹持器。用
5、 D-H 方法建立杆件坐标系。因为无论转动关节 0如何旋转,关节 0的长度对机械手臂的运动均无影响,所以令基准坐标系 0与坐标系 1原点重合,取在关节 0与关节 1的交点处。原点 O2位于关节 2的中心位置,Z2轴沿着关节 2的方向向右,X2轴垂直于关节 1和关节 2构成的平面指向上方。坐标系 3与坐标系 4原点重合。X3轴垂直于关节 2和关节 3构成的平面指向上方,X4轴垂直于夹持器所在平面,如图 1所示。图 1?机械手臂坐标示意图本文的基坐标系是 0坐标系,如果0A1表示第一个连杆坐标系相对于基坐标系的齐次变换矩阵,1A2表示第二个连杆坐标系相对于第一个坐标系的齐次变换矩阵,那么齐次变换矩
6、阵0A2可由0A1和1A2乘积给出:0A2=0A11A2。同理,齐次变换矩阵0A4为:0A4=0A11A22A33A4 2。当机械手臂全部的连杆坐标系建立好以后,就能够列出各连杆的常量参数,各连杆 D-H 参数及关节变量,如表 1所示。表 1?机械手臂 D-H参数表iai?idi?i?i范围10-90?0?1(-?/2,?/2)2a200?2-3?/4,-?/43a300?3-?/4,?/64000?4 0,?/2)表 1中机械手臂 D-H 参数的含义如下:ai定义为从 Zi-1轴到 Zi轴的距离,沿 Xi轴的指向为正;?i定义为从 Zi-1轴到 Zi轴的转角,绕Xi轴的正向转动为正;di定义
7、为从 Xi-1轴到 Xi轴的距离,沿 Zi-1轴的指向为正;i定义为从 Xi-1轴到 Xi轴的转角,绕Zi-1轴的正向转动为正。从图 1和 D-H 参数表,可以看出,系 i-1可以经过以下连续的相对运动变到系i:第 1步,沿Zi-1轴移动 di;第 2步,绕 Zi-1轴转动?i;第 3步,沿Xi轴移动 ai;第 4步,绕Xi轴转动?i。这样矩阵i-1Ai就成为关节变量?i的函数,一旦求得这些参数,就能确定 4个变换矩阵i-1Ai的值。由连续相对运动时齐次变换矩阵的求法,可得i-1Ai=Transz(di)Rotz(?i)T ransx(ai)Rotx(?i)=ci-c?isis?isiaici
8、sic?ici-s?iciaisi0s?ic?idi0001(1)式中:T ransx(di)为沿 Z轴移动 di;Rotz(?i)为绕 Z轴转动?T ransx(ai)为沿 X 轴移动 ai;Rotx(?i)为绕X 轴转动?i,si=sin?i,ci=cos?i,s?i=sin?i,c?i=cos?由式(1),可以得出机械手臂正运动学方程为0A4=0A11A22A33A4=c1c234-c1s234-s1a3c1c23+a2c1c2s1c234-s1s234c1a3s1c23+a2s1c2-s234-c2340-a3s23-a2s20001=nxoxaxpxnyoyaypynzozazpz0
9、001(2)式中:n=nx,ny,nzT、o=ox,oy,ozT、a=ax,ay,azT和 p=px,py,pzT分别为机器人末端坐标系 4相对于基坐标系 0位姿的 4个列向量;si j=sin(?i+?j),ci j=cos(?j),si j k=sin(?j+?k),ci j k=cos(?k),(i=2,3,4;j=2,3,4)。式(2)表示该机械手臂的正运动学变换矩阵,描述了连杆末端坐标系 4相对于基坐标系 0的位姿,是机械臂运动分析和综合的基础。2?逆向运动学分析机器人的逆向运动学,描述的就是由机器人末端的位置和姿态到机器人关节坐标系的坐标之间的映射关系。通过给出的机械手臂运动学方程
10、求解出各个关节角?765机 械 科 学 与 技 术第 28卷(1)求解?1由式(2)ax=-sin?1,ay=cos?1,可得axay=-sin?1cos?1=-tan?1,?1=-arctanaxay(3)?(2)求解?2,?3由式(2)pz=-a3sin(?2+?3)-a2sin?2,px=a3cos?1cos(?3)+a2cos?2,ay=cos?1,可得cos?3=pxay2+(pz)2-a23-a22a2ax(4)?由式(2),得到1A-20A-10A4=2A33A4(5)?将等式两端展开,得到1A-20A-10A4=1A-20A-1nxoxaxpxnyoyaypynzozazpz0
11、001=t11t12t13pxcos?2+pysin?2-pzsin?2-k2t21t22t23-pxcos?1sin?2-pysin?2-pzcos?2t31t32t33-pxsin?1+pycos?10001(6)2A4=2A33A4=cos(?3+?4)-sin(?4)0a3cos?3sin(?4)cos(?4)0a3sin?300100001(7)?由于只关心矩阵中的第四列,所以在式(6)中未列出具体的计算结果。利用式(6)和式(7)对应元素相等,则得cos?2px+sin?2py-sin?2pz-k2=a3cos?3(8)-cos?2px-sin?2py-cos?2pz=a3sin?
12、3(9)?因为,角度?2的范围是 -3?/4,则有sin?2=-1-cos2?2(10)?再将式(4)和式(10)代入式(8)可得cos?2(aypx-axpy)+1-cos2?2pz=a3cos?3+a2(11)?令 k1=aypx-axpy,k2=pz,k3=a3cos?3+a2,可得 cos?2的解式为cos?2=-12k2k1(k2?k22+4k21-4k1k3)-k3k1(12)故有?2=arccos-12k2k1(k2?k22+4k21-4k1k3)-k3k1-?(13)?由式(2)已知,ax=-sin?1,将式(13)代入式(9)中,并除以式(8)有tan?3=(axpy-ayp
13、x)sin?2(aypx-axpy)cos?2-a2(14)?由此可求得?3为?3=arctan(axpy-aypx)sin?2-a2(15)?(3)求解?4由式(2)nz=-sin(?4),oz=-cos(?4),可得nzoz=tan(?4)(16)所以,则有?4=arctannzoz-?2-?3(17)?通过上述分析,即可获得机械手臂 3个关节变量?3和?4的大小。机械臂雅可比矩阵求解机械手的笛卡尔空间运动速度与关节空间运动速度之间的变换,称之为雅可比矩阵(Jacobian m a-trix)。雅可比矩阵是关节空间速度向笛卡尔空间速度的传动比。即V=X?=J(q)q?(18)式中:q?为关
14、节速度矢量;X?是操作速度矢量。雅可比矩阵的行数等于机器人在操作空间运动的维数,列数等于关节数。因此,四自由度机械手臂雅可比矩阵是一个 6?4矩阵。构造雅可比矩阵的方法有矢量积法和微分变换法,本文用微分变换法推算该机器人的雅可比矩阵。雅可比矩阵第 i列的求解为Ji=(p?n)z,(p?o)z,(p?a)z,nz,oz,azT(19)766第 6期吴?四自由度机械手臂运动学分析及雅可比矩阵求解式中:n,o和 a矢量构成第 i连杆坐标系到机器人末端坐标系之间变换的旋转变换矩阵,p是第 i连杆坐标系到机器人末端坐标系之间变换的位置矢量。由于3A4=cos?4-sin?400sin?4cos?4000
15、0100001(20)1A4=1A22A33A4=cos(?4)0a3cos(?2sin(?4)0a3sin(?3)+a2sin?200100001?(21)?对于第 4个关节,式(20)就是该关节的连杆坐标系到末端坐标系的变换矩阵。将式(20)中的 n,o,a和 p 代入式(19),得到雅可比的列矢量 J4。J4=0,0,0,0,0,1T(22)?对于第 3个关节,将式(7)中的 n,o,a和 p代入式(19),得到雅可比的列矢量 J3。J3=0,0,0,0,0,1T(23)?对于第 2个关节,将式(21)中的 n,o,a和 p 代入式(19),得到雅可比的列矢量 J2。J2=0,0,0,0
16、,0,1T(24)?对于第 1个关节,将式(2)中的 n,o,a和 p代入式(19),得到雅可比的列矢量 J1。J1=sin(?4)(a3sin(?2)cos(?2)0-sin(?4)-cos(?4)0(25)故J=J1,J2,J3,J4?即为本机器人的雅可比矩阵。已知关节的速度即式(18)中的 q?,可以借助雅可比矩阵 J 求出操作空间的速度,即式(18)中的 X?。3?结论(1)这种优化的建立坐标系方法,减少了数学建模中的计算参数,简化了四自由度机械手臂正、逆运动学方程求解过程。(2)以此为依据,求出了雅可比矩阵,为仿真机械手臂的运动特性和设计控制器奠定了基础。(3)这种分析方法也适用于其
17、他类似的更多自由度关节型串联机械臂正、逆运动学求解以及雅可比矩阵求解。参考文献 1?谭民,徐德,侯增广等.先进机器人控制 M.高等教育出版社,2007 2?Lu G M,Sun L N,Zhang B.K inematics analysis of the 5-DOFrehabilitant mechanical ar m J.Journal ofHarbin Institute ofTechnology(New Series),2006,38(5):698 701 3?X ie J,Q iangW Y.Inverse kinematics problem for 6-DOF spaceman
18、ipulator based on the thery of screws A.Proceedings oft he 2007 IEEE InternationalConference on Robotics and B-iom i m etics C,Sanya.2007:1659 1663 4?Bilge G,et al.M echanical and electrical design of a four-whee-ldrive,four-whee-l steermobile manipulatorw ith PA-10 ar m J.Institute of Electrical an
19、d E lectron ics Engineers ComputerSociety,2007,2:1777 1782 5?D anevit J,Hartenberg R S.A kinematic notation for lower-pairmechanisns based on matrices J.AS M E Journal ofAppliedM echan ics,1995,22:215 221 6?宁柯军,杨汝清.一种串联机械臂通用位置逆解方法研究 J.机械科学与技术,2005,24(9):1066 1068 7?Siegwart R,Nourbakehsh I R.Introdu
20、ction to Autono m ousM obile RobotsM.The M IIPress,M assachusetts Institute ofTechnology,Cambridge,2004 8?Clark D.M ichael Owings:Building Robot Drive TrainsM.The M cGraw-H illCompanies,Inc,2003 9?霍伟.机器人动力学与控制 M.高等教育出版社,2005 10?姜哲,金奕江等译.人工智能一种现代方法(第二版)M.人民邮电出版社,2004 11?侯国柱,孔庆忠.一关节型机械臂运动学分析及雅可比矩阵求解 J.宁夏工程技术,2007,6(1):14 17767
copyright@ 2008-2023 冰点文库 网站版权所有
经营许可证编号:鄂ICP备19020893号-2