声学载荷到底是什么？——Abaqus 结构-声学耦合学习笔记（五）

继续学习 Abaqus Structural-Acoustic Analysis。前面主要看了声学材料和声学边界条件：空气区域需要定义密度和体积模量，声波进入模型以后，还要通过非反射边界尽量“自然离开”有限计算区域。

这一部分继续往下看一个更基础的问题：在声学有限元中，什么叫边界条件？什么叫载荷？为什么在结构分析里很熟悉的“位移—力”关系，到了声学分析中会变成“声压—体积加速度”关系？以下仍然是个人学习记录，理解不一定完全准确，欢迎交流指正。

从结构静力学里的      说起

图1展示的是一个很熟悉的结构有限元关系式：    。其中，     是材料刚度矩阵，     是位移向量，     是载荷向量。这个公式本身并不复杂，但放在“Acoustic loads and boundary conditions”这一节里，其实是为了说明一个更一般的概念：边界条件和载荷的区别，取决于基本求解变量。

在结构静力学中，基本求解变量是位移     。如果直接规定某个节点或某个面的位移，例如固定端     ，这就是边界条件。与位移对应的共轭量是力     ，所以集中力、面力或体力就是结构载荷。平时做结构分析时，这件事很自然，以至于很少专门思考。

但到了声学分析中，这个直觉不能简单照搬。因为声学单元的基本变量已经不是结构位移，而是声压。也就是说，在声学问题里，判断一个设置是边界条件还是载荷，不能再沿用“位移—力”的习惯，而要重新看它对应的是声压本身，还是声压的共轭量。

这也是图1的意义。它不是为了重新讲结构力学，而是先用最熟悉的      建立参照：基本变量被直接指定，就是边界条件；与基本变量共轭的量被施加，就是载荷。

声学问题中的基本变量变成了声压

图2把这个逻辑转到声学问题中。对于声学介质，可以写出关系式     。这里      是声压，     是介质密度，     是介质质点加速度。严格说，这是一个向量关系式，因为      是声压梯度向量，     也是加速度向量。在三维直角坐标中，它可以理解为     、    、    。

这个式子说明，声压的空间变化会驱动空气质点产生加速度。压力高的一侧和压力低的一侧之间存在压力梯度，空气质点就会被推动。也就是说，声学问题里同时存在声压波动和介质质点运动，它们不是孤立的。

Abaqus 声学单元的基本自由度是节点声压，也就是 DOF 8。因此，直接指定声压就是声学边界条件。例如前面消声器入口处的 prescribed acoustic pressure，本质上就是对基本求解变量      的规定。

那么声学载荷是什么？这一页给出的答案是：声压      的共轭量可以理解为体积加速度。也就是说，声学载荷不是我们在结构分析中熟悉的普通机械力，而是与声学介质局部压缩、膨胀和质点加速度相关的激励。

图2中还特别强调，集中载荷可以用来定义声学介质边界节点上的向内体积加速度。正的集中载荷作用在 DOF 8 上，表示向内的体积加速度。这里很容易误解，因为在结构单元里，*Cload 通常表示节点力；但在声学单元里，同样是集中载荷，作用到 DOF 8 后，物理意义已经变成声学介质边界上的体积加速度输入。

可以用扬声器振膜来理解。振膜向空气一侧运动时，会压缩附近空气，使局部声压升高，并激发声波向外传播。从声学有限元角度看，这类边界运动可以等效为对声学介质输入体积加速度，而不是简单地给空气节点施加一个结构意义上的力。

incident wave：一列已经传播过来的声波

图2还提到，声学单元的分布载荷可以通过 incident wave loading 施加。这里的 incident wave 通常译为入射波，指的是从声源传播过来、尚未与目标结构发生相互作用的声波。

入射波可以用不同物理量描述。如果从介质质点运动出发，可以写成位移形式，例如     。但在声学有限元中，由于基本变量是声压，更自然的写法通常是声压形式，例如     。两种写法描述的是同一列声波，只是选用的变量不同。位移、速度、加速度和声压之间可以通过声学关系联系起来。

incident wave loading 的意思，并不是在某一个节点上施加一个力，而是给定一列具有幅值、相位和传播方向的入射声场。这个声场作用到结构附近时，会撞击结构表面，产生反射、绕射、散射，也可能从结构周围继续传播过去。原文中 “impinging on and passing over a structure” 描述的正是这个过程：声波从声源传来，作用到结构上，并绕过结构继续传播。

从量纲上也可以理解为什么入射波可以作为声学分布载荷。声压      的量纲是 Pa，也就是 N/m²；作用在结构表面时，它本身就是分布压力。另一方面，声压梯度      的量纲是 Pa/m，也就是 N/m³，与      的量纲一致，对应单位体积上的惯性作用。因此，一列已知入射声场可以在有限元方程中转化为分布式声学激励。

一点学习体会

这两页内容让我感觉，声学载荷和结构载荷看起来在软件操作上有相似之处，但背后的物理含义并不一样。结构分析中，我们习惯了“位移—力”的关系；声学分析中，基本变量变成了声压，对应的载荷也要从声压场、压力梯度和介质质点运动来理解。

所以，指定声压是边界条件，DOF 8 上的集中声学载荷可以表示向内体积加速度，而 incident wave loading 则用于描述一列外部声波作为分布声学激励作用到模型中。理解这一点后，再看声学模型里的 prescribed acoustic pressure、volumetric acceleration 和 incident wave，就不会觉得它们只是几个软件选项，而是声学有限元变量体系中的不同入口。

以上内容是个人学习 Abaqus 结构-声学耦合分析过程中的整理，属于学习记录性质。部分理解可能还比较粗浅，欢迎相关领域老师和同行交流指正。