网格点数少看起来万能，但肯定也有缺点吧？

优点：仅边界离散化、无限区域的自然处理、高精度缺点：非线性问题难以适用、生成密矩阵（计算复杂度 $O(N^2)$）高速BEM：通过FMM（高速多极展开法）可降至 $O(N \log N)$

密矩阵听起来很麻烦。但如果要做声学或电磁波分析，BEM似乎是绕不过去的。

确实如此。在实际工作中也会使用FEM-BEM耦合法。根据问题的特性选择最适合的方法，这是实务中的智慧。

BEM — CAE用语解释

Q: 网格点数少看起来万能，但肯定也有缺点吧？

优点：仅边界离散化、无限区域的自然处理、高精度 缺点：非线性问题难以适用、生成密矩阵（计算复杂度 $O(N^2)$） 高速BEM：通过FMM（高速多极展开法）可降至 $O(N \log N)$

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

CAE visualization for bem - technical simulation diagram

BEM（边界元法 / Boundary Element Method）

🧑‍🎓

老师，我经常听说FEM，但"BEM"是什么？和网格划分的位置有什么不同吗？

BEM的理论基础

BEM的基本概念

🧑‍🎓

边界元法（BEM）和FEM在根本上有什么区别？说"只离散边界"，但内部的应力怎么求呢？

🎓

根本区别在于使用满足支配方程的"基本解"。FEM用近似函数填充整个区域，而BEM则将积分方程转化，只求解边界上的未知量。例如，对于三维拉普拉斯问题，支配方程

$$ \nabla^2 u = 0 $$

转化为边界Γ上的积分方程

$$ c(\mathbf{x}) u(\mathbf{x}) = \int_\Gamma \left[ G(\mathbf{x}, \mathbf{y}) \frac{\partial u}{\partial n} - u(\mathbf{y}) \frac{\partial G}{\partial n} \right] d\Gamma(\mathbf{y}) $$

其中

$$ G(\mathbf{x}, \mathbf{y}) = \frac{1}{4\pi |\mathbf{x} - \mathbf{y}|} $$

是基本解。求出边界上的u和∂u/∂n后，将内部点x代入上式，可直接计算区域内任意点的u。应力也是如此。

🧑‍🎓

基本解总是能用吗？对于非均质材料或非线性问题，能适用吗？

🎓

很好的问题。基本解只能对线性"常系数"微分方程进行解析求导。即，均质线性弹性体、稳态热传导、声学等是典型情况。当材料性质随位置变化（非均质）或出现塑性变形等材料非线性时，标准BEM难以适用。这种情况需要用领域分割法（Domain Decomposition）将问题分解为均质区域，或使用DR-BEM（双倒数边界元法）等扩展技术引入体积积分项。

🧑‍🎓

听说BEM能处理无限远处的条件。具体怎么建模呢？FEM的话需要把网格延伸到远处吧。

🎓

确实，这是BEM的重大优势之一。基本解本身具有在无限远处收敛到零的性质，因此在建立边界积分方程时，无需明确考虑无限远处的边界条件。例如，用CFD模拟飞机的风洞试验时，FEM/CFD需要在计算域外侧建立"远场边界"网格，并设定一致流等边界条件。而BEM（或面板法）只需对机身表面进行面板离散，通过基本解的叠加自动满足无限远处的一致流条件。这样可大幅简化网格生成。

BEM的数值计算手法

离散化与特异积分

🧑‍🎓

听说边界要分割成单元。和FEM的形状函数不一样吗？还有，对积分方程进行数值积分时，基本解在特异点（x=y）处发散怎么办？

🎓

单元和形状函数的思想和FEM相似，但未知量不同。BEM是将边界分割成"面板"或"单元"，在每个单元上用节点值和形状函数（线性或二次）对物理量（势函数u和流量∂u/∂n）进行近似。特异积分是个难题。被积函数中的

$$ G(\mathbf{x}, \mathbf{y}) = \frac{1}{4\pi r} $$

在r→0时会像1/r那样发散。为了处理这个问题，在特异点周围通过极坐标变换进行解析积分，或用特异性消除变量变换（如Telles变换），然后进行标准高斯积分。实现时，根据单元是否包含特异点来切换积分例程。

🧑‍🎓

离散化后最终会得到线性方程

$$ [H]{u} = [G]{q} $$

吧。这个系数矩阵[H]和[G]有什么性质？和FEM的刚度矩阵比较怎样？

🎓

FEM的刚度矩阵是稀疏的（带状矩阵）且对称正定，而BEM的系数矩阵是"密"的，一般来说是非对称的。这是因为基本解

$$ G(\mathbf{x}_i, \mathbf{y}_j) $$

即使离散节点i和j距离较远也不为零，表现所有节点间的相互作用。但在某些问题中也可以对矩阵进行对称化处理。这个密矩阵的生成需要计算成本，但由于未知量只有边界上的节点数，整体规模往往较小。求解联立一次方程采用直接法LU分解，或对大规模问题采用高速多极法（FMM）。

BEM的实务应用

适用的问题与前处理

🧑‍🎓

请告诉我在实际设计工作中什么时候应该用BEM。还有，模型制作时需要特别注意什么？

🎓

BEM的优势发挥在"无限或半无限区域""薄壳结构""应力集中部的详细分析"。具体例子：

1. **声学分析**：汽车室内噪声、扬声器辐射声场。可自然处理无限远处的辐射条件。 2. **电磁场分析**：天线辐射模式、EMC（电磁兼容性）问题。 3. **断裂力学**：裂纹尖端应力强度因子的计算。只需离散裂纹面，能精确捕捉特异场。 4. **势流**：船舶和飞机周围的流场（面板法）。模型制作时，需准备**仅包含边界的CAD数据（表面模型）**，然后进行网格化。角部和曲率大的部分需细化网格。另外，不同材料接触的"界面"存在于两侧，会增加节点数，需注意。

🧑‍🎓

边界条件的设置中，和FEM相比有容易出错的地方吗？

🎓

是的。FEM直接在节点或单元面上指定迪利克雷条件（固定）和诺伊曼条件（荷载），而BEM的边界条件与"基本解的类型"密切相关。例如，在稳态热传导中既有"温度u已知"的迪利克雷边界，又有"热流量q已知"的诺伊曼边界混存时，需要将离散化方程

$$ [H]{u} = [G]{q} $$

中的{u}和{q}的已知部分和未知部分重新整理到方程左右两边，这称为"系统矩阵的组装"。索引错误会导致完全错误的结果。商用软件会自动处理，但自编程序需特别小心。另外，无限远处的条件是基本解隐含的，无需明确设置。

BEM软件比较

主要BEM软件与其特征

🧑‍🎓

请告诉我市售的CAE软件中搭载BEM的有哪些。Ansys和Abaqus也有吗？

🎓

Ansys和Abaqus等通用结构/流体软件的主求解器是FEM/FVM。但针对特定物理领域有BEM模块的软件，或专门的BEM软件存在。

* **COMSOL Multiphysics**："AC/DC模块"的静电、静磁场分析一部分可进行BEM与FEM的耦合。"声学模块"中也可用BEM替代无限单元进行音场分析。 * **Altair HyperWorks**："HyperSizer"等工具组合中某些可能包含BEM功能，但不是主流。 * **BEM专用软件**: * **BEASY**：腐蚀、防腐设计、断裂力学专用的老牌BEM软件。 * **Actran**：在声学振动分析中很强，支持FEM与BEM混合分析。在汽车航空工业中应用广泛。 * **Wavenology**：专门用于电磁场和声学分析的BEM求解器。与"万能"的FEM软件不同，BEM软件往往专注于"声学""电磁场""断裂"等特定专业领域。

🧑‍🎓

听说可以把BEM和FEM连成一起用"混合分析"。什么时候、怎么样用？

🎓

这是很实用的方法。比如"汽车室内噪声分析"。发动机和路面的振动传通过复杂的结构。这个结构振动传播可以用FEM很好地计算。但从振动的车身辐射出的"声"传向无限空间，BEM最适合。这种情况下的工作流是：

1. 用FEM（如Abaqus/Standard）进行车身结构振动分析，得到振动速度。 2. 将振动速度数据作为边界条件输入BEM软件（如Actran）。 3. 用BEM离散车体外表面，计算辐射声场。这样，FEM负责"区域内的复杂现象"，BEM负责"无限区域的辐射"，各司其职发挥所长。这是混合分析的理念。数据交换需进行节点信息和物理量的映射，商用软件间通常提供专用接口。

BEM故障排除

常见的数值错误与对策

🧑‍🎓

听说BEM的计算结果有时在角部或尖端处出现振荡或明显异常（特别是超大值）。原因和对策是什么？

🎓

这是典型问题，主要有两个原因。

1. **几何特异点**：CAD模型的锐角（棱）或尖端。在这些点法线方向不唯一，边界积分方程的系数

$$ c(\mathbf{x}) $$

确定变得不稳定。对策是在CAD上给角加微小倒角（圆角），或对包含特异点的单元采用特殊处理（双节点法等）进行离散化。 2. **网格不当**：单元形状极度细长（纵横比很大）或相邻单元尺寸剧烈变化时，数值积分精度恶化，结果出现振荡。BEM中网格质量检查比FEM更重要，单元纵横比应保持在10以下，相邻单元尺寸比应在2倍以内。

🧑‍🎓

计算收敛了，但守恒律（比如闭合曲面的热流量总和为零）没有满足。应该查什么？

🎓

首先怀疑数值积分，特别是"特异积分"的精度不足。前面讲过的特异点处理不完全，周围的积分误差会集中。对策如下：

* **增加积分点数**：靠近特异点的单元中，不要用标准高斯积分点（如4点），而要增到8点或16点。 * **重新审视单元划分**：物理量急剧变化的区域（集中荷载作用点附近等）要细化网格。 * **后处理结果**：商用软件有时提供"结果平滑化"选项，用更高精度的积分公式重新计算守恒量。另外，在整个闭合曲面上检验守恒律是验证结果可信度的有效方法。允许误差因问题而异，但热流量总和在内部发热的1%以内基本满意。

🧑‍🎓

大规模模型中系数矩阵的生成和联立一次方程的求解花时间太长。有什么高速化方法吗？

🎓

BEM的瓶颈正是那个"密矩阵"。对策如下：

* **采用高速算法**：**高速多极法（FMM: Fast Multipole Method）** 是决定性技术。通过近似处理遥远单元间的相互作用，将计算复杂度从O(N^2)降至O(N log N)或O(N)。大规模声学分析软件多数都实现了这个算法。 * **并行计算（HPC）**：矩阵生成和求解在多核/GPU上进行分散处理。FMM本质上适合并行化。 * **聚类处理**：物理上遥远的单元组相互作用可用低秩近似（如ACA：自适应叉近似）表示。选择商用软件时，要检查是否支持这些高速化技术（如Actran的FMM求解器或GPU支持）。单元数超过1万时，没有这些技术是无法得到实用计算时间的。