学径XUEJING · 个人学习平台通知设置

课程库 Calculus 2 S6 · 二阶微分方程 › 第 6

弹簧振子(下)—— 过阻尼、临界阻尼与"越用力越慢"

普通课40 分钟

本课导读

上一课的阻尼还算温柔(c2<4mkc^2 < 4mk),物体好歹能摆几下。本课把阻尼拧到底:过阻尼(泡在蜂蜜里的弹簧,慢吞吞蠕回原位)和临界阻尼(恰好压住振荡的完美刻度,闭门器和仪表指针的秘方)。压轴是一个漂亮的反直觉事实:阻尼加得越猛,回到平衡位置反而越慢——工程师为什么锱铢必较地追求"临界",答案全在特征根的位置上。预计 40 分钟。

说明

前置工具箱:第 3 课情况一/二公式(两实根、重根);十字相乘解特征方程;指数衰减快慢比较(e2te^{-2t}et/2e^{-t/2} 死得快);上一课的判别式四档表。

🔧 30 秒回顾 · 本课要用的旧手艺(熟练可跳过)

① 情况一/二速查(第 3 课):两实根 → c1er1t+c2er2tc_1 e^{r_1 t} + c_2 e^{r_2 t};重根 → (c1+c2t)ert(c_1 + c_2 t)e^{rt}。微例:r2+52r+1=0r^2 + \tfrac{5}{2}r + 1 = 0 十字相乘 (r+12)(r+2)=0(r + \tfrac{1}{2})(r + 2) = 0——例题 1 原式。⚠️ 翻车点:重根写成两个独立指数项(影子解事故)。

② 指数衰减赛跑erte^{rt}r<0r < 0)中 r|r| 越大死得越快;tt \to \infty最慢的那项说了算。微例:43et/213e2t\tfrac{4}{3}e^{-t/2} - \tfrac{1}{3}e^{-2t} 的长期行为由 et/2e^{-t/2} 主宰。⚠️ 翻车点:误以为系数大的项主宰长期行为——指数面前系数是弟弟。

③ 求函数最大值:驻点法 x(t)=0x'(t) = 0。微例:x=ete2tx = e^{-t} - e^{-2t},令 x=et+2e2t=0x' = -e^{-t} + 2e^{-2t} = 0et=2e^t = 2t=ln2t = \ln 2——练习 3 原式。⚠️ 翻车点:解 et=2e2te^{-t} = 2e^{-2t} 时两边除以 e2te^{-2t} 的方向搞反。

过阻尼:蜂蜜里的弹簧

阻尼压倒弹性:c2>4mkc^2 > 4mk。特征方程 mr2+cr+k=0mr^2 + cr + k = 0 给出两个不同实根

r1,2=c±c24mk2mr_{1,2} = \frac{-c \pm \sqrt{c^2 - 4mk}}{2m}

关键观察:两个根都是负数。理由:c24mk<c2=c\sqrt{c^2 - 4mk} < \sqrt{c^2} = c,所以就连"加号根"r1=c+c24mk2mr_1 = \frac{-c + \sqrt{c^2 - 4mk}}{2m} 的分子也仍是负的。通解:

x=c1er1t+c2er2t,r2<r1<0x = c_1 e^{r_1 t} + c_2 e^{r_2 t}, \qquad r_2 < r_1 < 0

两个衰减指数的叠加,没有任何三角函数——不振荡,不过冲,只是无声无息地爬回平衡位置。这就是过阻尼(overdamped)。画面:把弹簧连物体整个泡进蜂蜜,拉开松手,它慢悠悠蠕回去,绝不多晃半下。

过阻尼最多穿越平衡位置一次。想看理由的点开:

🔍 选看 · 为什么最多过零一次

x(t)=c1er1t+c2er2t=0x(t) = c_1 e^{r_1 t} + c_2 e^{r_2 t} = 0c1,c20c_1, c_2 \neq 0,否则更没戏)。移项:

e(r1r2)t=c2c1e^{(r_1 - r_2)t} = -\frac{c_2}{c_1}

左边是严格单调的指数函数,任何水平线最多与它相交一次——所以方程最多一个解,即轨迹最多穿越 x=0x = 0 一次(初速度朝平衡位置猛冲时可能冲过头一次,然后从另一侧慢慢爬回)。同样的论证对临界阻尼的 (c1+c2t)ert=0(c_1 + c_2 t)e^{rt} = 0 也成立:线性函数 c1+c2tc_1 + c_2 t 最多一个零点。

例题 1 · 过阻尼全流程

系统 x+52x+x=0x'' + \tfrac{5}{2}x' + x = 0m=1,c=52,k=1m = 1, c = \tfrac{5}{2}, k = 1),x(0)=1x(0) = 1x(0)=0x'(0) = 0:从静止释放。求运动规律并分析长期行为。

第 1 步 · 判型与通解

Δ=2544=94>0\Delta = \tfrac{25}{4} - 4 = \tfrac{9}{4} > 0:过阻尼。十字相乘(🔧 工具箱 ①):(r+12)(r+2)=0\left(r + \tfrac{1}{2}\right)(r + 2) = 0r=12,2r = -\tfrac{1}{2}, -2。通解:

x=c1et/2+c2e2tx = c_1 e^{-t/2} + c_2 e^{-2t}
第 2 步 · 初值与长期行为

x(0)=c1+c2=1x(0) = c_1 + c_2 = 1x=12c1et/22c2e2tx' = -\tfrac{1}{2}c_1 e^{-t/2} - 2c_2 e^{-2t}x(0)=12c12c2=0c1=4c2x'(0) = -\tfrac{1}{2}c_1 - 2c_2 = 0 \Rightarrow c_1 = -4c_2。联立:4c2+c2=1c2=13, c1=43-4c_2 + c_2 = 1 \Rightarrow c_2 = -\tfrac{1}{3},\ c_1 = \tfrac{4}{3}

x=43et/213e2tx = \frac{4}{3}e^{-t/2} - \frac{1}{3}e^{-2t}

长期行为(🔧 工具箱 ②):e2te^{-2t} 项很快烧完,剩 x43et/2x \approx \tfrac{4}{3}e^{-t/2}——慢根 r1=12r_1 = -\tfrac{1}{2} 主宰归途。这个"慢根说了算"是本课压轴反转的伏笔,记住它。

临界阻尼:恰到好处的刻度

cc 从过阻尼往回拧,拧到 c2=4mkc^2 = 4mk 的瞬间:判别式归零,重根

r=c2m=km=ω0r = -\frac{c}{2m} = -\sqrt{\frac{k}{m}} = -\omega_0

(代入 c=2mkc = 2\sqrt{mk} 化简所得——临界时衰减率恰好等于固有频率,一个优雅的巧合。)通解(🔧 工具箱 ①重根版):

x=(c1+c2t)eω0tx = (c_1 + c_2 t)\,e^{-\omega_0 t}

这就是临界阻尼(critical damping):仍不振荡(线性因子撑不起摆动,最多过零一次),但它是不振荡阵营里回位最快的——阻尼再小一丝就开始晃,再大一丝就开始拖。临界阻尼系数记作:

ccr=2mkc_{\text{cr}} = 2\sqrt{mk}

例题 2 · 临界阻尼全流程

系统 x+2x+x=0x'' + 2x' + x = 0m=1,k=1m = 1, k = 1c=2=ccrc = 2 = c_{\text{cr}}),同样 x(0)=1x(0) = 1x(0)=0x'(0) = 0。求运动规律,并与例题 1 比谁先到家。

解法展开

(r+1)2=0(r + 1)^2 = 0,重根 r=1r = -1,通解 x=(c1+c2t)etx = (c_1 + c_2 t)e^{-t}

x(0)=c1=1x(0) = c_1 = 1。求导:x=(c2c1c2t)etx' = (c_2 - c_1 - c_2 t)e^{-t}x(0)=c21=0c2=1x'(0) = c_2 - 1 = 0 \Rightarrow c_2 = 1

x=(1+t)etx = (1 + t)\,e^{-t}

比赛结果:临界系统按 ete^{-t} 收敛,例题 1 的过阻尼系统按 et/2e^{-t/2} 收敛——临界更快。数值感受:t=6t = 6 时临界的 (1+6)e60.017(1+6)e^{-6} \approx 0.017,过阻尼的 43e30.066\tfrac{4}{3}e^{-3} \approx 0.066,差近 4 倍且差距持续拉大。同样的 m,km, k、更小的 cc,回家反而更快——怎么回事?往下看。

反直觉时刻:阻尼越猛,回家越慢

直觉说:阻尼是刹车,刹车踩得越狠停得越快。**错在哪?**看特征根在数轴上的走位。固定 m=k=1m = k = 1(则 ccr=2c_{\text{cr}} = 2),把 cc 从 2 往上拧:

cc特征根主宰归途的慢根
22(临界)1-1(重根)1-1
52\tfrac{5}{2}12, 2-\tfrac{1}{2},\ -212-\tfrac{1}{2}
550.21, 4.79\approx -0.21,\ -4.790.21\approx -0.21
10100.10, 9.90\approx -0.10,\ -9.900.10\approx -0.10

cc 越大,两根越"劈叉":快根冲向 -\infty(那部分运动瞬间压灭),慢根却越贴近 0——衰减近乎停摆。物理画面:蜂蜜太稠,物体连"回去"这个动作本身都被拖住了。收敛速度由慢根决定(工具箱 ②:最慢的项说了算),所以整体回位速度随 cc 增大而变慢

临界阻尼正是两根合并、"慢根最不慢"的黄金点——这就是工程师的执念所在:

  • 闭门器:临界(或略过阻尼)——门快速关上且不哐当反弹
  • 汽车悬架:调在临界附近(实际略欠阻尼)——过坎后车身迅速稳住,颠一下但不晃三晃
  • 指针式仪表:指针直奔读数不来回摆,等它晃停黄花菜都凉了
  • 高层建筑阻尼器:目标同样是快速吸收晃动能量而不引入新的振荡

✍️ 练习

练习 1 · 判型三连

固定 m=1m = 1c=6c = 6,判断 kk 取以下值时系统的阻尼类型,能写通解骨架更好:(a) k=5k = 5 (b) k=9k = 9 (c) k=13k = 13

✅ 参考答案

判别式 Δ=364k\Delta = 36 - 4k

(a) Δ=16>0\Delta = 16 > 0过阻尼r2+6r+5=(r+1)(r+5)=0r^2 + 6r + 5 = (r+1)(r+5) = 0x=c1et+c2e5tx = c_1 e^{-t} + c_2 e^{-5t}。 (b) Δ=0\Delta = 0临界阻尼(r+3)2=0(r+3)^2 = 0x=(c1+c2t)e3tx = (c_1 + c_2 t)e^{-3t}。 (c) Δ=16<0\Delta = -16 < 0欠阻尼r=3±2ir = -3 \pm 2ix=e3t(c1cos2t+c2sin2t)x = e^{-3t}(c_1\cos 2t + c_2\sin 2t)

同一根弹簧,kk 从软到硬把系统从"蠕行"推进到"振荡"——判别式一眼定乾坤。

练习 2 · 给工程师算刻度

质量 m=2m = 2、弹性系数 k=50k = 50 的减震系统,求临界阻尼系数 ccrc_{\text{cr}},并说明 c=25c = 25 时系统处于哪种状态。

✅ 参考答案

ccr=2mk=2100=20c_{\text{cr}} = 2\sqrt{mk} = 2\sqrt{100} = 20c=25>20c = 25 > 20:过阻尼——减震"过头",回位偏慢;工程上多半会把阻尼往 20 附近调。

练习 3 · 冲过头的一次

系统 x+3x+2x=0x'' + 3x' + 2x = 0x(0)=0x(0) = 0x(0)=1x'(0) = 1(在平衡位置被向正方向踢一脚)。求运动规律、最大位移及其发生时刻。

✅ 参考答案

(r+1)(r+2)=0(r+1)(r+2) = 0,过阻尼,x=c1et+c2e2tx = c_1 e^{-t} + c_2 e^{-2t}。初值:c1+c2=0c_1 + c_2 = 0c12c2=1-c_1 - 2c_2 = 1c1=1,c2=1c_1 = 1, c_2 = -1

x=ete2tx = e^{-t} - e^{-2t}

驻点(🔧 工具箱 ③):x=et+2e2t=0et=2t=ln20.69x' = -e^{-t} + 2e^{-2t} = 0 \Rightarrow e^t = 2 \Rightarrow t = \ln 2 \approx 0.69。最大位移 x(ln2)=1214=14x(\ln 2) = \tfrac{1}{2} - \tfrac{1}{4} = \tfrac{1}{4}。此后单调爬回平衡,不再过零——"最多穿越一次"这里体现为:出发时在零点,之后再也不回穿。

📝 随堂测验

随堂测验0 / 5 题正确

1. 过阻尼系统($c^2 > 4mk$)的两个特征根?

2. 临界阻尼的判别条件与解的形状是?

3. “阻尼越大回位越慢”的数学原因是?

4. 过阻尼或临界阻尼的系统,穿越平衡位置的次数最多是?

5. 闭门器和仪表指针追求临界阻尼(而非重阻尼)的原因是?

本课小结

  • 过阻尼c2>4mkc^2 > 4mk):两个负实根,解 c1er1t+c2er2tc_1 e^{r_1 t} + c_2 e^{r_2 t}——纯衰减无振荡;长期行为由慢根主宰
  • 临界阻尼c=ccr=2mkc = c_{\text{cr}} = 2\sqrt{mk}):重根 ω0-\omega_0,解 (c1+c2t)eω0t(c_1 + c_2 t)e^{-\omega_0 t}——不振荡阵营的回位冠军
  • 两者都最多穿越平衡位置一次(单调函数与水平线至多一交)
  • 反直觉核心:cc 越大两根越劈叉,慢根越贴近 0 ⟹ 阻尼过猛反而回家更慢;工程师因此死磕临界刻度(闭门器、悬架、仪表、楼宇阻尼器)
  • 至此齐次弹簧四档命运集齐:无阻尼永动、欠阻尼减幅摆、临界最速归、过阻尼蠕行

下一课给弹簧接上外力——非齐次方程(第 4 课的武器)重返战场,主角是那个能晃塌大桥、也能荡高秋千的现象:共振

划选正文任意文字可高亮、批注或加入复习卡

讨论

载入中…