CN108061660A - 基于线振动测量的卫星发动机在轨推力实时标定方法 - Google Patents
基于线振动测量的卫星发动机在轨推力实时标定方法 Download PDFInfo
- Publication number
- CN108061660A CN108061660A CN201710993325.8A CN201710993325A CN108061660A CN 108061660 A CN108061660 A CN 108061660A CN 201710993325 A CN201710993325 A CN 201710993325A CN 108061660 A CN108061660 A CN 108061660A
- Authority
- CN
- China
- Prior art keywords
- orbit
- thrust
- linearly coupled
- satellite
- response
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M15/00—Testing of engines
- G01M15/14—Testing gas-turbine engines or jet-propulsion engines
Landscapes
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Combustion & Propulsion (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
- Testing Of Engines (AREA)
Abstract
本发明提供了一种基于线振动测量的卫星发动机在轨推力实时标定方法,所述基于线振动测量的卫星发动机在轨推力实时标定方法采用了地面线振动测试数据进行响应仿真有限元模型的修正,并标定推力与线振动间的传递函数;仿真结果数据结合在轨测试数据进行干扰剔除,获得发动机在轨推力带来的线振动响应;并利用推力与线振动间的传递函数进行发动机在轨推力反演,实时标定卫星发动机在轨推力。本发明可广泛应用于各类航天器的各类发动机/推力器在轨推力的实时标定,利于航天器在轨进行轨道控制和姿态控制的策略修正与优化。
Description
技术领域
本发明涉及一种推力实时标定方法,具体地,涉及一种基于线振动测量的卫星发动机在轨推力实时标定方法。
背景技术
卫星在轨状态下,因气液状态变化使得发动机推力发生变化,使得卫星在轨轨道保持及变轨操作时存在较大风险,特别是在推进剂消耗而导致卫星质量特性发生变化时,发动机推力变化使得卫星在轨操控存在极大的风险,甚至关乎到高精度卫星的成败。然而,现阶段卫星发动机推力在轨尚无有效的标定方法,使得卫星变轨和轨道保持时存在一定风险,特别是卫星惯量与推力比较小,且卫星需要具备较长时间轨道机动或较大幅度的变轨操控时,卫星发动机在轨推力实时标定极为重要,是必须要重视的关乎卫星成败的重要问题。
因此,如何准确和高效的进行卫星发动机在轨推力实时标定,无论对卫星还是其它航天器均是极为重要的研究课题。
发明内容
针对现有技术中的缺陷,本发明的目的是提供一种基于线振动测量的卫星发动机在轨推力实时标定方法,其可广泛应用于各类航天器的各类发动机/推力器在轨推力的实时标定,利于航天器在轨进行轨道控制和姿态控制的策略修正与优化。
根据本发明的一个方面,提供一种基于线振动测量的卫星发动机在轨推力实时标定方法,其特征在于,所述基于线振动测量的卫星发动机在轨推力实时标定方法采用了地面线振动测试数据进行响应仿真有限元模型的修正,并标定推力与线振动间的传递函数;仿真结果数据结合在轨测试数据进行干扰剔除,获得发动机在轨推力带来的线振动响应;并利用推力与线振动间的传递函数进行发动机在轨推力反演,实时标定卫星发动机在轨推力。
优选地,所述基于线振动测量的卫星发动机在轨推力实时标定方法包括以下步骤:
步骤一,根据卫星构型及推力器布局,结合地面试验数据和有限元理论建立卫星结构有限元仿真建模;
步骤二,依据发动机推力热标试验数据,以及发动机推力主传递路径上设置线振动传感器,进行发动机推力地面模拟响应测试;
步骤三,依据发动机推力地面模拟响应测试结果,优化地面试验状态仿真模型,获得与校准推力器推力与线振动响应之间的传递函数,并建立在轨状态仿真模型进行在轨不同工况下的推力器点火响应仿真;
步骤四,设置与地面试验相同空间位置的线振动传感器进行在轨试验,将在轨不同工况下推力器点火时的线振动响应和卫星姿控数据传输至地面,利用该数据对比推力器点火响应仿真结果数据,剔除卫星附件、结构本体和姿态晃动,得到推力器点火在主传递路径上的响应;
步骤五,根据推力器点火响应仿真结果及推力器推力与线振动响应之间的传递函数,利用在轨线振动响应反演出推力器推力及其变化曲线,实现在轨实时标定卫星发动机推力的方法。
优选地,所述地面线振动测试数据为发动机推力地面模拟响应测试结果,是以发动机推力热标试验数据为输入,以发动机推力在星体上的主传递路径上设置的线振动传感器为测量部件,进行振动和冲击试验获得线振动传感器的测试数据。
优选地,所述推力器推力与线振动响应间的传递函数为不同推力输入条件与线振动传感器不同数值之间的对应关系。
优选地,所述剔除卫星附件、结构本体和姿态晃动等干扰是线振动传感器在轨实测数据结合仿真模型仿真结果,剥离卫星附件、结构本体和姿态晃动的频率对应的响应。
与现有技术相比,本发明具有如下的有益效果:本发明可广泛应用于各类航天器的各类发动机/推力器在轨推力的实时标定,利于航天器在轨进行轨道控制和姿态控制的策略修正与优化。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:
图1为本发明基于线振动测量的卫星发动机在轨推力实时标定方法的流程图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进。这些都属于本发明的保护范围。
如图1所示,本发明基于线振动测量的卫星发动机在轨推力实时标定方法采用了地面线振动测试数据进行响应仿真有限元模型的修正,并标定推力与线振动间的传递函数;仿真结果数据结合在轨测试数据进行干扰剔除,获得发动机在轨推力带来的线振动响应;并利用推力与线振动间的传递函数进行发动机在轨推力反演,实时标定卫星发动机在轨推力。
本发明基于线振动测量的卫星发动机在轨推力实时标定方法包括以下步骤:
步骤一,根据卫星构型及推力器布局,结合地面试验数据和有限元理论建立卫星结构有限元仿真建模;
步骤二,依据发动机推力热标试验数据,以及发动机推力主传递路径上设置线振动传感器,进行发动机推力地面模拟响应测试;
步骤三,依据发动机推力地面模拟响应测试结果,优化地面试验状态仿真模型,获得与校准推力器推力与线振动响应之间的传递函数,并建立在轨状态仿真模型进行在轨不同工况下的推力器点火响应仿真;
步骤四,设置与地面试验相同空间位置的线振动传感器进行在轨试验,将在轨不同工况下推力器点火时的线振动响应和卫星姿控数据传输至地面,利用该数据对比推力器点火响应仿真结果数据,剔除卫星附件、结构本体和姿态晃动等干扰,得到推力器点火在主传递路径上的响应;
步骤五,根据推力器点火响应仿真结果及推力器推力与线振动响应之间的传递函数,利用在轨线振动响应反演出推力器推力及其变化曲线,实现在轨实时标定卫星发动机推力的方法。
所述地面线振动测试数据为发动机推力地面模拟响应测试结果,是以发动机推力热标试验数据为输入,以发动机推力在星体上的主传递路径上设置的线振动传感器为测量部件,进行振动和冲击试验获得线振动传感器的测试数据。
所述推力器推力与线振动响应间的传递函数为不同推力输入条件与线振动传感器不同数值之间的对应关系。
所述剔除卫星附件、结构本体和姿态晃动等干扰是线振动传感器在轨实测数据结合仿真模型仿真结果,剥离卫星附件、结构本体和姿态晃动的频率对应的响应,即为发动机在轨推力带来的响应。
本发明从工程实际角度出发,采用了地面线振动测试数据进行响应仿真模型修正,仿真结果数据结合在轨测试数据进行干扰剔除,提出了一种利用线振动测量数据进行实时标定卫星发动机在轨推力的方法。所述方法包括:根据卫星构型及推力器布局,结合地面试验数据和有限元理论建立卫星结构有限元仿真建模;依据发动机推力热标试验数据,以及发动机推力主传递路径上设置线振动传感器,进行发动机推力地面模拟响应测试;依据发动机推力地面模拟响应测试结果,优化地面试验状态仿真模型,获得与校准推力器推力与线振动响应之间的传递函数,并建立在轨状态仿真模型进行在轨不同工况下的推力器点火响应仿真;将在轨不同工况下推力器点火响应和卫星姿控数据传输至地面,利用该数据对比推力器点火响应仿真结果数据,剔除卫星附件、结构本体和姿态晃动等干扰,得到推力器点火在主传递路径上的响应;根据推力器点火响应仿真及推力器推力与线振动响应之间的传递函数,利用在轨线振动响应反演出推力器推力及其变化曲线,实现在轨实时标定卫星发动机推力的方法。本发明可广泛应用于各类航天器的各类发动机/推力器在轨推力的实时标定,利于航天器在轨进行轨道控制和姿态控制的策略修正与优化。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变形或修改,这并不影响本发明的实质内容。
Claims (5)
1.一种基于线振动测量的卫星发动机在轨推力实时标定方法,其特征在于,所述基于线振动测量的卫星发动机在轨推力实时标定方法采用了地面线振动测试数据进行响应仿真有限元模型的修正,并标定推力与线振动间的传递函数;仿真结果数据结合在轨测试数据进行干扰剔除,获得发动机在轨推力带来的线振动响应;并利用推力与线振动间的传递函数进行发动机在轨推力反演,实时标定卫星发动机在轨推力。
2.根据权利要求1所述的基于线振动测量的卫星发动机在轨推力实时标定方法,其特征在于,所述基于线振动测量的卫星发动机在轨推力实时标定方法包括以下步骤:
步骤一,根据卫星构型及推力器布局,结合地面试验数据和有限元理论建立卫星结构有限元仿真建模;
步骤二,依据发动机推力热标试验数据,以及发动机推力主传递路径上设置线振动传感器,进行发动机推力地面模拟响应测试;
步骤三,依据发动机推力地面模拟响应测试结果,优化地面试验状态仿真模型,获得与校准推力器推力与线振动响应之间的传递函数,并建立在轨状态仿真模型进行在轨不同工况下的推力器点火响应仿真;
步骤四,设置与地面试验相同空间位置的线振动传感器进行在轨试验,将在轨不同工况下推力器点火时的线振动响应和卫星姿控数据传输至地面,利用该数据对比推力器点火响应仿真结果数据,剔除卫星附件、结构本体和姿态晃动,得到推力器点火在主传递路径上的响应;
步骤五,根据推力器点火响应仿真结果及推力器推力与线振动响应之间的传递函数,利用在轨线振动响应反演出推力器推力及其变化曲线,实现在轨实时标定卫星发动机推力的方法。
3.根据权利要求2所述的基于线振动测量的卫星发动机在轨推力实时标定方法,其特征在于,所述地面线振动测试数据为发动机推力地面模拟响应测试结果,是以发动机推力热标试验数据为输入,以发动机推力在星体上的主传递路径上设置的线振动传感器为测量部件,进行振动和冲击试验获得线振动传感器的测试数据。
4.根据权利要求2所述的基于线振动测量的卫星发动机在轨推力实时标定方法,其特征在于,所述推力器推力与线振动响应间的传递函数为不同推力输入条件与线振动传感器不同数值之间的对应关系。
5.根据权利要求2所述的基于线振动测量的卫星发动机在轨推力实时标定方法,其特征在于,所述剔除卫星附件、结构本体和姿态晃动等干扰是线振动传感器在轨实测数据结合仿真模型仿真结果,剥离卫星附件、结构本体和姿态晃动的频率对应的响应。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710993325.8A CN108061660B (zh) | 2017-10-23 | 2017-10-23 | 基于线振动测量的卫星发动机在轨推力实时标定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710993325.8A CN108061660B (zh) | 2017-10-23 | 2017-10-23 | 基于线振动测量的卫星发动机在轨推力实时标定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108061660A true CN108061660A (zh) | 2018-05-22 |
CN108061660B CN108061660B (zh) | 2019-09-17 |
Family
ID=62136968
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710993325.8A Active CN108061660B (zh) | 2017-10-23 | 2017-10-23 | 基于线振动测量的卫星发动机在轨推力实时标定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108061660B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109018433A (zh) * | 2018-06-28 | 2018-12-18 | 上海卫星工程研究所 | 基于角位移传感器的推力器推力在轨高精度标定方法 |
CN112208799A (zh) * | 2020-10-22 | 2021-01-12 | 上海卫星工程研究所 | 深空探测器轨控发动机推力方向偏差在轨标定方法和系统 |
CN112393835A (zh) * | 2020-11-03 | 2021-02-23 | 西北工业大学深圳研究院 | 一种基于扩展卡尔曼滤波的小卫星在轨推力标定方法 |
CN112520070A (zh) * | 2020-12-07 | 2021-03-19 | 上海卫星工程研究所 | 深空探测器推力矢量实时修正方法和系统 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN201277910Y (zh) * | 2008-10-22 | 2009-07-22 | 沈阳黎明航空发动机(集团)有限责任公司 | 一种航空发动机试车参数测量系统 |
CN103488812A (zh) * | 2013-08-12 | 2014-01-01 | 上海卫星工程研究所 | 卫星在轨微振动仿真建模修正方法 |
CN103712759A (zh) * | 2014-01-06 | 2014-04-09 | 北京卫星环境工程研究所 | 航天器整星微振动力学环境地面试验方法 |
CN104833466A (zh) * | 2015-04-30 | 2015-08-12 | 北京航空航天大学 | 一种航天器地面测试与在轨微振动力学环境映射方法 |
CN104932509A (zh) * | 2015-05-15 | 2015-09-23 | 上海新跃仪表厂 | 挠性航天器主动振动抑制地面测试系统 |
CN105446347A (zh) * | 2015-11-30 | 2016-03-30 | 上海卫星工程研究所 | 针对卫星太阳电池阵的在轨模态辨识系统及方法 |
CN107202660A (zh) * | 2017-06-06 | 2017-09-26 | 西安航天动力试验技术研究所 | 4‑25n姿控发动机真空热环境稳态推力校准测量装置 |
-
2017
- 2017-10-23 CN CN201710993325.8A patent/CN108061660B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN201277910Y (zh) * | 2008-10-22 | 2009-07-22 | 沈阳黎明航空发动机(集团)有限责任公司 | 一种航空发动机试车参数测量系统 |
CN103488812A (zh) * | 2013-08-12 | 2014-01-01 | 上海卫星工程研究所 | 卫星在轨微振动仿真建模修正方法 |
CN103712759A (zh) * | 2014-01-06 | 2014-04-09 | 北京卫星环境工程研究所 | 航天器整星微振动力学环境地面试验方法 |
CN104833466A (zh) * | 2015-04-30 | 2015-08-12 | 北京航空航天大学 | 一种航天器地面测试与在轨微振动力学环境映射方法 |
CN104932509A (zh) * | 2015-05-15 | 2015-09-23 | 上海新跃仪表厂 | 挠性航天器主动振动抑制地面测试系统 |
CN105446347A (zh) * | 2015-11-30 | 2016-03-30 | 上海卫星工程研究所 | 针对卫星太阳电池阵的在轨模态辨识系统及方法 |
CN107202660A (zh) * | 2017-06-06 | 2017-09-26 | 西安航天动力试验技术研究所 | 4‑25n姿控发动机真空热环境稳态推力校准测量装置 |
Non-Patent Citations (2)
Title |
---|
PETER J.WIKTOR: "On-orbit thruster calibration", 《JOURNAL OF GUIDANCE,CONTROL,AND DYNAMICS》 * |
张红军 等: "基于MME/KF的电推进器推力在轨标定算法", 《航天器环境工程》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109018433A (zh) * | 2018-06-28 | 2018-12-18 | 上海卫星工程研究所 | 基于角位移传感器的推力器推力在轨高精度标定方法 |
CN112208799A (zh) * | 2020-10-22 | 2021-01-12 | 上海卫星工程研究所 | 深空探测器轨控发动机推力方向偏差在轨标定方法和系统 |
CN112393835A (zh) * | 2020-11-03 | 2021-02-23 | 西北工业大学深圳研究院 | 一种基于扩展卡尔曼滤波的小卫星在轨推力标定方法 |
CN112520070A (zh) * | 2020-12-07 | 2021-03-19 | 上海卫星工程研究所 | 深空探测器推力矢量实时修正方法和系统 |
CN112520070B (zh) * | 2020-12-07 | 2022-03-29 | 上海卫星工程研究所 | 深空探测器推力矢量实时修正方法和系统 |
Also Published As
Publication number | Publication date |
---|---|
CN108061660B (zh) | 2019-09-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108061660B (zh) | 基于线振动测量的卫星发动机在轨推力实时标定方法 | |
Boyce et al. | SCRAMSPACE: scramjet-based access-to-space systems | |
Flandro et al. | Basic aerodynamics: incompressible flow | |
Lee et al. | Numerical prediction of scientific balloon trajectories while considering various uncertainties | |
Lugo et al. | Launch vehicle ascent trajectory simulation using the Program to Optimize Simulated Trajectories II (POST2) | |
Kreider | Mathematical modeling of high altitude balloon performance | |
Bardera et al. | Experimental and numerical characterization of the flow around the Mars 2020 rover | |
Chowdhury et al. | An integrated six degree-of-freedom trajectory simulator for hybrid sounding rockets | |
Cellier | Bond graphs: The right choice for educating students in modeling continuous-time physical systems | |
Carlson | Assessment of Using Ideal Gas for Predicting Boattail Flow at Cryogenic Temperatures | |
Kleb et al. | CFD: A Castle in the Sand? | |
CN103941593A (zh) | 低轨卫星姿态仿真方法 | |
Aboelezz et al. | Integrated online wind tunnel experiments and assessment for contingency scenario: case study at the british university in egypt (BUE) | |
CN108760365A (zh) | 一种软着陆试验中探测器受力状态模拟方法 | |
Wercinski et al. | Aerocapture Technology Development Needs for Outer Planet Exploration | |
Jiao et al. | Predictor–corrector method for scramjet inlet air mass flow rate measurement | |
Parish et al. | Rapid prototyping using the LabVIEW environment | |
Kharitonov et al. | Investigation of Aerogasdynamics of Reentry Vehicles in a New Hypersonic Wind Tunnel at ITAM | |
Costa et al. | Atmospheric reentry modeling and simulation | |
Dooling | An eight-parameter function for simulating model rocket engine thrust curves | |
Pinier | New Aerodynamic Data Dispersion Method with Application to Launch Vehicle Design | |
Cassaro et al. | Overview of the multipurpose aircraft simulation laboratory experience | |
Mehta | Air-breathing aerospace plane development essential: hypersonic propulsion flight tests | |
Desautel et al. | Development and integration of modern laboratories in aerospace education | |
LEROY | Huygens Close-out activities Final technical report |
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 |