一种计算圆形限制性三体问题周期轨道的新方法.pdfVIP

中国科学: 技术科学 2011 年 第41 卷 第8 期: 1078 ~ 1083 SCIENCE CHINA PRESS 论 文 一种计算圆形限制性三体问题周期轨道的新方法 * 张汉清 , 李言俊, 张科 西北工业大学航天学院, 西安 710072 * E-mail: hanqing.leben@ 收稿日期: 2010-12-23; 接受日期: 2011-04-15 国家自然科学基金(批准号:和国家重点基础研究发展计划(“863”计划)(批准号: G9KY1004)资助项目 摘要 周期轨道是理解圆形限制性三体问题动力系统的关键, 在深空探测实践中 关键词 具有重要的应用价值, 现有的周期轨道计算方法通常需要以高阶近似解析解作为迭 圆型限制性三体问题 代初值, 计算过程复杂且能得到的周期轨道种类较少. 本文利用圆形限制性三体问题 平动点 周期轨道 的对称性, 通过构造相空间内的一类流函数, 将位于对称平面上的状态再次映射到该 数值方法 对称面上, 得到了计算周期轨道的构造流函数方法, 该方法不需要以近似解析解作为 共振轨道 迭代初值, 也不需要计算状态转移矩阵, 从而便于编程实现, 其优点是不受圆形限制 性三体问题非线性的影响, 并可以系统地计算一大类具有x-z 平面对称性的周期轨道. 作为算例, 在平面情况下计算了 Lyapunov 轨道等周期轨道, 并进一步将构造流函数 方法与微分校正方法结合, 设计了能在不同共振频率上跳转的共振引力加速轨道, 最 后, 在三维空间情况下计算了晕轨道和多种三维周期轨道, 证明了该方法的有效性. 圆形限制性三体问题(CRTBP)动力学模型对于 共线平动点附近晕轨道的 Richardson 三阶解析展开 深空探测具有重要意义, 该模型下存在的共线平动 式[7], Gomez 和 Jorba 等人[8,9]利用 Lindstedt- Poincare 点区域具有特殊的优良性质[1,2], 位于该区域内的探 方法和 Fourier 级数展开得到了计算晕轨道的高阶展 测器不受大气阻力影响, 没有空间碎片、原子氧等侵 开式. Farquhar, Kamel, Richardson 和Gomez 所给出的 袭, 空间物资运输方便, 且非常适于进一步的深空探 方法虽然能得到晕轨道的近似解析解, 但其精度无 测, 是建立太空基地和星际航行港的最佳选择. 而周 法满足实际任务的要求, Breakwell 和 Brown 提出的 期轨道是理解 CRTBP 动力系统的关键, 在深空探测 方法能获得更精确的晕轨道数值解, 但需要以近似 实践中也通常作为探测器的工作轨道[3,4], 具有重要 解析解作为迭代初值, 并需要同时积分包含 42 个方 的应用价值. 程的微分方程组而得到状态转移矩阵, 因此计算量 Farquhar 和 Kamel[5] 研究了太阳引力下周期及 较大. 拟周期轨道的高阶解析式(Farquhar-Kamel 展开式), 本文通过构造一类流函数并利用 CRTBP 的对称 Breakwell 和 Brown[6]