CN113779788A - 一种大入轨偏差情况下组合体分离轨道确定方法和系统 - Google Patents

一种大入轨偏差情况下组合体分离轨道确定方法和系统 Download PDF

Info

Publication number
CN113779788A
CN113779788A CN202111028656.0A CN202111028656A CN113779788A CN 113779788 A CN113779788 A CN 113779788A CN 202111028656 A CN202111028656 A CN 202111028656A CN 113779788 A CN113779788 A CN 113779788A
Authority
CN
China
Prior art keywords
orbit
separation
determining
combined
deviation
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
Application number
CN202111028656.0A
Other languages
English (en)
Other versions
CN113779788B (zh
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 CN202111028656.0A priority Critical patent/CN113779788B/zh
Publication of CN113779788A publication Critical patent/CN113779788A/zh
Application granted granted Critical
Publication of CN113779788B publication Critical patent/CN113779788B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • 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
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Automation & Control Theory (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Navigation (AREA)

Abstract

本发明提供了一种大入轨偏差情况下组合体分离轨道确定方法和系统，所述方法包括如下步骤：步骤S1：根据星箭分离异常状态，确定大偏差入轨点轨道参数及组合体推进舱与轨道转移相关的平台参数；步骤S2：建立描述航天器轨道转移运动学模型；步骤S3：利用交轨共点方法对分离轨道点火点位置、速度增量及燃料消耗进行计算；步骤S4：以期望参数调整至目标值为原则，组合体推进舱最大可用燃料量为边界，积分迭代实现对分离轨道参数的确定。本发明解决了发生较大入轨偏差，纠偏燃料消耗导致原有转移任务无法完成情况下，组合体航天器最优分离轨道确定方法，为异常状况下航天任务补救提供支撑，具有一定的工程实用性。

Description

一种大入轨偏差情况下组合体分离轨道确定方法和系统
技术领域
本发明涉及宇航飞行器轨道运动学的技术领域，具体地，涉及一种大入轨偏差情况下组合体分离轨道确定方法和系统。
背景技术
组合体航天器通常包含多个舱段，并且每个舱段都携带推进系统，这也使得发生较大入轨偏差，纠偏燃料消耗导致原有轨道转移任务无法完成情况下，可以通过推进舱与任务舱段提前分离实现异常状况下航天任务抢救。
国内外目前针对组合体航天器的研究多是集中于模型预测与稳定控制领域，如名为“基于深度学习的组合体航天器模型预测控制”期刊论文中借助深度学习在多参数寻优上的优势，提出一种基于卷积神经网络的模型预测控制算法，实现组合体航天器多场景下姿态控制率重构；名为“组合体航天器的姿态无模型自适应控制”期刊论文中提出了一种基于无模型自适应控制方法，解决了转动惯量参数未知组合体航天器姿态精确控制问题；公开号为CN110789738A的中国专利公开了一种纳星-失效航天器组合体组态运动模型与控制模型，设计了分布式模型预测控制器。面向大入轨偏差的组合体航天器最优分离轨道参数确定方法未见公开的研究成果。
综上所述，需要针对大入轨偏差情况下组合体最优分离轨道确定方法进行设计优化。
发明内容
针对现有技术中的缺陷，本发明的目的是提供一种大入轨偏差情况下组合体分离轨道确定方法和系统。
根据本发明提供的一种大入轨偏差情况下组合体分离轨道确定方法，所述方法包括如下步骤：
步骤S1：根据星箭分离异常状态，确定大偏差入轨点轨道参数及组合体推进舱与轨道转移相关的平台参数；
步骤S2：建立描述航天器轨道转移运动学模型；
步骤S3：利用交轨共点方法对分离轨道点火点位置、速度增量及燃料消耗进行计算；
步骤S4：以期望参数调整至目标值为原则，组合体推进舱最大可用燃料量为边界，积分迭代实现对分离轨道参数的确定。
优选地，所述步骤S1包括如下步骤：
步骤S1.1：组合体航天器包含多个舱段，每个舱段携带推进系统，大入轨偏差情况下推进舱与任务舱需要提前分离，分离轨道的确定首先需要明确异常入轨时刻的组合体航天器轨道状态，即采用轨道瞬时六根数形式进行描述，即半长轴a0、偏心率e0、倾角i0、升交点赤经Ω0、近地点幅角ω0和平近点角M0；
步骤S1.2：明确组合体推进舱相关平台参数，包括星箭分离异常时刻重量m、最大可用燃料剩余量m′、发动机推力F和发动机比冲Isp。
优选地，所述步骤S2中建立J2地球引力摄动影响下的轨道转移运动学模型：
Figure BDA0003244292330000021
其中，
Figure BDA0003244292330000022
分别为点火时刻组合体航天器地心惯性系位置、速度矢量；
r为点火时刻组合体航天器位置矢量
Figure BDA0003244292330000023
的模值；
(x,y,z)为点火时刻航天器地心惯性系三轴位置分量；
Figure BDA0003244292330000024
为点火时刻组合体航天器地心惯性系三轴速度分量；
(α,δ)为组合体航天器推力方向的赤经、赤纬；
μ＝3.986005×1014m3·s-2为地球引力常数；
Re＝6378.14×103m为地球赤道半径；
g0＝9.80665m·s-2为地球重力加速度。
优选地，所述步骤S3包括如下步骤：
步骤S3.1：转移轨道点火点位置计算，目标升交点赤经Ωs、目标倾角is，计算得到目标轨道平面正法向矢量
Figure BDA0003244292330000031
为：
Figure BDA0003244292330000032
在星箭分离轨道上积分寻找t0时刻，确定点火点的位置速度
Figure BDA0003244292330000033
Figure BDA0003244292330000034
步骤S3.2：转移轨道速度增量计算，点火点P0的位置速度
Figure BDA0003244292330000035
目标速度大小为
Figure BDA0003244292330000036
依照如下方向构建相对坐标系Fxyz：
Figure BDA0003244292330000037
目标速度与相对坐标系Fxyz的轴向夹角为θ，得到目标速度在惯性坐标系下的表示
Figure BDA0003244292330000038
Figure BDA0003244292330000039
Fxyz'为构建的相对坐标系，以夹角θ为变量，
Figure BDA00032442923300000310
取极小为目标，[0°,360°]遍历搜索确定轨道转移速度增量Δv：
Figure BDA00032442923300000311
步骤S3.3：轨道转移耗肼量计算，计算得到轨道转移最小速度增量对应的燃料消耗Δm：
Figure BDA00032442923300000312
优选地，所述步骤S4包括如下步骤：
步骤S4.1：确定大入轨偏差条件下的期望调整参数序列Orb′，序列包含半长轴a′、倾角i′与升交点赤经Ω′中某一项或两项参数；
步骤S4.2：给定大入轨偏差条件下的分离轨道参数序列Orb，分离轨道参数序列包含期望调整参数与自由参数；
步骤S4.3：根据步骤S3的计算过程，以推进舱可用燃料消耗殆尽为目标，积分迭代确定分离轨道参数序列中的自由参数，确定输出分离轨道参数。
本发明还提供一种大入轨偏差情况下组合体分离轨道确定系统，所述系统包括如下模块：
模块M1：根据星箭分离异常状态，确定大偏差入轨点轨道参数及组合体推进舱与轨道转移相关的平台参数；
模块M2：建立描述航天器轨道转移运动学模型；
模块M3：利用交轨共点方法对分离轨道点火点位置、速度增量及燃料消耗进行计算；
模块M4：以期望参数调整至目标值为原则，组合体推进舱最大可用燃料量为边界，积分迭代实现对分离轨道参数的确定。
优选地，所述模块M1包括如下模块：
模块M1.1：组合体航天器包含多个舱段，每个舱段携带推进系统，大入轨偏差情况下推进舱与任务舱需要提前分离，分离轨道的确定首先需要明确异常入轨时刻的组合体航天器轨道状态，即采用轨道瞬时六根数形式进行描述，即半长轴a0、偏心率e0、倾角i0、升交点赤经Ω0、近地点幅角ω0和平近点角M0；
模块M1.2：明确组合体推进舱相关平台参数，包括星箭分离异常时刻重量m、最大可用燃料剩余量m′、发动机推力F和发动机比冲Isp。
优选地，所述模块M2中建立J2地球引力摄动影响下的轨道转移运动学模型：
Figure BDA0003244292330000041
其中，
Figure BDA0003244292330000042
分别为点火时刻组合体航天器地心惯性系位置、速度矢量；
r为点火时刻组合体航天器位置矢量
Figure BDA0003244292330000043
的模值；
(x,y,z)为点火时刻航天器地心惯性系三轴位置分量；
Figure BDA0003244292330000051
为点火时刻组合体航天器地心惯性系三轴速度分量；
(α,δ)为组合体航天器推力方向的赤经、赤纬；
μ＝3.986005×1014m3·s-2为地球引力常数；
Re＝6378.14×103m为地球赤道半径；
g0＝9.80665m·s-2为地球重力加速度。
优选地，所述模块M3包括如下模块：
模块M3.1：转移轨道点火点位置计算，目标升交点赤经Ωs、目标倾角is，计算得到目标轨道平面正法向矢量
Figure BDA0003244292330000052
为：
Figure BDA0003244292330000053
在星箭分离轨道上积分寻找t0时刻，确定点火点的位置速度
Figure BDA0003244292330000054
Figure BDA0003244292330000055
模块M3.2：转移轨道速度增量计算，点火点P0的位置速度
Figure BDA0003244292330000056
目标速度大小为
Figure BDA0003244292330000057
依照如下方向构建相对坐标系Fxyz：
Figure BDA0003244292330000058
目标速度与相对坐标系Fxyz的轴向夹角为θ，得到目标速度在惯性坐标系下的表示
Figure BDA0003244292330000059
Figure BDA00032442923300000510
以夹角θ为变量，
Figure BDA00032442923300000511
取极小为目标，[0°,360°]遍历搜索确定轨道转移速度增量Δv：
Figure BDA00032442923300000512
模块M3.3：轨道转移耗肼量计算，计算得到轨道转移最小速度增量对应的燃料消耗Δm：
Figure BDA0003244292330000061
优选地，所述模块M4包括如下模块：
模块M4.1：确定大入轨偏差条件下的期望调整参数序列Orb′，序列包含半长轴a′、倾角i′与升交点赤经Ω′中某一项或两项参数；
模块M4.2：给定大入轨偏差条件下的分离轨道参数序列Orb，分离轨道参数序列包含期望调整参数与自由参数；
模块M4.3：根据模块M3的计算过程，以推进舱可用燃料消耗殆尽为目标，积分迭代确定分离轨道参数序列中的自由参数，确定输出分离轨道参数。
与现有技术相比，本发明具有如下的有益效果：
本发明解决了发生较大入轨偏差，纠偏燃料消耗导致原有转移任务无法完成情况下，组合体航天器最优分离轨道确定方法，为异常状况下航天任务补救提供支撑，具有一定的工程实用性。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述，本发明的其它特征、目的和优点将会变得更明显：
图1为本发明的原理框图；
图2为本发明的组合体航天器(两舱段)构型示意图；
图3本发明的组合体航天器大入轨偏差轨道示意图；
图4本发明的仿真图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明，但不以任何形式限制本发明。应当指出的是，对本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变化和改进。这些都属于本发明的保护范围。
参照图1和图2，本发明提供一种大入轨偏差情况下组合体分离轨道确定方法，包括如下步骤：
步骤S1：根据星箭分离异常状态，确定大偏差入轨点轨道参数及组合体推进仓与轨道转移相关的平台参数。
步骤S1.1：组合体航天器通常包含多个舱段，并且每个舱段都携带推进系统，这也使得发生较大入轨偏差，纠偏燃料消耗导致原有转移任务无法完成情况下，可以通过推进舱与任务舱段提前分离实现异常状况下航天任务抢救。为了确定组合体最优分离轨道参数，首先需要明确异常星箭分离时刻的轨道状态。不失一般性，采用轨道瞬时六根数形式进行描述，即半长轴a0、偏心率e0、倾角i0、升交点赤经Ω0、近地点幅角ω0和平近点角M0。
步骤S1.2：此外，还需要明确组合体推进仓相关平台参数，一般包括星箭分离异常时刻重量m、最大可用燃料剩余量m′、发动机推力F、发动机比冲Isp。
步骤S2：建立描述航天器轨道转移运动学模型。
步骤S2中建立J2地球引力摄动影响下的轨道转移运动学模型：
Figure BDA0003244292330000071
其中，
Figure BDA0003244292330000072
分别为点火时刻组合体航天器地心惯性系位置、速度矢量；
r为点火时刻组合体航天器位置矢量
Figure BDA0003244292330000073
的模值；
(x,y,z)为点火时刻航天器地心惯性系三轴位置分量；
Figure BDA0003244292330000074
为点火时刻组合体航天器地心惯性系三轴速度分量；
(α,δ)为组合体航天器推力方向的赤经、赤纬；
μ＝3.986005×1014m3·s-2为地球引力常数；
Re＝6378.14×103m为地球赤道半径；
g0＝9.80665m·s-2为地球重力加速度。
步骤S3：利用交轨共点方法对虚拟分离轨道点火点位置、速度增量及燃料消耗进行估计。
步骤S3.1：转移轨道点火点位置计算，即：
假设目标升交点赤经Ωs、目标倾角is，则可以计算得到目标轨道平面正法向矢量
Figure BDA0003244292330000081
为：
Figure BDA0003244292330000082
通过在星箭分离轨道上积分寻找合适的t0时刻，使得如下条件满足，可以相应确定点火点P0的位置速度
Figure BDA0003244292330000083
Figure BDA0003244292330000084
步骤S3.2：转移轨道速度增量计算，即：
已知点火点P0的位置速度
Figure BDA0003244292330000085
假设由活力公式得到的目标速度大小为
Figure BDA0003244292330000086
则依照如下方向构建相对坐标系Fxyz：
Figure BDA0003244292330000087
假设目标速度与相对坐标系Fxyz的x轴向夹角为θ，则可以得到目标速度在惯性坐标系下的表示
Figure BDA0003244292330000088
Figure BDA0003244292330000089
以夹角θ为变量，
Figure BDA00032442923300000810
取极小为目标，[0°,360°]遍历搜索即可确定轨道转移速度增量Δv：
Figure BDA00032442923300000811
步骤S3.3：轨道转移耗肼量计算，即：
根据火箭公式，计算得到轨道转移最小速度增量对应的燃料消耗Δm：
Figure BDA00032442923300000812
步骤S4：以期望参数优先调整至目标值为原则，组合体推进舱最大可用燃料量为边界，积分迭代实现对分离轨道参数寻优确定。
步骤S4.1：确定大入轨偏差条件下的期望优先调整参数序列Orb′，不失一般性通常序列包含半长轴a′、倾角i′与升交点赤经Ω′中某一项或两项参数；
步骤S4.2：给定大入轨偏差条件下的分离轨道参数序列Orb，该序列需要包含期望优先调整参数与自由参数，如Orb＝(a,i′,Ω′)；
步骤S4.3：根据步骤S3计算过程，以推进仓可用燃料消耗殆尽为目标，积分迭代确定分离轨道参数序列中的自由参数，从而确定输出最优分离轨道参数。
参照图3，本实施例中，假设由于运载器故障导至入轨点参数发生较大偏，依照步骤S1确定分离轨道参数：
Orb0＝(a0,i0,e0,Ω0，ω0,M0)＝(20863.14km,30°,0.6861,353.76°,179.3°,1.39°)
确定组合体推进仓与轨道转移相关的平台参数为：入轨时刻重量m0＝3900kg、可用燃料剩余量m'0＝1780kg、发动机推力F0＝350N、发动机比冲Isp＝315s。
确定期望优先调整参数序列，假设Orb'＝(i′＝0°,Ω′＝264°)。则以推进仓燃料可用消耗量m’0为可达边界，以i′＝0°,Ω′＝264°为优先确定的目标轨道参数，以轨道半场轴a为寻优参数，利用步骤S3计算方法，结合步骤S2的转移轨道运动学模型，逐次计算不同目标轨道Orb＝(a，i′，Ω′)的燃料消耗，最终利用步骤S4所述的积分迭代程序确定最优分离轨道参数为：
Orbn＝(an，in,en,Ωn,ωn,Mn)＝(32971.07km,0°,0.06414,264.31°，257.40°,193.37°)
参照图4，由确定的分离轨道参数易见：分离轨道倾角in＝0°、升交点赤经Ωn＝264.31°满足期望参数调整目标值；轨道转移燃料消耗Δm＝1779.99kg，满足推进仓燃料可达边界。由此可见该分离轨道可以为组合体分离后的其他舱段提供最优的变轨转移条件。
本发明还提供一种大入轨偏差情况下组合体分离轨道确定系统，所述系统包括如下模块：
模块M1：根据星箭分离异常状态，确定大偏差入轨点轨道参数及组合体推进舱与轨道转移相关的平台参数；模块M1.1：组合体航天器包含多个舱段，每个舱段携带推进系统，大入轨偏差情况下推进舱与任务舱需要提前分离，分离轨道的确定首先需要明确异常入轨时刻的组合体航天器轨道状态，即采用轨道瞬时六根数形式进行描述，即半长轴a0、偏心率e0、倾角i0、升交点赤经Ω0、近地点幅角ω0和平近点角M0；模块M1.2：明确组合体推进舱相关平台参数，包括星箭分离异常时刻重量m、最大可用燃料剩余量m′、发动机推力F和发动机比冲Isp。
模块M2：建立描述航天器轨道转移运动学模型；所述模块M2中建立J2地球引力摄动影响下的轨道转移运动学模型：
Figure BDA0003244292330000101
其中，
Figure BDA0003244292330000102
分别为点火时刻组合体航天器地心惯性系位置、速度矢量；
r为点火时刻组合体航天器位置矢量
Figure BDA0003244292330000103
的模值；
(x,y,z)为点火时刻航天器地心惯性系三轴位置分量；
Figure BDA0003244292330000104
为点火时刻组合体航天器地心惯性系三轴速度分量；
(α,δ)为组合体航天器推力方向的赤经、赤纬；
μ＝3.986005×1014m3·s-2为地球引力常数；
Re＝6378.14×103m为地球赤道半径；
g0＝9.80665m·s-2为地球重力加速度。
模块M3：利用交轨共点方法对分离轨道点火点位置、速度增量及燃料消耗进行计算；模块M3.1：转移轨道点火点位置计算，目标升交点赤经Ωs、目标倾角is，计算得到目标轨道平面正法向矢量
Figure BDA0003244292330000105
为：
Figure BDA0003244292330000106
在星箭分离轨道上积分寻找t0时刻，确定点火点的位置速度
Figure BDA0003244292330000107
Figure BDA0003244292330000111
模块M3.2：转移轨道速度增量计算，点火点P0的位置速度
Figure BDA0003244292330000112
目标速度大小为
Figure BDA0003244292330000113
依照如下方向构建相对坐标系Fxyz：
Figure BDA0003244292330000114
目标速度与相对坐标系Fxyz的轴向夹角为θ，得到目标速度在惯性坐标系下的表示
Figure BDA0003244292330000115
Figure BDA0003244292330000116
以夹角θ为变量，
Figure BDA0003244292330000117
取极小为目标，[0°,360°]遍历搜索确定轨道转移速度增量Δv：
Figure BDA0003244292330000118
模块M3.3：轨道转移耗肼量计算，计算得到轨道转移最小速度增量对应的燃料消耗Δm：
Figure BDA0003244292330000119
模块M4：以期望参数调整至目标值为原则，组合体推进舱最大可用燃料量为边界，积分迭代实现对分离轨道参数的确定。模块M4.1：确定大入轨偏差条件下的期望调整参数序列Orb′，序列包含半长轴a′、倾角i′与升交点赤经Ω′中某一项或两项参数；模块M4.2：给定大入轨偏差条件下的分离轨道参数序列Orb，分离轨道参数序列包含期望调整参数与自由参数；模块M4.3：根据模块M3的计算过程，以推进舱可用燃料消耗殆尽为目标，积分迭代确定分离轨道参数序列中的自由参数，确定输出分离轨道参数。
本发明解决了发生较大入轨偏差，纠偏燃料消耗导致原有转移任务无法完成情况下，组合体航天器最优分离轨道确定方法，为异常状况下航天任务补救提供支撑，具有一定的工程实用性。
本领域技术人员知道，除了以纯计算机可读程序代码方式实现本发明提供的系统及其各个装置、模块、单元以外，完全可以通过将方法步骤进行逻辑编程来使得本发明提供的系统及其各个装置、模块、单元以逻辑门、开关、专用集成电路、可编程逻辑控制器以及嵌入式微控制器等的形式来实现相同功能。所以，本发明提供的系统及其各项装置、模块、单元可以被认为是一种硬件部件，而对其内包括的用于实现各种功能的装置、模块、单元也可以视为硬件部件内的结构；也可以将用于实现各种功能的装置、模块、单元视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
以上对本发明的具体实施例进行了描述。需要理解的是，本发明并不局限于上述特定实施方式，本领域技术人员可以在权利要求的范围内做出各种变化或修改，这并不影响本发明的实质内容。在不冲突的情况下，本申请的实施例和实施例中的特征可以任意相互组合。

Claims (10)

1.一种大入轨偏差情况下组合体分离轨道确定方法，其特征在于，所述方法包括如下步骤：
步骤S1：根据星箭分离异常状态，确定大偏差入轨点轨道参数及组合体推进舱与轨道转移相关的平台参数；
步骤S2：建立描述航天器轨道转移运动学模型；
步骤S3：利用交轨共点方法对分离轨道点火点位置、速度增量及燃料消耗进行计算；
步骤S4：以期望参数调整至目标值为原则，组合体推进舱最大可用燃料量为边界，积分迭代实现对分离轨道参数的确定。
2.根据权利要求1所述的一种大入轨偏差情况下组合体分离轨道确定方法，其特征在于，所述步骤S1包括如下步骤：
步骤S1.1：组合体航天器包含多个舱段，每个舱段携带推进系统，大入轨偏差情况下推进舱与任务舱需要提前分离，分离轨道的确定首先需要明确异常入轨时刻的组合体航天器轨道状态，即采用轨道瞬时六根数形式进行描述，即半长轴a0、偏心率e0、倾角i0、升交点赤经Ω0、近地点幅角ω0和平近点角M0；
步骤S1.2：明确组合体推进舱相关平台参数，包括星箭分离异常时刻重量m、最大可用燃料剩余量m′、发动机推力F和发动机比冲Isp。
3.根据权利要求1所述的一种大入轨偏差情况下组合体分离轨道确定方法，其特征在于，所述步骤S2中建立J2地球引力摄动影响下的轨道转移运动学模型：
Figure FDA0003244292320000011
其中，
Figure FDA0003244292320000012
分别为点火时刻组合体航天器地心惯性系位置、速度矢量；
r为点火时刻组合体航天器位置矢量
Figure FDA0003244292320000021
的模值；
(x,y,z)为点火时刻航天器地心惯性系三轴位置分量；
Figure FDA0003244292320000022
为点火时刻组合体航天器地心惯性系三轴速度分量；
(α,δ)为组合体航天器推力方向的赤经、赤纬；
μ＝3.986005×1014m3·s-2为地球引力常数；
Re＝6378.14×103m为地球赤道半径；
g0＝9.80665m·s-2为地球重力加速度。
4.根据权利要求1所述的一种大入轨偏差情况下组合体分离轨道确定方法，其特征在于，所述步骤S3包括如下步骤：
步骤S3.1：转移轨道点火点位置计算，目标升交点赤经Ωs、目标倾角is，计算得到目标轨道平面正法向矢量
Figure FDA0003244292320000023
为：
Figure FDA0003244292320000024
在星箭分离轨道上积分寻找t0时刻，确定点火点的位置速度
Figure FDA0003244292320000025
Figure FDA0003244292320000026
步骤S3.2：转移轨道速度增量计算，点火点P0的位置速度
Figure FDA0003244292320000027
目标速度大小为
Figure FDA0003244292320000028
依照如下方向构建相对坐标系Fxyz：
Figure FDA0003244292320000029
目标速度与相对坐标系Fxyz的轴向夹角为θ，得到目标速度在惯性坐标系下的表示
Figure FDA00032442923200000210
Figure FDA00032442923200000211
Fxyz'为构建的相对坐标系，以夹角θ为变量，
Figure FDA00032442923200000212
取极小为目标，[0°,360°]遍历搜索确定轨道转移速度增量Δv：
Figure FDA00032442923200000213
步骤S3.3：轨道转移耗肼量计算，计算得到轨道转移最小速度增量对应的燃料消耗Δm：
Figure FDA0003244292320000031
5.根据权利要求1所述的一种大入轨偏差情况下组合体分离轨道确定方法，其特征在于，所述步骤S4包括如下步骤：
步骤S4.1：确定大入轨偏差条件下的期望调整参数序列Orb′，序列包含半长轴a′、倾角i′与升交点赤经Ω′中某一项或两项参数；
步骤S4.2：给定大入轨偏差条件下的分离轨道参数序列Orb，分离轨道参数序列包含期望调整参数与自由参数；
步骤S4.3：根据步骤S3的计算过程，以推进舱可用燃料消耗殆尽为目标，积分迭代确定分离轨道参数序列中的自由参数，确定输出分离轨道参数。
6.一种大入轨偏差情况下组合体分离轨道确定系统，其特征在于，所述系统包括如下模块：
模块M1：根据星箭分离异常状态，确定大偏差入轨点轨道参数及组合体推进舱与轨道转移相关的平台参数；
模块M2：建立描述航天器轨道转移运动学模型；
模块M3：利用交轨共点方法对分离轨道点火点位置、速度增量及燃料消耗进行计算；
模块M4：以期望参数调整至目标值为原则，组合体推进舱最大可用燃料量为边界，积分迭代实现对分离轨道参数的确定。
7.根据权利要求6所述的一种大入轨偏差情况下组合体分离轨道确定系统，其特征在于，所述模块M1包括如下模块：
模块M1.1：组合体航天器包含多个舱段，每个舱段携带推进系统，大入轨偏差情况下推进舱与任务舱需要提前分离，分离轨道的确定首先需要明确异常入轨时刻的组合体航天器轨道状态，即采用轨道瞬时六根数形式进行描述，即半长轴a0、偏心率e0、倾角i0、升交点赤经Ω0、近地点幅角ω0和平近点角M0；
模块M1.2：明确组合体推进舱相关平台参数，包括星箭分离异常时刻重量m、最大可用燃料剩余量m′、发动机推力F和发动机比冲Isp。
8.根据权利要求6所述的一种大入轨偏差情况下组合体分离轨道确定系统，其特征在于，所述模块M2中建立J2地球引力摄动影响下的轨道转移运动学模型：
Figure FDA0003244292320000041
其中，
Figure FDA0003244292320000042
分别为点火时刻组合体航天器地心惯性系位置、速度矢量；
r为点火时刻组合体航天器位置矢量
Figure FDA0003244292320000043
的模值；
(x,y,z)为点火时刻航天器地心惯性系三轴位置分量；
Figure FDA0003244292320000044
为点火时刻组合体航天器地心惯性系三轴速度分量；
(α,δ)为组合体航天器推力方向的赤经、赤纬；
μ＝3.986005×1014m3·s-2为地球引力常数；
Re＝6378.14×103m为地球赤道半径；
g0＝9.80665m·s-2为地球重力加速度。
9.根据权利要求6所述的一种大入轨偏差情况下组合体分离轨道确定系统，其特征在于，所述模块M3包括如下模块：
模块M3.1：转移轨道点火点位置计算，目标升交点赤经Ωs、目标倾角is，计算得到目标轨道平面正法向矢量
Figure FDA0003244292320000045
为：
Figure FDA0003244292320000046
在星箭分离轨道上积分寻找t0时刻，确定点火点的位置速度
Figure FDA0003244292320000047
Figure FDA0003244292320000048
模块M3.2：转移轨道速度增量计算，点火点P0的位置速度
Figure FDA0003244292320000049
目标速度大小为
Figure FDA00032442923200000410
依照如下方向构建相对坐标系Fxyz：
Figure FDA0003244292320000051
目标速度与相对坐标系Fxyz的轴向夹角为θ，得到目标速度在惯性坐标系下的表示
Figure FDA0003244292320000052
Figure FDA0003244292320000053
以夹角θ为变量，
Figure FDA0003244292320000054
取极小为目标，[0°,360°]遍历搜索确定轨道转移速度增量Δv：
Figure FDA0003244292320000055
模块M3.3：轨道转移耗肼量计算，计算得到轨道转移最小速度增量对应的燃料消耗Δm：
Figure FDA0003244292320000056
10.根据权利要求6所述的一种大入轨偏差情况下组合体分离轨道确定系统，其特征在于，所述模块M4包括如下模块：
模块M4.1：确定大入轨偏差条件下的期望调整参数序列Orb′，序列包含半长轴a′、倾角i′与升交点赤经Ω′中某一项或两项参数；
模块M4.2：给定大入轨偏差条件下的分离轨道参数序列Orb，分离轨道参数序列包含期望调整参数与自由参数；
模块M4.3：根据模块M3的计算过程，以推进舱可用燃料消耗殆尽为目标，积分迭代确定分离轨道参数序列中的自由参数，确定输出分离轨道参数。
CN202111028656.0A 2021-09-02 2021-09-02 一种大入轨偏差情况下组合体分离轨道确定方法和系统 Active CN113779788B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111028656.0A CN113779788B (zh) 2021-09-02 2021-09-02 一种大入轨偏差情况下组合体分离轨道确定方法和系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111028656.0A CN113779788B (zh) 2021-09-02 2021-09-02 一种大入轨偏差情况下组合体分离轨道确定方法和系统

Publications (2)

Publication Number Publication Date
CN113779788A true CN113779788A (zh) 2021-12-10
CN113779788B CN113779788B (zh) 2023-12-12

Family

ID=78840797

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111028656.0A Active CN113779788B (zh) 2021-09-02 2021-09-02 一种大入轨偏差情况下组合体分离轨道确定方法和系统

Country Status (1)

Country Link
CN (1) CN113779788B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100250031A1 (en) * 2009-03-27 2010-09-30 The Charles Stark Draper Laboratory, Inc. Propulsive guidance for atmospheric skip entry trajectories
CN106570285A (zh) * 2016-11-09 2017-04-19 中国人民解放军装备学院 一种基于状态传递矩阵解析解的J2摄动Lambert问题求解方法
CN106989761A (zh) * 2017-05-25 2017-07-28 北京航天自动控制研究所 一种基于自适应滤波的空间飞行器制导工具在轨标定方法
CN112257343A (zh) * 2020-10-22 2021-01-22 上海卫星工程研究所 一种高精度地面轨迹重复轨道优化方法及系统

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100250031A1 (en) * 2009-03-27 2010-09-30 The Charles Stark Draper Laboratory, Inc. Propulsive guidance for atmospheric skip entry trajectories
CN106570285A (zh) * 2016-11-09 2017-04-19 中国人民解放军装备学院 一种基于状态传递矩阵解析解的J2摄动Lambert问题求解方法
CN106989761A (zh) * 2017-05-25 2017-07-28 北京航天自动控制研究所 一种基于自适应滤波的空间飞行器制导工具在轨标定方法
CN112257343A (zh) * 2020-10-22 2021-01-22 上海卫星工程研究所 一种高精度地面轨迹重复轨道优化方法及系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
李恒年 等: "有连续推力控制的卫星轨道确定算法", 系统工程与电子技术, vol. 32, no. 9, pages 1957 - 1960 *
贺波勇 等: "载人登月转移轨道偏差传播机理分析与稳健性设计", 物理学报, vol. 62, no. 19, pages 1 - 7 *

Also Published As

Publication number Publication date
CN113779788B (zh) 2023-12-12

Similar Documents

Publication Publication Date Title
CN112455723B (zh) 一种火箭推力下降故障下基于rbfnn的救援轨道决策方法
CN109080854B (zh) 航天器返回预定落点的大椭圆轨道变轨规划方法
Weiland Aerodynamic data of space vehicles
CN109911249B (zh) 低推重比飞行器的星际转移有限推力入轨迭代制导方法
CN112395689B (zh) 基于凸优化的火箭故障后在线重构方法
CN110096726B (zh) 基于月球借力的geo卫星应急转移轨道快速优化设计方法
CN112329136A (zh) 一种基于平衡飞行理论的运载火箭在线飞行程序重构方法
CN113189870A (zh) 一种火箭推力下降故障下椭圆救援轨道的轨迹重规划方法
CN115258196A (zh) 低轨卫星星座组网电推进变轨策略优化方法和系统
CN115743619A (zh) 超低轨道卫星星下点轨迹机动与保持方法、装置及介质
Heiligers et al. Design of optimal Earth pole-sitter transfers using low-thrust propulsion
Baranov et al. Ballistic aspects of large-size space debris flyby at low Earth near-circular orbits
CN113779788A (zh) 一种大入轨偏差情况下组合体分离轨道确定方法和系统
Xie et al. Guidance, navigation, and control for spacecraft rendezvous and docking: theory and methods
Petukhov Quasioptimal control with feedback for multiorbit low-thrust transfer between noncoplanar elliptic and circular orbits
CN112009727B (zh) 平动点轨道最优小推力转移分段设计方法
Wilken et al. Mission design for point-to-point passenger transport with reusable launch vehicles
CN112393648A (zh) 一种用于火箭推力故障模式下自主控制的平衡飞行理论方法
CN113741551A (zh) 一种基于代理模型的全过程轨迹优化方法及装置
Pontani et al. Lunar orbit dynamics and maneuvers for Lunisat missions
JP6542581B2 (ja) 宇宙機とその軌道面変更方法
Kumpel et al. A Conceptual Design for the Space Launch Capability of the Peacekeeper ICBM
CN113761670A (zh) 一种运载火箭的平衡飞行理论及在线入轨能力评估方法
Sell et al. Powered flight design and performance summary for the Mars science laboratory mission
CN111596678A (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