简介：本文系统阐述声学多物理场仿真的核心定义、技术原理、应用场景及实践要点。通过拆解仿真流程中的关键环节，帮助工程师与科研人员理解如何借助多物理场耦合分析解决实际声学问题，涵盖从基础建模到高级应用的完整知识框架，为声学设计与研究提供可落地的技术参考。
声学多物理场仿真是一种基于数值计算方法，在统一计算框架内同时求解声场与其他物理场（如结构振动、流体流动、电磁场、热场等）相互作用的分析技术。其核心目标不是孤立地计算声波传播，而是还原真实物理环境中声波与周围介质、结构及其他能量形式的耦合关系。
从技术视角看，该技术以有限元法、边界元法或有限体积法等离散化方法为基础，将连续的偏微分方程转化为可求解的代数方程组。用户通过定义几何模型、材料属性、边界条件与激励源，即可获得声压分布、声强矢量、振动位移、频率响应等关键结果。
从应用视角看，声学多物理场仿真解决的是“设计验证前移”的问题——在物理样机制作之前，通过数字样机预测产品声学性能，从而降低研发成本、缩短迭代周期。它特别适用于难以通过解析公式求解的复杂几何结构、非均匀介质以及强耦合问题。
传统声学设计高度依赖经验公式与样机测试。然而，现代工业产品日益复杂，例如扬声器单元与腔体结构耦合、汽车排气系统与气流相互作用、医疗超声探头与组织接触等场景，单纯依靠实验试错不仅成本高昂，而且难以获得内部场的全面数据。声学多物理场仿真的出现，正是为了应对以下核心挑战：
从行业需求来看，消费电子、汽车工业、建筑声学、医疗设备等领域对声学舒适度与功能性能的要求持续提升，相关技术人才需要掌握这一数字化设计工具，以应对从元件级到系统级的声学设计任务。
理解声学多物理场仿真，需要先认识其能力构成的五个层次：
仿真起始于几何准备。用户可以在仿真环境中直接构建简单结构（如管道、腔体），也可以从外部CAD系统导入复杂装配体。此阶段的关键问题包括：是否需要简化圆角、倒角等细节；如何将几何域划分为声学域（流体域）与结构域（固体域）；以及是否需要设置无限元域或完美匹配层（用于模拟开放空间）。
声学仿真需要精确的材料参数。对于流体介质（如空气、水），需指定密度、声速及动态粘度（考虑粘性损耗时）。对于多孔材料，则需要流阻率、孔隙率、曲折因子等参数。结构材料则需要弹性模量、泊松比、损耗因子等。多物理场耦合场景还需定义压电材料（用于换能器）或热粘性流体属性（用于微尺度声学）。
声学仿真涉及多种物理场接口，常见包括：
网格质量直接影响计算精度与资源消耗。声学问题的一般准则是每个波长至少包含5-10个网格单元，而热粘性边界层则需要更细的边界层网格。求解器选择同样关键：频域问题通常使用直接或迭代线性求解器；瞬态问题则需要考虑时间步长与稳定性条件（CFL条件）。
仿真输出包括声压级分布、声强流线、振动位移云图、模态振型、频谱响应曲线等。后处理不仅是云图展示，更重要的是提取设计所需的关键指标，如传输损耗（TL）、插入损耗（IL）、指向性图案、声功率、共振频率等。
声学多物理场仿真遵循“物理建模—数学离散—数值求解—结果验证”的基本工作流。以最常见的声-固耦合问题为例，其核心过程如下：
物理控制方程：
有限元离散：将计算域划分为有限数量的单元，在每个单元内用形函数近似未知场，将所有单元方程组集为整体矩阵方程。对于频域问题，最终求解形如\（(K - \omega^2 M + i \omega C) x = F\）的线性方程组，其中K为刚度矩阵、M为质量矩阵、C为阻尼矩阵，x为待求节点变量。
求解与验证：求解完成后，需要检查收敛性（如网格无关性验证）、能量守恒性（如总声功率是否一致），并将仿真结果与理论解或实验数据对比，确认建模假设的合理性。
声学多物理场仿真已在多个工业与科研领域形成成熟应用模式，以下按方向进行归纳：
扬声器、麦克风、耳机等设备中，声学仿真用于分析振膜振动模式与声场辐射特性。通过耦合结构力学与压力声学，设计者可优化音圈与振膜结构参数，改善频率响应曲线与指向性，同时分析腔体泄漏、倒相孔调谐等非线性效应。
在汽车、飞机、家电产品中，仿真帮助工程师定位主要噪声源（如发动机、轮胎、风噪），评估吸声材料布置方案，计算隔声量（STL）。通过参数化扫描，可快速获得不同阻尼材料厚度与覆盖面积对车内噪声的影响规律。
医用超声探头与工业无损检测设备广泛采用压电换能器。仿真同时求解压电方程、结构振动方程与声波传播方程，预测换能器的发射/接收灵敏度、声束聚焦特性以及在多层介质（如人体组织-骨骼界面）中的反射与透射信号，为探头阵列设计提供依据。
超材料依靠周期性微结构实现自然材料不具备的声学特性，如负等效密度、负等效模量。仿真在此领域主要进行两类任务：一是对单个单元进行特征频率与参数反演分析，提取等效参数；二是对有限尺寸超材料板进行全波仿真，验证其隔声、聚焦或隐身效果。
在微流控芯片中，利用驻波声场产生的声辐射力实现对微粒、细胞、液滴的无接触操控。仿真需同时求解声场与流场（或考虑粘性损耗），计算声辐射力分布与粒子运动轨迹。此类仿真要求极高网格精度，常需要使用热粘性声学接口，以捕捉边界层内的声能耗散。
虽然建筑声学通常采用射线追踪法或统计能量法，但在分析低频模态分布、墙体隔声性能及耦合空间声场时，也需要有限元/边界元声学仿真。典型应用包括：会议厅的混响时间优化，音乐厅的早期反射声分布分析等。
声学多物理场仿真与几种常见声学分析技术之间存在明确边界，理解这些区别有助于选型：
| 分析方法 | 适用频率 | 核心特点 | 典型场景 |
|---|---|---|---|
| 多物理场有限元仿真 | 低频-中频 | 求解精确，考虑结构与声场双向耦合，可处理复杂几何与材料非均匀性 | 扬声器单元设计、换能器优化、小尺度隔声结构 |
| 射线声学法（几何声学） | 高频 | 忽略波动性，以声线追踪能量传播路径，计算速度快 | 大型厅堂混响、厂房噪声分布 |
| 统计能量分析（SEA） | 高频 | 将系统划分为多个子系统，用能量流方程描述子系统间传递 | 航空航天舱室宽频噪声、船舶舱室噪声 |
| 边界元法（BEM） | 中频 | 仅需边界离散，自动满足无穷远辐射条件，适合开放空间外场问题 | 消声器插入损失、发动机外辐射噪声 |
多物理场仿真的独特优势在于处理“耦合”与“细节”——当振动与声场相互作用不可忽略、几何精细度影响性能、或需要观察内部场分布时，它是唯一可行的数值工具。而射线法与SEA在极高频宽频问题中计算效率更高，但无法提供场分布细节。
开展声学多物理场仿真时，以下几方面问题需要重点关注：
求解最高频率决定了最大网格尺寸。例如，在空气中模拟20kHz声波（波长约17mm），要求最小单元不超过3.4mm；若同时考虑1mm厚壳体的弯曲波，则壳体厚度方向可能需要至少2层单元。高频+大几何模型将产生巨大计算量，此时需要考虑对称性简化、二维近似或并行计算。
选用完美匹配层（PML）模拟无反射边界时，PML厚度需足够且与内域网格匹配，否则会产生虚假反射。对于刚性壁面与阻抗边界，需要确认材料阻抗频率相关性是否被计入。在瞬态分析中，需特别关注初值条件（如初始声压与速度）的设置。
并非所有场景都需要完整双向耦合。例如，研究轻质面板背后空腔的声吸收时，面板振动可能对声场影响微弱，可近似为刚性边界；而设计压电驱动器激励板结构辐射噪声时，必须双向耦合才能捕捉机电-结构-声学全过程。合理简化模型可显著降低计算成本。
声速随温度与湿度变化，多孔材料参数也可能随压缩状态改变。在仿真前应有明确的参数获取渠道（实测或文献数据），并评估参数不确定性对结果的影响，必要时进行参数扫描或敏感性分析。
仿真预测必须经过验证方可作为设计依据。最低限度应进行网格收敛性检查与能量平衡检查；有条件时，应使用激光测振仪、传声器阵列、阻抗管等设备获取实验数据，与仿真结果对比频率响应峰值位置、幅值水平及模态振型。
声学多物理场仿真的本质，是用统一的数值框架逼近真实世界中的声-固-流-电-磁耦合行为。它的价值不在于替代测试，而在于帮助研究者打开“设计黑箱”——在样机之前优化性能、在实验之前指导方案、在解析失效时提供精确解。该技术覆盖从微尺度声操控到大型消声结构的广阔尺度范围，是声学工程师数字技能库中不可或缺的一环。
对于计划进入该领域的工程师或研究者，建议从简单频域声学分析入手，逐步拓展到结构耦合与瞬态问题，并在实际项目中持续积累边界条件设置与网格划分经验。技术工具的最终目标，始终是服务于对声学现象更深刻的理解、对声学器件更高效的设计，以及对声学问题更可靠的决策。