电磁散射边界元作业 1、边界元的研究背景与意义 (BEM)是随着有限元计算方法的发展而在七十年代被提出来的求解偏微分方程初边值问题的数值计算方法。由于它一般只需进行边界剖分，有降低一维进行数值计算的特性，而且所占计算机内存小、计算时间省、与有限元法相棍合能很好地解决工程实际问题等优点，目前这一方法已在各个工程领域获得广泛应用。特别是边界元法中的基本解的性质往往能恰当地反应出无穷远处的边界条件，使得用它来研究无限域问题具有其他数值方法所不能具有的一些独特的优点。 2、边界元的研究方法 边界元法是在经典的边界积分方程法的基础上吸取了有限元离散技术而发展起来的一种偏微分的数值解法。它把区域内微分方程的边值问题归化为边界上的积分方程，然后在边界上离散求解。其基础在于边界归化，如果边界归化的途径不同，可以从同一边值问题得到几种不同形式的边界积分方程，因而导致不同的边界元方法。目前边界归化一般有以下几种方法:采用 Green 公式归化的直接边界元法;采用位势理论归化的间接边界元法;还有用 Green函数代替基本解的自然边界归化，称为自然边界元法。 边界元法也有一定的局限性。首先边界元方法必须事先明确偏微分方程的基本解，这使得它往往局限于求解常系数或某些特殊变系数线性偏微分方程。对于一般的偏微分方程，如果能够采用一些特殊的技术把它们的边值问题或初边值问题化为积分方程，这些积分方程往往是奇异积分方程，甚至是通常意义下不可积的奇异性积分方程，以至于只能在广义函数的意义上去理解这些积分。这一方面在技术上带来了数值积分的困难，另一方面也使对它的数学分析复杂化。 3、实例计算 1．已知正方形柱的Ⅰ，Ⅲ边界的,ⅡⅣ边界的，求Ⅰ，Ⅲ边界的和ⅡⅣ边界的。 参考文献：《边界元法基础》上海交大出版社王元淳 Page20-24 参考资料分析了H，K 矩阵元素的求法，其中对角元素 为边界元素的长度。非对角元素，其中为P（i）点到P（j）点的距离，为P（i）点到含 P（j）点边界单元的垂直距离。求解出 H，K 矩阵后利用求出未知边界条件 MATLAB 程序： % BEM.m % 本程序用边界元方法求解正方形柱体内电位分布 %*********************************初始化*********************************** clc; clear; clear all; close all; t1=cpu time; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %% 1.常数定义 a=6; % 正方形长 N=3; %...
