两维无源场域边界元法及其改进.pdfVIP

维普资讯 征西电力 第ls卷 I，’i年 第3期 (1。s)盯 两维无源场域边界元法及其改进 江西工业大学 雷 良饮 内謇摄鼻l本文首先介绍 自7O年代来兴起的边界元法，然后，作者针对两维场边界元海 的现行方案，作出了改进，即边界法向场强采用线性插值，但方程数不增加，实俯表碉对 凸 形边界可行。 引 盲 式中， 和 均为定义于域S中连续并有二阶 采用边界元法 (BEM )求 解 电磁场有 导数的任意标量函数。 它的优越性。首先，它可直接求出边界上的 2．现在，我们令 (1)式中的 为无蔼【 法向场强，同时，和有限元法比较，它的方 场域S中的电位 函数，因而在S申 程数 目大为减少，有助于降维、节约内存、 V = 0 …… (2) 缩短机时。但是，在边界元法中，对于边界 再令 =一L I— I …… (3) 单元上的法向场强，一般难于采 用 插 值 函 式中， 为坐标原点至S域中任一定点Pt的 数， 因为在节点处的法线方向变化不连续。 距离矢量，而r为坐标原点至 S域中积分甩 作者提出的法 向场强线性插值彦，考虑了节 动点的距离矢量 。 点处法向场强的不连续，对第一类边界问题 显然，对动点坐标进行微分运算，必有 特别有效。实例表 明，作者方法的精度 比一 1}=一2)zb(r—rk) …… (4·a) 般方法大有提高。 关键词 ·电磁场数值计算 边界元法 而』6(一一rk)ds=』。¨一一rk)dS=1 一 、 理论基碡 … … (●·b) 1．设有一乎面场域S，其边界 为 闭合线 式中，8可理解为包围定点P但半 径无限小 L。据矢量分析中的格林第二定理，垣 有 公 的匿。 式 l 将等式(2)．(3)、 (4·a)爰 (4·b) 代入式 (1)，得 2= (L 一 1(．加acp)+击 (L I-~kI)dL ：￡… (5) ■ 因此，若 已知边界L上的电位 值 甲及 电 fs(vt一甲vz)ds 场强度的法向分量En=一-鲁，便可求出 域内任一定点P政 的电位值 k。 (》吾一鲁)dL…_【1) 3．一般典型的边界向题，仅知边界L上 维普资讯 埠 (硒 ) 让西电力 第T5卷 1年 革a期 的电位分布甲(第一类边界问题 )或 仅知L 2．边界元法的第一步是在边界 L上设定 上的法 向场强Ea(第二类边界 问题 )，固 许多的节点，一般逆时针编号，顺次为节点 此无 1丧直接 由 (5)式求 出 k。 l