CN109018433A - 基于角位移传感器的推力器推力在轨高精度标定方法 - Google Patents

基于角位移传感器的推力器推力在轨高精度标定方法 Download PDF

Info

Publication number
CN109018433A
CN109018433A CN201810684259.0A CN201810684259A CN109018433A CN 109018433 A CN109018433 A CN 109018433A CN 201810684259 A CN201810684259 A CN 201810684259A CN 109018433 A CN109018433 A CN 109018433A
Authority
CN
China
Prior art keywords
satellite
thruster
angular displacement
thrust
orbit
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.)
Pending
Application number
CN201810684259.0A
Other languages
English (en)
Inventor
俞洁
吕建民
刘伟
铁琳
杨立峰
陆国平
梁奕瑾
黄杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Institute of Satellite Engineering
Original Assignee
Shanghai Institute of Satellite Engineering
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shanghai Institute of Satellite Engineering filed Critical Shanghai Institute of Satellite Engineering
Priority to CN201810684259.0A priority Critical patent/CN109018433A/zh
Publication of CN109018433A publication Critical patent/CN109018433A/zh
Pending legal-status Critical Current

Links

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00Cosmonautic vehicles
    • B64G1/10Artificial satellites; Systems of such satellites; Interplanetary vehicles
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00Cosmonautic vehicles
    • B64G1/22Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
    • B64G1/66Arrangements or adaptations of apparatus or instruments, not otherwise provided for

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Physics & Mathematics (AREA)
  • Astronomy & Astrophysics (AREA)
  • General Physics & Mathematics (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

本发明提供了一种基于角位移传感器的推力器推力在轨高精度标定方法,包括以下步骤:步骤一,使用高精度角位移传感器采集推力器在轨工作时的卫星三轴角位移数据,包含开启和关闭全过程,并下传至地面处理;步骤二,根据高精度角位移传感器下传卫星三轴角位移数据计算卫星三轴姿态角速度、角加速度;步骤三,通过建立卫星姿态动力学方程,将发射前测得的推力器在卫星上的安装位置、角度及卫星的转动惯量作为已知量,根据步骤二获得的卫星三轴角速度及角加速度可以计算得出推力器的在轨推力。本发明可以进行推力器在轨推力高精度标定,保证高精度遥感卫星在轨使用过程中卫星姿态控制精度和成像质量。

Description

基于角位移传感器的推力器推力在轨高精度标定方法
技术领域
本发明涉及一种标定方法,具体地,涉及一种基于角位移传感器的推力器推力在轨高精度标定方法。
背景技术
对于卫星等空间飞行器而言,推力器为主要的姿态控制部件,其在轨推力精度直接影响卫星姿态控制效果。在轨已多次出现因推力器推力值估计错误造成卫星姿态控制异常,甚至发生卫星在轨解体的事件;由于推力器开启与关闭依靠电磁阀的张合实现,在开启和关闭过程中推力将产生非线性波动,这种现象对高精度遥感卫星姿态控制及遥感成像造成不利影响,且该波动由于持续时间较短(小于10ms),常规姿态敏感器(测量带宽小于20Hz)很难进行测量,需使用高精度高采样频率的传感器进行精确标定(角位移传感器测量频率范围2Hz~400Hz,角位移传感器测量精度优于0.006″);电磁阀的张合实现推力器开启与关闭,其张合过程中的动态数据是寿命预估的重要参数,只有使用高精度高采用频率传感器才能对其进行长期监视,提高寿命预估精度。因此,采用一种基于角位移传感器的推力器推力在轨高精度标定方法,显得至关重要。
发明内容
针对现有技术中的缺陷,本发明的目的是提供一种基于角位移传感器的推力器推力在轨高精度标定方法,其实现推力器在轨推力高精度标定,保证高精度遥感卫星在轨使用过程中卫星姿态控制精度和成像质量;通过对推力器开启及关闭过程中的推力变化分析,可以进行推力器工作性能退化评估,预测推力器工作寿命。
根据本发明的一个方面,提供一种基于角位移传感器的推力器推力在轨高精度标定方法,其特征在于,包括以下步骤:
步骤一,使用高精度角位移传感器采集推力器在轨工作时的卫星三轴角位移数据,包含开启和关闭全过程,并下传至地面处理;
步骤二,根据高精度角位移传感器下传卫星三轴角位移数据计算卫星三轴姿态角速度、角加速度;
步骤三,通过建立卫星姿态动力学方程,将发射前测得的推力器在卫星上的安装位置、角度及卫星的转动惯量作为已知量,根据步骤二获得的卫星三轴角速度及角加速度可以计算得出推力器的在轨推力。
优选地,步骤三包括以下步骤:
步骤三十一:建立卫星姿态动力学方程,如下式:
其中,J为卫星转动惯量,通过地面测试获得,如下式:
为卫星角加速度矢量,通过角位移传感器在轨测量获得,如下式:
为推力器距离卫星质心的位移,通过地面测试获得,在卫星三轴的分量如下式:
为推力器的推力矢量,其在卫星三轴的推力分离与安装角度有关,其三轴分量如下:
F为推力器推力值,所求目标值;α、β为推力器在卫星上的安装角度;
步骤三十二:根据角位移传感器在轨测量数据,进行推力器推力值解算;得到以下公式:
从而可以获得推力器含开启和关闭全过程推力值F(i)。
与现有技术相比,本发明具有如下的有益效果:本发明实现推力器在轨推力高精度标定,保证高精度遥感卫星在轨使用过程中卫星姿态控制精度和成像质量;通过对推力器开启及关闭过程中的推力变化分析,可以进行推力器工作性能退化评估,预测推力器工作寿命。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:
图1为本发明基于角位移传感器的推力器推力在轨高精度标定方法的流程图。
图2为推力器推力值与三轴之间的关系示意图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进。这些都属于本发明的保护范围。
如图1所示,本发明基于角位移传感器的推力器推力在轨高精度标定方法包括以下步骤:
步骤一,使用高精度角位移传感器采集推力器在轨工作时的卫星三轴角位移数据包含开启和关闭全过程,并下传至地面处理;
步骤二,根据高精度角位移传感器下传卫星三轴角位移数据计算卫星三轴姿态角速度角加速度
使用差分法可以获得卫星三轴姿态角速度、角加速度,如下式(1)和(2):
其中,为角位移传感器下传的第i个采样点数据;Δt为两个采样点之间的间隔。
步骤三,通过建立卫星姿态动力学方程,将发射前测得的推力器在卫星上的安装位置、角度及卫星的转动惯量作为已知量,根据步骤二获得的卫星三轴角速度及角加速度可以计算得出推力器的在轨推力;步骤三包括以下步骤:
步骤三十一:建立卫星姿态动力学方程,如下式(3):
其中,J为卫星转动惯量,通过地面测试获得,如下式(4):
为卫星角加速度矢量,通过角位移传感器在轨测量获得,如下式(5):
为推力器距离卫星质心的位移,通过地面测试获得,在卫星三轴的分量如下式(6):
为推力器的推力矢量,其在卫星三轴的推力分离与安装角度有关,其三轴分量如下(7):
F为推力器推力值,所求目标值;α、β为推力器在卫星上的安装角度,与三轴间关系如图2所示。
步骤三十二:根据角位移传感器在轨测量数据,进行推力器推力值解算。
将式(2)、(4)~(7)代入式(3),可以得到以下公式(8):
从而可以获得推力器含开启和关闭全过程推力值F(i)。
本发明可以进行推力器在轨推力高精度标定,保证高精度遥感卫星在轨使用过程中卫星姿态控制精度和成像质量;通过对推力器开启及关闭过程中的推力变化分析,可以进行推力器工作性能退化评估,预测推力器工作寿命。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变形或修改,这并不影响本发明的实质内容。

Claims (2)

1.一种基于角位移传感器的推力器推力在轨高精度标定方法,其特征在于,包括以下步骤:
步骤一,使用高精度角位移传感器采集推力器在轨工作时的卫星三轴角位移数据,包含开启和关闭全过程,并下传至地面处理;
步骤二,根据高精度角位移传感器下传卫星三轴角位移数据计算卫星三轴姿态角速度、角加速度;
步骤三,通过建立卫星姿态动力学方程,将发射前测得的推力器在卫星上的安装位置、角度及卫星的转动惯量作为已知量,根据步骤二获得的卫星三轴角速度及角加速度可以计算得出推力器的在轨推力。
2.根据权利要求1所述的基于角位移传感器的推力器推力在轨高精度标定方法,其特征在于,步骤三包括以下步骤:
步骤三十一:建立卫星姿态动力学方程,如下式:
其中,J为卫星转动惯量,通过地面测试获得,如下式:
为卫星角加速度矢量,通过角位移传感器在轨测量获得,如下式:
为推力器距离卫星质心的位移,通过地面测试获得,在卫星三轴的分量如下式:
为推力器的推力矢量,其在卫星三轴的推力分离与安装角度有关,其三轴分量如下:
F为推力器推力值,所求目标值;α、β为推力器在卫星上的安装角度;
步骤三十二:根据角位移传感器在轨测量数据,进行推力器推力值解算;得到以下公式:
从而可以获得推力器含开启和关闭全过程推力值F(i)。
CN201810684259.0A 2018-06-28 2018-06-28 基于角位移传感器的推力器推力在轨高精度标定方法 Pending CN109018433A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810684259.0A CN109018433A (zh) 2018-06-28 2018-06-28 基于角位移传感器的推力器推力在轨高精度标定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810684259.0A CN109018433A (zh) 2018-06-28 2018-06-28 基于角位移传感器的推力器推力在轨高精度标定方法

Publications (1)

Publication Number Publication Date
CN109018433A true CN109018433A (zh) 2018-12-18

Family

ID=65522001

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810684259.0A Pending CN109018433A (zh) 2018-06-28 2018-06-28 基于角位移传感器的推力器推力在轨高精度标定方法

Country Status (1)

Country Link
CN (1) CN109018433A (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109668579A (zh) * 2019-01-23 2019-04-23 张过 基于角位移传感器的星载推扫光学载荷高频误差补偿方法
CN112208794A (zh) * 2020-10-22 2021-01-12 上海卫星工程研究所 深空探测器的质量在轨测量方法、系统及介质
CN112660419A (zh) * 2020-12-24 2021-04-16 中国人民解放军61540部队 一种重力卫星质心在轨标定方法及系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103940431A (zh) * 2014-04-11 2014-07-23 北京空间飞行器总体设计部 基于gnss精密定轨的圆轨道切向小推力在轨标定方法
US20160202131A1 (en) * 2013-08-26 2016-07-14 University Of Florida Research Foundation Inc. Method and Apparatus for Measuring Thrust
CN108061660A (zh) * 2017-10-23 2018-05-22 上海卫星工程研究所 基于线振动测量的卫星发动机在轨推力实时标定方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20160202131A1 (en) * 2013-08-26 2016-07-14 University Of Florida Research Foundation Inc. Method and Apparatus for Measuring Thrust
CN103940431A (zh) * 2014-04-11 2014-07-23 北京空间飞行器总体设计部 基于gnss精密定轨的圆轨道切向小推力在轨标定方法
CN108061660A (zh) * 2017-10-23 2018-05-22 上海卫星工程研究所 基于线振动测量的卫星发动机在轨推力实时标定方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张红军,等: "基于MME/KF的电推进器推力在轨标定算法", 《航天器环境工程》 *
霍红庆,等: "MHD角速度传感器在航天任务中的应用", 《真空与低温》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109668579A (zh) * 2019-01-23 2019-04-23 张过 基于角位移传感器的星载推扫光学载荷高频误差补偿方法
CN112208794A (zh) * 2020-10-22 2021-01-12 上海卫星工程研究所 深空探测器的质量在轨测量方法、系统及介质
CN112208794B (zh) * 2020-10-22 2022-07-01 上海卫星工程研究所 深空探测器的质量在轨测量方法、系统及介质
CN112660419A (zh) * 2020-12-24 2021-04-16 中国人民解放军61540部队 一种重力卫星质心在轨标定方法及系统
CN112660419B (zh) * 2020-12-24 2022-08-30 中国人民解放军61540部队 一种重力卫星质心在轨标定方法及系统

Similar Documents

Publication Publication Date Title
CN109018433A (zh) 基于角位移传感器的推力器推力在轨高精度标定方法
CN111680870B (zh) 目标运动轨迹质量综合评估方法
Weigel et al. Thermodynamic correction of particle concentrations measured by underwing probes on fast-flying aircraft
Rivest A decentred predictor for circular-circular regression
WO2021102669A1 (zh) 超低轨道卫星轨道自主维持方法
Garcia et al. Unscented Kalman filter for spacecraft attitude estimation using quaternions and euler angles
Bonin et al. Comparison and application of wind retrieval algorithms for small unmanned aerial systems
CN112857306B (zh) 一种视频卫星任意视向点的连续太阳高度角确定方法
Toth et al. Thermal recoil force, telemetry, and the Pioneer anomaly
CN104501833A (zh) 一种基准不确定情况下加速度计组合误差系数标定方法
Ely et al. Method of data reduction and uncertainty estimation for platinum‐group element data using inductively coupled plasma‐mass spectrometry
McFarquhar et al. Airborne instrumentation needs for climate and atmospheric research
CN110687024A (zh) 一种基于图像的爆破粉尘量的测量方法
Svanidze et al. Angular distribution of extensive air showers by TEL array under GELATICA experiment
US9803978B2 (en) Laser rod surface elevation table device and method
CN109443786A (zh) 姿轨控发动机燃气阀响应时间及羽流温度测量装置及方法
Cui et al. An improved regression analysing method for multi-position calibration test of gyrowheel system
King Supercooled Liquid Water Content Instrument Analysis and Winter 2014 Data With Comparisons to the NASA Icing Remote Sensing System and Pilot Reports
Greer et al. A Saltation Sensor for the Martian Aqueous Habitat Reconnaissance Suite (MAHRS)
Hegazy et al. INS Stochastic Noise Impact on Circular Error Probability of Ballistic Missiles
Zheng et al. Simulation of Two-Dimensional Temperature Field via Kriging Method Based on Limited Conditioning Points in an Arc Dam
CN105786770B (zh) 一种基于光谱特征的人造卫星指向分析方法
Gayet Calibration of Johnson–Williams and PMS ASSP probes in a wind tunnel
Patino et al. Estimation of the aerodynamical parameters of an experimental airship
Song et al. Polar observations and model predictions during May 4, 1998, magnetopause, magnetosheath, and bow shock crossings

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20181218