贾尼别科夫效应的稳定性判据推导（网球拍定理 / 中间轴不稳定性）

0. 摘要

对无外力矩自由刚体，若其主惯量满足$$I_1 < I_2 < I_3$$I1​<I2​<I3​

则绕 最小主惯量轴 $I_1$I1​ 与 最大主惯量轴 $I_3$I3​ 的匀速自旋是线性稳定的；绕 中间主惯量轴 $I_2$I2​ 的匀速自旋是线性不稳定的。
不稳定性的直接原因来自欧拉方程线性化后出现的“正实指数”增长模态。

---

1. 模型与假设

1.1 假设

• 刚体为理想刚体，质量分布固定
• 在惯性空间中无外力矩：$\vec{\tau}=\vec{0}$τ=0
• 在刚体主轴系（principal axes）中分析
• 只做线性稳定性分析（小扰动）

1.2 变量与符号

• $I_1, I_2, I_3$I1​,I2​,I3​：刚体关于三条主轴的主转动惯量（常数）
• $\boldsymbol{\omega}=(\omega_1,\omega_2,\omega_3)$ω=(ω1​,ω2​,ω3​)：角速度在主轴系的分量
• $\boldsymbol{L}=(L_1,L_2,L_3)$L=(L1​,L2​,L3​)：角动量在主轴系的分量
• 上点 $\dot{(\ )}$( )˙​：对时间求导
• $\Omega$Ω：匀速自旋的主轴角速度标量（常数）
• $\delta\omega_i$δωi​：小扰动量

---

2. 基本方程：欧拉刚体方程（无外力矩）

在主轴系中，惯量张量对角化：$\mathbf{I}=\mathrm{diag}(I_1,I_2,I_3)$I=diag(I1​,I2​,I3​)，且$$L_i = I_i \omega_i$$Li​=Ii​ωi​

无外力矩的欧拉方程为：$$I_1 \dot{\omega}_1 = (I_2 – I_3)\,\omega_2\omega_3$$I1​ω˙1​=(I2​−I3​)ω2​ω3​ $$I_2 \dot{\omega}_2 = (I_3 – I_1)\,\omega_3\omega_1$$I2​ω˙2​=(I3​−I1​)ω3​ω1​ $$I_3 \dot{\omega}_3 = (I_1 – I_2)\,\omega_1\omega_2$$I3​ω˙3​=(I1​−I2​)ω1​ω2​

---

3. 平衡解：绕某主轴的匀速自旋

考虑绕第 $k$k 条主轴匀速自旋。以绕 $I_2$I2​ 轴为例，取平衡解：$$\omega_2=\Omega,\quad \omega_1=0,\quad \omega_3=0$$ω2​=Ω,ω1​=0,ω3​=0

代入欧拉方程右端可得 $\dot{\omega}_i=0$ω˙i​=0，因此这是一个平衡（匀速）解。

同理，绕 $I_1$I1​ 或 $I_3$I3​ 的匀速自旋也分别给出平衡解。

---

4. 线性化：对小扰动建立判别式

以绕 $I_2$I2​ 轴自旋为例，在平衡解附近引入小扰动：$$\omega_2=\Omega+\delta\omega_2,\quad \omega_1=\delta\omega_1,\quad \omega_3=\delta\omega_3$$ω2​=Ω+δω2​,ω1​=δω1​,ω3​=δω3​

线性化原则：只保留一阶小量，忽略二阶项（如 $\delta\omega_1\delta\omega_3$δω1​δω3​）。

4.1 对 $\dot{\omega}_1$ω˙1​ 线性化

$$I_1 \dot{\omega}_1 = (I_2 – I_3)\omega_2\omega_3 \approx (I_2 – I_3)\,(\Omega)\,(\delta\omega_3)$$I1​ω˙1​=(I2​−I3​)ω2​ω3​≈(I2​−I3​)(Ω)(δω3​)

因此：$$\dot{\delta\omega_1}=\frac{I_2-I_3}{I_1}\,\Omega\,\delta\omega_3$$δω1​˙​=I1​I2​−I3​​Ωδω3​

4.2 对 $\dot{\omega}_3$ω˙3​ 线性化

$$I_3 \dot{\omega}_3 = (I_1 – I_2)\omega_1\omega_2 \approx (I_1 – I_2)\,(\delta\omega_1)\,(\Omega)$$I3​ω˙3​=(I1​−I2​)ω1​ω2​≈(I1​−I2​)(δω1​)(Ω)

因此：$$\dot{\delta\omega_3}=\frac{I_1-I_2}{I_3}\,\Omega\,\delta\omega_1$$δω3​˙​=I3​I1​−I2​​Ωδω1​

4.3 $\delta\omega_2$δω2​ 的一阶行为

$$I_2 \dot{\omega}_2=(I_3-I_1)\omega_3\omega_1 \approx (I_3-I_1)\delta\omega_3\delta\omega_1 \approx 0$$I2​ω˙2​=(I3​−I1​)ω3​ω1​≈(I3​−I1​)δω3​δω1​≈0

所以一阶近似下：$$\dot{\delta\omega_2}\approx 0$$δω2​˙​≈0

这表示：线性层面上 $\omega_2$ω2​ 保持常数（扰动主要耦合在 $\delta\omega_1,\delta\omega_3$δω1​,δω3​ 平面内）。

---

5. 得到二阶方程与增长率（稳定性判据的核心）

对 $\dot{\delta\omega_1}$δω1​˙​ 再求导，并代入 $\dot{\delta\omega_3}$δω3​˙​：$$\ddot{\delta\omega_1} =\frac{I_2-I_3}{I_1}\Omega\,\dot{\delta\omega_3} =\frac{I_2-I_3}{I_1}\Omega\left(\frac{I_1-I_2}{I_3}\Omega\,\delta\omega_1\right)$$δω1​¨​=I1​I2​−I3​​Ωδω3​˙​=I1​I2​−I3​​Ω(I3​I1​−I2​​Ωδω1​)

整理得：$$\ddot{\delta\omega_1} =\underbrace{\left(\frac{(I_2-I_3)(I_1-I_2)}{I_1 I_3}\Omega^2\right)}_{\lambda}\,\delta\omega_1$$δω1​¨​=λ(I1​I3​(I2​−I3​)(I1​−I2​)​Ω2)​​δω1​

同理可得 $\ddot{\delta\omega_3}=\lambda\,\delta\omega_3$δω3​¨​=λδω3​。

5.1 判别逻辑

• 若 $\lambda<0$λ<0：解为正弦/余弦振荡（有界），线性稳定
• 若 $\lambda>0$λ>0：解包含 $e^{\sqrt{\lambda}t}$eλ​t 指数增长项（无界），线性不稳定

因此稳定性取决于：$$\lambda \propto (I_2-I_3)(I_1-I_2)$$λ∝(I2​−I3​)(I1​−I2​)

---

6. 结论：中间轴不稳定（完整判据）

假设排序 $I_1<I_2<I_3$I1​<I2​<I3​，则：

• $I_2-I_3 < 0$I2​−I3​<0
• $I_1-I_2 < 0$I1​−I2​<0
• 两者乘积 $(I_2-I_3)(I_1-I_2) > 0$(I2​−I3​)(I1​−I2​)>0

从而：$$\lambda>0 \quad\Rightarrow\quad \text{绕 } I_2 \text{ 轴自旋线性不稳定}$$λ>0⇒绕 I2​ 轴自旋线性不稳定

这就是中间轴不稳定性，对应实验上看到的周期性翻转（贾尼别科夫效应 / 网球拍定理）。

---

7. 对称推导：绕 $I_1$I1​ 与 $I_3$I3​ 的稳定性

7.1 绕 $I_1$I1​ 轴自旋

平衡：$\omega_1=\Omega,\omega_2=0,\omega_3=0$ω1​=Ω,ω2​=0,ω3​=0。线性化得到：$$\ddot{\delta\omega_2} =\left(\frac{(I_3-I_1)(I_1-I_2)}{I_2 I_3}\Omega^2\right)\delta\omega_2$$δω2​¨​=(I2​I3​(I3​−I1​)(I1​−I2​)​Ω2)δω2​

在 $I_1<I_2<I_3$I1​<I2​<I3​ 下：

• $I_3-I_1>0$I3​−I1​>0
• $I_1-I_2<0$I1​−I2​<0
  乘积为负，故 $\lambda<0$λ<0 ⇒ 稳定。

7.2 绕 $I_3$I3​ 轴自旋

平衡：$\omega_3=\Omega,\omega_1=0,\omega_2=0$ω3​=Ω,ω1​=0,ω2​=0。线性化得到：$$\ddot{\delta\omega_1} =\left(\frac{(I_2-I_3)(I_3-I_1)}{I_1 I_2}\Omega^2\right)\delta\omega_1$$δω1​¨​=(I1​I2​(I2​−I3​)(I3​−I1​)​Ω2)δω1​

在 $I_1<I_2<I_3$I1​<I2​<I3​ 下：

• $I_2-I_3<0$I2​−I3​<0
• $I_3-I_1>0$I3​−I1​>0
  乘积为负，故 $\lambda<0$λ<0 ⇒ 稳定。

---

8. 最终稳定性判据（可直接引用）

令 $I_1<I_2<I_3$I1​<I2​<I3​，自由刚体绕主轴匀速自旋的线性稳定性为：

• 绕 $I_1$I1​（最小主惯量轴）：稳定
• 绕 $I_2$I2​（中间主惯量轴）：不稳定
• 绕 $I_3$I3​（最大主惯量轴）：稳定

更一般地，绕第 $k$k 轴自旋时，线性化二阶系数的符号由$$(I_k – I_i)(I_k – I_j)$$(Ik​−Ii​)(Ik​−Ij​)

决定（$\{i,j,k\}=\{1,2,3\}${i,j,k}={1,2,3}）。若该乘积为正则不稳定，为负则稳定；中间惯量轴会使两项同号从而乘积为正，导致指数增长。

---

9. 物理含义的一句话版本

自由刚体同时受两条守恒约束（角动量与动能），在主轴自旋附近的相空间结构使得最大/最小惯量轴对应“能量极值”（稳定），而中间惯量轴对应“鞍点”（不稳定），小扰动会被放大，进而出现翻转。