边界元法（BEM）声学分析
边界元法（BEM）声学理论基础
声学BEM是什么
老师，声学BEM是什么？
将Helmholtz方程转换为边界积分方程，仅用表面网格计算声场的方法。不需要体积网格，所以能够自然处理外部声场（无限域）是最大的优点。
支配方程
声学场的Helmholtz方程：
$k = \omega/c$：波数。用Green定理将其转换为边界积分形式：
$G$：自由空间Green函数：
$c(\mathbf{x})$：边界上为$1/2$，域内为$1$，域外为$0$。
与FEM相比有什么优势？
总结
BEM起源于1960年代的弹性理论
边界元法的数学基础由斯坦福大学的Jaswon（1963年）和Sympson奠定。他们发表了弹性问题的积分方程离散化方法，但当时并未考虑声学应用。声学BEM的应用转用始于十年后，Schenck在1968年作为IBM技术报告发表的CHEFS方法成为先驱，随后影响了NASTRAN和商用求解器。
边界元法（BEM）声学数值计算手法
声学BEM的离散化
怎样在数值上求解边界积分方程？
用三角形或四边形单元离散表面。声压$p$和法向速度$v_n$用节点值近似。
$[H]$, $[G]$：影响矩阵。各单元间Green函数的积分值。
奇异积分的处理
BEM实现中技术上最困难的是奇异积分（当$r \to 0$时Green函数发散）的处理：
- 弱奇异性（$1/r$）：用极坐标变换正则化
- 强奇异性（$1/r^2$）：用Cauchy主值定义
- 超奇异积分：Hadamard有限部分积分
非唯一性问题
在外部BEM中，闭合边界的内部固有振动频率会导致解不唯一。对策：
- CHIEF法：在内部点配置附加方程（超定系统）
- Burton-Miller法：将通常的BIE与法向导数BIE线性组合。在理论上完全
Burton-Miller法是标准方法。需要处理超奇异积分，但能确保去除非唯一性。
FMM（快速多极法）
你提到密矩阵限制了大规模问题的求解…
FMM（Fast Multipole Method）可以解决这个问题。将远方的单元群组合作近似计算：
- 内存：$O(N^2) \to O(N)$
- 计算量：$O(N^2) \to O(N\log N)$
百万级单元规模的BEM变得实用。Actran、COMSOL等已实装FMM。
总结
FMM将BEM计算量从O(N²)降至O(N log N)
传统BEM对节点数N需要$O(N^2)$的内存和计算量，大规模模型实际上不可行。1987年由Greengard和Rokhlin发表的快速多极法（Fast Multipole Method）带来了革命，使汽车车身规模（数十万节点）的声学分析变得现实。NuancesVirtualLab Acoustics和LMS Sysnoise 5.x等工具都是较早实现FMM-BEM的产品。
边界元法（BEM）声学实务应用
声学BEM实务
发动机辐射音、轮胎噪声、变压器噪声、排气管的声学辐射是典型应用。
分析流程
1. 结构振动分析 — 用FEM获得表面振动速度$v_n$
2. 创建BEM表面网格 — 从FEM网格提取或独立创建
3. 设置边界条件 — 设置$v_n$（Neumann条件）
4. BEM求解 — 计算表面声压$p$
5. 声场评估 — 计算任意观测点的声压、辐射功率
实务检查清单
BEM网格指南
| 最高频率 [Hz] | 音速340m/s下的波长 [m] | 单元尺寸上限 [mm] |
|---|---|---|
| 500 | 0.68 | 113 |
| 1000 | 0.34 | 57 |
| 2000 | 0.17 | 28 |
| 5000 | 0.068 | 11 |
| 10000 | 0.034 | 5.7 |
5kHz以上网格需要非常细密。
是的。高频下BEM的网格成本急剧增加，应考虑切换到SEA。
车室内声学需要合理选择BEM和FEM
汽车车室（舱室）噪声分析中，1kHz以下的低频用FEM有利，1kHz以上的高频用BEM有利。Audi AG在2000年代初发表的基准研究表明，BEM-FEM耦合在1500Hz频带达到了实验值3dB以内的精度。但实务中吸音材料复素阻抗值的设置是影响精度的最大课题。
边界元法（BEM）声学软件比较
工具
Actran看起来压倒性地处于优势…
汽车NVH领域Actran确实是事实标准。航空航天也广泛应用。COMSOL适合小规模研究和验证。Siemens系与Nastran用户兼容性高。
Sysnoise收购强化了Actran的声学BEM功能
声学BEM老牌求解器LMS Sysnoise（比利时LMS International制造）在2000年代确立了行业标准地位，但2012年被Siemens收购并整合到Simcenter Testlab品牌。与此同时，FFT公司（现Siemens旗下）的Actran通过声学有限元和边界元的混合方式进行差异化，在航空发动机短舱透射损失分析方面有优势。OptiStruct和NASTRAN也在后来增加了BEM模块。
边界元法（BEM）声学先端研究
先端研究
IGA-BEM太厉害了。完全不需要网格吗？
直接使用CAD的NURBS表面，所以不存在"网格生成"这个工序。对曲面多的形状（汽车车身、飞机机身）特别有威力。但仍处于研究阶段，商用实装有限。
水下声学BEM在潜艇设计中曾是机密
水下辐射音的边界元分析在1970-80年代的美国海军DARPA计划中快速发展。潜艇螺旋桨辐射音低减是军事机密，技术转让到民间延迟到了冷战结束后的1990年代。现在ABAQUS/Acoustics和Actran等标准提供水下声学BEM，也被用于评估海上风电基础的水下噪声。
边界元法（BEM）声学故障排除
声学BEM故障
| 症状 | 原因 | 对策 |
|---|---|---|
| 特定频率声压发散 | 非唯一性（内部共鸣） | 激活Burton-Miller法 |
| 声压结果全面反转 | 法向方向指向内侧 | 修正所有单元的法向向外 |
| 结果与振动分析不符 | FEM→BEM速度映射错误 | 确认网格对应关系。验证插值精度 |
| 计算缓慢/内存不足 | 密矩阵直接求解 | 激活FMM。检查网格数量 |
| 近场声压不精确 | 观测点在单元面上或极近 | 将观测点移至单元尺寸的2-3倍外 |
| 高频结果不稳定 | 网格过粗（<6单元/波长） | 细分网格。或切换到SEA |
法向方向搞反了影响这么大？
法向反转会改变边界积分方程的符号。声压正负反向，Burton-Miller法也无法正常工作。BEM模型最初必须检查法向。
没有CHIEF法BEM在固有频率会破坏
声学BEM存在"不正则频率（irregular frequency）"问题。在结构的特定固有频率处，BEM方程变为奇异矩阵，解发散。Schenck（1968年）提出的CHIEF（Combined Helmholtz Integral Equation Formulation）法通过在内部点配置附加方程来回避此问题，现在成为许多商用求解器的标准方法。但CHIEF点位置不当又会再次产生不稳定，需要注意。
更详细
错误