📚 课程目录第4章二阶线性方程与振动模型

4. 二阶线性方程与振动模型

许多工程振动问题最终都归结为同一个数学对象：常系数二阶线性微分方程。弹簧-质量系统和RLC串联电路虽然物理本质完全不同，却共享完全相同的方程结构与解的行为。理解这种"同构"关系，是本节的核心目标。

先约定参照系：取质量块静止平衡位置为原点，向下为正方向，位移记为 x(t)。设质量 m、弹簧刚度 k、阻尼系数 c，外加激励力 F(t)。对质量块受力分析：弹簧恢复力为 −kx，阻尼力为 −c·ẋ，外力为 F(t)。由牛顿第二定律得

m·ẍ + c·ẋ + k·x = F(t)

两边除以 m，写成标准形式 ẍ + 2ζω_n·ẋ + ω_n²·x = f(t)，其中无阻尼固有频率 ω_n = √(k/m)，阻尼比 ζ = c/(2√(km))。这两个参数完全决定了系统的动态特性。

RLC串联电路以电容电压 u_C 为未知量，由基尔霍夫电压定律得 L·u̇_C″ 形式：LC·ü_C + RC·u̇_C + u_C = u_s(t)。对照可见：L↔m，R↔c，1/C↔k，u_s↔F。这种一一对应意味着：弹簧系统的任何结论都可以直接翻译到电路上。

齐次方程 m·ẍ + c·ẋ + kx = 0 对应特征方程 mr² + cr + k = 0，其根为

r₁,₂ = [−c ± √(c² − 4mk)] / (2m)

判别式 Δ = c² − 4mk 的符号决定三种情形，这正是阻尼比 ζ 与 1 比较的几何含义：

过阻尼（ζ > 1，Δ > 0）：两个不等实根，通解 x = C₁e^r₁t + C₂e^r₂t，无振荡，缓慢回到平衡。
临界阻尼（ζ = 1，Δ = 0）：重根 r = −ω_n，通解 x = (C₁ + C₂t)e^−ω_nt，最快无振荡回归。
欠阻尼（0 < ζ < 1，Δ < 0）：共轭复根 −ζω_n ± iω_d，通解 x = e^−ζω_nt(C₁cosω_dt + C₂sinω_dt)，其中阻尼振荡频率 ω_d = ω_n√(1−ζ²)。

非齐次方程的通解 = 齐次通解 + 一个特解。特解由激励形式决定：若 F(t) = F₀cos(Ωt)，则设特解 x_p = Acos(Ωt) + Bsin(Ωt)，代入比较系数即可求出 A、B。齐次部分随 t 增大而衰减（瞬态），特解持续存在（稳态），这是受迫振动的核心结构。

以一组教学示例参数说明：m = 1 kg，k = 100 N/m，则 ω_n = 10 rad/s。改变阻尼 c 观察行为变化：

c (N·s/m)	ζ	特征根	行为
0	0	±10i	等幅振荡，周期 0.628 s
4	0.2	−2 ± 9.80i	欠阻尼，衰减振荡
20	1.0	−10（重根）	临界阻尼，最快无振荡回归
40	2.0	−2.68, −37.32	过阻尼，缓慢回归

受迫情形下，当激励频率 Ω 接近 ω_n 时振幅急剧增大，即共振。稳态振幅为 F₀/√((k−mΩ²)² + (cΩ)²)。ζ 越小，共振峰越高越尖；ζ 增大则峰被压低并左移。工程上常通过增大阻尼或使工作频率远离 ω_n 来避开共振区。

建模检查清单：① 明确正方向与平衡位置；② 写出各力（或各电压）的符号表达式；③ 化为标准形式并算出 ω_n、ζ；④ 由判别式判断阻尼类型；⑤ 用初始条件 x(0)、ẋ(0) 定常数；⑥ 受迫时区分瞬态与稳态。

数值求解时，可将二阶方程降为一阶方程组：令 y₁ = x，y₂ = ẋ，则 ẏ₁ = y₂，ẏ₂ = (F(t) − c·y₂ − k·y₁)/m。这样即可交给任意一阶 ODE 求解器处理，为后续数值方法章节做好准备。