发明名称 一种测量星上有效载荷三轴角位移高频微振动的方法
摘要 一种测量星上有效载荷三轴角位移高频微振动的方法,本发明通过测量线加速度,利用了高频信号经过两次积分后产生的放大效应,从而以较低代价实现高精度测量,分析表明,对于50Hz幅值为10-6度的角位移,在1m测量基线的情况下,采用10-4g精度的加速度计即可,大大降低了仪器了精度,降低了测量成本;本发明提供了整段时间、时间窗口内两个前提下振动量的测量方法。
申请公布号 CN102023051B 申请公布日期 2012.06.13
申请号 CN201010297966.8 申请日期 2010.09.29
申请人 北京控制工程研究所 发明人 黎康;刘一武;杨孟飞;牟小刚;张勇智;张芸香
分类号 G01H17/00(2006.01)I;G01C1/00(2006.01)I;G01P15/18(2006.01)I 主分类号 G01H17/00(2006.01)I
代理机构 中国航天科技专利中心 11009 代理人 杨虹
主权项 1.一种测量星上有效载荷三轴角位移高频微振动的方法,其特征在于通过以下步骤实现:第一步,安装三轴加速度计组件三轴加速度计组件安装在有效载荷的安装支架平台上,即在有效载荷的前后两个表面各安装一套三轴加速度计组件;每套三轴加速度计组件包括三个用来测量有效载荷三个轴振动量的加速度计,三个加速度计的敏感轴两两正交于一点,其中两个加速度计与有效载荷的中心轴垂直;测量有效载荷同一方向振动的两个加速度计敏感轴相互平行,即两套三轴加速度计组件中加速度计的敏感轴两两平行,加速度计敏感轴交点的连线与有效载荷的中心轴平行;第二步,采集第一步中加速度计的测量得到的有效载荷前后两个平面的线加速度信号,利用公式(1)得到有效载荷的角加速度a(t),<maths num="0001"><![CDATA[<math><mrow><mi>a</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>a</mi><mn>1</mn></msub><mo>-</mo><msub><mi>a</mi><mn>2</mn></msub></mrow><mi>L</mi></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,a<sub>1</sub>、a<sub>2</sub>为两个敏感轴相互平行的加速度计测量得到的线加速度信号,L为测量基线长度,即两个敏感轴相互平行的加速度计之间的距离,t为测量时间;第三步,将第二步得到有效载荷的角加速度a(t)转化为频域功率谱密度函数A(f),利用公式(2)得到有效载荷的角速度功率密度<img file="FSB00000635855300012.GIF" wi="123" he="58" />利用公式(3)得到有效载荷角位移功率谱密度θ(f),<maths num="0002"><![CDATA[<math><mrow><mover><mi>&theta;</mi><mo>&CenterDot;</mo></mover><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>A</mi><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow></mrow><mrow><mn>2</mn><mi>&pi;f</mi></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><mi>&theta;</mi><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>A</mi><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow></mrow><msup><mrow><mo>(</mo><mn>2</mn><mi>&pi;f</mi><mo>)</mo></mrow><mn>2</mn></msup></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths>f为线加速度信号的频率;第四步,利用公式(4)得到有效载荷一个轴方向上的在测量时间内的角度振动量Φ<sub>rms</sub>,Φ<sup>2</sup><sub>rms</sub>=∑θ<sup>2</sup>(f)·Δf    (4);第五步,计算某个时间窗口T内有效载荷一个轴方向上的角度振动量,A5.1、利用公式(5)得到低于临界频率f<sub>0</sub>范围的角度振动量Δθ<sub>0</sub>,<maths num="0004"><![CDATA[<math><mrow><msub><mi>&Delta;&theta;</mi><mn>0</mn></msub><mo>=</mo><munder><mi>&Sigma;</mi><mi>i</mi></munder><mfrac><mi>T</mi><mn>2</mn></mfrac><msub><mover><mi>&theta;</mi><mo>&CenterDot;</mo></mover><mrow><mn>0</mn><mo>,</mo><mi>i</mi></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中<img file="FSB00000635855300022.GIF" wi="54" he="65" />为利用公式(2)得到的低于临界频率f<sub>0</sub>的各个抖动频率的角速度幅值,i∈[1,2,…I],I是低于临界频率f<sub>0</sub>的抖动频率的角速度幅值的总数,临界频率f<sub>0</sub>=1/(πT);A5.2、利用公式(6)得到高于临界频率f0范围的角度振动量Δθ<sub>1</sub>,<maths num="0005"><![CDATA[<math><mrow><msub><mi>&Delta;&theta;</mi><mn>1</mn></msub><mo>=</mo><msqrt><msubsup><mo>&Integral;</mo><msub><mi>f</mi><mn>0</mn></msub><mrow><mo>+</mo><mo>&infin;</mo></mrow></msubsup><msup><msub><mi>&theta;</mi><mn>1</mn></msub><mn>2</mn></msup><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mi>df</mi></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中θ<sub>1</sub>为利用公式(3)得到的高于临界频率f<sub>0</sub>的各个抖动频率的角度幅值;A5.3、利用公式(7)得到在时间窗口T内的有效载荷的角度振动量Δθ,Δθ=Δθ<sub>0</sub>+Δθ<sub>1</sub>    (7);第六步,重复第二步~第五步,得到有效载荷其余两个轴的角度振动量。
地址 100080 北京市2729信箱