固体力学中的边界积分方程及其边界元法综述

计算固体力学

读书报告

固体力学中的边界积分方程及其边界元法

综述

Review of the Boundary Integral Equation and Boundary Element Method in Solid Mechanics

土木工程系 2014年03月17日

评语

目录

摘要 (2)

A BSTRACT (2)

一、引言 (3)

1)什么是边界元法[1] (3)

2)积分方程和边界元法的发展历史[2] (3)

二、边界元法[5] (4)

1)概述 (4)

2)基本解 (4)

3)拉普拉斯（Laplace）积分方程 (5)

4)拉普拉斯（Laplace）边界积分方程 (6)

5)拉普拉斯（Laplace）积分方程离散化与解法 (6)

6)泊松（Poisson）边界积分方程 (7)

三、结束语 (8)

参考文献 (9)

摘要

本文综述了边界元法的历史、现状及发展，并对积分方程和边界元法的原理进行了简单推导。边界元法是在经典的积分方程的基础上，吸收了有限元法的离散技术而发展起来的计算方法，具有计算简单、适应性强、精度高的优点。它以边界积分方程为数学基础，同时采用了与有限元法相似的划分单元离散技术，通过将边界离散为边界元，将边界积分方程离散为代数方程组，再用数值方法求解代数方程组，从而得到原问题边界积分方程的解。用传统的有限单元法求解不可压缩材料会遇到严重困难，但是用边界元法求解这类材料不会有任何问题。近年来随着将快速多级算法引入边界元法，使边界元法的计算效率和解题规模都有了几个数量级的提高。

关键词：边界元法积分方程边界离散快速多级算法

Abstract

This paper reviews the history, current situation and development of the boundary element method and deduced the integral equation. The boundary element method is based on the integral equation and absorbed the discrete technology of finite element method. It has the advantages of simple calculation, strong adaptability and high accuracy. It is based on the boundary integral equation, though boundary discretization discrete boundary integral equations into algebraic equations, and then by the numerical method solving algebraic equations, thus obtain the original problem solution of boundary integral equations. The solution of nearly or exactly incompressible material problems presents serious difficulties and errors when using the conventional displacement-based finite element method, because the general stress-strain equations of elasticity contain terms that become infinite as Poisson’s ratio reaches 0.5, while the boundary element method accommodates such problems without any difficulty due to the nature of the integral equations used in the analysis. In recent years, the fast multi-pole boundary element method has received much attention because some large-scale engineering design and analysis problems were analyzed faster using boundary element method than with finite element method. This new trend suggests future prospects for boundary element method applications.

Keywords：Boundary Element Method; Integral Equation; Boundary Discretization Method; Fast Multipole Algorithm

一、引言

1)什么是边界元法[1]

边界元法是在经典的积分方程的基础上，吸收了有限元法的离散技术而发展起来的计算方法。从计算格式形成的全过程来看，关键问题有两个。一个是问题的边界化，即将给定区域上的定解问题化为可以只考虑边界的问题。这一步的关键是格林公式，这是边界元法的基石。边界化的结果使问题降维，如果是各维尺度相近的大型问题，代数方程组的未知数按指数规律减少，这无疑将大大减少准备工作、存储量与机时。有限元法要将全部区域及边界离散，并要求将全部节点纳入方程进行计算，这是因为这些节点值包含在同一个封闭方程组中。边界元法的封闭方程组中只有边界结点上的未知量，解完方程组再根据需要有目的地求内部值，因此减少了计算的盲目性。

第二个关键是边界的离散化。单就离散技术本身而言，与有限元法没有很大的不同。由于只对边界离散，因此计算误差只来源于边界，这就使我们将减少误差的全部注意力放在边界上。有由于区域内由解析公式计算，这就提高了计算精度，并达到内部值的某些特殊要求，如连续性、可微性等。

计算格式形成的两个关键已经体现了计算简单、适应性强、精度高的优点，这正是这种方法的生命力所在。

2)积分方程和边界元法的发展历史[2]

1828年George Green就对位势问题提出了三个Green等式，其中包括解的积分表示式，即将问题的解表示为单层势和双层势与边界变量及其法向导数乘积的边界积分。1903年由Erik Ivar Fredholm奠基了积分方程理论。1926年Erich Trefftz提出了区别于Ritz法的积分方程边界解法，即Trefftz法。1963年M.A.Jaswon将边界积分方程直接法用于位势问题。1967年F.J.Rizzo在Quarterly Applied Mathematics发表了关于弹性静力学问题直接法边界积分方程方法的论文[3]。C.A.Brebbia于1978年出版了书籍《工程师用边界元法》，P.K.Nanerjee和R.Butterfield于1981年出版了《工程科学中的边界元法》[4]。

我国的边界元法研究起步于20世纪70年代末，清华大学的杜庆华在推动我国边界元法研究方面起了重要作用。清华大学张楚汉运用动力学边界元法与断裂力学原理提出重力坝地震断裂与拱坝裂缝扩展模型，在边坡与地下工程研究方面提出了时域边界元与离散元耦合模型。

边界元法由于形成的代数方程组的系数矩阵是满阵，因此对于求解规模有很大的制约，传统的边界元法难以处理工程实际中复杂的大规模计算问题。由于20世纪90年代开始将计