第 4–5 章 · 刚体运动学与动力学

刚体 = 任意两点距离不变的质点系，自由度骤降到 6（质心 3 + 姿态 3）。姿态的描述（欧拉角）与惯性的编码（惯量张量）是刚体力学的两大基建，本章把这两块砖砌好，再拿陀螺当试金石。

欧拉角与角速度

姿态三个角：$\phi$（进动，绕固定 $z$）、$\theta$（章动，绕节线）、$\psi$（自转，绕本体对称轴）。角速度在本体坐标系分解：

写成三个欧拉角速度的线性组合（具体系数表是必背基本功）。为什么本体坐标系？——惯量张量在本体系中是常数矩阵，动力学方程大幅简化；代价是坐标系数表本身是时间函数，这个"扭"由欧拉方程的附加项吸收。

惯量张量与主轴

$\mathbf{I}$ 是实对称矩阵 ⟹ 谱定理保证存在主轴系（对角化）。对称陀螺 $I_1 = I_2 \ne I_3$；球陀螺三个全等。平行轴与转向轴定理是手工计算的省钱工具。角动量与角速度一般不共线（$\mathbf{L} = \mathbf{I}\boldsymbol{\omega}$ 是矩阵乘法）——这就是刚体反直觉行为的根源。

欧拉方程与两大现象

本体系中的转动方程（$\mathbf{I}\dot{\boldsymbol{\omega}} + \boldsymbol{\omega}\times\mathbf{I}\boldsymbol{\omega} = \mathbf{N}$）：

• 自由对称陀螺：$\boldsymbol{\omega}$ 绕对称轴规则进动（无外力也进动！）——地球的极移、自由陀螺的定向；
• 重对称陀螺（拉格朗日陀螺）：拉格朗日量 $L = \frac{I_1}{2}(\dot\theta^2 + \dot\phi^2\sin^2\theta) + \frac{I_3}{2}(\dot\psi + \dot\phi\cos\theta)^2 - mgl\cos\theta$，两个循环坐标（$\phi,\psi$）给两个守恒量，有效势方法解出慢进动与快进动——回转仪导航的原理。
• 网球拍定理（中间轴不稳定）：绕主轴的自由转动中，最大与最小惯量轴稳定、中间轴不稳定——欧拉方程线性化的直接推论，太空翻转的手感来源。

老师的提醒

• 欧拉角速度分解表（$\omega_1, \omega_2, \omega_3$ 用 $\dot\phi, \dot\theta, \dot\psi$ 表出）必须亲手推一遍——所有陀螺习题的第一步。
• "角动量守恒但角速度改变"不再矛盾：$\mathbf{I}$ 在空间系里随姿态变，$\mathbf{L}$ 不变 ⟹ $\boldsymbol\omega$ 得转。

自检

• 写出对称陀螺的欧拉方程并解释自由进动。
• 为什么主轴系让惯量张量对角？与谱定理的关系？
• 网球拍定理的稳定性分析思路？

更深一层：惯量张量是"方向的惯性签名"

惯量张量 $\mathbf{I}$ 的深层身份：把角速度翻译成角动量的线性机器，且这个机器随你在体内选的方向不同而不同。主轴系是它的"出厂说明书"——三个特征值（主惯量 $I_1 \le I_2 \le I_3$）是机器的三个"方向灵敏度"。三条工程含义：

• 旋转稳定性排序：绕最大或最小主轴的自由旋转稳定，绕中间轴不稳定（网球拍定理）——欧拉方程线性化的直接推论；
• 能量-角动量预算：同样角动量下，绕大惯量轴转（扁陀螺）储存的能量更少（$T = \frac{L^2}{2I}$）——耗散系统（带液体的卫星）会自发漂移到绕最大惯量轴转（Explorer-1 卫星的姿态翻转事故——柔性天线耗散能量）；
• 对称陀螺：$I_1 = I_2$ 时进动退化规则——回转仪定向、自行车轮陀螺的原理。

推导细读：欧拉方程的"本体系智慧"

欧拉方程在本体（随体）坐标系里写：$I_1\dot\omega_1 + (I_3 - I_2)\omega_2\omega_3 = N_1$（循环）。为什么在本体系？惯量张量在体内是常数、在空间里随姿态变——常系数方程好解，代价是多出 $\boldsymbol\omega\times\mathbf{I}\boldsymbol\omega$ 项（坐标系的"欺骗项"）。这个"坐标系换得便宜"的权衡是刚体动力学的核心设计决策，也是流体力学随体导数、连续介质力学的随体框架的同一思想：数学描述应该"住"在最自然的坐标系里。

自由对称陀螺的解直接可读：$\omega_3$ 常数（对称轴方向分量）、$\omega_1, \omega_2$ 以频率 $\Omega = \frac{I_3 - I_1}{I_1}\omega_3$ 绕轴转——无外力也进动。视觉版：抛起的瓷砖/手机绕中间轴翻滚（稳定性翻转的实拍视频在网上极多）。

在真实世界里：刚体动力学的工程出口

• 卫星姿控：自旋稳定（绕最大惯量轴自旋，像陀螺保持定向）与三轴稳定（反作用轮 + 磁力矩器）的选择题；
• 手机 IMU：姿态估计（四元数 + 陀螺仪积分）的算法内核就是刚体运动学——欧拉角奇异（万向锁）问题由四元数解决；
• 航天器翻滚救援：失效卫星绕中间轴翻滚的救援时机计算——欧拉方程的周期性 Analysis 直接输出"何时角速度最小适合抓捕"。

专题深入：惯量张量的"方向惯性签名"完整版

惯量张量 $\mathbf{I}$ 的深层身份值得展开：把角速度翻译成角动量的线性机器，且这个机器随你在体内选的方向不同而不同。主轴系是它的"出厂说明书"——三个特征值（主惯量 $I_1 \le I_2 \le I_3$）是机器的三个"方向灵敏度"。三条工程含义：

• 旋转稳定性排序：绕最大或最小主轴的自由旋转稳定，绕中间轴不稳定（网球拍定理）——欧拉方程线性化的直接推论；
• 能量-角动量预算：同样角动量下，绕大惯量轴转（扁陀螺）储存的能量更少（$T = \frac{L^2}{2I}$）——耗散系统（带液体的卫星）会自发漂移到绕最大惯量轴转（Explorer-1 卫星的姿态翻转事故——柔性天线耗散能量）；
• 对称陀螺：$I_1 = I_2$ 时进动退化规则——回转仪定向、自行车轮陀螺的原理。

数学细节展开：欧拉方程的"本体系智慧"完整版

欧拉方程在本体（随体）坐标系里写：$I_1\dot\omega_1 + (I_3 - I_2)\omega_2\omega_3 = N_1$（循环）。为什么在本体系？惯量张量在体内是常数、在空间里随姿态变——常系数方程好解，代价是多出 $\boldsymbol\omega\times\mathbf{I}\boldsymbol\omega$ 项（坐标系的"欺骗项"）。这个"坐标系换得便宜"的权衡是刚体动力学的核心设计决策，也是流体力学随体导数、连续介质力学的随体框架的同一思想：数学描述应该"住"在最自然的坐标系里。

自由对称陀螺的解直接可读：$\omega_3$ 常数（对称轴方向分量）、$\omega_1, \omega_2$ 以频率 $\Omega = \frac{I_3 - I_1}{I_1}\omega_3$ 绕轴转——无外力也进动。视觉版：抛起的瓷砖/手机绕中间轴翻滚（稳定性翻转的实拍视频在网上极多）。

工程案例细节：刚体动力学的工程出口（工程版）

• 卫星姿控：自旋稳定（绕最大惯量轴自旋，像陀螺保持定向）与三轴稳定（反作用轮 + 磁力矩器）的选择题；
• 手机 IMU：姿态估计（四元数 + 陀螺仪积分）的算法内核就是刚体运动学——欧拉角奇异（万向锁）问题由四元数解决；
• 航天器翻滚救援：失效卫星绕中间轴翻滚的救援时机计算——欧拉方程的周期性 Analysis 直接输出"何时角速度最小适合抓捕"。