材料力学数值方法：边界元法(BEM)：BEM的离散化技术.pdfVIP

材料力学数值方法：边界元法(BEM)：BEM的离散化技术

1引言

1.1边界元法(BEM)概述

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

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

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

在于它只需要对问题的边界进行离散化，而不是整个域，这在处理无限域或半

无限域问题时尤其有效，可以显著减少计算资源的需求。

BEM的基本思想是将偏微分方程转化为边界积分方程，然后对边界上的未

知量进行数值求解。这一转化过程通常涉及到格林函数（Green’sfunction）的

使用，格林函数描述了在边界上施加单位力或单位热流时，整个域内的响应。

通过格林函数，可以将域内的问题转化为边界上的问题，从而简化计算。

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

边界元法（BEM）与有限元法（FEM）在处理材料力学问题时，各有其特

点和适用范围。FEM是一种基于域内离散化的数值方法，它将整个问题域划分

为许多小的单元，然后在每个单元内近似解的分布，通过单元间的连接来构建

整个域的解。相比之下，BEM仅对问题的边界进行离散化，这在处理无限域或

边界条件复杂的问题时，可以大大减少计算量和模型的复杂度。

1.2.1代码示例：使用BEM求解二维弹性力学问题

下面是一个使用Python和scipy库来实现BEM求解二维弹性力学问题的简

化示例。在这个例子中，我们将求解一个受集中力作用的无限域中的半圆形边

界上的位移。

importnumpyasnp

fromscipy.integrateimportquad

fromscipy.specialimportjv

#定义格林函数

defgreen_function(r,theta,r0,theta0):

#r和r0是从原点到点和源点的距离

#theta和theta0是点和源点的极角

#这里使用的是二维弹性力学的格林函数

G=1/(2*np.pi)*np.log(r/r0)*np.cos(theta-theta0)

returnG

1

#定义边界上的点

boundary_points=np.array([(1,0),(np.cos(np.pi/4),np.sin(np.pi/4)),(0,1)])

#定义集中力的位置和大小

force_location=(0,0)

force_magnitude=1

#计算边界上的位移

defcalculate_displacement(point):

r=np.sqrt(point[0]**2+point[1]**2)

theta=np.arctan2(point[1],point[0])

displacement=0

forbpinboundary_points:

r0=np.sqrt((bp[0]-force_location[0])**2+(bp[1]-force_location[1])**2)

theta0=np.arctan2(bp[1]-force_location[1],bp[0]-force_location[0])

displacement+=force_magnitude*quad(lambdat:green_function(r,theta,r0,t),0,2*np.pi)

[0]

returndisplacement

#输出边界上的位移

displacements=[calculate_displacement(bp)forbpinboundary_points]

print(边界上的位移:,displacements)

1.2.2解释

在这个示例中，我们首先定义了一个格林函数green_function，它描述了在

二维弹性力学中，任意点与源点之间的相互作用。然后，我们定义了边界上的

几个点boundary_points，以及集中力的作用位置和大小force_location和

force_magnitude。

calculate_displacement函数用于计算边界上每个点的位移。它通过积