CN108100307A - 一种用于复杂约束下低能量小天体精确探测轨道转移方法 - Google Patents

一种用于复杂约束下低能量小天体精确探测轨道转移方法 Download PDF

Info

Publication number
CN108100307A
CN108100307A CN201711264944.XA CN201711264944A CN108100307A CN 108100307 A CN108100307 A CN 108100307A CN 201711264944 A CN201711264944 A CN 201711264944A CN 108100307 A CN108100307 A CN 108100307A
Authority
CN
China
Prior art keywords
orbit
constraints
transfer
quasi
detector
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
CN201711264944.XA
Other languages
English (en)
Other versions
CN108100307B (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.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
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 Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN201711264944.XA priority Critical patent/CN108100307B/zh
Publication of CN108100307A publication Critical patent/CN108100307A/zh
Application granted granted Critical
Publication of CN108100307B publication Critical patent/CN108100307B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00Cosmonautic vehicles
    • B64G1/10Artificial satellites; Systems of such satellites; Interplanetary vehicles
    • B64G1/105Space science
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00Cosmonautic vehicles
    • B64G1/10Artificial satellites; Systems of such satellites; Interplanetary vehicles
    • B64G1/105Space science
    • B64G1/1064Space science specifically adapted for interplanetary, solar or interstellar exploration
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00Cosmonautic vehicles
    • B64G1/22Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
    • B64G1/24Guiding or controlling apparatus, e.g. for attitude control
    • B64G1/242Orbits and trajectories

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Astronomy & Astrophysics (AREA)
  • General Physics & Mathematics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Sustainable Development (AREA)
  • Measuring Magnetic Variables (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开的一种用于复杂约束下低能量小天体精确探测轨道转移方法,属于航空航天技术领域。本发明首先确定探测轨道设计任务所需满足的多种复杂非一致强耦合约束,建立多种复杂非一致强耦合约束与轨道设计参数的映射关系;在质心旋转坐标系下建立探测器动力学方程;通过建立的线性化探测器动力学方程提供初值,采用非线性降维方法和二阶微分修正得到星历模型下精确的拟周期轨道;基于星历模型下精确的拟周期轨道,采用拟流形扰动法优化获得转移轨道初值;针对多种复杂非一致强耦合约束对得到的转移轨道初值进行修正,获得精确的低能量转移轨道。本发明具有效率高、收敛性好、转移所需能量小的优点。

Description

一种用于复杂约束下低能量小天体精确探测轨道转移方法
技术领域
本发明涉及一种用于复杂约束下低能量小天体精确探测轨道转移方法,尤其涉及适用于考虑燃料约束,测控约束和飞行机动状态约束的小天体低能量转移轨道实现方法,属于航空航天技术领域。
背景技术
小天体富含的矿物质可成为人类开采的重要来源,其蕴含的有机物也将为生命起源提供直接信息与见解。因此小天体探测已经在多个航天大国的深空探测任务中取得了不同程度的进展,已成为未来世界航天发展的新方向。
在已发展的关于探测器小行星探测轨道设计方法中在先技术[1](参见:乔栋,崔祜涛,崔平远.小行星探测最优两脉冲交会轨道设计与分析[J].宇航学报2005,26(3):362-367)给出了小天体探测中任意两个非共面非共轴椭圆轨道之间的最优两脉冲转移方法。在先技术[2](参见:李九天,罗亚中,唐国金.小行星探测多脉冲交会轨道多目标优化[J].国防科技大学学报,2011,33(3):5-9.)基于Lambert交会算法建立了包含地球逃逸轨道和日心转移轨道的多脉冲交会轨道,并以燃料消耗最小和转移时间最短为优化变量,从而对小行星探测轨道进行优化设计。在先技术[3](参加:陈杨,宝音贺西,李俊峰.我国小行星探测目标分析与电推进轨道设计[J].中国科学物理学力学天文学,2011,(09):1104-1111)采用间接法求解燃料最优控制问题最终得到小天体探测小推力的最优轨迹。对于上述小天体探测任务轨道设计方法均从地球出发,不适用与平衡点出发的探测轨道设计,且仅能对少量约束进行考虑,不能确保满足工程实践任务中的不同的强耦合约束。
发明内容
本发明公开的一种用于复杂约束下低能量小天体精确探测轨道转移方法要解决的技术问题是,以满足多种复杂非一致强耦合约束为前提,实现多约束低能量小天体精确探测轨道转移,具有效率高、收敛性好、转移所需能量小的优点。所述的多种复杂非一致强耦合约束包括考虑燃料约束,测控约束和飞行机动状态约束。
本发明的目的是通过下述技术方案实现的。
本发明公开的一种用于复杂约束下低能量小天体精确探测轨道转移方法,首先确定探测轨道设计任务所需满足的多种复杂非一致强耦合约束,建立多种复杂非一致强耦合约束与轨道设计参数的映射关系。在质心旋转坐标系下建立探测器动力学方程。通过建立的线性化探测器动力学方程提供初值,采用非线性降维方法和二阶微分修正得到星历模型下精确的拟周期轨道。基于星历模型下精确的拟周期轨道,采用拟流形扰动法优化获得转移轨道初值。针对多种复杂非一致强耦合约束对得到的转移轨道初值进行修正,获得精确的低能量转移轨道。
本发明公开的一种用于复杂约束下低能量小天体精确探测轨道转移方法,包括如下步骤:
步骤一:确定探测轨道设计任务所需满足的多种复杂非一致强耦合约束,建立多种复杂非一致强耦合约束与轨道设计参数的映射关系。
确定探测轨道设计任务所需满足的多种复杂非一致强耦合约束,所述的多种复杂非一致强耦合约束包括考虑燃料约束、测控约束和飞行机动状态约束。首先考虑燃料约束,即所提供的速度增量有限,通过燃料质量确定最大速度增量。其次需要满足深空探测测控约束,即要求飞行关键节点探测器距离小于通信有效距离,通信无遮挡、地面可测控。最后需要满足飞行机动状态约束,即转移时间小于任务时间。所述的飞行关键节点包括变轨、修正、飞越。
步骤二:在质心旋转坐标系下建立探测器动力学方程。
定义太阳为质量为M1主天体P1,地球为质量为M2的主天体P2,探测器P3质量为M3。将P1、P2与P3视为质点,且P1、P2与P3的质量满足关系M1>M2>>M3。约束太阳与地球的运动为圆运动,则建立起圆型限制性三体问题。选择太阳与地球构成的系统的质心为原点建立质心旋转坐标系。X-Y坐标平面即两个主天体相对运动平面,X轴方向由太阳P1质心指向地球P2质心,Z轴方向与系统角速度方向相同,Y轴与X轴、Z轴构成右手坐标系。
选取特征长度为主天体质心间距离,特征质量为主天体质量之和M1+M2。无量纲化后的探测器P3在质心旋转坐标系的动力学方程为,
其中:μ表示太阳的引力常数,xyz表示探测器P3在质心旋转坐标系Oxyz下的位置矢量,表示探测器P3在质心旋转坐标系Oxyz下的速度矢量,表示探测器P3在质心旋转坐标系Oxyz下的加速度矢量。
步骤三:通过步骤二建立的线性化探测器动力学方程提供初值,通过非线性降维方法和二阶微分修正得到星历模型下精确的拟周期轨道。
对步骤二建立的线性化探测器动力学方程求解可知,探测器在平衡点附近运动方程的通解为,
其中:ξ、η与ζ为扰动坐标,φ、描述每个具体的初始点,Ai(i=1,2,x,z)为由初始条件确定的积分常数,A1、A2为双曲振幅,Ax、Az表示拟周期轨道平面内与平面内的振幅。c、k是与轨道的雅各比常数C相关的常数,分别表示为,
三个特征频率ω、ν与λ分别为,
由于圆型限制性三体模型忽略系统偏心率以及摄动因素,使得该模型下的拟周期轨道设计结果不满足工程约束,且存在不收敛的情况,因而采用非线性降维方法,在相空间设置若干截面,并通过方程(2)给出的线性化运动提供初值给拟周期轨道,确定轨道穿过截面的交点作为节点,采用二级微分修正方法改变交点在截面上的位置,实现满足始末端约束且节点位置与速度连续的轨道。
对探测器在拟周期轨道上的位置进行修正时,首先需根据近似解析解提供所需修正点的状态矢量与时间的初值X0,X1,X2与t0,t1,t2;然后令控制变量为初始点的速度分量约束变量为由初始点轨道递推后的终端状态的位置矢量终止条件为固定的积分时间。
根据式(6)求解微分修正关系式,反复迭代至误差位置在容许范围内。重复此步骤至所有轨迹都位置连续。
对速度进行修正时,通过改变修正点的位置矢量使连接点处的速度跳变Δv接近于零。对于n段轨迹,有n+1个连接点x0,x1,……xn。相应的微分修正关系式为,
即实现通过非线性降维方法和二阶微分修正得到星历模型下精确的拟周期轨道。
步骤四:基于步骤三得到的星历模型下精确的拟周期轨道,采用拟流形扰动法优化获得转移轨道初值。
单值矩阵M为状态转移矩阵Φ在一个轨道周期T处的取值,即:
M=Φ(t=T) (9)
单值矩阵的特征值对应的特征向量可决定拟周期轨道附近的运动特征,其稳定与不稳定特征向量可决定受到初始扰动时的发散或收敛。拟流形扰动法利用单值矩阵M不稳定特征向量作为扰动量施加的方向,使得探测器远离初始轨道并到达目标小行星附近。施加的扰动速度矢量ΔV表达式为,
ΔV=λu (10)
其中:λ速度扰动量大小,u为不稳定特征向量速度分量的单位向量。针对不同的探测器从拟周期轨道分离的时刻t0,以速度扰动量大小λ为优化变量得到转移轨道与目标小天体相距最近的距离D与相距最近的时刻tf,从而为转移轨道设计初值提供依据。
即探测器从拟周期轨道分离的时刻t0、速度扰动量大小λ为采用拟流形扰动法优化获得转移轨道初值。
步骤五:针对步骤一确定的多种复杂非一致强耦合约束对步骤四得到的转移轨道初值进行修正,获得精确的低能量转移轨道。
综合考虑多种复杂非一致强耦合约束,根据有限的速度增量大小对速度扰动量大小λ进行约束;考虑深空测控约束,以步骤四得到的转移轨道初值为基础,设置满足飞行机动状态约束的节点范围,利用二级微分修正,设计不同机动时刻的小行星精确转移轨道,最后考虑飞行机动状态约束,筛选出符合给定时间范围内且燃料消耗最省的转移轨道,即以满足多种复杂非一致强耦合约束为前提,实现低能量小天体精确探测轨道转移。
有益效果:
1、本发明公开的一种用于复杂约束下低能量小天体精确探测轨道转移方法,采用非线性降维方法和微分修正方法,能够快速的找到非线性强扰动三体模型下精确稳定的拟周期轨道,具有效率高、收敛性好的优点。
2、本发明公开的一种用于复杂约束下低能量小天体精确探测轨道转移方法,采用拟流形扰动设计转移轨道,利用动力学特性进行小天体转移,转移所需能量小。
3、本发明公开的一种用于复杂约束下低能量小天体精确探测轨道转移方法,通过确定探测轨道设计任务所需满足多种复杂非一致强耦合约束并建立多种复杂非一致强耦合约束与轨道设计参数的映射关系,实现对转移轨道的筛选,即以满足多种复杂非一致强耦合约束为前提,实现低能量小天体精确探测轨道转移,具有较好的工程实用性。
附图说明:
图1是本发明一种用于复杂约束下低能量小天体精确探测轨道转移方法的流程图;
图2是本发明一种用于复杂约束下低能量小天体精确探测轨道转移方法的质心旋转坐标系的示意图;
图3是本发明一种用于复杂约束下低能量小天体精确探测轨道转移方法的星历模型下的日地L2点拟周期轨道;
图4是本发明一种用于复杂约束下低能量小天体精确探测轨道转移方法的飞越目标小行星的转移轨道设计方案。
具体实施方式
为了更好地说明本发明的目的和优点,下面选择2001AV43小行星作为目标星,探测器初始轨道为日地L2点拟周期轨道,进行低能量小天体转移轨道设计,对本发明做出详细解释。
实施例1:
如图1所示,本实例公开的一种用于复杂约束下低能量小天体精确探测轨道转移方法,包括如下步骤:
步骤一:确定探测轨道设计任务所需满足的多种复杂非一致强耦合约束,建立多种复杂非一致强耦合约束与轨道设计参数的映射关系。
确定探测轨道设计任务所需满足的多种复杂非一致强耦合约束,所述的多种复杂非一致强耦合约束包括考虑燃料约束、测控约束和飞行机动状态约束。首先考虑燃料约束,即所提供的速度增量有限,通过燃料质量确定最大速度增量,在实施例中给定变轨能力约束为200m/s,即速度扰动量-200m/s≤λ≤200m/s。其次需要满足深空探测测控约束,即要求飞行关键节点探测器距离小于通信有效距离,通信无遮挡、地面可测控。最后需要满足飞行机动状态约束,即转移时间小于任务时间。所述的飞行关键节点包括变轨、修正、飞越。
步骤二:在质心旋转坐标系下建立探测器动力学方程。
定义太阳为质量为M1主天体P1,地球为质量为M2的主天体P2,探测器P3质量为M3。将P1、P2与P3视为质点,且P1、P2与P3的质量满足关系M1>M2>>M3。约束太阳与地球的运动为圆运动,则建立起圆型限制性三体问题。选择太阳与地球构成的系统的质心为原点建立质心旋转坐标系。X-Y坐标平面即两个主天体相对运动平面,X轴方向由太阳P1质心指向地球P2质心,Z轴方向与系统角速度方向相同,Y轴与X轴、Z轴构成右手坐标系。
选取特征长度为主天体质心间距离,特征质量为主天体质量之和M1+M2。无量纲化后的探测器P3在质心旋转坐标系的动力学方程为,
其中:μ表示太阳的引力常数,xyz表示探测器P3在质心旋转坐标系Oxyz下的位置矢量,表示探测器P3在质心旋转坐标系Oxyz下的速度矢量,表示探测器P3在质心旋转坐标系Oxyz下的加速度矢量。
步骤三:通过步骤二建立的线性化探测器动力学方程提供初值,通过非线性降维方法和二阶微分修正得到星历模型下精确的拟周期轨道。
对步骤二建立的线性化探测器动力学方程求解可知,探测器在平衡点附近运动方程的通解为,
其中:ξ、η与ζ为扰动坐标,φ、描述每个具体的初始点,Ai(i=1,2,x,z)为由初始条件确定的积分常数,A1、A2为双曲振幅,Ax、Az表示拟周期轨道平面内与平面内的振幅。c、k是与轨道的雅各比常数C相关的常数,分别表示为,
三个特征频率ω、ν与λ分别为,
由于圆型限制性三体模型忽略系统偏心率以及摄动因素,使得该模型下的拟周期轨道设计结果不满足工程约束,且存在不收敛的情况,因而采用非线性降维方法,在相空间设置若干截面,并通过方程(2)给出的线性化运动提供初值给拟周期轨道,确定轨道穿过截面的交点作为节点,采用二级微分修正方法改变交点在截面上的位置,实现满足始末端约束且节点位置与速度连续的轨道。
对探测器在拟周期轨道上的位置进行修正时,首先需根据近似解析解提供所需修正点的状态矢量与时间的初值X0,X1,X2与t0,t1,t2;然后令控制变量为初始点的速度分量约束变量为由初始点轨道递推后的终端状态的位置矢量终止条件为固定的积分时间。
根据式(6)求解微分修正关系式,反复迭代至误差位置在容许范围内。重复此步骤至所有轨迹都位置连续。
对速度进行修正时,通过改变修正点的位置矢量使连接点处的速度跳变Δv接近于零。对于n段轨迹,有n+1个连接点x0,x1,……xn。相应的微分修正关系式为,
即实现通过非线性降维方法和二阶微分修正得到星历模型下精确的拟周期轨道。轨道的初始时刻为2033年3月1日,Z向振幅为345750km,修正后星历模型下拟周期轨道如图3所示。
步骤四:基于步骤三得到的星历模型下精确的拟周期轨道,采用拟流形扰动法优化获得转移轨道初值。
单值矩阵M为状态转移矩阵Φ在一个轨道周期T处的取值,即:
M=Φ(t=T) (9)
单值矩阵的特征值对应的特征向量可决定拟周期轨道附近的运动特征,其稳定与不稳定特征向量可决定受到初始扰动时的发散或收敛。拟流形扰动法利用单值矩阵M不稳定特征向量作为扰动量施加的方向,使得探测器远离初始轨道并到达目标小行星2001AV43附近。施加的扰动速度矢量ΔV表达式为,
ΔV=λu (10)
其中:λ速度扰动量大小,u为不稳定特征向量速度分量的单位向量。针对不同的探测器从拟周期轨道分离的时刻t0,以速度扰动量大小λ为优化变量得到转移轨道与目标小天体2001AV43相距最近的距离D与相距最近的时刻tf,从而为转移轨道设计初值提供依据。
即探测器从拟周期轨道分离的时刻t0、速度扰动量大小λ为采用拟流形扰动法优化获得转移轨道初值。
步骤五:针对步骤一确定的约束对选定的初值对步骤五得到的轨道初值进行修正,获得精确的低能量转移轨道。
综合考虑非一致强耦合约束,根据有限的速度增量大小对速度扰动量大小λ进行约束;考虑深空测控约束,以设计初值为基础,设置满足约束的节点范围,利用二级微分修正,设计不同机动时刻的小行星精确转移轨道,最后考虑时间约束,筛选出符合给定时间范围内且燃料消耗最省的转移轨道,即为满足多约束的低能量小行星探测精确轨道设计。图4给出了满足多约束的小行星Toutatis的转移轨道设计方案。其中转移时间为133.9天,所需速度增量为190m/s,小行星飞越时间为2038年8月32日。
以上所述的具体描述,对发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例,用于解释本发明,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (6)

1.一种用于复杂约束下低能量小天体精确探测轨道转移方法,其特征在于:包括如下步骤,
步骤一:确定探测轨道设计任务所需满足的多种复杂非一致强耦合约束,建立多种复杂非一致强耦合约束与轨道设计参数的映射关系;
步骤二:在质心旋转坐标系下建立探测器动力学方程;
步骤三:通过步骤二建立的线性化探测器动力学方程提供初值,通过非线性降维方法和二阶微分修正得到星历模型下精确的拟周期轨道;
步骤四:基于步骤三得到的星历模型下精确的拟周期轨道,采用拟流形扰动法优化获得转移轨道初值;
步骤五:针对步骤一确定的多种复杂非一致强耦合约束对步骤四得到的转移轨道初值进行修正,获得精确的低能量转移轨道。
2.如权利要求1所述的一种用于复杂约束下低能量小天体精确探测轨道转移方法,其特征在于:步骤一的具体实现方法为,
确定探测轨道设计任务所需满足的多种复杂非一致强耦合约束,所述的多种复杂非一致强耦合约束包括考虑燃料约束、测控约束和飞行机动状态约束;首先考虑燃料约束,即所提供的速度增量有限,通过燃料质量确定最大速度增量;其次需要满足深空探测测控约束,即要求飞行关键节点探测器距离小于通信有效距离,通信无遮挡、地面可测控;最后需要满足飞行机动状态约束,即转移时间小于任务时间;所述的飞行关键节点包括变轨、修正、飞越。
3.如权利要求2所述的一种用于复杂约束下低能量小天体精确探测轨道转移方法,其特征在于:步骤二的具体实现方法为,
定义太阳为质量为M1主天体P1,地球为质量为M2的主天体P2,探测器P3质量为M3;将P1、P2与P3视为质点,且P1、P2与P3的质量满足关系M1>M2>>M3;约束太阳与地球的运动为圆运动,则建立起圆型限制性三体问题;选择太阳与地球构成的系统的质心为原点建立质心旋转坐标系;X-Y坐标平面即两个主天体相对运动平面,X轴方向由太阳P1质心指向地球P2质心,Z轴方向与系统角速度方向相同,Y轴与X轴、Z轴构成右手坐标系;
选取特征长度为主天体质心间距离,特征质量为主天体质量之和M1+M2;无量纲化后的探测器P3在质心旋转坐标系的动力学方程为,
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mover> <mi>x</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mo>-</mo> <mn>2</mn> <mover> <mi>y</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>-</mo> <mi>x</mi> <mo>=</mo> <mfrac> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>&amp;mu;</mi> <mo>)</mo> <mo>(</mo> <mi>x</mi> <mo>+</mo> <mi>&amp;mu;</mi> <mo>)</mo> </mrow> <msup> <mrow> <mo>&amp;lsqb;</mo> <msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>+</mo> <mi>&amp;mu;</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>z</mi> <mn>2</mn> </msup> <mo>&amp;rsqb;</mo> </mrow> <mrow> <mn>3</mn> <mo>/</mo> <mn>2</mn> </mrow> </msup> </mfrac> <mo>-</mo> <mfrac> <mrow> <mi>&amp;mu;</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>-</mo> <mn>1</mn> <mo>+</mo> <mi>&amp;mu;</mi> <mo>)</mo> </mrow> </mrow> <msup> <mrow> <mo>&amp;lsqb;</mo> <msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>-</mo> <mn>1</mn> <mo>+</mo> <mi>&amp;mu;</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>z</mi> <mn>2</mn> </msup> <mo>&amp;rsqb;</mo> </mrow> <mrow> <mn>3</mn> <mo>/</mo> <mn>2</mn> </mrow> </msup> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mover> <mi>y</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mo>+</mo> <mn>2</mn> <mover> <mi>x</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>-</mo> <mi>y</mi> <mo>=</mo> <mfrac> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>&amp;mu;</mi> <mo>)</mo> <mi>y</mi> </mrow> <msup> <mrow> <mo>&amp;lsqb;</mo> <msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>+</mo> <mi>&amp;mu;</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>z</mi> <mn>2</mn> </msup> <mo>&amp;rsqb;</mo> </mrow> <mrow> <mn>3</mn> <mo>/</mo> <mn>2</mn> </mrow> </msup> </mfrac> <mo>-</mo> <mfrac> <mrow> <mi>&amp;mu;</mi> <mi>y</mi> </mrow> <msup> <mrow> <mo>&amp;lsqb;</mo> <msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>-</mo> <mn>1</mn> <mo>+</mo> <mi>&amp;mu;</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>z</mi> <mn>2</mn> </msup> <mo>&amp;rsqb;</mo> </mrow> <mrow> <mn>3</mn> <mo>/</mo> <mn>2</mn> </mrow> </msup> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mover> <mi>z</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mo>=</mo> <mo>-</mo> <mfrac> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>&amp;mu;</mi> <mo>)</mo> <mi>z</mi> </mrow> <msup> <mrow> <mo>&amp;lsqb;</mo> <msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>+</mo> <mi>&amp;mu;</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>z</mi> <mn>2</mn> </msup> <mo>&amp;rsqb;</mo> </mrow> <mrow> <mn>3</mn> <mo>/</mo> <mn>2</mn> </mrow> </msup> </mfrac> <mo>-</mo> <mfrac> <mrow> <mi>&amp;mu;</mi> <mi>z</mi> </mrow> <msup> <mrow> <mo>&amp;lsqb;</mo> <msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>-</mo> <mn>1</mn> <mo>+</mo> <mi>&amp;mu;</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> <mo>+</mo> <msup> <mi>z</mi> <mn>2</mn> </msup> <mo>&amp;rsqb;</mo> </mrow> <mrow> <mn>3</mn> <mo>/</mo> <mn>2</mn> </mrow> </msup> </mfrac> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
其中:μ表示太阳的引力常数,xyz表示探测器P3在质心旋转坐标系Oxyz下的位置矢量,表示探测器P3在质心旋转坐标系Oxyz下的速度矢量,表示探测器P3在质心旋转坐标系Oxyz下的加速度矢量。
4.如权利要求3所述的一种用于复杂约束下低能量小天体精确探测轨道转移方法,其特征在于:步骤三的具体实现方法为,
对步骤二建立的线性化探测器动力学方程求解可知,探测器在平衡点附近运动方程的通解为,
其中:ξ、η与ζ为扰动坐标,φ、描述每个具体的初始点,Ai(i=1,2,x,z)为由初始条件确定的积分常数,A1、A2为双曲振幅,Ax、Az表示拟周期轨道平面内与平面内的振幅;c、k是与轨道的雅各比常数C相关的常数,分别表示为,
<mrow> <mi>c</mi> <mo>=</mo> <mfrac> <mrow> <msup> <mi>&amp;lambda;</mi> <mn>2</mn> </msup> <mo>-</mo> <mn>1</mn> <mo>-</mo> <mn>2</mn> <mi>C</mi> </mrow> <mrow> <mn>2</mn> <mi>&amp;lambda;</mi> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mi>k</mi> <mo>=</mo> <mfrac> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <msup> <mi>&amp;omega;</mi> <mn>2</mn> </msup> <mo>+</mo> <mn>1</mn> <mo>+</mo> <mn>2</mn> <mi>C</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mn>2</mn> <mi>&amp;omega;</mi> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
三个特征频率ω、ν与λ分别为,
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mi>&amp;omega;</mi> <mo>=</mo> <msqrt> <mfrac> <mrow> <mn>2</mn> <mo>-</mo> <mi>C</mi> <mo>+</mo> <msqrt> <mrow> <mn>9</mn> <msup> <mi>C</mi> <mn>2</mn> </msup> <mo>-</mo> <mn>8</mn> <mi>C</mi> </mrow> </msqrt> </mrow> <mn>2</mn> </mfrac> </msqrt> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>d</mi> <mo>=</mo> <msqrt> <mfrac> <mrow> <mi>C</mi> <mo>-</mo> <mn>2</mn> <mo>+</mo> <msqrt> <mrow> <mn>9</mn> <msup> <mi>C</mi> <mn>2</mn> </msup> <mo>-</mo> <mn>8</mn> <mi>C</mi> </mrow> </msqrt> </mrow> <mn>2</mn> </mfrac> </msqrt> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>v</mi> <mo>=</mo> <msqrt> <mi>C</mi> </msqrt> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>
由于圆型限制性三体模型忽略系统偏心率以及摄动因素,使得该模型下的拟周期轨道设计结果不满足工程约束,且存在不收敛的情况,因而采用非线性降维方法,在相空间设置若干截面,并通过方程(2)给出的线性化运动提供初值给拟周期轨道,确定轨道穿过截面的交点作为节点,采用二级微分修正方法改变交点在截面上的位置,实现满足始末端约束且节点位置与速度连续的轨道;
对探测器在拟周期轨道上的位置进行修正时,首先需根据近似解析解提供所需修正点的状态矢量与时间的初值X0,X1,X2与t0,t1,t2;然后令控制变量为初始点的速度分量约束变量为由初始点轨道递推后的终端状态的位置矢量[δx1 *,δy1 *,δz1 *],终止条件为固定的积分时间;
<mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msup> <msub> <mi>&amp;delta;x</mi> <mn>1</mn> </msub> <mo>*</mo> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msup> <msub> <mi>&amp;delta;y</mi> <mn>1</mn> </msub> <mo>*</mo> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msup> <msub> <mi>&amp;delta;z</mi> <mn>1</mn> </msub> <mo>*</mo> </msup> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <msup> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>*</mo> </msup> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mover> <mi>x</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>0</mn> </msub> </mrow> </mfrac> </mtd> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <msup> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>*</mo> </msup> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mover> <mi>y</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>0</mn> </msub> </mrow> </mfrac> </mtd> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <msup> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>*</mo> </msup> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mover> <mi>z</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>0</mn> </msub> </mrow> </mfrac> </mtd> <mtd> <mrow> <msup> <msub> <mover> <mi>x</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>1</mn> </msub> <mo>*</mo> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <msup> <msub> <mi>y</mi> <mn>1</mn> </msub> <mo>*</mo> </msup> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mover> <mi>x</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>0</mn> </msub> </mrow> </mfrac> </mtd> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <msup> <msub> <mi>y</mi> <mn>1</mn> </msub> <mo>*</mo> </msup> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mover> <mi>y</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>0</mn> </msub> </mrow> </mfrac> </mtd> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <msup> <msub> <mi>y</mi> <mn>1</mn> </msub> <mo>*</mo> </msup> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mover> <mi>z</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>0</mn> </msub> </mrow> </mfrac> </mtd> <mtd> <mrow> <msup> <msub> <mover> <mi>y</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>1</mn> </msub> <mo>*</mo> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <msup> <msub> <mi>z</mi> <mn>1</mn> </msub> <mo>*</mo> </msup> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mover> <mi>x</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>0</mn> </msub> </mrow> </mfrac> </mtd> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <msup> <msub> <mi>z</mi> <mn>1</mn> </msub> <mo>*</mo> </msup> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mover> <mi>y</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>0</mn> </msub> </mrow> </mfrac> </mtd> <mtd> <mfrac> <mrow> <mo>&amp;part;</mo> <msup> <msub> <mi>z</mi> <mn>1</mn> </msub> <mo>*</mo> </msup> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mover> <mi>z</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>0</mn> </msub> </mrow> </mfrac> </mtd> <mtd> <mrow> <msup> <msub> <mover> <mi>z</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>1</mn> </msub> <mo>*</mo> </msup> </mrow> </mtd> </mtr> </mtable> </mfenced> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mi>&amp;delta;</mi> <msub> <mover> <mi>x</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>0</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>&amp;delta;</mi> <msub> <mover> <mi>y</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>0</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>&amp;delta;</mi> <msub> <mover> <mi>z</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>0</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>&amp;delta;</mi> <mrow> <mo>(</mo> <mrow> <msub> <mi>t</mi> <mi>f</mi> </msub> <mo>-</mo> <msub> <mi>t</mi> <mn>0</mn> </msub> </mrow> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>
根据式(6)求解微分修正关系式,反复迭代至误差位置在容许范围内;重复此步骤至所有轨迹都位置连续;
对速度进行修正时,通过改变修正点的位置矢量使连接点处的速度跳变Δv接近于零;对于n段轨迹,有n+1个连接点x0,x1,……xn;相应的微分修正关系式为,
<mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mi>&amp;delta;</mi> <mi>&amp;Delta;</mi> <mi>v</mi> </mtd> </mtr> <mtr> <mtd> <mo>...</mo> </mtd> </mtr> <mtr> <mtd> <mi>&amp;delta;</mi> <mi>&amp;Delta;</mi> <mi>v</mi> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mi>M</mi> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>r</mi> <mn>0</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>t</mi> <mn>0</mn> </msub> </mtd> </mtr> <mtr> <mtd> <mo>...</mo> </mtd> </mtr> <mtr> <mtd> <msub> <mi>r</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>t</mi> <mn>1</mn> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>
即实现通过非线性降维方法和二阶微分修正得到星历模型下精确的拟周期轨道。
5.如权利要求4所述的一种用于复杂约束下低能量小天体精确探测轨道转移方法,其特征在于:步骤四的具体实现方法为,
单值矩阵M为状态转移矩阵φ在一个轨道周期T处的取值,即:
M=Φ(t=T) (9)
单值矩阵的特征值对应的特征向量决定拟周期轨道附近的运动特征,其稳定与不稳定特征向量可决定受到初始扰动时的发散或收敛;拟流形扰动法利用单值矩阵M不稳定特征向量作为扰动量施加的方向,使得探测器远离初始轨道并到达目标小行星附近;施加的扰动速度矢量ΔV表达式为,
ΔV=λu (10)
其中:λ速度扰动量大小,u为不稳定特征向量速度分量的单位向量;针对不同的探测器从拟周期轨道分离的时刻t0,以速度扰动量大小λ为优化变量得到转移轨道与目标小天体相距最近的距离D与相距最近的时刻tf,从而为转移轨道设计初值提供依据;
即探测器从拟周期轨道分离的时刻t0、速度扰动量大小λ为采用拟流形扰动法优化获得转移轨道初值。
6.如权利要求5所述的一种用于复杂约束下低能量小天体精确探测轨道转移方法,其特征在于:步骤五的具体实现方法为,
综合考虑多种复杂非一致强耦合约束,根据有限的速度增量大小对速度扰动量大小λ进行约束;考虑深空测控约束,以步骤四得到的转移轨道初值为基础,设置满足飞行机动状态约束的节点范围,利用二级微分修正,设计不同机动时刻的小行星精确转移轨道,最后考虑飞行机动状态约束,筛选出符合给定时间范围内且燃料消耗最省的转移轨道,即以满足多种复杂非一致强耦合约束为前提,实现低能量小天体精确探测轨道转移。
CN201711264944.XA 2017-12-05 2017-12-05 一种用于复杂约束下低能量小天体精确探测轨道转移方法 Active CN108100307B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711264944.XA CN108100307B (zh) 2017-12-05 2017-12-05 一种用于复杂约束下低能量小天体精确探测轨道转移方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711264944.XA CN108100307B (zh) 2017-12-05 2017-12-05 一种用于复杂约束下低能量小天体精确探测轨道转移方法

Publications (2)

Publication Number Publication Date
CN108100307A true CN108100307A (zh) 2018-06-01
CN108100307B CN108100307B (zh) 2018-12-21

Family

ID=62208108

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711264944.XA Active CN108100307B (zh) 2017-12-05 2017-12-05 一种用于复杂约束下低能量小天体精确探测轨道转移方法

Country Status (1)

Country Link
CN (1) CN108100307B (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110736469A (zh) * 2019-11-06 2020-01-31 北京理工大学 基于日地旋转坐标系的小行星探测精确轨道转移方法
CN110736470A (zh) * 2019-11-06 2020-01-31 北京理工大学 一种不规则形状小天体附近连续推力轨道混合搜索方法
CN110889504A (zh) * 2019-12-05 2020-03-17 北京理工大学 一种基于能量分布评估的多目标小行星交会序列搜索方法
CN111605736A (zh) * 2020-04-29 2020-09-01 北京理工大学 地月l2点转移轨道最优误差修正点选择方法
CN111737814A (zh) * 2020-05-25 2020-10-02 北京理工大学 一种利用地月旁近的平衡点周期轨道逃逸方法
CN111832957A (zh) * 2020-07-21 2020-10-27 中国人民解放军32035部队 一种适用于解体初期解体时刻的选优分析方法
CN112009727A (zh) * 2020-08-21 2020-12-01 北京空间技术研制试验中心 平动点轨道最优小推力转移分段设计方法
CN112660426A (zh) * 2020-12-15 2021-04-16 北京航天自动控制研究所 一种火箭软着陆制导方法
CN113247309A (zh) * 2021-05-20 2021-08-13 中国科学院国家空间科学中心 基于碰撞带图谱的转移轨道族初值搜索方法及系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103112600A (zh) * 2013-03-04 2013-05-22 北京理工大学 一种星际转移轨道设计方法
CN103863584A (zh) * 2014-03-26 2014-06-18 南京航空航天大学 一种小行星变轨方法
CN104369875A (zh) * 2014-10-31 2015-02-25 中国运载火箭技术研究院 基于非线性轨道计算的航天器制导控制方法及系统
CN105329464A (zh) * 2015-09-25 2016-02-17 北京理工大学 一种基于平衡点周期轨道的行星低能量捕获轨道方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103112600A (zh) * 2013-03-04 2013-05-22 北京理工大学 一种星际转移轨道设计方法
CN103863584A (zh) * 2014-03-26 2014-06-18 南京航空航天大学 一种小行星变轨方法
CN104369875A (zh) * 2014-10-31 2015-02-25 中国运载火箭技术研究院 基于非线性轨道计算的航天器制导控制方法及系统
CN105329464A (zh) * 2015-09-25 2016-02-17 北京理工大学 一种基于平衡点周期轨道的行星低能量捕获轨道方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
乔栋: "嫦娥二号卫星飞越Toutatis小行星转移轨道设计", 《中国科学:科学技术》 *

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110736470A (zh) * 2019-11-06 2020-01-31 北京理工大学 一种不规则形状小天体附近连续推力轨道混合搜索方法
CN110736469A (zh) * 2019-11-06 2020-01-31 北京理工大学 基于日地旋转坐标系的小行星探测精确轨道转移方法
CN110889504B (zh) * 2019-12-05 2022-10-21 北京理工大学 一种基于能量分布评估的多目标小行星交会序列搜索方法
CN110889504A (zh) * 2019-12-05 2020-03-17 北京理工大学 一种基于能量分布评估的多目标小行星交会序列搜索方法
CN111605736A (zh) * 2020-04-29 2020-09-01 北京理工大学 地月l2点转移轨道最优误差修正点选择方法
CN111605736B (zh) * 2020-04-29 2021-06-22 北京理工大学 地月l2点转移轨道最优误差修正点选择方法
CN111737814A (zh) * 2020-05-25 2020-10-02 北京理工大学 一种利用地月旁近的平衡点周期轨道逃逸方法
CN111832957A (zh) * 2020-07-21 2020-10-27 中国人民解放军32035部队 一种适用于解体初期解体时刻的选优分析方法
CN111832957B (zh) * 2020-07-21 2024-04-26 中国人民解放军32035部队 一种适用于解体初期解体时刻的选优分析方法
CN112009727A (zh) * 2020-08-21 2020-12-01 北京空间技术研制试验中心 平动点轨道最优小推力转移分段设计方法
CN112660426B (zh) * 2020-12-15 2021-09-14 北京航天自动控制研究所 一种火箭软着陆制导方法
CN112660426A (zh) * 2020-12-15 2021-04-16 北京航天自动控制研究所 一种火箭软着陆制导方法
CN113247309B (zh) * 2021-05-20 2022-04-05 中国科学院国家空间科学中心 基于碰撞带图谱的转移轨道族初值搜索方法及系统
CN113247309A (zh) * 2021-05-20 2021-08-13 中国科学院国家空间科学中心 基于碰撞带图谱的转移轨道族初值搜索方法及系统

Also Published As

Publication number Publication date
CN108100307B (zh) 2018-12-21

Similar Documents

Publication Publication Date Title
CN108100307A (zh) 一种用于复杂约束下低能量小天体精确探测轨道转移方法
CN104527994B (zh) 异面交叉快变轨道固定时间稳定姿态指向跟踪控制方法
CN105631095B (zh) 一种等间隔发射的多约束地月转移轨道簇搜索方法
CN109592079A (zh) 一种限定时间的航天器共面交会变轨策略确定方法
CN105574261B (zh) 一种月球借力约束的地月平动点转移轨道设计方法
CN104345738B (zh) 一种绳系释放稳定控制方法及电动力绳系末级离轨稳定控制方法
CN105912819B (zh) 一种地月l1拉格朗日点转移轨道的快速设计方法
CN102819266B (zh) 一种拟周期j2不变相对轨道编队飞行控制方法
CN105718659B (zh) 一种高面质比航天器轨道动力学分析方法
CN111552003B (zh) 基于球卫星编队的小行星引力场全自主测量系统及方法
CN105930305B (zh) 一种三脉冲交会接近制导方法
CN104309822A (zh) 一种基于参数优化的航天器单脉冲水滴形绕飞轨迹悬停控制方法
Arlulkar et al. Solution based on dynamical approach for multiple-revolution lambert problem
CN108082538B (zh) 一种考虑始末约束的多体系统低能量捕获轨道方法
CN109240340A (zh) 一种基于拟周期轨道的洛伦兹力多星编队构型方法
He et al. Solution set calculation of the Sun-perturbed optimal two-impulse trans-lunar orbits using continuation theory
Shen et al. Point return orbit design and characteristics analysis for manned lunar mission
Kang et al. Nanosat formation flying design for SNIPE mission
Samsam et al. Multi-impulse shape-based trajectory optimization for target chasing in on-orbit servicing missions
Farrés et al. Periodic motion for an imperfect solar sail near an asteroid
Ploen et al. Dynamics of earth orbiting formations
Wang et al. Optimal four-impulse rendezvous between coplanar elliptical orbits
CN112009727B (zh) 平动点轨道最优小推力转移分段设计方法
CN115113638A (zh) 一种燃料最优主动漂移三维成像轨道控制方法
Korolev Problems of spacecraft multi-impulse trajectories modelling

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