第8章_边界元法概念.ppt

（2）线性单元 1）主对角线元素 采用线性单元离散曲线边界，相当于用N条直线段去逼近曲线。设定节点为线段的端点，显然节点处曲线不再光滑，如图8-8所示。根据式（8-51）和式（8-52），对角线元素为 和 式中，ci=1-θ/2π，θ是节点i所张的平面角。由于单元中r0与en相互垂直，经过与式（8-54）和式（8-56）相类似的推导，可以得出 和 2）非对角线元素 此时 根据式（8-47）和式（8-48），并注意到l=Lj（1+ξ）/2，所以有 同理可得 和 式中，r是ξ的函数， 如上所述，应用高斯求积公式，即可算得以上各系数矩阵的元素值。 8.5 典型算例 仍以长直接地金属槽中电位和电场强度分布的计算为例，具体阐述二维边界元法的实施及其计算精度的评价。该金属槽如图8-9所示，侧壁和底面电位为零，顶盖电位的相对值为100。基于场分布的对称性，此二维电场的求解区域可进而压缩而归结为D/2。这样，对应于本问题的数学模型为 采用8.4节中线性单元的离散化模式，沿计算场域边界按逆时针方向编写单元的节点号，如图8-10所示。然后，基于线性插值函数，建立系数矩阵［H］和［G］，由所得边界元方程（8-38）解出边界上的未知量后，即可计算场域内任一点的电位与场强。以上采用线性单元边界元法的数值求解过程，可应用列于附录8.1的线性单元边界元法的通用计算程序（BEM2D）来实现。为此，按该程序输入原始给定数据的子程序（INPUT）的规定顺序与格式，依次输入以下数据： 在场域边界剖分完全相同的前提下，表8-1和表8-2给出当a=20，b=17时，分别应用有限元法、边界元法和解析法计算典型场点处电位及电场强度的分布，并列出所得数值解的相对误差。由表8-1、表8-2计算结果可以看出，边界元法和有限元法对于位函数的计算，都有令人满意的计算程度。但是，在计算场强时，很明显，边界元法较之于有限元法有更高的计算精度，这是因为采用边界元法计算时，如8.4.1节所述，场强与位函数有同阶的计算精度，而用有限元法求场强时，则需再进行一次微分运算，从而导致计算精度的降低。 第 8 章 边 界 元 法 第 8 章 边 界 元 法 本章基于加权余量法，阐述了构造边界元法的数学基础———边界积分方程。由此，引入微分方程基本解和格林公式，进一步导出了对应于边界上未知量为场量φ、A、E、和H的直接边界积分方程。 以数值求解边界积分方程为目的，本章介绍了两种最基本的边界元法的构造模式：常数单元和线性单元的计算模式。并借助于高斯求积公式给出各系数矩阵的元素值。 本章最后给出二维边界元法典型应用的示例，讨论了方法实施的全过程，以及可供参考使用的计算程序，并通过本方法与其他方法计算结果的对比，进一步展示了本方法的特点。 8.1 概述 边界元法（Boundary Element Method，简称BEM）是近20余年来发展形成的一种数值计算方法。该方法的工程应用起始于弹性力学，现进而应用于流体力学、热力学、电磁工程、土木工程等诸多领域，并已从线性、静态问题延拓到非线性、时变问题的研究范畴。 边界元法是把边值问题等价地转化为边界积分方程问题，然后利用有限元离散技术所构造的一种方法，其主要特点是： 1）降低问题求解的空间维数。本方法将给定场域的边值问题通过包围该场域边界面上的边界积分方程来表示，从而降低了问题求解的空间维数。也就是说，三维问题可利用边界表面积分降维为二维问题；而二维问题则利用边界的线积分降维为一维问题。因此，有限元离散仅对应于二维曲面单元或一维曲线单元，使方法的构造大为简化。 2）方程组阶数降低，输入数据量减少。如前所述，待求量将仅限于边界节点，这不仅简化了问题的前处理过程，而且大幅度降低了待求离散方程组的阶数。 3）计算精度高。本方法直接求解的是边界广义场源的分布。根据不同的问题，广义场源可以是位势、场源或等效场源。场域中任一点的场量将通过线性叠加各离散的广义场源的作用而求得，毋需再经微分运算。此外，由于只对边界离散，离散化误差仅仅来源于边界。所以边界元法较之有限元法，可望有较高的计算精度。 4）易于处理开域问题。本方法只对有限场域或无限场域的有限边界进行离散化处理并求解，因此特别适用于开域问题。 然而，边界元法与有限元法相比较，其明显的不足之处是： 1）系数矩阵为非对称性的满阵。显然，这就引发了应用计算机求解大型离散方程组的困难，从而约束了边界元方程组的阶数。 2）系数矩阵元素值需经数值积分处理，故系数矩阵的建立需要较多的计算机时。 3）不易处理多种媒质共存的问题。 8.2 基础知识 8.2.1 格林公式 设V为空间中