第1章：声学基础回顾——波动方程、声压与声强、声阻抗、声波在介质中的传播特性
各位同学，大家好。我是你们这门课的主讲。说实话，每次开新课，我都要纠结一下：基础讲多深？讲太浅，后面做超材料优化时你们会懵；讲太深，又怕把你们吓跑。
我的经验是：基础不牢，地动山摇。尤其是声学超材料，说白了就是在跟声波“斗智斗勇”。你不了解声波怎么走、怎么反射、怎么折射，你怎么去设计结构控制它？
所以，今天这一章，咱们把声学最核心的几个概念过一遍。我尽量用我这些年踩过的坑，帮你们把抽象的东西具象化。
本章核心目标：理解声波的本质，掌握波动方程、声压、声强、声阻抗这四个基石。它们是后续所有仿真和优化的“语言”。
1.1 波动方程——声学的“牛顿定律”
我个人习惯，把波动方程看作是声学世界的牛顿第二定律。你想想看，F=ma 描述了力怎么让物体动起来；波动方程则描述了压力扰动怎么在介质中传播。
标准的三维波动方程长这样：
∇²p = (1/c²) · ∂²p/∂t²其中 ∇² 是拉普拉斯算子，p 是声压，c 是声速，t 是时间。
这个方程怎么来的？其实是从三个物理定律推导出来的：
- 运动方程（牛顿第二定律）：压力梯度让质点加速
- 连续性方程（质量守恒）：流入的质量等于密度的变化
- 状态方程（热力学关系）：压力和密度通过声速关联
我在项目中遇到过一件事：有次做消声器仿真，怎么算结果都不对。折腾了两天，最后发现是介质假设错了——我把空气当成了不可压缩流体。嗯，波动方程的前提就是介质可压缩，这个坑你们别踩。
我的小技巧：做仿真时，先用手算一个简单的一维情况验证。比如平面波在管道中传播，波动方程退化成 d²p/dx² = (1/c²)d²p/dt²，解起来很直观。这一步能帮你快速发现模型设置有没有问题。
1.2 声压与声强——耳朵听到的和仪器测到的
声压，说白了就是声音引起的压力变化。大气压是 101325 Pa，你说话时产生的声压可能只有 0.01 Pa。但就是这点微小的变化，我们的耳朵能感知到。
声压用 p 表示，单位是帕斯卡（Pa）。有效声压 p_rms 是均方根值：
p_rms = p_peak / √2 （对于正弦波）声强呢？它表示单位面积上通过的声功率。我经常跟新人说：声压是“多响”，声强是“多有力”。
I = p² / (ρc)其中 ρ 是介质密度，c 是声速。对于空气，ρc ≈ 415 Rayl（瑞利）。
这里有个容易混淆的地方：声压是标量，声强是矢量。声强有方向！
注意：在近场测量时，声压和声强的关系会变得复杂。我曾经在测量一个微型扬声器时，用声压推算声强，结果差了 3dB。后来才发现，近场存在“无功声强”——能量在来回振荡，并没有向外传播。这个在超材料设计中也要小心。
实际工程中，我们更常用声压级（SPL）和声强级（SIL）：
| 物理量 | 定义 | 参考值 |
|---|---|---|
| 声压级 SPL | Lp = 20·log₁₀(p/p_ref) | p_ref = 2×10⁻⁵ Pa |
| 声强级 SIL | LI = 10·log₁₀(I/I_ref) | I_ref = 10⁻¹² W/m² |
1.3 声阻抗——介质的“脾气”
声阻抗 Z，定义为声压与质点速度的比值：
Z = p / v对于平面波在自由场中传播，声阻抗就是特性阻抗：
Z₀ = ρc为什么说它是介质的“脾气”？你想想看：
- 阻抗大的介质（比如水，ρc ≈ 1.5×10⁶ Rayl），同样的声压只能产生很小的质点速度——介质“不配合”
- 阻抗小的介质（比如空气，ρc ≈ 415 Rayl），同样的声压能产生很大的质点速度——介质“很配合”
这个特性在超材料设计中至关重要。我记得有一次做水下声学隐身衣，核心思路就是通过超材料结构，让空气和水之间的阻抗渐变，从而减少反射。说白了，就是给声波铺一条“平滑的路”。
关键点：声阻抗不匹配 → 反射；声阻抗匹配 → 透射。超材料阵列优化的核心目标之一，就是实现“阻抗渐变”或“阻抗定制”。
1.4 声波在介质中的传播特性
声波在介质中传播时，会遇到四种基本现象：反射、折射、衍射、衰减。
1.4.1 反射与折射
当声波遇到两种介质的界面时，一部分能量反射回去，一部分透射过去。反射和透射的比例由声阻抗决定：
反射系数 R = (Z₂ - Z₁) / (Z₂ + Z₁)
透射系数 T = 2Z₂ / (Z₂ + Z₁)注意，这里说的是声压的反射和透射系数。能量上，反射率 = R²，透射率 = 1 - R²（无损耗时）。
我曾经做过一个项目，要在钢板上贴一层吸声材料。一开始随便选了一种泡沫，结果反射率还是很高。后来一算，泡沫的阻抗和空气差太多，根本匹配不上。换了阻抗渐变的多层结构，效果立竿见影。
1.4.2 衍射
声波遇到障碍物时，会绕过障碍物继续传播。这就是衍射。衍射的强弱取决于障碍物尺寸与波长的比值：
- 障碍物尺寸 ≪ 波长：声波几乎不受影响，完全绕过去
- 障碍物尺寸 ≈ 波长：衍射明显，障碍物后方出现“声影区”
- 障碍物尺寸 ≫ 波长：衍射很弱，障碍物后方形成明显的声影区
这个特性在超材料设计中经常被利用。比如，我们可以用亚波长结构（尺寸小于波长）来控制声波，利用的就是衍射效应。
1.4.3 衰减
声波在介质中传播时，能量会逐渐衰减。衰减的原因主要有两个：
- 黏滞损耗：介质黏性导致质点运动时产生摩擦
- 热传导损耗：声波压缩和膨胀时，热量在介质中传导
衰减系数 α 的单位是 Np/m（奈培每米）或 dB/m。对于空气，衰减与频率的平方成正比：
α ≈ 1.8 × 10⁻¹¹ · f² (dB/m, 20°C, 50%湿度)这意味着高频声波衰减很快。我做过一个实验，20kHz 的超声波在空气中传播 10 米，衰减超过 30dB。所以，做超材料设计时，频率范围决定了你的结构尺寸和材料选择。
避坑指南：我曾经在仿真中忽略了空气的黏滞损耗，结果实验测出来的吸声系数比仿真低了 20%。后来在 COMSOL 中加入了“热黏性声学”模块，才和实验对上。记住：高频仿真一定要考虑衰减。
本章小结
好了，这一章的内容就到这里。我们回顾了四个核心概念：
- 波动方程：描述声波传播的数学工具
- 声压与声强：衡量声音“大小”和“能量”的物理量
- 声阻抗：决定声波在界面处行为的“脾气”参数
- 传播特性：反射、折射、衍射、衰减——声波在介质中的“行为模式”
这些概念就像声学世界的“字母表”。后面我们做超材料阵列优化，所有的设计、仿真、实验，都离不开它们。下一章，我们会进入超材料的世界，看看怎么用这些基础知识去“操控”声波。
有什么问题，欢迎随时交流。我的微信是 deep3321，公众号是“蓝海资料掘金营”。咱们下章见。