CN112698055B - 加速度计在精密离心机上的参数标定方法 - Google Patents

加速度计在精密离心机上的参数标定方法 Download PDF

Info

Publication number
CN112698055B
CN112698055B CN202110313204.0A CN202110313204A CN112698055B CN 112698055 B CN112698055 B CN 112698055B CN 202110313204 A CN202110313204 A CN 202110313204A CN 112698055 B CN112698055 B CN 112698055B
Authority
CN
China
Prior art keywords
coordinate system
accelerometer
error
axis
shaft
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110313204.0A
Other languages
English (en)
Other versions
CN112698055A (zh
Inventor
王常虹
夏红伟
刘庆博
任顺清
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenrui Technology Beijing Co ltd
Harbin Institute of Technology
Original Assignee
Shenrui Technology Beijing Co ltd
Harbin Institute of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shenrui Technology Beijing Co ltd, Harbin Institute of Technology filed Critical Shenrui Technology Beijing Co ltd
Priority to CN202110313204.0A priority Critical patent/CN112698055B/zh
Publication of CN112698055A publication Critical patent/CN112698055A/zh
Application granted granted Critical
Publication of CN112698055B publication Critical patent/CN112698055B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P21/00—Testing or calibrating of apparatus or devices covered by the preceding groups

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Length Measuring Devices With Unspecified Measuring Means (AREA)
  • Centrifugal Separators (AREA)

Abstract

本发明公开一种加速度计在精密离心机上的参数标定方法，包括：获取精密离心机的各静态误差以及动态误差，并根据精密离心机的结构建立坐标系，以及根据所述各静态误差以及动态误差计算所述坐标系下的位姿误差；驱动精密离心机的主轴以匀角速率旋转，以产生向心加速度标定加速度计，基于所述坐标系下的位姿误差计算向心加速度、重力加速度和Coriolis加速度的比力分配，以确定加速度计误差模型；对加速度计在三种不同安装方式下的六个对称位置的指示输出，利用加减消元的方法标定加速度计误差模型表达式中的高阶项误差系数。本发明可有效提高石英加速度计高阶误差模型系数的标定精度。

Description

加速度计在精密离心机上的参数标定方法
技术领域
本发明涉及离心机标定领域，具体涉及一种加速度计在精密离心机上的参数标定方法。
背景技术
文献“加速度计精密离心机试验的优化设计”分析了加速度计在精密离心机测试时的实际量测噪声特性，在此基础上指出传统的优化设计方法，即饱和D最优试验设计，存在工程适用性问题。然后为了改善饱和D最优试验设计的适用性，并且考虑到试验代价和精度的折中关系，提出了D最优改进试验设计方案。该方案将饱和D最优试验谱点作为基本谱点，在基本谱点之间均匀***其他谱点来降低输入加速度偏差的影响，并通过加权的方法来分配基本谱点和新增谱点的测度，权值的选取依据实际的噪声特性。虽然文献“加速度计精密离心机试验的优化设计”对石英加速度计在精密离心机上进行了具体的标定试验，但没有考虑离心机误差对误差模型系数标定精度的影响，这可能会引入额外的标定误差，并且文献中所辨识的加速度计的误差模型系数较少。
文献“精密离心机误差对石英加速度计误差标定精度分析”分析了离心机各个误差源，用齐次变换法精确地计算了产生的向心加速度，给出了向心加速度、重力加速度和哥氏加速度在加速度计坐标系下的分量，推导了被试加速度计输入加速度的精确表达式。采用了10位置测试方法来辨识误差模型的高阶系数，着重讨论了误差模型系数的计算值与离心机误差之间的关系。但是二次项误差系数
Figure 103085DEST_PATH_IMAGE001
和三次项误差系数
Figure 489067DEST_PATH_IMAGE002
、
Figure 940908DEST_PATH_IMAGE003
未得到辨识，并且需要已知动态和静态误差对辨识结果进行修正和补偿，无法规避离心机的各项误差。
发明内容
有鉴于此，本发明提供一种加速度计在精密离心机上的参数标定方法，包括：
获取精密离心机的各静态误差以及动态误差，并根据精密离心机的结构建立坐标系，以及根据所述各静态误差以及动态误差计算所述坐标系下的位姿误差；
驱动精密离心机的主轴以匀角速率旋转，以产生向心加速度标定加速度计，基于所述坐标系下的位姿误差计算向心加速度、重力加速度和Coriolis加速度的比力分配，以确定加速度计误差模型；
对加速度计在三种不同安装方式下的六个对称位置的指示输出，利用加减消元的方法标定加速度计误差模型表达式中的高阶项误差系数。
本发明一种加速度计在精密离心机上的参数标定方法，在分析精密离心机各项动、静态误差源的基础上，给出了离心机输入比力的精确表达式；结合加速度计误差模型，利用加减消元的方法标定加速度计误差模型表达式中的高阶项误差系数，在离心机误差稳定的情况下，监测和补偿动态失准角和动态半径，就可以完全消除离心机的动态误差和静态误差，可有效提高石英加速度计高阶误差模型系数的标定精度。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其它的附图。
图1为本发明精密离心机结构示意图。
图2为本发明精密离心机各个坐标系示意图。
图3为本发明加速度计3种不同安装方式下的6个对称位置组合。
具体实施方式
下面结合附图对本发明实施例进行详细描述。
需说明的是，在不冲突的情况下，以下实施例及实施例中的特征可以相互组合；并且，基于本公开中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本公开保护的范围。
需要说明的是，下文描述在所附权利要求书的范围内的实施例的各种方面。应显而易见，本文中所描述的方面可体现于广泛多种形式中，且本文中所描述的任何特定结构及/或功能仅为说明性的。基于本公开，所属领域的技术人员应了解，本文中所描述的一个方面可与任何其它方面独立地实施，且可以各种方式组合这些方面中的两者或两者以上。举例来说，可使用本文中所阐述的任何数目个方面来实施设备及/或实践方法。另外，可使用除了本文中所阐述的方面中的一或多者之外的其它结构及/或功能性实施此设备及/或实践此方法。
如图1所示，精密离心机有主轴、水平轴和方位轴3个轴系，3个轴系均有精密位置功能，水平轴轴端安装有360齿多齿分度盘，可以
Figure 601696DEST_PATH_IMAGE004
的精度定位到360个位置，主轴轴系与方位轴轴系均有精密角速率功能，当主轴以
Figure 884910DEST_PATH_IMAGE005
的匀角速率旋转时，在工作半径
Figure 674749DEST_PATH_IMAGE006
处，将产生
Figure 145045DEST_PATH_IMAGE007
的向心加速度。
离心机的静态误差源主要包括主轴轴线的二维铅垂度误差
Figure 812787DEST_PATH_IMAGE008
水平轴轴线与主轴轴线的垂直度
Figure 747245DEST_PATH_IMAGE009
相交度
Figure 209450DEST_PATH_IMAGE010
水平轴轴线与方位轴轴线的垂直度
Figure 167042DEST_PATH_IMAGE011
相交度
Figure 372895DEST_PATH_IMAGE012
以及方位轴的初始零位误差
Figure 161859DEST_PATH_IMAGE013
安装惯性仪表的工作基面对方位轴轴线的垂直度
Figure 60545DEST_PATH_IMAGE014
加速度计安装基面姿态误差
Figure 505433DEST_PATH_IMAGE015
偏心误差
Figure 750863DEST_PATH_IMAGE016
以及初始对零误差
Figure 659913DEST_PATH_IMAGE017
主轴、水平轴和方位轴三个轴的角位置误差分别为
Figure 463921DEST_PATH_IMAGE018
等。图1和图2标出了离心机结构简图以及建立的相应坐标系。
离心机的动态误差源主要包括主轴径向回转误差
Figure 661684DEST_PATH_IMAGE019
轴向窜动
Figure 943761DEST_PATH_IMAGE020
及倾角回转误差
Figure 910580DEST_PATH_IMAGE021
动态半径误差
Figure 947806DEST_PATH_IMAGE022
动态失准角
Figure 101707DEST_PATH_IMAGE023
水平轴径向回转误差
Figure 453054DEST_PATH_IMAGE024
轴向窜动
Figure 274379DEST_PATH_IMAGE025
以及倾角回转误差
Figure 184304DEST_PATH_IMAGE026
方位轴径向回转误差
Figure 887818DEST_PATH_IMAGE027
轴向窜动
Figure 511698DEST_PATH_IMAGE028
倾角回转误差
Figure 453109DEST_PATH_IMAGE029
等。
为了方便研究半径误差的影响，将半径的静态误差与动态误差综合，
Figure 35400DEST_PATH_IMAGE030
其中，
Figure 226210DEST_PATH_IMAGE031
为静态半径标称值，是由计量部门标定出的已知量，但半径的静态测试误差
Figure 919359DEST_PATH_IMAGE032
是未知量，
Figure 184118DEST_PATH_IMAGE033
为用双频激光干涉仪监测的离心机在运行状态下的实际工作半径相对于离心机静态半径的变化量，是主轴角速率
Figure 734048DEST_PATH_IMAGE034
的函数。
下面将建立如下坐标系：
(1)地理坐标系
Figure 615417DEST_PATH_IMAGE035
轴水平指东，
Figure 846678DEST_PATH_IMAGE036
轴水平指北，
Figure 998567DEST_PATH_IMAGE037
轴指天，构成右手坐标系。
(2)主轴轴套坐标系
Figure 453819DEST_PATH_IMAGE038
主轴轴套坐标系相对于地理坐标系的位姿为
Figure 822484DEST_PATH_IMAGE039
(3)主轴坐标系
Figure 591856DEST_PATH_IMAGE040
主轴坐标系相对于主轴轴套坐标系的位姿为
Figure 362366DEST_PATH_IMAGE041
其中
Figure 988520DEST_PATH_IMAGE042
表示主轴旋转的角度。
(4)水平轴轴套坐标系
Figure 844480DEST_PATH_IMAGE043
水平轴轴套坐标系相对于主轴坐标系的位姿为
Figure DEST_PATH_IMAGE044
(5)水平轴坐标系
Figure 886386DEST_PATH_IMAGE045
水平轴坐标系相对于水平轴轴套坐标系的位姿为
Figure 42561DEST_PATH_IMAGE046
其中
Figure 308457DEST_PATH_IMAGE047
表示水平轴旋转的角度。
(6)方位轴轴套坐标系
Figure 884669DEST_PATH_IMAGE048
方位轴轴套坐标系相对于水平轴坐标系的位姿为
Figure 261424DEST_PATH_IMAGE049
(7)方位轴坐标系
Figure 740947DEST_PATH_IMAGE050
方位轴坐标系相对于方位轴轴套坐标系的位姿为
Figure 177745DEST_PATH_IMAGE051
其中
Figure 539456DEST_PATH_IMAGE052
表示方位轴旋转的角度。
(8)工作基面坐标系
Figure 719901DEST_PATH_IMAGE053
工作基面坐标系相对于方位轴坐标系的位姿为
Figure 788351DEST_PATH_IMAGE054
其中
Figure 192788DEST_PATH_IMAGE055
为
Figure 245058DEST_PATH_IMAGE056
点相对
Figure 698036DEST_PATH_IMAGE057
点位移。
(9)加速度计坐标系
Figure 653616DEST_PATH_IMAGE058
加速度计坐标系相对于工作基面坐标系的位姿为
Figure 432216DEST_PATH_IMAGE059
其中
Figure 768519DEST_PATH_IMAGE060
为
Figure 290767DEST_PATH_IMAGE061
点相对
Figure 333810DEST_PATH_IMAGE062
点位移。
以上离心机的各个位姿误差均视为小位移和小角度。加速度计坐标系相对于地理坐标系的位姿为
Figure 283311DEST_PATH_IMAGE063
其中
Figure 44594DEST_PATH_IMAGE064
表示加速度计坐标系与地理坐标系之间的姿态变换矩阵，
Figure 370533DEST_PATH_IMAGE065
为加速度计坐标系与地理坐标系的相对位移矢量。
加速度计坐标系相对于主轴坐标系的位姿为
Figure 330399DEST_PATH_IMAGE066
其中
Figure 185222DEST_PATH_IMAGE067
表示加速度计坐标系与主轴坐标系之间的姿态变换矩阵。
加速度计坐标系原点在主轴坐标系下表示为
Figure 197915DEST_PATH_IMAGE068
忽略二阶小量后可得，
Figure DEST_PATH_IMAGE069
Figure 327545DEST_PATH_IMAGE070
Figure 345180DEST_PATH_IMAGE071
和
Figure 636484DEST_PATH_IMAGE072
将在后文用来计算加速度计坐标系原点的精确向心加速度。
具体地，本实施例所涉及的一种加速度计在精密离心机上的参数标定方法，石英加速度计输入比力的计算过程为：
精密离心机当主轴以匀角速率旋转产生的向心加速度标定加速度计时，加速度计的比力输入有3个来源，即向心加速度、重力加速度、Coriolis加速度，可以得出各个加速度来源的比力分配为：
（1）重力加速度产生的比力在被测加速度计三个轴上的分配
设重力加速度在被测加速度计输入轴、摆轴和输出轴上的分量分别为
Figure 372359DEST_PATH_IMAGE073
Figure 40101DEST_PATH_IMAGE074
重力加速度产生的比力在地理坐标系下表示为
Figure 708979DEST_PATH_IMAGE075
，则在加速度计坐标系下表示为
Figure 436764DEST_PATH_IMAGE076
（2）向心加速度在被测加速度计三个轴上的分配
根据上面所述，加速度计坐标原点处的向心加速度在主轴坐标系下表示为
Figure 394356DEST_PATH_IMAGE077
，它被测加速度计输入轴、摆轴和输出轴上的分量分别为
Figure 89955DEST_PATH_IMAGE078
根据式(10)可得：
Figure 82182DEST_PATH_IMAGE079
（3）地球自转产生的Coriolis加速度分量
在加速度计原点处由地球自转角速率产生的Coriolis加速度很小，由离心机位姿误差引起的计算误差要小得多，可以忽略不计，因此考虑Coriolis加速度的标称值即可。此时，Coriolis加速度表达式为：
Figure 777605DEST_PATH_IMAGE080
其中
Figure 222493DEST_PATH_IMAGE081
为当地纬度。
综上可得加速度计三个轴上的精密比力为
Figure 966458DEST_PATH_IMAGE082
因为回转误差项对于比力的影响呈正弦、余弦形式的变化，采用整周积分时可以忽略，因含
Figure 78771DEST_PATH_IMAGE083
和
Figure 679516DEST_PATH_IMAGE084
的整周积分为零，亦可忽略，经计算得
Figure 877280DEST_PATH_IMAGE085
Figure 424936DEST_PATH_IMAGE086
Figure 126175DEST_PATH_IMAGE087
计算了加速度计的精确的比力输入，下面将用12位置法标定加速度计，可通过3种安装方式，利用公式（16）计算比力输入，再设计相应的试验方法。
具体地，本实施例所涉及的一种加速度计在精密离心机上的参数标定方法，石英加速度计高阶误差系数的具体计算过程为：
石英加速度计误差模型表达式采用如下形式：
Figure 599620DEST_PATH_IMAGE088
其中，
Figure 81417DEST_PATH_IMAGE089
为加速度计输出值，单位：V；
Figure 432764DEST_PATH_IMAGE090
为加速度计的输出当量，单位：g；
Figure 254089DEST_PATH_IMAGE091
为标度因数，单位：V/g；
Figure 399900DEST_PATH_IMAGE092
分别为加速度计输入轴，摆轴和输出轴上的加速度分量，单位：g；
Figure 306676DEST_PATH_IMAGE093
为零偏，单位：g；
Figure 992872DEST_PATH_IMAGE094
为交叉轴敏感度，单位：rad；
Figure 934283DEST_PATH_IMAGE095
为二阶非线性系数，单位：
Figure 250995DEST_PATH_IMAGE096
；
Figure 379488DEST_PATH_IMAGE097
为奇异二次项系数，单位：
Figure 574102DEST_PATH_IMAGE098
；
Figure 166758DEST_PATH_IMAGE099
为三阶非线性系数，单位：
Figure 919950DEST_PATH_IMAGE100
；
Figure 801318DEST_PATH_IMAGE101
为交叉耦合系数，单位：
Figure 767000DEST_PATH_IMAGE102
；
Figure 417425DEST_PATH_IMAGE103
随机误差，单位：g。
本发明主要针对石英加速度计高阶误差模型系数的测试与标定方法，因此，误差模型系数中的常值项和一次项当作已知量。本发明将采用6个对称位置组合来标定石英加速度计误差模型表达式中的高阶项误差系数。
通过图3所示6个对称位置组合来辨识石英加速度计的高阶误差模型系数，其中
Figure 138256DEST_PATH_IMAGE104
表示向心加速度矢量。图中所示的各个安装位置所对应的能标定的加速度计误差模型系数如表1所示。
表1对称位置组合与石英加速度计可辨识高阶误差模型系数的关系
Figure 506920DEST_PATH_IMAGE105
图3中总共采用3种安装方式，成对位置1-2，3-4，7-8是第1种安装方式，此时加速度计的输出轴始终与离心机的方位轴轴线一致，离心机的水平轴始终处于
Figure 541873DEST_PATH_IMAGE106
位置，方位轴处于如表1所示的6个位置，可获得3对成对位置。5-6，9-10位置为第2种安装方式，此时加速度计的输入轴始终与离心机的方位轴轴线一致，水平轴处于
Figure 781224DEST_PATH_IMAGE107
位置，方位轴处于4个位置可获得2对成对位置。11-12位置为第3种安装方式，此时加速度计的摆轴与离心机的方位轴轴线方向相反，水平轴始终处于
Figure 374754DEST_PATH_IMAGE108
位置，方位轴处于
Figure 27452DEST_PATH_IMAGE109
这2个位置。
根据式（16），可以得出第1~12安装位置对应的实际加速度计各轴的比力输入，具体计算时，输入轴上的比力精确到一阶小量，摆轴与输出轴上的比力只计算标称值，一阶小量也忽略，因为与这两个轴的输入比力相关的系数也是小量。公式（16）中的
Figure 334937DEST_PATH_IMAGE110
是已知量，用于计算加速度计的指示输出
Figure 959953DEST_PATH_IMAGE111
取到一阶小量，与其它系数相关的
Figure 225850DEST_PATH_IMAGE112
取标称值即可。要标定到加速度计的3阶误差模型系数，至少每个成对位置需要4个比力输入，即要求主轴运行于4个不同角速率
Figure 100265DEST_PATH_IMAGE113
，并采集加速度计的输出的整周均值。为方便起见，12位置采用统一的结构矩阵如公式（16）所示，当然也可以增加更多的角速率点进行测试。
Figure 742599DEST_PATH_IMAGE114
位置1上石英加速度计输入轴、摆轴和输出轴的比力分别为：
Figure 222121DEST_PATH_IMAGE115
其中
Figure 393340DEST_PATH_IMAGE116
均以g为单位，以下表达式相同。
将式(18)代入到式（16），石英加速度计在位置1的指示输出为：
Figure 958313DEST_PATH_IMAGE117
位置2上石英加速度计输入轴、摆轴和输出轴的比力分别为：
Figure 669917DEST_PATH_IMAGE118
将式(20)代入式（16），石英加速度计在位置2的指示输出为：
Figure 770991DEST_PATH_IMAGE119
分别将式(19)和式(21)相加和相减得到：
Figure 113110DEST_PATH_IMAGE120
对于式（22）为加速度的常数项，一次项和二次项组成。对于式（23）为常数项、一次项、二次项和三次项组成。综合以上分析，采用主轴4个速率点进行测试，可以辨识
Figure 165380DEST_PATH_IMAGE121
。
Figure 680675DEST_PATH_IMAGE122
其中“
Figure 869211DEST_PATH_IMAGE123
”表示此项理论上为零或者因为是很多位姿误差项的合成，而又没有必要写出。
式(24)写成矩阵形式为
Figure 647811DEST_PATH_IMAGE124
根据最小二乘可得：
Figure 984114DEST_PATH_IMAGE125
在式(24)中，辨识
Figure 506363DEST_PATH_IMAGE126
项规避了离心机误差
Figure 814984DEST_PATH_IMAGE127
，从而提高了
Figure 498906DEST_PATH_IMAGE128
项的标定精度。
根据式(23)可得
Figure 56927DEST_PATH_IMAGE129
其中
Figure 881401DEST_PATH_IMAGE130
Figure 778950DEST_PATH_IMAGE131
根据最小二乘可得：
Figure 899353DEST_PATH_IMAGE132
在观测向量
Figure 210248DEST_PATH_IMAGE133
中补偿了动态误差项
Figure 74299DEST_PATH_IMAGE134
，
Figure 91934DEST_PATH_IMAGE135
和Coriolis加速度项，在误差系数向量中加入了离心机的位姿误差项
Figure 383238DEST_PATH_IMAGE136
，自动补偿了静态半径测试误差
Figure 915850DEST_PATH_IMAGE137
以及回转误差项等，消除了离心机误差和Coriolis加速度的影响，从而提高了
Figure 318013DEST_PATH_IMAGE138
和
Figure 190154DEST_PATH_IMAGE139
项的标定精度。
位置3和4上石英加速度计输入轴、摆轴和输出轴的比力分别为：
Figure 419403DEST_PATH_IMAGE140
将式（28）和式（29）分别代入到式（16）中，计算出石英加速度计的指示输出
Figure 173733DEST_PATH_IMAGE141
和
Figure 379586DEST_PATH_IMAGE142
，并进行加减消元运算，得到以下表达式：
Figure 371813DEST_PATH_IMAGE143
Figure 4919DEST_PATH_IMAGE144
根据式（30）可得：
Figure 512124DEST_PATH_IMAGE145
其中
Figure 521668DEST_PATH_IMAGE146
根据式（31）可得：
Figure 102822DEST_PATH_IMAGE147
其中
Figure 172410DEST_PATH_IMAGE148
补偿掉动态失准角所产生的附加加速度后，可以辨识出
Figure 901331DEST_PATH_IMAGE149
和
Figure 714566DEST_PATH_IMAGE150
项。
位置5和6上石英加速度计输入轴、摆轴和输出轴的比力分别为：
Figure 914341DEST_PATH_IMAGE151
将式（34）和式（35）分别代入到式（16）中，计算出石英加速度计的指示输出
Figure 889251DEST_PATH_IMAGE152
和
Figure 371048DEST_PATH_IMAGE153
，并进行加减消元运算，得到以下表达式：
Figure 722394DEST_PATH_IMAGE154
根据式（36）可得：
Figure 278141DEST_PATH_IMAGE155
其中
Figure 423951DEST_PATH_IMAGE156
根据式(37)可得：
Figure 127465DEST_PATH_IMAGE157
其中
Figure 16924DEST_PATH_IMAGE158
同样补偿掉动态失准角所产生的附加加速度后，可以辨识出
Figure 692756DEST_PATH_IMAGE159
和
Figure 71784DEST_PATH_IMAGE160
项。
位置7和8上石英加速度计输入轴、摆轴和输出轴的比力分别为：
Figure 465857DEST_PATH_IMAGE161
Figure 406610DEST_PATH_IMAGE162
将式（39）和式（40）分别代入到式（16）中，计算出石英加速度计的指示输出
Figure 202528DEST_PATH_IMAGE163
，并进行加减消元运算，得到以下表达式：
Figure 955721DEST_PATH_IMAGE164
根据式（42）可得：
Figure 368247DEST_PATH_IMAGE165
其中
Figure 599509DEST_PATH_IMAGE166
精确辨识出
Figure 249933DEST_PATH_IMAGE167
系数后，再减去以前辨识的
Figure 970764DEST_PATH_IMAGE168
，可以辨识
Figure 339428DEST_PATH_IMAGE169
误差模型系数。
根据式（43）可得：
Figure 374381DEST_PATH_IMAGE170
其中
Figure 613732DEST_PATH_IMAGE171
位置9和10上石英加速度计输入轴、摆轴和输出轴的比力分别为：
Figure 207262DEST_PATH_IMAGE172
将式（45）和式（46）分别代入到式（16）中，计算出石英加速度计的指示输出
Figure 594381DEST_PATH_IMAGE173
，并进行加减消元运算，得到以下表达式：
Figure 167445DEST_PATH_IMAGE174
Figure 792461DEST_PATH_IMAGE175
根据式（48）可得：
Figure 58358DEST_PATH_IMAGE176
其中
Figure 932773DEST_PATH_IMAGE177
补偿动态失准角的影响之后，辨识
Figure 309527DEST_PATH_IMAGE178
后，再减去
Figure 523471DEST_PATH_IMAGE179
即可得
Figure 757006DEST_PATH_IMAGE180
。
根据式（49）可得：
Figure 321980DEST_PATH_IMAGE181
其中
Figure 236846DEST_PATH_IMAGE182
位置11和12上石英加速度计输入轴、摆轴和输出轴的比力分别为：
Figure 337920DEST_PATH_IMAGE183
将式(51)和式(52)分别代入到式(16)中，计算出石英加速度计的指示输出
Figure 742356DEST_PATH_IMAGE184
和
Figure 529047DEST_PATH_IMAGE185
，并进行加减消元运算，得到以下表达式：
Figure 247604DEST_PATH_IMAGE186
Figure 701719DEST_PATH_IMAGE187
根据式（54）可得：
Figure 480319DEST_PATH_IMAGE188
其中
Figure 551043DEST_PATH_IMAGE189
辨识出
Figure 807712DEST_PATH_IMAGE190
后，再减去辨识出的
Figure 116334DEST_PATH_IMAGE191
即可得到
Figure 862573DEST_PATH_IMAGE192
项。
根据式（55）可得：
Figure 889435DEST_PATH_IMAGE193
其中
Figure 182751DEST_PATH_IMAGE194
综合前面推出的公式，可得石英加速度计高阶误差项的标定结果如下：
Figure 611458DEST_PATH_IMAGE195
可归纳出石英加速度计高阶误差模型系数的表达式为
Figure 528598DEST_PATH_IMAGE196
如图2所示，本实施例所涉及的一种加速度计在精密离心机上的参数标定方法，令
Figure 777177DEST_PATH_IMAGE197
，可得误差模型系数
Figure 641228DEST_PATH_IMAGE198
项的表达式为：
Figure 658863DEST_PATH_IMAGE199
其中
Figure 950167DEST_PATH_IMAGE200
表示矩阵
Figure 482779DEST_PATH_IMAGE201
第
Figure 884942DEST_PATH_IMAGE202
行，第
Figure 22662DEST_PATH_IMAGE203
列的元素。假设石英加速度计的指示输出独立且精度相等，其不确定度为
Figure 281605DEST_PATH_IMAGE204
，那么
Figure 6241DEST_PATH_IMAGE205
项的不确定度为
Figure 212094DEST_PATH_IMAGE206
假设离心机提供5g，10g，15g和20g的向心加速度，石英加速度计的输出的不确定度为
Figure 204321DEST_PATH_IMAGE207
，动态失准角不确定度
Figure 634165DEST_PATH_IMAGE208
，动态半径误差的不确定度
Figure 344632DEST_PATH_IMAGE209
。经计算可得，石英加速度计的二次项和交叉二次项的不确定度分别为
Figure 88597DEST_PATH_IMAGE210
Figure 935331DEST_PATH_IMAGE212
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以权利要求的保护范围为准。

Claims (6)

1.一种加速度计在精密离心机上的参数标定方法，其特征在于，包括：
获取精密离心机的各静态误差以及动态误差，并根据精密离心机的结构建立坐标系，以及根据所述各静态误差以及动态误差计算所述坐标系下的位姿误差；
驱动精密离心机的主轴以匀角速率旋转，以产生向心加速度标定加速度计，基于所述坐标系下的位姿误差计算向心加速度、重力加速度和Coriolis加速度的比力分配，以确定加速度计误差模型；
对加速度计在三种不同安装方式下的六个对称位置的指示输出，利用加减消元的方法标定加速度计误差模型表达式中的高阶项误差系数。
2.根据权利要求1所述的加速度计在精密离心机上的参数标定方法，其特征在于，所述精密离心机包括主轴、水平轴和方位轴；
所述精密离心机的静态误差包括主轴轴线的二维铅垂度误差
Figure 619406DEST_PATH_IMAGE001
水平轴轴线与主轴轴线的垂直度
Figure 940666DEST_PATH_IMAGE002
相交度
Figure 429416DEST_PATH_IMAGE003
水平轴轴线与方位轴轴线的垂直度
Figure 759904DEST_PATH_IMAGE004
相交度
Figure 752130DEST_PATH_IMAGE005
以及方位轴的初始零位误差
Figure 916396DEST_PATH_IMAGE006
安装惯性仪表的工作基面对方位轴轴线的垂直度
Figure 751496DEST_PATH_IMAGE007
加速度计安装基面姿态误差
Figure 761041DEST_PATH_IMAGE008
偏心误差
Figure 997987DEST_PATH_IMAGE009
以及初始对零误差
Figure 801995DEST_PATH_IMAGE010
主轴、水平轴和方位轴三个轴的角位置误差分别为
Figure 265337DEST_PATH_IMAGE011
所述精密离心机的动态误差包括主轴径向回转误差
Figure 937627DEST_PATH_IMAGE012
轴向窜动
Figure 904446DEST_PATH_IMAGE013
及倾角回转误差
Figure 535148DEST_PATH_IMAGE014
动态半径误差
Figure 954628DEST_PATH_IMAGE015
动态失准角
Figure 571554DEST_PATH_IMAGE016
水平轴径向回转误差
Figure 517513DEST_PATH_IMAGE017
轴向窜动
Figure 928903DEST_PATH_IMAGE018
以及倾角回转误差
Figure 225892DEST_PATH_IMAGE019
方位轴径向回转误差
Figure 115350DEST_PATH_IMAGE020
，轴向窜动
Figure 181395DEST_PATH_IMAGE021
倾角回转误差
Figure 29266DEST_PATH_IMAGE022
；
半径误差
Figure 423338DEST_PATH_IMAGE023
其中，
Figure 241121DEST_PATH_IMAGE024
为静态半径标称值，是由计量部门标定出的已知量，半径的静态测试误差
Figure 37039DEST_PATH_IMAGE025
是未知量，
Figure 914865DEST_PATH_IMAGE026
为用双频激光干涉仪监测的精密离心机在运行状态下的实际工作半径相对于离心机静态半径的变化量，是主轴角速率
Figure 61813DEST_PATH_IMAGE027
的函数。
3.根据权利要求2所述的加速度计在精密离心机上的参数标定方法，其特征在于，所述根据精密离心机的结构建立坐标系，以及根据所述各静态误差以及动态误差计算所述坐标系下的位姿误差包括：
建立地理坐标系
Figure 558653DEST_PATH_IMAGE028
轴水平指东，
Figure 333711DEST_PATH_IMAGE029
轴水平指北，
Figure 257805DEST_PATH_IMAGE030
轴指天，构成右手坐标系；
建立主轴轴套坐标系
Figure 751103DEST_PATH_IMAGE031
得到主轴轴套坐标系相对于地理坐标系的位姿；
建立主轴坐标系
Figure 786055DEST_PATH_IMAGE032
得到主轴坐标系相对于主轴轴套坐标系的位姿；
建立水平轴轴套坐标系
Figure 822144DEST_PATH_IMAGE033
得到水平轴轴套坐标系相对于主轴坐标系的位姿；
建立水平轴坐标系
Figure 776194DEST_PATH_IMAGE034
得到水平轴坐标系相对于水平轴轴套坐标系的位姿；
建立方位轴轴套坐标系
Figure 632154DEST_PATH_IMAGE035
得到方位轴轴套坐标系相对于水平轴坐标系的位姿；
建立方位轴坐标系
Figure 861010DEST_PATH_IMAGE036
得到方位轴坐标系相对于方位轴轴套坐标系的位姿；
建立工作基面坐标系
Figure 220447DEST_PATH_IMAGE037
得到工作基面坐标系相对于方位轴坐标系的位姿；
建立加速度计坐标系
Figure 610978DEST_PATH_IMAGE038
得到加速度计坐标系相对于工作基面坐标系的位姿、加速度计坐标系相对于地理坐标系的位姿、加速度计坐标系相对于主轴坐标系的位姿。
4.根据权利要求3所述的加速度计在精密离心机上的参数标定方法，其特征在于，
主轴轴套坐标系相对于地理坐标系的位姿为
Figure 954234DEST_PATH_IMAGE039
主轴坐标系相对于主轴轴套坐标系的位姿为
Figure 596568DEST_PATH_IMAGE040
其中
Figure 466304DEST_PATH_IMAGE041
表示主轴旋转的角度；
水平轴轴套坐标系相对于主轴坐标系的位姿为
Figure 903102DEST_PATH_IMAGE042
水平轴坐标系相对于水平轴轴套坐标系的位姿为
Figure 616146DEST_PATH_IMAGE043
其中
Figure 62171DEST_PATH_IMAGE044
表示水平轴旋转的角度；
方位轴轴套坐标系相对于水平轴坐标系的位姿为
Figure 396201DEST_PATH_IMAGE045
方位轴坐标系相对于方位轴轴套坐标系的位姿为
Figure 394112DEST_PATH_IMAGE046
其中
Figure 446382DEST_PATH_IMAGE047
表示方位轴旋转的角度；
工作基面坐标系相对于方位轴坐标系的位姿为
Figure 430519DEST_PATH_IMAGE048
其中
Figure 9268DEST_PATH_IMAGE049
为
Figure 53447DEST_PATH_IMAGE050
点相对
Figure 717646DEST_PATH_IMAGE051
点位移；
加速度计坐标系相对于工作基面坐标系的位姿为
Figure 239895DEST_PATH_IMAGE052
其中
Figure 814095DEST_PATH_IMAGE053
为
Figure 888231DEST_PATH_IMAGE054
点相对
Figure 915093DEST_PATH_IMAGE055
点位移；
加速度计坐标系相对于地理坐标系的位姿为
Figure 365665DEST_PATH_IMAGE056
其中
Figure 528794DEST_PATH_IMAGE057
表示加速度计坐标系与地理坐标系之间的姿态变换矩阵，
Figure 914775DEST_PATH_IMAGE058
为加速度计坐标系与地理坐标系的相对位移矢量；
加速度计坐标系相对于主轴坐标系的位姿为
Figure 553567DEST_PATH_IMAGE059
其中
Figure 417618DEST_PATH_IMAGE060
表示加速度计坐标系与主轴坐标系之间的姿态变换矩阵；
加速度计坐标系原点在主轴坐标系下表示为
Figure 825466DEST_PATH_IMAGE061
忽略二阶小量后可得，
Figure 382349DEST_PATH_IMAGE062
Figure 383803DEST_PATH_IMAGE063
5.根据权利要求1所述的加速度计在精密离心机上的参数标定方法，其特征在于，所述加速度计误差模型的表达式为：
Figure 176178DEST_PATH_IMAGE064
其中，
Figure 48320DEST_PATH_IMAGE065
为加速度计输出值；
Figure 900738DEST_PATH_IMAGE066
为加速度计的输出当量；
Figure 389488DEST_PATH_IMAGE067
为标度因数；
Figure 719975DEST_PATH_IMAGE068
分别为加速度计输入轴，摆轴和输出轴上的加速度分量；
Figure 977781DEST_PATH_IMAGE069
为零偏；
Figure 876467DEST_PATH_IMAGE070
为交叉轴敏感度；
Figure 445989DEST_PATH_IMAGE071
为二阶非线性系数；
Figure 721112DEST_PATH_IMAGE072
为奇异二次项系数；
Figure 958059DEST_PATH_IMAGE073
为三阶非线性系数；
Figure 762066DEST_PATH_IMAGE074
为交叉耦合系数；
Figure 225409DEST_PATH_IMAGE075
—随机误差；
所述高阶项误差系数包括所述二阶非线性系数、奇异二次项系数、三阶非线性系数以及交叉耦合系数。
6.根据权利要求5所述的加速度计在精密离心机上的参数标定方法，其特征在于，所述对加速度计在三种不同安装方式下的六个对称位置的指示输出，利用加减消元的方法标定加速度计误差模型表达式中的高阶项误差系数包括：
当加速度计的输出轴始终与离心机的方位轴轴线一致，离心机的水平轴始终处于
Figure 897699DEST_PATH_IMAGE076
位置，可获得3对成对位置，分别为：位置1与位置2、位置3与位置4、位置5与位置6；当加速度计的输入轴始终与离心机的方位轴轴线一致，水平轴处于
Figure 864518DEST_PATH_IMAGE077
位置，可获得2对成对位置，分别为：位置7与位置8、位置9与位置10；当加速度计的摆轴与离心机的方位轴轴线方向相反，水平轴始终处于
Figure 229640DEST_PATH_IMAGE078
位置，方位轴处于
Figure 180278DEST_PATH_IMAGE079
2个位置，分别为：位置11与位置12；
12个位置采用统一的结构矩阵如式(17)所示，
Figure 531625DEST_PATH_IMAGE080
根据位置1和位置2上加速度计输入轴、摆轴和输出轴的比力辨识出
Figure 477585DEST_PATH_IMAGE081
项和
Figure 154554DEST_PATH_IMAGE082
项；
根据位置3和位置4上加速度计输入轴、摆轴和输出轴的比力辨识出
Figure 185963DEST_PATH_IMAGE083
和
Figure 341001DEST_PATH_IMAGE084
项；
根据位置5和位置6上加速度计输入轴、摆轴和输出轴的比力辨识出
Figure 282412DEST_PATH_IMAGE085
项和
Figure 723758DEST_PATH_IMAGE086
项；
根据位置7和位置8上加速度计输入轴、摆轴和输出轴的比力辨识出
Figure 117830DEST_PATH_IMAGE087
项；
根据位置9和位置10上加速度计输入轴、摆轴和输出轴的比力辨识出
Figure 201193DEST_PATH_IMAGE088
项；
根据位置11和位置12上加速度计输入轴、摆轴和输出轴的比力辨识出
Figure 997111DEST_PATH_IMAGE089
项。
CN202110313204.0A 2021-03-24 2021-03-24 加速度计在精密离心机上的参数标定方法 Active CN112698055B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110313204.0A CN112698055B (zh) 2021-03-24 2021-03-24 加速度计在精密离心机上的参数标定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110313204.0A CN112698055B (zh) 2021-03-24 2021-03-24 加速度计在精密离心机上的参数标定方法

Publications (2)

Publication Number Publication Date
CN112698055A CN112698055A (zh) 2021-04-23
CN112698055B true CN112698055B (zh) 2021-06-25

Family

ID=75515561

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110313204.0A Active CN112698055B (zh) 2021-03-24 2021-03-24 加速度计在精密离心机上的参数标定方法

Country Status (1)

Country Link
CN (1) CN112698055B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113865583B (zh) * 2021-07-20 2024-02-09 北京航天控制仪器研究所 一种加速度计组合动态安装偏差矩阵确定及补偿方法
CN113865585B (zh) * 2021-09-07 2023-08-29 北京航天控制仪器研究所 一种陀螺仪组合高阶误差系数分离与补偿方法和***
CN113804221B (zh) * 2021-10-14 2023-10-03 天津科技大学 一种基于模观测法的离心机加速度计组合标定方法
CN114034885B (zh) * 2021-11-11 2024-04-26 哈尔滨工业大学 一种基于全误差分析的陀螺加速度计在双轴离心机上的测试方法
CN114018292A (zh) * 2021-11-11 2022-02-08 中国航空工业集团公司北京长城计量测试技术研究所 双轴离心机控制方法、装置、计算机设备和存储介质
CN118032013B (zh) * 2024-04-11 2024-06-18 伸瑞科技（北京）有限公司 正交双加速度计在分度头上的标定精度验证方法和***

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102221372B (zh) * 2011-03-25 2012-10-10 北京航空航天大学 使用离心机和转台对惯性测量单元进行误差标定的方法
US9207256B2 (en) * 2012-02-14 2015-12-08 Infineon Technologies Ag Auto-calibration of acceleration sensors
CN102841218A (zh) * 2012-08-21 2012-12-26 哈尔滨工业大学 一种基于双轴离心机的陀螺加速度计测试方法
CN108917787B (zh) * 2018-04-20 2021-04-13 北京航天控制仪器研究所 一种mems陀螺仪标度因数的加速度敏感性补偿方法
CN109813343B (zh) * 2019-03-21 2021-06-08 哈尔滨工业大学 一种离心机初始对准误差的测量方法
CN111879335A (zh) * 2019-09-20 2020-11-03 天津科技大学 一种基于离心机的多位置陀螺仪漂移系数的标定方法
CN111781400B (zh) * 2020-07-10 2021-08-10 哈尔滨工业大学 一种标定加速度计高阶误差系数的方法
CN111947683B (zh) * 2020-07-17 2022-07-05 北京航天控制仪器研究所 一种精密离心机半径误差离线测量与在线补偿方法和装置

Also Published As

Publication number Publication date
CN112698055A (zh) 2021-04-23

Similar Documents

Publication Publication Date Title
CN112698055B (zh) 加速度计在精密离心机上的参数标定方法
CN106017507B (zh) 一种用于中低精度的光纤惯组快速标定方法
CN110006450B (zh) 一种激光捷联惯导***在卧式三轴转台上的标定方法
CN101290326B (zh) 石英挠性加速度计测量组件的参数辨识标定方法
CN107655493B (zh) 一种光纤陀螺sins六位置***级标定方法
CN108592952A (zh) 基于杆臂补偿与正反倍速率同时标定多mimu误差的方法
CN113156166A (zh) 石英加速度计在精密离心机上的对称融消测试方法
CN109974749B (zh) 一种三轴转台综合指向误差的评定方法
CN107677292B (zh) 基于重力场模型的垂线偏差补偿方法
Ren et al. Calibration method of accelerometer’s high-order error model coefficients on precision centrifuge
CN101629830A (zh) 一种三轴一体高精度光纤陀螺的标定方法及装置
Wang et al. An improve hybrid calibration scheme for strapdown inertial navigation system
CN113701747B (zh) 一种基于离心机激励的惯性测量***姿态角误差分离方法
Ermakov et al. Angular velocity estimation of rotary table bench using aggregate information from the sensors of different physical nature
Sun et al. Sequential calibration method of nonlinear errors of PIGA on counter-rotating platform centrifuge
Yingbo et al. Calibration method of quartz accelerometer on dynamic centrifuge
CN115979311B (zh) Piga交叉二次项系数标定方法、***、设备及介质
RU2717566C1 (ru) Способ определения погрешностей инерциального блока чувствительных элементов на двухосном поворотном столе
CN113945230B (zh) 一种惯性器件的高阶误差系数的辨识方法
CN115876225A (zh) 基于二自由度转台的mems imu标定方法及***
CN109871658A (zh) 用于导弹弹头转动惯量和惯性积测量的多姿态最优估计方法
CN114034885B (zh) 一种基于全误差分析的陀螺加速度计在双轴离心机上的测试方法
CN102707092A (zh) 一种基于角速率台的单波束激光测速仪标定方法
CN108716925A (zh) 一种九轴传感器的标定方法及装置
CN113865585B (zh) 一种陀螺仪组合高阶误差系数分离与补偿方法和***

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant