材料力学数值方法：边界元法(BEM)：BEM在断裂力学中的应用.pdfVIP

材料力学数值方法：边界元法(BEM)：BEM在断裂力学中的

应用

1材料力学数值方法：边界元法(BEM)在断裂力学中的应用

1.1引言

1.1.1边界元法(BEM)概述

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

于解决偏微分方程问题。与有限元法（FEM）相比，BEM主要关注于问题的边

界条件，将问题域的内部信息转化为边界上的积分方程，从而大大减少了计算

的自由度，提高了计算效率。在材料力学领域，BEM被广泛应用于弹性、塑性、

断裂力学等问题的求解。

1.1.2断裂力学基础理论

断裂力学是研究材料在裂纹存在下行为的学科，主要关注裂纹的扩展条件

和裂纹尖端的应力强度因子（StressIntensityFactor,SIF）。SIF是衡量裂纹尖端应

力集中程度的重要参数，其值的大小决定了裂纹是否会发生扩展。在断裂力学

中，线弹性断裂力学（LinearElasticFractureMechanics,LEFM）是最基础的理论，

它假设材料在裂纹尖端附近处于线弹性状态，裂纹扩展路径是确定的。

1.2BEM在断裂力学中的应用

在断裂力学中，BEM的应用主要集中在计算裂纹尖端的应力强度因子上。

由于BEM仅需要在边界上进行离散，对于裂纹问题，这种方法可以精确地捕捉

裂纹尖端的细节，而无需在裂纹尖端附近进行密集的网格划分，这在FEM中是

必要的但计算成本较高。

1.2.1应用原理

BEM将弹性力学的基本方程转化为边界上的积分方程。对于二维弹性问题，

基本方程可以表示为：

+=0

,

其中，是应力张量，是体积力。在BEM中，通过格林函数（Green’s

function）将上述方程转化为边界上的积分方程：

1

​​∂

=+

∂

​

+

其中，是位移，是格林函数，是边界，是问题域，是边界上的外

法向量。

1.2.2应用内容

在断裂力学中，BEM被用于计算裂纹尖端的应力强度因子。具体步骤如下：

1.建立模型：定义裂纹的几何形状和边界条件。

2.离散边界：将裂纹边界和问题域的边界离散为一系列的单元。

3.求解积分方程：对于每个单元，求解上述的边界积分方程，得到

单元上的应力和位移。

4.计算SIF：利用特定的公式，从单元上的应力和位移计算裂纹尖端

的应力强度因子。

1.2.3示例

2

假设我们有一个二维的裂纹问题，裂纹位于无限大平面的中心，长度为。

我们使用BEM来计算裂纹尖端的应力强度因子。

1.2.3.1数据样例

裂纹的几何参数和材料属性如下：

裂纹长度2=1.0m