对称面为什么可以“不设声压”？——Abaqus 结构-声学耦合学习笔记（七）
继续学习 Abaqus Structural-Acoustic Analysis。前面已经梳理了声学边界条件的三种基本形式：指定声压、自然边界和阻抗边界。这次看到的两个例子，一个是卡车驾驶室的半对称声学模型，一个是消声器入口的指定声压，正好把前面比较抽象的第一类、第二类边界条件落到了具体模型中。尤其是消声器入口设置中的
REAL，还涉及稳态谐波分析中一个非常重要的概念——复声压。
半个驾驶室为什么能代表完整声场
图1展示的是一个卡车驾驶室的半对称模型。结构域只建立了一半，包围驾驶室的声学域同样只保留了一半，中间切开的平面就是 acoustic symmetry plane，也就是声学对称面。在这个平面上采用 Natural boundary condition，也就是自然边界条件。
前面已经知道，声学介质满足动量关系 。沿边界法向投影后，可以写成：
自然边界条件对应 ，因此有 。对于真实的静止刚性壁面，这意味着空气不能沿法向穿过壁面；而对于图1中的声学对称面，物理意义稍有不同。假设完整声场关于中间平面对称，那么对称面两边相同位置的声压应该完全相等。声压从一侧跨过对称面以后呈镜像变化，因此在正中央位置沿法向的斜率恰好为零，也就满足 。
这让我重新理解了“自然边界条件”这件事：相同的数学边界条件，并不一定代表相同的物理边界。 静止刚性墙和声学对称面都满足 ，但前者是真的存在一堵不可穿透的墙，后者只是人为把完整模型切成了一半，另外半边声场可以通过镜像恢复。
当然，能不能使用这种半模型，不能只看几何外形是否对称。结构、材料、约束和激励也必须满足相应的对称性。例如左右车身完全一样，但如果只在左侧施加一个不对称激励，那么结构振动和声场已经不再关于中心面对称，这时再使用半模型就不合适了。
消声器入口为什么直接规定声压
图2又回到了前面的消声器声辐射算例。
这里选中了内部空气入口 INLET，并直接规定节点声压 DOF 8：
INLET, 8, 8, 1.0E-6
声学单元的基本求解变量就是声压 ，所以这种设置本质上就是直接规定：
这对应一般偏微分方程中的第一类 Dirichlet 边界条件。它和结构分析中直接规定位移 是同样的逻辑。
这个模型采用 mm–tonne–s 单位体系，压力单位是 MPa，因此 MPa 正好等于 1 Pa。也就是说，消声器内部空气入口处规定的是 1 Pa 的声压。
这里很容易把它理解成“在入口施加了一个声学力”，实际上并不是。声压本身就是声学单元的基本自由度，因此它属于边界条件。入口声压规定以后，入口和内部空气之间形成声压的空间梯度，再根据 推动空气质点运动，声波才进一步向消声器内部传播。
所以这一过程可以简单理解成：先规定入口压力状态，再由压力梯度驱动介质运动。
为什么这里又出现了 REAL
图2中真正值得多想一步的是关键词：
*BOUNDARY, REAL
为什么一个声压边界还需要特别说明 REAL？
首先要明确，真实世界中的声压始终是实数。例如麦克风测到的 ，在任何一个具体时刻都是一个真实的压力值。这里出现复数，并不是说物理声压变成了复数，而是因为当前求解的是稳态谐波问题。
假设某个位置的声压是一个单频简谐信号，可以写成：
这里其实包含三个信息：幅值 、相位 和频率 。为了让后续动力方程处理起来更方便，通常把它写成复数形式：
其中 称为复声压幅值。这个复数保存了两个信息：它的模 表示声压幅值，它的辐角 表示相位。
这里还有一个很关键的问题：频率并没有包含在 这个复数里。 频率由 单独描述，或者写成普通频率 ，二者满足 。
所以稳态声场更完整的表达应该是：
其中 描述“以多快的频率振动”，而 描述“在这个频率下，这个位置振多大、相位是多少”。
这也是稳态谐波分析和一般时域分析很不一样的地方。稳态分析已经预先假设系统最终只剩下与激励相同频率的周期响应，所以没有必要一步一步追踪每个时刻的声压，只需要在指定频率下求复幅值 即可。
1 Pa、REAL，到底代表了什么
当前消声器算例后面只分析 170 Hz，因此入口的完整物理条件其实可以理解成一个 170 Hz 的谐波声压。
关键词中给出 REAL=1.0E-6 MPa，而没有给虚部，因此对应：
它的幅值是 1 Pa，相位是 。因此这个入口可以看作整个稳态声学问题的相位参考。
恢复到真实时域声压，就是：
这里的 1 Pa 和 来自复声压边界，而 170 Hz 来自后面的稳态动力分析步。
声波进入消声器以后，虽然所有位置仍然以 170 Hz 振动，但由于传播距离、反射、结构振动以及声学共振等作用，不同位置的幅值和相位都会发生变化。例如某个位置可能得到 Pa。这个复数不是另一种频率，而是表示同样在 170 Hz 下，这里的声压幅值和相位已经与入口不同。
这也是为什么复数形式在声学问题中 特别方便。对于一个指定频率，整个复杂声场最终都可以表示成一个空间分布的复数场：
一个复数同时保存幅值和相位，而频率则作为独立的分析坐标存在。如果继续计算 100 Hz、150 Hz、200 Hz 等多个频率，就会分别得到这些频率下不同的复声压场，最后也就形成我们熟悉的声压频率响应。
一点学习体会
这两页看起来分别只是“对称面”和“入口声压”两个简单设置，但放到前面学过的边界条件框架里以后，逻辑其实非常清楚。
声学对称面对应第二类边界条件，控制的是声压沿法向的变化，满足 ；消声器入口则对应第一类边界条件，直接规定基本变量 。
而 REAL 又进一步提醒我：进入稳态谐波分析以后，求解对象已经不是完整的 时间历程，而是在每个指定频率下求一个复声压 。频率决定“振多快”，复数决定“振多大以及相位是多少”。
把这一点弄清楚以后，再看到稳态声学分析中的 Real、Imaginary、Magnitude 和 Phase，就会发现它们其实只是在用不同方式描述同一个谐波声场。
以上内容是个人学习 Abaqus 结构-声学耦合分析过程中的整理，属于学习记录性质。部分理解可能还比较粗浅，欢迎相关领域老师和同行交流指正。
