边界单元法基础(直接法讲义)教程分析.doc

边界单元法基础（直接法） 一、概述 近年来在边界法方面人们发表了大量的文章和著作。这些方法是以不同的名称而提出来的，如“边界积分方程方法”“边界积分解”，等等。这种方法的数值解形式是把所考虑的域的边界划分为一系列的单元。边界单元法简称BEM是七十年代兴起的一种新的计算方法。它将边界上的广义位移和广义力作为独立变量且同时用满足场方程的奇异函数（源函数）作为加权函数。所以，它是一种特殊格式的加权余量法。 边界元法只需将求解域的边界划分成单元，故使求解问题的维数降低，如三维问题可转变成二维问题求解。二维问题可化为一维问题。因而，输入数据大为减少，计算时间缩短。由于它只对边界离散，故离散误差仅为来源于边界，而域内变量可由解析式的离散形式直接求得。因此，提高了计算精度。求域内变量时，只须改变其数量和坐标位置即可。二、泊松方程的边界单元法 1．积分方程的建立和基本解 为了说明边界单元法的积分方程是如何由加权余量法推导得来的，我们以泊松方程为例来阐述其全部求解过程，这对了解其它问题的边界元法求解是有益的。考虑势函数φ。它在域内满足微分方程，即 （在Ω上） （9-1） 在边界上满足边界条件，即 （在上） （在上） （9-2） 如图9-1所示，求解域Ω的总边界。 图9-1 位势问题的域和边界 可以证明，对于泊松方程或拉普拉斯方程（），一般加权余量表示式为 （9-3） 式中，W是权函数，在边界法中可令W，称为相应于方程（9-1）的基本解，后面描述它。 对式（9-3）左边第二项进行分部积分，即 （9-4） 式中右边第一项的被积函数形式称爱因斯坦求和约定，即 把式（9-4）右端第一项再次分部积分，得 所以 把上式代入式（9-3）或去（），考虑到，把左右两边界积分合并时，得 （9-5） 这就是求解泊松方程的积分方程。如果f = 0，即为求解拉普拉斯方程的积分方程。 下面来导出相应该方程的基本解。假如有一单位势作用在物体上的P点，那么基本解就是满足下列方程 ▽2φ*（p,q）+δ（p-q）=0 （9-6） （9-7） 式中，是狄拉克δ函数，它具有如下性质： （9-8） r为p点到q点之间的距离。方程（9-6）或（9-7）的解称为拉普拉斯方程的基本解。 在式（9-7）中，当时，则为 （9-9） 它的解为 （9-10） 式中，C１和Ｃ2为积分常数，可由r = 0时的条件确定。由于 （9-11） 应用格林公式 （9-12） 由图9-2（b）所示，以p为圆心，以ε为半径作一小圆，将式（9-10）代入式（9-12）并使时取极限，于是得 由式（9-11）得 （9-13） 由于是势函数，C1可取作零，所以 （9-14） 式（9-14）即为二维拉普拉斯问题的基本解。 同理可求得三维拉普拉斯问题的基本解为 （9-15） 现在根据式（9-5）来建立积分方程。由于 于是得 （9-16） 式中； 考虑时，可写成 式（9-16）就是域内任一点p的函数值与边界的积分关系，除了微分方程非齐次项f的积分项外，其余都是边界积分表达式，故称积分方程。 图9-3 用半球包围的边界点 下面我们来建立边界解的积分表达式。如果在式（9-16）中把内点p取到边界上P点（三维举例），如图9-3所示。这时，为了避免奇性，需要在P点附近稍加改变，即以P点为中心，以为半径作一小半球包含P点。这样P点仍为内点。此时把边界分为两部分：一部分是鼓起的半球部分，用表示；另一部分即。假如点P取在上（上也一样分析），对于式（9-16）在边上的积分分为两部分，则 （9-17） 式（9-17）右端第二工面当时取极限，并把三维基本解代入时，则对光滑表面有 当时，而式（9-17）右端第一项仍为上积分，即 式（9-16）右端第一项在上积分同样分为两部分，即和。而 而 如果P点取在上也可推出同样的结论。所以，把P点取在边界上时，则式（9-16）为 （9-