三维引力辅助解析的分析方法研究.pdfVIP

三维引力辅助解析分析方法研究1 *,2 * 贾建华 ，王琪 * （北京航空航天大学航空科学与工程学院， 北京 100191） 摘 要：针对星际探测中的引力辅助技术，基于二体问题，导出了三维引力辅助前后飞行器速度、能量、角动量和轨道倾角 与引力辅助参数之间的解析公式。将解析公式计算的结果与基于圆型限制性三体问题的数值仿真结果进行了对比，各项绝对 误差小于 0.01，除极少数情况外，相对误差均小于0.7% 。分析结果表明本文导出的解析公式与实际情况对比具有较高的精 度，可以满足引力辅助轨道初步设计的需要。 关键词：三维引力辅助；解析公式；二体问题；限制性三体问题 引 言 引力辅助技术 Gravity-Assist ，GA ；也称为Swing-by 、Fly-by ）可以有效降低星际探测任务所需的发射 能量和总的速度增量，因而引起了广泛关注并已经成为深空探测轨道设计的一项热门研究。从历届国际深 [1] 空轨迹优化竞赛的结果看，排名前7位的队伍全部采用了引力辅助技术 表明该项技术在深空探测中的 ， [2] 地位非常重要。为此，中国已将其列为需要重点发展的星际航行节能技术之一 。 目前，对引力辅助机理的研究主要有两种方法：基于二体问题的解析法和基于限制性三体问题的数值 积分方法。1988 年Broucke 基于二体问题对二维引力辅助前后飞行器的速度、能量和角动量进行了分析， [3] 揭开了人类对引力辅助机理认识的新篇章 。此后，一系列的文献[4-7]都是采用基于限制性三体问题的数 值积分方法对引力辅助机理进行了研究。基于限制性三体问题的数值积分方法虽然求解精度较高，但是难 以找出有关参数的函数关系，不易揭示各量间的内在联系。基于二体问题的解析法，虽然精度稍差，但不 需进行数值计算，就可以给出各量间的关系，而且物理含义明确，特别是对于初步设计是十分方便的工具。 因此，本文采用二体模型，导出三维引力前后飞行器速度、能量、角动量和轨道倾角与引力辅助参数 之间的解析公式。以飞行器利用木星引力辅助为例，将解析公式计算的结果与基于圆型限制性三体问题的 数值仿真结果进行了对比，验证了解析公式的正确性。 1 问题描述 设主要天体M 1 和次要天体M 2 ，以开普勒轨道绕其公共质心C 运动；第三体M 3 （本文中为飞行器） 的质量远远小于两个主天体的质量，在两个主天体的作用下运动，而不影响两个主天体的运动。 本文基于圆锥曲线拼凑模型，即引力辅助变轨过程不是瞬间完成的。此时就需要考虑M 2 的影响球边 1 国家自然科学基金 资助项目 2 E-mail ：01座机电话号码01座机电话号码@ase.buaa.edu.cn 界的限制。M 3 的运动可分为三段：（1）M 3 处于M 2 的影响球之外时，M 3 绕M 1 运动，忽略M 2 的影响。 （2）当M 3 进入M 2 的影响球内时，此时M 3 绕M 2 运动，而忽略M 1 的影响。（3）M 3 离开M 2 的影响球 时，M 3 再次绕M 1 运动，忽略M 2 的影响。 第二段即为本文所要研究的引力辅助阶段。如图 1所示，飞行器从A点进入M 2 的影响球，经过近拱 点 飞向B 点离开M 的影响球。 P 2 图1 三维引力辅助示意图 三维引力辅助可以通过五个自由参数来描述：1 ：M 在 点时相对于M 的速度 点为M 在 V P P p 3 2 3 ? M 的影响球内绕M 运动轨道的近拱点 ；2 ， ：指定M 在三维空间绕M 的运行轨道的近拱点的 2 2 ? 3 2 方向；3 ：M 在近拱点 与M 之间的距离；4 ：速度 同过 点平行于 的平面 平面1 与 r P ? V P x ?y p 3 2 p ? ? ? 过 点和向量 并且垂直于 的平面 平面2 的交线的夹角。? 和 的取值范围为：??[0 ,360 ， P V r 、? ? p p ? ?[?90 ,90 ] ，? ?[0 ,180 ] 。? ? ? ? 本文的目的是：导出引力辅助前后，即飞行器在A 、B 点处的速度、能量、角动量与上述五个引力辅 助参数之间的解析公式，进而导出能量变化和轨道倾角变化的解析公式。 2 三维引力辅助解析公式 为了讨论上的方便，本文建立四个坐标系，如图2 所示。 坐标系1：惯性坐标系xyz 。 以M 1 和M 2 的公共质心为原点；x ?y 平面位于M 1 和M 2 的运动平面； x 轴为t 0 时刻M 和M 的连线。以M 到达 点时刻为t 0 时刻。 P 1 2 3 2 坐标系2：旋转坐标系x y z 。 以M 和M 的公共质心为原点；x ?y 平面位于M 和M