CN109085654B - 一种旋转加速度计重力梯度仪数字建模仿真方法 - Google Patents
一种旋转加速度计重力梯度仪数字建模仿真方法 Download PDFInfo
- Publication number
- CN109085654B CN109085654B CN201810599648.3A CN201810599648A CN109085654B CN 109085654 B CN109085654 B CN 109085654B CN 201810599648 A CN201810599648 A CN 201810599648A CN 109085654 B CN109085654 B CN 109085654B
- Authority
- CN
- China
- Prior art keywords
- accelerometer
- coordinate system
- vector
- axis
- time
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V7/00—Measuring gravitational fields or waves; Gravimetric prospecting or detecting
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Force Measurement Appropriate To Specific Purposes (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开一种旋转加速度计重力梯度仪数字建模仿真方法,基于重力梯度仪的测量原理与结构,综合考虑重力梯度仪内部加速度计的径向安装误差、高度安装误差、输入轴失准误差、电路增益不匹配误差、线性标度系数不匹配误差、零偏不匹配误差、加速度计二阶误差系数、噪声、环境物体自梯度等因素,能够模拟真实旋转加速度计重力梯度仪。它能解决重力梯度仪研发中,一些关键技术(比如自梯度标定、补偿,加速计二阶误差系数在线补偿,加速度计安装误差在线补偿,加速度计线性标度系数在线补偿,故障诊断等)的实验验证问题。
Description
技术领域
本发明涉及一种旋转加速度计重力梯度仪数字建模仿真方法,属于精密测量技术领域。
背景技术
重力梯度仪是一种高精度万有引力梯度检测仪,能够测量仪器所在位置点的万有引力 梯度值。航空重力梯度勘探能高精度、高效率的收集重力梯度信息,测量的重力梯度数据 可直接用于分析勘探区域的地质密度、地质结构、重建重力场等。目前国外已经投入商业 应用的重力梯度仪主要为旋转加速度计重力梯度仪,在研的具有应用前景的重力梯度仪有 冷原子重力梯度仪、超导重力梯度仪、MEMS重力梯度仪等,我国的重力梯度仪样机正处 于研制中。
目前旋转加速度计重力梯度仪信号仿真方法,过于理想化,没有综合考虑重力梯度仪 内部加速度计的安装误差、加速度计二阶非线性误差系数、加速度计电流转电压增益不匹 配、载体运动等因素,不能模拟真实重力梯度仪。
发明内容
技术问题:本发明提供一种能模拟真实的重力梯度仪,各项参数、指标可控,解决了 重力梯度仪研发中各项技术的实验验证问题的旋转加速度计重力梯度仪数字建模仿真方 法。
技术方案:本发明的旋转加速度计重力梯度仪数字建模仿真方法,包括以下步骤:
1)设置重力梯度仪的参数:
设置重力梯度仪加速度计安装圆盘的旋转角频率Ω;
设置四只加速度计的安装参数及输出模型参数:设置加速度计A1的径向安装距离R1, 加速度计A1的初相位角β1z,加速度计A1的高度角β1x,加速度计A1的输入轴失准角 加速度计A1的零偏k10,加速度计A1的线性标度系数k11,加速度计A1的二阶误差 系数k12,k14,k15,k16,k17,k18,加速度计A1的电流转电压增益k1V/I;设置加速度计 A2的径向安装距离R2,加速度计A2的初相位角β2z,加速度计A2的高度角β2x,加速度 计A2的输入轴失准角加速度计A2的零偏k20,加速度计A2的线性标度系数k21, 加速度计A2的二阶误差系数k22,k24,k25,k26,k27,k28,加速度计A2的电流转电压 增益k2V/I;设置加速度计A3的径向安装距离R3,加速度计A3的初相位角β3z,加速度计 A3的高度角β3x,加速度计A3的输入轴失准角加速度计A3的零偏k30,加速 度计A3的线性标度系数k31,加速度计A3的二阶误差系数k32,k34,k35,k36,k37,k38, 加速度计A3的电流转电压增益k3V/I;设置加速度计A4的径向安装距离R4,加速度计A4 的初相位角β4z,加速度计A4的高度角β4x,加速度计A4的输入轴失准角加 速度计A4的零偏k40,加速度计A4的线性标度系数k41,加速度计A4的二阶误差系数k42, k44,k45,k46,k47,k48,加速度计A4的电流转电压增益k4V/I;
设置加速度计的噪声系数α,b,ωT;
设置载体加速度向量fcmm,角速度向量ω,角加速度向量ωa;
设置采样频率fs,仿真总时间Ttotal;
2)根据步骤1)中加速度计的安装参数,计算t时刻加速度计A1在重力梯度仪测量坐标系中的位置向量计算t时刻加速度计A2在重力梯度仪测量坐标系中的位置向量计算t时刻加速度计A3在重力梯度仪测量坐标系中的位置向量计算t时刻加速 度计A4在重力梯度仪测量坐标系中的位置向量
根据步骤1)中加速度计的安装参数,计算t时刻加速度计A1的标称安装坐标系的x轴单位方向向量τ1x,y轴单位方向向量τ1y,z轴单位方向向量τ1z;计算t时刻加速度计 A2的标称安装坐标系的x轴单位方向向量τ2x,y轴单位方向向量τ2y,z轴单位方向向量 τ2z;计算t时刻加速度计A3的标称安装坐标系的x轴单位方向向量τ3x,y轴单位方向向 量τ3y,z轴单位方向向量τ3z;计算t时刻加速度计A4的标称安装坐标系的x轴单位方向 向量τ4x,y轴单位方向向量τ4y,z轴单位方向向量τ4z;
根据步骤1)中加速度计的安装参数,计算t时刻加速度计A1的标称坐标系到加速度 计A1的测量坐标系的变换矩阵C1;计算t时刻加速度计A2的标称坐标系到加速度计A2的测量坐标系的变换矩阵C2;计算t时刻加速度计A3的标称坐标系到加速度计A3的测 量坐标系的变换矩阵C3;计算t时刻加速度计A4的标称坐标系到加速度计A4的测量坐 标系的变换矩阵C4;
根据以下功率谱密度模型计算t时刻加速度计A1的噪声f1noise、加速度计A2的噪声f2noise、加速度计A3的噪声f3noise、加速度计A4的噪声f4noise:
Φ(f)noise=αf-b+ωT
4)根据下式分别计算t时刻加速度计A1的比力向量f1、t时刻加速度计A2的比力向量f2、t时刻加速度计A3的比力向量f3、t时刻加速度计A4的比力向量f4:
式中,×为向量叉乘;G为万有引力常数;
5)根据下式分别计算t时刻四只加速度计在标称安装坐标系的x方向的比力、y方向 的比力、z方向的比力:
式中f1x,f1y,f1z分别是加速度计A1在其标称安装坐标系的x方向,y方向,z方向的比 力;f2x,f2y,f2z分别是加速度计A2在其标称安装坐标系的x方向,y方向,z方向的比力;f3x,f3y,f3z分别是加速度计A3在其标称安装坐标系的x方向,y方向,z方向的比力; f4x,f4y,f4z分别是加速度计A4在其标称安装坐标系的x方向,y方向,z方向的比力;式 中.为向量点积;
6)根据下式分别计算t时刻四只加速度计在加速度计测量坐标系的输入轴的比力,输 出轴的比力,摆轴的比力:
式中f1i,f1o,f1p分别是加速度计A1在其测量坐标系的输入轴、输出轴、摆轴的比力; f2i,f2o,f2p分别是加速度计A2在其测量坐标系的输入轴、输出轴、摆轴的比力,f3i,f3o,f3p分别是加速度计A3在其测量坐标系的输入轴、输出轴、摆轴的比力,f4i,f4o,f4p分别是加速度计A4在其测量坐标系的输入轴、输出轴、摆轴的比力;
7)根据下式分别计算t时刻加速度计A1的输出电压V1、加速度计A2的输出电压V2、加速度计A3的输出电压V3、加速度计A4的输出电压V4:
8)根据下式计算t时刻重力梯度仪解调前的信号Gout:
Gout=V1+V2-V3-V4
9)令t=t+1/fs,然后判断更新后的时间t是否小于仿真总时间Ttotal,如是,则返回步 骤2),否则进入步骤10);
10)对信号Gout抗混叠滤波、降采样、带通滤波,然后进行正交幅值解调,恢复重力梯度信号。
进一步,本发明方法中,步骤1)中,检测质量既可以是密度均匀、形状规则的物体,也可以是密度、形状不规则的物体;检测质量的运动状态,可以是静止,也可以是运动;
进一步,本发明方法中,步骤1)到5)中向量均为3维向量,矩阵为3×3方阵;
进一步,本发明方法中,步骤2)和步骤5)中的加速度计标称安装坐标系的原点在加速度计的质心,加速度计标称安装坐标系的x轴定义为该点切线的方向,y轴定义为旋 转圆盘的中心点到加速度计安装点的径向方向,z轴平行旋转圆盘轴线,方向垂直向下。
进一步,本发明方法中,步骤2)和步骤6)中的加速度计测量坐标系的原点在加速度计的质心,加速度计测量坐标系的x轴定义为加速度计输入轴方向,y轴定义为加速度 计输出轴的方向,z轴定义为加速度计摆轴的方向。
进一步,本发明方法中,步骤1)和步骤2)中的重力梯度仪测量坐标系的原点位于旋转圆盘中心,重力梯度仪测量坐标系的x轴定义为旋转圆盘中心点到加速度计A1的初 始位置,y轴定义为旋转圆盘的中心点到加速度计A3的初始位置,z轴定义为旋转圆盘轴 线方向。
在旋转加速度计重力梯度仪的研发中,一些关键技术比如自梯度的标定和补偿,误差 增益系数的标定,加速度计安装误差在线补偿,加速度计二阶误差系数在线补偿,加速度 计一阶标度系数在线补偿,故障诊断等,这些关键技术需要大量实验论证,验证实验需要 一台更高精度的旋转加速度计重力梯度仪作对照,来判断实验效果。然而,实际中无法获 得高精度的旋转加速度计重力梯度仪。为了解决该问题,本发明根据旋转加速度计重力梯 度仪的原理、结构,综合考虑了重力梯度仪内部加速度计的安装误差、二阶非线性误差系 数、噪声、电路增益等因素,建立了虚拟旋转加速度计重力梯度仪,它能模拟真实的重力 梯度仪,而且它的各项参数、指标可控,它解决了重力梯度仪研发中各项技术的实验验证 问题。
有益效果:本发明与现有技术相比,具有以下优点:
目前有关文献中的旋转加速度计重力梯度仪信号仿真方法,过于理想化,没有综合考 虑重力梯度仪内部加速度计的安装误差、加速度计二阶非线性误差系数、加速度计电流转 电压增益不匹配、载体运动等因素,不能模拟真实重力梯度仪。本发明根据旋转加速度计 重力梯度仪的原理、结构,综合考虑了重力梯度仪内部加速度计的安装误差、二阶非线性 误差系数、噪声、电路增益等因素,建立了虚拟旋转加速度计重力梯度仪,它能模拟真实 的重力梯度仪,而且它的各项参数、指标可控,它解决了重力梯度仪研发中各项技术的实 验验证问题。
附图说明
图1加速度计安装误差示意图。
图2旋转加速度计重力梯度仪数字模型原理图。
图3为时间长度为重力梯度仪旋转圆盘旋转一周的加速度计输出电压数据
图4为实验条件下重力梯度仪的输出。
具体实施方式
下面结合实施例和说明书附图对本发明作进一步地说明。
图1是旋转加速度计重力梯度仪加速度计安装误差示意图。图中坐标系xmymzm是重力 梯度仪测量坐标系,它的x轴指向加速度计A1的初始安装位置,y轴指向加速度计A3的初始安装位置,z轴和圆盘旋转轴重合。旋转加速度计重力梯度仪要求四只加速度计的质心处于同一个平面,加速度计的质心到圆盘中心的距离相等,且加速度计的输出轴要与加速度计安装位置点的旋转圆盘切线方向一致。图1中加速度计A1,A2,A3,A4均无安装 误差,这四只加速度计的质点处于同一平面,这个平面称为参考平面,四只加速度计到圆 盘中心的距离相等,该距离称为径向标称距离。加速度计是加速度计A1存在安装误差 时的情况,图1以加速度计A1为例描述加速度计安装误差。加速度计的安装误差包括安 装位置点误差,输入轴失准误差。加速度计的安装位置由3个参数确定:径向距离,高度 角,初始相位角。径向距离,加速度计的径向距离等于标称径向距离加上安装误差;加速 度计的高度角,描述加速度计的质心偏离参考平面的程度;加速度计的初相位角,是加速 度计的实际安装位置点在参考平面的投影到加速度计理论安装点的圆心角,它能描述加速 度计安装位置的偏离程度。图1中加速度计A1 0的径向距离为R1,高度角为β1x,初始相位 角β1z。输入轴失准角是加速度计输入轴偏离安装位置点圆盘切线方向,图1中A1 0的标称 安装坐标系为x1y1z1坐标系,它的原点在加速度计的质心,它的x1轴方向和圆盘切线方向一 致,它的y1轴从圆盘中心指向加速度计安装位置,z1轴同圆盘旋转轴线平行,方向向下; 图1中A1 0的加速度计测量坐标系为aiaoap;输入轴失准角是加速计输入轴ai偏离x1方向的 误差,可以用两个角度描述。因此,重力梯度仪旋转圆盘上的每只加速度计安装 误差可以用5个安装参数来描述。同时,每只加速度计的输出模型有9个参数,这些参数 分别属于零偏、线性标度系数、二阶非线性误差系数、电流转电压增益系数。
本发明的旋转加速度计重力梯度仪数字建模仿真方法,包括以下步骤:
1)设置重力梯度仪的参数:
设置重力梯度仪加速度计安装圆盘的旋转角频率Ω=1.571rad/s;
设置四只加速度计的安装参数及输出模型参数:四只加速度计的安装参数及输出模型 参数均列出在下表中:
表中Rj表示加速度计的径向安装距离,βjz表示加速度计的初相位角,βjx表示加速度 计的高度角,表示加速度计的输入轴失准角,kj0表示加速度计的零偏,kj1表示加速度计的线性标度系数,kj2,kj4,kj5,kj6,kj7,kj8表示加速度计的二阶误差系数,kjV/I表示加速度计的电流转电压增益;
设置N=1个检测质量的参数:检测质量的质量500kg,检测质量在重力梯度仪测量坐 标系中的位置向量设置加速度计的噪声系数b=1.7, 设置载体加速度向量fcmm=[0.001,0.004,0.002]g,角速度向量 ω=[100,300,600]arcsec/s,角加速度向量ωa=[0,0,0];设置采样频率fs=10Hz,仿真总时间 Ttotal=1600秒;
2)根据步骤1)中加速度计的安装参数,计算t时刻加速度计A1在重力梯度仪测量坐标系中的位置向量计算t时刻加速度计A2在重力梯度仪测量坐标系中的位置向量计算t时刻加速度计A3在重力梯度仪测量坐标系中的位置向量计算t时刻加速 度计A4在重力梯度仪测量坐标系中的位置向量
根据步骤1)中加速度计的安装参数,计算t时刻加速度计A1的标称安装坐标系的x轴单位方向向量τ1x,y轴单位方向向量τ1y,z轴单位方向向量τ1z;计算t时刻加速度计 A2的标称安装坐标系的x轴单位方向向量τ2x,y轴单位方向向量τ2y,z轴单位方向向量 τ2z;计算t时刻加速度计A3的标称安装坐标系的x轴单位方向向量τ3x,y轴单位方向向 量τ3y,z轴单位方向向量τ3z;计算t时刻加速度计A4的标称安装坐标系的x轴单位方向 向量τ4x,y轴单位方向向量τ4y,z轴单位方向向量τ4z;
根据步骤1)中加速度计的安装参数,计算t时刻加速度计A1的标称坐标系到加速度 计A1的测量坐标系的变换矩阵C1;计算t时刻加速度计A2的标称坐标系到加速度计A2的测量坐标系的变换矩阵C2;计算t时刻加速度计A3的标称坐标系到加速度计A3的测 量坐标系的变换矩阵C3;计算t时刻加速度计A4的标称坐标系到加速度计A4的测量坐 标系的变换矩阵C4;根据以下功率谱密度模型:
Φ(f)noise=αf-b+ωT
计算t时刻加速度计A1的噪声f1noise,加速度计A2的噪声f2noise,加速度计A3的噪声f3noise,加速度计A4的噪声f4noise,参数α,b,ωT是步骤1)中设置的加速度计噪声参数;
4)根据下式分别计算t时刻加速度计A1的比力向量f1,计算t时刻加速度计A2的比力向量f2、t时刻加速度计A3的比力向量f3、t时刻加速度计A4的比力向量f4:
式中,×为向量叉乘;G为万有引力常数
5)根据下式分别计算t时刻四只加速度计在标称安装坐标系的x方向的比力,y方向 的比力,z方向的比力:
式中f1x,f1y,f1z分别是加速度计A1在其标称安装坐标系的x方向,y方向,z方向的比 力;f2x,f2y,f2z分别是加速度计A2在其标称安装坐标系的x方向,y方向,z方向的比力;f3x,f3y,f3z分别是加速度计A3在其标称安装坐标系的x方向,y方向,z方向的比力; f4x,f4y,f4z分别是加速度计A4在其标称安装坐标系的x方向,y方向,z方向的比力;式 中[·]为向量点积;
6)计算t时刻四只加速度计在加速度计测量坐标系的输入轴的比力,输出轴的比力, 摆轴的比力:
式中f1i,f1o,f1p分别是加速度计A1在其测量坐标系的输入轴,输出轴,摆轴的比力; f2i,f2o,f2p分别是加速度计A2在其测量坐标系的输入轴,输出轴,摆轴的比力;f3i,f3o,f3p分别是加速度计A3在其测量坐标系的输入轴,输出轴,摆轴的比力;f4i,f4o,f4p分别是加速度计A4在其测量坐标系的输入轴,输出轴,摆轴的比力;
7)计算t时刻加速度计A1的输出电压V1,加速度计A2的输出电压V2,加速度计A3 的输出电压V3,加速度计A4的输出电压V4:
8)计算t时刻重力梯度仪解调前的信号Gout:
Gout=V1+V2-V3-V4
9)令t=t+1/fs,然后判断更新后的时间t是否小于仿真总时间Ttotal,如是,则返回步 骤2),否则进入步骤10);
步骤1)设置了重力梯度仪的参数、检测质量的参数,其中重力梯度仪的参数中,加速度计的安装误差、加速度计模型参数,决定了载体加速度、角速度、角加速度的误差传 递特性;通过步骤2)到步骤9),可以得到该实验条件下,重力梯度仪四只加速度计的输 出电压数据,为了清晰的显示该数据,选取了时间长度为重力梯度仪旋转圆盘旋转一周的 加速度计输出数据,如图3所示。
10)如图2所示,对信号Gout抗混叠滤波、降采样、带通滤波,然后进行正交幅值解调,恢复重力梯度信号。图4是该实验条件下,重力梯度仪的输出。从图4可以看出,inline通道的输出高达1929Eo,cross通道高达-744Eo。该实验中,检测质量的质量为500Kg, 坐标为(0.3,0.5,0.6),理论上在该检测质量的激励下,重力梯度仪的Inline通道的输出应为39.0713Eo,cross通道的输出应为-36.6293Eo。由于本实验中,设置了加速度计的标度系数不匹配,加速度计二阶误差系数,加速度计安装误差,因此载体的加速度,角速度,传递 到重力梯度仪的输出,使重力梯度仪的输出远大于检测质量激励的输出,该实验现象和重 力梯度仪的特性相一致。
Claims (5)
1.一种旋转加速度计重力梯度仪数字建模仿真方法,其特征在于,该方法包括以下步骤:
1)设置重力梯度仪的参数:
设置重力梯度仪加速度计安装圆盘的旋转角频率Ω;
设置四只加速度计的安装参数及输出模型参数:设置加速度计A1的径向安装距离R1,加速度计A1的初相位角β1z,加速度计A1的高度角β1x,加速度计A1的输入轴失准角θ1y、θ1z,加速度计A1的零偏k10,加速度计A1的线性标度系数k11,加速度计A1的二阶误差系数k12,k14,k15,k16,k17,k18,加速度计A1的电流转电压增益k1V/I;设置加速度计A2的径向安装距离R2,加速度计A2的初相位角β2z,加速度计A2的高度角β2x,加速度计A2的输入轴失准角θ2y、θ2z,加速度计A2的零偏k20,加速度计A2的线性标度系数k21,加速度计A2的二阶误差系数k22,k24,k25,k26,k27,k28,加速度计A2的电流转电压增益k2V/I;设置加速度计A3的径向安装距离R3,加速度计A3的初相位角β3z,加速度计A3的高度角β3x,加速度计A3的输入轴失准角θ3y、θ3z,加速度计A3的零偏k30,加速度计A3的线性标度系数k31,加速度计A3的二阶误差系数k32,k34,k35,k36,k37,k38,加速度计A3的电流转电压增益k3V/I;设置加速度计A4的径向安装距离R4,加速度计A4的初相位角β4z,加速度计A4的高度角β4x,加速度计A4的输入轴失准角θ4y、θ4z,加速度计A4的零偏k40,加速度计A4的线性标度系数k41,加速度计A4的二阶误差系数k42,k44,k45,k46,k47,k48,加速度计A4的电流转电压增益k4V/I;
设置加速度计的噪声系数α,b,ωT;
设置载体加速度向量fcmm,角速度向量ω,角加速度向量ωa;
设置采样频率fs,仿真总时间Ttotal;
2)根据步骤1)中加速度计的安装参数,计算t时刻加速度计A1在重力梯度仪测量坐标系中的位置向量计算t时刻加速度计A2在重力梯度仪测量坐标系中的位置向量计算t时刻加速度计A3在重力梯度仪测量坐标系中的位置向量计算t时刻加速度计A4在重力梯度仪测量坐标系中的位置向量
根据步骤1)中加速度计的安装参数,计算t时刻加速度计A1的标称安装坐标系的x轴单位方向向量τ1x,y轴单位方向向量τ1y,z轴单位方向向量τ1z;计算t时刻加速度计A2的标称安装坐标系的x轴单位方向向量τ2x,y轴单位方向向量τ2y,z轴单位方向向量τ2z;计算t时刻加速度计A3的标称安装坐标系的x轴单位方向向量τ3x,y轴单位方向向量τ3y,z轴单位方向向量τ3z;计算t时刻加速度计A4的标称安装坐标系的x轴单位方向向量τ4x,y轴单位方向向量τ4y,z轴单位方向向量τ4z;
根据步骤1)中加速度计的安装参数,计算t时刻加速度计A1的标称坐标系到加速度计A1的测量坐标系的变换矩阵C1;计算t时刻加速度计A2的标称坐标系到加速度计A2的测量坐标系的变换矩阵C2;计算t时刻加速度计A3的标称坐标系到加速度计A3的测量坐标系的变换矩阵C3;计算t时刻加速度计A4的标称坐标系到加速度计A4的测量坐标系的变换矩阵C4;
根据以下功率谱密度模型计算t时刻加速度计A1的噪声f1noise、加速度计A2的噪声f2noise、加速度计A3的噪声f3noise、加速度计A4的噪声f4noise:
Φ(f)noise=αf-b+ωT
式中,Φ(f)noise表示加速度计的噪声功率谱密度函数,f表示频率;
4)根据下式分别计算t时刻加速度计A1的比力向量f1、t时刻加速度计A2的比力向量f2、t时刻加速度计A3的比力向量f3、t时刻加速度计A4的比力向量f4:
式中,×为向量叉乘;G为万有引力常数;
5)根据下式分别计算t时刻四只加速度计在标称安装坐标系的x方向的比力、y方向的比力、z方向的比力:
式中f1x,f1y,f1z分别是加速度计A1在其标称安装坐标系的x方向,y方向,z方向的比力;f2x,f2y,f2z分别是加速度计A2在其标称安装坐标系的x方向,y方向,z方向的比力;f3x,f3y,f3z分别是加速度计A3在其标称安装坐标系的x方向,y方向,z方向的比力;f4x,f4y,f4z分别是加速度计A4在其标称安装坐标系的x方向,y方向,z方向的比力;式中·为向量点积;
6)根据下式分别计算t时刻四只加速度计在加速度计测量坐标系的输入轴的比力,输出轴的比力,摆轴的比力:
式中f1i,f1o,f1p分别是加速度计A1在其测量坐标系的输入轴、输出轴、摆轴的比力;f2i,f2o,f2p分别是加速度计A2在其测量坐标系的输入轴、输出轴、摆轴的比力,f3i,f3o,f3p分别是加速度计A3在其测量坐标系的输入轴、输出轴、摆轴的比力,f4i,f4o,f4p分别是加速度计A4在其测量坐标系的输入轴、输出轴、摆轴的比力;
7)根据下式分别计算t时刻加速度计A1的输出电压V1、加速度计A2的输出电压V2、加速度计A3的输出电压V3、加速度计A4的输出电压V4:
8)根据下式计算t时刻重力梯度仪解调前的信号Gout:
Gout=V1+V2-V3-V4
9)令t=t+1/fs,然后判断更新后的时间t是否小于仿真总时间Ttotal,如是,则返回步骤2),否则进入步骤10);
10)对信号Gout抗混叠滤波、降采样、带通滤波,然后进行正交幅值解调,恢复重力梯度信号。
2.根据权利要求1所述的一种旋转加速度计重力梯度仪数字建模仿真方法,其特征在于,所述步骤1)到5)中向量均为3维向量,矩阵为3×3方阵。
3.根据权利要求1所述的一种旋转加速度计重力梯度仪数字建模仿真方法,其特征在于,所述步骤2)和步骤5)中的加速度计标称安装坐标系的原点在加速度计的质心,加速度计标称安装坐标系的x轴定义为该点切线的方向,y轴定义为旋转圆盘的中心点到加速度计安装点的径向方向,z轴平行旋转圆盘轴线,方向垂直向下。
4.根据权利要求1、2或3所述的一种旋转加速度计重力梯度仪数字建模仿真方法,其特征在于,所述步骤2)和步骤6)中的加速度计测量坐标系的原点在加速度计的质心,加速度计测量坐标系的x轴定义为加速度计输入轴方向,y轴定义为加速度计输出轴的方向,z轴定义为加速度计摆轴的方向。
5.根据权利要求1、2或3所述的一种旋转加速度计重力梯度仪数字建模仿真方法,其特征在于,所述步骤1)和步骤2)中的重力梯度仪测量坐标系的原点位于旋转圆盘中心,重力梯度仪测量坐标系的x轴定义为旋转圆盘中心点到加速度计A1的初始位置,y轴定义为旋转圆盘的中心点到加速度计A3的初始位置,z轴定义为旋转圆盘轴线方向。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810599648.3A CN109085654B (zh) | 2018-06-11 | 2018-06-11 | 一种旋转加速度计重力梯度仪数字建模仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810599648.3A CN109085654B (zh) | 2018-06-11 | 2018-06-11 | 一种旋转加速度计重力梯度仪数字建模仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109085654A CN109085654A (zh) | 2018-12-25 |
CN109085654B true CN109085654B (zh) | 2020-01-07 |
Family
ID=64839930
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810599648.3A Active CN109085654B (zh) | 2018-06-11 | 2018-06-11 | 一种旋转加速度计重力梯度仪数字建模仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109085654B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109709628B (zh) * | 2019-02-15 | 2020-08-14 | 东南大学 | 一种旋转加速度计重力梯度仪标定方法 |
CN114089441B (zh) * | 2021-11-26 | 2022-08-30 | 华中科技大学 | 一种重力梯度仪测量系统数值仿真方法 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5357802A (en) * | 1993-05-07 | 1994-10-25 | Textron, Incorporated | Rotating accelerometer gradiometer |
US5922951A (en) * | 1997-06-11 | 1999-07-13 | The Broken Hill Proprietary Company Ltd. | Gravity gradiometer |
AU2003217957B2 (en) * | 2002-03-06 | 2008-07-31 | Lockheed Martin Corporation | Gravity gradiometer and method for calculating a gravity tensor with increased accuracy |
AUPS114702A0 (en) * | 2002-03-18 | 2002-04-18 | Bhp Billiton Innovation Pty Ltd | Enhancement of sensors for airborne operation |
CN102608668B (zh) * | 2011-12-19 | 2014-07-09 | 华中科技大学 | 重力梯度测量系统及其测量方法 |
US9519076B2 (en) * | 2014-02-20 | 2016-12-13 | Lockheed Martin Corporation | De-centralized control architecture for improved sensitivity of accelerometer-based gravity gradiometers |
CN104374401A (zh) * | 2014-10-15 | 2015-02-25 | 哈尔滨工程大学 | 一种捷联惯导初始对准中重力扰动的补偿方法 |
CN104459826B (zh) * | 2014-11-03 | 2017-09-26 | 东南大学 | 旋转加速度计重力梯度仪重力梯度信号仿真方法 |
CN104570142B (zh) * | 2015-01-05 | 2017-02-22 | 中国船舶重工集团公司第七0七研究所 | 一种重力梯度仪重力梯度测量信号的解调方法 |
CN105044795A (zh) * | 2015-05-29 | 2015-11-11 | 中国船舶重工集团公司第七0七研究所 | 利用光纤陀螺减小旋转加速度计重力梯度仪测量误差的方法 |
CN105717553B (zh) * | 2016-01-29 | 2018-01-16 | 东南大学 | 一种旋转加速度计重力梯度仪标定方法 |
-
2018
- 2018-06-11 CN CN201810599648.3A patent/CN109085654B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109085654A (zh) | 2018-12-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109766812B (zh) | 一种旋转加速度计重力梯度仪运动误差事后补偿方法 | |
CN109709628B (zh) | 一种旋转加速度计重力梯度仪标定方法 | |
CN104459826B (zh) | 旋转加速度计重力梯度仪重力梯度信号仿真方法 | |
CN102620892A (zh) | 一种转动部件的动平衡测试方法 | |
CN103323625B (zh) | 一种mems-imu中加速度计动态环境下的误差标定补偿方法 | |
CN109085654B (zh) | 一种旋转加速度计重力梯度仪数字建模仿真方法 | |
CN104155054A (zh) | 一种基于气浮扭摆台的转动惯量的频域检测方法 | |
CN107870371A (zh) | 一种动基座重力梯度仪自梯度补偿方法 | |
CN106840100A (zh) | 一种数字式倾角传感器及测量方法 | |
CN104597520B (zh) | 一种重力梯度仪重力梯度敏感器输入激励产生方法 | |
CN111650664B (zh) | 一种航空重力梯度仪实时重力梯度解调方法及装置 | |
CN113885098A (zh) | 一种重力敏感器低频频率响应误差在线建模及补偿方法 | |
CN110068876A (zh) | 基于载体自振动航空重力梯度仪运动误差补偿方法 | |
CN102778335A (zh) | 一种基于等效初相矢的各向异性转子全息动平衡方法 | |
CN105478245B (zh) | 基于主轴振动检测的双自由度精密离心机副轴动不平衡量辨识方法 | |
WO2020042696A1 (zh) | 一种动基座旋转加速度计重力梯度仪误差补偿装置及方法 | |
CN111879335A (zh) | 一种基于离心机的多位置陀螺仪漂移系数的标定方法 | |
CN108120439B (zh) | 一种三分量感应线圈姿态测量方法及装置 | |
CN108931824A (zh) | 一种动基座旋转加速度计重力梯度仪误差增益系数标定方法 | |
CN107576992A (zh) | 一种重力梯度仪自标定方法及离心梯度补偿方法 | |
CN112345199A (zh) | 一种暂冲式高速风洞迎角传感器振动影响修正方法 | |
CN105716626B (zh) | 一种悬浮类陀螺仪的定子旋转调制误差补偿方法 | |
Yu et al. | Mathematic model and error analysis of moving-base rotating accelerometer gravity gradiometer | |
CN114324977B (zh) | 一种基于全误差分析的在离心机上对惯性导航的加速度计标定的方法 | |
CN105842481A (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 |