音响-结构相互作用
概述
老师，音响-结构相互作用就是说音和结构互相影响吗？
音响-结构相互作用的理论基础
正是如此。当结构物体振动时，会推动周围的空气或流体产生声波。反过来，声压会作为载荷作用在结构上，引起振动。这种双向耦合就是音响-结构相互作用（ASI）。
具体在哪些场景会成为问题呢？
汽车车室内噪声是典型例子。发动机的振动激励车体面板，面板压缩和膨胀车室内的空气，乘客听到的就是这个声音。其他还有飞机机身透射噪声、建筑物隔音设计、潜水艇声呐罩设计等。
支配方程
结构侧和音响侧分别用什么方程描述呢？
结构侧通常用弹性体的运动方程。
右边的 $[A]\{p\}$ 是声压对结构的载荷项，$[A]$ 是耦合矩阵。音响侧是Helmholtz方程的时间域版本，即波动方程。
这两个方程在哪里连接呢？
在耦合界面的边界条件处连接。从音响侧看，界面法向速度必须与结构的变位速度一致。
写成矩阵形式，耦合系统整体变成非对称方程组。
非对称性看起来很麻烦。
确实如此。Nastran中使用FLUID单元（CAABSF、CHEXA等）会自动内部组织这个非对称耦合系统。Abaqus中为音响-结构耦合提供了*TIE约束。
周频域的定式化
振动噪声通常用频率讨论，对吧？
假设调和响应 $\{u\} = \{\hat{u}\}e^{i\omega t}$、$\{p\} = \{\hat{p}\}e^{i\omega t}$，可以变换到频率域。
这样可以逐频率求解，高效地得到频率响应函数（FRF）。
这样的定式化能清楚看出哪个模式在哪个频率有问题，对吧。
正是。特别是结构的固有频率与音响腔体的固有频率接近时，会发生耦合共振，声压会被放大。汽车车室的"轰鸣"问题就是这种典型情况。
实务上的要点
分析时需要注意什么？
首先是音响网格的单元尺寸。对于最高分析频率 $f_{max}$，要保证每个单元内至少有6个节点在波长内。
在空气中 $c = 343$ m/s，$f_{max} = 1000$ Hz 时，$h \leq 57$ mm。二阶单元可以略粗一些，但一阶单元应该坚守这个规则。
看来低阶单元需要非常细的网格呢。整体情况逐渐清楚了。
"声音能改变结构的形状"——音响辐射压现象
学习音响-结构相互作用理论时，很多人对"声压对结构施加力"的事实感到震惊。其实超声波清洗机就是这个原理的应用，28~40 kHz的声波在金属部件表面产生微小气泡（空化），物理上剥除污垢。声压级可超过150 dB。人类听不到的声音，却能改变金属形状——音响耦合理论就是要用方程式描述这种"看不见的力"的本质。
音响-结构相互作用的数值计算方法
FEM的音响-结构耦合离散化
音响-结构耦合用FEM怎么处理？结构和音响的单元不同吧。
好问题。结构侧用通常的壳单元或实体单元，音响侧用有压力自由度的音响流体单元。耦合界面把结构的变位自由度和音响的压力自由度结合起来。
Nastran里怎么做？
Nastran用FLUID单元（CHEXA、CPENTA、CTETRA的音响版本），在结构面定义ACMODAL面。用SOL 108（频率响应）或SOL 111（模态频率响应）求解耦合系统。
FEM-BEM耦合法
外部辐射问题用FEM会困扰无限域，对吧？
正确。内部问题（如车室内的闭空间）只需FEM，但外部辐射问题要和BEM（边界元法）结合。用FEM求解结构振动，把辐射面的速度作为BEM的边界条件。
BEM的Kirchhoff-Helmholtz积分方程如下。
这里 $G = \frac{e^{ikR}}{4\pi R}$ 是三维自由空间的Green函数。
BEM不需要体网格，对外部问题很适合呢。
是的。Siemens Simcenter 3D（旧LMS Virtual.Lab）和FFT ASTRA是FEM-BEM耦合的代表工具。HEAD acoustics的ARTEMIS也是基于边界元。
统计能量分析（SEA）
高频时FEM的网格不会膨胀吗？
说得好。高频域（汽车行业大概500 Hz以上），模态密度极高，FEM的确定性方法不再实用。这时用统计能量分析（SEA）。
SEA的基本方程是功率平衡。
这里 $\eta_{ij}$ 是耦合损失系数，$n_i$ 是子系统 $i$ 的模态密度，$E_i$ 是能量。
所以要按频率段选择不同方法，对吗？
对。低频用FEM，高频用SEA，中频用混合FE-SEA（Simcenter 3D的Hybrid FE-SEA或Wave支持）。这是振动音响分析实务的基础。
| 频率段 | 方法 | 代表工具 |
|---|---|---|
| 低频（~500 Hz） | FEM / FEM-BEM | Nastran, Abaqus, COMSOL |
| 中频 | 混合FE-SEA | Simcenter 3D, Wave6 |
| 高频（500 Hz~） | SEA | VA One (ESI), AutoSEA |
时间域的显式解法
碰撞声、瞬间音怎么处理？
要在时间域求解。LS-DYNA支持音响-结构耦合的显式解法，用于安全气囊展开声、车门关闭声的瞬态分析。Abaqus/Explicit也能与音响单元耦合。
每个方法的适用范围都清楚了。振动音响的世界真有意思啊。
FEM与边界元法——音响耦合分析的"搭档"选择
选择音响-结构相互作用的数值解法时，现场工程师必然会纠结："只用FEM求闭合系统，还是加BEM来处理外部音场？"汽车车室噪声这种封闭空间的情况下FEM很有利，但算车外辐射噪声时需要处理无限远边界，这时BEM出场。但BEM会产生稠密矩阵，模型变大时计算成本陡增。近年高速多极法（FMM）的出现缓解了这个问题，但求解器选定仍然是"各有利弊"的局面。
音响-结构相互作用的实务应用
一般流程是这样的。
1. 建立结构模型 -- 从CAD出发构建结构FEM模型。准确建模面板连接和密封部分
2. 建立音响腔体模型 -- 用音响单元网格化车室等闭空间
3. 定义耦合界面 -- 设置结构面和音响面的接口（Nastran用ACMODAL，Abaqus用TIE）
4. 设置激励条件 -- 设定发动机挂架力、路面激励力等输入
5. 执行频率响应分析 -- 在目标频率范围内计算FRF
6. 后处理 -- 音压级分布、贡献分析、面板贡献分解
第3步的耦合界面定义是关键，听说过。
正是。界面的法向必须一致，不一致的话能量传递不会正确计算。Nastran中把wet surface的法向统一是标准检查点。
网格设计的诀窍
结构和音响的网格尺寸考虑方法不一样吧？
音响网格用波长准则 $h \leq c/(6 f_{max})$。结构网格用弯曲波长准则 $h \leq \lambda_b / 6$。弯曲波长用下式求得。
这里 $D$ 是弯曲刚度，$t$ 是板厚。通常结构侧需要更细的网格。
结构和音响的网格不一致时怎么办？
就变成非匹配网格耦合。Nastran的ACMODAL和Abaqus的TIE都支持非匹配网格。但网格差异太大的话耦合精度会下降，界面附近最好控制在2:1左右的比率。
面板贡献分析
怎样找出哪个面板最吵？
用面板贡献分析（Panel Contribution Analysis）。把接收点的音压分解成各面板的贡献。面板 $j$ 的贡献可以写成。
由于包含位相和振幅，有些面板的贡献可能相互抵消。Siemens Simcenter 3D的贡献分析功能很有名。
改进措施的优先顺序就能定下来了。
常见失败与对策
| 症状 | 原因 | 对策 |
|---|---|---|
| 特定频率出现异常的音压峰值 | 耦合共振（结构与音响固有频率一致） | 结构侧固有频率调整（加强度）或加制振材 |
| 实测和分析的频率不符 | 材料常数（尤其杨氏模量）误差 | 用实验模态同定的值 |
| 全频段音压过大 | 阻尼设定不足 | 适当设定结构阻尼（$\eta_s$=0.01~0.03）和音响吸声率 |
| 高频分析结果波动 | 网格过粗（波长分解度不足） | 检查"6个单元/波长"规则 |
阻尼参数设定看起来很难。
实务中结构阻尼系数按面板来设定，座椅和地毯等的吸声率作为阻抗边界条件输入。这个调整对分析精度影响很大。
车门关闭音是"音响-结构相互作用"的浓缩
高级车的车门关闭时不是"砰"的声音，而是"咚"——调整这种音质是一门学问。车门关闭音其实是金属面板的振动（结构）→车室空气振动（音响）→密封件减衰（材料）的相互作用的体现。实务中改变面板厚度0.1 mm就能改变音质，不通过试制直接从仿真就能达到目标声音，音响耦合分析在这里发挥了重要作用。有专门的"车门声设计工程师"这样的职位，足见其重要性。
音响-结构相互作用的软件比较
主要工具的对比如下。
| 工具名 | 开发商 | 擅长领域 | 音响方法 |
|---|---|---|---|
| MSC Nastran | Hexagon (旧MSC) | 汽车NVH | FEM音响、ACMODAL |
| Simcenter 3D | Siemens | NVH集成环境 | FEM-BEM、混合FE-SEA |
| Actran | FFT (Siemens子公司) | 音响专业 | FEM、BEM、无限单元 |
| COMSOL Multiphysics | COMSOL AB | 多物理耦合 | FEM音响、PML |
| Abaqus | Dassault Systemes | 通用非线性 | FEM音响-结构耦合 |
| VA One | ESI Group | 中高频SEA | SEA、FE-SEA |
Actran是第一次听说。
Actran是比利时FFT（Free Field Technologies）开发的音响专业求解器。2019年被Siemens收购。无限单元和自适应PML（Perfectly Matched Layer）功能很充实，外部辐射问题特别强。在航空发动机风扇噪声分析中广泛应用。
NVH业界的工具分工
汽车制造商一般怎样用这些工具？
典型的分工是这样的。
- 低频轰鸣: Nastran SOL 111（模态频率响应） + 音响腔体FEM
- 路面噪声: Nastran + 传递路径分析（TPA）
- 风噪: CFD（Fluent/STAR-CCM+） → Actran or Simcenter BEM
- 中高频: VA One（SEA） or Simcenter Hybrid FE-SEA
Nastran真的是NVH的行业标准呢。
是啊。用DMAP（Direct Matrix Abstraction Program）还能写自定义处理。但最近Simcenter 3D因为有完整的前后处理环境，势头很强。
开源选项
商业工具以外还有什么选择？
OpenFOAM的acousticFoam求解器能做LES+音响类比。Code_Aster+Salome-Meca也能做音响FEM。学术界有pyFEM这样的音响FEM的Python库。但商业工具那样的自动耦合接口需要自己实现。
生产效率上商业工具还是更有优势呢。工具选定的参考很有帮助。
Nastran与Abaqus——音响耦合工具选择的现实
"Nastran还是Abaqus"这个问题在音响-结构耦合分析中根深蒂固。Nastran以DMAP的定制灵活性和轻量化模型见长，在汽车和航空航天的结构部门压倒性的市场份额。而Abaqus在材料非线性和复杂接触条件处理上更强，复合材料和树脂部件多的产品更容易选它。专业音响工具如Actran或VirtualLab在精度方面有优势。但实务选择时，"既有公司许可"和"工程师技能资产"占9成的权重——工具性能差异远不如能否用对来得重要。
音响-结构相互作用的前沿研究
混合FE-SEA
低频FEM、高频SEA这样说了，中频的混合方法具体是怎样的？
混合FE-SEA是这样：要用决定论方法处理的部件（发动机缸体等刚度高、模态少的）作为FE子系统，统计方法处理的部件（面板等模态密度高）作为SEA子系统。
理论上由剑桥大学Langley等人在2005年左右体系化。假定随机点质量，用"扩散场互反定理（diffuse field reciprocity）"实现FE和SEA的一致耦合。
用什么工具？
Simcenter 3D的Hybrid FE-SEA功能和ESI Group的VA One是两大工具。Wave6（旧Wave、Dassault Systemes）也支持。
内饰音响建模
汽车内饰（座椅等）对音响有什么影响？
地毯、仪表板隔音等多层吸隔音材用Biot多孔弹性体理论。支配方程包含3个场（固体骨架变位、流体变位、声压），耦合非常复杂。
Actran的Trim Acoustic Module（TAM）和COMSOL的Poroacoustics Module对此有支持。实务中常用Johnson-Champoux-Allard模型等等效流体模型近似。
所以要测定多孔材料的流阻系数 $\sigma$ 这样的参数，对吧。
机器学习的融合
现在AI的应用怎样了？
有几个方向。一个是音响传递函数的代理模型化。从大量FEM分析结果学习传递函数到神经网络，设计变更时实时预测响应。
二是用PINN（Physics-Informed Neural Network）直接近似波动方程解。精度还有问题，但对逆问题（音源同定）很有前景。
三是与拓扑优化结合。用Adjoint法求声学功率最小的结构形状。COMSOL和OptiStruct都支持。
CAE的世界也在变化呢。不断追赶前沿很重要。
音响超材料——"负的质量密度"这奇异概念
音响-结构相互作用的前沿研究中注目的是"音响超材料"。在特定频率范围内等效呈现"负质量密度"或"负弹性率"的结构体，用完全不同的原理控制声音屏蔽和折射。比如10 cm见方的板状超材料性能超过面积相同的常规防音墙的情况也有报告。"负密度"这个概念有悖常理，但设计成局部共振结构（弹簧-质量系统的集合体）就能实现。作为新一代静音设计的变革技术，研究加速。
音响-结构相互作用的故障排除
一个个来看。
1. 耦合界面法向不一致
症状: 音压完全为零或期望值的数倍。
原因: 结构面的法向没有指向音响腔体内侧。Nastran中法向必须向外（音响侧），否则connectivity不对。
对策: 用前处理器（HyperMesh、ANSA等）可视化法向并统一。Nastran的BGPDT卡可检查。
这种花好多小时都找不到原因的问题呢。
2. 模态漏失
症状: 特定频段实测和分析差距很大。
原因: 模态频率响应（SOL 111等）用的模态数不足。必须包括到分析频率上限1.5倍为止的音响和结构的所有模态。
对策: 使用残差向量（residual vector）。Nastran的RESTVEC参数、Abaqus的RESIDUAL MODES。
3. 高频精度恶化
症状: 分析频率上限附近FRF出现波动，可信度低。
原因: 音响网格波长分解度不足。没有满足"6单元/波长"规则。
对策: 重新检讨网格密度或缩小目标频率范围。用二阶单元的话相同节点数下分解度更高。
4. 阻尼参数不匹配
结构阻尼和音响吸声是分别设定的吧？
症状: 共振峰比实测尖锐或宽阔。
原因: 结构损失系数 $\eta_s$ 与音响壁面阻抗与实际不符。
对策: 结构阻尼从实验模态同定（半值宽度法）。音响吸声率用残响室法或阻抗管的实测值。
| 参数 | 典型值 | 备注 |
|---|---|---|
| 钢板损失系数 | 0.001~0.003 | 无制振 |
| 有制振材钢板 | 0.01~0.05 | CLD处理后 |
| 车室吸声率（低频） | 0.05~0.15 | 含座椅、地毯 |
| 车室吸声率（高频） | 0.2~0.5 | 同上 |
5. Nastran典型错误
耦合分析特有的陷阱都明白了。实务中有错误就回头看这个清单。
"分析收敛了却与实测差10分贝"——音响耦合常见陷阱
音响-结构相互作用分析中，"计算完成了但实验值差10 dB以上"这样的困惑在现场很常见。最常见的原因是"边界条件设置失误"。实际零件是用螺栓固定的，但分析用完全固定，所以实际紧固刚度与分析不同导致固有值偏，共鸣峰频率偏。再加上内部损失系数（阻尼）设定用教科书值而不是实测值的话，峰值附近会产生很大误差。通常以为"模型不对＝网格粗"，但多数情况改善边界条件和阻尼值就解决了。
有帮助
更详细
错误