发明名称 基于带单轴反转台离心机的挠性陀螺比力敏感项非线性测试方法
摘要 本发明公开了一种基于带单轴反转台离心机的挠性陀螺比力敏感项非线性测试方法,属于挠性陀螺过载项测试技术领域。所述的测试方法包括挠性陀螺和光纤陀螺的安装、挠性陀螺和光纤陀螺加电预热、静态数据采集和过载条件下的数据采集,获得陀螺开机重复性的多组输出数据进行非线性计算处理,得到过载加速度下挠性陀螺一次比力敏感项矢量和的非线性度。本发明可以得到挠性陀螺静态漂移模型中包含的比力敏感误差项与环境过载加速度之间的非线性关系,在实际应用中对该非线性变化进行补偿,从而提高挠性陀螺的测量精度和挠性惯组的导航精度,具有非常重要的实用价值。
申请公布号 CN102636185A 申请公布日期 2012.08.15
申请号 CN201210093684.5 申请日期 2012.03.31
申请人 北京航空航天大学 发明人 芦佳振;张春熹;李保国;叶勉
分类号 G01C25/00(2006.01)I 主分类号 G01C25/00(2006.01)I
代理机构 北京永创新实专利事务所 11121 代理人 姜荣丽
主权项 1.基于带单轴反转台离心机的挠性陀螺比力敏感项非线性测试方法,其特征在于包括如下步骤:第一步:挠性陀螺和光纤陀螺的安装:测试试验前,采用水平校准仪器调整随动反转台的台面和离心机的台面与水平面平行,然后将挠性陀螺和单轴光纤陀螺通过转接工装安装在随动反转台上,使挠性陀螺的敏感轴X轴与敏感轴Y轴均与水平面平行,光纤陀螺敏感轴与挠性陀螺敏感轴X轴平行,挠性陀螺的自转轴与反转台的转动轴重合;挠性陀螺和光纤陀螺的供电及信号通过随动反转台和离心机的滑环连接到供电电源和数据采集计算机;第二步:挠性陀螺和光纤陀螺加电预热:通过离心机的控制界面将离心机和随动反转台分别控制在其零位位置并保持静止,然后,挠性陀螺与光纤陀螺加电,预热稳定10min;第三步:静态数据采集:打开数据采集计算机中的采集软件开始采集挠性陀螺和光纤陀螺的输出,采集时间5min,得到初始静态下挠性陀螺X轴和光纤陀螺的静态输出数据;第四步:使离心机以角速度ω<sub>t</sub>=ω<sub>0</sub>旋转,随动反转台以角速度ω<sub>s</sub>=ω<sub>t</sub>-40旋转,挠性陀螺受到的环境过载加速度峰值为:<maths num="0001"><![CDATA[<math><mrow><mi>a</mi><mo>=</mo><msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>&omega;</mi><mi>t</mi></msub><mo>&CenterDot;</mo><mi>&pi;</mi></mrow><mn>180</mn></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>&CenterDot;</mo><mi>R</mi><mo>/</mo><msub><mi>g</mi><mn>0</mn></msub></mrow></math>]]></maths>其中,R是挠性陀螺自转轴到离心机转动轴的距离,单位:m;g<sub>0</sub>是当地标准重力加速度,单位m/s<sup>2</sup>;角速度ω<sub>t</sub>、ω<sub>s</sub>的单位为°/s;第五步:过载条件下的数据采集:离心机和随动反转台转速稳定后,采集挠性陀螺和光纤陀螺的输出数据15min,得到过载条件中挠性陀螺X轴和光纤陀螺的过载输出数据;第六步:重复第四步和第五步,设定不同的离心机角速度ω<sub>t</sub>和随动反转台的角速度ω<sub>s</sub>,即ω<sub>t</sub>=ω<sub>i</sub>,ω<sub>s</sub>=ω<sub>i</sub>-40,i=1,2……n,获得不同环境过载加速度a<sub>i</sub>下的挠性陀螺输出数据;其中<maths num="0002"><![CDATA[<math><mrow><msub><mi>a</mi><mi>i</mi></msub><mo>=</mo><msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>&omega;</mi><mi>i</mi></msub><mo>&CenterDot;</mo><mi>&pi;</mi></mrow><mn>180</mn></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>&CenterDot;</mo><mi>R</mi><mo>/</mo><msub><mi>g</mi><mn>0</mn></msub><mo>;</mo></mrow></math>]]></maths>第七步:测试完成后,离心机与随动反转台关机,停止数据采集,系统断电;第八步:重复第二步至第七步,获得陀螺开机重复性的多组输出数据;第九步:非线性计算处理,具体步骤为:(a)对采集到的初始静态下挠性陀螺X轴和光纤陀螺的数据求取平均值;(b)对采集到的过载条件中挠性陀螺X轴和光纤陀螺的数据分别减去其初始静态数据的平均值;(c)计算挠性陀螺X轴输出残差:<img file="FDA0000149564780000021.GIF" wi="1639" he="113" />(d)对挠性陀螺X轴输出残差在不同过载加速度a<sub>i</sub>下的数据D<sub>ij</sub>(X),截取整周期数据,周期<img file="FDA0000149564780000022.GIF" wi="172" he="84" />截取的数据个数N<sub>i</sub>满足<img file="FDA0000149564780000023.GIF" wi="266" he="50" />其中t为采样周期,m为大于500的正整数,<img file="FDA0000149564780000024.GIF" wi="34" he="46" />为过载加速度a<sub>i</sub>下的数据采集时间,ω<sub>si</sub>为随动反转台的角速度值,i=1,2……n;(e)对截取的各过载加速度下的挠性陀螺X轴输出残差数据进行傅里叶级数运算,分别得到傅里叶级数一次余弦和一次正弦谐波项系数A<sub>1xi</sub>和B<sub>1xi</sub>;<maths num="0003"><![CDATA[<math><mrow><msub><mi>A</mi><mrow><mn>1</mn><mi>xi</mi></mrow></msub><mo>=</mo><mn>2</mn><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>i</mi></msub></munderover><msub><mi>D</mi><mi>ij</mi></msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mi>cos</mi><mrow><mo>(</mo><msub><mi>j&pi;&omega;</mi><mi>si</mi></msub><mi>t</mi><mo>/</mo><mn>180</mn><mo>)</mo></mrow><mo>)</mo></mrow><mo>/</mo><msub><mi>N</mi><mi>i</mi></msub><mo>,</mo><msub><mi>B</mi><mrow><mn>1</mn><mi>xi</mi></mrow></msub><mo>=</mo><mn>2</mn><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>i</mi></msub></munderover><msub><mi>D</mi><mi>ij</mi></msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mi>sin</mi><mrow><mo>(</mo><msub><mi>j&pi;&omega;</mi><mi>si</mi></msub><mi>t</mi><mo>/</mo><mn>180</mn><mo>)</mo></mrow><mo>)</mo></mrow><mo>/</mo><msub><mi>N</mi><mi>i</mi></msub></mrow></math>]]></maths>其中,D<sub>ij</sub>(X)为过载加速度a<sub>i</sub>时采集到的挠性陀螺X轴第j个输出数据,ω<sub>si</sub>为过载加速度a<sub>i</sub>时随动反转台的角速度值,j=1,2……N<sub>i</sub>,i=1,2……n,n是用户选取的过载加速度测试点的数量;t为采样周期,N<sub>i</sub>为该第i个过载加速度段的挠性陀螺X轴输出残差数据个数;(f)求取不同过载加速a<sub>i</sub>的下挠性陀螺比力敏感项的矢量和:<maths num="0004"><![CDATA[<math><mrow><msub><mi>K</mi><mi>i</mi></msub><mo>=</mo><mfrac><mn>1</mn><msub><mi>a</mi><mi>i</mi></msub></mfrac><msqrt><msubsup><mi>A</mi><mrow><mn>1</mn><mi>xi</mi></mrow><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>B</mi><mrow><mn>1</mn><mi>xi</mi></mrow><mn>2</mn></msubsup></msqrt><mo>,</mo></mrow></math>]]></maths>(g)由下式计算非线性:<maths num="0005"><![CDATA[<math><mrow><msub><mi>&zeta;</mi><mi>i</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>K</mi><mi>i</mi></msub><mo>-</mo><msub><mi>K</mi><mn>1</mn></msub></mrow><msub><mi>K</mi><mn>1</mn></msub></mfrac><mo>&times;</mo><mn>100</mn><mo>%</mo></mrow></math>]]></maths>其中,i=1,2……n,ζ<sub>i</sub>为过载加速度a<sub>i</sub>下挠性陀螺比力敏感项矢量和的非线性度,K<sub>i</sub>、K<sub>1</sub>分别为过载加速度a<sub>i</sub>和a<sub>1</sub>下的挠性陀螺比力敏感项矢量和。
地址 100191 北京市海淀区学院路37号
您可能感兴趣的专利