基于不变流形与引力辅助的低能量行星逃逸轨道设计方法技术

本发明专利技术涉及一种基于不变流形与引力辅助的低能量行星逃逸轨道设计方法，特别适合于利用动平衡点的低能量深空探测任务轨道设计，属于航天器轨道机动技术领域。本方法基于近拱点庞加莱映射的分段匹配，首先通过引入近拱点庞加莱映射得到不变流形在近拱点处的轨道状态集；然后，根据实现行星际转移所需逃逸双曲线超速要求，采用数值迭代方法计算轨道状态集中不变流形每一分支应施加的机动脉冲，通过机动脉冲大小对比可以确定燃料最省逃逸对应的不变流形分支和需要的机动脉冲；具有算法简单、计算效率高等优点，能够为精确设计提供合理可行的初值猜测。

【技术实现步骤摘要】

本专利技术涉及一种，特别适合于利用动平衡点的低能量深空探测任务轨道设计，属于航天器轨道机动

技术介绍
三体动力学系统中动平衡点及附近的周期轨道因其特殊的位置和动力学特性在空间探测任务中具有重要的应用价值，截止目前全球已成功实施了多次动平衡点轨道任务，例如ISEE-3、S0H0、Genesis等。近些年一些学者提出可将动平衡点作为深空探测中继站的构想，这为低能量实现深空探测提供了新的转移方式。由出发星附近动平衡点周期轨道向深空逃逸的轨道设计是该类型转移方式应用的一个重要基础，目前对这一问题的研究较少，且多集中于直接转移或利用不变流形直接逃逸出发行星引力场的方式。由于太阳一行星系统的质量比通常较小，单纯借用不变流形对逃逸探测器飞行轨迹进行改变的效果不明显，无法达到有效降低燃料消耗的作用，同时还大大延长了飞行时间。若同时利用不变流形和行星引力辅助则可以大幅度的降低探测器由周期轨道逃逸出发行星的燃料消耗，达到进行行星际转移所应具备的能量。然而，在三体动力学系统下，探测器在行星附近的轨道动力学非线性很强，采用传统的轨道动力学方法对逃逸轨道进行设计异常困难。如何选择合适的不变流形分支，计算机动点的需要施加的速度脉冲，即发展一种可靠有效的从周期轨道实现行星逃逸的轨道设计方法是当前需要解决的热点问题。在已发展的低能量行星逃逸轨道设计方法中，在先技术 (Alonso G. P. The design of system-to-system transfer arcs using invariantmanifolds in the multi-body problem. Purdue University, Indiana, USA. 2006.),针对日地系统Halo 轨道至日火系统Halo轨道间的脉冲转移轨道设计问题，给出了一种直接借助不变流形的转移轨道设计方法。在地球逃逸阶段，根据Halo轨道不变流形与火星轨道的接近程度选择合适的流形分支，然后根据日心轨道二体能量的变化选取出在不变流形上的机动位置和需要施加的速度脉冲。该方法只是针对于单纯借助不变流形逃逸地球引力场这一问题研究的， 无法适用于结合不变流形与行星引力辅助的低能量逃逸轨道。在先技术(Nakamiya M. , Yamakawa H. and Scheeres D. J. , Yoshikawa M.Interplanetary Transfers Between Halo Orbits:Connectivity Between Escape and Capture Trajectories . Journal ofGuidance,Control and Dynamics, 2010，33(3) : 803-813.)，针对日地Halo轨道向日火Halo轨道转移问题，采用了结合不变流形与行星引力辅助的地球逃逸和火星俘获策略实现转移。该技术中只是对这一类型的转移方式的性质进行了分析，对地球逃逸和火星俘获需要的燃料消耗进行了初步评估，并未给出具体的流形分支选择和速度脉冲计算方法，无法用于实际的低能量行星逃逸轨道的设计问题。
技术实现思路
本专利技术针对现有技术无法实现燃料最省、不变流形分支选择和机动脉冲计算等问题，提出了一种。本方法基于近拱点庞加莱映射的分段匹配，适合于基于不变流形与引力辅助的低能量行星逃逸轨道设计问题。首先通过引入近拱点庞加莱映射得到不变流形在近拱点处的轨道状态集；然后，根据实现行星际转移所需逃逸双曲线超速要求，采用数值迭代方法计算轨道状态集中不变流形每一分支应施加的机动脉冲，通过机动脉冲大小对比可以确定燃料最省逃逸对应的不变流形分支和需要的机动脉冲。本专利技术方法的技术方案是通过如下步骤实现的步骤I，行星周围近拱点计算由于在行星周围近拱点处改变轨道的能力最强，故逃逸点选在近拱点达到消耗较少的燃料完成逃逸的目的。为完整描述不变流形在行星周围的近拱点，引入近拱点庞加莱映射。旋转坐标系下，令权利要求1.，其特征在于包括以下步骤步骤I，逃逸点选在近拱点，采用庞加莱映射计算得到多个行星周围近拱点旋转坐标系下，令x=表示航天器相对于行星的位置矢量，μ为太阳-行星三体模型中的质量比，则航天器处于近拱点的约束表示为2.根据权利要求I所述的，其特征在于选取逃逸轨道时只考虑速度增量最小，则k1=l，k2=0。3.根据权利要求I所述的，其特征在于若航天器在近拱点处角动量为正，则η的取值范围为(η，3π/2)，对应逆行双曲线轨道；若航天器在近拱点处角动量为负，则η的取值范围为(η/2，π)，对应顺行双曲线轨道。4.根据权利要求I所述的，其特征在于所述速度坐标系的原点处于行星质心上，t轴指向双曲超速方向，在轨道面内将t轴顺时针旋转90°得到S轴；所述旋转坐标系的原点位于太阳和行星的共同质心上，X轴由太阳指向行星，y轴指向行星速度方向。全文摘要本专利技术涉及一种，特别适合于利用动平衡点的低能量深空探测任务轨道设计，属于航天器轨道机动
本方法基于近拱点庞加莱映射的分段匹配，首先通过引入近拱点庞加莱映射得到不变流形在近拱点处的轨道状态集；然后，根据实现行星际转移所需逃逸双曲线超速要求，采用数值迭代方法计算轨道状态集中不变流形每一分支应施加的机动脉冲，通过机动脉冲大小对比可以确定燃料最省逃逸对应的不变流形分支和需要的机动脉冲；具有算法简单、计算效率高等优点，能够为精确设计提供合理可行的初值猜测。文档编号B64G99/00GK102923324SQ20121049939公开日2013年2月13日 申请日期2012年11月29日 优先权日2012年11月29日专利技术者尚海滨, 崔平远, 吴伟仁, 王帅, 赵遵辉, 窦强 申请人:北京理工大学本文档来自技高网...

【技术保护点】
基于不变流形与引力辅助的低能量行星逃逸轨道设计方法，其特征在于：包括以下步骤：步骤1，逃逸点选在近拱点，采用庞加莱映射计算得到多个行星周围近拱点：旋转坐标系下，令χ=[x？1+μ,y]表示航天器相对于行星的位置矢量，μ为太阳？行星三体模型中的质量比，则航天器处于近拱点的约束表示为：χχ·T=0,χ·χ·T+χχ··T≥0在航天器所处Lyapunov轨道上均匀取点，计算各点接近地球的不稳定流形所对应的摄动；给定逃逸轨道的最大时间，对施加摄动后各个点的状态进行时间积分，当满足近拱点的约束时记录该点为近拱点；步骤2，计算近拱点逃逸所需速度增量根据近拱点状态与所需双曲超速，将不变流形近拱点的位置和速度矢量转换至速度坐标系中，建立速度坐标系下的逃逸双曲轨道所满足的非线性方程：g=||r1||-μP(e2-1)||v∞||2(1+ecosf)=0其中：r1为探测器在逃逸位置，v∞为最终达到的双曲超速；μP为行星的引力常数，e为双曲线轨道的偏心率，令逃逸双曲线近拱点处的位置矢量和速度矢量分别为r0和v0，双曲超速与r0和r1的夹角分别为η和α，取值范围均为[0,2π]，f=η？α为r1在双曲线轨道上对应的真近点角；通过牛顿迭代法求得η，迭代过程为η=η-(dgdη)-1g其中:dgdη=μP||v∞||2(cosfdedη-esinfdfdη)-2ededη||v∞||4(1+ecosf)2dedη=sinηcos2η,dfdη=1在取值范围内任意选取η的初值，当迭代的结果不在η的取值范围内时更改η的取值，重复迭代过程直至满足任务设计所需的精度条 件；利用牛顿迭代得到η后，求解探测器速度v1v1=μPa(e2-1)[-sinf·p+(cosf+e)·q]其中a=μ/||v∞||2，p和q分别为双曲线近拱点处位置和速度方向单位矢量；令r1处初始的速度为vi，则逃逸的速度增量Δv为Δv=v1？vi再将速度脉冲Δv转换至旋转坐标系中；步骤3，选取逃逸轨道采用步骤2的方法计算步骤1得到的所有近拱点处逃逸所需的速度增量，选取速度增量的大小和到达近拱点的时间t作为逃逸轨道选取的目标函数：f=k1|Δv|+k2t→min其中k1和k2为参数的权重；利用目标函数选取最佳的逃逸近拱点及相应的速度脉冲，完成基于不变流形与引力辅助的低能量行星逃逸轨道设计。FDA00002496847300016.jpg...

【技术特征摘要】

【专利技术属性】
技术研发人员：尚海滨，崔平远，吴伟仁，王帅，赵遵辉，窦强，
申请(专利权)人：北京理工大学，
类型：发明
国别省市：