Page 24 - 《中国医疗器械杂志》2026年第2期
P. 24

Chinese Journal of Medical Instrumentation                                         2026年 第50卷 第2期

                                                     医   疗   机   器   人



                                                                               T
              杆扭转角4个参数确定各个连杆之间的运动学关                             式中:[p , p , p ] 代表末端工具坐标系相对于原点
                                                                              z
                                                                           y
                                                                        x
              系。如图4所示建立六自由度机械臂连杆坐标系,                            的位置向量;[n, o, a]代表末端工具坐标系相对于
              确定的DH参数如表1所示。                                     世界坐标系的旋转矩阵。

                                                                                                         [8]
                                                     X
                                    a 2                             对机械臂进行控制时,通常对逆运动学 更感
                                                     Y          兴趣,即已知机械臂的连杆参数,给定末端工具坐
                                                     Z
                  d 1   1                                       标系的期望位置和姿态,找到相应的关节变量,使
                                               d 4
                        0            2               4  d 5     末端工具坐标系能够到达期望位姿。
                                                                    本文所研究的机械臂后3个关节轴交于一点,
                                               3      6
                                        a 3
                                                    5           满足Pieper准则,故逆运动学解析解存在。逆运动
                                                                学求解思路如下:首先,将关节变量逆向解耦为前
                                                       d 6
                      图4   机械臂3D模型和连杆坐标系分布                      3个关节变量和后3个关节变量;然后,利用末端位
               Fig.4  3D model of the robotic arm and distribution of link coordinate
                                   frames                       置求出前3个关节变量的值,利用末端朝向和已知


                             表1   机械臂DH参数                       关节角求出后3个关节变量的值,即可得到机械臂
                        Tab.1  DH parameters of the robotic arm
                                                                的8组完整解析解,其逆解流程如图5所示。
                 关节i      θ i /rad  d i /m   a i /m   α i /rad
                                                                  输入末端     求解腕心     求解关节      求解关节     输出可能
                  1        θ 1     0.152      0        π/2        位姿 T 6    位置       1~3位置    4~6姿态     关节角
                                                                     0
                  2        θ 2      0       −0.425     0
                                                                               图5   机械臂逆解流程
                  3        θ 3      0       −0.395     0
                                                                       Fig.5  Inverse solution flowchart of the robotic arm

                  4        θ 4     0.102      0        π/2
                                                                    首先计算前3个关节的腕心位置,腕心是后3个
                  5        θ 5     0.102      0       −π/2
                                                                关节轴线的交点,根据正运动学可得腕心位置:
                  6        θ 6     0.100      0        0
                                                                                                
              注:关节角θ i 、连杆长度a i 、连杆偏距d i 、连杆扭转角α i 。                     0      wc x       p x       a x    
                                                                               
                                                                                       
                                                                                                 
                                                                                                 
                                                                               
                                                                                       
                                                                               
                                                                                                 
                                                                                       
                                                                                                  (5)
                                                                               
                                                                          P wc =  wc y  =  p y  −d 6  a y 
                                                                                       
                                                                                                
                  根据上文内容建立关节1~6的位姿变换关系                                           wc z    p z      a z
                                                                       0
              0                                                 式中: P 为腕心位置;d =0.100 m为关节6偏移
               T 。方法是通过矩阵变换相邻关节的DH变换矩阵                                   wc              6
                6
                                                                             T
              链乘:                                               量;[a , a , a ] 是末端的z轴方向,来自旋转矩阵第
                                                                         y
                                                                            z
                                                                      x
                     0    0    1    2   3    4    5             3列。
                      T 6 = T 1 × T 2 × T 3 × T 4 × T 5 × T 6  (1)
                  由机械臂DH参数可以得到连杆坐标系{i}相对                            如图6所示,从基坐标系俯视xy平面,用解耦
              于坐标系{i−1}的变换矩阵:                                   法求解关节1的θ :
                                                                               1
                                                                  θ 1 = arctan2(wc y ,wc x )或θ 1 = arctan2(wc y ,wc x )+π (6)
                                                       
                      cosθ i −sinθ i cosα i sinθ i sinα i a i cosθ i 
                                                       
                                                       
                                                       
               i−1     sinθ i cosθ i cosα i −cosθ i sinα i a i sinθ i    (2)  此时产生2个候选解左侧构型和右侧构型,需
                     
                                                        
                 T i = 
                                                       
                       0     sinα i     cosα i     d i  
                                                              根据实际情况选择。
                                                       
                                                       
                        0       0          0        1
                                                                          y                    腕心
                  由变换矩阵定义可知,连杆坐标系{0}相对于
                                                                         wc y
              世界坐标系{B}的变换矩阵和工具坐标系{7}相对于
              连杆坐标系{6}的变换矩阵分别为:                                                 关节2
                                         1 0 0    0  
                                          
                     0 0 1 −d 0                       
                                                     
                                                   
               B     0 −1 0   0      6    0 0 −1 −d 7    (3)
                    
                                          
                                                        
                    
                T 0 =           , T 7 =            
                     1 0 0    0         0 1 0    0   
                                  
                    
                                          
                                                     
                                                     
                      0 0 0    1            0 0 0    1
                  通过变换矩阵链乘可以得到末端工具坐标系                                            θ 1
              {7}相对于世界坐标系{B}的变换矩阵:                                       O    基座                wc x  x
                                                                           图6   腕心几何关系俯视图
                                       n x  o x  a x  p x 
                      [            ]                              Fig.6  Top view of wrist center geometric relationships
                        n  o  a  p                   
                 B                     n y  o y  a y  p y  
                                      
                  T 7 =             =                  (4)
                        0  0  0  1       n z  o z  a z  p z     如图7所示,在竖直平面内建立二连杆模型,
                                                     
                                                     
                                        0   0   0   1           求解关节2、3,计算腕心到基座的距离:
                                                             138
   19   20   21   22   23   24   25   26   27   28   29