跳到主要内容
Sandbox Physics

M056 · 自由刚体 / 中间轴不稳定

网球拍定理

分别绕三个主轴释放不对称刚体。旋转三维视角、设置微小横向自转,比较半周翻转与稳定旋转。跟随体轴角速度,并在角动量球与能量椭球的交线上观察运动。

交互模型网球拍定理
实际经过时间Pending\text{Pending}
第一主惯量Pending\text{Pending}
第二主惯量Pending\text{Pending}
第三主惯量Pending\text{Pending}
第一体轴角速度Pending\text{Pending}
第二体轴角速度Pending\text{Pending}
第三体轴角速度Pending\text{Pending}
所选轴偏离角Pending\text{Pending}
本体转动动能Pending\text{Pending}
总角动量大小Pending\text{Pending}
设置状态Pending\text{Pending}
局部增长或摆动速率Pending\text{Pending}
最大归一化账本误差Pending\text{Pending}
两步长姿态角差Pending\text{Pending}
世界角动量误差Pending\text{Pending}

物理教程

网球拍定理

背景在相同形状和扰动下比较三种主轴释放。

为什么值得理解自由物体翻转时,能量和世界角动量依然可以守恒。

先抓重点

核心问题
哪个主轴会放大微小横向自转?
一句话直觉
中间惯量产生局部指数增长。非线性运动反复翻转体轴,但不改变守恒的世界矢量。

核心数学模型

主惯量

I1=m(b2+c2)12,I2=m(a2+c2)12,I3=m(a2+b2)12I_1=\frac{m(b^2+c^2)}{12},\quad I_2=\frac{m(a^2+c^2)}{12},\quad I_3=\frac{m(a^2+b^2)}{12}

尺寸变化可重排惯量大小;编号始终固定在刚体上。

欧拉方程

I1ω˙1=(I2−I3)ω2ω3and cyclic permutationsI_1\dot\omega_1=(I_2-I_3)\omega_2\omega_3\quad\text{and cyclic permutations}

没有力矩不代表体轴角速度恒定;体坐标系本身在旋转。

世界角动量

L=RIωb,L˙=0\mathbf L=R I\boldsymbol\omega_b,\qquad\dot{\mathbf L}=\mathbf 0

惯性系矢量固定,但在运动体坐标中的分量改变。

角动量球与能量椭球

∑iLb,i2=∣L∣2,∑iLb,i2Ii=2E\sum_i L_{b,i}^2=|\mathbf L|^2,\qquad\sum_i\frac{L_{b,i}^2}{I_i}=2E

两个约束都使用体坐标角动量而非角速度;青色轨迹位于交线。

中间轴局部增长

λ=∣Ω∣(Ij−Imin⁡)(Imax⁡−Ij)Imin⁡Imax⁡\lambda=|\Omega|\sqrt{\frac{(I_j-I_{\min})(I_{\max}-I_j)}{I_{\min}I_{\max}}}

仅在所选轴惯量严格居中时有此正增长率;它是无穷小估计,不是翻转频率。

稳定轴振荡

ν=∣Ω∣(Ij−Ik)(Ij−Iℓ)IkIℓ\nu=|\Omega|\sqrt{\frac{(I_j-I_k)(I_j-I_\ell)}{I_k I_\ell}}

最小与最大惯量时乘积为正。重复惯量消去指数增长预测,并不意味着耗散收敛。

常见难点

无力矩就不会翻转

容易误解没有外力矩时,角速度必须固定。

正确理解世界角动量固定;各向异性惯量张量随物体旋转,会改变角速度。

第二轴总是中间轴

容易误解轴编号决定稳定性。

正确理解编辑尺寸后读取惯量;编号固定在刚体上。

动手实验

  1. 01

    比较三个轴

    选择最小、中间和最大惯量预设,再回看完整记录。

    应该观察到: 中间轴释放会翻转所选轴;其余两种偏离保持很小。
  2. 02

    扰动精确解

    比较精确无扰动与极小扰动预设;拖动橙色幅度手柄。

    应该观察到: 无扰动时不稳定解仍可精确保持;更小扰动推迟离开该状态。
  3. 03

    跟随几何交线

    在翻转附近暂停,观察体轴角速度与角动量几何。

    应该观察到: 体坐标分量沿两张守恒曲面运动,世界角动量保持固定。
  4. 04

    改变几何

    令两个或三个尺寸相等,再比较质量缩放与步长误差。

    应该观察到: 重复惯量改变局部增长定律;质量缩放改变能量和角动量而不改变姿态。