并联机构的模型图如图1所示。由图1可知,并联机构共有3条支链,其中底边两条支链为相同的UPS支链,顶边为RPR支链,RPR支链中两个R副空间位置相互垂直。为了方便分析该机构的正反解,规定定平台左端顶点为A
,逆时针分布得到A
、A
;动平台左端顶点为B
,逆时针分布得到B
、B
。其中,A
、A
、A
分别为U
、U
、R
的质心点;B
、B
、B
分别为S
、S
、R
的质心点。以定平台底边中心为定坐标系原点建立坐标系o
-xyz
,x
轴与定平台底边重合,x
轴方向指向第二条支链底端U
的质心处,z
轴方向与定平台方向垂直,y
轴方向由右手螺旋定则确定;以动平台底边中心为动坐标系原点建立坐标系o
-x
y
z
,x
轴与动平台底边重合,方向指向第二条支链顶端S
的质心处,z
轴方向与动平台方向垂直,y
轴与动平台高重合,方向指向动平台顶点,x
轴根据右手螺旋定则确定。2-UPS/RPR机构简图如图2所示。![](https://img.fx361.cc/images/2022/1111/e1392993a128a2d52c9e04d2c10773724bd38f9f.webp)
图1 2-UPS/RPR机构模型图2 2-UPS/RPR机构简图
1.2 自由度运算
给定螺旋$
=(L
,M
,N
;P
,Q
,R
)与螺旋$
=(L
,M
,N
;P
,Q
,R
),则两螺旋的互易积表示为:$
○$
=L
P
+M
Q
+N
R
+L
P
+M
Q
+N
R
,(1)
若两螺旋的互易积为0,则两螺旋互为反螺旋。在并联机构中,各分支运动螺旋数目与其对应的约束螺旋数目之和为6。对于2-UPS/RPR并联机构,先求得3条支链的运动螺旋系,然后根据互易积理论可得到与其对应的约束螺旋系,接着将3条支链的约束螺旋系合并后得到动平台的约束螺旋系,最后对其求反螺旋系,从而得到动平台的运动螺旋系。
第一条UPS支链共有6个运动螺旋,组成的运动螺旋系可表示为:
![](https://img.fx361.cc/images/2022/1111/af25f2658bc76957026b69c4f1fe271b4d8ee37d.webp)
(2)
式中,l
、m
、n
中的ij
表示第i
条支链的第j
个运动副;l
表示该运动副位置矢量的方向余弦;X
1、Y
1、Z
1为$
的位置矢量。从式(2)中还可看出,UPS支链的6个运动螺旋线性无关,因此该运动螺旋系不存在与之对应的约束螺旋系,即该支链对动平台无约束力与约束力偶。第二条支链的构型与第一条支链构型相同,且相对于机构几何中心点呈对称分布,因此第二条支链也无约束螺旋系。第三条RPR支链在定坐标系中的运动螺旋为:
![](https://img.fx361.cc/images/2022/1111/4865f344200f02c7f0aaaaaa0a4682f2fdeb6534.webp)
(3)
根据互易积公式,RPR支链的约束螺旋系有3个约束螺旋:
![](https://img.fx361.cc/images/2022/1111/86f92b9cd38155bb1c78dc4d0aa673bd182e80ed.webp)
(4)
根据3条支链的约束螺旋系得到动平台的运动螺旋系为:
![](https://img.fx361.cc/images/2022/1111/eb94d3b936e62714d3bc3677b14c894a143d7ad8.webp)
(5)
![](https://img.fx361.cc/images/2022/1111/1784433cca13082ef680042b30ac210de30f3ce2.webp)
如果机构中运动螺旋数目超过6,则超出部分为并联冗余约束,所以过去常用的Grübler-Kutzbach(以下简称G-K)公式无法对所有机构求得正确的结果,因此对G-K公式加以修正,修正后G-K公式如式(6)所示。
![](https://img.fx361.cc/images/2022/1111/57b9f04b5a3a8f2443241a309aa3edf1ba01680b.webp)
(6)
式中,dof
为机构自由度;m
为刚体自由度;N
为构件数量(定平台也看作为一个构件);J
为关节的数目;f
为第i
个关节的自由度数;ζ
为机构中全部过约束的总数;υ
表示并联冗余约束;本机构属于单环机构,因此无并联冗余。对于2-UPS/RPR并联机构,m
与λ
之和为6,公共约束λ
=0,即m
=6-λ
=6,构件数量N
=11,关节数目J
=12。其中转动副的自由度为1,移动副的自由度为1。虎克副的自由度为2,冗余自由度υ
=0,过约束自由度ζ
=0。将其代入修正的G-K公式可得:dof
=6×(11-12-1)+(6×2+3)+0-0=3,由G-K公式计算得出2-UPS/RPR并联机构的自由度为3,与螺旋理论计算出的结果一致。
2 机构的运动分析
2.1 位置反解
在动平台位姿确定后求其他构件运动状态为位置反解。将动平台中心点的位置设为P
,因机构只能绕x
、y
轴转动,沿z
轴移动,所以3个姿态角中绕z
轴转动的角度γ
为0,设绕x
轴转动的角度为α
,绕y
轴转动的角度为β
,根据以上参数可求得支链上的3个驱动副的位移距离。![](https://img.fx361.cc/images/2022/1111/7355617ae07c5ad417c5416918d6d0330e90ecb3.webp)
![](https://img.fx361.cc/images/2022/1111/1841820fa5f72036ad71f9a31ff38f1bc75e3964.webp)
(7)
将动坐标系原点O
表示为定坐标系中的坐标矢量P
:P
=(x
,y
,z
),(8)
![](https://img.fx361.cc/images/2022/1111/9fd9d8947c6daae4cdcaaec2c03b75e3c416ae92.webp)
![](https://img.fx361.cc/images/2022/1111/b2bfa72ea447b37b8b45be9c47640563a0d2e1ae.webp)
(9)
式中,s
表示sin
;c
表示cos
。动平台在固定坐标系中的坐标矢量为OB
,将OA
与OB
的矢量差长度设为l
,可得:![](https://img.fx361.cc/images/2022/1111/4afa1f2d8c836e24b4a81ad3a44dc894456ddff4.webp)
(10)
l
=|OB
-OA
|,(11)
代入计算得:
![](https://img.fx361.cc/images/2022/1111/ea07a21099f4db928efe9b7f1df6cdfe73c8d280.webp)
(12)
2.2 位置正解
位置正解即根据驱动值求解动平台中心点位置参数,该方法是位置反解的一种逆运用。本例为一般构型,采用方法为数值分析法,将位置反解方程式整理得到动平台坐标系原点位置的求解方程:
![](https://img.fx361.cc/images/2022/1111/6460d1eae50f908e0041699f1691962d12a07cf0.webp)
(13)
该方程为多元非齐次线性方程,常规求根公式无法求得该方程的解析解,但该方程在单根附近平方收敛,因此采用牛顿-拉夫逊迭代法对函数进行更新迭代,从而求得近似解,牛顿迭代公式如下:
![](https://img.fx361.cc/images/2022/1111/bbd8bfe9a06b5fadfa0486ad0f7f2f0a8bc61c35.webp)
(14)
已知3个移动副伸缩量l
、l
、l
,初始向量T
,经过不断更新迭代后得到最终向量T
,同时计算得到3条支链移动副移动的距离,该方法可求得动平台最终的位姿变化。3 速度雅可比矩阵分析
选取并联机构3条支链的移动副作为驱动向动平台输入速度,输出速度由雅可比矩阵来映射到动平台上。将位置逆解分别对α
、β
、z
求一阶导数,求得的参数代入雅可比矩阵中,雅可比矩阵J
为:![](https://img.fx361.cc/images/2022/1111/e197c6ea22121ca5dbaff9692cd8805ea5e565e8.webp)
(15)
其中,
![](https://img.fx361.cc/images/2022/1111/1633beed2806e97dde058f58acda63b707a9d906.webp)
(16)
对位置逆解方程两侧同时求导并加以整理得到:
![](https://img.fx361.cc/images/2022/1111/e4eef1b35bb6dd57cecdca8ca644bba321ee9492.webp)
(17)
![](https://img.fx361.cc/images/2022/1111/d0d45be9756900905162c16a6d9f5728c621d596.webp)
4 运动学仿真
为了分析2-UPS/RPR机构在输入3个驱动函数下动平台变化的规律,在3条支链上的3个移动副设置驱动函数进行驱动,驱动函数如下:
![](https://img.fx361.cc/images/2022/1111/67a1849f10c63df8bc98ca38c4fb4c501681a7f2.webp)
(18)
设定仿真时间为30 s,步数为500,机构开始运动后得到动平台运动仿真曲线云图。位移、角速度、速度、加速度变化曲线图分别如图3、图4、图5、图6所示。由图3、图4可知,Adams仿真云图呈周期性变化,周期为12.5 s,动平台在x
轴方向位移变化较小,角速度变化较大;在y
轴方向位移变化较大,角速度变化较小。由图5、图6可知,动平台在z
轴方向速度与加速度变化高于x
轴、y
轴。因为z
轴为移动,另外两轴为转动,而移动副行程较长,从位置反解中也能体现出角度变化率比杆长变化率小。总体来看,动平台参数变化曲线光滑连续,中间无断点与突变,表明机构能平稳运行。综上,2-UPS/RPR并联机构的运动性能良好。![](https://img.fx361.cc/images/2022/1111/d25864cbedbf8dee44307e69fcda801d157f2271.webp)
图3 位移变化曲线图
![](https://img.fx361.cc/images/2022/1111/ad4cdb7571d4d5f357d999fa0ef85cc71a848930.webp)
图4 角速度变化曲线图
![](https://img.fx361.cc/images/2022/1111/c2c0604c349e2dd7535945fff8843a5d46233324.webp)
图5 速度变化曲线图
![](https://img.fx361.cc/images/2022/1111/a9ac7923bd0d255616fc2c1ebf3ed9152b2004ab.webp)
图6 加速度变化曲线图
5 结论
基于螺旋理论求得2-UPS/RPR并联机构存在绕x
轴、y
轴转动与沿z
轴移动自由度,表明该机构可调节不同角度对不规则曲面产品进行喷涂。运用封闭矢量法与坐标转换法求得并联机构反解方程与正解方程,对其运用微分法求得雅可比矩阵,从而得到其运动学特性。利用Adams动态仿真得出机构的运动变化曲线图,从曲线良好的运动学性能可看出该机构在进行喷涂工作时平稳流畅,具有一定的实用价值。