1. 课程导入与数值模拟概览

振动与声波传播是工程中极为常见的一类物理过程：机床的颤振、汽车车身的低频轰鸣、管道中的气动噪声、扬声器振膜向空气辐射的声场，本质上都可以归结为“扰动在弹性介质或流体介质中随时间与空间传播”。要定量回答“振动有多大”“噪声在何处最强”“改变结构后声场如何变化”，仅靠解析公式往往力不从心，因为真实结构几何复杂、材料非均匀、边界条件多样。数值模拟正是为这类问题提供定量答案的核心手段。

一、从物理量到控制方程

振动与声波传播的数学描述建立在几个基本量之上：位移 u、速度 v、加速度 a、应力 σ、应变 ε，以及声学中的声压 p、质点振速 v 与介质密度 ρ。结构振动通常由弹性动力学方程描述，声波在流体中则由波动方程或亥姆霍兹方程描述。以简谐假设 p(x,t)=P(x)e^{jωt} 代入，时域波动方程转化为频域亥姆霍兹方程：

∇²P + k²P = 0,   k = ω/c = 2πf/c

其中 k 为波数，c 为声速，f 为频率。这个方程看似简单，但只有在极少数规则几何与边界下才有解析解。数值模拟的任务，就是把连续的偏微分方程离散为有限维代数方程组，再借助计算机求解。

二、数值模拟在振动声学中的角色

数值模拟并非替代理论，而是把理论“落地”为可计算的工程工具。它主要承担四类工作：其一，预测——在设计阶段估计结构的固有频率、振型、频响函数与辐射声压级；其二，诊断——对已有异常振动或噪声问题定位主要贡献源与传递路径；其三，优化——比较不同板厚、加强筋布置、吸声材料参数对声学指标的影响；其四，外推——在试验成本高或工况危险的场合，用模型探索极端条件。

需要强调的是，数值结果的可信度取决于三件事：物理模型是否抓住了主要机理、离散网格是否足以分辨最短波长、边界条件与材料参数是否接近真实。工程上常用“每波长至少 6～10 个单元”作为声学网格的经验准则。例如空气中 c≈343 m/s，在 f=1000 Hz 时波长约 0.343 m，若按每波长 8 个单元计，单元尺寸应控制在 0.04 m 量级；频率升高到 4000 Hz，单元尺寸需缩小到约 0.01 m，计算量随之显著增加。这一“频率—波长—网格”的制约关系，是贯穿本课程的重要主线。

三、典型工程场景

在汽车行业，车身板件在路面激励下的低频结构声是典型问题：先用有限元计算白车身模态，再把振型作为边界条件输入声学边界元模型，得到车内驾驶员耳旁的声压级频谱，进而判断哪个板件贡献最大。在电子设备中，风扇与硬盘引起的振动通过支架传到外壳并辐射噪声，工程师会建立“激励源—传递路径—辐射面”的链路模型，逐段评估隔振垫刚度与阻尼的影响。在建筑声学中，楼板撞击声的预测需要同时考虑结构振动与房间声场的耦合。这些场景的共同点是：单一物理场不足以描述问题，必须把结构域与声域连接起来。

四、整体学习路径与工具链

本课程按“基础—方法—实践—综合”推进。基础部分建立波动方程、边界条件与离散化概念；方法部分分别展开有限元、边界元、时域差分与统计能量法的适用条件与实现要点；实践部分围绕典型算例完成建模、求解、后处理与误差判断；综合部分处理流固耦合与多物理场问题。工具链上，几何建模与网格划分负责前处理，求解器完成离散方程求解，后处理把场量转化为频谱、云图与贡献量排序，脚本则用于参数扫描与批量计算。

五、一个可动手的入门示例

为建立直观感受，可先做一个一维杆纵向振动的小算例。取长度 L=1 m、弹性模量 E=210 GPa、密度 ρ=7850 kg/m³ 的钢杆，两端自由。其一阶纵向固有频率解析解为 f₁=c/(2L)，其中 c=√(E/ρ)≈5172 m/s，故 f₁≈2586 Hz。若用有限元把杆划分为 N 个线性单元，单元长度 h=L/N，则一阶频率的离散结果会随 N 增大而单调逼近解析值。建议读者取 N=5、10、20、40 分别计算，记录频率误差，观察“网格加密—误差下降”的趋势。这一过程虽简单，却完整体现了数值模拟的核心逻辑：用离散模型逼近连续物理，并用收敛性判断结果是否可信。示例中的材料与几何参数为教学假设，实际工程应以真实数据替换。

当读者能够独立完成上述算例，并解释误差来源时，就已经具备了进入后续各章方法学习的基础。