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

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

Info

Publication number
CN113779788B
CN113779788B CN202111028656.0A CN202111028656A CN113779788B CN 113779788 B CN113779788 B CN 113779788B CN 202111028656 A CN202111028656 A CN 202111028656A CN 113779788 B CN113779788 B CN 113779788B
Authority
CN
China
Prior art keywords
orbit
track
parameters
spacecraft
determining
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.)
Active
Application number
CN202111028656.0A
Other languages
English (en)
Other versions
CN113779788A (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

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地球引力摄动影响下的轨道转移运动学模型：
其中，分别为点火时刻组合体航天器地心惯性系位置、速度矢量；
r为点火时刻组合体航天器位置矢量的模值；
(x,y,z)为点火时刻航天器地心惯性系三轴位置分量；
为点火时刻组合体航天器地心惯性系三轴速度分量；
(α,δ)为组合体航天器推力方向的赤经、赤纬；
μ＝3.986005×1014m3·s-2为地球引力常数；
Re＝6378.14×103m为地球赤道半径；
g0＝9.80665m·s-2为地球重力加速度。
优选地，所述步骤S3包括如下步骤：
步骤S3.1：转移轨道点火点位置计算，目标升交点赤经Ωs、目标倾角is，计算得到目标轨道平面正法向矢量为：
在星箭分离轨道上积分寻找t0时刻，确定点火点的位置速度
步骤S3.2：转移轨道速度增量计算，点火点P0的位置速度目标速度大小为依照如下方向构建相对坐标系Fxyz：
目标速度与相对坐标系Fxyz的轴向夹角为θ，得到目标速度在惯性坐标系下的表示
Fxyz'为构建的相对坐标系，以夹角θ为变量，取极小为目标，[0°,360°]遍历搜索确定轨道转移速度增量Δv：
步骤S3.3：轨道转移耗肼量计算，计算得到轨道转移最小速度增量对应的燃料消耗Δm：
优选地，所述步骤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地球引力摄动影响下的轨道转移运动学模型：
其中，分别为点火时刻组合体航天器地心惯性系位置、速度矢量；
r为点火时刻组合体航天器位置矢量的模值；
(x,y,z)为点火时刻航天器地心惯性系三轴位置分量；
为点火时刻组合体航天器地心惯性系三轴速度分量；
(α,δ)为组合体航天器推力方向的赤经、赤纬；
μ＝3.986005×1014m3·s-2为地球引力常数；
Re＝6378.14×103m为地球赤道半径；
g0＝9.80665m·s-2为地球重力加速度。
优选地，所述模块M3包括如下模块：
模块M3.1：转移轨道点火点位置计算，目标升交点赤经Ωs、目标倾角is，计算得到目标轨道平面正法向矢量为：
在星箭分离轨道上积分寻找t0时刻，确定点火点的位置速度
模块M3.2：转移轨道速度增量计算，点火点P0的位置速度目标速度大小为依照如下方向构建相对坐标系Fxyz：
目标速度与相对坐标系Fxyz的轴向夹角为θ，得到目标速度在惯性坐标系下的表示
以夹角θ为变量，取极小为目标，[0°,360°]遍历搜索确定轨道转移速度增量Δv：
模块M3.3：轨道转移耗肼量计算，计算得到轨道转移最小速度增量对应的燃料消耗Δm：
优选地，所述模块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地球引力摄动影响下的轨道转移运动学模型：
其中，分别为点火时刻组合体航天器地心惯性系位置、速度矢量；
r为点火时刻组合体航天器位置矢量的模值；
(x,y,z)为点火时刻航天器地心惯性系三轴位置分量；
为点火时刻组合体航天器地心惯性系三轴速度分量；
(α,δ)为组合体航天器推力方向的赤经、赤纬；
μ＝3.986005×1014m3·s-2为地球引力常数；
Re＝6378.14×103m为地球赤道半径；
g0＝9.80665m·s-2为地球重力加速度。
步骤S3：利用交轨共点方法对虚拟分离轨道点火点位置、速度增量及燃料消耗进行估计。
步骤S3.1：转移轨道点火点位置计算，即：
假设目标升交点赤经Ωs、目标倾角is，则可以计算得到目标轨道平面正法向矢量为：
通过在星箭分离轨道上积分寻找合适的t0时刻，使得如下条件满足，可以相应确定点火点P0的位置速度
步骤S3.2：转移轨道速度增量计算，即：
已知点火点P0的位置速度假设由活力公式得到的目标速度大小为/>则依照如下方向构建相对坐标系Fxyz：
假设目标速度与相对坐标系Fxyz的x轴向夹角为θ，则可以得到目标速度在惯性坐标系下的表示
以夹角θ为变量，取极小为目标，[0°,360°]遍历搜索即可确定轨道转移速度增量Δv：
步骤S3.3：轨道转移耗肼量计算，即：
根据火箭公式，计算得到轨道转移最小速度增量对应的燃料消耗Δm：
步骤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地球引力摄动影响下的轨道转移运动学模型：
其中，分别为点火时刻组合体航天器地心惯性系位置、速度矢量；
r为点火时刻组合体航天器位置矢量的模值；
(x,y,z)为点火时刻航天器地心惯性系三轴位置分量；
为点火时刻组合体航天器地心惯性系三轴速度分量；
(α,δ)为组合体航天器推力方向的赤经、赤纬；
μ＝3.986005×1014m3·s-2为地球引力常数；
Re＝6378.14×103m为地球赤道半径；
g0＝9.80665m·s-2为地球重力加速度。
模块M3：利用交轨共点方法对分离轨道点火点位置、速度增量及燃料消耗进行计算；模块M3.1：转移轨道点火点位置计算，目标升交点赤经Ωs、目标倾角is，计算得到目标轨道平面正法向矢量为：
在星箭分离轨道上积分寻找t0时刻，确定点火点的位置速度
模块M3.2：转移轨道速度增量计算，点火点P0的位置速度目标速度大小为依照如下方向构建相对坐标系Fxyz：
目标速度与相对坐标系Fxyz的轴向夹角为θ，得到目标速度在惯性坐标系下的表示
以夹角θ为变量，取极小为目标，[0°,360°]遍历搜索确定轨道转移速度增量Δv：
模块M3.3：轨道转移耗肼量计算，计算得到轨道转移最小速度增量对应的燃料消耗Δm：
模块M4：以期望参数调整至目标值为原则，组合体推进舱最大可用燃料量为边界，积分迭代实现对分离轨道参数的确定。模块M4.1：确定大入轨偏差条件下的期望调整参数序列Orb′，序列包含半长轴a′、倾角i′与升交点赤经Ω′中某一项或两项参数；模块M4.2：给定大入轨偏差条件下的分离轨道参数序列Orb，分离轨道参数序列包含期望调整参数与自由参数；模块M4.3：根据模块M3的计算过程，以推进舱可用燃料消耗殆尽为目标，积分迭代确定分离轨道参数序列中的自由参数，确定输出分离轨道参数。
本发明解决了发生较大入轨偏差，纠偏燃料消耗导致原有转移任务无法完成情况下，组合体航天器最优分离轨道确定方法，为异常状况下航天任务补救提供支撑，具有一定的工程实用性。
本领域技术人员知道，除了以纯计算机可读程序代码方式实现本发明提供的系统及其各个装置、模块、单元以外，完全可以通过将方法步骤进行逻辑编程来使得本发明提供的系统及其各个装置、模块、单元以逻辑门、开关、专用集成电路、可编程逻辑控制器以及嵌入式微控制器等的形式来实现相同功能。所以，本发明提供的系统及其各项装置、模块、单元可以被认为是一种硬件部件，而对其内包括的用于实现各种功能的装置、模块、单元也可以视为硬件部件内的结构；也可以将用于实现各种功能的装置、模块、单元视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
以上对本发明的具体实施例进行了描述。需要理解的是，本发明并不局限于上述特定实施方式，本领域技术人员可以在权利要求的范围内做出各种变化或修改，这并不影响本发明的实质内容。在不冲突的情况下，本申请的实施例和实施例中的特征可以任意相互组合。

Claims (8)

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

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
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 上海卫星工程研究所 一种高精度地面轨迹重复轨道优化方法及系统

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8489258B2 (en) * 2009-03-27 2013-07-16 The Charles Stark Draper Laboratory, Inc. Propulsive guidance for atmospheric skip entry trajectories

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
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
有连续推力控制的卫星轨道确定算法;李恒年 等;系统工程与电子技术;第32卷(第9期);第1957-1960页 *
载人登月转移轨道偏差传播机理分析与稳健性设计;贺波勇 等;物理学报;第62卷(第19期);第1-7页 *

Also Published As

Publication number Publication date
CN113779788A (zh) 2021-12-10

Similar Documents

Publication Publication Date Title
CN112455723B (zh) 一种火箭推力下降故障下基于rbfnn的救援轨道决策方法
CN112257343B (zh) 一种高精度地面轨迹重复轨道优化方法及系统
CN109080854B (zh) 航天器返回预定落点的大椭圆轨道变轨规划方法
CN109911249B (zh) 低推重比飞行器的星际转移有限推力入轨迭代制导方法
CN110989644A (zh) 一种考虑目标点多终端约束的飞行器轨迹规划方法
Pasand et al. A study of spacecraft reaction thruster configurations for attitude control system
CN110244767A (zh) 一种采用有限元方法的编队构型重构优化
Zhu et al. Rotating object specific tracking based on orbit-attitude coordinated adaptive control
CN115258196A (zh) 低轨卫星星座组网电推进变轨策略优化方法和系统
CN113779788B (zh) 一种大入轨偏差情况下组合体分离轨道确定方法和系统
Fan et al. Quick-response attitude takeover control using multiple servicing spacecraft based on inertia properties identification
CN112231831A (zh) 固体运载火箭强耦合条件下终端多约束能量管理方法
Tao et al. Attitude maneuvering and vibration reducing control of flexible spacecraft subject to actuator saturation and misalignment
Misra et al. Asteroid landing guidance design in the framework of coupled orbit-attitude spacecraft dynamics
Lian et al. RBF network based adaptive sliding mode control for solar sails
Petukhov Quasioptimal control with feedback for multiorbit low-thrust transfer between noncoplanar elliptic and circular orbits
Losa et al. Electric station keeping of geostationary satellites: a differential inclusion approach
CN111596678B (zh) 一种航天器逼近空间翻滚目标的相对位置自适应控制方法
Pontani et al. Lunar orbit dynamics and maneuvers for Lunisat missions
Moradi et al. Attitude tracking control using an online identification and a linear quadratic regulator–based strategy in the presence of orbital eccentricity
Omagari et al. CMG configuration and control for rapid attitude maneuver of small spacecraft
Sedelnikov et al. Dynamic characteristics modeling of rotary platform installed on board of a small spacecraft
CN114872934B (zh) 基于无损凸化和Chebyshev正交配点法的月面飞行器轨迹优化方法
CN113247310B (zh) 一种适用于卫星可连续姿态机动次数的估算方法及系统
CN116039970B (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