声学分析的基本方程是什么？

线性声学的基础源于连续方程、动量守恒和状态方程，通过假设无粘非等熵流，将压力和密度分解为平均值和微小扰动，线性化后得到赫尔姆霍兹方程（∇²p + k²p = 0），其中p是声压，k是波数。

音速是如何确定的？

音速由媒质的体积弹性率和密度决定。对于气体，c₀ = √(γRT)（γ:比热比，R:气体常数，T:绝对温度）。实务中常用值为空气(20℃)343m/s、水约1482m/s、钢材约5900m/s。

声学分析中处理的声压大小是多少？

声压p相对于大气压P₀(101325Pa)非常小，会话水平约为0.02Pa~2Pa。其比值p/P₀约为10⁻⁷，极小，因此方程线性化成为可能。

CAE声学分析中常用的软件是什么？

Ansys Mechanical的Acoustics模块等被广泛使用。这些工具基于赫尔姆霍兹方程进行数值分析，应用于噪声预测和声学设计。

声学分析 — CAE术语解释

分类：术语集 | 2026-01-15

CAE visualization for acoustic analysis - technical simulation diagram

声学分析

🧑‍🎓

声学分析就是"声音的分析"吧。比如汽车车内的静音性和扬声器的指向性？用FEM可以计算声音，真是令人惊讶。

声学的理论基础

声学波动的基本方程

🧑‍🎓

声学分析最基本的控制方程是什么？是如何从流体运动方程导出的？

🎓

线性声学的基础源于连续方程、动量守恒和状态方程。假设无粘非等熵流，将压力和密度分解为平均值和微小扰动，线性化后得到如下赫尔姆霍兹方程：

$$ \nabla^2 p + k^2 p = 0 $$

其中，

$$ p $$

是声压（扰动分量），

$$ k = \omega / c_0 $$

是波数，

$$ c_0 $$

是音速（空气中20℃约为343 m/s），

$$ \omega $$

是角频率。

🧑‍🎓

"线性化"具体是怎样的操作？为什么可以忽略非线性项？

🎓

很好的疑问。例如，将压力分解为

$$ P = P_0 + p $$

。其中

$$ P_0 $$

是平均压力（101325 Pa），

$$ p $$

是声压（通常会话水平约0.02 Pa～2 Pa）。声压相对于平均压力极小（

$$ p / P_0 \sim 10^{-7} $$

）。因此，运动方程中的对流项

$$ (\mathbf{v} \cdot \nabla) \mathbf{v} $$

和状态方程的高阶项都可以忽略。这就是线性近似的物理根据。当声压级超过约140 dB（

$$ p \sim 200 Pa $$

）时，非线性效应不再可忽略。

🧑‍🎓

音速

$$ c_0 $$

随媒质的变化规律是什么？分析中常用的值是多少？

🎓

音速由媒质的体积弹性率和密度决定。气体中

$$ c_0 = \sqrt{\gamma R T} $$

（

$$ \gamma $$

是比热比，

$$ R $$

是气体常数，

$$ T $$

是绝对温度）。实务中常用的默认值为空气（20℃）343 m/s、水（20℃）约1482 m/s、钢材约5900 m/s。Ansys Mechanical的Acoustics模块和COMSOL的Pressure Acoustics模块都预设了这些材料库值。

声学的数值计算手法

FEM的离散化和边界条件

🧑‍🎓

用有限元法（FEM）求解赫尔姆霍兹方程时，如何进行离散化？特别是含波数k的项如何处理？

🎓

用伽辽金法推导弱形式。在控制方程两边乘以权函数

$$ w $$

进行积分，分部积分后得到下列形式：

$$ \int_{\Omega} (\nabla w \cdot \nabla p - k^2 w p) d\Omega - \int_{\Gamma} w \frac{\partial p}{\partial n} d\Gamma = 0 $$

第二项的边界积分是引入边界条件的关键。用形状函数近似节点变量p，离散化后最终得到周频率依赖的线性方程

$$ (\mathbf{K} - \omega^2 \mathbf{M}) \mathbf{p} = \mathbf{f} $$

其中

$$ \mathbf{M} $$

是相当于"质量矩阵"的矩阵，在声学中习惯上这样称呼。

🧑‍🎓

边界积分中的

$$ \frac{\partial p}{\partial n} $$

具体如何设置？"完全反射"和"吸收"如何区分？

🎓

这就是边界条件的本质。硬墙（完全反射）时，法向粒子速度为0，即

$$ \frac{\partial p}{\partial n} = 0 $$

（诺伊曼条件）。而贴着吸音材料的墙或模拟无限远的"吸收边界条件"，使用局部关系式

$$ \frac{\partial p}{\partial n} = -ikp $$

（阻抗条件或一阶吸收边界）等。为获得更高精度的吸收，使用PML（完全匹配层）。Siemens Simcenter 3D的Acoustics模块可以将吸音材料特性作为周频率依赖的复阻抗直接输入。

🧑‍🎓

网格尺寸如何确定？需要根据频率改变吗？

🎓

声学FEM需要能充分分辨波长的网格。经验规则是，针对解析的最高频率

$$ f_{max} $$

，该频率下的波长

$$ \lambda_{min} = c_0 / f_{max} $$

需要每波长6～10个以上单元。例如，空气中解析到5000 Hz，

$$ \lambda_{min} \approx 0.0686 m $$

，单元尺寸应小于约6.9～11.4 mm。Abaqus Acoustics指南推荐"1波长用6次单元1个单元，2次单元3个单元"。网格过粗会产生"数值分散"，使计算声速变慢。

声学的实务应用

室内声学分析的工作流程

🧑‍🎓

如何用CAE评估会议室等室内的声学特性（残响时间等）？工作步骤是什么？

🎓

典型工作流程如下：1) 用3D CAD制作室内形状。2) 按材料分别设置声阻抗（吸音率）。使用JIS A 1409或ISO 354规范的实测值。3) 网格生成。按照上述规则，例如500 Hz解析则单元尺寸～70mm以下。4) 定义音源。可使用点音源或面振动源。5) 执行频率响应分析。一般使用125、250、500、1000、2000、4000 Hz的1/1倍频程中心频率。6) 后处理评估声压分布、特定点的频率响应、及残响时间（从萨宾公式或能量衰减曲线计算）。

🧑‍🎓

设置吸音率时需要注意什么？可以直接输入"吸音率0.5"吗？

🎓

绝对不可以。吸音率是实数，但声阻抗一般是复数。虽然存在从吸音率反演阻抗的近似公式，但特别在低频时误差很大。实务中应该直接输入材料厂家提供的"复特征阻抗"或"复反射系数"的频率数据。例如吸音材料厂家的数据表中，125-4000Hz范围内的实部和虚部已列出。COMSOL的"阻抗边界条件"可以将这个实部和虚部用表格输入。

🧑‍🎓

分析结果的残响时间与实测不符时，首先应该怀疑什么？

🎓

按下列检查清单逐一确认：1) **材料数据**：阻抗值是否在正确的频率输入。特别是低频（125Hz）的数据不确定性较高。2) **网格**：是否针对最高频率足够细密。粗网格会过大评估衰减。3) **泄漏**：门隙等使声能逸出的路径是否建模。4) **家具影响**：会议室的桌椅会散射和吸收声。这些可能需要作为等效吸音面积添加。5) **求解器设置**：衰减（如果有）的设置或求解频率点数是否充分。

声学的软件比较

主要软件的特点和适用范围

🧑‍🎓

用Ansys、COMSOL、Abaqus做声学分析，分别有什么特点和优势领域？

🎓

各有特色。**Ansys**主要是"Mechanical + ACT Acoustics"或"Harmonic Acoustics"。结构振动连成（FSI）强，在汽车车内噪声和家电振动噪声方面很强。大规模模型的并行计算性能优异。**COMSOL Multiphysics**的"Pressure Acoustics"模块为核心，多物理场连成是本领。可在同一环境灵活设置声-结构-热-流耦合。适合研发。**Abaqus**用"Abaqus/Standard"的声学单元。与结构分析无缝集成，可处理轮胎路面噪声和水中声学等含固液大变形复杂接触问题。

🧑‍🎓

开源免费软件（OSS）如何？能实用吗？

🎓

有限但可行。**CalculiX**和**Code_Aster**等FEM求解器有声学单元。前处理、后处理用**Salome-Meca**、**Gmsh**、**ParaView**组合。但商用软件的高级吸收边界条件（PML）、复杂周频率依赖阻抗设置、高效周频扫频求解器需要自己实现或找现有宏。教学和小规模各向同性模型验证有用。在实务产品开发中，商用软件因支持和已验证功能更有优势。

🧑‍🎓

还听说"用BEM（边界元法）做声学分析"。与FEM有什么区别？

🎓

根本区别是，FEM对全域进行网格分割，而BEM只对边界进行网格化。这样能自然地表现无限远声场（辐射声），是最大优点。适合汽车外部辐射噪声和扬声器远场声响分析。代表性商用BEM求解器有**Altair Acoustics**（原LMS Virtual.Lab）、**FFT Actran**（FEM/BEM两者可用）、**COMET**等。缺点是系数矩阵为稠密矩阵，大规模模型计算成本高。近年来高速BEM（快速多极子BEM等）在克服这一问题。

声学的故障处理

常见收敛错误和对策

🧑‍🎓

频率响应分析中，特定频率（如275Hz）解发散或出现异常大声压值。原因是什么？

🎓

几乎肯定是"固有模式共振"。声学腔室有自己的固有振动数（声学模式），像鼓膜一样。如果减衰没有充分建模，该频率响应理论上会趋向无穷。首先，执行该频率附近的固有值解析，确认固有频率是否一致。对策有3种：1) 在模型中加入现实的减衰（材料内部损失和边界吸收）。2) 进行结构-流体耦合分析（现实的墙不是完全刚体）。3) 评估结果时回避共振频率，或进行平均。Ansys可设置"Modal Damping"或"Rayleigh Damping"有效。

🧑‍🎓

细化网格后反而结果出错，求解器报"Matrix is singular"。为什么？

🎓

这是"数值锁定"或"虚假模式"的典型症状。声学FEM多用仅压力为变量的定式（压力公式）。这类似于非压缩流FEM的"体积锁定"，网格过细时单元层面会过度满足

$$ \nabla \cdot \mathbf{u} = 0 $$

（非压缩条件），刚性矩阵趋向奇异。对策是1) 降低单元阶数（2次到1次）。2) 用混合定式（压力和位移势为变量）的软件。3) 避免网格纵横比过极端。Abaqus可用"AC3D8"等线性声学单元来规避。

🧑‍🎓

设了PML（完全匹配层），边界仍有明显反射。原因是什么？

🎓

PML设置有误。检查以下项目：1) **层厚度**：一般需要目标最低波长的1/2～1波长厚度。500Hz下空气中约需34cm以上。2) **衰减系数分布**：从内部域到PML边界应平滑递增（二次函数等）。陡峭会反射。3) **PML配置**：需垂直于波传播方向。复杂形状或散射体应球形或圆柱形配置。4) **网格**：PML区内也需至少3～5层单元。检查COMSOL"scaling coordinate system"或Ansys"PML单元类型"设置是否未改变默认值。

🧑‍🎓

结构-声学耦合分析中，流体（空气）网格和结构网格的接口如何处理？需要一致吗？

🎓

不一定要一致，但要用一致方法耦合。主流有2种方法：1) **一致网格**：接口共享节点。最简单精度高，但网格生成灵活性受限。2) **非一致网格**：分别独立网格，用约束方程（MPC）或拉格朗日乘数法课加"压力和法向位移连续条件"。Ansys"FSI接口"和Abaqus"TIE"约束即此。用非一致网格时，若声学网格比结构网格粗很多，则高频能量传递会不正确，所以接口附近声学网格应与结构网格同级或更细。