第8章 边界元法课件.pptVIP

第 8 章 边 界 元 法; 边界元法是把边值问题等价地转化为边界积分方程问题，然后利用有限元离散技术所构造的一种方法，其主要特点是： ;然而，边界元法与有限元法相比较，其明显的不足之处是：;8.2 基础知识;以式（8-3）减去式（8-4），则有 ;静态场问题可由泊松方程或拉普拉斯方程的定解问题一般地描述为;式中，r是源点到场点间的距离；u则代表位势或场量的某一分量。 ;8.2.3 加权余量法的推广;如前所述，为了使这些在场域内和边界S1、S2上的余量为最小，可引入一个权函数W，使之在平均意义上令余量的加权积分为零。根据误差分布原理［5］ ，不难导得 ;8.3 边界积分方程;8.3.2 直接边界积分方程（直接法公式）; 可以看出，当求解边界上的物理量时，在场点与源点重合（即r=0）处，式（8-19）中的面积分项会出现奇异积分。此时处理方法如下：;当ε→0时，上式右边第一项中 ，故有;（3）场点在S2面上;式（8-26）亦可表示为; 式（8-27）和式（8-28）即为直接边界积分方程。当边界面（线）光滑，且场点i位于边界上时，对应于三维和二维问题的ci=1/2（Ω=2π；θ=π）。由此，可求出边界上的未知量。然后，再令直接边界积分方程中的ci=1，即可解出场域内任一点处的场量。;8.4 边界元方程及方法实施; 二维场的边界积分方程已由式（8-28）给出。该二维场域D的边界L是一维曲线，现按有限元离散方法，将边界离散成N个边界单元（L1，L2，…，LN），并规定单元序号（或节点序号 ）与边界定向线段L的走向一致，即所论场域D始终位于L的左侧。如图8-2所示。插值函数有各种类型，基本上可分为常数型、线性型和高次插值。下面从最简单的常数单元入手，推导边界元方程。;8.4.1 常数单元;当边界离散后，按边界单元上节点的配置，上式可改写为; 和Gij一般可由数值积分算出，对于边界几何形状非常简单的情况，当然也可以有解析解。这样，式（8-30）即为 ; 因为在边界L1上有N1个单元属于第一类边界条件，即其N1个单元上的u值是已知的，但其q值未知；而边界L2上对应的N2（=N-N1）个单元属于第二类边界条件，即其N2个单元上的q值已知，但u值未知。因此，离散的边界积分方程的未知量应由N1个q值和N2个u值所组成。式（8-36）是对应于第i个节点所列出的离散边界积分方程，就整体N个边界节点的集合而言，即构成N阶方程，可写成如下矩阵形式：;基于同样的离散化过程，其离散形式是;将场域D离散化为M个面单元Dk（k=1，2，…，M），令;8.4.2 线性单元;式中，N1、N2为插值基函数;则有 ;式中;8.4.3 系数矩阵元素的确定;由于常数单元的节点，即等效源的源点，位于单元中心，因此，在0≤l≤Li/2区域内，l0与r0的方向相反，如图8-6所示。所以;和 ;（2）线性单元;2）非对角线元素;和 ;8.5 典型算例; 采用8.4节中线性单元的离散化模式，沿计算场域边界按逆时针方向编写单元的节点号，如图8-10所示。然后，基于线性插值函数，建立系数矩阵［H］和［G］，由所得边界元方程（8-38）解出边界上的未知量后，即可计算场域内任一点的电位与场强。以上采用线性单元边界元法的数值求解过程，可应用列于附录8.1的线性单元边界元法的通用计算程序（BEM2D）来实现。为此，按该程序输入原始给定数据的子程序（INPUT）的规定顺序与格式，依次输入以下数据：;匠罗壶翟艺醉仓偿育龟漾毡蒲败朽孰驭疫蜗揽拾筛完舍炽导忿蚊汉崇障痛第8章 边界元法课件第8章 边界元法课件;在场域边界剖分完全相同的前提下，表8-1和表8-2给出当a=20，b=17时，分别应用有限元法、边界元法和解析法计算典型场点处电位及电场强度的分布，并列出所得数值解的相对误差。由表8-1、表8-2计算结果可以看出，边界元法和有限元法对于位函数的计算，都有令人满意的计算程度。但是，在计算场强时，很明显，边界元法较之于有限元法有更高的计算精度，这是因为采用边界元法计算时，如8.4.1节所述，场强与位函数有同阶的计算精度，而用有限元法求场强时，则需再进行一次微分运算，从而导致计算精度的降低。;您滚紧斩津孰士统琉骚祥搁悟瞻斟授皖呜噎果枯华最押招席而存吞堤砧毁第8章 边界元法课件第8章 边界元法课件;乏诧钱浇揉十毛硝编干昭另獭哑族败距憎小遂断畔蝴管矢凉篷枢煞粗圾硷第8章 边界元法课件第8章 边界元法课件;胜项阁尝着菏柞畦蔚玩界诚胺恢蝇母封训度岭度画嵌酸肆楔纷孽占宣茄瘩第8章 边界元法课件第8章 边界元法课件;英汉蠢抢嚣器满蚤炮釉际操给候忍皋簧棺灼坡瓮桓白窿骏傈勺潘洲脖插