CN109018433A - 基于角位移传感器的推力器推力在轨高精度标定方法 - Google Patents
基于角位移传感器的推力器推力在轨高精度标定方法 Download PDFInfo
- 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
Links
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/10—Artificial satellites; Systems of such satellites; Interplanetary vehicles
-
- 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/66—Arrangements 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)。
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)
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)
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 | 上海卫星工程研究所 | 基于线振动测量的卫星发动机在轨推力实时标定方法 |
-
2018
- 2018-06-28 CN CN201810684259.0A patent/CN109018433A/zh active Pending
Patent Citations (3)
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)
Title |
---|
张红军,等: "基于MME/KF的电推进器推力在轨标定算法", 《航天器环境工程》 * |
霍红庆,等: "MHD角速度传感器在航天任务中的应用", 《真空与低温》 * |
Cited By (5)
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 |