噪声预测仿真 — FEM-BEM耦合到SEA

说白了就两步。首先用FEM做结构振动分析，求出表面速度分布。然后用BEM（边界元法）计算这个振动表面辐射的声压分布。

最典型的就是汽车路面噪声。轮胎→悬架→车体面板，振动一层层传上来，最后在车室内形成噪声。这种"结构→声学"的连锁反应就是FEM-BEM耦合的标准应用。还有压缩机噪音、家电马达声、建筑机械的环境噪声，都是这个套路。

对头。低频（～500 Hz）就用FEM-BEM，中高频（500 Hz～）就得换SEA（统计能量法），空气动力噪声有FW-H法。不同频率不同法宝，这就是声学工程的诀窍。把这个全局图景先装在脑子里，后面就好理解了。

起点是波动方程。音压 $p(\mathbf{x}, t)$ 满足：

其中 $c$ 是声速（空气中约343 m/s）。把它变换到频域，角频率为 $\omega$ 的定常声场就满足Helmholtz方程：

完全对。经验法则是1波长至少6个单元（二阶单元可以是3个）。例如1000 Hz时 $\lambda \approx 0.34$ m，所以单元尺寸要在57 mm以下。频率越高网格数量暴增，这就是声学分析的"宿命"。

无限刚壁中的振动面（baffle面）辐射可以用Rayleigh积分直接求解。面的法向速度 $v_n(\mathbf{y})$ 已知，任意点 $\mathbf{x}$ 的音压就是：

这里 $r = |\mathbf{x} - \mathbf{y}|$ 是音源点到接收点的距离，$\rho_0$ 是空气密度。$e^{-jkr}/r$ 就是格林函数，代表点源的球面波。扬声器辐射模式计算用这个就够了。

这就要用Kirchhoff-Helmholtz积分方程和BEM了。那才是通用方法，处理任意闭曲面的辐射。

对Helmholtz方程用格林定理，外部区域（辐射问题）的Kirchhoff-Helmholtz积分方程是：

$G(\mathbf{x},\mathbf{y}) = \dfrac{e^{-jk|\mathbf{x}-\mathbf{y}|}}{4\pi|\mathbf{x}-\mathbf{y}|}$ 是3维自由空间的格林函数。$c(\mathbf{x})$ 在边界上取1/2，外部取1。由于 $\partial p/\partial n = j\omega\rho_0 v_n$，只要从结构FEM得到表面速度 $v_n$，代进去就能解出表面音压 $p$。然后用同样的积分再算外部任意点的音压。

完全同意。外部辐射问题用FEM的话，得把周围空气也网格化，还要设吸收边界（PML之类的）。BEM只需表面网格，远场辐射条件自动满足。汽车发动机这样的外部辐射问题上BEM优势太明显了。

主要两个指标。首先是声压级（SPL），接收点的音压与基准值的比较：

然后是声功率级（$L_W$），代表声源本身放出的能量，与接收位置无关：

声功率 $W$ 是把振动面的声强积分：

就是这个意思。标准里ISO 3744和ISO 3745规定了声功率的测量方法。仿真结果跟测试对标的时候就按这套框架来。

原理是一样的，但噪声预测有独特的重点。目的是求出表面法向速度 $v_n(\omega)$。结构动力学方程是：

位移 $\mathbf{u}$ 算出来后，表面速度就是 $v_n = j\omega\, \mathbf{u} \cdot \mathbf{n}$。关键是阻尼模型的选择。Rayleigh阻尼（$\mathbf{C} = \alpha\mathbf{M} + \beta\mathbf{K}$）方便但广带宽内阻尼比会飘。汽车NVH分析通常用实测模态阻尼来逐频段设置才可靠。

完全对。阻尼比真实值的2倍就能让共振峰的SPL变6 dB。所以结构FEM要跟试验模态对标，MAC值要≥0.8，这是声学预测信度的根基。

把表面分成 $N$ 个单元离散化，积分方程变成矩阵形式：

$\mathbf{p}$ 是节点音压向量，$\mathbf{q} = \partial p/\partial n = j\omega\rho_0 v_n$ 是法向压力梯度。结构FEM已经给出 $v_n$，所以 $\mathbf{q}$ 是已知的。解这个方程得到表面音压 $\mathbf{p}$，再用积分就能算出外部任意点的音压。

BEM矩阵是稠密矩阵（满矩阵）。所有节点之间都有相互作用，所以是 $N \times N$ 的满矩阵。内存是 $O(N^2)$，直接法求解是 $O(N^3)$。这就是大规模问题的BEM瓶颈，后面的FMM-BEM就是来破解这个的。

单向耦合（弱耦合）最简单，实务中大部分用的就这个。结构FEM算表面速度，直接拿去给BEM当边界条件。忽略了音压反作用力。

双向耦合（强耦合）是当音响负荷显著影响结构时需要。比如薄钢板在水里振动。这时耦合方程是：

$\mathbf{L}$ 是音压→结构力的耦合矩阵。空气中的辐射问题几乎都用单向就够了。空气密度比结构小3个数量级，音响负荷的影响可以忽视。水下声学或密闭车室才需要双向耦合。

FMM（快速多极子法）是革命性的。把远场单元的贡献分组计算，矩阵-向量乘法从 $O(N^2)$ 降到 $O(N\log N)$。内存也变成 $O(N)$。这样百万级节点的问题就解得了，而传统BEM不可能。

完全是。Siemens Simcenter 3D、FFT ACTRAN、ESI VA One这些商用工具都内置FMM-BEM。OpenBEM、FastBEM这样的开源也有，但大型产业问题上商业工具的稳定性和前后处理还是领先。

SEA（统计能量分析）登场。把结构和声场分割成"子系统"，用能量平衡方程描述各子系统间的能量流：

子系统间的功率流用耦合损失因子 $\eta_{ij}$ 表示：$P_{ij} = \omega\eta_{ij}(E_i/n_i - E_j/n_j)$。$n_i$ 是模态密度，$E_i$ 是总能。不追踪单个模式，而是统计平均能量，所以几千Hz高频也只需几十个子系统。

SEA出不来空间分布，只有平均能量。"车室内这个座位比那个座位吵"这样的问题回答不了。只能说"这个面板向这个空间流多少能量"这种趋势。精度就看耦合损失因子的准度，实测同定最靠谱。

空气动力噪声用Ffowcs Williams-Hawkings（FW-H）方程。非定常CFD（通常LES或DES）算出流体压力波动，用声学比拟法推远场音压。

第一项是单极子（厚度噪声），第二项是偶极子（负荷噪声）。"ret"表示延迟时间求值。汽车侧镜风切音、飞机起落架噪声、风力机叶片噪声都是这个方法的经典应用。Ansys Fluent和OpenFOAM都集成了FW-H求解器。

就是。LES足够时间窗的非定常数据才能保证统计频谱稳定。汽车侧镜空力噪声分析，数百核×数日的计算常见。但FW-H的好处是CFD网格和声学网格分离，接收点位置改变不用重新算CFD，直接改变远场评估位置就行。

FEM-BEM耦合外部辐射问题的标准流程：

车体各面板（顶棚、地板、车门等）对接收点SPL的贡献分解。"哪个面板处理最有效"来优化对策成本。BEM积分按面板分组算。

铁律是"6单元/波长"。最大分析频率 $f_{\max}$ 下：

实務"FEM-BEM2000 Hz、以上SEA"使分多。FEM-SEA最近的。

超关键。车室内噪声预测中，吸音材阻抗错了就能让中高频SPL差5～10 dB。多孔材料用Biot理论或Johnson-Champoux-Allard模型，需要流阻、孔隙率、迂曲度等参数。有的厂商提供材料库，最好自己测。

外部BEM用的Kirchhoff-Helmholtz方程在物体内部的固有周波数处矩阵奇异。物理上跟外部音场无关但数值上解被污染。Burton-Miller定式（法向微分方程和原方程组合）完全消除。商用BEM通常默认开启，不用操心。

噪声预测（NVH: Noise, Vibration, Harshness）需要专门声学求解器。主要工具：

Simcenter 3D（）、Actran（MSC傘下）、VA One（ESI Japan）、COMSOL（直販）日本語技術体制。費「困電話相談」实务死活的重要。

答案就是混合FEM-SEA。长波长振动用FEM决定论处理，短波长用SEA概率处理。同一模型里"低频FEM精度"+"高频SEA效率"共存。

核心思想是用"扩散场的互易性"结合确定FEM区间和统计SEA区间。Langley & Cotoni (2004) 的定式是基础，Simcenter 3D和VA One实装。刚强骨架（FEM）+柔软面板（SEA）混合的车身特别有用。

代理模型已接近实用。少数FEM-BEM结果训练NN再优化，比跑几百case计算经济。PINN还主要学术论文，但潜力巨大。"5年后CAE工程师必学"的预测也有。

现在噪声预测目标从"减dB"转向"设计音质"。电动车没内燃机遮蔽，马达高频纯音和逆变器开关噪声突出。心理声学指标登场：

CAE中，BEM或SEA算出接收点音压时间序列或频谱后，后处理软件转换成心理声学指标。Simcenter Testlab、HEAD acoustics ArtemiS、ANSYS Sound分析工具都有这功能。电动车开发中"行驶时锐度≤1.5 acum"这样的目标就开始出现了。

10 dB偏差不稀奇，原因三大类：

铁规则："结构差音响必差"。相反结构FEM与实测对标好的话，BEM SPL通常3～5 dB精度以内吻合。

典型的网格不足症状。6单元/波长不足时数值色散误差暴增，相位积累失控，SPL乱飘。对策：

外部BEM的非唯一性问题（虚拟固有频率）。物体内部音响固有频率处矩阵条件数恶化，解失效。对策两个：

最新Simcenter 3DActranBurton-Miller標準装備。通常心配无用。古自作BEM確認要。"奇妙"出Burton-Miller有効化疑。

SEA准度关键由耦合损失因子（CLF）决定。常见问题：

黄金对策是功率注入法（PIM）实测同定CLF。逐个子系统输入振动，能量响应反演。费时但SEA结果信度大幅提升。

噪声预測结构FEM的出来全。相関始。结构合音響自然。、周波数帯方法切替判断力实务鍛、最初FEM-BEM的低周波問題経験積勧。分聞。