CN102927995A - 一种在五个陀螺仪配置下的一致性故障诊断方法 - Google Patents
一种在五个陀螺仪配置下的一致性故障诊断方法 Download PDFInfo
- Publication number
- CN102927995A CN102927995A CN2012104092075A CN201210409207A CN102927995A CN 102927995 A CN102927995 A CN 102927995A CN 2012104092075 A CN2012104092075 A CN 2012104092075A CN 201210409207 A CN201210409207 A CN 201210409207A CN 102927995 A CN102927995 A CN 102927995A
- Authority
- CN
- China
- Prior art keywords
- delta
- theta
- gyroscope
- formula
- centerdot
- 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.)
- Granted
Links
Images
Landscapes
- Gyroscopes (AREA)
Abstract
Description
技术领域
本发明属于惯性导航领域,涉及一种故障诊断方法。
背景技术
为了提高捷联惯组的可靠性,目前有效的办法是采用冗余技术,其中需要对冗余结构中的同一信息测量进行故障诊断。第二炮兵工程学院的单斌和缪栋2004年第三期的《导弹与航天运载技术》杂志上,发表文献《捷联惯性测量组合的冗余设计和优化》,该文实现5个单自由度陀螺配置下的冗余优化设计,该配置采用的是圆锥体方式,即各表均匀放置在一倒圆锥体的圆周壁上,其平面夹角α=72°,半锥顶角θ=54.75°。虽然可以容忍二度故障(即允许两个表损坏),但由于不存在着两两垂直的正交表,在无故障的情况下,由于斜置表要参与导航运算,必须先投影到正交轴上,这样也增加了算法的复杂度。同时受到斜置表安装工艺的限制,投影系数的精度会受影响,从而使得导航精度也受影响。
北京航天自动控制研究所的巩庆海、吕新广等人提出了《一种在加速度表八表配置下的故障诊断方法》(专利号:ZL200910243093.x),该专利采用的冗余结构为:两套惯组配置相同,分别各有四只加速度表,其中三只正交安装(x1、y1、z1表),一只斜置安装(s1表),斜置加速度表可以通过投影到正交轴的方式对正交加速度表输出的一致性进行检验。该方法虽然在无故障的情况下,无需投影轴加表参与导航计算,导航精度也不受影响。但是由于采用了两套惯组进行故障诊断,不但对惯组的整体安装精度要求提高,而且也增加成本。
发明内容
本发明所要解决的技术问题是提供一种既能够降低成本,又不损失导航精度和可靠性的一致性故障诊断方法,
本发明进一步所要解决的技术问题是提供一种既能够降低成本,又不损失导航精度和可靠性的一致性故障诊断方法,该故障诊断方法能够提供无故障的角速度增量,用于导航计算。
本发明包括如下技术方案:
一种在五个陀螺仪配置下的一致性故障诊断方法,所述五个陀螺仪包括三个正交安装的陀螺仪Wx、Wy、Wz,和两个斜置的陀螺仪Ws、Wt;五个陀螺仪安装在一个本体结构内;所述一致性故障诊断方法包括如下步骤:
(1)构建角速度一致性故障判别公式,
角速度一致性故障判别公式包括如下五个公式:
公式1:
公式2:
公式3:
公式4:
公式5:
Es_x、Es_y、Es_z为陀螺仪Ws在正交轴的投影系数;
Et_x、Et_y、Et_z为陀螺仪Wt在正交轴的投影系数;
εGxy、εGyz、εGxz、εGs、εGt为一致性故障门限值;
Δt3为滚动累加周期;
(3)将第(2)步的结果代入公式1~5,判断公式1~5是否成立;并通过如下方法获得故障诊断结果;
当5个公式都成立时,则诊断为无故障;当5个公式都不成立时,则诊断为2个陀螺仪及以上故障;当公式1成立,其它公式不成立时,则诊断为陀螺仪Wt故障;当公式2成立,其它公式不成立时,则诊断为陀螺仪Ws故障;当公式3成立,其它公式不成立时,则诊断为陀螺仪Wz故障;当公式4成立,其它公式不成立时,则诊断为陀螺仪Wy故障;当公式5成立,其它公式不成立时,则诊断为陀螺仪Wx故障。
上述故障门限值根据陀螺仪的测量误差σ来确定,取6σ~9σ。
Δt3=25·τ,τ为控制周期。
在第(3)步后还包括如下步骤:
根据故障诊断结果选择陀螺仪进行导航计算;
当无故障、或者陀螺仪Ws故障、或者陀螺仪Wt故障时,采用三个正交安装的陀螺仪Wx、Wy、Wz输出的角速度增量进行导航计算;
当三个正交安装的陀螺仪Wx、Wy、Wz中的任一个故障时,用剩余两个无故障的正交安装的陀螺仪输出其所在坐标轴的角速度增量,然后利用四个无故障的陀螺仪输出的角速度增量获得故障陀螺仪所在坐标轴的角速度增量。
当陀螺仪Wy故障时,根据陀螺仪Ws、Wt、Wx、Wz输出的角速度增量获得陀螺仪Wy所在坐标轴的角速度增量
本发明与现有技术相比具有如下优点:
(1)本发明采用五个陀螺仪进行冗余配置后安装在一个本体结构内,构成一套惯组;只需在一套惯组内实现一致性故障诊断;采用上述冗余配置方式,不但降低成本,而且降低惯组整体的安装精度的要求。
(2)本发明在五个陀螺仪配置下,可以运用测量信息和公式组进行一致性故障诊断、切除故障信息、使用正确的陀螺信息进行导航计算,算法简单实用,能够检测出所有一度故障,可靠性高。
(3)本发明优选正交轴的陀螺信息进行导航计算,这是因为正交轴陀螺信息直接敏感三轴角速率信息,无需投影计算,可以降低由于投影系数误差导致的导航误差。
附图说明
图1为本发明的陀螺五表冗余结构配置图;
图2为本发明的故障诊断方法流程图。
具体实施方式
下面就结合附图对本发明做进一步介绍。
一、陀螺五表冗余结构
本发明的冗余结构采用五个陀螺仪组成一套惯组的形式进行配置。如图1所示,所述五个陀螺仪包括三个正交安装的陀螺仪Wx、Wy、Wz,和两个斜置的陀螺仪Ws、Wt;五个陀螺仪安装在一个本体结构内,从而构成一套惯组。五个陀螺仪固定在同一个方位安装基准面上。陀螺仪Wx、Wy、Wz的测量轴正向为箭体坐标系O-X1Y1Z1的坐标轴方向,陀螺仪Ws、Wt测量轴具体指向可以根据情况具体确定。其中,O-X1Y1Z1为箭体坐标系,OX1为纵轴,OY1为法向轴,OZ1为横向轴。陀螺仪Wx、Wy、Wz主要用于导航计算,斜置的陀螺仪Ws、Wt主要用于故障诊断和冗余计算。
二、一致性故障诊断方法
利用上述五个陀螺仪形成的冗余信息进行一致性故障诊断。如图2所示,本发明的故障诊断方法具体包括如下步骤:
(1)构建角速度一致性故障判别公式
根据空间投影原理确定如下角速度一致性故障判别公式:
公式1:
公式2:
公式3:
公式4:
公式5:
Es_x、Es_y、Es_z为陀螺仪Ws在正交轴的投影系数,其根据陀螺仪Ws的安装角度确定,为已知量;
Et_x、Et_y、Et_z为陀螺仪Wt在正交轴的投影系数,其根据陀螺仪Wt的安装角度确定,为已知量;
εGxy、εGyz、εGxz、εGs、εGt为一致性故障门限值;上述故障门限值根据陀螺仪的测量误差σ来确定,一般取6σ~9σ。例如,按9σ取值,确定出一致性故障门限值:εGxT=εGyT=εGzT=εGs=εGt=0.0901°。采用上述一致性故障门限值,可以保证其陀螺信息不会误切除。
Δt3为滚动累加周期。为降低高频干扰信号对故障诊断的影响,对角速度一致性判别进行滚动累加,累加时间段从开20ms中断以后开始,为防止长时间滚动累加误差造成影响,在允许起飞时刻进行一次清零,之后重新开始滚动累加。Δt3根据数学仿真等情况综合确定,例如累加周期可以设置为Δt3=25·τ(τ为控制周期20ms)。
(3)将第(2)步的结果代入公式1~5进行计算,并将计算结果按照角速度一致性故障判别定位表进行判别定位,得到故障判别定位结果。
所述角速度一致性判别定位表如下:
表1角速度一致性故障判别定位表
表中“√”表示公式成立;
“×”表示公式不成立;
对应角速度一致性故障判别定位表,对故障进行判别和定位;具体方法如下:当5个公式都成立时,则诊断为无故障;当5个公式都不成立时,则诊断为2个陀螺仪及以上故障;当公式1成立,其它公式不成立时,则诊断为陀螺仪Wt故障;当公式2成立,其它公式不成立时,则诊断为陀螺仪Ws故障;当公式3成立,其它公式不成立时,则诊断为陀螺仪Wz故障;当公式4成立,其它公式不成立时,则诊断为陀螺仪Wy故障;当公式5成立,其它公式不成立时,则诊断为陀螺仪Wx故障。
可以根据上述一致性故障判别定位结果进行导航计算。由于斜置表存在较大的安装误差,在导航计算过程中,优选三个正交安装的陀螺仪Wx、Wy、Wz作为导航计算的陀螺,但任意一路陀螺出现故障后,进行故障切除,并使用如下的方法来选择陀螺仪进行导航计算:
表2导航计算表
表2中,Wx1、Wy1,Wz1的值是通过计算得到。
Wy1对应的角速度增量值可通过陀螺仪Ws、Wt、Wx、Wz输出的角速度增量计算得到,利用如下公式:
本发明未详细说明部分属本领域技术人员公知常识。
Claims (5)
1.一种在五个陀螺仪配置下的一致性故障诊断方法,所述五个陀螺仪包括三个正交安装的陀螺仪Wx、Wy、Wz,和两个斜置的陀螺仪Ws、Wt;五个陀螺仪安装在一个本体结构内;其特征在于,所述一致性故障诊断方法包括如下步骤:
(1)构建角速度一致性故障判别公式,
角速度一致性故障判别公式包括如下五个公式:
公式1:
公式2:
公式3:
公式4:
公式5:
Es_x、Es_y、Es_z为陀螺仪Ws在正交轴的投影系数;
Et_x、Et_y、Et_z为陀螺仪Wt在正交轴的投影系数;
εGxy、εGyz、εGxz、εGs、εGt为一致性故障门限值;
Δt3为滚动累加周期;
(3)将第(2)步的结果代入公式1~5,判断公式1~5是否成立;并通过如下方法获得故障诊断结果;
当5个公式都成立时,则诊断为无故障;当5个公式都不成立时,则诊断为2个陀螺仪及以上故障;当公式1成立,其它公式不成立时,则诊断为陀螺仪Wt故障;当公式2成立,其它公式不成立时,则诊断为陀螺仪Ws故障;当公式3成立,其它公式不成立时,则诊断为陀螺仪Wz故障;当公式4成立,其它公式不成立时,则诊断为陀螺仪Wy故障;当公式5成立,其它公式不成立时,则诊断为陀螺仪Wx故障。
2.如权利要求1所述的一种在五个陀螺仪配置下的一致性故障诊断方法,其特征在于:上述故障门限值根据陀螺仪的测量误差σ来确定,例如取6σ~9σ。
3.如权利要求1所述的一种在五个陀螺仪配置下的一致性故障诊断方法,其特征在于:Δt3=25·τ,τ为控制周期。
4.如权利要求1所述的一种在五个陀螺仪配置下的一致性故障诊断方法,其特征在于,在第(3)步后还包括如下步骤:
根据故障诊断结果选择陀螺仪进行导航计算;
当无故障、或者陀螺仪Ws故障、或者陀螺仪Wt故障时,采用三个正交安装的陀螺仪Wx、Wy、Wz输出的角速度增量进行导航计算;
当三个正交安装的陀螺仪Wx、Wy、Wz中的任一个故障时,用剩余两个无故障的正交安装的陀螺仪输出其所在坐标轴的角速度增量,然后利用四个无故障的陀螺仪输出的角速度增量获得故障陀螺仪所在坐标轴的角速度增量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210409207.5A CN102927995B (zh) | 2012-10-24 | 2012-10-24 | 一种在五个陀螺仪配置下的一致性故障诊断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210409207.5A CN102927995B (zh) | 2012-10-24 | 2012-10-24 | 一种在五个陀螺仪配置下的一致性故障诊断方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102927995A true CN102927995A (zh) | 2013-02-13 |
CN102927995B CN102927995B (zh) | 2015-04-22 |
Family
ID=47642853
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210409207.5A Active CN102927995B (zh) | 2012-10-24 | 2012-10-24 | 一种在五个陀螺仪配置下的一致性故障诊断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102927995B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103234553A (zh) * | 2013-03-29 | 2013-08-07 | 北京控制工程研究所 | 一种陀螺测量***的故障诊断方法 |
CN103759740A (zh) * | 2014-01-03 | 2014-04-30 | 北京控制工程研究所 | 一种巡视器静态时陀螺故障自主诊断方法 |
CN103940451A (zh) * | 2014-04-30 | 2014-07-23 | 北京控制工程研究所 | 基于零空间向量自主优化选择的冗余陀螺故障定位方法 |
CN104061933A (zh) * | 2014-07-18 | 2014-09-24 | 上海新跃仪表厂 | 一种用于运载火箭的捷联惯性导航*** |
CN104898678A (zh) * | 2015-03-30 | 2015-09-09 | 北京航天自动控制研究所 | 一种运载火箭飞行控制的冗余诊断方法 |
CN107544466A (zh) * | 2017-09-15 | 2018-01-05 | 北京控制工程研究所 | 一种单框架控制力矩陀螺低速框架故障诊断方法 |
CN109781101A (zh) * | 2019-03-15 | 2019-05-21 | 重庆零壹空间航天科技有限公司 | 运载火箭惯性测量***及其组件的冗余方法 |
CN112945277A (zh) * | 2021-02-09 | 2021-06-11 | 北京航天自动控制研究所 | 一种在陀螺双五表配置下的故障判别方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005189083A (ja) * | 2003-12-25 | 2005-07-14 | Sumitomo Precision Prod Co Ltd | 自動車用角速度または自動車用角度センサ、およびそれらの故障検出方法 |
CN101672657A (zh) * | 2009-09-30 | 2010-03-17 | 北京航天自动控制研究所 | 基于多余度惯性测量信息的故障判断方法 |
CN101738499A (zh) * | 2009-12-24 | 2010-06-16 | 北京航天自动控制研究所 | 一种在加速度表八表配置下的故障诊断方法 |
-
2012
- 2012-10-24 CN CN201210409207.5A patent/CN102927995B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005189083A (ja) * | 2003-12-25 | 2005-07-14 | Sumitomo Precision Prod Co Ltd | 自動車用角速度または自動車用角度センサ、およびそれらの故障検出方法 |
CN101672657A (zh) * | 2009-09-30 | 2010-03-17 | 北京航天自动控制研究所 | 基于多余度惯性测量信息的故障判断方法 |
CN101738499A (zh) * | 2009-12-24 | 2010-06-16 | 北京航天自动控制研究所 | 一种在加速度表八表配置下的故障诊断方法 |
Non-Patent Citations (3)
Title |
---|
徐景硕,程传金,邸亚洲: ""冗余配置陀螺仪故障诊断方法研究"", 《科技信息》 * |
潘鸿飞,任尚云,袁立群: ""捷联惯导***故障诊断与处理技术研究"", 《上海航天》 * |
贾鹏,张洪钺: ""基于冗余惯性组件故障诊断方法的比较研究"", 《***仿真学报》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103234553A (zh) * | 2013-03-29 | 2013-08-07 | 北京控制工程研究所 | 一种陀螺测量***的故障诊断方法 |
CN103234553B (zh) * | 2013-03-29 | 2015-08-19 | 北京控制工程研究所 | 一种陀螺测量***的故障诊断方法 |
CN103759740A (zh) * | 2014-01-03 | 2014-04-30 | 北京控制工程研究所 | 一种巡视器静态时陀螺故障自主诊断方法 |
CN103759740B (zh) * | 2014-01-03 | 2015-07-08 | 北京控制工程研究所 | 一种巡视器静态时陀螺故障自主诊断方法 |
CN103940451A (zh) * | 2014-04-30 | 2014-07-23 | 北京控制工程研究所 | 基于零空间向量自主优化选择的冗余陀螺故障定位方法 |
CN104061933A (zh) * | 2014-07-18 | 2014-09-24 | 上海新跃仪表厂 | 一种用于运载火箭的捷联惯性导航*** |
CN104898678A (zh) * | 2015-03-30 | 2015-09-09 | 北京航天自动控制研究所 | 一种运载火箭飞行控制的冗余诊断方法 |
CN104898678B (zh) * | 2015-03-30 | 2017-06-06 | 北京航天自动控制研究所 | 一种运载火箭飞行控制的冗余诊断方法 |
CN107544466A (zh) * | 2017-09-15 | 2018-01-05 | 北京控制工程研究所 | 一种单框架控制力矩陀螺低速框架故障诊断方法 |
CN107544466B (zh) * | 2017-09-15 | 2019-08-09 | 北京控制工程研究所 | 一种单框架控制力矩陀螺低速框架故障诊断方法 |
CN109781101A (zh) * | 2019-03-15 | 2019-05-21 | 重庆零壹空间航天科技有限公司 | 运载火箭惯性测量***及其组件的冗余方法 |
CN112945277A (zh) * | 2021-02-09 | 2021-06-11 | 北京航天自动控制研究所 | 一种在陀螺双五表配置下的故障判别方法 |
CN112945277B (zh) * | 2021-02-09 | 2024-06-11 | 北京航天自动控制研究所 | 一种在陀螺双五表配置下的故障判别方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102927995B (zh) | 2015-04-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102927995B (zh) | 一种在五个陀螺仪配置下的一致性故障诊断方法 | |
CN103471615B (zh) | 一种双冗余惯导***快速故障检测方法 | |
CN110411438B (zh) | 一种基于多星敏感器的自适应组合确定卫星姿态角的方法 | |
CN109955872B (zh) | 一种用于高速铁路列车的实时定位方法及装置 | |
CN103616030A (zh) | 基于捷联惯导解算和零速校正的自主导航***定位方法 | |
CN107421534A (zh) | 一种冗余式捷联惯导***多故障隔离方法 | |
CN101858748A (zh) | 高空长航无人机的多传感器容错自主导航方法 | |
CN104750999B (zh) | 基于轨道面的地基探测设备过境计算目标及时段筛选方法 | |
CN102176159A (zh) | 一种基于状态观测器和等价空间的卫星姿控***故障诊断装置及方法 | |
CN101738499B (zh) | 一种在加速度表八表配置下的故障诊断方法 | |
CN112650281B (zh) | 多传感器三余度***、控制方法、无人机、介质及终端 | |
CN103162689A (zh) | 辅助车载定位***及车辆的辅助定位方法 | |
CN110567457B (zh) | 一种基于冗余的惯导自检测*** | |
CN107525505A (zh) | 机车轮对空转及滑行检测方法及*** | |
CN1932444A (zh) | 适用于高速旋转体的姿态测量方法 | |
CN102608641A (zh) | 单轴陀螺仪和单轴加速度计的车载组合导航***及方法 | |
CN103453904B (zh) | 一种惯性测量单元的冗余配置结构 | |
CN109975851A (zh) | 一种列车线路故障点精确定位方法与*** | |
CN112945229B (zh) | 一种双十表惯组主从冗余方法 | |
CN117405109A (zh) | 基于四元数球面线性加权的三套惯导***姿态表决方法 | |
CN109677508A (zh) | 一种车辆运动数据获取方法、装置、设备和存储介质 | |
CN104764464A (zh) | 一种利用全量信息进行飞行器冗余诊断的方法 | |
CN104833529A (zh) | 基于惯性导航汽车运动性能测试***及其测试方法 | |
CN104122412A (zh) | 一种基于北斗二代速度信息的加速度计标定方法 | |
CN111060098A (zh) | 一种矿井综采工作面用惯性导航定位定向*** |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |