常微分方程（十七）：物理与工程应用

微分方程不是纯数学游戏——它是理解物理世界的语言。 从天体运动到电路响应，从荡来荡去的单摆到桥缆背后的旋涡脱落，所有动力系统都在"讲" ODE。

本章是一次精心安排的巡礼。每个例子都会调用我们 1-16 章建立起来的全部工具：相平面、特征值、拉普拉斯变换、模态分析、守恒律、数值积分、控制。所有案例都不是玩具——它们是真实的工作物理，写得尽量紧凑，让结构清晰可见。

本章要点

非线性单摆：小角度线性化 vs 完整方程，周期与振幅，分界线（separatrix）
RLC 电路：三种阻尼情形，谐振，Q 因子，调谐
开普勒轨道：万有引力、偏心率、能量与角动量守恒、Bertrand 定理
多自由度结构振动：模态分析、谐振、调谐质量阻尼器（Tacoma 大桥的教训）
流体力学：Poiseuille 流、雷诺数、旋涡脱落、Stokes 沉降

前置知识

第 1-3 章 ── 一阶/二阶线性理论
第 6-7 章 ── 线性方程组与相图
第 8 章 ── 非线性稳定性与能量法
第 11 章 ── 数值方法（贯穿始终）

1. 非线性单摆 ── 非线性 ODE 的"Hello World"

刚性长度 $L$、重力场 $g$ 中的质点单摆满足

$$ \ddot\theta + \frac{g}{L}\sin\theta = 0. $$

小角度下用 $\sin\theta \approx \theta$ 替换，得到简谐振动，周期 $T_0 = 2\pi\sqrt{L/g}$，与振幅无关。这是伽利略当年利用的等时性奇迹。

完整方程并不等时：周期依赖振幅，精确表达式涉及第一类完全椭圆积分

$$ T(\theta_0) = \frac{4}{\omega_0}\,K\!\bigl(\sin(\theta_0/2)\bigr), \qquad \omega_0 = \sqrt{g/L}. $$

单摆：大角度 vs 线性化时间响应、周期-振幅曲线、全局相图（小幅摆动 / 分界线 / 翻转）。

左上：从同一个大初始角 $\theta_0 = 1.5$ rad 出发，线性方程和完整方程在前半个周期内一致，之后就开始漂——线性模型先丢相位，再丢幅值。右上：数值积分得到的 $T(\theta_0)$；0.5 rad 以下小角度等时几乎完美，但 $\theta_0 \to \pi$ 时周期发散，因为达到了不稳定平衡。下：完整 $(\theta, \dot\theta)$ 相图。分界线（红色）以内是小幅摆动（libration），以外是越过顶点的转动（rotation），$2\pi$ 周期重复。

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
from scipy.integrate import solve_ivp
import numpy as np

g, L = 9.81, 1.0

def pendulum(t, s):
    theta, omega = s
    return [omega, -(g/L)*np.sin(theta)]

# 不同初角速度：小幅摆动 vs 翻转
for v0 in [1.5, 4.0, 6.5, 7.5]:    # 后两者越顶
    sol = solve_ivp(pendulum, (0, 10), [0.0, v0],
                    t_eval=np.linspace(0, 10, 2000),
                    rtol=1e-10, atol=1e-12)
    # sol.y 中可见闭合环或 theta 单调

为什么重要。 单摆是最简单的非平凡保守非线性系统，相图已包含第 7-8 章所有概念：稳定 / 不稳定平衡、分界线、能量-周期关系、可积性。所有其他振子——分子振动、Josephson 结、等离子体波——都借用这个骨架。

2. RLC 电路 ── 同一个方程，换成铜线

串联 $R$-$L$-$C$ 回路在电压 $V(t)$ 驱动下，电容上电荷 $q$ 满足

$$ L\ddot q + R\dot q + \frac{1}{C}\,q = V(t). $$

对照机械的弹簧-质量-阻尼 $m\ddot x + c\dot x + kx = F(t)$。完全是同一个方程——所以阻尼比、固有频率的所有直觉都可以直接搬过来。

定义 $\omega_n = 1/\sqrt{LC}$（固有频率），$\zeta = (R/2)\sqrt{C/L}$（阻尼比）。从 $V$ 到 $q$ 的传递函数

$$ H(s) = \frac{1/L}{s^2 + (R/L)\,s + 1/(LC)}, $$

带来的还是第 16 章那三种阶跃响应：欠阻尼、临界阻尼、过阻尼。

左上：经典串联 RLC。右上：三种阻尼比下的阶跃响应——欠阻尼绕稳态值振荡，临界阻尼最快无超调，过阻尼缓慢。左下：$|H(j\omega)|$ 在 $\zeta$ 小时出现尖锐谐振峰。右下：$R$ 增大时 Q 因子 $Q = (1/R)\sqrt{L/C}$ 下降，3-dB 带宽 $\Delta\omega = R/L$ 变宽——选择性与速度的根本权衡。

无阻尼振子在共振频率驱动下幅值线性增长——1940 年的 Tacoma Narrows 大桥正是把这个写进了历史（详见第 4 节）。

3. 行星轨道 ── 牛顿遇见开普勒

牛顿万有引力在二维下

$$ \ddot{\mathbf r} \;=\; -\frac{GM\,\mathbf r}{|\mathbf r|^3}. $$

守恒量是总能量 $E = \tfrac12|\mathbf v|^2 - GM/|\mathbf r|$ 与角动量 $\mathbf L = \mathbf r \times \mathbf v$。它们 不需要 解析求解 ODE 即可推出开普勒三定律。

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
from scipy.integrate import solve_ivp
import numpy as np

GM = 1.0
def kepler(t, s):
    x, y, vx, vy = s
    r = np.hypot(x, y)
    return [vx, vy, -GM*x/r**3, -GM*y/r**3]

# 圆、椭圆、近抛物线
for v0 in [1.0, 0.85, 0.55]:
    sol = solve_ivp(kepler, (0, 15), [1, 0, 0, v0],
                    t_eval=np.linspace(0, 15, 4000),
                    method='DOP853', rtol=1e-10)

不同初速度下的开普勒轨道；$1/r^2$ vs $1/r^3$ 对比；能量与角动量的数值守恒。

左：从 $(1,0)$ 出发，切向速度 $v_y$ 不同——轨道在圆与狭长椭圆之间过渡，正是开普勒观察到的家族。右上：思想实验——把 $1/r^2$ 换成 $1/r^3$，轨道不再闭合，螺旋坠入中心（与 Bertrand 定理 一致：所有球对称势中只有 $1/r$ 与谐振势能让所有束缚轨道闭合）。右下：高阶数值积分在多个公转后仍把两个守恒量保持在 $10^{-7}$ 量级。

Bertrand 定理 是个安静而震撼的事实：所有球对称势中，仅 $\propto 1/r$ 和 $\propto r^2$ 两种能让所有束缚轨道闭合。这正是 为什么 我们的太阳系有稳定的闭合行星轨道存在。先有数学，后有天文。

4. 结构振动 ── 楼宇、桥梁、调谐质量阻尼器

一栋多层楼大致可看作由刚性弹簧（梁柱）连接的质量串。两层楼，质量 $m_1, m_2$，刚度 $k_1, k_2$：

$$ M\ddot{\mathbf x} + C\dot{\mathbf x} + K\mathbf x = \mathbf F(t), \qquad M = \begin{pmatrix} m_1 & 0 \\ 0 & m_2 \end{pmatrix}, \; K = \begin{pmatrix} k_1+k_2 & -k_2 \\ -k_2 & k_2 \end{pmatrix}. $$

模态分析。 特征值问题 $K\,\boldsymbol\phi = \omega^2 M\,\boldsymbol\phi$ 给出固有频率 $\omega_i$ 与模态形状 $\boldsymbol\phi_i$。把 $\mathbf x = \sum_i q_i(t)\,\boldsymbol\phi_i$ 代入，方程组解耦为 $n$ 个独立的单自由度振子——多自由度问题化为 $n$ 个标量问题。

左：两个模态形状——模态 1 同相（整楼一起摇晃），模态 2 反相（楼层相对运动）。中列：从 2 楼初始位移激发的时间响应（上）、含小阻尼的频率响应有两个干净的谐振峰（下）。右列：上——共振频率驱动且无阻尼时幅值线性增长，正是 Tacoma Narrows 的经典图景；下——加上一个针对该模态的小型 调谐质量阻尼器（TMD），单一高峰被拆成两个低峰，是摩天大楼防风抗震的标准技巧（台北 101 顶部 660 吨的金球）。

这个套路是普适的：无论系统是小提琴弦、机翼、涡轮叶片还是千禧桥，步骤都是 组装 $M, C, K$ → 解特征值 → 用模态解耦 → 在需要处加阻尼。

5. 流体力学 ── 从稳态管流到旋涡脱落

稳态管流（Poiseuille）

不可压粘性流在长直圆管中，稳态 Navier-Stokes 退化为径向 ODE

$$ \frac{1}{r}\frac{d}{dr}\!\Bigl(r\,\frac{du}{dr}\Bigr) = -\frac{1}{\mu}\frac{dp}{dz}, $$

边界条件 $u(R) = 0$（无滑移）、$u(0)$ 有限。解为抛物线轮廓 $u(r) = u_{\max}(1 - r^2/R^2)$，体积流率 $Q = \pi R^4 \Delta p / (8\mu L)$ 即 Hagen-Poiseuille 律。流率对半径 四次方 的依赖，解释了为什么动脉扩张能高效降低血压：半径加倍，同样压差下流量增 16 倍。

雷诺数与层流-湍流转捩

无量纲数 $\mathrm{Re} = \rho U L / \mu$ 衡量惯性 / 粘性效应之比。$\mathrm{Re} \lesssim 2300$ 时层流，线性 $Q \propto \Delta p$ 成立；越过该值，湍涡耗散额外能量，同样压差下流量更小。

旋涡脱落（Strouhal 数）

钝体绕流后，对称涡以频率 $f$ 交替脱落，使得 Strouhal 数 $\mathrm{St} = f D / U$ 在很宽的 Re 范围内基本恒定（约 0.21）。当此频率与结构模态匹配——又是 Tacoma 大桥的故事。

沉降小球（Stokes 阻力）

致密小球在粘性流体中下落

$$ m\dot v = (\rho_p - \rho_f) V g - 6\pi\mu R\,v, $$

是一阶线性 ODE，终端速度 $v_t = (\rho_p - \rho_f) V g / (6\pi\mu R)$，时间常数 $\tau = m/(6\pi\mu R)$。

流体力学组合：Poiseuille 抛物线、层流-湍流流率标度、Strouhal 数随 Re 变化、Stokes 沉降。

左上：Hagen-Poiseuille 抛物速度剖面。右上：$Q$-$\Delta p$ 关系展示层流的线性如何在湍流出现后转为 $\sqrt{\Delta p}$ 标度。左下：Strouhal 数在巨大 Re 范围上的曲线——亚临界区平台 $\approx 0.21$ 正是为何缆绳和烟囱以可预测频率振荡的原因。右下：三种半径的小球在重力下指数地达到终端速度，$v_t \propto R^2$。

6. 共同的语法

本章每个领域都说同一种五步语法：

写出牛顿 / 基尔霍夫 / Navier-Stokes 律 对一个无穷小元。
化为 ODE：靠对称性、量纲分析或模态投影。
找平衡点 并线性化；得到固有频率和阻尼比。
加上激励，找谐振；决定加阻尼或 改变激励频率。
数值积分 处理非线性主导的部分（大幅、湍流、多体）。

掌握这个循环，物理与工程的相当大一片江山就被串起来了。

领域	控制 ODE	关键无量纲数
力学	$\ddot\theta + (g/L)\sin\theta = 0$	–
电学	$L\ddot q + R\dot q + q/C = V(t)$	$Q = (1/R)\sqrt{L/C}$
轨道	$\ddot{\mathbf r} = -GM\mathbf r/	\mathbf r
结构	$M\ddot x + C\dot x + Kx = F$	模态阻尼 $\zeta_i$
流体（管）	$\nabla \cdot \boldsymbol\sigma = 0$ → 径向 ODE	Re
流体（尾迹）	– (经验)	St $\approx 0.21$

小结

微分方程是 物理学的通用语言。前 16 章你建好了语言体系；本章证明了几个经典模型——单摆、RLC、开普勒、多自由度楼宇、管流——能把你从本科物理直接带到工程实践。它们都没用到新的数学，只是在 识别结构 后挑了正确的工具。

下一章是终章，把视野推到经典菜单之外：Neural ODE、随机与分数阶方程，以及 ODE 与现代机器学习之间的连接。

参考资料

Kreyszig, Advanced Engineering Mathematics, Wiley (2011).
Taylor, Classical Mechanics, University Science Books (2005).
Goldstein, Poole & Safko, Classical Mechanics, Pearson (2002).
Den Hartog, Mechanical Vibrations, Dover (1985).
Acheson, Elementary Fluid Dynamics, Oxford (1990).

系列导航


上一章	第十六章：控制理论基础
当前	第十七章：物理与工程应用
下一章	第十八章：前沿专题与总结