2. 振动与声波传播的物理基础

数值模拟的对象是物理场，而物理场的行为由控制方程和定解条件共同决定。在动手写差分格式或调用求解器之前，必须先把“系统如何振动、波如何传播、边界与初始状态如何约束解”这几件事在物理层面理清楚。本节围绕单自由度与多自由度振动、一维波动方程、声压与质点速度的关系，以及边界条件和初始条件的物理含义展开，为后续离散化与编程实现建立参照系。

一、单自由度振动：一切振动问题的起点

最简单的振动系统是质量—弹簧—阻尼器构成的单自由度系统。取质量块偏离平衡位置的位移 u(t) 为未知量，规定向右为正，则其运动方程为：

m·ü + c·u̇ + k·u = F(t)

其中 m 为质量，c 为阻尼系数，k 为刚度，F(t) 为外激励力。三个参数组合出两个关键量：无阻尼固有角频率 ω_n = √(k/m)，以及阻尼比 ζ = c/(2√(km))。阻尼比决定了自由振动的形态：ζ < 1 为欠阻尼，响应是幅值按 e^−ζω_nt 衰减的振荡；ζ = 1 为临界阻尼，最快回到平衡而不振荡；ζ > 1 为过阻尼，缓慢单调衰减。

举个可手算的示例：取 m = 1 kg，k = 400 N/m，c = 4 N·s/m，则 ω_n = 20 rad/s，ζ = 4/(2×20) = 0.1，属于弱阻尼。若初始位移 u(0) = 0.01 m、初速度为零，自由振动周期约 2π/ω_d ≈ 0.315 s（ω_d = ω_n√(1−ζ²) ≈ 19.9 rad/s），每个周期幅值衰减到前一周期的 e^{−2πζ/√(1−ζ²)} ≈ 0.53 倍。这个“周期 + 衰减率”正是数值积分结果最直接的校验依据：若你的时程曲线周期明显偏离 0.315 s，或衰减速度与 0.53 不符，说明时间步长、阻尼处理或积分格式存在问题。

当外激励为简谐力 F = F₀sin(ωt) 时，稳态振幅由动力放大系数决定，在 ω ≈ ω_n 附近出现共振峰，峰值高度受 ζ 控制。这提醒我们：数值模拟中若激励频率接近系统固有频率，必须保证时间步长足够小以分辨共振响应，否则会严重低估或高估振幅。

二、多自由度振动：模态是理解复杂系统的钥匙

实际结构很少是单自由度。以两个质量块串联弹簧为例，其方程写成矩阵形式为 Mü + Cu̇ + Ku = F。无阻尼自由振动对应广义特征值问题 Kφ = ω²Mφ，解出各阶固有频率和振型向量 φ。振型的物理意义是：系统按某一阶模态振动时，各自由度位移保持固定比例同步运动，形成“驻波式”的整体形态。

多自由度系统的响应可以分解为各阶模态响应的叠加。这一性质对数值模拟极为重要：它意味着网格或自由度增加后，高频模态的周期越来越短，若时间积分格式是显式的，稳定性条件（如 CFL 条件）将由最高阶频率决定，步长必须相应缩小。因此建模时应有意识地判断“我关心的频段到哪里”，避免为不必要的高频模态付出计算代价。

三、波动方程：从离散振子到连续介质

把连续介质看作无限密集的质点链，单自由度方程就过渡为波动方程。以一维杆的纵向位移 u(x,t) 为例：

∂²u/∂t² = c²·∂²u/∂x²,   c = √(E/ρ)

其中 E 为弹性模量，ρ 为密度，c 为波速。例如钢材 E ≈ 2.1×10¹¹ Pa、ρ ≈ 7800 kg/m³，纵波波速约 5190 m/s。波速是数值模拟中最关键的参数之一：它直接决定显式格式的稳定步长 Δt ≤ Δx/c。若空间步长取 0.01 m，则时间步长需小于约 1.9×10⁻⁶ s，这个量级必须在写代码前就估算清楚。

波动方程的通解为 u = f(x−ct) + g(x+ct)，即右行波与左行波的叠加。这一图像帮助我们理解边界的作用：波传到边界时，若边界不能自由通过，就会产生反射波，入射波与反射波叠加形成驻波或复杂干涉图样。

四、声压与质点速度：声场的两个基本变量

在流体声学中，描述声场常用声压 p 与质点速度 v。二者通过动量方程耦合：

ρ₀·∂v/∂t = −∂p/∂x

对平面简谐波，二者满足 p = ρ₀c·v，比例系数 ρ₀c 称为特性阻抗。空气中 ρ₀ ≈ 1.2 kg/m³、c ≈ 343 m/s，特性阻抗约 412 Pa·s/m。这个关系在数值实现中非常实用：很多声学求解器只求解声压一个变量，质点速度通过上述关系或动量方程后处理得到；反之若采用速度—压力混合格式，则需保证两个变量在时间推进上同步耦合，否则会出现非物理的漂移。

需要注意声压是标量（对流体而言），而质点速度是矢量，方向与波传播方向一致。在设置激励时，若给定的是声压边界，实际注入的是压力扰动；若给定的是速度边界（如振速边界），则对应的是质点运动。两者不能混用，否则幅值会相差一个特性阻抗因子。

五、边界条件与初始条件：定解的两类约束

波动方程本身有无穷多解，必须由边界条件和初始条件共同确定唯一解。二者的角色不同，不能互相替代。

边界条件约束的是空间边界上的行为，常见三类：

• 狄利克雷条件（固定边界）：给定边界位移或声压，如 u(0,t) = 0，对应刚性壁面或固定端，波在此全反射且相位翻转。
• 诺伊曼条件（自由边界）：给定边界梯度为零，如 ∂u/∂x = 0，对应自由端或声学软边界，反射波不翻转相位。
• 阻抗/吸收边界：给定压力与速度的比例关系，模拟部分吸收或辐射到无限域，是避免虚假反射的关键手段。

初始条件约束的是 t = 0 时刻整个场的状态，对二阶方程需要两个：初始位移 u(x,0) 与初始速度 ∂u/∂t(x,0)。只给位移不给速度，或反之，方程无法唯一推进。常见的错误是把“初始静止”理解为只设位移为零而忘记速度也为零，或在脉冲激励问题中把激励错误地放进初始条件而非边界/源项。