M055 · 旋转刚体受力矩

陀螺进动实验台

为什么旋转飞轮会进动，而不只是向下倒？ 因为重力对飞轮的力矩指向侧面，不是向下。力矩朝哪个方向，角动量就往哪个方向变。飞轮的自转角动量很大，受到这个侧向力矩，只会绕着竖直方向转，不会倒下去。整个过程中，能量、竖直方向的角动量和沿轴的自转都保持不变。

飞轮由一个固定支点卡住，可以直接释放，也可以从慢或快的稳态进动开始，还能反转自转、关掉重力。在三维中跟着角动量和力矩箭头看，量一量章动，再和慢进动估算对比。

回放时刻 $\text{Pending}$

轴线倾角 $\text{Pending}$

瞬时进动速率 $\text{Pending}$

章动速率 $\text{Pending}$

实际轴向自转 $\text{Pending}$

慢进动近似 $\text{Pending}$

精确慢分支 $\text{Pending}$

实际设置初态 $\text{Pending}$

转动动能 $\text{Pending}$

重力势能 $\text{Pending}$

总机械能 $\text{Pending}$

竖直角动量 $\text{Pending}$

轴向角动量 $\text{Pending}$

重力力矩 $\text{Pending}$

支点反力大小 $\text{Pending}$

竖直支点反力 $\text{Pending}$

最大归一化能量误差 $\text{Pending}$

最大归一化角动量误差 $\text{Pending}$

两步长轴线差异 $\text{Pending}$

竖直自转阈值 $\text{Pending}$

自转与竖直阈值之比 $\text{Pending}$

设置、测量，再验证

设置初始角动量

显示物理矢量每种箭头的长度各自归一化；真实大小看读数。刚体质量

m

1.00\,\mathrm{kg}

刚体半径

R

0.20\,\mathrm{m}

支点至质心距离

\ell

0.22\,\mathrm{m}

实际轴向自转

\omega_3

100.00\,\mathrm{rad/s}

初始轴线倾角

\theta_0

75.00\,\mathrm{{}^\circ}

普通释放进动速率

\dot\phi_0

0.00\,\mathrm{rad/s}

初始章动速率

\dot\theta_0

0.00\,\mathrm{rad/s}

初态设置

s

0.00\,

记录时长

T

10.00\,\mathrm{s}

均匀重力

g

9.81\,\mathrm{m/s^2}

回放比例

t/T

0.00\,\%

拖动下方的橙色倾角手柄。初态设置决定起点：0 是普通释放，1 是慢速稳态进动分支，2 是快速分支。章动扰动作用在慢分支上；快分支始终零扰动，按解析解运动。如果不存在实数稳态，就按普通释放的参数设置，并在状态栏里注明。

飞轮是均匀薄圆盘，装在无质量的轴上；距离从固定支点量到盘心。画面里的盘厚、支架和轮毂只是示意，不计质量。模型是均匀重力场中一点固定的理想轴对称刚体，直接积分完整方程，不用小倾角或快自转近似。这里的自转指角速度沿转轴的分量，不同于欧拉自转角的变化率。轴承能推也能拉，所以这不是桌面上自由滑动或弹跳的尖端。摩擦、空气阻力、形变、与地面或支架的碰撞以及自转衰减都不考虑。完整记录图里包括提前算好的后续样本。转得快时，屏幕帧率会造成混叠，画面看起来不对；可以暂停，用精细单步查看相位。直立稳定阈值只适用于竖直旋转附近的无穷小倾斜。

物理教程

陀螺进动实验台

背景飞轮一边自转，一边由一个固定支点卡住，三维视角可以随意转动。你可以直接释放它，让它从慢或快的稳态进动开始，也可以反转自转、关掉重力。跟着角动量和力矩箭头看，量一量章动，再把完整动力学的结果和慢进动估算对比。

为什么值得理解仪器里的飞轮和边转边晃的陀螺，遵循同一条刚体定律。只盯着轮子转，会错过真正要紧的运动：转轴本身在动。

先抓重点

核心问题: 为什么旋转飞轮会进动，而不只是向下倒？
一句话直觉: 重力产生侧向力矩，改变空间中的角动量；能量、角动量的竖直分量和它沿转轴的分量仍然守恒。

核心数学模型

对支点的转动惯量

I_3=\frac12mR^2,\qquad I_\perp=\frac14mR^2+m\ell^2

飞轮是均匀薄圆盘，装在无质量的轴上；距离从固定支点量到盘心。画面里的盘厚、支架和轮毂只是示意，不计质量。

矢量形式的运动方程

\dot{\mathbf L}=mg\ell\,\mathbf e_z\times\mathbf n,\qquad\dot{\mathbf n}=\frac{\mathbf L\times\mathbf n}{I_\perp}

所有量都对固定支点计算。对称轴和空间中的角动量一起变化，自转角用四元数记录。

角速度与轴向自转

\boldsymbol\omega=\frac{\mathbf L}{I_\perp}+\left(\frac1{I_3}-\frac1{I_\perp}\right)(\mathbf L\cdot\mathbf n)\mathbf n,\qquad\omega_3=\dot\psi+\dot\phi\cos\theta

沿轴的实际自转角速度，不等于欧拉自转角的变化率。两者一混，精确稳态进动的条件就会算错。

两条稳态进动分支

I_\perp\cos\theta\,\Omega^2-I_3\omega_3\Omega+mg\ell=0

只有初始章动为零时，实数根才给出恒定倾角。加一个扰动，这个解就被打破。转轴水平时，方程变成一次方程。

慢进动近似

\Omega_{\mathrm{slow}}\simeq\frac{mg\ell}{I_3\omega_3}

这个估算略去了横向转动惯量那一项。进动比自转慢得多时它才准，自转为零时无定义。

三个守恒量

E=\frac{I_\perp}{2}(\dot\theta^2+\dot\phi^2\sin^2\theta)+\frac{I_3\omega_3^2}{2}+mg\ell\cos\theta,\qquad L_z=\mathrm{const},\quad\mathbf L\cdot\mathbf n=\mathrm{const}

盘心与支点等高时势能为零。最大误差除以系统本身的非零能量和角动量尺度，所以哪怕某个守恒量本身为零，误差仍有意义。

直立稳定阈值

|I_3\omega_3|>2\sqrt{I_\perp mg\ell}

这是竖直附近小倾角的线性稳定条件。低于阈值时，严格竖直的旋转仍是精确解；给一点倾角，才看得出它不稳定。

支点对轴的作用力

\mathbf R=m\ell\ddot{\mathbf n}+mg\mathbf e_z

理想轴承既能向上推，也能向下拉。如果尖端只是自由搁在桌面上，就得另外建立单边接触和碰撞模型。

常见难点

角动量总是沿着转轴

容易误解不管转多快，角动量都和转轴同向。

正确理解只要转轴在动，横向转动惯量也会贡献角动量。在低自转或快分支下，比较蓝色箭头和转轴。

竖着不倒就是稳定

容易误解陀螺竖着转、一直不倒，就说明它稳定。

正确理解低于阈值时，不受扰动的竖直解依然存在。比较从小倾角开始的运动，看它对扰动怎么反应。

支点和桌面上的尖端一样

容易误解普通桌面也能把尖端往下拉。

正确理解这里的支点由轴承固定。桌面拉不住尖端；滑动、脱离和碰撞都不在模型里。

动手实验

01
设置并转动视角
选择普通释放，依次看全景、俯视和侧视，再拖动橙色倾角手柄。
你会看到： 转轴的初始方向一改，装置画面和记录的数值都跟着变。
02
比较稳态分支
分别选慢分支和快分支，章动扰动保持为零，再反转自转。
你会看到： 倾角保持不变，是因为初态正好设在这个解上；释放时进动速度不对，一般就会章动。
03
检验近似和直立陀螺
降低自转，看近似还准不准；圆锥陀螺可以比较稳定侧和失稳侧的小倾角预设。
你会看到： 慢进动近似可能失效；直立稳定阈值要靠扰动才看得出来。
04
检查数值误差
回看完整记录，看守恒误差和两种步长的差别；用精细相位单步检查转子标记。
你会看到： 数值误差小，说明模型算得准；但模型略去的摩擦和桌面接触，它验证不了。