一种动态离心机加速度场的高精度测算方法技术

本发明专利技术涉及一种动态离心机加速度场的高精度测算方法，属于力学环境试验与测量技术领域，在动态离心机可视为刚体的前提下，根据刚体上各点加速度关于位置线性变化的原理，利用动态离心机上多个点的加速度实测值可解算出该线性变化的具体系数，从而得到动态离心机上试件敏感位置处的加速度测算值；再对此计算过程对应的数学公式进行偏导分析，可得到该测算值的不确定度。利用本方法可设计出结构简单、便于安装、操作方便、精度高的加速度场测算系统，适用于单轴/双轴/三轴等各类动态离心机上加速度场的高精度测算。加速度场的高精度测算。加速度场的高精度测算。

【技术实现步骤摘要】
一种动态离心机加速度场的高精度测算方法


[0001]本专利技术属于力学环境试验与测量
，具体涉及一种动态离心机加速度场的高精度测算方法。

技术介绍

[0002]加表等惯性器件是飞行器上的重要组件，其精度对这些飞行器的导航精度具有重要影响。因此，过去人们利用精密离心机创造高精度的加速度环境，从而对这些惯性器件进行考核。考虑到精密离心机无法模拟实际的飞行加速度环境——三维动态加速度环境，中国工程物理研究院总体工程研究所研制了三轴动态离心机，为上述惯性器件的考核提供了更加逼真的飞行加速度环境。目前，该设备能够以
±
0.4g的加载精度在设备控制点处产生三维动态加速度环境。然而，一方面该精度还不够高，无法对高精度加表进行精度验证；另一方面，试验中加表等试件一般无法准确安装在设备控制点处，考虑到动态离心机上加速度场具有一定梯度，各点加速度不同，因此控制点处的加速度值不能直接用以表示试件处的加速度值。因此，亟需发展动态离心机上加速度场的高精度测算方法，从而得到试件处加速度的高精度测算值，同时给出该测算值的不确定度，从而验证加表等惯性器件的精度性能。此外，发展本方法后，也可对动态离心机控制点处的加速度进行高精度测算，从而实现设备的校准，确保动态离心试验设备提供的加速度环境条件量值真实可靠。
[0003]目前，现有动态离心机上加速度场测算方法的大体思路如下：利用多刚体运动学给出加速度计算的数学模型，然后测量动态离心机各轴的运动参数、尺寸以及试件安装位置、姿态，最后带入该数学模型计算试件处的加速度。由于这种方式依赖于动态离心机的各种参数，因此将其简称为“参数法”。
[0004]采用参数法时，为了对加速度进行高精度测算，需考虑离心机主轴动态轴偏、离心机主轴与外框转轴的平行度误差、外框转轴与内框转轴的垂直度误差、动态离心机的动态半径等因素，试件处加速度表达式的函数形式如下。
[0005]a
M
＝f(ω1,α1,q2,ω2,α2,q3,ω3,α3,ρ
0M
,A
sn
,R,ΔR,Δθ1,Δθ2,Δθ3)
ꢀꢀ
(1)
[0006]式中，ω1,α1分别表示主轴的角速度和角加速度，q2,ω2,α2分别表示外框的转角、角速度和角加速度，q3,ω3,α3分别表示内框的转角、角速度和角加速度，ρ
0M
,A
sn
分别表示试件敏感位置相对内框中心点的位置、姿态，R表示内框中心点距离主轴的静态距离(即静态半径)，ΔR为离心机运动过程中R的改变量，Δθ1,Δθ2,Δθ3分别表示动态离心机运转过程中主轴、外框转轴、内框转轴的动态姿态偏转量，决定了离心机主轴动态轴偏、离心机主轴与外框转轴的平行度误差、外框转轴与内框转轴的垂直度误差等误差量。可见采用参数法对试件处的加速度进行高精度测算时，面临如下难题：
[0007]1)数学建模异常复杂：参数多，多达22个，这会导致更加复杂的不确定度分析公式；
[0008]2)动态离心机运转过程中的动态半径ΔR，各轴的动态姿态偏转Δθ1,Δθ2,Δθ3均很难进行高精度的定量分析，且由于缺乏测量基准而很难测量；
[0009]3)试件安装位置ρ
0M
、姿态A
sn
均由于缺乏测量基准而难以得到高精度测量值；
[0010]4)各轴角速度、角加速度的高精度测量难度高：考虑到离心机编码器的安装误差，导致编码器的刻度具有一定的不均匀性，导致各轴角速度、角加速度的测量误差难以精确分析。

技术实现思路

[0011]本专利技术目的在于提供一种动态离心机加速度场的高精度测算方法，用于解决上述现有技术中存在的技术问题，即，采用参数法对动态离心机上加速度场进行高精度测算时，其不确定度分析同时面临理论建模异常复杂、工程实施难度极高两大难题。
[0012]为实现上述目的，本专利技术的技术方案是：
[0013]在动态离心机可视为刚体的前提下，根据刚体上各点加速度关于位置线性变化的原理，利用动态离心机上多个点的加速度实测值可解算出该线性变化的具体系数，从而得到动态离心机上试件敏感位置处的加速度测算值；再对此计算过程对应的数学公式进行偏导分析，可得到该测算值的不确定度。
[0014]本方法的特点：考虑到各类动态离心机均可视为刚体，因此本测算方法同时适用于单轴/双轴/三轴等各类动态离心机上加速度场的高精度测算。
[0015]进一步的，动态离心机上各点加速度关于位置ρ
0M
线性变化，可归结为如下公式：a
M
(t)＝a0(t)+K(t)ρ
0M
(2)式中a
M
为点M处的加速度，a0为点0处的加速度，ρ
0M
为0到M的矢径，K∈R3×3为刚体的加速度场分布系数矩阵；因此，利用动态离心机上多个点处的加速度实测值即可解算出K，进而将K代入式(2)即可得到动态离心机上任意一点M处的加速度测算值。
[0016]进一步的，根据动态离心机上多个点的加速度实测值解算K；解算K的核心算法如下：
[0017][0018][0019][0020]式中，a
nx
,a
ny
,a
nz
(n＝0,i,j,k)分别表示加速度计(0、i、j、k)的测量结果；ρ
0i
,ρ
0j
,ρ
0k
分别表示加速度计0敏感位置到加速度计i、j、k敏感位置的位置矢量在方块坐标系f(图1中的坐标系o
‑
xyz)中的投影；上标“T”表示矩阵转置；上标
“‑
1”表示矩阵的逆，k
11
,...,k
33
为K的元素；
[0021]然后根据试件坐标系s(表示试件敏感元件的位置、姿态)相对f系的位置、姿态即
可得到试件敏感位置处的加速度测算值如下：
[0022][0023]式中，A
sf
表示f系到s系的方向余弦矩阵，a0表示加速度计0的测量结果，ρ
0M
表示试件敏感位置M相对加速度计0敏感位置的位置矢量，即投影到f系；可见加速度测算值由如下测量值决定：
[0024]1)加速度计0、i、j、k的测量值，a
0x
、a
0y
、a
0z
、a
ix
、a
iy
、a
iz
、a
jx
、a
jy
、a
jz
、a
kx
、a
ky
、a
kz
；
[0025]2)加速度计0、i、j、k相对f系的位置坐标，x0、y0、z0、x
i
、y
i
、z
i
、x
j
、y
j
、z
j
...

【技术保护点】

【技术特征摘要】
1.一种动态离心机加速度场的高精度测算方法，其特征在于，利用加速度场分布系数矩阵K描述动态离心机上加速度场关于位置的分布规律，然后根据某个基点处的加速度实测值推算试件敏感位置处的加速度。2.如权利要求1所述的一种动态离心机加速度场的高精度测算方法，其特征在于，动态离心机上各点加速度关于位置ρ
0M
线性变化，可归结为如下公式：a
M
(t)＝a0(t)+K(t)ρ
0M
ꢀꢀ
(2)式中a
M
为点M处的加速度，a0为点0处的加速度，ρ
0M
为0到M的矢径，K∈R3×3为刚体的加速度场分布系数矩阵；因此，利用动态离心机上多个点处的加速度实测值即可解算出K，进而将K代入式(2)即可得到动态离心机上任意一点M处的加速度测算值。3.如权利要求2所述的一种动态离心机加速度场的高精度测算方法，其特征在于，根据动态离心机上多个点的加速度实测值解算K；解算K的核心算法如下：动态离心机上多个点的加速度实测值解算K；解算K的核心算法如下：动态离心机上多个点的加速度实测值解算K；解算K的核心算法如下：式中，a
nx
,a
ny
,a
nz
(n＝0,i,j,k)分别表示加速度计(0、i、j、k)的测量结果；ρ
0i
,ρ
0j
,ρ
0k
分别表示加速度计0敏感位置到加速度计i、j、k敏感位置的位置矢量在方块坐标系f中的投影；上标“T”表示矩阵转置；上标
“‑
1”表示矩阵的逆，k
11
,...,k
33
为K的元素。4.如权利要求3所述的一种动态离心机加速度场的高精度测算方法，其特征在于，求得K后，试件敏感位置处的加速度测算值计算方法如下：根据试件坐标系s相对f系的位置、姿态即可得到试件敏感位置处的加速度测算值为：式中，A
sf
表示f系到s系的方向余弦矩阵，a0表示加速度计0的测量结果，ρ
0M
表示试件敏感位置M相对加速度计0敏感位置的位置矢量，即投影到f系；可见加速度测算值由如下测量值决定：1)加速度计0、i、j、k的测量值，a
0x
、a
0y
、a
0z
、a
ix
、a
iy
、a
iz
、a
jx
、a
jy
、a
jz
、a

【专利技术属性】
技术研发人员：周邦召，刘谦，杨敏，陈胜来，毛勇建，陈颖，牛宝良，袁彪，欧峰，张帅，吴静，张东峰，
申请(专利权)人：中国工程物理研究院总体工程研究所，
类型：发明
国别省市：