经典力学/陀螺进动章动与稳定性

经典力学/陀螺进动章动与稳定性

基本说明

本章介绍刚体陀螺在外力矩作用下的典型运动：进动、章动（nutation）以及与之相关的稳定性判据。我们以“对称重陀螺”（一主惯量轴对称，且质心不在支点）为主线，结合 欧拉方程 与有效势方法，说明：

• 为什么会出现匀速进动（steady precession）；
• 为什么会叠加章动（进动角速度随时间起伏）；
• 什么条件下“几乎无章动”的进动是稳定的（小扰动不发散）。

相关背景可先读：经典力学/欧拉方程与对称陀螺。

前置知识与符号

• 刚体在体坐标系中的主惯量：$I_1=I_2,I_3$（对称陀螺，3轴为对称轴）。
• 欧拉角：$(\varphi ,\theta ,\psi )$，其中 $\theta$ 为陀螺对称轴与竖直方向夹角（倾角）。
• 角速度在体轴分解：$𝝎$。
• 重陀螺：质心到支点距离 $l$，重力势能 $V(\theta )=mgl\cos \theta$（竖直向上为 $\theta =0$ 的约定可按课程统一；若取向下为零，公式会等价改写）。

参见：陀螺仪、刚体动力学、欧拉角。

物理图像：为什么会“转着转着开始绕圈”

当陀螺自转角动量很大时，重力对支点产生的力矩 $𝝉=𝐫\times m𝐠$ 不会立刻把陀螺“压倒”，而是主要改变角动量方向，使其缓慢绕竖直轴旋转——这就是进动。如果初始条件不恰好落在“匀速进动解”上，系统会在匀速进动的基础上出现周期性摆动：$\theta (t)$ 在两个极值之间往返，这就是章动。

动力学推导（对称重陀螺）

角速度与动能

对称陀螺用欧拉角表示角速度，可写成： ${\omega }_3=\dot{\psi }+\dot{\varphi }\cos \theta$ 且 $T=\frac{1}{2}{I}_1({\dot{\theta }}^2+{\dot{\varphi }}^2{\sin }^2\theta )+\frac{1}{2}{I}_3(\dot{\psi }+\dot{\varphi }\cos \theta {)}^2.$

势能取 $V(\theta )=mgl\cos \theta .$

两个守恒量（循环坐标）

由于 $\varphi$ 与 $\psi$ 为循环坐标，对应广义动量守恒：

• 自转轴角动量分量

$p_{\psi }=\frac{\partial L}{\partial \dot{\psi }}=I_3(\dot{\psi }+\dot{\varphi }\cos \theta )\equiv \mathrm{\Lambda }$ 常记作 $L_3$）。

• 空间竖直方向的角动量投影（或与 $p_{\varphi }$ 等价的常数）

$p_{\varphi }=\frac{\partial L}{\partial \dot{\varphi }}=I_1\dot{\varphi }{\sin }^2\theta +I_3(\dot{\psi }+\dot{\varphi }\cos \theta )\cos \theta \equiv \mathrm{\Pi }.$

这些结果与对称陀螺欧拉方程分析一致。

有效势与章动条件

将守恒量代回，可把 $\theta$ 的运动写成一维“粒子”形式： $\frac{1}{2}{I}_1{\dot{\theta }}^2+U_{\mathrm{e}\mathrm{f}\mathrm{f}}(\theta )=E^{\prime },$ 其中有效势（形式可写成） $U_{\mathrm{e}\mathrm{f}\mathrm{f}}(\theta )=mgl\cos \theta +\frac{(\mathrm{\Pi }-\mathrm{\Lambda }\cos \theta {)}^2}{2I_1{\sin }^2\theta }+\frac{{\mathrm{\Lambda }}^2}{2I_3}.$

于是：

• 若 $E^{\prime }$ 与 $U_{\mathrm{e}\mathrm{f}\mathrm{f}}(\theta )$ 只有一个接触点（切点），则 $\theta$ 为常数：无章动的匀速进动。
• 若有两个转折点 ${\theta }_{min},{\theta }_{max}$，则 $\theta$ 在其间往返：出现章动。

章动的“快慢”受 $I_1$、自转强度（$\mathrm{\Lambda }$）、以及倾角附近有效势曲率控制。

匀速进动（steady precession）

匀速进动的条件

匀速进动意味着 $\dot{\theta }=0$ 且 $\ddot{\theta }=0$。在有效势图像中等价于： $\frac{d{U}_{\mathrm{e}\mathrm{f}\mathrm{f}}}{d\theta }=0$ 且该点为极小值（稳定）或极大值（不稳定）。

对重陀螺的经典近似（自转很快、倾角不太大）给出常见匀速进动角速度量级： ${\mathrm{\Omega }}_p\sim \frac{\tau }{L}=\frac{mgl\sin \theta }{I_3{\omega }_3}.$ 这就是“力矩改变角动量方向”的直观比值。

快进动与慢进动（两支解）

更精确地解匀速进动条件时，常会得到关于 $\dot{\varphi }$ 的二次方程，从而出现两支解（常被称为“慢进动”和“快进动”分支）。慢进动分支在高速自转时更常见；快进动分支在某些初始条件/能量范围也可出现，但对扰动更敏感（具体稳定性取决于有效势局部性质）。

章动（nutation）

章动从何而来

如果初始倾角、进动角速度与自转角速度不满足匀速进动的“切触条件”，系统会在有效势井中往返，表现为：

• $\theta (t)$ 周期性变化；
• $\dot{\varphi }(t)$ 与 $\dot{\psi }(t)$ 随之调制；
• 轨迹在空间中呈现“花边”式进动（进动叠加轻微起伏）。

这与 nutation 的一般定义一致。

“睡眠陀螺”（sleeping top）与小章动

当陀螺几乎竖直（$\theta$ 很小）并高速自转时，$U_{\mathrm{e}\mathrm{f}\mathrm{f}}$ 在 $\theta \approx 0$ 附近可形成很陡的势阱，小扰动仅产生小幅振荡（小章动），视觉上像“几乎不抖”的稳定进动/定轴。

稳定性：什么叫“稳定进动”？

本节讨论的是：在某个匀速进动解附近加入小扰动，$\theta$ 是否保持小幅振荡（稳定），还是逐渐偏离（不稳定）。

有效势判据（最常用）

若匀速进动对应的 $\theta ={\theta }_0$ 满足 ${\frac{d{U}_{\mathrm{e}\mathrm{f}\mathrm{f}}}{d\theta }|}_{{\theta }_0}=0$ 且 ${\frac{d^2{U}_{\mathrm{e}\mathrm{f}\mathrm{f}}}{d{\theta }^2}|}_{{\theta }_0}>0,$ 则在能量不变、守恒量固定条件下，该匀速进动对小扰动是稳定的；小扰动导致的小章动频率近似与 ${\omega }_{\mathrm{n}\mathrm{u}\mathrm{t}}^2\approx \frac{1}{I_1}{\frac{d^2{U}_{\mathrm{e}\mathrm{f}\mathrm{f}}}{d{\theta }^2}|}_{{\theta }_0}$ 同阶（精确形式取决于变量选取与是否做小角近似）。

与欧拉方程观点的对应

从欧拉方程出发，也能得到等价的线性化稳定性结论：在稳定情形下，某些扰动模为振荡而非指数增长。可对照阅读相关讲义/课程资料中的“对称欧拉陀螺转动的稳定性”段落。

常见误区

• “看起来残影很圆”不等于稳定：章动幅度可能很小但仍存在。
• “残差/摆动小”与“参数扰动后仍小”是两件事；稳定性要看小扰动演化。
• 仅凭 ${\mathrm{\Omega }}_p\sim \tau /L$ 的量级估计无法判定章动是否发生；是否章动取决于能量与守恒量是否落在“切触”条件上。

实用小结（把公式变成直觉）

• 高速自转使角动量大：力矩主要改变方向 → 产生进动。
• 初始条件不匹配匀速解：倾角在有效势井中往返 → 产生章动。
• 匀速进动是否稳定：看 $U_{\mathrm{e}\mathrm{f}\mathrm{f}}$ 的极值是“碗底”还是“山顶”。（碗底稳定，山顶不稳定。）

练习

1. 对称重陀螺的拉格朗日量写成 $L=T-V$，指出 $\varphi$、$\psi$ 为循环坐标，并写出对应守恒量。
2. 在给定 $\mathrm{\Lambda }$ 与 $\mathrm{\Pi }$ 的条件下，画出（定性）$U_{\mathrm{e}\mathrm{f}\mathrm{f}}(\theta )$ 的形状，说明何时出现两个转折点。
3. 在小角近似 $\theta \ll 1$ 下，对 $U_{\mathrm{e}\mathrm{f}\mathrm{f}}(\theta )$ 做泰勒展开到二阶，给出小章动频率的量级估计。

参见

• 经典力学/欧拉方程与对称陀螺
• 陀螺仪
• 进动
• 章动
• 欧拉角
• 刚体