张国英, 刘冠峰, 管贻生

(1.广东工业大学 机电工程学院, 广州 510006； 2.广东技术师范大学 机电学院, 广州 510665)

两转一移(2R1T)并联机构因其自由度少、结构简单、制造成本低、适用于多种场合等优势，已引起国内外学者的广泛关注. 自Hunt[1]于1983年提出一种构型为3-RPS并联机构以来，2R1T并联机构的构型综合近年来在国内也掀起一股研究热潮[2-9]. 这类机构目前已成功应用于五轴混联机床的转动中心[10]、姿态调整装置[11]以及坐标测量机[12]等.

速度分析是并联机构运动性能分析和动力学分析的基础，属于微分运动的范畴，常用的分析方法有闭环矢量法、旋量法、影响系数法[13-14]和网络分析法[15]等. 前三种方法实质是求解速度映射即雅克比矩阵，一般通过对位置约束方程求导来获得速度映射关系，但在位置约束方程本身就很复杂的情况下，对其进行一阶求导就更加困难. 网络分析法以图论为基础，适合多环强耦合并联机构的速度和力分析，对于少自由度并联机构的速度分析通用性不强.

为避开雅克比矩阵和位置导数的求解，本文在文献[9]的基础上，利用Riemann对称空间理论对一种自由度类球面并联机构进行了瞬时速度分析. 从机构子链旋量系的对称性出发，建立各子链关节运动的约束，从而简化机构动平台瞬时速度的计算. 根据计算和仿真的对比结果来看，角速度误差控制在-0.004～0.006 rad/s以内，线速度误差控制在-0.01～0.015 mm/s以内. 采用对称空间理论方法计算瞬时速度用时5.187 s.

1 机构描述及其坐标系建立

机构由等半径的动平台和基座以及3个相同的支链组成. 在初始位形下，机构的几何简图如图1所示. 3个支链呈120°均匀分布于基座和动平台之间. 与基座相连的转动副记为Bi(i=1,2,3)，它们的轴线相交于基座中心OB，与动平台相连的转动副记为Di(i=1,2,3)，它们的轴线相交于动平台中心OD；每条支链上的球副记为Si(i=1,2,3)，其空间布置为SiOB⊥OBBi,SiOD⊥ODDi,‖SiOD‖=‖SiOB‖. 任意位形下，动平台和基座关于3个球副所组成的平面S1S2S3对称，动平台可以绕对称面内的任意直线发生连续转动，也可以沿两平台中心连线方向进行连续平移.ODS1S2S3OB为一个六面体，平面S1S2S3将此六面体分为上下对称的两部分，OB和OD始终关于S1,S2,S3组成的平面对称.

图1 并联机构的几何简图

建立基坐标系OBxyz(记为{FW})和动坐标系ODuvw(记为{FT})，如图2所示. 初始位形下，基座与动平台平行，各转动副轴线也相互平行.OBOD⊥S1S2S3交S1S2S3于点M，‖OBM‖=‖ODM‖,‖OBSi‖=‖ODSi‖=l.

图2 机构坐标系

2 位置正解

在基坐标系{FW}中，S1,S2,S3的位置由驱动副Bi的转角αia(i=1,2,3)和l决定，可写为

Si=Ri[-lcαia,0,lsαia]T,

(1)

其中

对称平面方程ax+by+cz+d=0由各球副的位置坐标确定，满足下式:

(2)

根据点关于平面对称点的求解原理，容易得到动平台中心OD在{FW}下的位置坐标：

(3)

其中，A=a2+b2+c2. 展开式(3)即可求得机构的位置正解为

(4)

其中

si表示sinαia，ci表示cosαia，s12表示sin(α1a+α2a)，s13表示sin(α1a+α3a).

3 瞬时速度分析

3.1 Riemann对称空间理论

定义1[16]设p是Riemann流形M中的一点，如果存在M的一个变换σp：M→M满足下面3个条件：

1)σp是对合的，即σp2=id(id表示恒等变换)但σp≠id；

2)σp是保长变换(即等距变换)；

3)p是σp的孤立不动点，即有p的邻域U，使得∀q∈U,σp(q)=pq-1p.

则称M对于p成中心对称，σp叫作关于p的中心对称.

定义2如果一个Riemann流形M对于M中任意一点p，q都成中心对称，则称M为Riemann对称空间(简称对称空间).

Riemann对称空间的2个典型例子就是欧几里德空间和球体[17]. 在机构学中存在一种符合Riemann对称的并联机构，它是一类具有对称几何结构支链的并联机构. 每条支链可分为相连的两个子支链(C+,C-)，如图3所示. 它们的运动副的空间分布满足一种特殊的对称，即可以通过符号相反的一对运动{φ,-φ}使得两个子链重合为一个子链.

图3 对称并联机构的子支链图

如果重合子链的旋量系为ξk,…,ξ1，则原有支链的2个子支链的旋量系为

容易得到子链{C++C-}的运动学正解为

f(α)=eξ1(α1+φ1)…eξk(αk+φk)eξk(αk+φk)…eξ1(α1+φ1)

如果在刚体运动群SE(3)上定义一个反射映射：σp(q)=pq-1p，则上述子链的正向运动学可以写为

f(α)=σeξ1(α1+φ1)(σeξ2(α2+φ2))(…σeξk(αk+φk)(e)…).

考虑与f像空间对应的SE(3)的子流形Q. 由于f是多个反射的组合，Q必须是SE(3)的对称子空间.

定理1如果TeQ=ξk,…,ξ1⊂SE(3)，则f的像空间包含Q的关于e的一个邻域.

定理2如果并联机构的所有支链的约束空间的和正好是TeQ的正交补，则并联机构的任务空间包含Q的关于e的一个邻域. 并且每条子链的C+，C-转过的角度正好相反.

本机构任意时刻的运动都可以想象为在一个球体(或椭球体)内进行. 动平台与基座始终对称，对称面为过S1,S2,S3的圆面(或椭圆面)，如图4(a)所示. 由于球副S可以等效3个转动副R串联，各等效转动副的交点为球副的中心，故将本机构中各支链的球副Si(i=1,2,3)局部等效为3个轴线相交于球副中心的转动副R2-,R2+,R3+(R3-)，如图4(b)所示，R2-的轴线沿SiOB方向，R2+的轴线沿SiOD方向，R3-(R3+)的轴线位于对称面S1S2S3内且与连接两平台的转动副轴线同向. 这样，机构的每条支链上下两部分可以相互折叠后重合，成为在对称平面里的同一个RRR结构.

(a) (b)

至此可以得出结论：本机构是一个具有Riemann对称空间性质的机构，其每条支链经过等效处理后，支链的上下2条子链也是对称的，因此，上子链与动平台的夹角和下子链与基座的夹角相互为负. 即转动副Bi,Di的转角关系为：αip=-αia(i=1,2,3)，其中αia为平面OBBiSi与基座的夹角，αip为平面ODDiSi与动平台的夹角. 以B3,D3的转角关系α3p=-α3a为例，如图5所示.

图5 转动副B3, D3的转角

3.2 等效瞬时速度空间

构造运动旋量ξi(i=1,2,3)与每个关节对应，它表示第i个关节的旋量运动，此时除第i个关节外其他所有关节均固定于αja=0的位置. 在研究转动关节时，其运动旋量ξi可表示为

(5)

其中，wi∈R3表示旋转轴线上单位矢量，qi∈R3表示轴线上任一点. 对移动关节有

(6)

下面分别以机构的3种特殊运动为例，推导其运动旋量表达式.

图6 第一种运动：固定S1, S2，让S3自由运动

Fig.6 Motion pattern Ⅰ: fixS1,S2and freeS3

该运动情形下的运动旋量为

(7)

(8)

(9)

运动旋量系{ξ1,ξ2,ξ3}等效为绕对称面S1S2S3内任意一条直线的转动(2个旋转自由度)，以及沿对称面法线方向的平移(1个平移自由度)，其维数等于机构的自由度数，很好地解释了类球面并联机构2R1T的自由度性质，物理意义明确.

由于动平台的运动是由3个驱动关节共同作用引起，其运动可以由旋距为零的速度旋量来描述:

(10)

(11)

式中，ξ1,ξ2,ξ3分别为驱动关节B1,B2,B3在基坐标系下的等效运动旋量，ξ的集合就构成了机构末端在基坐标系下的速度空间.

4 计算与仿真

选取机构主要参数为：l=100 mm. 不失一般性，令驱动关节的角位移分别为

(12)

根据式(1)～(11)对动平台的瞬时速度(包括角速度和线速度)进行ADAMS仿真并与Mathematica计算结果进行对比，结果如图7所示. 由图7可知，仿真结果与理论计算几乎完全吻合，说明采用对称空间理论分析机构的瞬时速度是正确、可行的.

(a)角速度的仿真值与计算值对比

(b)线速度的仿真值与计算值对比

Fig.7 Instantaneous velocity computed and simulated by symmetric space theory

5 结 论

对一种对称结构的3自由度类球面并联机构进行了位置分析. 由于位置解析式较为复杂，若通过对其进行一阶求导或推导其雅克比矩阵的方式来获得机构的瞬时速度将是十分困难的事情. 为简化计算，利用对称空间理论，从机构子链旋量系的对称性出发，建立各子链关节运动的约束. 由于机构有三个驱动关节，每次锁定2个驱动关节，只让1个关节运动，这样就可以得到3种特殊的运动及其运动旋量，这些运动旋量和对应的驱动关节速度的线性组合构成了机构的瞬时速度空间. 这个空间的基正好解释了机构的自由度性质为两旋转一平移，物理意义明确，计算方法简单.