声波到达边界以后会发生什么？——Abaqus 结构-声学耦合学习笔记（六）
继续学习 Abaqus Structural-Acoustic Analysis。前面已经梳理了声学载荷的基本含义：声学单元以声压为基本求解变量，指定声压属于边界条件，而体积加速度和入射波则对应不同形式的声学激励。
这一部分继续往下看声学边界条件。声波传播到边界以后，是被反射、被吸收，还是继续离开计算区域，取决于边界条件的定义。把这个问题和偏微分方程中常见的三类边界条件联系起来，理解会更清楚。以下仍然是个人学习记录，理解不一定完全准确，欢迎交流指正。
从一般偏微分方程的三类边界条件说起
设区域内的未知函数为 ，边界为 。常见边界条件可以分为三类。
第一类边界条件直接规定未知函数本身：
这就是 Dirichlet 边界条件。第二类边界条件规定未知函数沿边界法向的导数：
这就是 Neumann 边界条件。第三类边界条件同时包含未知函数及其法向导数：
这就是 Robin 边界条件，也常称为混合边界条件。
声学有限元中的基本求解变量是声压 ，因此图1中的 Prescribed acoustic pressure、Natural boundary condition 和 Impedance，实际上可以分别放到第一类、第二类和第三类边界条件中理解。
声学介质满足局部动量关系 。将其投影到边界法向上，有：
其中， 是介质质点沿边界法向的加速度。这样一来，三类声学边界条件就可以统一理解为：规定声压 、规定法向运动 ，或者规定声压与边界运动之间的关系。
Prescribed acoustic pressure：第一类边界条件
Prescribed acoustic pressure 是指定声压边界条件，对应 Dirichlet 条件。它直接给出边界上的声压：
例如在消声器入口处，可以规定一个随时间变化的压力脉动。此时边界声压不再由方程求解，而是由用户直接给定。这和结构分析中直接规定位移 是同一种思路。
从 的角度看，指定声压并不是直接给空气一个加速度，而是先规定边界压力。边界与内部之间形成声压梯度后，附近空气质点才会被推动。
所以，第一类边界条件回答的是：这里的声压是多少。
Natural boundary condition：第二类边界条件
Natural boundary condition 对应 Neumann 条件。它不直接规定声压，而是规定声压沿边界法向的变化。
对于静止刚性壁面，空气不能沿法向穿过边界，因此 。代入 ，可得：
这并不表示边界声压为零，而是表示声压沿法向没有梯度。声波到达刚性壁面后不能继续穿过，通常会被反射回来。
声学对称面也满足同样的数学形式。虽然对称面和刚性壁面的物理含义不同，但它们都对应法向声压梯度为零。
这一点很值得注意：声学边界上什么都不设置，并不表示开放空间。自然边界通常是默认条件，因此“什么都不设”往往意味着刚性反射。如果把本应开放的外边界留空，声波就可能被人为反射回来，形成不真实的驻波和声压放大。
所以，第二类边界条件回答的是：这里的空气沿法向怎样运动。
Impedance：第三类边界条件
Impedance，也就是声学阻抗，对应第三类或混合边界条件。它既不单独规定声压，也不单独规定法向运动，而是建立二者之间的关系。
对于最简单的阻性边界，可以写成 ，其中 是声学阻抗， 是介质质点的法向速度。又因为 ，而边界法向上的动量关系为：
将 代入，可得：
整理后为：
这个式子虽然是时域形式，但已经能看出阻抗边界的“混合”特征：它把声压随时间的变化与声压沿边界法向的空间变化联系起来。也就是说，边界既没有被完全固定，也没有直接指定声压，而是要求声压和边界运动按照一定关系共同变化。
如果阻抗很大，边界几乎不动，效果接近刚性反射；如果阻抗与介质传播特性匹配，声波更容易离开边界；如果阻抗中包含耗散作用，一部分声能会在边界处被吸收。
图1中提到，阻抗有两个主要用途。一类用于外部声学区域，使声能通过非反射边界离开有限网格；另一类用于结构—声学界面，描述声学介质与固体之间的能量损失。对于平面波，可以采用平面非反射边界；对于圆柱、椭圆等简单区域，还可以采用与几何形状相匹配的处理。
图中的三类边界应该怎样理解
图1虽然只是一个分类式概览，但实际上给出了一个很实用的判断思路。
如果边界代表静止刚性壁面或对称面，就采用自然边界，对应第二类条件；如果边界上的声压已经知道，就采用指定声压，对应第一类条件；如果希望声波离开有限区域，或者需要考虑界面吸声和耗散，就采用阻抗边界，对应第三类条件。
也可以把三类边界理解成声波到达边界后的三种反应：自然边界让声波反射；指定声压边界直接规定边界压力状态；阻抗边界则允许声压和边界运动相互匹配，从而产生反射、吸收或向外传播。
图中还强调，同一片区域通常只应施加一种边界条件，结构—声学界面除外。因为三类条件约束的是同一个声压场，如果同一位置既直接规定声压，又规定阻抗关系，就可能造成重复或矛盾约束。结构—声学界面比较特殊，它不仅要传递结构运动和流体压力，还可能通过界面阻抗描述能量损失。
一点学习体会
把声学边界条件和偏微分方程的三类边界条件对应起来以后，很多设置就不再显得零散。
指定声压，本质上是直接规定基本变量 ；自然边界规定的是声压法向梯度，等价于限制法向介质运动；阻抗边界则把声压与边界运动联系起来。它们分别对应第一类、第二类和第三类边界条件。
这部分学习对我最大的提醒是：声学边界条件不是简单决定“边界在哪里”，而是在决定声波到达边界以后会发生什么。尤其对于外部声场，如果把开放边界误设成自然边界，本应离开的声波就会被反射回来。边界条件一旦判断错误，后面的声压结果再精细，也可能失去真实意义。
以上内容是个人学习 Abaqus 结构-声学耦合分析过程中的整理，属于学习记录性质。部分理解可能还比较粗浅，欢迎相关领域老师和同行交流指正。
