CN107131949B - 一种利用星载加速度计间接测量微振动物理量参数的方法 - Google Patents
一种利用星载加速度计间接测量微振动物理量参数的方法 Download PDFInfo
- Publication number
- CN107131949B CN107131949B CN201710245123.5A CN201710245123A CN107131949B CN 107131949 B CN107131949 B CN 107131949B CN 201710245123 A CN201710245123 A CN 201710245123A CN 107131949 B CN107131949 B CN 107131949B
- Authority
- CN
- China
- Prior art keywords
- acceleration
- displacement
- time domain
- angular
- measurand
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种利用星载加速度计间接测量微振动物理量参数的方法,用于对被测对象进行微振动物理量参数的测量;首先在被测对象上选定一个安装平面,确定坐标系;然后在安装平面上选择p和q两点设置加速度计,加速度计的敏感轴与p和q连线相垂直、并在同一平面内;2个加速度计敏感到的加速度以数字信号下传到地面,地面处理后得两组微振动线加速度时域数据,分别进行FFT,并计算出线速度和线位移的幅度谱,线速度和线位移的幅度谱分别进行IFFT,得到线速度和线位移的时域数据;求出被测对象绕p,q两点确定的转轴线旋转的相对角加速度时域数据,进行FFT,得并计算出角速度和角位移的幅度谱,分别进行IFFT得到角速度和角位移的时域数据。
Description
技术领域
本发明属于非电量计量测试技术领域,尤其涉及一种利用星载加速度计间接测量微振动物理量参数的方法。
背景技术
卫星空间资源有着严格的限制,以尽可能少的资源实现任务目标是工程设计的基本原则。采用高精度加速度计测量航天器微振动是通用技术手段,获取微振动速度和位移通常通过对加速度一次积分和二次积分获得,由此则带来积分常数误差,如图3、图4和图5所示,存在明显图像畸变;同理,对时域角加速度进行一次积分和二次积分也可以得到时域角速度和角位移,存在的问题也是同样的。
发明内容
有鉴于此,本发明提供了一种利用星载加速度计间接测量微振动物理量参数的方法,本发明仅适用两个加速度计就能进行微振动物理量的测量,且采用时频域转换的方式进行计算,大大减小了误差。
为了达到上述目的,本发明的技术方案为:一种利用星载加速度计间接测量微振动物理量参数的方法,用于对被测对象进行微振动物理量参数的测量,其中微振动物理量参数包括被测对象的线速度、线位移、角速度和角位移;
步骤如下:
步骤(1)在被测对象上选定一个安装平面,确定坐标系;
步骤(2)在所述安装平面上选择p和q两点分别设置加速度计,p和q两点相距L,加速度计的敏感轴与p和q两点间的连线相垂直、并在同一平面内;
步骤(3)2个加速度计敏感所在环境的加速度,输出两组加速度模拟信号,两组加速度模拟信号经处理后以数字信号的形式下传到地面;
步骤(4)对下传的两组加速度数字信号解析后进行预处理,预处理包括:零偏修正、物理量变换、添加时间序列,得两组微振动线加速度时域数据:ap(kΔt)和aq(kΔt);
其中k为采样序号,k=0,l,2,…,N-1;Δt为采样时间间隔;N为样本量;
步骤(5)分别对测量得到的两组微振动线加速度时域数据进行快速傅里叶变换FFT,得到加速度幅度谱,进而计算出线速度和线位移的幅度谱,分别对线速度和线位移的幅度谱进行快速傅里叶逆变换IFFT,得到线速度和线位移的时域数据;
步骤(6)求出时间T=kΔt时被测对象绕p,q两点确定的转轴线旋转的相对角加速度时域数据β(kΔt):
步骤(7)对步骤(6)计算得到的角加速度时域数据β(kΔt)进行FFT,得到角加速度幅度谱,进而计算出角速度和角位移的幅度谱,分别对角速度和角位移的幅度谱进行IFFT得到角速度和角位移的时域数据。
进一步地,步骤(1)中所选安装平面具备被测对象的角振动特征。
进一步地,两个加速度计的测量方向平行,不平行度不大于0.5°。
进一步地,步骤(5)中,针对线加速度时域数据a(kΔt)的处理具体分为如下步骤:
其中a(kΔt)可以为ap(kΔt)和aq(kΔt);
Step501、针对线加速度时域数据a(kΔt)进行FFT得到线加速度幅度谱:
式中:j为虚数单位;i为加速度幅度值序列的序号;Δf为频率间隔,单位为Hz;单位为rad;X(iΔf)为被测对象频率在频率f=iΔf处的加速度FFT函数;
Step502、依据所述线加速度幅度谱,计算出线速度幅度谱和线位移幅度谱:
A(iΔf)为被测对象频率在f=iΔf处的加速度幅度值,为X(iΔf)的共轭复数;Ph(iΔf)为X(iΔf)在频率f=iΔf处的相位;Re(X(iΔf)为X(iΔf)在频率f=iΔf处的实部;Im(X(iΔf)为X(iΔf)在频率f=iΔf处的虚部;
式中:V(iΔf)为被测对象在频率f=iΔf处的速度幅度值;S(iΔf)为被测对象在频率f=iΔf处的位移幅度值;Phv(iΔf)为V(iΔf)在频率f=iΔf处的相位;Phs(iΔf)为S(iΔf)在频率f=iΔf处的相位;
Step503、分别使用iΔf、V(iΔf)、Phv(iΔf)和iΔf、S(iΔf)、Phs(iΔf)进行IFFT,可得到速度和位移的时域数据v(kΔt)和s(kΔt),:
k=0,1,2,...,N-1
所述样本量N为2的整数次幂。
有益效果:
本发明通过设置加速度计的位置,并以两个加速度计检测得到的加速度数据进行反演计算,就能够得到全参数(速度、位移、角加速度以及角速度,角位移)的时频域的数据,并且在计算的过程中摒弃了传统的时域积分算法,而采用时频域转换的方式进行计算,大大减小了误差,且由于仅使用了两个加速度计,实现简单。
附图说明
图1为本发明流程图;
图2为采用加速度计测量刚体角加速度的原理示意图;
图3为微振动加速度测量时域结果;
图4为速度量时域结果;
图5为位移量时域结果;
图6为IFFT运算后得到的时域速度结果;
图7为IFFT运算后得到的时域位移结果。
具体实施方式
下面结合附图并举实施例,对本发明进行详细描述。
实施例1、一种利用星载加速度计间接测量微振动物理量参数的方法,用于对被测对象进行微振动物理量参数的测量,其中微振动物理量参数包括被测对象的线速度、线位移、角速度和角位移;
如图2所示,本发明方法流程如下:
步骤(1)在被测对象上选定一个安装平面,确定坐标系;被测对象可以为卫星或其有效载荷,可视为刚体;所选安装平面代表被测对象的角振动特征。
步骤(2)在安装平面上选择p和q两点,p和q两点可以随机选择分别设置加速度计,p和q两点相距L,加速度计的敏感轴与p和q两点间的连线相垂直、并在同一平面内。安装方式如图2所示。
两个加速度计的测量方向平行,不平行度不大于0.5°。
步骤(3)2个加速度计敏感所在环境的加速度,输出两组加速度模拟信号,两组加速度模拟信号经处理后以数字信号的形式下传到地面。
步骤(4)对下传的两组加速度数字信号解析后进行预处理,预处理包括:零偏修正、物理量变换、添加时间序列,得两组微振动线加速度时域数据:ap(kΔt)和aq(kΔt);如图3为微振动加速度测量时域结果。
其中k为采样序号,k=0,l,2,…,N-1;Δt为采样时间间隔,s;N为样本量;
步骤(5)分别对测量得到的两组加速度时域数据进行快速傅里叶变换FFT,得到加速度幅度谱,进而计算出线速度和线位移的幅度谱,分别对线速度和线位移的幅度谱进行快速傅里叶逆变换IFFT,得到线速度和线位移的时域数据;
针对线加速度时域数据a(kΔt)的处理具体分为如下步骤:
其中a(kΔt)可以为ap(kΔt)和aq(kΔt);
Step501、针对线加速度时域数据a(kΔt)进行FFT得到线加速度幅度谱:
式中:j为虚数单位;i为加速度幅度值序列的序号;Δf为频率间隔,单位为Hz;单位为rad;X(iΔf)为被测对象频率在频率f=iΔf处的加速度FFT函数;
Step502、依据线加速度幅度谱,计算出线速度幅度谱和线位移幅度谱:
A(iΔf)为被测对象频率在f=iΔf处的加速度幅度值,为X(iΔf)的共轭复数;Ph(iΔf)为X(iΔf)在频率f=iΔf处的相位;Re(X(iΔf)为X(iΔf)在频率f=iΔf处的实部;Im(X(iΔf)为X(iΔf)在频率f=iΔf处的虚部;如图4为速度量时域结果。
式中:V(iΔf)为被测对象在频率f=iΔf处的速度幅度值;S(iΔf)为被测对象在频率f=iΔf处的位移幅度值;Phv(iΔf)为V(iΔf)在频率f=iΔf处的相位;Phs(iΔf)为S(iΔf)在频率f=iΔf处的相位;图5为位移量时域结果。
Step503、分别使用iΔf、V(iΔf)、Phv(iΔf)和iΔf、S(iΔf)、Phs(iΔf)进行IFFT,可得到速度和位移的时域数据v(kΔt)和s(kΔt),:
k=0,1,2,...,N-1
样本量N为2的整数次幂。
步骤(6)求出时间T=kΔt时被测对象绕p,q两点确定的转轴线旋转的相对角加速度时域数据β(kΔt):
当p、q两点间的线加速度相反时,转轴线位于p、q两点之间;当p、q的线加速度相同时,如图2所示,转轴线位于p、q两点连线的延长线上。
步骤(7)对计算得到的角加速度时域数据采用步骤(5)的方法进行FFT,得到角加速度幅度谱,进而计算出角速度和角位移的幅度谱,分别对角速度和角位移的幅度谱进行IFFT得到角速度和角位移的时域数据。
其中图6为IFFT运算后得到的时域速度结果;图7为IFFT运算后得到的时域位移结果。
综上,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (3)
1.一种利用星载加速度计间接测量微振动物理量参数的方法,其特征在于,用于对被测对象进行微振动物理量参数的测量,其中微振动物理量参数包括被测对象的线速度、线位移、角速度和角位移;
步骤如下:
步骤(1)在被测对象上选定一个安装平面,确定坐标系;
步骤(2)在所述安装平面上选择p和q两点分别设置加速度计,p和q两点相距L,加速度计的敏感轴与p和q两点间的连线相垂直、并在同一平面内;
步骤(3)2个加速度计敏感所在环境的加速度,输出两组加速度模拟信号,两组加速度模拟信号经处理后以数字信号的形式下传到地面;
步骤(4)对下传的两组加速度数字信号解析后进行预处理,预处理包括:零偏修正、物理量变换、添加时间序列,得两组微振动线加速度时域数据:ap(kΔt)和aq(kΔt);
其中k为采样序号,k=0,l,2,…,N-1;Δt为采样时间间隔;N为样本量;
步骤(5)分别对测量得到的两组微振动线加速度时域数据进行快速傅里叶变换FFT,得到加速度幅度谱,进而计算出线速度和线位移的幅度谱,分别对线速度和线位移的幅度谱进行快速傅里叶逆变换IFFT,得到线速度和线位移的时域数据;
步骤(6)求出时间T=kΔt时被测对象绕p,q两点确定的转轴线旋转的相对角加速度时域数据β(kΔt):
步骤(7)对步骤(6)计算得到的相对角加速度时域数据β(kΔt)进行FFT,得到角加速度幅度谱,进而计算出角速度和角位移的幅度谱,分别对角速度和角位移的幅度谱进行IFFT得到角速度和角位移的时域数据;
所述步骤(5)中,针对线加速度时域数据a(kΔt)的处理具体分为如下步骤:
其中a(kΔt)可以为ap(kΔt)和aq(kΔt);
Step501、针对线加速度时域数据a(kΔt)进行FFT得到线加速度幅度谱:
式中:j为虚数单位;i为加速度幅度值序列的序号;Δf为频率间隔,单位为Hz;X(iΔf)为被测对象频率在频率f=iΔf处的加速度FFT函数;
Step502、依据所述线加速度幅度谱,计算出线速度幅度谱和线位移幅度谱:
A(iΔf)为被测对象频率在f=iΔf处的加速度幅度值,为X(iΔf)的共轭复数;Ph(iΔf)为X(iΔf)在频率f=iΔf处的相位;Re(X(iΔf)为X(iΔf)在频率f=iΔf处的实部;Im(X(iΔf)为X(iΔf)在频率f=iΔf处的虚部;
式中:V(iΔf)为被测对象在频率f=iΔf处的速度幅度值;S(iΔf)为被测对象在频率f=iΔf处的位移幅度值;Phv(iΔf)为V(iΔf)在频率f=iΔf处的相位;Phs(iΔf)为S(iΔf)在频率f=iΔf处的相位;
Step503、分别使用iΔf、V(iΔf)、Phv(iΔf)和iΔf、S(iΔf)、Phs(iΔf)进行IFFT,可得到速度和位移的时域数据v(kΔt)和s(kΔt),:
所述样本量N为2的整数次幂。
2.如权利要求1所述的一种利用星载加速度计间接测量微振动物理量参数的方法,其特征在于,所述步骤(1)中所选安装平面具备被测对象的角振动特征。
3.如权利要求1所述的一种利用星载加速度计间接测量微振动物理量参数的方法,其特征在于,所述两个加速度计的测量方向平行,不平行度不大于0.5°。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710245123.5A CN107131949B (zh) | 2017-04-14 | 2017-04-14 | 一种利用星载加速度计间接测量微振动物理量参数的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710245123.5A CN107131949B (zh) | 2017-04-14 | 2017-04-14 | 一种利用星载加速度计间接测量微振动物理量参数的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107131949A CN107131949A (zh) | 2017-09-05 |
CN107131949B true CN107131949B (zh) | 2019-08-23 |
Family
ID=59715578
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710245123.5A Active CN107131949B (zh) | 2017-04-14 | 2017-04-14 | 一种利用星载加速度计间接测量微振动物理量参数的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107131949B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114754959B (zh) * | 2022-04-27 | 2023-07-18 | 清华大学 | 一种激光通信链路微振动试验的角位移信号生成方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102023051A (zh) * | 2010-09-29 | 2011-04-20 | 北京控制工程研究所 | 一种测量星上有效载荷三轴角位移高频微振动的方法 |
CN102878989A (zh) * | 2012-09-21 | 2013-01-16 | 航天东方红卫星有限公司 | 一种利用星载线加速度计进行三轴角振动测量的方法 |
CN103134947A (zh) * | 2012-12-21 | 2013-06-05 | 中国飞行试验研究院 | 航空发动机整机振动总量速度有效值数值计算方法 |
CN104406680A (zh) * | 2014-11-29 | 2015-03-11 | 国网山西省电力公司晋城供电公司 | 一种电力变压器表面测点振动加速度信号特征的提取方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH1130566A (ja) * | 1997-07-10 | 1999-02-02 | Masaaki Okuma | 振動特性解析装置 |
JP3398582B2 (ja) * | 1997-10-08 | 2003-04-21 | 富士通株式会社 | ディスク装置の筐体の振動解析方法及び解析装置 |
-
2017
- 2017-04-14 CN CN201710245123.5A patent/CN107131949B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102023051A (zh) * | 2010-09-29 | 2011-04-20 | 北京控制工程研究所 | 一种测量星上有效载荷三轴角位移高频微振动的方法 |
CN102878989A (zh) * | 2012-09-21 | 2013-01-16 | 航天东方红卫星有限公司 | 一种利用星载线加速度计进行三轴角振动测量的方法 |
CN103134947A (zh) * | 2012-12-21 | 2013-06-05 | 中国飞行试验研究院 | 航空发动机整机振动总量速度有效值数值计算方法 |
CN104406680A (zh) * | 2014-11-29 | 2015-03-11 | 国网山西省电力公司晋城供电公司 | 一种电力变压器表面测点振动加速度信号特征的提取方法 |
Non-Patent Citations (2)
Title |
---|
"基于振动加速度测量的振动速度和位移信号识别方法探讨";顾名坤 等;《机械科学与技术》;20110415;第30卷(第4期);第522-526页 |
"石英挠性加速度计测量航天器微振动的方法";赵伟 等;《航天器环境工程》;20160415;第33卷(第2期);第199-205页 |
Also Published As
Publication number | Publication date |
---|---|
CN107131949A (zh) | 2017-09-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104897962B (zh) | 基于互素感知的单频信号短样本高精度测频方法及其装置 | |
CN104697623B (zh) | 变速激励下叶片异步振动参数识别方法 | |
CN109856689B (zh) | 一种超导航磁梯度张量数据抑噪处理方法和*** | |
CN101216512A (zh) | 一种非正弦周期信号实时高精度检测方法 | |
CN104155054B (zh) | 一种基于气浮扭摆台的转动惯量的频域检测方法 | |
CN102661782B (zh) | 一种超低频振动参数量值的快速测量方法 | |
CN108875706A (zh) | 基于滑动平均与能量归集的海洋结构时频分析方法 | |
CN111624671B (zh) | 旋转加速度计重力梯度仪重力梯度解调相位角确定方法及装置 | |
CN107064559A (zh) | 一种基于角摇摆运动的sins加速度计频率特性测试方法 | |
CN111650664B (zh) | 一种航空重力梯度仪实时重力梯度解调方法及装置 | |
CN102269803B (zh) | 基于时间延迟的离散频谱低频成分的校正方法 | |
CN106546400A (zh) | 一种多输入多输出非高斯随机振动试验***及试验方法 | |
CN102878989A (zh) | 一种利用星载线加速度计进行三轴角振动测量的方法 | |
CN113885098A (zh) | 一种重力敏感器低频频率响应误差在线建模及补偿方法 | |
CN101738203A (zh) | 挠性陀螺仪静态漂移零次和一次加速度相关项误差模型最优位置标定方法 | |
CN107741312B (zh) | 一种非均布式结构振动场工程测试方法 | |
CN107131949B (zh) | 一种利用星载加速度计间接测量微振动物理量参数的方法 | |
Bekchanov et al. | Mathematical model of vibrating air pump unit | |
CN102735937A (zh) | 信号相位差测量的方法 | |
CN105699686B (zh) | 基于esmd的气固两相流颗粒速度检测方法 | |
CN105043667B (zh) | 转子不平衡振动信号幅值、相位实时计算方法 | |
CN109446552A (zh) | 多轴相关随机激励下结构疲劳寿命时域计算方法 | |
CN117686232A (zh) | 一种燃气轮机振动基频实时提取方法、装置及存储介质 | |
CN103438899A (zh) | 补偿惯性测量***在运动中产生的误差的方法及*** | |
CN107218996B (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 |