文章快速检索  
  高级检索
捕获轨迹系统并联机构地面标定方法
谢峰1,2, 洪冠新1, 张晨凯2, 魏忠武2, 马汉东2     
1. 北京航空航天大学 航空科学与工程学院, 北京 100083;
2. 中国航天空气动力技术研究院, 北京 100074
摘要: 捕获轨迹系统(CTS)是一种先进的预测外挂物投放轨迹的试验系统,普遍采用六自由度(6-DoF)串联机构作为其运动机构,串联机构因惯性力大和关节累积误差大使其定位精准度不足。相比串联机构,并联机构具有惯性力小和关节误差不累积等优点。采用6-PTRT并联机构作为CTS试验系统的六自由度运动机构,在空间受限的风洞环境中对CTS并联机构进行地面标定:提出动平台位姿的测量和计算方法,建立包含直线驱动平台安装夹角修正的标定模型,并基于非线性最小二乘法辨识结构参数。辨识后CTS并联机构的位移定位准度优于0.1 mm,姿态定位准度优于0.05°,最后以CTS并联机构和常规攻角机构进行8#标模的对比风洞试验。风洞试验结果表明,CTS并联机构的风载定位准度满足测力试验精准度要求。
关键词: 捕获轨迹系统    六自由度并联机构    地面标定方法    定位准度    风洞试验    
Ground calibration method of captive trajectory system parallel mechanism
XIE Feng1,2, HONG Guanxin1, ZHANG Chenkai2, WEI Zhongwu2, MA Handong2     
1. School of Aeronautic Science and Engineering, Beihang University, Beijing 100083, China;
2. China Academy of Aerospace Aerodynamics, Beijing 100074, China
Abstract: Captive Trajectory System (CTS) is an advanced test system for predicting the release trajectory of the external storages, which commonly uses the six-degree-of-freedom (6-DOF) series mechanism as the motion mechanism. As the large inertia force and the accumulated joint error, the positioning precision of the series mechanism is not enough. Compared with the series mechanism, the parallel mechanism has the advantages of small inertia force and non-accumulation of joint errors. In this paper, 6-PTRT parallel mechanism is employed as the 6-DOF motion mechanism for the CTS. The ground calibration method of the CTS parallel mechanism is studied in the space-constrained wind tunnel environment:the method of measuring and calculating the position and posture of the moving platform is presented, the calibration model which included the clamped angel repairman of the linear drive platform is established, the structural parameters are identified base on the method of nonlinear least squares. After identification, the position accuracy of the CTS parallel mechanism is better than 0.1 mm and the posture accuracy is better than 0.05°. Finally the CTS parallel mechanism and the regular attack angle mechanism are compared in the wind tunnel test. The wind tunnel test results show that the position and posture accuracy of the CTS mechanism under wind-load meets the requirements of the force wind tunnel test.
Keywords: captive trajectory system    6-DOF parallel mechanism    ground calibration method    position and posture accuracy    wind tunnel test    

捕获轨迹系统[1-4](Captive Trajectory System, CTS)是一种先进的预测外挂物投放轨迹的试验系统。CTS试验系统包含一套六自由度运动机构,用来控制外挂物模型在风洞中进行线位移运动和角位移运动。目前CTS试验系统普遍采用六自由度串联机构作为其运动机构,六自由度串联机构使用独立的电机和装置控制每个自由度,各自由度串联组合实现外挂物模型的六自由度运动[5-6]。六自由度串联机构运动关节逐级叠加,位移和角度误差逐渐累积,使得外挂物模型定位准度不足。

在精确定位的工业机器人、医用机器人、数控机床等领域,并联机构得到广泛应用。并联机构是一种以若干并联杆件连接静平台和动平台的空间机构。并联机构运动部件的质量明显减小,且各支链并联连接支撑起动平台,各支链结构误差不会对动平台位姿造成累积影响,因而相比串联机构而言并联机构具有惯性力小、关节误差不累积等优点[7]。常见的六自由度并联机构的结构形式有:6-SPS、6-PSS和6-PTRT等。

关于并联机构的定位准度,零件的加工与装配误差、构件的变形、运动副间隙等误差可占到全部定位误差的60%~70%[7]。在不改变现有并联机构任何结构和硬件尺寸条件下,开展运动学标定是提高并联机构定位准度的一种有效手段。

机构学研究者针对并联机构的标定方法开展了研究。樊锐等[8]基于拆分标定方法对6-UPS并联机构进行了标定。曾磊[9]针对3-RRR并联机构,提出了考虑减速机背隙误差的运动学标定误差模型。戴智武等[10]针对Delta并联机器人的误差模型,采用了粒子群算法辨识误差参数。李寅翔[11]分析了2-UPR-RPU并联机构的几何误差源,以其零点误差模型进行标定。马志强等[12]基于3-PPPS并联机构的位姿误差与几何误差之间的映射关系建立了运动学逆解方程,对其误差参数进行了标定。张立杰等[13]基于D-H参数建立了球面5R并联机构的运动学反解模型,并对其误差参数进行了标定。皮阳军等[14]提出了一种基于关节力传感器的并联机构的运动学标定方法,通过关节驱动力构造辨识模型,实现了结构参数的标定。张建中[15]通过构造Stewart平台姿态的雅克比矩阵,利用倾角仪测量并用遗传算法对参数进行辨识,完成了对Stewart平台的标定工作。孙皓[16]针对一种新型5-UPS/PRPUU并联机床,采用逐次固定每根杆的方法,标定每个支链铰点坐标和杆长。刘文涛等[17]提出了一种利用正解自动构造平台定位误差矢量和平台结构误差Jacobian矩阵的“鸡尾酒法”。Mansour等[18]以Hexaglide型并联机构为研究对象,通过百分表测量相同位置不同姿态下的关节变量向量和定位误差完成了标定。David[19]针对六自由度并联机构提出了基于末端误差最小子集检测信息的运动学标定算法。

本文将6-PTRT并联机构引入CTS试验系统,作为CTS试验系统的运动机构。为了避免CTS并联机构静平台的安装误差对标定结果产生影响,在风洞试验段内安装CTS并联机构后,采用可移动的ROMER绝对臂测量机对机构进行标定,标定完成后CTS并联机构不予拆除或重新安装,直接开展风洞试验。在空间受限的风洞试验段环境中,CTS并联机构标定需要解决的难题有:动平台相对风洞坐标系的空间位姿测量与计算、CTS并联机构标定模型的建立以及CTS并联机构结构参数的辨识。本文针对风洞环境中的CTS并联机构,提出CTS并联机构动平台空间位姿的测量和计算方法,在CTS并联机构的运动学反解条件中加入直线驱动平台与风洞轴线夹角修正以构造标定模型,采用非线性最小二乘法辨识结构参数,最后对辨识后的并联机构定位准度进行校核和风洞试验验证。

1 CTS并联机构的地面标定方法

图 1所示,CTS并联机构包含6条支链,每条支链的组成形式是:移动副-虎克铰-旋转副-虎克铰。CTS并联机构的直线驱动平台内埋于风洞试验段的两侧壁内,6根支撑杆通过虎克铰分别与直线驱动平台和动平台连接,由交流伺服电机驱动直线驱动平台在导轨上滑动实现并联机构运动。

图 1 CTS并联机构简图 Fig. 1 Diagram of CTS parallel mechanism

CTS并联机构的地面标定方法是通过实测动平台位姿辨识真实结构参数,以提高CTS并联机构定位准度的方法,其内容包含动平台位姿的测量和计算、标定模型的建立和结构参数辨识3个部分。CTS并联机构地面标定的流程图如图 2所示。

图 2 标定流程图 Fig. 2 Flowchart of calibration process
1.1 动平台位姿的测量和计算

动平台空间位姿$ \boldsymbol{x}=\left[\begin{array}{llll} {x} & {y} & {z} & {\psi} & {\vartheta} \end{array}\right]$既是CTS并联机构结构参数辨识的输入量,又是CTS并联机构定位准度的量化指标。

动平台相对风洞的姿态角以3个欧拉角($ \psi、\vartheta、\gamma $)表示,风洞坐标系$O_{b}-X_{b} Y_{b} Z_{b} $可依次旋转$ \psi、\vartheta、\gamma $转换到动平台坐标系$ O_{a}-X_{a} Y_{a} Z_{a}$,如图 3所示。

图 3 风洞坐标系与动平台坐标系转换 Fig. 3 Transformation between wind tunnel and moving platform coordinate systems

设计专门的十字校准架作为动平台位姿的测量平台,十字校准架与动平台刚性连接,以保证十字校准架与动平台的相对位置精度,且通过地面加载实验表明,十字校准架与模型的质量差异对CTS并联机构的定位精度影响可以忽略。十字校准架上有4个标准测量孔,如图 4所示。以孔1作为动平台的位移参考点,十字校准架的姿态角作为动平台的姿态角。通过测量孔1~孔4的坐标值计算动平台的位移和姿态角。

图 4 十字校准架 Fig. 4 Cross-shape calibration frame

考虑风洞洞壁的空间限制,动平台位姿的测量设备采用瑞士海克斯康公司的便携式ROMER绝对臂测量机(如图 5所示),通过手持触发方式实现三维空间坐标的测量。ROMER绝对臂测量机的测量范围为2.0 m,空间长度和单点重复测量精度分别为±0.023 mm和±0.016 mm。

图 5 ROMER绝对臂测量机 Fig. 5 ROMER absolute arm

通过ROMER绝对臂测量机测量风洞3个垂直平面的法向方向构造动平台位姿的参考坐标系,即风洞坐标系$O_{b}-X_{b} Y_{b} Z_{b} $。通过ROMER绝对臂测量机点触十字校准架的标准测量孔,获取标准测量孔1~孔4在风洞坐标系下的坐标,并以此构造动平台坐标系$O_{a}-X_{a} Y_{a} Z_{a} $:原点位于孔1,$O_{a}X_{a} $轴由孔1指向孔2,$O_{a}Z_{a} $轴由孔4指向孔3,$O_{a}Y_{a} $方向按照右手直角坐标系确定。

通过孔1~孔4的坐标值,可以方便获得$O_{a}X_{a} $轴与$O_{b}Y_{b} $轴的夹角${\lambda _\vartheta } $$O_{a}X_{a} $轴与$O_{b}Z_{b} $的夹角$ {\lambda _\psi }$$O_{a}Z_{a} $轴与$O_{b}Y_{b} $轴的夹角$ {\lambda _\gamma }$,如图 6所示。基于夹角${\lambda _\theta } $$ {\lambda _\psi }$$ {\lambda _\gamma }$计算姿态角$\vartheta 、\psi、\gamma $的方法为

图 6 两坐标系夹角 Fig. 6 Included angles of two coordinate systems

1) 俯仰角$\vartheta $

根据$O_{a}X_{a} $轴和$O_{b}Y_{b} $轴的方向余弦关系:

$ \cos {\lambda _\vartheta } = \sin \vartheta $

从而得到

$ \vartheta = \arcsin \left( {\cos {\lambda _\theta }} \right) $ (1)
 

2) 偏航角$\psi $

根据$O_{a}X_{a} $轴与$O_{b}Z_{b} $轴的方向余弦关系:

$ \cos {\lambda _\psi } = - \cos \vartheta \sin \psi $

从而得到

$ \psi = - \arcsin \left( {\cos {\lambda _\psi }/\cos \vartheta } \right) $ (2)
 

3) 滚转角$ \gamma $

根据$O_{a}Z_{a} $轴与$O_{b}Y_{b} $轴的方向余弦关系:

$ \cos {\lambda _\gamma } = - \cos \vartheta \sin \gamma $

从而得到

$ \gamma = - \arcsin \left( {\cos {\lambda _\gamma }/ cos \vartheta } \right) $ (3)
 
1.2 标定模型的建立

CTS并联机构的标定模型基于直线驱动平台移动距离和动平台位姿的运动学关系,是CTS并联机构各部件空间位置在每个动平台位姿下的数学关系,其建立涉及到风洞坐标系下支撑杆长、静平台虎克铰中心坐标和动平台虎克铰中心坐标的获取。

CTS并联机构的直线驱动平台移动距离$\mathit{\boldsymbol{d}} = \left[ {\begin{array}{*{20}{l}} {{d_1}}&{{d_2}}&{{d_3}}&{{d_4}}&{{d_5}}&{{d_6}} \end{array}} \right] $和动平台位姿$ \mathit{\boldsymbol{x}}=\left[ {\begin{array}{*{20}{c}} x&y&z&\psi &\vartheta &\gamma \end{array}} \right]$的运动学关系由CTS并联机构的运动学反解条件决定[20]。CTS并联机构的反解条件为:已知动平台的位姿,存在唯一的直线驱动平台移动距离满足CTS并联机构结构的几何关系。对于CTS并联机构的任一条支链,其结构矢量图如图 7所示,在风洞坐标系中可描述为

$ \mathit{\boldsymbol{L}}_i^b = {\mathit{\boldsymbol{t}}^b}\mathit{\boldsymbol{ + A}}_i^b - \mathit{\boldsymbol{B}}_i^b $
图 7 单一支链矢量图 Fig. 7 Vector diagram of single chain

标量形式为

$ L_i^2 = {\left( {{\mathit{\boldsymbol{t}}^b} + \mathit{\boldsymbol{A}}_i^b - \mathit{\boldsymbol{B}}_i^b} \right)^{\rm{T}}}\left( {{\mathit{\boldsymbol{t}}^b} + \mathit{\boldsymbol{A}}_i^b - \mathit{\boldsymbol{B}}_i^b} \right) $ (4)
 

式中:tb为风洞坐标系下,风洞坐标系原点Ob到动坐标系原点Oa的向量;Aib为风洞坐标系下,Oa到支链i的动平台虎克铰中心ai的向量;Bib为风洞坐标系下,Ob到支链i的静平台虎克铰中心bi的向量;Lib为风洞坐标系下,支链i的静平台虎克铰中心bi到动平台虎克铰中心ai的向量;Li为支链i的静平台虎克铰中心bi到动平台虎克铰中心ai的距离。

确定式(4)中各变量的数值,即建立了CTS并联机构的标定模型。tb可以从动平台位移参考点坐标直接得到,Li可以从支撑杆杆长、静虎克铰和动虎克铰的几何尺寸得到,Aib是动平台姿态角的函数,Bib是直线驱动平台移动距离的函数。AibBib将通过下面的计算方法得到。

1) Aib的计算方法

已知Aia是动平台虎克铰中心在动坐标系下的坐标表达,由动平台刚体几何尺寸决定。将Aia从动坐标系转换到风洞坐标系可以计算得到Aib,转换矩阵由动平台姿态角组成。所以Aib可表达成

$ \mathit{\boldsymbol{A}}_i^b = \mathit{\boldsymbol{R}} \cdot \mathit{\boldsymbol{A}}_i^a $ (5)
 

式中:

$ \mathit{\boldsymbol{R}} = \left[ {\begin{array}{*{20}{c}} {{\rm{c}}\vartheta {\rm{c}}\psi }&{ - {\rm{s}}\vartheta {\rm{c}}\psi {\rm{c}}\gamma + {\rm{s}}\psi {\rm{s}}\gamma }&{{\rm{s}}\vartheta {\rm{c}}\psi {\rm{s}}\gamma + {\rm{s}}\psi {\rm{c}}\gamma }\\ {{\rm{s}}\vartheta }&{{\rm{c}}\vartheta {\rm{c}}\gamma }&{ - {\rm{c}}\vartheta {\rm{s}}\gamma }\\ { - {\rm{c}}\vartheta s\psi }&{{\rm{s}}\vartheta {\rm{s}}\psi {\rm{c}}\gamma + {\rm{c}}\psi {\rm{s}}\gamma }&{ - {\rm{s}}\vartheta {\rm{s}}\psi {\rm{s}}\gamma + {\rm{c}}\psi {\rm{c}}\gamma } \end{array}} \right] $

式中:sin、cos分别简记为s和c。

2) Bib的计算方法

静平台虎克铰固连在直线驱动平台上,随着直线驱动平台一起运动,由于安装误差,直线驱动平台与风洞坐标系的ObXb轴不可避免地存在夹角,这里采用线性关系表达静平台虎克铰中心在风洞坐标系下的坐标与直线驱动平台移动距离的关系,则Bib向量可以写成

$ \mathit{\boldsymbol{B}}_i^b = \left[ {\begin{array}{*{20}{l}} {{b_{ix}}}\\ {{b_{iy}}}\\ {{b_{ix}}} \end{array}} \right] = \left[ {\begin{array}{*{20}{l}} {{k_x}}\\ {{k_y}}\\ {{k_z}} \end{array}} \right]{d_i} + \left[ {\begin{array}{*{20}{c}} {{s_x}}\\ {{s_y}}\\ {{s_z}} \end{array}} \right] = \mathit{\boldsymbol{K}}{d_i} + \mathit{\boldsymbol{S}} $ (6)
 

式中:di为直线驱动平台移动距离;K为线性方程系数组成的数组;S为线性方程截距组成的数组。

结合式(4)、式(5)和式(6),CTS并联机构标定模型可表示为

$ \begin{array}{l} L_i^2 \equiv {\left( {{\mathit{\boldsymbol{t}}^b} \times \mathit{\boldsymbol{RA}}_i^a - \mathit{\boldsymbol{K}}{d_i} - \mathit{\boldsymbol{S}}} \right)^{\rm{T}}} \cdot \\ \left( {{\mathit{\boldsymbol{t}}^b} \times \mathit{\boldsymbol{RA}}_i^a - \mathit{\boldsymbol{K}}{d_i} - \mathit{\boldsymbol{S}}} \right) \end{array} $ (7)
 
1.3 结构参数辨识

由于存在加工和安装误差,名义结构参数和真实结构参数间不可避免地存在差异。结构参数辨识是基于标定模型辨识真实结构参数的过程。根据标定模型,当实测位姿对应的直线驱动平台移动距离di与名义直线驱动平台移动距离dni相等时,实测位姿将与名义位姿相等。引入目标函数:

$ I = \mathop {\min }\limits_{\left\{ a \right\}} \sum\limits_{i = 1}^6 {{{\left( {{d_i} - {d_{ni}}} \right)}^2}} $ (8)
 

结构参数辨识即辨识CTS并联机构的结构参数{a},使得目标函数值最小。

影响CTS并联机构定位准度的结构参数有:静平台虎克铰中心位置、动平台虎克铰中心位置和支撑杆的长度。在CTS并联机构的任一支链中,静平台端虎克铰中心位置由其在风洞坐标系下的3个坐标值$\left(b_{i x}, b_{i y}, b_{i z}\right) $决定,由式(6)可知,$\left(b_{i x}, b_{i y}, b_{i z}\right) $是直线驱动平台移动距离的函数,可选取直线驱动平台移动距离的误差dei作为1个结构参数;动平台端虎克铰中心位置由其在动坐标下的3个坐标值$ \left(a_{i x}, a_{i y}, a_{i z}\right)$决定,即3个结构参数;杆长Li,即1个结构参数。因此每条支链有5个结构参数,6条支链共有30个结构参数。

动平台位姿作为结构参数辨识的输入量,用于标定的位姿选取遵循以下原则:

1) 在CTS试验所需的位姿空间中(如表 1所示),随机均匀地选取标定位姿,每个位姿均包含6个自由度。

表 1 CTS试验空间 Table 1 Workspace of CTS test
位姿空间 范围
x/mm 0~300
y/mm -300~300
z/mm -300~300
$\vartheta $ /(°) -20~20
ψ/(°) -20~20
γ/(°) 0

2) CTS并联机构共有6条支链,对于30个结构参数,至少需要5个位姿,利用运动雅克比矩阵$ \boldsymbol{J}=\partial \boldsymbol{d} / \partial \boldsymbol{x}$的条件数对奇异性位姿进行甄别和剔除。

3) 保证所选取的位姿都在ROMER绝对臂测量机的可测空间内。

根据以上位姿选取原则,选取130组位姿,选用的标定点的空间情况如图 8所示。

图 8 标定点的位置空间分布 Fig. 8 Spatial distribution of calibration points

在测量和计算标定点位姿的基础上,以CTS并联机构的结构参数设计值作为迭代初值,如表 2所示,采用非线性最小二乘法辨识结构参数。辨识后的结构参数见表 3

表 2 结构参数迭代初值 Table 2 Initial value of structural parameters  
mm
支链 Li aix aiy aiz dei
1 1 000 -133.8 68.6 -50 0
2 900 -133.8 -18 -100 0
3 1 100 -133.8 -104.6 -50 0
4 1 100 -133.8 -104.6 50 0
5 900 -133.8 -18 100 0
6 1 000 -133.8 68.6 50 0
表 3 辨识后结构参数值 Table 3 Post-calibration value of parameters 
mm
支链 Li aix aiy aiz dei
1 999.3 -134.3 69.1 -49.1 -0.4
2 899.5 -134.0 -16.9 -100.0 -1.1
3 1 099.0 -133.9 -104.1 -50.7 0.9
4 1 100.262 -133.686 -105.202 49.1 -2.3
5 899.960 -133.859 -19.168 99.8 -0.8
6 1 000.027 -134.101 67.958 51.0 -1.3

参数辨识完成后,以辨识后的结构参数值控制CTS并联机构运动,进行定位准度的校核和风洞试验验证。

2 定位准度的校核和风洞试验验证

位姿定位准度的评价指标主要有:平均绝对误差(Mean of Absolute Error, MAE)和误差均方根(Root Mean Square, RMS)[21]。MAE可以估算可能的误差范围,RMS则能反映误差波动的幅值信息。MAE、RMS的表达式为

$ {\rm{MAE}} = \frac{1}{n}\sum\limits_i^n {{\rm{ABS}}\left( {{{\hat r}_i} - {r_i}} \right)} $ (9)
 
$ {\rm{RMS}} = \sqrt {\frac{1}{n}\sum\limits_i^n {{{\left( {{e_i} - {{\bar e}_i}} \right)}^2}} } $ (10)
 

式中:$\hat{r}_{i} $为名义位姿;ri为实测位姿;n为校验点数;ei为单点误差;${{\bar e}_i} $为平均误差。

校验点的选取原则与标定点的选取原则相同,选取10组位姿,对CTS并联机构的位姿定位准度进行辨识前后的对比。

辨识前后定位准度的MAE指标对比结果如表 4所示,RMS指标对比结果如表 5所示。

表 4 辨识前后的定位准度MAE对比 Table 4 Comparison of MAE for pre-calibration and post-calibration
参数 辨识前MAE 辨识后MAE 比值/%
x/mm 1.022 0.046 4.5
y/mm 0.608 0.081 13.4
z/mm 0.175 0.054 31.2
$\vartheta $ /(°) 0.064 0.028 43.9
ψ/(°) 0.743 0.030 4.1
γ/(°) 1.361 0.022 1.6
注:比值=(辨识后指标/辨识前指标),表 5同。
表 5 辨识前后的定位准度RMS对比 Table 5 Comparison of RMS for pre-calibration and post-calibration
参数 辨识前RMS 辨识后RMS 比值/%
x/mm 0.626 0.027 4.3
y/mm 0.413 0.013 3.3
z/mm 0.137 0.031 22.5
$\vartheta $ /(°) 0.059 0.007 12.4
ψ/(°) 0.087 0.013 15.6
γ/(°) 0.152 0.010 1.6

表 4可以看出,辨识后的定位误差范围小于辨识前的定位误差范围,尤其是xyψγ这4个自由度;从表 5可以看出,辨识前位姿定位误差波动明显,辨识后6个自由度的位姿定位误差都仅有小幅度的波动,表明各自由度的定位准度稳定。

辨识后CTS并联机构的定位准度与CTS试验准度要求对比如表 6所示,可以看出,辨识后并联机构定位准度满足CTS试验准度要求。

表 6 CTS并联机构准度与CTS试验要求准度对比 Table 6 Comparison of position accuracy between CTS test requirement and the CTS parallel mechanism
参数 试验要求准度 并联机构准度
x/mm ±0.1 -0.087~0.061
y/mm ±0.1 -0.095~-0.060
z/mm ±0.1 -0.087~0.078
$ \vartheta $ /(°) ±0.05 -0.035~0.016
ψ/(°) ±0.05 -0.042~0.050
γ/(°) ±0.05 -0.039~0.021

上述对CTS并联机构的定位准度校核是在地面无风载情况下进行的,为了考核CTS并联机构在风载情况下的定位准度,在FD12风洞[1]中进行了并联机构与常规攻角机构的对比风洞试验,在相同测控系统和试验模型条件下,风洞试验结果的精准度由机构的定位准度决定。通常以标准模型按阶梯攻角运动的测力风洞试验考核运动机构的风载定位准度,采用的常规攻角机构刚度极大且已通过定位准度标定校验,将其得出的气动数据作为真值。试验模型为8#标模(见图 9),试验马赫数Ma=0.8,试验攻角序列:-10°~10°。并联机构动平台端的8#标模在试验段均匀流场中的3处位置进行了风洞测力试验,用以验证并联机构在不同位置的定位准度。

图 9 8#标模照 Fig. 9 8# standard model

试验后,法向力系数CN和俯仰力矩系数MZ曲线对比如图 10所示,图上#1、#2和#3分别表示模型质心位于试验段迎风横截面中心、中心正上方和正下方200 mm位置。试验数据表明,并联机构支撑与常规攻角机构支撑的8#标模的CN均方差σCN小于0.005,MZ均方差σMZ小于0.005,满足了测力试验精准度的要求[1]

图 10 并联机构和常规攻角机构的风洞试验数据对比 Fig. 10 Comparison of wind tunnel test data between parallel mechanism and regular attack angle mechanism
3 结论

1) 6-PTRT并联构型可以用于CTS试验的运动机构,并能有效提高CTS试验的定位准度,以其支撑8#标模的风洞试验数据满足测力风洞试验精准度要求。

2) 提出的针对空间受限风洞环境的动平台空间位姿的测量和计算方法,可以准确高效地得到动平台空间位姿,建立的包含直线驱动平台与风洞轴线夹角修正的标定模型,能准确地描述风洞安装环境下的CTS并联机构的运动学关系。

3) 结构参数辨识后的CTS并联机构位姿误差范围和误差波动幅值均显著减小,位移准度优于±0.1 mm,姿态角准度优于±0.05°。

致谢

感谢秦永明研究员、王强研究员、张江研究员、董金刚高级工程师在本文工作开展中给予的指导与帮助!

参考文献
[1] 李周复. 风洞试验手册[M]. 北京: 航空工业出版社, 2015: 554-581.
LI Z F. Wind tunnel handbook[M]. Beijing: Aviation Industry Press, 2015: 554-581. (in Chinese)
[2] 黄叙辉, 罗新福, 于志松. FL-24风洞新型捕获轨迹系统设计与发展[J]. 空气动力学报, 2008, 26(2): 146-149.
HUANG X H, LUO X F, YU Z S. Design & development of a new captive trajectory simulation system in FL-24 wind tunnel[J]. Acta Aerodynamic Sinica, 2008, 26(2): 146-149. (in Chinese)
Cited By in Cnki (15) | Click to display the text
[3] MOSHE Z, MOTY S. The use of the captive trajectory system for computation of trajectories to the impact point: AIAA-92-4031[R]. Reston, VA: AIAA, 1992.
[4] MICHAEL J H. 战术导弹空气动力学[M]. 北京: 宇航出版社, 1999: 578-662.
MICHAEL J H. Tactical missile aerodynamics:General topics[M]. Beijing: Astronavigation Press, 1999: 578-662. (in Chinese)
[5] 黄叙辉, 庞旭东, 宋斌. 1.2 m跨超声速风洞新型捕获轨迹系统研制[J]. 实验流体力学, 2008, 22(2): 95-98.
HUANG X H, PANG X D, SONG B. Development of a new captive trajectory simulation system in the 1.2 m transonic wind tunnel[J]. Journal of Experiments in Fluid Mechanics, 2008, 22(2): 95-98. (in Chinese)
Cited By in Cnki | Click to display the text
[6] LOMBARDI G. Use of a captive trajectory system in a wind tunnel: AGARD CP 570[R]. 1996.
[7] 张曙, HEISEL U. 并联运动机床[M]. 北京: 机械工业出版社, 2003: 18-34.
ZHANG S, HEISEL U. Parallel kinematics machine tool[M]. Beijing: Machinery Industry Press, 2003: 18-34. (in Chinese)
[8] 樊锐, 王进, 侯立果. 基于拆分测量的6-UPS并联机构精度的标定[J]. 黑龙江科技大学学报, 2018, 28(4): 426-429.
FAN R, WANG J, HOU L G. Calibration of 6-UPS parallel mechanism based on split measurement[J]. Journal of Heilongjiang University of Science & Technology, 2018, 28(4): 426-429. (in Chinese)
Cited By in Cnki | Click to display the text
[9] 曾磊.平面3-RRR并联机构标定方法与控制系统研究[D].广州: 华南理工大学, 2017: 16-25.
ZENG L. Research on the calibration method and control system of planar 3-RRR parallel mechanism[D]. Guangzhou: South China University of Technology, 2017: 16-25(in Chinese).
Cited By in Cnki | Click to display the text
[10] 戴智武, 刘超, 盛鑫军, 等. Delta并联机器人机构参数误差分析及标定方法研究[J]. 机电一体化, 2016(3): 8-12.
DAI Z W, LIU C, SHENG X J, et al. The error analysis and calibration of delta parallel robot[J]. Mechatronics, 2016(3): 8-12. (in Chinese)
Cited By in Cnki | Click to display the text
[11] 李寅翔. 2-UPR-RPU并联机构运动学标定研究[D].杭州: 浙江理工大学, 2015: 23-50.
LI Y X. Kinematic calibration of a 2-UPR-RPU parallel manipulator[D]. Hangzhou: Zhejiang Sci-Tech University, 2015: 23-50(in Chinese).
Cited By in Cnki | Click to display the text
[12] 马志强, 李泷杲, 刑宏文, 等. 3-PPPS并联机翼调姿机构运动学标定[J]. 计算机集成制造系统, 2015, 21(9): 2378-2383.
MA Z Q, LI L G, XING H W, et al. Kinematic calibration of 3-PPPS parallel wing posture adjustment mechanism[J]. Computer Integrated Manufacturing Systems, 2015, 21(9): 2378-2383. (in Chinese)
Cited By in Cnki | Click to display the text
[13] 张立杰, 李立泉, 王艮川. 基于三坐标测量机的球面5R并联机构运动学标定研究[J]. 中国机械工程, 2013, 24(22): 2997-3002.
ZHANG L J, LI L Q, WANG G C. Research on kinematic calibration of spherical 5R parallel manipulator based on coordinate measuring machine[J]. China Mechanical Engineering, 2013, 24(22): 2997-3002. (in Chinese)
Cited By in Cnki (3) | Click to display the text
[14] 皮阳军, 王宣银, 胡玉梅. 基于关节力传感器的并联六自由度机构标定方法[J]. 农业机械学报, 2012, 43(10): 215-218.
PI Y J, WANG X Y, HU Y M. Calibration of 6-DOF parallel mechanism using joint force sensors[J]. Transactions of the Chinese Society for Agricultural Machinery, 2012, 43(10): 215-218. (in Chinese)
Cited By in Cnki (4) | Click to display the text
[15] 张建中.基于倾角测量位姿的六自由度平台运动学标定研究[D].哈尔滨: 哈尔滨工业大学, 2007: 33-44.
ZHANG J Z. Kinematics calibration of stewart platforms based on configurations measured with two inclinometers[D]. Harbin: Harbin Institute of Technology, 2007: 33-44(in Chinese).
[16] 孙皓. 5-UPS/PRPU并联机器人精度设计与运动学参数辨识[D].天津: 天津理工大学, 2009: 23-33.
SUN H. Research on the accuracy design and kinematics parameters calibration of 5-UPS/PRPU parallel manipulator[D]. Tianjin: Tianjin University of Technology, 2009: 23-33(in Chinese).
Cited By in Cnki (3) | Click to display the text
[17] 刘文涛, 唐德威, 王知行. Stewart平台机构标定的鸡尾酒法[J]. 机械工程学报, 2004, 40(12): 48-52.
LIU W T, TANG D W, WANG Z X. Cocktail method for stewart platform calibration[J]. Journal of Mechanical Engineering, 2004, 40(12): 48-52. (in Chinese)
Cited By in Cnki (41) | Click to display the text
[18] MANSOUR A, HODJAT P, ARIA A, et al. Experimental kinematic calibration of parallel manipulators using a relative position error measurement system[J]. Robotics and Computer Integrated Manufacturing, 2010, 26(6): 799-804.
Click to display the text
[19] DAVID D. Kinematic calibration of the gough platform[J]. Robotica, 2003, 21(6): 677-690.
Click to display the text
[20] 李鹭扬. 6-SPS型并联机床若干关键理论研究[D].南京: 南京航空航天大学, 2011: 18-20.
LI L Y. Research on some key theoretical issues of 6-SPS parallel kinematic machines[D]. Nanjing: Nanjing University of Aeronautics and Astronautics, 2011: 18-20(in Chinese).
[21] 盛骤, 谢式千, 潘承毅. 概率论与数理统计[M]. 北京: 高等教育出版社, 2008: 90-113.
SHENG Z, XIE S Q, PAN C Y. Probability theory and mathematical statistics[M]. Beijing: Higher Education Press, 2008: 90-113. (in Chinese)
http://dx.doi.org/10.7527/S1000-6893.2019.23175
中国航空学会和北京航空航天大学主办。
0

文章信息

谢峰, 洪冠新, 张晨凯, 魏忠武, 马汉东
XIE Feng, HONG Guanxin, ZHANG Chenkai, WEI Zhongwu, MA Handong
捕获轨迹系统并联机构地面标定方法
Ground calibration method of captive trajectory system parallel mechanism
航空学报, 2020, 41(1): 423175.
Acta Aeronautica et Astronautica Sinica, 2020, 41(1): 423175.
http://dx.doi.org/10.7527/S1000-6893.2019.23175

文章历史

收稿日期: 2019-05-21
退修日期: 2019-06-12
录用日期: 2019-09-10
网络出版时间: 2019-10-12 15:25

相关文章

工作空间