Lagrange三角平动点邻近的相空间结构.doc

天 文 学 报 第 40 卷 第 3 期 1999 年 8 月 Vol . 40 , No . 3 Aug. , 1999 AC TA AS TRONOM ICA SIN ICA L agrange 三角平动点邻近的相空间结构 周礼勇 孙义燧 周济林 ( 南京大学天文系 南京 210093) 摘要 构造描述一类特殊平面圆型限制性三体问题的一个映射 ,并用这个映射讨论了该 类三体问题 L agrange 三角平动点邻近的相空间结构以及它的稳定性 ,发现当两个主天体的 质量比 μ 0 . 02165 时 ,除去 μ = 0 . 01440 的例外情况 ,三角平动点被不变曲线包围 ,是稳定 的 ,这与理论结果相符 ,由此可解释特罗央群和希腊群小行星的稳定存在 . 该结果也验证了映 射方法的可靠性和有效性. 关键词 Lagrange 平动点 , 小行星 , 映射方法 引言 1 三体运动最简单的模型是圆型限制性三体问题. 人们在理论上研究太阳系中小行星 的运动时 ,经常近似地采用这个模型 ,其中把太阳和木星作为两个主天体 ,这是因为木星 轨道靠近小行星轨道 ,木星的质量又远远大于其他行星质量 ,而且太阳 - 木星轨道的偏心 率很小1 . 此模型的运动方程存在 L agrange 三角平动解 ,而实际上此处恰好存在着两群 小行星 ,即特罗央群和希腊群 . 所以研究圆型限制性三体问题三角平动解的稳定性有着实 际意义. 众所周知 ,用映射方法研究动力系统的长期演化在理论上比较简单 ,在数值探索 时能够节省大量的计算时间 ,并且减少计算过程中的误差积累 . 本文研究两个主天体质量 比很小的平面圆型限制性三体问题 ,利用 Hadjidemet rio u2 提出的方法构造了一个映射 , 通过该映射对三角平衡解邻近的相空间结构以及它的稳定性作了讨论 ,并与理论结果作 了比较. 2 运动方程 设质量为 m 的木星和质量为 M 的太阳围绕它们的质心以圆型轨道绕转 ,质量很小 的第三体即小行星在太阳 - 木星轨道确定的平面上运动 . 小行星的运动由圆型限制性三 体问题的方程描述3 ,4 . 我们定义质量 、长度 、时间的单位如下 : 3 1/ 2 G ( M + m ) a S - J / M + m = 1 , aS - J = 1 , = 1 , 1998 - 07 - 03 收到原稿 ,1998 - 09 - 09 收到修改稿  其中 aS - J 指木星 - 太阳相对运动轨道的半长径. 在这种单位制下 , 万有引力常数 G = 1 , 木星 - 太阳相对运动角速度 = 1 , 令 μ = m / ( M + m ) , 则木星质量为 μ, 太阳质量为 1 - μ, 质心与太阳之间的距离为 s1 = μ, 质心与木星之间的距离为 s2 = 1 - μ, 我们建立以 太阳 - 木星的质心为原点 , X 轴位于它们两 者连线上 ,以角速度 = 1 旋转的旋转参考系 XO Y ,如图 1 所示. 此时 , 描述小行星运动的 哈密顿函数写成极坐标形式为 : H = H0 + μH1 . ( 1) 其中 , P2 1 - μ 1 φ P2 H0 = + - - Pφ , ( 2) r r2 2 r 图 1 旋转坐标系 Fig. 1 The rotaing f rame 1 - μ 1 - μ μ μH1 = - - r r1 r2 = 1 - μ - 1 - μ μ - r r2 + μ2 + 2 rμco sφ 1 r2 + ( 1 - μ) 2 - 2 r ( 1 - μ) co sφ co sφ - 2 = μ r2 μ + O . ( 3) r2 + 1 - 2 rco sφ (3) 式中 r1 和 r2 如图 1 所示分别代表小行星与太阳 、木星之间的距离 ,φ 是如图所示的 角度. 因为 μ是小量 ( 对于太阳 - 木星系统 , 量级为 10 - 3 ) , 所以 ( 1) 式中的 μH1 可以认为 是摄动项 ; 而 H0 则是描述小行星围绕主天体质心作开普勒运动的未受摄的哈密顿函数. 我们引入未受摄哈密顿函数的作用量 ———角变量 I r 、Iφ、θr 、θφ ,这些新的共轭变量与 原共轭变量 Pr 、Pφ、r 、φ 之间的关系是 : 1 - μ ( 4) I r = - Pφ + , 2 + 2 ( 1 - μ) Pφ P2 - - r r2 r ( 5) Iφ = Pφ , ( 1 - μ) 2 r2 1 2 θr = - Iφ + 2 ( 1 - μ)