埃农--黑尔斯模型.pdfVIP

埃农--黑尔斯模型 谢保国 罗克章 付雷 我们已经知道，仅有两个星体相互作用的两体问题，牛顿力学能够给出完全的解析解答， 然而天文学家研究的星团是包含有许多星体的星系，即使是只有三个星体的三体问题，就已 经无法得到严格的解析解，甚至令人更为失望。为研究星团中的星体运动，法国天文学家埃 农（M.Henon ）和美国的研究生黑尔斯（C.Heiles ）研究了一个简化的星体势模型。本文在 这里想再现这一问题的研究和计算，并打算对混沌问题作进一步的讨论。 星体势模型是： 2 2 1 2 1 3 V x y +x y +( y − ) （1） 2 3 式中第一项弹性的二维谐振势，后面两项则是三角对称型的扰动势。后来，有人在研究分子 的振动能级时，也引进了这一形式的电势模型。这个势场的等势线，如图1 所示。显然，这 个势场给出的不是有心力，质点在该势场中运动时不存在角动量守恒积分。但由于势函数的 存在，质点所受的力依旧是保守力，机械能仍然是守恒的。质点（星体）在这样的力场中如 何运动？除能量积分外，是否还可能存在其他积分？这正是这个课题探求的目的。 取质点的质量 m=1,则质点的运动微分方程为： ⎧ ∂V − − + x x (1 2y ) ⎪⎪ ∂x ⎨ （2 ） ⎪ ∂V 2 − − − y y y ( x1) ⎪⎩ ∂y 2 2 2 3 0.5 (x +y )+x y−y /3−0.5 = 0 6 4 2 y 0 −2 −4 −6 −6 −4 −2 0 2 4 6 x 图表 1-势能曲线 由于在z 方向没有力的作用，只要初速度没有z 分量，质点将保持在Oxy 平面内运动。 埃农--黑尔斯对非线性运动方程（2 ）作出了著名的计算研究。我们也对此方方程组也编 制了程序。采用龙格—库塔方法，在各种不同的初始条件下求解，并画出相应的庞加莱面上 的相迹点。相迹点在庞加莱面上的分布与质点的能量有密切关系。实际上，由（1）式可见， 能量的大小是与扰动项的大小相关的；质点的能量E 可作为非线性扰动强度的参量。 图2 是质点在不同能量情形下，各种不同初始条件下的相轨迹在庞加莱面上分别留下的相迹 点，其中（a ）,(b)和（c ）三幅图分别对应于三个不同质点总机械能E 的数值，而每一幅图 形又包含了若干个初始条件下的相轨迹交截点，现分述如下： （a ）E ＝1/12 很明显，庞加莱面上的相点可分别连成明确的闭合曲线，每一闭合曲线代表一个环面与庞加 莱面的交截线。这表明在这一能量状况下每一条相轨迹都缠绕在各自的不变环面上，不变环 面的存在，