机器人正逆运动学

1、机器人正逆运动学课件第1页，共59页，2022年，5月20日，2点18分，星期四D-H表示法学习目标：1. 理解D-H法原理 2. 学会用D-H法对机器人建模学习重点：1. 给关节指定参考坐标系 2. 制定D-H参数表 3. 利用参数表计算转移矩阵第2页，共59页，2022年，5月20日，2点18分，星期四背景简介： 1955年，Denavit和Hartenberg(迪纳维特和哈坦伯格)提出了这一方法，后成为表示机器人以及对机器人建模的标准方法，应用广泛。 总体思想： 首先给每个关节指定坐标系，然后确定从一个关节到下一个关节进行变化的步骤，这体现在两个相邻参考坐标系之间的变化，将所有变化结合起

2、来，就确定了末端关节与基座之间的总变化，从而建立运动学方程，进一步对其求解。第3页，共59页，2022年，5月20日，2点18分，星期四1.第一个关节指定为关节n,第二个关节为n+1,其余关节以此类推。坐标系的确定2.Z轴确定规则：如果关节是旋转的，Z轴位于按右手规则旋转的方向，转角 为关节变量。如果关节是滑动的，Z轴为沿直线运动的方向，连杆长度d为关节变量。关节n处Z轴下标为n-1。第4页，共59页，2022年，5月20日，2点18分，星期四3.X轴确定规则情况1：两关节Z轴既不平行也不相交取两Z轴公垂线方向作为X轴方向，命名规则同Z轴。情况2：两关节Z轴平行此时，两Z轴之间有无数条公垂线，

3、可挑选与前一关节的公垂线共线的一条公垂线。情况3：两关节Z轴相交取两条Z轴的叉积方向作为X轴。4.Y轴确定原则取X轴、Z轴叉积方向作为Y轴方向。（右手）5.变量选择原则用n+1角表示Xn到Xn+1绕Zn轴的旋转角;dn+1表示从Xn到Xn+1沿Zn测量的距离;an+1表示关节偏移，an+1是从Zn到Zn+1沿Xn+1测量的距离;角表示关节扭转, n+1是从Zn到Zn+1绕Xn+1旋转的角度。 通常情况下，只有和d是关节变量。第5页，共59页，2022年，5月20日，2点18分，星期四斯坦福机器人斯坦福机器人开始的两个关节是旋转的，第三个关节是滑动的，最后三个腕关节全是旋转关节第6页，共59页，

4、2022年，5月20日，2点18分，星期四第7页，共59页，2022年，5月20日，2点18分，星期四例1：Stanford机器人运动学方程第8页，共59页，2022年，5月20日，2点18分，星期四A1A2A3A4A5A6d1z1x1y1O1d2z2x2y2O2z3y3x3O3y4z4x4O4z5y5x5O5d3z6x6y6O6d6z0y0 x0O0为右手坐标系原点Oi： Ai与Ai+1关节轴线的交点zi轴：与Ai+1关节轴重合，指向任意 xi轴： Zi和Zi-1构成的面的法线yi轴：按右手定则 ai沿 xi 轴， zi-1 轴与 xi 轴交点到Oi 的距离i 绕 xi 轴，由 zi-1 转

5、向zidi 沿 zi-1 轴，zi-1 轴和 xi 交点至Oi 1 坐标 系原点的距离i 绕 zi-1 轴，由 xi-1转向 xi关节1坐标系0关节2坐标系1关节3坐标系2连杆0连杆1连杆2连杆3连杆4连杆5关节4坐标系3关节5坐标系4关节6坐标系5第9页，共59页，2022年，5月20日，2点18分，星期四解：第10页，共59页，2022年，5月20日，2点18分，星期四第11页，共59页，2022年，5月20日，2点18分，星期四第12页，共59页，2022年，5月20日，2点18分，星期四第13页，共59页，2022年，5月20日，2点18分，星期四123456关节变量都是例2、PUMA

6、560运动学方程（六个自由度，全部是旋转关节）第14页，共59页，2022年，5月20日，2点18分，星期四PUMA560机器人的连杆及关节编号第15页，共59页，2022年，5月20日，2点18分，星期四为右手坐标系，Yi轴：按右手定则 Zi轴：与Ai+1关节轴重合，指向任意 Xi轴： Zi和Zi-1构成的面的法线， 或连杆i两端轴线Ai 与Ai+1的公垂线(即： Zi和Zi-1的公垂线)原点Oi： Ai与Ai+1关节轴线的交点，或Zi与Xi的交点ai沿 xi 轴， zi-1 轴与 xi 轴交点到Oi 的距离i 绕 xi 轴，由 zi-1 转向zidi 沿 zi-1 轴，zi-1 轴和 xi

7、 交点至Oi 1 坐标 系原点的距离i 绕 zi-1 轴，由 xi-1转向 xiA1A2A3A4A5A6O1O0第16页，共59页，2022年，5月20日，2点18分，星期四第17页，共59页，2022年，5月20日，2点18分，星期四第18页，共59页，2022年，5月20日，2点18分，星期四 对下图所示简单机器人，根据D-H法，建立必要坐标系及参数表。例 3第19页，共59页，2022年，5月20日，2点18分，星期四第一步：根据D-H法建立坐标系的规则建立坐标系 第20页，共59页，2022年，5月20日，2点18分，星期四第二步：将做好的坐标系简化为我们熟悉的线图形式第21页，共59

8、页，2022年，5月20日，2点18分，星期四第三步：根据建立好的坐标系，确定各参数，并写入D-H参数表#da1009020030040-90500906000第22页，共59页，2022年，5月20日，2点18分，星期四#da1009020030040-90500906000第23页，共59页，2022年，5月20日，2点18分，星期四第四步：将参数代入A矩阵，可得到第24页，共59页，2022年，5月20日，2点18分，星期四第5步 求出总变化矩阵第25页，共59页，2022年，5月20日，2点18分，星期四 依次写出从基坐标系到手爪坐标系之间相邻两坐标系的齐次变换矩阵，它们依次连乘的结果

9、就是末端执行器（手爪）在基坐标系中的空间描述，即已知q1,q2,qn，求 ，称为运动学正解；已知 ，求q1,q2,qn，称为运动学反解。上式称为运动方程。综上：第26页，共59页，2022年，5月20日，2点18分，星期四正解反解第27页，共59页，2022年，5月20日，2点18分，星期四1.5 机器人的逆运动学解 给定机器人终端位姿，求各关节变量，称求机器人运动学逆解。让我们通过下面这道例题来了解一下机器人逆运动学求解的一般步骤。前面例子最后方程为：求逆运动学方程的解第28页，共59页，2022年，5月20日，2点18分，星期四 根据第3行第4列元素对应相等可得到依次用 左乘上面两个矩阵，

10、得到：第29页，共59页，2022年，5月20日，2点18分，星期四根据1，4元素和2，4元素，可得到：将上面两个方程两边平方相加，并利用和差化积公式得到第30页，共59页，2022年，5月20日，2点18分，星期四已知于是可得到：依次类推，分别在方程2.19两边左乘A1A4的逆，可得到第31页，共59页，2022年，5月20日，2点18分，星期四接下来再一次利用式由于C12=C1C2-S1S2以及S12=S1C2+C1S2，最后得到：第32页，共59页，2022年，5月20日，2点18分，星期四最后用A5的逆左乘式2.67，再利用2,1元素和2,2元素，得到：123456关节变量都是第33页

11、，共59页，2022年，5月20日，2点18分，星期四2.10 机器人的运动学编程 在实际应用中，对运动学的求解是相当繁琐和耗时的，因此需要用计算机编程来实现。并且应尽量避免使用矩阵求逆或高斯消去法等相对繁琐的算法。正确的算法是：第34页，共59页，2022年，5月20日，2点18分，星期四2.11 设计项目 利用本书中所介绍的四自由度机器人，结合本章所学的知识进行四自由度机器人的正逆运动学分析。 SCARA型机器人的运动学模型的建立，包括机器人运动学方程的表示，以及运动学正解、逆解等，这些是研究机器人控制的重要基础，也是开放式机器人系统轨迹规划的重要基础。为了描述SCARA型机器人各连杆之间

12、的数学关系，采用D-H法。SCARA型机器人操作臂可以看作是一个开式运动链。它是由一系列连杆通过转动或移动关节串联而成的。为了研究操作臂各连杆之间的位移关系，可在每个连杆上固接一个坐标系，然后描述这些坐标系之间的关系。第35页，共59页，2022年，5月20日，2点18分，星期四SCARA（Selective Compliance Assembly Robot Arm装配机器人臂）机器人坐标系的建立 1.SCARA机器人坐标系建立原则根据D-H坐标系建立方法，SCARA机器人的每个关节坐标系的建立可参照以下的三原则(1) 轴沿着第n个关节的运动轴;基坐标系的选择为:当第一关节变量为零时，零坐标

13、系与一坐标系重合。(2) 轴垂直于 轴并指向离开 轴的方向。(3) 轴的方向按右手定则确定。 2.构件参数的确定根据D-H构件坐标系表示法，构件本身的结构参数 、 和相对位置参数 、 可由以下的方法确定:(1) 为绕 轴(按右手定则)由 轴到 轴的关节角。(2) 为沿 轴，将 轴平移至 轴的距离。(3) 为沿 轴从 量至 轴的距离。(4) 为绕 轴(按右手定则)由 轴到 轴的偏转角。第36页，共59页，2022年，5月20日，2点18分，星期四 3.变换矩阵的建立全部的连杆规定坐标系之后，就可以按照下列的顺序来建立相邻两连杆n-1和n之间的相对关系:(1)绕 轴转 角。(2)沿 轴移动 。(3

14、)绕 轴转 角。(4)沿 轴移动 。这种关系可由表示连杆n对连杆n-1相对位置齐次变换 来表征。即：展开上式得 第37页，共59页，2022年，5月20日，2点18分，星期四由于 描述第n个连杆相对于第n-1连杆的位姿，对于SCARA教学机器人(四个自由度)，机器人的末端装置即为连杆4的坐标系，它与基座的关系为: 第38页，共59页，2022年，5月20日，2点18分，星期四如上图坐标系，可写出连杆n相对于n-1变换矩阵 : 其中： 以下相同。相应连杆初始位置及参数列于表2.4，表中 、 为关节变量。 构件1000102001030010400010第39页，共59页，2022年，5月20日，

15、2点18分，星期四 各连杆变换矩阵相乘，可得到SCARA机器人末端执行器的位姿方程(正运动学方程)为下 式它表示了SCARA手臂变换矩阵 ，它描述了末端连杆坐标系4相对基坐标系0的位姿 。SCARA机器人的正运动学分析 第40页，共59页，2022年，5月20日，2点18分，星期四 SCARA机器人的逆运动学分析 1.求关节变量 为了分离变量，对方程的两边同时左乘 ,得:即：第41页，共59页，2022年，5月20日，2点18分，星期四 左右矩阵中的第一行第四个元素(1.4)，第二行第四个元素(2.4)分别相等。即: 由以上两式联立可得: 式中： 第42页，共59页，2022年，5月20日，2

16、点18分，星期四2 求关节变量 由式(2.87)可得: 式中：第43页，共59页，2022年，5月20日，2点18分，星期四3 求关节变量 再令左右矩阵中的第三行第四个元素(3.4)相等，可得: 4 求关节变量 再令左右矩阵中的第一行第一个元素、第二行第一个元素(1.1，2.1)分别相等，即：由上两式可求得: 第44页，共59页，2022年，5月20日，2点18分，星期四 至此，机器人的所有运动学逆解都已求出。在逆解的求解过程中只进行了一次矩阵逆乘，从而使计算过程大为简化，从 的表达式中可以看出它有两个解，所以SCARA机器人应该存在两组解。运动学分析提供了机器人运动规划和轨迹控制的理论基础。

17、第45页，共59页，2022年，5月20日，2点18分，星期四对机器人相关概念的补充退化：当机器人失去一个自由度，并因此不按所期望的状态运动时即称为退化。退化发生条件：1.机器人达到物理极限，不能进一步运动2.两个相似关节共线不灵巧区域：能对机器人定位不定姿的区域称为不灵巧区域。D-H法的局限性：无法表示关于y轴的运动。第46页，共59页，2022年，5月20日，2点18分，星期四退化状态下的机器人第47页，共59页，2022年，5月20日，2点18分，星期四总 结1 用矩阵表示点，向量，坐标系及变换的方法2 正逆运动学方程的建立3 用D-H法建立坐标系及变化方程4 正逆运动学方程的求解第48

18、页，共59页，2022年，5月20日，2点18分，星期四9.2 机器人杆件，关节和它们的参数 9.2.1 杆件与关节操作机由一串用转动或平移（棱柱形）关节连接的刚体（杆件）组成每一对关节杆件构成一个自由度，因此N个自由度的操作机就有N对关节杆件。0号杆件（一般不把它当作机器人的一部分）固联在机座上，通常在这里建立一个固定参考坐标系，最后一个杆件与工具相连关节和杆件均由底座向外顺序排列，每个杆件最多和另外两个杆件相联，不构成闭环。 关节杆件末端操作手机座两自由度关节第49页，共59页，2022年，5月20日，2点18分，星期四关节：一般说来，两个杆件间是用低副相联的只可能有6种低副关节：旋转（转

19、动）、棱柱（移动）、圆柱形、球形、螺旋和平面，其中只有旋转和棱柱形关节是串联机器人操作机常见的，各种低副形状如下图所示：旋转棱柱形柱形球形螺旋形平面第50页，共59页，2022年，5月20日，2点18分，星期四9.2.2 杆件参数的设定 条件关节串联每个杆件最多与2个杆件相连，如Ai与Ai-1和 Ai+1相连。第 i 关节的关节轴 Ai 位于2个杆件相连接处，如图所示， i -1关节和 i +1关节也各有一个关节轴 Ai-1 和 Ai+1。AiAi+1Ai-1第51页，共59页，2022年，5月20日，2点18分，星期四 杆件参数的定义 、 、 和 li 和 li-1 在 Ai 轴 线上的交点

20、之间 的距离。AiAi+1Ai-1 li 和 li-1 之间的夹 角，按右手定则 由li-1 转向 li。 由运动学的观点来看，杆件保持其两端关节间的形态不变，这种形态由两个参数决定：杆件长度 li 和杆件扭转角 。杆件的相对位置关系，由另外两个参数决定：杆件的距离 di 和杆件的回转角 。 li 关节 Ai 轴和 Ai+1 轴线公法线的长度。 关节i 轴线与i+1轴线在垂直于li 平面内的夹角。 第52页，共59页，2022年，5月20日，2点18分，星期四 上述4个参数，就确定了杆件的结构形态和相邻杆件相对位置关系。在转动关节中，li, i, di是固定值，i是变量。在移动关节中，li,

21、i , i是固定值， d i 是变量。第53页，共59页，2022年，5月20日，2点18分，星期四9.3 机器人关节坐标系的建立对于每个杆件都可以在关节轴处建立一个正规的笛卡儿坐标系（xi, yi, zi），（i=1, 2, , n），n是自由度数，再加上基座坐标系，一共有（n+1）个坐标系。基座坐标系 O0定义为0号坐标系（x0, y0, z0）,它也是机器人的惯性坐标系，0号坐标系在基座上的位置和方向可任选，但z0轴线必须与关节1的轴线重合，位置和方向可任选；最后一个坐标系（n关节），可以设在手的任意部位，但必须保证 zn与zn-1 垂直。第54页，共59页，2022年，5月20日，2点

22、18分，星期四机器人关节坐标系的建立主要是为了描述机器人各杆件和终端之间的相对运动，对建立运动方程和动力学研究是基础性的工作。为了描述机器人各杆件和终端之间转动或移动关系，Denavit和Hartenberg于1955年提出了一种为运动链中每个杆件建立附体坐标系的矩阵方法（D-H方法） ，建立原则如下： 9.3.1 D-H关节坐标系建立原则右手坐标系原点Oi：设在li与Ai+1轴线的交点上 Zi轴： 与Ai+1关节轴重合，指向任意 Xi轴： 与公法线Li重合，指向沿Li由Ai轴线指向Ai+1轴线 Yi轴： 按右手定则 第55页，共59页，2022年，5月20日，2点18分，星期四 9.3.2 关节坐标系的建立方法AiAi+1Ai-1原点Oi：设在li与Ai+1轴线的交点上 zi轴：与Ai+1关节轴重