六自由度运动平台的仿真研究

六自由度运动平台的仿真研究

天津工程机械研究院

杨永立

摘要：本文分析了六自由度运动平台分别采用球铰链和万向节铰链进行连接时的自由度，运用欧拉角、旋转变换的方法推导出位置反解方程，介绍了数值迭代法进行位置正解的过程。

关键词：并联，局部自由度，位置反解，位置正解。 1. 简介

运动平台按结构形式可分为串联和并联两大类。与串联形式相比，并联形式具有刚度大、承载能力强、结构简单、运动负荷小、能实现包括横移、纵移、升沉等多个自由度运动等特点。同时，串联形式的优点也很明显，其具有运动空间大，测量精度高，运动、受力分析相对简单、控制、测量的实现相对容易，且每个自由度都能独立运动等特点。

六自由度运动平台（如图1所示）是由六条油缸通过万向节铰链（或球铰链）将上、下两个平台连接而成，下平台固定在基础上，借助六条油缸的伸缩运动，完成上平台在三维空间六个自由度（X，Y，Z，α，β，γ）的运动，从而可以模拟出各种空间运动姿态。 2. 自由度的确定

若在三维空间有n个完全不受约束的物体，任选其中一个作为固定参照物，因每个物体相对参照物都有6个运动自由度，则n个物体相对参照物共有6(n-1)个运动自由度。若在所有物体之间用运动副联接起来组成机构，设第i个运动副的约束为ui（1到5之间的整数），如果运动副的总数为g，则机构的自由度M为：

1

M?6(n?1)??ui

i?1g利用上述公式计算一下如图1所示运动平台(采用球铰链)的自由度数。将油缸分解为缸筒和活塞杆，则总的构件数n=14，油缸与上下平台之间的连接为12个球铰链（约束为3），缸筒和活塞杆构成6个既可以相对移动，又可以相对转动的运动副（约束为4），则平台的自由度M为：

M?6(n?1)??ui=6 (14－1)－(3×12＋4×6)=18

i?1g计算结果出人意料，平台似乎无法只通过六条油缸进行驱动。但是，如果保持上平台和缸筒固定不动，由球铰链的特性可知，活塞杆仍然可以相对其轴线转动；同理，缸筒也具有同样的效应。实践证明，这种转动并不影响上平台的空间运动姿态，因此属于局部自由度。

在六自由度运动平台的实际设计中，由于球铰链的刚度差，结构不稳定，所以一般采用万向节铰链（如图2所示，约束为4）来代替图1中的球铰链，则自由度M为：

M?6(n?1)??ui=6 (14－1)－(4×12＋4×6)=6

i?1g

3. 六自由度运动平台空间姿态的解算

图2 万向节铰链

要实现对平台空间姿态的控制和测量，必须掌握它两个方向上的解算方法，即位置反解和位置正解。 3.1 位置反解（逆向解）：

已知输出件的位置和姿态，求解输入件的位置称为机构的位置反解。在运动平台的实际应用当中，用户所给定的一般都是平台的六个空间姿态参数X，Y，Z，α，β，γ，然而要实现对平台的控制，需要的是六条油缸的长度L1、L2?L6，这正好是已知输出求输入，属于位置反解。也就是说，要实现对平台空间姿态的控制，就必需推导出平台的位置反解方程。

如图1所示，在上平台建立动坐标系o-xyz，在下平台建立静坐标系O-XYZ，

2

那么，上平台的运动可分解为随o-xyz坐标原点o沿O-XYZ三个坐标轴方向上的平移（X、Y、Z），以及绕坐标轴的旋转（α，β，γ）。为了避免发生角度间的“耦合”，一般采用欧拉角来描述刚体的旋转状态，而欧拉角的定义又随旋转次序的不同而不同。本文将欧拉角定义为依次绕z轴旋转γ，绕y轴旋转β，绕x轴旋转α。下平台各铰点A1、A2、?A6的坐标（XA1,YA1,ZA1）、（XA2,YA2,ZA2）?（XA6,YA6,ZA6）和上平台各铰点a1、a2、?a6的动坐标（xa1,ya1,za1）、（xa2,ya2,za2）?（xa6,ya6,za6）为已知，只要求出对应姿态参数X，Y，Z，α，β，γ的上平台各铰点的静坐标（Xa1,Ya1,Za1）、（Xa2,Ya2,Za2）?（Xa6,Ya6,Za6），运用两点间距离公式便可以求出L1、L2?L6。

以A1和a1为例来计算与其相连的油缸的长度L1。在如上所述对运动进行分解的情况下，静坐标（Xa1,Ya1,Za1）和动坐标（xa1,ya1,za1）有如下变换公式：

?Xa1??xa1??X??Y???T???y???Y? ?a1??a1??????Za1???za1????Z??其中[T]是关于α、β、γ的旋转变换矩阵，公式中只有矩阵[T]未知，下面就来求该旋转变换矩阵。

a 绕z轴旋转γ b 绕y?轴旋转β c 绕x??轴旋转α

图3 欧拉角坐标系

根据本文定义的欧拉角，建立如图3所示的四个坐标系o?xyz、o??x?y?z?、

o???x??y??z??、o????x???y???z???。首先绕z轴旋转γ，其变换关系如下：

x=x?cosγ-y?sinγ

3

y=x?sinγ+y?cosγ z=z? 写成矩阵形式为：

?x??cos??y???sin??????z????0?sin?cos?00??x???x???y??=[C]??y?? 0????????1????z????z???绕y?轴旋转β，其变换关系为：

?x???cos??y????0?????z??????sin?0sin???x????x????y???=[B]??y??? 10????????0cos?????z?????z????绕x??轴旋转α，其变换关系为：

0?x????1?y?????0cos??????z??????0sin???x?????x?????y????=[A]??y???? ?sin?????????cos?????z??????z?????0将上述三个变换公式合并，得：

?x??x?????y???C???B???A???y???? ???????z???z?????由上述变换过程可知，在只有旋转，没有平移的情况下，（Xa1,Ya1, Za1）等同于(x, y, z)，同为静坐标，而（xa1, ya1, za1）等同于(x???,y???,z???)，同为动坐标，所以，旋转变换矩阵[T]= [C]×[B]×[A]=

?cos?cos??cos?sin?????sin??cos?sin??sin?sin?cos?cos?cos??sin?sin?sin?sin?cos?sin?sin??cos?sin?cos???sin?cos??cos?sin?sin???

?cos?cos??将矩阵[T]代入前面的变换公式，即可求出a1的静坐标（Xa1, Ya1, Za1），从而可以求出油缸的长度L1为：

L1?(Xa1?XA1)2?(Ya1?YA1)2?(Za1?ZA1)2

同理可得：

L2?(Xa2?XA2)2?(Ya2?YA2)2?(Za2?ZA2)2

4

? ? ? ? ?

L6?(Xa6?XA6)2?(Ya6?YA6)2?(Za6?ZA6)2

3.2 位置正解（顺向解）：

已知机构输入件的位置，求解机构输出件的位置和姿态称为机构的位置正解。由于六自由度运动平台是并联机构，直接测量平台的六个自由度的空间姿态相当困难，但可以通过位移传感器测出每条油缸的长度，再经过位置正解间接求出平台的空间姿态。到目前为止，还没有直接的正解方程式，只能采用数值迭代的方法，利用计算机快速运算的特点来逼近求解平台姿态。目前所提出的迭代求解的方法很多，本文所采用的是牛顿法，其基本原理就是将非线性方程组变成线性方程组，求出近似解，然后在此近似解基础上进一步迭代，逐步逼近非线性方程组真解。由位置反解方程组可得：

f1(X，Y，Z，α，β，γ)=L12 – [( Xa1–XA1)2+( Ya1–YA1)2+( Za1–ZA1)2]=0 f2(X，Y，Z，α，β，γ)=L22 – [( Xa2–XA2)2+( Ya2–YA2)2+( Za2–ZA2)2]=0 ? ? ? ? ? ? ? ? f6(X，Y，Z，α，β，γ)=L62 – [( Xa6–XA6)2+( Ya6–YA6)2+( Za6–ZA6)2]=0 此时，六条油缸长度L1、L2 ?L6已知，需要求解平台的空间姿态X，Y，Z，α，β，γ。设X=(X，Y，Z，α，β，γ)＇，F(X)=( f1()，f2()，f3()，f4()，f5()，f6())＇，则位置正解的迭带公式为：

?(Xk) –1 F (Xk) Xk+1= Xk –F移项，得到六个方程六个未知数的线性方程组：

?(Xk)( Xk+1 – Xk)+ F (Xk)=0 F求解方程组便可以得到Xk+1。选择适当的初始点X0（如：X0= (0,0,0,0,0,0)）和终止条件（如：|Xk+1- Xk|

本文运用SolidWorks软件建立了六自由度运动平台的三维模型（如图4所示），并构造了与图3对应的欧拉角坐标系（如图5所示），分别输入油缸的长度和平台的空间姿态参数，经验证，本文所述正解和反解方法的计算结果均符合要求。

5

图4 六自由度运动平台 图5 欧拉角坐标系

参考文献：

[1]北京亿美博科技有限公司：杨世祥 杨 涛，北京大学：徐悦桐，《大型数字式六自由度运动平台的开发》，发表在《工业数字化先锋》。

[2]秦皇岛燕山大学：黄真、孔令富、方跃法，《并联机器人机构学理论及控制》，北京：机械工业出版社，1997。

[3]林成森，《数值计算方法》，北京：科学出版社，2001。

6

图4 六自由度运动平台 图5 欧拉角坐标系

参考文献：

[1]北京亿美博科技有限公司：杨世祥 杨 涛，北京大学：徐悦桐，《大型数字式六自由度运动平台的开发》，发表在《工业数字化先锋》。

[2]秦皇岛燕山大学：黄真、孔令富、方跃法，《并联机器人机构学理论及控制》，北京：机械工业出版社，1997。

[3]林成森，《数值计算方法》，北京：科学出版社，2001。

6

本文来源：https://www.bwwdw.com/article/jn7p.html