第42届复赛·理论第六题 精讲(零跳步版)

弹簧-双球转速仪 · 转动参考系 · 角动量的矢量性 · 45分 · 每次叉乘、每次积分、每次泰勒展开全部写全
本题不配图(请对照原卷图 6a)。全题只有一台机器:一根被弹簧"拽着"的斜杆,套在竖直转轴上。
三个物理看点依次登场:①离心力与弹簧的博弈定平衡;②角动量不沿转轴 → 轴承受力(动不平衡);③冲击后 \(L_z\) 守恒 → 滑冰式变速振动

〇、背景:瓦特的蒸汽机如何"自己踩油门"

0.1 飞球调速器——工业革命的自动控制

1788 年,瓦特给蒸汽机装上了一个小装置:离心调速器(flyball governor)。一根竖直轴由蒸汽机带动旋转,轴顶用铰链挂着两根杆,杆端各挂一个重球。机器转得越快,球被"甩"得越高;球通过一套杠杆连着进汽阀门——球升高,阀门关小,机器减速;球下落,阀门开大,机器加速。

这是人类历史上最早的负反馈自动控制系统之一:不需要人盯着,机器自己把转速稳定在设定值附近。后来麦克斯韦还专门写了论文《论调速器》(1868),开创了控制理论这门学科。

0.2 本题装置 = 弹簧版调速器

本题的转速仪是同一个思想的"竞赛加强版":

所以这道 45 分的大题,本质上是用竞赛语言把"调速器 + 动平衡 + 角动量守恒"三件事串成了一条龙。


一、读题:把这台机器拆成物理模型

1.1 逐个零件过一遍(每个"光滑"都要翻译成功课)

1.2 坐标系与两个关键位置矢量(后面一切计算的弹药库)

地面系 O-xyz:O 在轴心,z 竖直向上,x 水平(图 6a 中向右),y 垂直纸面向里。三者满足右手系:

$$ \boldsymbol e_x\times\boldsymbol e_y=\boldsymbol e_z,\qquad \boldsymbol e_y\times\boldsymbol e_z=\boldsymbol e_x,\qquad \boldsymbol e_z\times\boldsymbol e_x=\boldsymbol e_y . $$

本题最常用的时刻是图 6a 所示瞬间:杆正好转到 x-z 平面内,且 a 在 z 轴左侧(\(-x\) 一侧)。定义沿杆指向 a 的单位矢量 \(\hat{\boldsymbol n}\),则此刻

$$ \hat{\boldsymbol n}=-\sin\theta\,\boldsymbol e_x+\cos\theta\,\boldsymbol e_z, \qquad \boldsymbol r_a=l\hat{\boldsymbol n},\qquad \boldsymbol r_b=-l\hat{\boldsymbol n}. $$
图 6a 里 a 在左、b 在右,所以 \(\boldsymbol r_a\) 的 x 分量是 负的。这个正负号在第 1.2 问叉乘时会连环出现,定错一个,后面 \(\boldsymbol L\) 的 x 分量符号全反。官方解答⑦式的 \((\cos\theta_0\,\boldsymbol e_x+\sin\theta_0\,\boldsymbol e_z)\) 方向正是"a 在 \(-x\) 侧"的结果。

1.3 本题有"两个转动",务必从第一秒就分清

第(1)问:\(\theta\) 被锁死在 \(\theta_0\),只有整体自转,是定轴转动问题(放到转动参考系里就是静力学)。第(2)问:冲击把 \(\theta\) 解放出来,\(\theta\) 与自转 两个自由度耦合,靠 \(L_z\) 守恒 + 能量守恒两把锁管住。全题结构一目了然:

$$ \text{(1.1) 静平衡}\;\longrightarrow\;\text{(1.2) 转起来的 }\boldsymbol L\text{ 与轴承力}\;\longrightarrow\;\text{(2) 解锁 }\theta\text{ 后的耦合振动}. $$

1.4 符号速查表

符号含义
\(2h\)轴承 A、B 间距,O 为中点
\(2l,\,M\)杆长、杆质量(匀质,线密度 \(M/2l\))
\(m_1,m_2\)球 a、b 质量(第 2 问起 \(m_1=m_2=m\))
\(k\)盘簧扭转系数,回复力矩 \(=k\theta\)
\(\theta,\ \theta_0,\ \theta_M\)Oa 与 \(+z\) 夹角;平衡值;冲击后最大值
\(\omega,\ \omega_0,\ \omega_d\)整体自转角速度;平衡自转角速度;\(\theta\) 小幅振动角频率
\(\Omega_0\)冲击后瞬间的 \(\dot\theta\)
\(\delta_0=\theta_M-\theta_0\)振幅(小量,2.3 问只保留它的一阶)

二、(1.1) 平衡角速度:离心力矩与弹簧的博弈

工具 1:转动参考系与惯性离心力(为什么这题可以"静止"着做)

跟着系统以 \(\omega_0\) 绕 z 轴一起转的参考系里,杆、球、轴全部静止——动力学问题降级成静力学问题,代价是每个质量 \(m\) 要补上两个"假想力":

物理图像:离心力把 a、b 两球(以及杆的每一小段)往"远离 z 轴"的方向甩,效果都是把杆往 \(\theta\) 增大(更平)的方向压;弹簧和重力差把杆往回拉。三者在 \(\theta_0\) 处打平。

工具 2:力矩正负号约定——"以 \(\theta\) 增大为正"到底是哪个矢量方向

官方解答①式开头注明"以 \(\theta\) 增大方向为正"。我们要把这个约定翻译成矢量语言,否则每个力矩的正负都没法检验。在图 6a 瞬间,若 \(\theta\) 增大,\(\hat{\boldsymbol n}\) 怎么动?对

$$ \hat{\boldsymbol n}=-\sin\theta\,\boldsymbol e_x+\cos\theta\,\boldsymbol e_z $$

求导(\(\boldsymbol e_x,\boldsymbol e_z\) 不动):

$$ \frac{d\hat{\boldsymbol n}}{d\theta}=-\cos\theta\,\boldsymbol e_x-\sin\theta\,\boldsymbol e_z . $$

另一方面,若杆以角速度 \(\boldsymbol\Omega=\dot\theta\,\boldsymbol u\) 转动(\(\boldsymbol u\) 为待定方向),则 \(\dfrac{d\hat{\boldsymbol n}}{dt}=\boldsymbol\Omega\times\hat{\boldsymbol n}\),即

$$ \frac{d\hat{\boldsymbol n}}{d\theta}=\boldsymbol u\times\hat{\boldsymbol n}. $$

试 \(\boldsymbol u=-\boldsymbol e_y\):

$$ (-\boldsymbol e_y)\times(-\sin\theta\,\boldsymbol e_x+\cos\theta\,\boldsymbol e_z) =\sin\theta\,(\boldsymbol e_y\times\boldsymbol e_x)-\cos\theta\,(\boldsymbol e_y\times\boldsymbol e_z) =-\sin\theta\,\boldsymbol e_z-\cos\theta\,\boldsymbol e_x \;\checkmark $$

完全一致!所以"\(\theta\) 增大"的转动方向 \(=-\boldsymbol e_y\)。约定:力矩矢量 \(\boldsymbol M\) 在 \(-\boldsymbol e_y\) 上的投影(即 \(-M_y\))为正,称为"使 \(\theta\) 增大"的力矩。下面每算一个力矩,我们都先算出矢量,再读它的正负。

2.1 四个力矩,一个一个慢慢算

① 盘簧力矩

题设:回复力矩、大小 \(k\theta_0\)、"回复"= 把 \(\theta\) 往小推。所以在我们的约定下直接写

$$ M_{\text{簧}}=-k\theta_0 . $$

② 球 a 的重力矩

力:\(-m_1g\,\boldsymbol e_z\),作用点 \(\boldsymbol r_a=l(-\sin\theta_0\,\boldsymbol e_x+\cos\theta_0\,\boldsymbol e_z)\)。力矩 \(=\boldsymbol r_a\times(-m_1g\,\boldsymbol e_z)\),逐项分配:

$$ \boldsymbol r_a\times(-m_1g\,\boldsymbol e_z) =l(-\sin\theta_0)\,\boldsymbol e_x\times(-m_1g\,\boldsymbol e_z)+l\cos\theta_0\,\boldsymbol e_z\times(-m_1g\,\boldsymbol e_z). $$

第二项 \(\boldsymbol e_z\times\boldsymbol e_z=0\),没了。第一项里 \(\boldsymbol e_x\times\boldsymbol e_z=-(\boldsymbol e_z\times\boldsymbol e_x)=-\boldsymbol e_y\),所以

$$ \boldsymbol r_a\times(-m_1g\,\boldsymbol e_z)=l\sin\theta_0\cdot m_1g\cdot(-\boldsymbol e_y)=-m_1gl\sin\theta_0\,\boldsymbol e_y . $$

它在 \(-\boldsymbol e_y\) 方向,按约定为正:

$$ M_{\text{重}a}=+m_1gl\sin\theta_0 . $$

物理检验:a 在斜上方、偏在 z 轴一侧,重力把它竖直往下拉,杆只会倒得更平——\(\theta\) 增大。✓

③ 球 b 的重力矩

作用点 \(\boldsymbol r_b=-\boldsymbol r_a=l(\sin\theta_0\,\boldsymbol e_x-\cos\theta_0\,\boldsymbol e_z)\),力 \(-m_2g\,\boldsymbol e_z\):

$$ \boldsymbol r_b\times(-m_2g\,\boldsymbol e_z)=l\sin\theta_0(-m_2g)(\boldsymbol e_x\times\boldsymbol e_z)=l\sin\theta_0(-m_2g)(-\boldsymbol e_y)=+m_2gl\sin\theta_0\,\boldsymbol e_y . $$

它在 \(+\boldsymbol e_y\) 方向,按约定为负:

$$ M_{\text{重}b}=-m_2gl\sin\theta_0 . $$

物理检验:b 在斜下方,重力把它继续往下拉 = 把 a 往上顶 = \(\theta\) 减小。✓ 若 \(m_1=m_2\),两球重力矩恰好抵消——这正是第 2 问"重力出局"的根源。

④ 两球的惯性离心力矩

在转动系中,a 到 z 轴的距离 \(\rho_a=l\sin\theta_0\),离心力大小 \(m_1\omega_0^2l\sin\theta_0\),方向:水平背离 z 轴。a 在 \(-x\) 侧,所以方向为 \(-\boldsymbol e_x\):

$$ \boldsymbol F_a^{\text{离}}=-m_1\omega_0^2l\sin\theta_0\,\boldsymbol e_x . $$

力矩:

$$ \boldsymbol r_a\times\boldsymbol F_a^{\text{离}} =l(-\sin\theta_0\,\boldsymbol e_x+\cos\theta_0\,\boldsymbol e_z)\times(-m_1\omega_0^2l\sin\theta_0\,\boldsymbol e_x). $$

\(\boldsymbol e_x\times\boldsymbol e_x=0\),只剩第二项:

$$ =l\cos\theta_0\cdot(-m_1\omega_0^2l\sin\theta_0)\,(\boldsymbol e_z\times\boldsymbol e_x) =-m_1\omega_0^2l^2\sin\theta_0\cos\theta_0\,\boldsymbol e_y . $$

按约定(\(-\boldsymbol e_y\) 为正)取正:

$$ M_{\text{离}a}=+m_1\omega_0^2l^2\sin\theta_0\cos\theta_0 . $$

b 完全对称:\(\boldsymbol F_b^{\text{离}}=+m_2\omega_0^2l\sin\theta_0\,\boldsymbol e_x\)(b 在 \(+x\) 侧,向外即 \(+x\)),

$$ \boldsymbol r_b\times\boldsymbol F_b^{\text{离}} =(-l\cos\theta_0\,\boldsymbol e_z)\times(m_2\omega_0^2l\sin\theta_0\,\boldsymbol e_x) =-m_2\omega_0^2l^2\sin\theta_0\cos\theta_0\,\boldsymbol e_y , $$ $$ M_{\text{离}b}=+m_2\omega_0^2l^2\sin\theta_0\cos\theta_0 . $$

注意一个容易糊涂的点:两球的离心力矩同号。因为"球在哪一侧"和"离心力往哪边甩"同时反号,叉乘完负负得正。直觉:甩 a 向外,杆变平;甩 b 向外,杆也变平——都是 \(\theta\) 增大。

⑤ 杆的惯性离心力矩(本题唯一要动积分的地方)

杆是连续体,切成小段。杆线密度 \(\lambda=\dfrac{M}{2l}\)。在 a 侧取距 O 为 \(r\)(\(0\le r\le l\))、长 \(dr\) 的一小段:

$$ \text{质量 }dm=\frac{M}{2l}\,dr,\qquad \text{到 z 轴距离 }r\sin\theta_0, $$ $$ \text{离心力大小 }dF=\omega_0^2r\sin\theta_0\cdot\frac{M}{2l}\,dr,\quad\text{方向与 }\boldsymbol F_a^{\text{离}}\text{ 相同}. $$

与④同理,它产生的力矩(按约定取正)

$$ dM_{\text{离杆}}=\underbrace{\omega_0^2r\sin\theta_0\cdot\frac{M}{2l}\,dr}_{dF}\cdot\underbrace{r\cos\theta_0}_{\text{力臂:力的作用线到 O 的垂直距离}} =\frac{M}{2l}\,\omega_0^2\sin\theta_0\cos\theta_0\,r^2\,dr . $$

力臂为什么是 \(r\cos\theta_0\):力是水平的,位矢沿杆长 \(r\),二者夹角的垂直分量 \(=\) 位矢在竖直方向的投影 \(r\cos\theta_0\)。这和④中 \(l\cos\theta_0\) 是同一个几何。

b 侧镜像的一段给出相同结果(同样负负得正),所以全杆:

$$ M_{\text{离杆}}=2\int_0^l \frac{M}{2l}\,\omega_0^2\sin\theta_0\cos\theta_0\,r^2\,dr =\frac{M}{l}\,\omega_0^2\sin\theta_0\cos\theta_0\int_0^l r^2\,dr . $$

积分 \(\displaystyle\int_0^l r^2\,dr=\left[\frac{r^3}{3}\right]_0^l=\frac{l^3}{3}\),代回:

$$ M_{\text{离杆}}=\frac{M}{l}\,\omega_0^2\sin\theta_0\cos\theta_0\cdot\frac{l^3}{3} =\frac{M}{3}\,\omega_0^2l^2\sin\theta_0\cos\theta_0 . $$

系数 \(\dfrac{M}{3}l^2\) 不是巧合:\(\displaystyle 2\int_0^l \frac{M}{2l}r^2dr=\frac{1}{3}Ml^2=\frac{1}{12}M(2l)^2\),正是杆绕中心垂直轴的转动惯量。离心力的总力矩 = "把全杆质量按转动惯量等效"的结果。

⑥ 汇总:①②式

把五个力矩按约定相加并令其为零(转动系中杆静止、力矩平衡):

$$ -k\theta_0+m_1gl\sin\theta_0-m_2gl\sin\theta_0+M_{\text{离}}=0 \tag{①} $$ $$ M_{\text{离}}=\left(m_1+m_2+\frac{M}{3}\right)\omega_0^2l^2\sin\theta_0\cos\theta_0 \tag{②} $$
官方解答②式第一行把杆的积分印成了 "\(-2\int_0^l\cdots\)",那是印刷笔误——只有取 "\(+2\int\)" 才能得到它自己第二行的结果 \(\left(m_1+m_2+\frac{M}{3}\right)\)。物理上也必须是加:杆的两半各自被往外甩,力矩同号叠加。考场上若发现自己的中间式与答案差个符号,先检查这类印刷坑。

2.2 解出 \(\omega_0\)(③式)

把②代入①,把含 \(\omega_0\) 的项移到右边:

$$ k\theta_0-(m_1-m_2)gl\sin\theta_0=\left(m_1+m_2+\frac{M}{3}\right)\omega_0^2l^2\sin\theta_0\cos\theta_0 , $$

两边除以 \(\left(m_1+m_2+\frac{M}{3}\right)l^2\sin\theta_0\cos\theta_0\),开根号:

$$ \boxed{\;\omega_0=\sqrt{\frac{k\theta_0-(m_1-m_2)gl\sin\theta_0}{\left(m_1+m_2+\dfrac{M}{3}\right)l^2\sin\theta_0\cos\theta_0}}\;} \tag{③} $$

量纲体检:分子 \(k\theta_0\) 是力矩(\(\mathrm{N\cdot m=kg\,m^2/s^2}\)),\(mgl\) 也是力矩,同纲可减 ✓;分母是"质量×长度²" = 转动惯量量纲。力矩/转动惯量 = 角加速度量纲 \(\mathrm{s^{-2}}\),开根号得 \(\mathrm{s^{-1}}\) ✓。

2.3 存在条件讨论:题目为什么让你"讨论"

③式根号里是一个分式。"实际存在此 \(\omega_0\)" ⟺ 根号内有意义且 \(\omega_0\) 为非零实数 ⟺ 分子、分母同号(题设已排除 \(\theta_0=0,\pi/2,\pi\) 这些让分母为零的退化情形)。讨论就是按 \(m_1\) 与 \(m_2\) 的大小、\(\theta_0\) 所在的区间,把"同号"翻译成不等式。

情形(i):\(m_1>m_2\)。记分子 \(N(\theta_0)=k\theta_0-(m_1-m_2)gl\sin\theta_0\)。

情形(ii):\(m_1\le m_2\)。此时 \(N=k\theta_0+(m_2-m_1)gl\sin\theta_0>0\) 恒成立(两项皆非负,且 \(k\theta_0>0\)),所以只可能要求分母为正:

$$ 0<\theta_0<\frac{\pi}{2} \tag{⑥} $$

直觉:b 不比 a 轻时,重力矩与弹簧一起把杆往竖直拉,离心力必须能赢,平衡只可能出现在上半区。

(1.1) 收官:平衡角速度③ + 存在条件④⑤⑥,共 7 分(①式 2 分,②③④⑤⑥各 1 分)。核心手法一句话:换转动参考系 → 静力矩平衡 → 根号非负分区间讨论。

三、(1.2) 角动量 \(\boldsymbol L\) 与轴承力:本题的"矢量陷阱"

工具 3:叉乘速成(本题一切矢量计算的乘法表)

右手系三条基矢的乘法表(循环序为正,反序添负号):

$$ \boldsymbol e_x\times\boldsymbol e_y=\boldsymbol e_z,\quad \boldsymbol e_y\times\boldsymbol e_z=\boldsymbol e_x,\quad \boldsymbol e_z\times\boldsymbol e_x=\boldsymbol e_y; \qquad \boldsymbol b\times\boldsymbol a=-\boldsymbol a\times\boldsymbol b;\qquad \boldsymbol a\times\boldsymbol a=0 . $$

配合分配律 \(\boldsymbol a\times(\boldsymbol b+\boldsymbol c)=\boldsymbol a\times\boldsymbol b+\boldsymbol a\times\boldsymbol c\),任何"基矢组合 × 基矢组合"都能逐项算到底。本题所有叉乘我们都只用这张表,不用行列式。

3.1 球 a 的角动量(⑦式,每步写全)

某时刻杆在 x-z 平面内。整体以 \(\boldsymbol\omega_0=\omega_0\boldsymbol e_z\) 自转,a 的速度 = 角速度叉位矢:

$$ \boldsymbol v_a=\boldsymbol\omega_0\times\boldsymbol r_a =\omega_0\boldsymbol e_z\times l(-\sin\theta_0\,\boldsymbol e_x+\cos\theta_0\,\boldsymbol e_z) =-\omega_0l\sin\theta_0\,(\boldsymbol e_z\times\boldsymbol e_x)+0 =-\omega_0l\sin\theta_0\,\boldsymbol e_y . $$

检验:a 在 \(-x\) 侧、绕 z 轴逆时针(俯视)转,此刻应朝 \(-y\)(纸面外)运动 ✓。

角动量 \(\boldsymbol L_a=\boldsymbol r_a\times m_1\boldsymbol v_a\):

$$ \boldsymbol L_a=l(-\sin\theta_0\,\boldsymbol e_x+\cos\theta_0\,\boldsymbol e_z)\times m_1(-\omega_0l\sin\theta_0\,\boldsymbol e_y) $$ $$ =m_1\omega_0l^2\sin\theta_0\Big[\sin\theta_0\,(\boldsymbol e_x\times\boldsymbol e_y)-\cos\theta_0\,(\boldsymbol e_z\times\boldsymbol e_y)\Big] =m_1\omega_0l^2\sin\theta_0\Big[\sin\theta_0\,\boldsymbol e_z+\cos\theta_0\,\boldsymbol e_x\Big]. $$ $$ \boldsymbol L_a=\boldsymbol r_1\times m_1\boldsymbol v_1=m_1l^2\sin\theta_0\,\omega_0\,(\cos\theta_0\,\boldsymbol e_x+\sin\theta_0\,\boldsymbol e_z) \tag{⑦} $$

3.2 球 b 的角动量(⑧式):方向居然和 a 相同

\(\boldsymbol r_b=-\boldsymbol r_a\),所以 \(\boldsymbol v_b=\boldsymbol\omega_0\times\boldsymbol r_b=+\omega_0l\sin\theta_0\,\boldsymbol e_y=-\boldsymbol v_a\):

$$ \boldsymbol L_b=\boldsymbol r_b\times m_2\boldsymbol v_b=(-\boldsymbol r_a)\times m_2(-\omega_0l\sin\theta_0\,\boldsymbol e_y)=\frac{m_2}{m_1}\,\boldsymbol L_a\ \text{的结构}, $$ $$ \boldsymbol L_b=\boldsymbol r_2\times m_2\boldsymbol v_2=m_2l^2\sin\theta_0\,\omega_0\,(\cos\theta_0\,\boldsymbol e_x+\sin\theta_0\,\boldsymbol e_z) \tag{⑧} $$

两个负号(位置反、速度反)再次抵消,a、b 的角动量方向完全一致。记住这个方向

$$ \boldsymbol u_L\equiv\cos\theta_0\,\boldsymbol e_x+\sin\theta_0\,\boldsymbol e_z , $$

它在 x-z 平面内、与 z 轴夹角 \(\dfrac{\pi}{2}-\theta_0\)——不沿转轴。这是本小问全部物理的源头。

3.3 杆的角动量(⑨式,积分再来一遍)

与 2.1⑤ 同样的切法:a 侧 \(r\) 处小段 \(dm=\dfrac{M}{2l}dr\),位置 \(r\hat{\boldsymbol n}\),速度 \(=-\omega_0r\sin\theta_0\,\boldsymbol e_y\),它贡献的角动量与⑦同结构(只是把 \(l\to r\)):

$$ d\boldsymbol L=\omega_0r^2\sin\theta_0\,(\cos\theta_0\,\boldsymbol e_x+\sin\theta_0\,\boldsymbol e_z)\,dm . $$

b 侧镜像小段贡献相同,全杆:

$$ \boldsymbol L_M=2\int_0^l\frac{M}{2l}\,\omega_0r^2\sin\theta_0\,(\cos\theta_0\,\boldsymbol e_x+\sin\theta_0\,\boldsymbol e_z)\,dr =\frac{M}{l}\,\omega_0\sin\theta_0\,(\cdots)\int_0^l r^2dr, $$ $$ \boldsymbol L_M=\frac{M}{3}\,l^2\sin\theta_0\,\omega_0\,(\cos\theta_0\,\boldsymbol e_x+\sin\theta_0\,\boldsymbol e_z) \tag{⑨} $$

——又是那个 \(\dfrac{M}{3}l^2\)。同一个转动惯量,既管离心力矩(②式),又管角动量(⑨式),因为它们都是"质量×\(r^2\)"的矩。

3.4 系统总角动量(⑩式)与它的"第二种写法"

⑦⑧⑨方向全同,直接相加。记等效转动惯量因子:

$$ \boldsymbol L=\boldsymbol L_a+\boldsymbol L_b+\boldsymbol L_M =\left(m_1+m_2+\frac{M}{3}\right)l^2\sin\theta_0\,\omega_0\,(\cos\theta_0\,\boldsymbol e_x+\sin\theta_0\,\boldsymbol e_z) \tag{⑩} $$

把③式的 \(\omega_0\) 代入,把系数收进根号(注意 \(\left(m_1+m_2+\frac{M}{3}\right)l^2\sin\theta_0\) 整体平方后进根号,与根号内分母约掉一个因子):

$$ \boldsymbol L=\sqrt{\left(m_1+m_2+\frac{M}{3}\right)l^2\sin\theta_0\,\frac{k\theta_0-(m_1-m_2)gl\sin\theta_0}{\cos\theta_0}}\;(\cos\theta_0\,\boldsymbol e_x+\sin\theta_0\,\boldsymbol e_z). $$
评分标准原话:⑦⑧⑨⑩式各 1 分,"若只给出角动量的某个分量,⑦⑧⑨⑩式不得分"。即必须写出矢量形式(两个分量都在),只写 \(L_z\) 或只写大小都算白给。\(\boldsymbol L\) 有 x 分量这一点本身就是采分点。

3.5 震撼事实:\(\boldsymbol L\) 不沿转轴,而且它在"扫圈"

\(\boldsymbol L\) 的大小不变,方向却跟着杆一起以 \(\omega_0\) 绕 z 轴旋转(x 分量其实是"此刻杆所在方位"的标签,下一时刻杆转到别的方位,\(\boldsymbol L\) 的水平分量跟着转)。一个方向在变的角动量意味着 \(\dfrac{d\boldsymbol L}{dt}\neq0\) ——由角动量定理,系统必须持续受到外力矩。重力提供不了这么多(且方向不对),这份力矩只能由轴承 A、B 来出。这就是"斜杆转动必振动机"的根本原因,工程上叫动不平衡。

工具 4:常模旋转矢量的求导——\(\dfrac{d\boldsymbol L}{dt}=\boldsymbol\omega\times\boldsymbol L\)

一个长度不变、以角速度 \(\boldsymbol\omega\) 整体旋转的矢量 \(\boldsymbol L\),其矢端在空间画圆。设 \(\boldsymbol L\) 与 \(\boldsymbol\omega\) 夹角 \(\alpha\),则矢端圆周半径 \(=L\sin\alpha\),线速度大小 \(=\omega L\sin\alpha\),方向沿切线。而

$$ |\boldsymbol\omega\times\boldsymbol L|=\omega L\sin\alpha,\qquad \text{方向}=\text{切线方向}, $$

所以 \(\dfrac{d\boldsymbol L}{dt}=\boldsymbol\omega\times\boldsymbol L\)。记忆法:旋转矢量的"速度"= 角速度叉它自己(和 \(\boldsymbol v=\boldsymbol\omega\times\boldsymbol r\) 一模一样——位矢就是跟着刚体转的特殊矢量)。本题 \(\boldsymbol\omega=\omega_0\boldsymbol e_z\)。

由⑩,\(\boldsymbol L=L_x\boldsymbol e_x+L_z\boldsymbol e_z\),其中

$$ L_x=\left(m_1+m_2+\frac{M}{3}\right)l^2\omega_0\sin\theta_0\cos\theta_0 . $$

用工具 4:

$$ \frac{d\boldsymbol L}{dt}=\omega_0\boldsymbol e_z\times(L_x\boldsymbol e_x+L_z\boldsymbol e_z) =\omega_0L_x(\boldsymbol e_z\times\boldsymbol e_x)+0=\omega_0L_x\,\boldsymbol e_y , $$ $$ \frac{d\boldsymbol L}{dt}=\left(m_1+m_2+\frac{M}{3}\right)\omega_0^2l^2\sin\theta_0\cos\theta_0\,\boldsymbol e_y \tag{⑯} $$

(官方把⑯排在受力分析之后,我们这里按逻辑先算好它,等下直接用。)

3.6 质心运动定理:先拿到两个轴承力的"和"(⑪–⑮)

第一步:质心在哪。三部分的 z、x 坐标(沿 \(\hat{\boldsymbol n}\) 方向的带符号距离):a 在 \(+l\),b 在 \(-l\),杆质心在 O(0):

$$ r_C=\frac{m_1l+m_2(-l)+M\cdot 0}{m_1+m_2+M}=\frac{m_1-m_2}{m_1+m_2+M}\,l \tag{⑪} $$

\(r_C>0\) 表示质心偏 a 侧。它是"沿杆方向"的带符号距离,质心的实际位置矢量为 \(r_C\hat{\boldsymbol n}\)。

第二步:质心加速度。质心绕 z 轴做匀速圆周运动,半径 \(=r_C\sin\theta_0\)(位置矢量 \(r_C\hat{\boldsymbol n}\) 的水平投影),角速度 \(\omega_0\)。向心加速度指向 z 轴;质心在 \(-x\) 侧(设 \(r_C>0\)),指向轴即 \(+\boldsymbol e_x\):

$$ \boldsymbol a_C=\omega_0^2r_C\sin\theta_0\,\boldsymbol e_x \tag{⑫} $$

若 \(r_C<0\)(\(m_1

第三步:受力清单。系统在竖直方向:总重力向下,地面在轴底的支撑力向上,二者抵消(题设"系统重力与地面对转轴的支持力抵消")。水平方向只剩轴承力:

$$ \boldsymbol F_A=F_{Ax}\boldsymbol e_x+F_{Ay}\boldsymbol e_y,\qquad \boldsymbol F_B=F_{Bx}\boldsymbol e_x+F_{By}\boldsymbol e_y . $$

(光滑轴承的力过轴线,没有 z 分量;此刻杆在 x-z 平面,我们允许轴承力有 y 分量——待会儿方程会证明它是零。)

由质心运动定理(系统总外力 = 总质量 × 质心加速度):

$$ \boldsymbol F_A+\boldsymbol F_B=(m_1+m_2+M)\,\boldsymbol a_C \tag{⑬} $$

写分量。x 方向,用⑫⑪:

$$ F_{Ax}+F_{Bx}=(m_1+m_2+M)\,\omega_0^2r_C\sin\theta_0=(m_1-m_2)\,\omega_0^2l\sin\theta_0 \tag{⑭} $$

y 方向(\(\boldsymbol a_C\) 无 y 分量):

$$ F_{Ay}+F_{By}=0 \tag{⑮} $$

3.7 角动量定理:再拿到两个轴承力的"差"(⑰–⑳)

工具 5:质点系对定点的角动量定理

\(\displaystyle\sum\boldsymbol M_{\text{外}}^{(O)}=\frac{d\boldsymbol L}{dt}\),其中力矩与角动量必须对同一定点(这里取 O)。系统的外力:重力(a、b、杆)、地面对轴底的支撑、轴承力 \(\boldsymbol F_A,\boldsymbol F_B\)、水平轴对杆的约束力(过 O,力矩为零)、弹簧(内力,杆与轴之间,不出现在系统外力清单里)。

各外力对 O 的力矩(y 分量才是主角):

① 重力矩 \(\boldsymbol M_G\):3.3 节算过 a、b 各自的重力矩矢量:a 贡献 \(-m_1gl\sin\theta_0\,\boldsymbol e_y\),b 贡献 \(+m_2gl\sin\theta_0\,\boldsymbol e_y\),杆的重力作用线过 O 力矩为零:

$$ \boldsymbol M_G=-(m_1-m_2)gl\sin\theta_0\,\boldsymbol e_y . $$

② 轴承力的力矩:A 在 O 正上方 \(h\boldsymbol e_z\) 处,B 在 \(-h\boldsymbol e_z\) 处:

$$ \boldsymbol M_A=h\boldsymbol e_z\times\boldsymbol F_A=hF_{Ax}(\boldsymbol e_z\times\boldsymbol e_x)+hF_{Ay}(\boldsymbol e_z\times\boldsymbol e_y)=hF_{Ax}\boldsymbol e_y-hF_{Ay}\boldsymbol e_x , $$ $$ \boldsymbol M_B=-h\boldsymbol e_z\times\boldsymbol F_B=-hF_{Bx}\boldsymbol e_y+hF_{By}\boldsymbol e_x . $$

③ 地面支撑力沿 z 轴、作用点在 z 轴上,力矩为零。

角动量定理(矢量式):

$$ \boldsymbol M_G+\boldsymbol M_A+\boldsymbol M_B=\frac{d\boldsymbol L}{dt} \tag{⑰} $$

x 分量:左边 \(=0-hF_{Ay}+hF_{By}\),右边⑯没有 x 分量 \(=0\):

$$ -hF_{Ay}+hF_{By}=0 \tag{⑱} $$

y 分量:左边 \(=-(m_1-m_2)gl\sin\theta_0+hF_{Ax}-hF_{Bx}\),右边 \(=\) ⑯:

$$ -(m_1-m_2)gl\sin\theta_0+hF_{Ax}-hF_{Bx}=\left(m_1+m_2+\frac{M}{3}\right)\omega_0^2l^2\sin\theta_0\cos\theta_0 \tag{⑲} $$

由⑮⑱联立:\(F_{Ay}+F_{By}=0\) 且 \(-F_{Ay}+F_{By}=0\),两式相加得 \(2F_{By}=0\):

$$ F_{Ay}=F_{By}=0 \tag{⑳} $$

轴承力此刻没有 y 分量——合理:需要的向心力(x 方向)和力矩(y 方向力偶)都由 x 方向的轴承力就能凑出来。

3.8 解方程组(㉑㉒):把四个式子拧成两个数

未知量只剩 \(F_{Ax},F_{Bx}\)。把⑭⑲写成"和与差"的标准形:

$$ F_{Ax}+F_{Bx}=\underbrace{(m_1-m_2)\omega_0^2l\sin\theta_0}_{\equiv\,P}, $$ $$ F_{Ax}-F_{Bx}=\frac{1}{h}\underbrace{\left[(m_1-m_2)gl\sin\theta_0+\left(m_1+m_2+\frac{M}{3}\right)\omega_0^2l^2\sin\theta_0\cos\theta_0\right]}_{\equiv\,R}. $$

两式相加除以 2 得 \(F_{Ax}\),相减除以 2 得 \(F_{Bx}\):

$$ F_{Ax}=\frac{P}{2}+\frac{R}{2h}=\frac{1}{2}(m_1-m_2)\omega_0^2l\sin\theta_0+\frac{1}{2h}(m_1-m_2)gl\sin\theta_0+\frac{1}{2h}\left(m_1+m_2+\frac{M}{3}\right)\omega_0^2l^2\sin\theta_0\cos\theta_0 , $$ $$ F_{Bx}=\frac{P}{2}-\frac{R}{2h}=\frac{1}{2}(m_1-m_2)\omega_0^2l\sin\theta_0-\frac{1}{2h}(m_1-m_2)gl\sin\theta_0-\frac{1}{2h}\left(m_1+m_2+\frac{M}{3}\right)\omega_0^2l^2\sin\theta_0\cos\theta_0 . $$

化简的关键一步:注意 \(R\) 里的第二项 \(\left(m_1+m_2+\frac{M}{3}\right)\omega_0^2l^2\sin\theta_0\cos\theta_0\),正是①②式右边的 \(M_{\text{离}}\),由①它等于 \(k\theta_0-(m_1-m_2)gl\sin\theta_0\)。于是

$$ R=(m_1-m_2)gl\sin\theta_0+\Big[k\theta_0-(m_1-m_2)gl\sin\theta_0\Big]=k\theta_0 . $$

重力项恰好消掉!再把 \(P\) 里的 \(\omega_0^2\) 用③代入:

$$ \frac{P}{2}=\frac{1}{2}(m_1-m_2)l\sin\theta_0\cdot\frac{k\theta_0-(m_1-m_2)gl\sin\theta_0}{\left(m_1+m_2+\dfrac{M}{3}\right)l^2\sin\theta_0\cos\theta_0} =\frac{1}{2}(m_1-m_2)\,\frac{k\theta_0-(m_1-m_2)gl\sin\theta_0}{\left(m_1+m_2+\dfrac{M}{3}\right)l\cos\theta_0}. $$

最终(官方㉑㉒的第二行形式):

$$ F_{Ax}=\frac{1}{2}(m_1-m_2)\,\frac{k\theta_0-(m_1-m_2)gl\sin\theta_0}{\left(m_1+m_2+\dfrac{M}{3}\right)l\cos\theta_0}+\frac{k\theta_0}{2h} \tag{㉑} $$ $$ F_{Bx}=\frac{1}{2}(m_1-m_2)\,\frac{k\theta_0-(m_1-m_2)gl\sin\theta_0}{\left(m_1+m_2+\dfrac{M}{3}\right)l\cos\theta_0}-\frac{k\theta_0}{2h} \tag{㉒} $$
物理拆读(比公式本身更值得讲):轴承力 = 两件完全不同来源的事叠加。第一项(A、B 同号、各一半):偏心质心做圆周运动需要的向心力,两个轴承平摊——\(m_1=m_2\) 时这项为零。第二项 \(\pm\dfrac{k\theta_0}{2h}\)(A、B 反号):一对力偶,力偶矩 \(=2h\cdot\dfrac{k\theta_0}{2h}=k\theta_0\),不多不少正好等于盘簧的扭矩——弹簧"拧"着杆,反作用就拧在轴上,轴承必须用一对力偶顶住。

3.9 解法二(官方附的"秒杀版"):两分钟推出同一个结果

转轴自身绕 z 轴匀速自转,它自己的角动量沿 z 轴且恒定,故作用在转轴上的外力矩之和为零。轴受什么力矩?轴承 A、B 的力(对 O 的力矩 \(\boldsymbol M_A,\boldsymbol M_B\))+ 盘簧的反作用扭矩 \(-\boldsymbol M_k\)(杆被弹簧以 \(k\theta_0\) 拧,弹簧就被杆以 \(k\theta_0\) 反拧,弹簧装在轴上)。水平轴对杆的支持力过 O 点、地面支撑沿 z 轴,力矩都是零。
$$ \boldsymbol M_A+\boldsymbol M_B-\boldsymbol M_k=0,\qquad \boldsymbol M_k=k\theta_0\,\boldsymbol e_y . \tag{⑯'⑰'} $$

写分量(用 3.7② 的展开):

$$ -hF_{Ay}+hF_{By}=0 \tag{⑱'} $$ $$ hF_{Ax}-hF_{Bx}-k\theta_0=0 \tag{⑲'} $$

⑲' 一句话就说完了 3.7 节⑯⑰⑲三个式子的事:轴承力偶之差 = 弹簧扭矩。再联立⑭(质心定理还是省不掉)即得㉑'㉒',与㉑㉒完全相同。

考场上用哪个?解法一(角动量定理)是通法,任何"斜转子"问题都能套;解法二依赖"弹簧扭矩已知"这个巧合,快但特殊性强。建议:主用解法一,用解法二验算——两条路得到同一组数,基本可以确信没错。

(1.2) 收官:\(\boldsymbol L\)(⑩,矢量形式!)+ \(\boldsymbol F_A,\boldsymbol F_B\)(㉑㉒ + ⑳),共 17 分。核心认知一句话:斜着转的转子,角动量不沿轴,\(d\boldsymbol L/dt\neq0\),轴承被迫出力——动不平衡。

四、(2) 冲击之后:\(\theta\) 解锁,两个自由度跳起舞来

从现在起 \(m_1=m_2=m\)。系统原处于第(1)问的稳定状态(\(\theta_0,\omega_0\))。\(t=0\) 时刻球 a 受到一个冲击:方向在"杆-z 平面"内、垂直于杆、使 \(\theta\) 增大,作用时间极短。此后 \(\theta\) 最大摆到 \(\theta_M\)(\(\theta_0<\theta_M<\pi\))。

工具 6:冲击近似三句话 + 为什么守恒的是 \(L_z\)

冲击时间 \(\Delta t\to0\)、冲力 \(F\to\infty\)、冲量有限。三条推论:

冲击力的冲量作用于 a,方向在杆-z 平面内 \(\Rightarrow\) 它对 O 的冲量矩 \(\boldsymbol r_a\times\boldsymbol J\) 垂直于杆-z 平面(沿水平轴方向),z 分量为零;O 点的约束冲量过 O 点,力矩为零;轴承对轴的冲量过 z 轴,对 z 轴力矩为零。所以冲击前后 \(L_z\) 守恒。又因冲击期间 \(\theta\) 不变(\(\sin\theta\) 不变),由 \(L_z\propto\omega\sin^2\theta\) 立刻知道:自转角速度在冲击瞬间不变,冲击只把 \(\dot\theta\) 从 0 拉到 \(\Omega_0\)。冲击后瞬间状态:\((\theta,\omega,\dot\theta)=(\theta_0,\omega_0,\Omega_0)\)。

4.0 先备弹药:两个自由度下的动能(㉕㉖)

此后任一时刻,每个质量元的速度都是两部分的矢量和,且这两部分互相垂直(一个在杆-z 平面内摆动,一个水平画圈):

$$ v^2=\underbrace{(r\dot\theta)^2}_{\text{摆动:在杆-z平面内,速率}=r|\dot\theta|} +\underbrace{(\omega r\sin\theta)^2}_{\text{自转:水平圆圈,半径 }r\sin\theta} . $$

距 O 为 \(r\) 的杆元:摆动平面内速率 \(=r|\dot\theta|\)(绕水平轴的圆周运动);到 z 轴距离 \(r\sin\theta\),自转线速率 \(=\omega r\sin\theta\)。球就是 \(r=l\) 的特例。两个方向垂直(一个在竖直平面内沿"经线",一个水平沿"纬线"),所以平方直接相加。

系统动能 = 摆动部分 + 自转部分:

$$ E_k=\frac{1}{2}I\dot\theta^2+\frac{1}{2}\cdot 2m(\omega l\sin\theta)^2+2\int_0^l\frac{M\,dr}{2l}\,(\omega r\sin\theta)^2 \tag{㉕} $$

其中自转部分 = 两球(\(\frac12\cdot 2m\omega^2l^2\sin^2\theta\))+ 杆(\(\frac{M}{l}\omega^2\sin^2\theta\cdot\frac{l^3}{3}=\frac{M}{3}\omega^2l^2\sin^2\theta\) 的一半……仔细点:\(\frac12\cdot\frac{M}{3}\omega^2l^2\sin^2\theta\));摆动部分里的 \(I\) 是系统对 O 点水平摆动轴的转动惯量:

$$ I=(m_1+m_2)l^2+\frac{1}{12}M(2l)^2=\left(2m+\frac{M}{3}\right)l^2 \tag{㉖} $$

合并(球、杆的自转系数与 \(I\) 完全相同——都正比于 \(\int r^2dm\)):

$$ E_k=\frac{1}{2}\left(2m+\frac{M}{3}\right)l^2\dot\theta^2+\frac{1}{2}\left(2m+\frac{M}{3}\right)l^2\omega^2\sin^2\theta . $$

势能:\(m_1=m_2\) 时两球一升一降,重力势能之和恒为常数(取为零);杆的重力势能不变;只剩盘簧 \(\frac12k\theta^2\)。机械能:

$$ E=\frac{1}{2}\left(2m+\frac{M}{3}\right)l^2\dot\theta^2+\frac{1}{2}\left(2m+\frac{M}{3}\right)l^2\omega^2\sin^2\theta+\frac{1}{2}k\theta^2 \tag{㉗} $$

为书写省事,以下记 \(\displaystyle I=\left(2m+\frac{M}{3}\right)l^2\)(官方㉖式同款)。

4.1 (2.1) 求冲击后瞬间的 \(\dot\theta=\Omega_0\)

第一步:冲击后 \(L_z\) 永远守恒。冲击之后系统只受重力(力矩水平)、弹簧力矩(沿水平轴方向)、轴承力(过 z 轴)——对 z 轴的力矩全为零,所以任意时刻:

$$ L_z=I\sin^2\theta\,\omega=I\sin^2\theta_0\,\omega_0 \tag{㉓} $$

\(L_z\) 的算法:由⑩式(把 \(m_1+m_2\to2m\)),\(\boldsymbol L\) 的 z 分量 \(=I\sin\theta\,\omega\cdot\sin\theta=I\omega\sin^2\theta\)。更初等的看法:每个质点对 z 轴的角动量 \(=m\rho\cdot\omega\rho=m\omega\rho^2\),\(\rho=r\sin\theta\),求和即得。

解出 \(\omega\) 作为 \(\theta\) 的函数(这步把两个自由度"降"成一个):

$$ \omega=\frac{\sin^2\theta_0}{\sin^2\theta}\,\omega_0 \tag{㉔} $$

第二步:两个特殊状态之间用机械能守恒。冲击刚结束的瞬间:\(\theta=\theta_0,\ \dot\theta=\Omega_0,\ \omega=\omega_0\)。摆到最大角的瞬间:\(\theta=\theta_M,\ \dot\theta=0\),由㉔ \(\omega=\omega_M\equiv\dfrac{\sin^2\theta_0}{\sin^2\theta_M}\omega_0\)。㉗式两边相等:

$$ \frac{1}{2}I\Omega_0^2+\frac{1}{2}I\omega_0^2\sin^2\theta_0+\frac{1}{2}k\theta_0^2 =\frac{1}{2}I\omega_M^2\sin^2\theta_M+\frac{1}{2}k\theta_M^2 \tag{㉘} $$

第三步:解出 \(\Omega_0\)(代数一步不跳)。把含 \(\Omega_0\) 的项孤立:

$$ \frac{1}{2}I\Omega_0^2=\frac{1}{2}I\Big[\omega_M^2\sin^2\theta_M-\omega_0^2\sin^2\theta_0\Big]+\frac{1}{2}k(\theta_M^2-\theta_0^2). $$

括号第一项用 \(\omega_M\) 的定义化简:

$$ \omega_M^2\sin^2\theta_M=\frac{\sin^4\theta_0}{\sin^4\theta_M}\omega_0^2\cdot\sin^2\theta_M=\omega_0^2\,\frac{\sin^4\theta_0}{\sin^2\theta_M}, $$

于是

$$ \omega_M^2\sin^2\theta_M-\omega_0^2\sin^2\theta_0 =\omega_0^2\sin^2\theta_0\left(\frac{\sin^2\theta_0}{\sin^2\theta_M}-1\right). $$

再把③在 \(m_1=m_2\) 下的形式 \(\displaystyle\omega_0^2=\frac{k\theta_0}{I\sin\theta_0\cos\theta_0}\) 代入:

$$ \frac{1}{2}I\omega_0^2\sin^2\theta_0\left(\frac{\sin^2\theta_0}{\sin^2\theta_M}-1\right) =\frac{1}{2}I\cdot\frac{k\theta_0}{I\sin\theta_0\cos\theta_0}\cdot\sin^2\theta_0\left(\cdots\right) =\frac{1}{2}k\theta_0\tan\theta_0\left(\frac{\sin^2\theta_0}{\sin^2\theta_M}-1\right). $$

两边乘 \(\dfrac{2}{I}\),开根号:

$$ \boxed{\;\Omega_0=\sqrt{\frac{k}{\left(2m+\dfrac{M}{3}\right)l^2}}\; \sqrt{\theta_0\tan\theta_0\left(\frac{\sin^2\theta_0}{\sin^2\theta_M}-1\right)+\theta_M^2-\theta_0^2}\;} \tag{㉙} $$
㉙根号里第一项是负的(\(\theta_M>\theta_0\Rightarrow\sin^2\theta_0/\sin^2\theta_M<1\),当 \(\theta_M\) 不超过 \(\pi-\theta_0\) 时),别慌:它与第二项(正)相加之和恒正——4.4 节我们会把它展开到 \(\delta_0\) 二阶,亲眼看到它等于 \(\omega_d^2\delta_0^2>0\)。考场上见到"负+正"的根号,先别推翻自己,这是能量守恒的正常长相。

4.2 (2.2) 小幅振动频率:有效势能法(解法一)

㉔已经把 \(\omega\) 消成了 \(\theta\) 的函数——把它塞回能量㉗,问题就坍缩成一个自由度 \(\theta\) 在"有效势能"里的一维运动。这招和中心力场里 \(U_{\text{eff}}=U+\dfrac{L^2}{2mr^2}\) 是同一个模子:守恒的角动量把"被消掉的自由度"变成一项排斥势。

把㉔代入㉗的自转项:

$$ \frac{1}{2}I\omega^2\sin^2\theta=\frac{1}{2}I\,\frac{\sin^4\theta_0}{\sin^4\theta}\,\omega_0^2\cdot\sin^2\theta =\frac{1}{2}I\omega_0^2\,\frac{\sin^4\theta_0}{\sin^2\theta}, $$ $$ E=\frac{1}{2}I\dot\theta^2+\underbrace{\frac{1}{2}I\omega_0^2\,\frac{\sin^4\theta_0}{\sin^2\theta}+\frac{1}{2}k\theta^2}_{\equiv\,U_{\text{eff}}(\theta)} \tag{㉚} $$

第一步:验证 \(\theta_0\) 是平衡点(\(U'_{\text{eff}}(\theta_0)=0\))。对 \(\sin^{-2}\theta\) 求导:

$$ \frac{d}{d\theta}\,\frac{1}{\sin^2\theta}=-2\cdot\frac{1}{\sin^3\theta}\cdot\cos\theta=-\frac{2\cos\theta}{\sin^3\theta}, $$ $$ U'_{\text{eff}}(\theta)=\frac{1}{2}I\omega_0^2\sin^4\theta_0\cdot\left(-\frac{2\cos\theta}{\sin^3\theta}\right)+k\theta =-I\omega_0^2\sin^4\theta_0\,\frac{\cos\theta}{\sin^3\theta}+k\theta . $$

代入 \(\theta=\theta_0\):

$$ U'_{\text{eff}}(\theta_0)=-I\omega_0^2\sin\theta_0\cos\theta_0+k\theta_0=0\quad\checkmark $$

(正是 \(m_1=m_2\) 的①②式)。意料之中——\(\theta_0\) 本来就是平衡位置,有效势能在此处取极值。

第二步:算 \(U''_{\text{eff}}(\theta_0)\)(决定频率的"等效劲度系数")。需要对 \(\dfrac{\cos\theta}{\sin^3\theta}\) 再求一次导,用乘积法则 \(\big(u\cdot v\big)'=u'v+uv'\),取 \(u=\cos\theta,\ v=\sin^{-3}\theta\):

$$ u'=-\sin\theta,\qquad v'=-3\sin^{-4}\theta\cdot\cos\theta, $$ $$ \frac{d}{d\theta}\,\frac{\cos\theta}{\sin^3\theta} =-\sin\theta\cdot\frac{1}{\sin^3\theta}+\cos\theta\cdot\left(-\frac{3\cos\theta}{\sin^4\theta}\right) =-\frac{1}{\sin^2\theta}-\frac{3\cos^2\theta}{\sin^4\theta}. $$

于是

$$ U''_{\text{eff}}(\theta)=-I\omega_0^2\sin^4\theta_0\left(-\frac{1}{\sin^2\theta}-\frac{3\cos^2\theta}{\sin^4\theta}\right)+k =I\omega_0^2\sin^4\theta_0\left(\frac{1}{\sin^2\theta}+\frac{3\cos^2\theta}{\sin^4\theta}\right)+k . $$

代入 \(\theta=\theta_0\)(注意 \(\sin^4\theta_0\) 与分母约掉):

$$ U''_{\text{eff}}(\theta_0)=I\omega_0^2\Big(\sin^2\theta_0+3\cos^2\theta_0\Big)+k =I\omega_0^2\Big(1+2\cos^2\theta_0\Big)+k . $$

最后一步用了 \(\sin^2+3\cos^2=(\sin^2+\cos^2)+2\cos^2=1+2\cos^2\theta_0\)。

第三步:泰勒展开能量。令 \(\theta=\theta_0+\delta\)(\(\delta\ll\theta_0\)),在 \(\theta_0\) 处展开到 \(\delta\) 二阶:

$$ U_{\text{eff}}(\theta_0+\delta)=U_{\text{eff}}(\theta_0)+\underbrace{U'_{\text{eff}}(\theta_0)}_{=0}\delta+\frac{1}{2}U''_{\text{eff}}(\theta_0)\delta^2+\cdots $$ $$ E=\frac{1}{2}I\dot\delta^2+\frac{1}{2}\Big[I\omega_0^2(2\cos^2\theta_0+1)+k\Big]\delta^2+\text{常量} $$

再用 \(\omega_0^2=\dfrac{k\theta_0}{I\sin\theta_0\cos\theta_0}\) 把方括号第一项改写:

$$ I\omega_0^2(2\cos^2\theta_0+1)=\frac{k\theta_0(2\cos^2\theta_0+1)}{\sin\theta_0\cos\theta_0}, $$ $$ E=\frac{1}{2}\left(2m+\frac{M}{3}\right)l^2\dot\delta^2+\frac{1}{2}k\left(1+\theta_0\,\frac{2\cos^2\theta_0+1}{\sin\theta_0\cos\theta_0}\right)\delta^2+\text{常量} \tag{㉛} $$

对比简谐振子的标准能量 \(E=\frac12m_{\text{eff}}\dot\delta^2+\frac12k_{\text{eff}}\delta^2\),读出 \(\omega_d^2=k_{\text{eff}}/m_{\text{eff}}\):

$$ \boxed{\;\omega_d=\sqrt{\frac{k}{\left(2m+\dfrac{M}{3}\right)l^2}\left(1+\theta_0\,\frac{2\cos^2\theta_0+1}{\sin\theta_0\cos\theta_0}\right)}\;} \tag{㉜} $$

由⑥式,\(m_1=m_2\) 属于情形(ii),\(0<\theta_0<\dfrac{\pi}{2}\),故 \(\sin\theta_0\cos\theta_0>0\),括号恒正,振动确实存在。物理上:离心力提供的"离心势垒"让弹簧变得更硬——\(\omega_d\) 永远大于裸弹簧频率 \(\sqrt{k/I}\),而且 \(\theta_0\) 越接近 \(\pi/2\)(转得越快),修正项越大。

4.3 (2.2) 解法二:转动定理直接怼(官方附法)

在随系统以角速度 \(\omega\) 同步转动的参考系中,杆绕 O 点水平轴做定轴转动。对这根轴用转动定理:回复力矩(弹簧)+ 离心力矩(②式推广到任意 \(\theta,\omega\),注意 \(m_1=m_2\) 时重力矩抵消)= \(I\ddot\theta\):

$$ -k\theta+\left(2m+\frac{M}{3}\right)\omega^2l^2\sin\theta\cos\theta=I\ddot\theta \tag{㉝} $$

用㉔消 \(\omega\)(\(\omega^2=\omega_0^2\sin^4\theta_0/\sin^4\theta\)):

$$ -k\theta+I\omega_0^2\sin^4\theta_0\,\frac{\cos\theta}{\sin^3\theta}=I\ddot\theta \tag{㉞'} $$

令 \(\theta=\theta_0+\delta\),左边在 \(\theta_0\) 处泰勒展开到 \(\delta\) 一阶:常数项 \(=0\)(平衡条件),一次项系数用 4.2 节第二步已算好的导数:

$$ \Big[-k-I\omega_0^2(2\cos^2\theta_0+1)\Big]\,\delta =-k\left(1+\theta_0\,\frac{2\cos^2\theta_0+1}{\sin\theta_0\cos\theta_0}\right)\delta=I\ddot\delta \tag{㉟'} $$

——同一个 \(\omega_d\),殊途同归。解法一(能量)省思考、解法二(力矩)省书写,考场上任选,另一个用来验算。

4.4 (2.3) 转速随时间的起伏:\(\omega(t)\) 的一阶展开

第一步:写出 \(\theta(t)\)。小振动下 \(\theta\) 在 \(\theta_0\) 附近简谐振荡;\(t=0\) 时 \(\theta=\theta_0\) 且 \(\dot\theta=\Omega_0>0\),所以取正弦相位;振幅 = 最大偏离 \(=\theta_M-\theta_0=\delta_0\):

$$ \theta(t)=\theta_0+\delta_0\sin\omega_dt\qquad(\delta_0\ll\theta_0) \tag{㉝} $$

第二步:把㉝代入㉔,展开到 \(\delta_0\) 一阶。核心是把 \(\sin^2(\theta_0+\delta)\)(其中 \(\delta=\delta_0\sin\omega_dt\))展开。先展开 \(\sin(\theta_0+\delta)\):

$$ \sin(\theta_0+\delta)=\sin\theta_0\cos\delta+\cos\theta_0\sin\delta\approx\sin\theta_0+\delta\cos\theta_0 $$

(用了 \(\cos\delta\approx1,\ \sin\delta\approx\delta\),都是一阶近似)。平方,扔掉 \(\delta^2\) 项:

$$ \sin^2(\theta_0+\delta)\approx\sin^2\theta_0+2\delta\sin\theta_0\cos\theta_0 =\sin^2\theta_0\,\Big(1+2\delta\cot\theta_0\Big). $$

再用近似公式 \(\dfrac{1}{1+x}\approx1-x\)(\(x=2\delta\cot\theta_0\ll1\)):

$$ \omega(t)=\frac{\sin^2\theta_0}{\sin^2(\theta_0+\delta)}\,\omega_0 =\frac{\omega_0}{1+2\delta\cot\theta_0} \approx\omega_0\Big(1-2\delta\cot\theta_0\Big). $$

代回 \(\delta=\delta_0\sin\omega_dt\):

$$ \boxed{\;\omega(t)=\frac{\sin^2\theta_0}{\sin^2\big(\theta_0+\delta_0\sin\omega_dt\big)}\,\omega_0 \;\approx\;\omega_0\Big(1-2\delta_0\cot\theta_0\,\sin\omega_dt\Big)\;} \tag{㉞} $$

其中

$$ \omega_d=\sqrt{\frac{k}{\left(2m+\dfrac{M}{3}\right)l^2}\left(1+\theta_0\,\frac{2\cos^2\theta_0+1}{\sin\theta_0\cos\theta_0}\right)}, \qquad \omega_0=\sqrt{\frac{k\theta_0}{\left(2m+\dfrac{M}{3}\right)l^2\sin\theta_0\cos\theta_0}} . $$
物理图像(讲给学生的高光时刻):\(\theta\) 增大 → 球和杆离 z 轴更远(\(\sin\theta\) 变大)→ 转动惯量变大 → \(L_z\) 守恒逼着 \(\omega\) 变小;\(\theta\) 回来,\(\omega\) 又涨回去。转速与摆角反相振荡,频率与摆角振动同为 \(\omega_d\),相对幅度 \(2\delta_0\cot\theta_0\)。这就是花样滑冰收臂加速、放臂减速的机器版——监测 \(\omega(t)\) 的微小波动,就能反推杆的振动,这正是"转速仪"作为传感器的另一面。

4.5 彩蛋(自洽性检验):\(\Omega_0\) 必须等于 \(\omega_d\delta_0\)

一个精讲课该做的闭环:㉙给的 \(\Omega_0\) 与㉝的初始条件 \(\dot\theta(0)=\delta_0\omega_d\cos0=\omega_d\delta_0\) 必须一致。把㉙按 \(\delta_0\) 展开验证。记 \(\theta_M=\theta_0+\delta_0\),逐块展开到 \(\delta_0^2\):

块 1(要带二阶,因为一阶会相互抵消):

$$ \sin^2\theta_M=\sin^2\theta_0+2\delta_0\sin\theta_0\cos\theta_0+\delta_0^2(\cos^2\theta_0-\sin^2\theta_0)+\cdots =\sin^2\theta_0\Big[1+2\delta_0\cot\theta_0+\delta_0^2(\cot^2\theta_0-1)\Big], $$ $$ \frac{\sin^2\theta_0}{\sin^2\theta_M}\approx1-2\delta_0\cot\theta_0-\delta_0^2(\cot^2\theta_0-1)+(2\delta_0\cot\theta_0)^2 =1-2\delta_0\cot\theta_0+\delta_0^2(3\cot^2\theta_0+1). $$

用了 \(\dfrac{1}{1+u}\approx1-u+u^2\),\(u=2\delta_0\cot\theta_0+\delta_0^2(\cot^2\theta_0-1)\),\(u^2\) 只保留首项平方。

块 2:\(\theta_M^2-\theta_0^2=2\theta_0\delta_0+\delta_0^2\)。

合成㉙的根号内(第一项乘上 \(\theta_0\tan\theta_0\)):

$$ \theta_0\tan\theta_0\Big[-2\delta_0\cot\theta_0+\delta_0^2(3\cot^2\theta_0+1)\Big]+2\theta_0\delta_0+\delta_0^2 $$ $$ =\underbrace{\Big[-2\theta_0\delta_0+2\theta_0\delta_0\Big]}_{\text{一阶精确抵消!}} +\delta_0^2\Big[\theta_0(3\cot\theta_0+\tan\theta_0)+1\Big]. $$

一阶项恒等于零——这就是为什么㉙看起像"\(\sqrt{\text{负}+\text{正}}\)"却在 \(\delta_0\) 小时恒正。再把 \(3\cot\theta_0+\tan\theta_0=\dfrac{3\cos^2\theta_0+\sin^2\theta_0}{\sin\theta_0\cos\theta_0}=\dfrac{1+2\cos^2\theta_0}{\sin\theta_0\cos\theta_0}\) 代入:

$$ \Omega_0^2=\frac{k}{I}\,\delta_0^2\left(1+\theta_0\,\frac{2\cos^2\theta_0+1}{\sin\theta_0\cos\theta_0}\right)=\omega_d^2\delta_0^2 \quad\Longrightarrow\quad \Omega_0=\omega_d\delta_0\ \checkmark $$

三套独立推导(能量守恒㉙、有效势能㉜、简谐运动㉝)在同一个数上握手——做到这一步,可以确信全题无懈可击。


五、复盘

5.1 全题逻辑地图

(1.1) 转动参考系(杆静止)→ 五个力矩(弹簧/重力×2/离心×球/离心×杆积分)→ 力矩平衡①② → \(\omega_0\)③ → 根号同号讨论④⑤⑥。

(1.2) 逐质点叉乘 \(\boldsymbol r\times m\boldsymbol v\) → ⑦⑧⑨同方向叠加 → \(\boldsymbol L\)⑩(不沿 z 轴)→ 旋转矢量求导 \(d\boldsymbol L/dt=\boldsymbol\omega_0\times\boldsymbol L\)⑯ → 质心定理给"和"⑭⑮、角动量定理给"差"⑱⑲ → ㉑㉒(巧解:轴上力矩平衡⑲'一步到位)。

(2.1) 冲击三近似 → \(L_z\) 守恒㉓ → \(\omega(\theta)\)㉔ → 动能㉕㉖+机械能㉗ → 两特殊态能量守恒㉘ → \(\Omega_0\)㉙。

(2.2) ㉔代入㉗得有效势能㉚ → \(U'=0\) 验证、\(U''\) 求出 → 泰勒到二阶㉛ → \(\omega_d\)㉜(或转动定理㉝–㉟')。

(2.3) 简谐解㉝ → 代回㉔,\(\sin^2\) 展开 + \(1/(1+x)\) 展开 → \(\omega(t)\)㉞。闭环:\(\Omega_0=\omega_d\delta_0\)。

5.2 本题最值得记住的五个"手法"

5.3 官方解答的一处印刷瑕疵

②式第一行的杆积分前印成 "\(-2\int_0^l\)",按其第二行结果 \(\left(m_1+m_2+\frac{M}{3}\right)\) 及物理符号分析,应为 "\(+2\int_0^l\)"。引用官方解答时以结果为准。

5.4 易错清单(考前自查)

5.5 评分标准对照(45 分)

问分值明细
(1.1)7①式 2 分;②③④⑤⑥式各 1 分
(1.2)17⑦⑧⑨⑩式各 1 分(只给分量不得分);⑭式 2 分;⑮式 1 分;⑯式 2 分;⑱式 1 分;⑲式 2 分;⑳式 1 分;㉑㉒式各 2 分(解法二带撇式同分)
(2.1)9㉓式 2 分;㉔式 1 分;㉕式 2 分;㉖式 1 分;㉘式 2 分;㉙式 1 分
(2.2)8㉚㉛式各 3 分;㉜式 2 分(解法二同分)
(2.3)4㉝㉞式各 2 分

分数分布的启示:这 45 分里,"会物理"(建模、守恒律选择)大约决定 1/3,"算得对"(叉乘、积分、泰勒展开、解方程组)决定 2/3。第 (2) 问 21 分几乎全是"守恒律 + 展开"的基本功——这正是本精讲把每一步代数都写出来的原因。