CN114861320A - 一种航天器姿控推力建模及定轨解算方法 - Google Patents
一种航天器姿控推力建模及定轨解算方法 Download PDFInfo
- Publication number
- CN114861320A CN114861320A CN202210554091.8A CN202210554091A CN114861320A CN 114861320 A CN114861320 A CN 114861320A CN 202210554091 A CN202210554091 A CN 202210554091A CN 114861320 A CN114861320 A CN 114861320A
- Authority
- CN
- China
- Prior art keywords
- spacecraft
- attitude control
- control thrust
- acceleration
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 53
- 230000001133 acceleration Effects 0.000 claims abstract description 49
- 238000004364 calculation method Methods 0.000 claims abstract description 18
- 230000008569 process Effects 0.000 claims abstract description 6
- 239000011159 matrix material Substances 0.000 claims description 23
- 230000006872 improvement Effects 0.000 claims description 15
- 230000007704 transition Effects 0.000 claims description 14
- 230000010354 integration Effects 0.000 claims description 10
- 230000008859 change Effects 0.000 claims description 5
- 230000009466 transformation Effects 0.000 claims description 5
- 238000009434 installation Methods 0.000 claims description 4
- 230000000694 effects Effects 0.000 abstract description 4
- 238000005259 measurement Methods 0.000 abstract description 2
- 238000000053 physical method Methods 0.000 abstract description 2
- 229910000831 Steel Inorganic materials 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000006386 neutralization reaction Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Automation & Control Theory (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明公开了一种航天器姿控推力建模及定轨解算方法,涉及航天测控领域。该方法根据航天器姿控发动机的相关参数,基于物理方法建立理论姿控推力模型,并将本体系姿控推力加速度作为待估参数进行定轨解算,解算过程中依据理论姿控推力模型添加合理的先验约束,提高航天器轨道及姿控推力加速度的解算精度。因此,本发明能够有效降低解算姿控推力加速度与其他解算参数之间的相关性,提高轨道解算结果的稳健性,特别在处理连续、不同方向的脉冲加速度时具有较好的解算效果,可进一步提高航天器轨道确定和预报的精度。
Description
技术领域
本发明涉及航天测控领域,具体涉及一种航天器姿控推力建模及定轨解算方法。
背景技术
航天器在轨飞行期间,为了维持指定的姿态需要提供额外的力矩以抵消环境力矩。常规的姿态控制方法包括发动机喷气模式和角动量输出模式。发动机喷气模式通常为力偶的,但由于存在误差,实际工作状态下仍会对航天器产生微小的加速度,从而影响航天器的飞行轨迹。角动量输出模式通常动量轮或控制力矩陀螺,当动量轮转速饱和或外力矩超出角动量输出装置能力时,必须依靠发动机喷气来辅助控制姿态,该喷气控制往往并非力偶的,同样会对航天器产生一定的加速度。
随着我国航天事业的发展,越来越多的航天器被送入太空,精细化、复杂化的航天任务实施对轨道力学模型精度的要求越来越高。然而,受复杂环境力矩、自身气体释放、发动机安装位置偏差、脉冲时间较短导致的发动机推力变化等诸多因素影响,依靠理论模型计算姿控推力加速度存在一定的系统误差和随机误差。特别对于深空探测任务,姿控推力模型误差已经成为制约轨道计算精度的主要误差源。
姿控发动机喷气一般为脉冲式的,产生速度增量往往较小,在实际的工程任务中,通常采用经验加速度模型进行吸收。但是,频繁的、随机的脉冲喷气给经验加速度解算带来一定挑战(解算值间存在较强的相关性),往往需要将不同时刻、不同方向的姿控推力加速度进行平均处理,导致解算结果与实际情况存在较大偏差。因此,亟需建立适用性更强的姿控推力理论模型及相关参数的解算方法,以满足未来载人登月、小行星探测等复杂艰巨航天任务对高精度轨道计算的要求。
发明内容
针对上述问题,本发明的目的在于:提供一种航天器姿控推力建模及定轨解算方法,能够降低解算姿控推力加速度与其他解算参数之间的相关性,提高轨道解算结果的稳健性,特别在处理连续、不同方向的脉冲加速度时具有较好的解算效果,可进一步提高航天器轨道确定和预报的精度。
为实现上述发明目的,本发明提供以下技术方案:
一种航天器姿控推力建模及定轨解算方法,其包括以下步骤:
S1:根据航天器的质量、发动机的安装位置、额定推力和开机时长,计算姿控推力在本体系的加速度,利用姿控推力在本体系中加速度的先验分布信息构造伪观测方程;
S2:在惯性系中计算姿控推力加速度,及其关于航天器状态参数的偏导数;
S3:基于航天器初始状态参数,采用数值积分的方法计算航天器星历和状态转移矩阵;
S4:根据航天器的星历计算其理论观测值,结合实际观测值和状态转移矩阵构造观测方程;
S5:利用最小二乘方法统一求解观测方程和伪观测方程,得到航天器初始状态参数的改进量,采用迭代改进方法,获取最终的初始状态参数。
根据一种具体的实施方式,本发明的航天器姿控推力建模及定轨解算方法的步骤S1中,首先构造航天器的状态参数,定义航天器在t时刻的状态参数为:
根据一种具体的实施方式,本发明的航天器姿控推力建模及定轨解算方法的步骤S1中,在航天器本体系中计算姿控发动机脉冲喷气产生的速度增量为:其中,为发动机推力,Δt脉冲喷气时长,m为航天器质量;速度增量除以时间即为加速度,则姿控推力在本体系的加速度为:
根据一种具体的实施方式,本发明的航天器姿控推力建模及定轨解算方法的步骤S1中,由于不同姿控发动机在本体系中产生的加速度方向相对固定,更有利于参数求解,因此将其作为解算模型参数。利用姿控推力在本体系中加速度的先验分布信息,构造解算参数的伪观测方程作为先验约束,包括:假设姿控推力加速度服从正态分布:
ai~N(μi,σi),i=x,y,z
其中,μi为均值,利用姿控发动机的额定推力计算;σi为标准差,取μi的10~30%;
根据上述先验信息构造关于本体系姿控推力加速度改进变量的伪观测方程为:
根据一种具体的实施方式,本发明的航天器姿控推力建模及定轨解算方法的步骤S2中,根据航天器姿态变换矩阵,计算姿控推力在惯性系中的加速度为:其中,M为本体系到惯性系的姿态变换矩阵;以及,计算姿控推力在惯性系中加速度关于航天器状态参数的偏导数:
根据一种具体的实施方式,本发明的航天器姿控推力建模及定轨解算方法的步骤S3中,构造航天器的运动方程和轨道变分方程:
对于m维观测量,则有观测方程为:
根据一种具体的实施方式,本发明的航天器姿控推力建模及定轨解算方法的
根据一种具体的实施方式,本发明的航天器姿控推力建模及定轨解算方法的步骤S5中,将伪观测方程和观测方程联立组成新的观测方程,
其中,
考虑观测量的权重信息,即两边都乘以权重系数矩阵,利用最小二乘方法对其进行求解得:
在迭代过程中,当连续两次损耗函数的相对改变量小于设定门限时,迭代收敛,获取的航天器状态参数即为最优解,从而确定最终的初始状态参数。
附图说明
图1,航天器姿控推力建模及定轨解算流程图。
具体实施方式
下面结合附图,对本发明的具体实施方式进行详细描述,但应当理解奔放的保护范围不受具体方式的限制。
如图1所示,本发明航天器姿控推力建模及定轨解算方法,其包括以下步骤:
S1:根据航天器的质量、发动机的安装位置、额定推力和开机时长,计算姿控推力在本体系的加速度,并利用该先验信息构造伪观测方程;
S2:在惯性系中计算姿控推力加速度,及其关于航天器状态参数的偏导数;
S3:基于航天器初始状态参数,采用数值积分的方法计算航天器星历和状态转移矩阵;
S4:根据航天器的星历计算其理论观测值,结合实际观测值和状态转移矩阵构造观测方程;
S5:利用最小二乘方法统一求解观测方程和伪观测方程,得到航天器初始状态参数的改进量,采用迭代改进方法,获取最终的初始状态参数。
可选地,S1中所述的姿控推力加速度在本体系中计算。由于脉冲喷气引起的航天器质量变化很小,忽略其影响。根据牛顿第二定律,发动机开机期间产生的速度增量为:
为了便于轨道积分,将速度增量在解算区间ΔT内平均,ΔT为整数个积分步长,得到该区间内的平均加速度:
可选地,假设姿控推力加速度服从正态分布:
ai~N(μi,σi),i=x,y,z
其中,μi为均值,利用姿控发动机的额定推力计算;σi为标准差,取μi的10~30%。
可选地,根据上述先验信息构造关于本体系姿控推力加速度改进变量的伪观测方程为:
可选地,S2中所述的姿控推力加速度在惯性系中计算。根据航天器的姿态模式,计算本体系到惯性系的姿态变换矩阵,以三轴对地稳定姿态为例,
再将本体系姿控推力加速度转换惯性坐标系:
可选地,S3中星历和状态转移矩阵计算方法为:构造航天器的运动方程和轨道变分方程:
可选地,S3中针对姿控推力加速度解算区间的数值积分采取分段处理方法:在姿控发动机开机和关机时刻开始重新开始积分,中断节点处采取内插的方式获取星历及状态转移矩阵。
定义状态转移矩阵为:
由S3可知,
则
对于m维观测量,则有观测方程为:
其中,Λ为权重系数矩阵,wi=1/σi(i=1,2,…,m)为观测量的权重系数。
可选地,步骤S5中,将伪观测方程和观测方程联立组成新的观测方程,
其中,
考虑观测量的权重信息,即两边都乘以权重系数矩阵W12,利用最小二乘方法对其进行求解得:
根据一种具体的实施方式,本发明的航天器姿控推力建模及定轨解算方法的步骤S5中,由于对观测模型的线性化将引入误差,单次求解的存在误差,因此采用迭代改进方法对初始状态参数进行修正:设第n次最小二乘求解得到的改进量为将其叠加到轨道初值上,作为下一次轨道改进的初值
在迭代过程中,当连续两次损耗函数的相对改变量小于设定门限时,迭代收敛,获取的航天器状态参数即为最优解,从而确定最终的初始状态参数。
本发明航天器姿控推力建模及定轨解算方法,根据航天器姿控发动机的相关参数,基于物理方法建立理论姿控推力模型,并将本体系姿控推力加速度作为待估参数进行定轨解算,解算过程中依据理论姿控推力模型添加合理的先验约束,提高航天器轨道及姿控推力加速度的解算精度。因此,本发明能够有效降低解算姿控推力加速度与其他解算参数之间的相关性,提高轨道解算结果的稳健性,特别在处理连续、不同方向的脉冲加速度时具有较好的解算效果,可进一步提高航天器轨道确定和预报的精度。
Claims (7)
1.一种航天器姿控推力建模及定轨解算方法,其特征在于,包括以下步骤:
S1:根据航天器的质量、发动机的安装位置、额定推力和开机时长,计算姿控推力在本体系的加速度,并利用姿控推力在本体系中加速度的先验分布信息构造伪观测方程;
S2:在惯性系中计算姿控推力加速度,及其关于航天器状态参数的偏导数;
S3:基于航天器初始状态参数,采用数值积分的方法计算航天器星历和状态转移矩阵;
S4:根据航天器的星历计算其理论观测值,结合实际观测值和状态转移矩阵构造观测方程;
S5:利用最小二乘方法统一求解观测方程和伪观测方程,得到航天器初始状态参数的改进量,采用迭代改进方法,获取最终的初始状态参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210554091.8A CN114861320B (zh) | 2022-05-19 | 2022-05-19 | 一种航天器姿控推力建模及定轨解算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210554091.8A CN114861320B (zh) | 2022-05-19 | 2022-05-19 | 一种航天器姿控推力建模及定轨解算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114861320A true CN114861320A (zh) | 2022-08-05 |
CN114861320B CN114861320B (zh) | 2023-02-10 |
Family
ID=82639865
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210554091.8A Active CN114861320B (zh) | 2022-05-19 | 2022-05-19 | 一种航天器姿控推力建模及定轨解算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114861320B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115906523A (zh) * | 2022-12-28 | 2023-04-04 | 北京航天飞行控制中心 | 用于轨道计算的待估参数优化方法及其装置 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010266468A (ja) * | 2010-08-30 | 2010-11-25 | Seiko Epson Corp | 位置算出方法及び位置算出装置 |
WO2012021760A2 (en) * | 2010-08-12 | 2012-02-16 | The Government Of The United States Of America As Represented By The Secretary Of The Navy | Improved orbit covariance estimation and analysis (ocean) system and method |
US20120084068A1 (en) * | 2010-09-30 | 2012-04-05 | University Of Utah Research Foundation | Incorporation of Uncertainty in Simulation Analysis |
CN106508029B (zh) * | 2010-12-16 | 2013-11-13 | 中国人民解放军国防科学技术大学 | 一种航天器动力学参数综合辨识方法 |
CN105867167A (zh) * | 2016-04-13 | 2016-08-17 | 中国人民解放军63680部队 | 跨平台的自旋稳定卫星飞行仿真系统 |
WO2017005052A1 (zh) * | 2015-07-09 | 2017-01-12 | 北京航空航天大学 | 航天器脉冲交会轨迹的梯度分割区间优化设计方法 |
CN113375677A (zh) * | 2021-08-12 | 2021-09-10 | 中国人民解放军国防科技大学 | 一种基于脉冲星观测的航天器定速方法 |
CN114063122A (zh) * | 2021-12-02 | 2022-02-18 | 武汉大学 | 电推进转移轨道航天器星载gnss在轨实时定轨方法 |
EP3960638A1 (en) * | 2020-08-24 | 2022-03-02 | Beijing Interstellar Glory Space Technology Co., Ltd. | Energy management method for aircraft, control method and aircraft |
CN114239131A (zh) * | 2021-11-24 | 2022-03-25 | 沈阳航天新光集团有限公司 | 姿控动力系统空间飞行力学环境预示方法 |
-
2022
- 2022-05-19 CN CN202210554091.8A patent/CN114861320B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012021760A2 (en) * | 2010-08-12 | 2012-02-16 | The Government Of The United States Of America As Represented By The Secretary Of The Navy | Improved orbit covariance estimation and analysis (ocean) system and method |
JP2010266468A (ja) * | 2010-08-30 | 2010-11-25 | Seiko Epson Corp | 位置算出方法及び位置算出装置 |
US20120084068A1 (en) * | 2010-09-30 | 2012-04-05 | University Of Utah Research Foundation | Incorporation of Uncertainty in Simulation Analysis |
CN106508029B (zh) * | 2010-12-16 | 2013-11-13 | 中国人民解放军国防科学技术大学 | 一种航天器动力学参数综合辨识方法 |
WO2017005052A1 (zh) * | 2015-07-09 | 2017-01-12 | 北京航空航天大学 | 航天器脉冲交会轨迹的梯度分割区间优化设计方法 |
CN105867167A (zh) * | 2016-04-13 | 2016-08-17 | 中国人民解放军63680部队 | 跨平台的自旋稳定卫星飞行仿真系统 |
EP3960638A1 (en) * | 2020-08-24 | 2022-03-02 | Beijing Interstellar Glory Space Technology Co., Ltd. | Energy management method for aircraft, control method and aircraft |
CN113375677A (zh) * | 2021-08-12 | 2021-09-10 | 中国人民解放军国防科技大学 | 一种基于脉冲星观测的航天器定速方法 |
CN114239131A (zh) * | 2021-11-24 | 2022-03-25 | 沈阳航天新光集团有限公司 | 姿控动力系统空间飞行力学环境预示方法 |
CN114063122A (zh) * | 2021-12-02 | 2022-02-18 | 武汉大学 | 电推进转移轨道航天器星载gnss在轨实时定轨方法 |
Non-Patent Citations (3)
Title |
---|
NI Z等: "On-orbit identification of time-varying moment of inertia for spacecraft based on a recursive subspace method", 《CONTROL CONFERENCE》 * |
张宇等: "长期姿控扰动情况下空间实验室轨道影响分析及建模", 《宇航学报》 * |
杨天社等: "航天器测控仿真方法研究", 《飞行器测控学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115906523A (zh) * | 2022-12-28 | 2023-04-04 | 北京航天飞行控制中心 | 用于轨道计算的待估参数优化方法及其装置 |
CN115906523B (zh) * | 2022-12-28 | 2024-04-30 | 北京航天飞行控制中心 | 用于轨道计算的待估参数优化方法及其装置 |
Also Published As
Publication number | Publication date |
---|---|
CN114861320B (zh) | 2023-02-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Olds et al. | IRVE-3 post-flight reconstruction | |
Searcy et al. | Magnetometer-only attitude determination using novel two-step Kalman filter approach | |
Steltzner et al. | Mars Science Laboratory entry, descent, and landing system | |
CN105573337B (zh) | 一种满足再入角和航程约束的离轨制动闭路制导方法 | |
CN104764467B (zh) | 空天飞行器惯性传感器误差在线自适应标定方法 | |
Woffinden et al. | Linear covariance analysis techniques to generate navigation and sensor requirements for the safe and precise landing integrated capabilities evolution (SPLICE) project | |
Geller et al. | Linear covariance analysis for powered lunar descent and landing | |
CN102819266B (zh) | 一种拟周期j2不变相对轨道编队飞行控制方法 | |
CN107792396B (zh) | 发动机工作时干扰力矩实时估计方法 | |
Abdelrahman et al. | Sigma-point Kalman filtering for spacecraft attitude and rate estimation using magnetometer measurements | |
Ma et al. | Magnetometer-only attitude and angular velocity filtering estimation for attitude changing spacecraft | |
CN114861320B (zh) | 一种航天器姿控推力建模及定轨解算方法 | |
CN114253291B (zh) | 基于线性伪谱模型预测控制的航天器编队制导方法和系统 | |
CN108516107B (zh) | 一种在线发动机推力和比冲估计方法及系统 | |
CN111638643B (zh) | 一种位移模式无拖曳控制动力学协调条件确定方法 | |
Li et al. | Adaptive sliding mode control for spacecraft rendezvous with unknown system parameters and input saturation | |
CN112817233B (zh) | 一种基于迭代学习控制的小天体探测器绕飞段轨道跟踪控制方法 | |
Somov et al. | Health checking autonomous attitude control system of Earth-observing miniature satellite in initial orientation modes | |
Challa et al. | Simultaneous determination of spacecraft attitude and rates using only a magnetometer | |
Vetrisano et al. | Navigating to the Moon along low-energy transfers | |
CN106326576B (zh) | 一种任意基准系下的整星偏置角动量的偏航估计方法 | |
CN110955255A (zh) | 基于cmg的高精度轨控姿态维持方法、系统及介质 | |
Yu et al. | Approximate analytical solutions for launch-vehicle ascent trajectory | |
Crain et al. | Mars entry navigation: atmospheric interface through parachute deploy | |
Lapin et al. | Autonomous stabilization of a spacecraft orbital orientation at the lack of angular velocity measurements |
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 |