CN108750145B - 一种磁力矩器极性在轨检测方法 - Google Patents
一种磁力矩器极性在轨检测方法 Download PDFInfo
- Publication number
- CN108750145B CN108750145B CN201810316841.1A CN201810316841A CN108750145B CN 108750145 B CN108750145 B CN 108750145B CN 201810316841 A CN201810316841 A CN 201810316841A CN 108750145 B CN108750145 B CN 108750145B
- Authority
- CN
- China
- Prior art keywords
- magnetic
- kinetic energy
- magnetic torquer
- axis
- polarity
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 20
- 230000006698 induction Effects 0.000 claims abstract description 9
- 238000009434 installation Methods 0.000 claims description 2
- 238000000034 method Methods 0.000 abstract description 15
- 238000005259 measurement Methods 0.000 description 5
- 230000006340 racemization Effects 0.000 description 3
- 230000010355 oscillation Effects 0.000 description 2
- CBENFWSGALASAD-UHFFFAOYSA-N Ozone Chemical compound [O-][O+]=O CBENFWSGALASAD-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/244—Spacecraft control systems
- B64G1/245—Attitude control algorithms for spacecraft attitude control
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/36—Guiding or controlling apparatus, e.g. for attitude control using sensors, e.g. sun-sensors, horizon sensors
- B64G1/366—Guiding or controlling apparatus, e.g. for attitude control using sensors, e.g. sun-sensors, horizon sensors using magnetometers
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Radar, Positioning & Navigation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Automation & Control Theory (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明提供了一种磁力矩器极性在轨检测方法,在检测过程中,每个控制周期,通过卫星本体坐标系中的星体角速度和磁感应强度计算磁力矩器待检测轴需要执行的单轴消旋指令;在检测过程中累加星体转动动能,在每一个判断周期结束时计算星体在该判断周期中的平均动能;设置判定磁力矩器极性正确的平均动能减小次数阈值和判定磁力矩器极性相反的平均动能增加次数阈值,对每个判断周期计算得到的平均动能减小、增大的情况进行累计,判断其次数是否达到相应的阈值,从而判断磁力矩器的极性。本发明实现了磁力矩器极性的在轨检测。
Description
技术领域
本发明涉及一种航天器故障诊断方法,用于对磁力矩器在轨进行极性检测。
背景技术
磁力矩器广泛应用于多种卫星的姿态控制中,并且具有较高的可靠性,但是磁力矩器可靠性高不代表不会出现故障。如全球臭氧探测器TOMS-EP在1996年7月2日由于磁力矩器极性接错导致空间任务中断。在已知的卫星故障情况分析中,对于磁力矩器的地面检测完成之后,后续的地面检测等操作也可能使得磁力矩器装反。目前对磁力矩器的极性检测通常在地面整星测试时进行,如通过指北针或者磁强计来测量磁力矩器的极性,即当磁力矩器工作时,其周围的磁场改变,磁强计的输出信号随之改变,通过磁强计的输出判断磁力矩器的输出极性([2]倪亦斌,陈德祥.某种新型磁强计在星载部件测试中的应用研究[J].中国空间科学技术,2000,20(3):66-70.)。
现有方法是在地面进行磁力矩器极性检测,地面检测时,卫星所处的磁场环境稳定,并且可以通过外部放置在磁力矩器附近的磁强计测量磁力矩器工作时的磁场变化。若将该方法应用到磁力矩器极性的在轨检测中,面临着两个问题,首先,虽然卫星姿态测量与控制系统中包含磁强计,但是为了避免磁力矩器干扰其测量卫星所处环境的磁场,不能将卫星姿态测量与控制系统中的磁强计安装在磁力矩器附近,故测量磁力矩器的输出需要额外安装磁强计;其次,卫星在轨运行时,因为轨道位置和姿态发生变化,卫星所处的磁场环境相应变化,使得磁强计的测量结果并不能真实反映磁力矩器工作引起的磁场变化。所以该地面磁力矩器极性检测的方法不适用于在轨的磁力矩器极性检测。
发明内容
为了克服现有技术的不足,本发明提供一种采用磁力矩器单轴消旋的方法,基于卫星星体转动动能的变化情况进行判断,在轨完成磁力矩器的极性检测。本发明使用卫星姿态测量与控制系统中已有的磁强计等器件,通过磁力矩器工作时卫星的动能变化情况,对磁力矩器极性进行有效判断,不需要直接测量磁力矩器工作时的磁场变化。
本发明解决其技术问题所采用的技术方案包括以下步骤:
(1)记卫星本体坐标系中星体角速度为ω,其X、Y、Z三轴的分量为ωx、ωy、ωz;卫星本体坐标系中磁感应强度为B,其X、Y、Z三轴的分量为Bx、By、Bz;正装的磁力矩器的X、Y、Z三轴与卫星本体坐标系X、Y、Z三轴对应平行,记磁力矩器X、Y、Z三轴能够提供的最大磁矩分别为mxmax、mymax、mzmax;
每个控制周期Tc,通过卫星星体角速度和磁感应强度计算磁力矩器待检测轴需要执行的单轴消旋指令;
对磁力矩器X轴进行检测时,磁力矩器Y轴、Z轴执行的磁矩指令为0,磁力矩器X轴执行的磁矩指令mx=-mxmaxsign(Byωz-Bzωy);
对磁力矩器Y轴进行检测时,磁力矩器X轴、Z轴执行的磁矩指令为0,磁力矩器Y轴执行的磁矩指令my=-mymaxsign(Bzωx-Bxωz);
对磁力矩器Z轴进行检测时,磁力矩器X轴、Y轴执行的磁矩指令为0,磁力矩器Z轴执行的磁矩指令mz=-mzmaxsign(Bxωy-Byωx);
联立上述方程,计算得到ωx、ωy、ωz,进而得到当前控制周期卫星本体坐标系中星体角速度ω;
(3)记平均动能减小次数计数值为Cdec,初值为0;设定判定磁力矩器极性正确的平均动能减小次数阈值Ndec;记平均动能增加次数计数值为Cinc,初值为0;设定判定磁力矩器极性相反的平均动能增加次数阈值Ninc;第一个判断周期结束,以首次计算得到的平均动能作为参考平均动能;
(4)重复步骤(1)和(2),从第二个判断周期开始,每个判断周期结束时进行一次磁力矩器极性判断,若当前判断周期的平均动能小于参考平均动能,使用当前判断周期的平均动能更新参考平均动能,Cinc清零,Cdec加1,当Cdec达到Ndec时,判定磁力矩器该轴极性正确;若当前判断周期的平均动能不小于参考平均动能,则Cdec清零,Cinc加1,当Cinc达到Ninc时,判定磁力矩器该轴极性相反。
本发明的有益效果是:通过采用磁力矩器单轴消旋的方法,基于星体平均转动动能的变化情况进行判断,实现了磁力矩器极性的在轨检测。磁力矩器极性的在轨检测为后续的在轨故障处理奠定基础,增强了航天器姿态控制系统的可靠性。
附图说明
图1为磁力矩器极性判断流程图。
图2为算例平均动能变化曲线图。
具体实施方式
下面结合附图和实施例对本发明进一步说明,本发明包括但不仅限于下述实施例。
本发明包括以下步骤:
(1)单轴消旋指令计算
考虑磁力矩器三轴与卫星本体坐标系三轴对应平行的情况,在检测过程中,每个控制周期,通过卫星本体坐标系中的星体角速度和磁感应强度计算磁力矩器待检测轴需要执行的单轴消旋指令。
(2)平均转动动能计算
考虑到星体转动惯量、陀螺测量的角速度存在误差等原因,每个控制周期计算得到的转动动能会呈震荡变化的趋势,故使用一段时间的转动动能平均值来代替瞬时转动动能进行判断,将该段时间称为判断周期。即在检测过程中累加星体转动动能,在每一个判断周期结束时计算星体在该判断周期中的平均动能。
(3)磁力矩器极性判断
设置判定磁力矩器极性正确的平均动能减小次数阈值和判定磁力矩器极性相反的平均动能增加次数阈值,对每个判断周期计算得到的平均动能减小、增大的情况进行累计,判断其次数是否达到相应的阈值,从而判断磁力矩器的极性。
本发明通过判断使用磁力矩器单轴消旋时,卫星星体转动动能变化的情况,进行磁力矩器的极性检测。在工程实际中,无法获得准确的星体转动惯量,陀螺测量的角速度也存在噪声和常值偏差等误差,因此计算得到的转动动能会呈震荡上升或震荡下降,故本发明使用一段时间内的转动动能平均值来代替瞬时转动动能以解决此问题。
综上,本发明提出的在轨磁力矩器极性检测方法为:
a.单轴消旋指令计算
记控制周期为Tc,卫星本体坐标系中星体角速度为ω,其X、Y、Z三轴的分量为ωx,ωy,ωz;卫星本体坐标系中磁感应强度为B,其X、Y、Z三轴的分量为Bx,By,Bz。在本发明中认为磁力矩器正装,并且其X、Y、Z三轴与卫星本体坐标系X、Y、Z三轴对应平行,记磁力矩器X、Y、Z三轴能够提供的最大磁矩分别为mxmax、mymax、mzmax。
检测过程中,每个控制周期,通过卫星星体角速度和磁感应强度计算磁力矩器待检测轴需要执行的单轴消旋指令。
对磁力矩器X轴进行检测时,磁力矩器Y轴、Z轴执行的磁矩指令为0,磁力矩器X轴执行的磁矩指令mx为:
mx=-mxmaxsign(Byωz-Bzωy) (1)
对磁力矩器Y轴进行检测时,磁力矩器X轴、Z轴执行的磁矩指令为0,磁力矩器Y轴执行的磁矩指令my为:
my=-mymaxsign(Bzωx-Bxωz) (2)
对磁力矩器Z轴进行检测时,磁力矩器X轴、Y轴执行的磁矩指令为0,磁力矩器Z轴执行的磁矩指令mz为:
mz=-mzmaxsign(Bxωy-Byωx) (3)
b.平均转动动能计算
因为工程实际中,星体转动惯量、陀螺测量的角速度存在误差等原因,每个控制周期计算得到的转动动能会呈震荡变化趋势,故本发明使用一段时间的转动动能平均值来代替瞬时转动动能以解决此问题。将计算平均动能的该段时间称为判断周期,用包含的控制周期数Nj表示。在每一个判断周期结束时计算星体在该判断周期中的平均动能Ek_mean。每一个控制周期的瞬时转动动能Eki计算如下:
式中,
I——星体转动惯量。
每个判断周期Nj中的平均动能Ek_mean计算如下:
c.磁力矩器极性判断
在进行磁力矩器极性判断时,记平均动能减小次数计数值为Cdec,初值为0;记判定磁力矩器极性正确的平均动能减小次数阈值为Ndec;记平均动能增加次数计数值为Cinc,初值为0;记判定磁力矩器极性相反的平均动能增加次数阈值为Ninc。
磁力矩器极性判断流程如图1所示。即第一个判断周期结束,首次计算得到平均动能,记录该平均动能做为参考平均动能,不进行判断。从第二个判断周期开始,每个判断周期结束时判断一次,判断方法为:若当前判断周期的平均动能小于参考平均动能,使用当前判断周期的平均动能更新参考平均动能,Cinc清零,Cdec加1,当Cdec达到Ndec时,判定磁力矩器该轴极性正确;若当前判断周期的平均动能不小于参考平均动能,Cdec清零,Cinc加1,当Cinc达到Ninc时,判定磁力矩器该轴极性相反。
使用上述方法,依次对磁力矩器的三轴进行检测,完成磁力矩器的极性检测。
下面选择磁力矩器X轴装反为例进行说明。本实例中,磁力矩器X轴最大磁矩mxmax=0.46Am2,卫星的转动惯量矩阵I=diag(0.153,0.151,0.097)kg·m2,星体初始角速度为(10,10,10)°/s。
a.单轴消旋指令计算
取控制周期Tc=1s,检测过程中,每个控制周期通过卫星星体角速度和磁感应强度按照式(1)计算磁力矩器X轴需要执行的单轴消旋指令。
b.平均转动动能计算
取判断周期Nj=300,在每一个控制周期按照式(4)计算瞬时转动动能Eki,在每一个判断周期结束时按照式(5)计算星体在该判断周期中的平均动能Ek_mean。
c.磁力矩器极性判断
取判定磁力矩器极性正确的平均动能减小次数阈值Ndec=5;取判定磁力矩器极性相反的平均动能增加次数阈值Ninc=5。
该实例中,对磁力矩器X轴检测时,平均动能Ek_mean的变化曲线如图2所示,可以看到,平均动能增加次数计数值Cinc达到了5次,即达到了判定阈值,故判定磁力矩器X轴极性相反。
Claims (1)
1.一种磁力矩器极性在轨检测方法,其特征在于包括下述步骤:
(1)记卫星本体坐标系中星体角速度为ω,其X、Y、Z三轴的分量为ωx、ωy、ωz;卫星本体坐标系中磁感应强度为B,其X、Y、Z三轴的分量为Bx、By、Bz;正装的磁力矩器的X、Y、Z三轴与卫星本体坐标系X、Y、Z三轴对应平行,记磁力矩器X、Y、Z三轴能够提供的最大磁矩分别为mxmax、mymax、mzmax;
每个控制周期Tc,通过卫星星体角速度和磁感应强度计算磁力矩器待检测轴需要执行的单轴消旋指令;
对磁力矩器X轴进行检测时,磁力矩器Y轴、Z轴执行的磁矩指令为0,磁力矩器X轴执行的磁矩指令mx=-mxmaxsign(Byωz-Bzωy);
对磁力矩器Y轴进行检测时,磁力矩器X轴、Z轴执行的磁矩指令为0,磁力矩器Y轴执行的磁矩指令my=-mymaxsign(Bzωx-Bxωz);
对磁力矩器Z轴进行检测时,磁力矩器X轴、Y轴执行的磁矩指令为0,磁力矩器Z轴执行的磁矩指令mz=-mzmaxsign(Bxωy-Byωx);
(3)记平均动能减小次数计数值为Cdec,初值为0;设定判定磁力矩器极性正确的平均动能减小次数阈值Ndec;记平均动能增加次数计数值为Cinc,初值为0;设定判定磁力矩器极性相反的平均动能增加次数阈值Ninc;第一个判断周期结束,以首次计算得到的平均动能作为参考平均动能;
(4)重复步骤(1)和(2),从第二个判断周期开始,每个判断周期结束时进行一次磁力矩器极性判断,若当前判断周期的平均动能小于参考平均动能,使用当前判断周期的平均动能更新参考平均动能,Cinc清零,Cdec加1,当Cdec达到Ndec时,判定磁力矩器该轴极性正确;若当前判断周期的平均动能不小于参考平均动能,则Cdec清零,Cinc加1,当Cinc达到Ninc时,判定磁力矩器该轴极性相反。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810316841.1A CN108750145B (zh) | 2018-04-10 | 2018-04-10 | 一种磁力矩器极性在轨检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810316841.1A CN108750145B (zh) | 2018-04-10 | 2018-04-10 | 一种磁力矩器极性在轨检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108750145A CN108750145A (zh) | 2018-11-06 |
CN108750145B true CN108750145B (zh) | 2021-07-13 |
Family
ID=63981666
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810316841.1A Active CN108750145B (zh) | 2018-04-10 | 2018-04-10 | 一种磁力矩器极性在轨检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108750145B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109050982A (zh) * | 2018-08-30 | 2018-12-21 | 西北工业大学 | 磁力矩器极性在轨检测方法 |
CN112230174B (zh) * | 2020-09-14 | 2022-07-29 | 中国空间技术研究院 | 一种适应中高轨星载强磁场输出设备的在轨健康确定方法 |
CN112526416B (zh) * | 2020-11-13 | 2023-04-14 | 航天东方红卫星有限公司 | 一种磁极性测试方法及装置 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6356814B1 (en) * | 1999-02-03 | 2002-03-12 | Microcosm, Inc. | Spacecraft magnetic torquer feedback system |
CN100559126C (zh) * | 2008-01-31 | 2009-11-11 | 北京控制工程研究所 | 一种星敏感器极性检验方法 |
CN102114917B (zh) * | 2009-12-31 | 2012-11-21 | 北京控制工程研究所 | 一种提高磁力矩器控制精度的处理方法 |
CN203870231U (zh) * | 2014-05-06 | 2014-10-08 | 北京控制工程研究所 | 卫星用手持磁极性测试仪 |
CN106218923B (zh) * | 2016-07-27 | 2018-10-26 | 中国科学院长春光学精密机械与物理研究所 | 一种磁力矩器的控制方法 |
-
2018
- 2018-04-10 CN CN201810316841.1A patent/CN108750145B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN108750145A (zh) | 2018-11-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108750145B (zh) | 一种磁力矩器极性在轨检测方法 | |
MXPA00006237A (es) | Metodo y aparato para determinar la ubicacion de caracteristicas de una tuberia. | |
CN107421534B (zh) | 一种冗余式捷联惯导系统多故障隔离方法 | |
CN110702104B (zh) | 一种基于车辆零速检测的惯性导航误差修正方法 | |
CN100565115C (zh) | 多位置捷联寻北系统方位效应的标定方法 | |
RU2558724C2 (ru) | Устройство диагностического комплекса для определения положения трубопровода и способ определения относительного перемещения трубопровода по результатам двух и более инспекционных пропусков диагностического комплекса для определения положения трубопровода | |
CN108061560B (zh) | 一种动中通天线的修正惯导方法及其构成的混合导航方法 | |
CN106153069B (zh) | 自主导航系统中的姿态修正装置和方法 | |
CN103712622B (zh) | 基于惯性测量单元旋转的陀螺漂移估计补偿方法及装置 | |
CN104132662A (zh) | 基于零速修正的闭环卡尔曼滤波惯性定位方法 | |
JP2004212382A (ja) | ロボット掃除機の回転位置誤差補正方法 | |
CN102927995B (zh) | 一种在五个陀螺仪配置下的一致性故障诊断方法 | |
CN112414413B (zh) | 一种基于相对角动量的仅测角机动检测及跟踪方法 | |
CN112459721B (zh) | 旋转导向钻井工具测控系统的故障诊断方法、装置及应用 | |
CN101738203B (zh) | 挠性陀螺仪静态漂移零次和一次加速度相关项误差模型最优位置标定方法 | |
WO2015143546A1 (en) | Absolute vector gravimeter and methods of measuring an absolute gravity vector | |
CN115290080A (zh) | 基于无偏有限脉冲响应滤波器的管道机器人定位方法 | |
CN111829503A (zh) | 一种光纤陀螺阈值测试方法及装置 | |
CN107219335B (zh) | 基于复连续小波变换的管道连接器检测方法 | |
CN109677508A (zh) | 一种车辆运动数据获取方法、装置、设备和存储介质 | |
US8701297B2 (en) | Method of determining heading by turning an inertial device | |
JP4509006B2 (ja) | 人工衛星の姿勢決定系 | |
CN108089434A (zh) | 一种基于磁强计的皮纳卫星姿态捕获方法 | |
CN111947668A (zh) | 基于在线估计的木星探测器测角/测距组合导航方法 | |
CN114370885A (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 | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20221012 Address after: East 1st Road, Science and Technology Park, Jiangshan Town, Yinzhou District, Ningbo City, Zhejiang Province, 315153 Patentee after: Soaring Aerospace Technology (Ningbo) Co.,Ltd. Address before: 710072 No. 127 Youyi West Road, Shaanxi, Xi'an Patentee before: Northwestern Polytechnical University |