结构力学数值方法：边界元法(BEM)：边界元法概论与历史发展.pdf

结构力学数值方法：边界元法(BEM)：边界元法概论与历史

发展

1边界元法概论

1.1边界元法的基本概念

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

于解决偏微分方程问题，特别是那些在工程和科学领域中常见的问题。与有限

元法（FEM）不同，BEM主要关注问题的边界条件，将整个问题域的积分转化

为边界上的积分，从而减少了问题的维数，使得计算更加高效。

在BEM中，问题的解通过边界上的未知量（如应力或位移）来表示，这些

未知量通过边界积分方程（BoundaryIntegralEquation,BIE）来求解。BIE是通过

格林函数（Green’sfunction）和问题的边界条件构建的，格林函数描述了在边

界上施加单位力或单位位移时，整个域内的响应。

1.1.1示例

假设我们有一个二维弹性问题，需要求解一个圆盘在边界上受到均匀压力

的情况。我们可以使用BEM来求解这个问题，首先定义边界上的未知量，然后

通过边界积分方程来求解这些未知量。

importnumpyasnp

fromscipy.integrateimportquad

#定义格林函数

defgreen_function(r,s):

return1/(2*np.pi*r)*np.log(np.abs(r-s))

#定义边界上的积分

defboundary_integral(r,boundary,pressure):

integral=0

forsinboundary:

integral+=quad(lambdat:green_function(r,s)*pressure(s)*t,0,1)[0]

returnintegral

#定义边界

boundary=np.linspace(0,2*np.pi,100)

#定义边界上的压力

pressure=lambdatheta:np.cos(theta)

1

#求解边界上的未知量

unknowns=[boundary_integral(r,boundary,pressure)forrinboundary]

这个例子中，我们定义了一个简单的格林函数和边界积分方程，然后使用

scipy.integrate.quad函数来数值求解积分。boundary和pressure分别定义了边

界和边界上的压力分布，最后通过循环求解边界上的未知量。

1.2边界元法与有限元法的比较

边界元法与有限元法在处理问题时有显著的不同。FEM将整个问题域离散

化，通过在每个节点上求解局部方程来得到全局解。而BEM只关注边界，将问

题域的积分转化为边界上的积分，这样可以减少计算量和存储需求，特别是在

处理无限域或半无限域问题时。

1.2.1优势

减少计算量：BEM只需要处理边界上的未知量，而FEM需要处理

整个域内的未知量。

无限域问题：BEM在处理无限域问题时不需要无限大的网格，而

FEM则需要通过特殊技术来模拟无限域。

高精度：BEM在处理某些类型的问题时，如弹性问题，可以提供

更高的精度。

1.2.2劣势

复杂边界条件：对于复杂的边界条件，BEM的实现可能比FEM更

复杂。

内部解：BEM在求解内部点的解时不如FEM直接和方便。

非线性问题：对于非线性问题，BEM的处理可能比FEM更复杂。

1.3边界元法的应用领域

边界元法在多个领域都有广泛的应用，包括但不限于：

结构力学：用于求解弹性、塑性、断裂力学等问题。

流体力学：用于求解流体动力学、声学等问题。

热传导：用于求解热传导、热辐射等问题。

电磁学：用于求解电磁场、电磁波散射等问题。

边界元法在处理无限域、半无限域、多连通域以及高维问题时特别有效，

这使得它在工程和科学计算中成为一个重要的工具。

以上内容详细介绍了边界元法的基本概念、与有限元法的比较以及其应