材料力学数值方法：边界元法(BEM)：BEM中的边界条件处理.pdfVIP

材料力学数值方法：边界元法(BEM)：BEM中的边界条件处

理

1材料力学数值方法：边界元法(BEM)：BEM中的边界条件

处理

1.1引言

1.1.1边界元法(BEM)简介

边界元法（BoundaryElementMethod,BEM）是一种数值方法，主要用于解

决偏微分方程问题，特别是在材料力学领域中，它被广泛应用于求解弹性力学、

热传导、流体力学等问题。与传统的有限元法（FEM）相比，BEM的主要优势

在于它将问题的求解域从整个区域缩减到边界上，从而大大减少了计算量和所

需的存储空间。BEM基于格林函数和边界积分方程，通过将边界条件转化为积

分方程，然后对边界进行离散化，最终通过数值方法求解。

1.1.2BEM与有限元法(FEM)的比较

边界元法（BEM）和有限元法（FEM）是两种常用的数值求解偏微分方程

的方法，它们在处理材料力学问题时各有优势和局限性。

计算域的处理：FEM需要对整个计算域进行离散化，而BEM仅需

对边界进行离散化，这使得BEM在处理无限域或半无限域问题时更为高

效。

存储需求：由于BEM仅处理边界，其矩阵规模通常远小于FEM，

因此存储需求较低。

计算效率：BEM在处理复杂边界条件时，如接触问题，可能比

FEM更有效率，因为它直接在边界上处理这些条件。

数值稳定性：BEM在处理某些问题时，如高频率声学问题，可能

比FEM更稳定，因为它避免了内部节点的离散化。

适用范围：FEM在处理非线性问题、多物理场耦合问题时更为灵

活，而BEM在处理线性问题时更为有效。

1.2BEM中的边界条件处理

在BEM中，边界条件的处理是关键步骤之一，它直接影响到求解的准确性

和效率。边界条件可以分为三种主要类型：Dirichlet边界条件、Neumann边界

条件和混合边界条件。

1

1.2.1Dirichlet边界条件

Dirichlet边界条件，也称为第一类边界条件，规定了边界上的位移或温度

等物理量的值。在BEM中，处理Dirichlet边界条件通常涉及到将边界条件直接

代入到边界积分方程中，从而形成一个混合的代数方程组。例如，在弹性力学

问题中，如果边界上的位移已知，那么在构建边界积分方程时，可以直接将这

些位移值代入，从而减少未知数的数量。

1.2.2Neumann边界条件

Neumann边界条件，也称为第二类边界条件，规定了边界上的应力或热流

等物理量的值。在BEM中，处理Neumann边界条件通常需要使用格林函数的

导数，将边界上的应力或热流转化为边界上的积分项。例如，在热传导问题中，

如果边界上的热流已知，那么可以通过格林函数的导数，将热流转化为边界上

的积分项，然后通过数值方法求解。

1.2.3混合边界条件

混合边界条件是Dirichlet和Neumann边界条件的组合，它在边界上的某些

部分规定了位移或温度，而在其他部分规定了应力或热流。在BEM中，处理混

合边界条件需要同时使用上述两种方法，形成一个包含Dirichlet和Neumann条

件的代数方程组，然后通过迭代或其他数值方法求解。

1.3示例：使用BEM求解弹性力学问题

假设我们有一个二维弹性力学问题，其中边界上既有位移边界条件

（Dirichlet），也有应力边界条件（Neumann）。我们将使用BEM来求解这个问

题。

1.3.1数据样例

考虑一个矩形区域，其边界条件如下：

左边界：=0（位移边界条件）

右边界：=100（应力边界条件）

上下边界：=0（位移边界条件）

其中，和分别表示x和y方向的位移，表示x方向的应力。

1