---
title: "角动量怎样把轨道形状与倾角连起来"
description: "从开普勒能量推导偏心率、倾角和角动量投影的关系，再估算远处伴星施加力矩的时间尺度。"
canonical_url: "https://sandboxphysics.com/zh/stories/secular-orbits/angular-momentum-and-slow-torques/"
language: "zh-CN"
page_type: "story"
alternate_languages:
  en: "https://sandboxphysics.com/stories/secular-orbits/angular-momentum-and-slow-torques/"
last_modified: "2026-10-28"
---

[从远处力矩走到 Kozai–Lidov 轨道循环](https://sandboxphysics.com/zh/stories/secular-orbits/)

AC09 · L1

# 角动量怎样把轨道形状与倾角连起来

Sandbox Physics · 2026-10-28

在圆外轨、双平均的近似下，内轨角动量沿外轨轴的分量保持不变。偏心率增大时，角动量的大小减小；为了保住这个分量，轨道平面也必须转动。能量决定椭圆的总体尺度，角动量则把形状和方向连在一起。

## 开始前需要什么

需要会用向量叉积，知道牛顿引力下的轨道能量和椭圆几何。本章会推导后两章要用的轨道关系。

## 先看行星沿轨道横向运动的多少

考虑质量可以忽略的行星。位置与速度都相对于中心恒星测量，坐标轴本身不转动。两者的叉积叫作比角动量，其中“比”表示已经除以行星质量。它垂直于瞬时轨道平面。

行星经过近心点时，径向速度恰好为零，速度完全沿横向。角动量大小便等于近心距乘以此处速度。把速度写成角动量除以近心距，再代入开普勒能量；最后使用椭圆近心距与半长轴的关系，就能得到下面的恒等式。

固定半长轴时，圆轨道拥有最大的角动量。偏心率达到零点八的椭圆，只剩同尺度圆轨道六成的角动量。行星在远处停留很久，随后迅速掠过近心点；叉积始终记录着这份横向运动。

$$
\boldsymbol h=\boldsymbol r\times\boldsymbol v,\qquad E=\frac{v^2}{2}-\frac{GM}{r}=-\frac{GM}{2a}
$$

位置、速度和能量均相对于中心恒星。能量等式对孤立的开普勒轨道成立。

$$
q=a(1-e),\quad v_q=\frac{h}{q},\quad \frac{h^2}{2q^2}-\frac{GM}{q}=-\frac{GM}{2a}\quad\Longrightarrow\quad h^2=GMa(1-e^2)
$$

在近心点代入速度，便得到偏心率与角动量的关系；这里没有把椭圆速度误当成圆轨道速度。

对应阅读：Tremaine & Yavetz, Why do Earth satellites stay up?, arXiv:1309.5244v1

## 保住沿外轨轴的那一部分

把内外两条轨道上的快速运动都平均掉后，圆外轨的引力作用具有绕轴对称性。让整个内轨椭圆绕外轨轴转动，平均相互作用能不会变化，对应的角动量分量便守恒。

用同半长轴圆轨道的角动量来归一化，所得向量长度就是一减偏心率平方的平方根。互倾角是这个向量与外轨轴的夹角。于是沿外轨轴的投影，等于向量长度乘以倾角余弦。

默认初始偏心率为零点零五、互倾角为六十五度，固定投影约为零点四二二零九。假如偏心率后来升至零点八，角动量长度只剩零点六，互倾角便要降到约四十五点三度。轨道变细长时，平面也随之转动。

单靠投影关系还能得到一个几何上限：向量的总长度不能小于任何分量的绝对值。默认参数下，它允许偏心率接近零点九零七。但实际平均演化的最大值只有约零点八三八九零九；完整的引力哈密顿量还施加了更强的限制，下一章会继续推导。

$$
\boldsymbol j=\frac{\boldsymbol h}{\sqrt{GMa}},\quad j=\sqrt{1-e^2},\quad j_\parallel=\boldsymbol j\cdot\hat{\boldsymbol n}_2=\sqrt{1-e^2}\cos i
$$

轴向由外轨单位法向量定义。这里的倾角是两条轨道之间的夹角，不是相对于观测者天空的倾角。

$$
j_\parallel=\sqrt{1-0.05^2}\cos65^\circ=0.422090,\qquad \cos i=\frac{0.422090}{0.600}
$$

此例采用顺行轨道，并保持圆外轨双平均模型的轴向投影。

$$
e\leq\sqrt{1-j_\parallel^2}
$$

这是必要的几何条件，尚未用尽平均哈密顿量的约束。

对应阅读：Tremaine & Yavetz, Why do Earth satellites stay up?, arXiv:1309.5244v1 · Liu, Muñoz & Lai, Short-range forces in triple systems, arXiv:1409.6717v2

## 变细长与向内迁移是两回事

这个保守平均模型保持半长轴不变。偏心率从零点零五升到零点八时，近心距从半长轴的零点九五倍减至零点二倍，远心距则从一点零五倍增至一点八倍。行星可以靠近恒星，而椭圆的总体尺度仍然相同。

若要让这些近距离经过形成短周期行星，还需要移走轨道能量。潮汐迁移设想中，承担这件事的是潮汐耗散。保守的相对论近心点进动只转动椭圆方向，并不耗散能量。判断迁移需要继续追踪能量与半长轴的长期变化。

直接三体积分中，瞬时密切半长轴也会有小幅变化。密切轨道是用当时的位置和速度拟合的开普勒椭圆，它与平均理论中固定的轨道作用量有所不同。读数据曲线时，要留意采用的是哪一种描述。

$$
q=a(1-e),\qquad Q=a(1+e)
$$

固定半长轴，近心点与远心点便向相反方向移动。

对应阅读：Liu, Muñoz & Lai, Short-range forces in triple systems, arXiv:1409.6717v2 · Gupta et al., A hot-Jupiter progenitor on a super-eccentric retrograde orbit, Nature 632 (2024)

## 估算远处力矩需要多久

伴星对内系统的差分加速度，正比于内轨尺度，反比于外轨距离的三次方。平均摄动势具有下面的尺度。用圆轨道角动量除以这个势，便得到特征世俗时间。

本组计算令引力常数、两颗恒星质量和初始内轨半长轴都为一。内轨平均运动角速度为一，周期约为六点二八个时间单位。外轨距离为十二时，伴星相对于中心恒星的周期约为一百八十四点六八七，世俗时间则为一千七百二十八。默认保存轨迹的一次偏心率循环大约需要五个世俗时间。

外轨相对运动的开普勒时钟使用两颗恒星的质量之和；测试粒子的内轨时钟使用中心恒星质量。忽略这项区别，快相位就会改变，两组直接积分也无法保持一致。扩大外轨距离时，摄动力矩减弱得比外轨周期增长快得多。

当两个公转周期都远短于角动量明显变化的时间时，平均通常更可靠。高偏心率状态下，一次外轨公转就可能显著改变内轨角动量，即使半长轴看起来仍相隔很远。第三章会用实际保存的轨迹检验这个近似。

$$
\Phi_0=\frac{Gm_2a_1^2}{a_2^3},\qquad t_K=\frac{\sqrt{GMa_1}}{\Phi_0}=\frac{1}{n_1}\frac{M}{m_2}\left(\frac{a_2}{a_1}\right)^3
$$

外轨为圆。这个定义固定了本组所有平均方程的时间归一化。

$$
n_1=\sqrt{\frac{GM}{a_1^3}},\quad n_2=\sqrt{\frac{G(M+m_2)}{a_2^3}},\quad \tau=\frac{t}{t_K}
$$

完整循环周期需要另行测量，以特征世俗时间为单位。

对应阅读：Tremaine & Yavetz, Why do Earth satellites stay up?, arXiv:1309.5244v1 · REBOUND 5.2.2 official Kozai example and IAS15 integrator

保存的数值例子与计算方法

合成动力学单位令引力常数、两颗恒星质量和初始内轨半长轴均为一。平均原始向量采用 RK4，最大世俗步长为零点零零一；直接参照为固定 REBOUND 5.2.2 IAS15 状态，保存间隔为零点零零五个世俗时间。

| 初始偏心率与互倾角 | $e_0=0.0500,\quad i_0=65.0^\circ$ |
| --- | --- |
| 外内半长轴比 | $a_2/a_1=12.0$ |
| 内轨周期、外轨周期、世俗时间 | $P_1=6.28319,\quad P_2=184.687,\quad t_K=1728$ |
| 代表性回放时刻 | $t/t_K=1.800$ |
| 偏心率与互倾角 | $e=0.503185,\quad i=60.762649^\circ$ |
| 近心距与轴向角动量 | $q/a_1=0.496815,\quad j_\parallel=0.422090$ |

代表状态取世俗时刻一点八。有限初始偏心率与近圆解析极限分别列明；可下载分析能复算峰值，并提供实际直接积分状态。

## [亲手改变轨道交换](https://sandboxphysics.com/zh/labs/kozai-lidov-cycles/)

1. 先打开低倾角预设，再切换偏心率交换预设，比较曲线振幅和内轨平面的变化。
2. 在偏心率接近峰值时暂停，核对角动量总长度减小而轴向分量保持稳定的关系。

## 练习与解释性答案

1. 偏心率为零点六时，近心距、远心距和角动量长度分别是多少？以相同半长轴的圆轨道为参照。

   展开解释

   近心距是半长轴的零点四倍，远心距是一点六倍；角动量是圆轨道的零点八倍。分别使用椭圆几何和本章从能量推导的恒等式即可得到。

2. 保持中心恒星质量和两个半长轴不变，把伴星质量加倍，所有时间尺度都会减半吗？

   展开解释

   世俗时间减半，因为摄动力矩尺度加倍。外轨周期取决于恒星总质量的平方根；从等质量双星出发，它会乘以三分之二的平方根。内轨开普勒周期保持不变。

3. 默认轨道保留同一个轴向投影，能否在平均演化中达到九十度互倾角？

   展开解释

   不能。默认的正投影约为零点四二二，而九十度会让投影为零。轨道翻转需要另一种动力学条件，或引入能够改变此处约束的物理项。

[上一篇：远处伴星怎样把轨道拉长](https://sandboxphysics.com/zh/stories/secular-orbits/why-a-distant-companion-stretches-an-orbit/) [下一篇：从四极势推导偏心率循环](https://sandboxphysics.com/zh/stories/secular-orbits/derive-the-quadrupole-cycle/)

## 来源与计算范围

1. [Tremaine & Yavetz, Why do Earth satellites stay up?, arXiv:1309.5244v1](https://arxiv.org/abs/1309.5244v1)

   固定为2013年9月20日版本，DOI10.1119/1.4874853。势函数18、向量方程23、条件24和第IV节给出四极模型及近圆轨道临界角极限。

2. [Liu, Muñoz & Lai, Short-range forces in triple systems, arXiv:1409.6717v2](https://arxiv.org/abs/1409.6717v2)

   固定为2015年1月7日版本，DOI10.1093/mnras/stu2396。阅读第2–4节、四极势7、八极势8、GR项32–33/42和图3–4、7。浏览器采用圆外轨、测试粒子四极模型，导读继续阅读更完整的模型。

3. [REBOUND 5.2.2 official Kozai example and IAS15 integrator](https://rebound.hanno-rein.de/c_examples/kozai/)

   官方例程用于选择独立直接积分参照。固定轨迹使用REBOUND 5.2.2、源提交c34134716cd484ed9cdb9f7d54ef8d859bdf547d、IAS15 精度参数为十的负十一次方，初值与停止条件另行列明。软件采用GPL-3.0-only，Lab没有复制发布者代码。

4. [Gupta et al., A hot-Jupiter progenitor on a super-eccentric retrograde orbit, Nature 632 (2024)](https://www.nature.com/articles/s41586-024-07688-3)

   正式版发表于2024年7月17日，DOI10.1038/s41586-024-07688-3，CC BY 4.0。图1–2为凌星和视向速度观测；扩展图5–6讨论依赖模型的初值与示意潮汐迁移史。核读了作者存档全文PMC11291287的XML，分清观测和演化模型。