电磁学发展与计算方法.DOC

电磁学发展与计算方法 摘要：介绍电磁学计算方法的研究进展和状态，对几种富有代表性的算法做了介绍，并比较了各自的优势和不足。同时也说明电磁学与电之间的相互关系。 关键词：电磁学、矩量法、有限元法 一 引言 现代人的生活，似乎离不开电。电灯、电话、电视、 HYPERLINK /map/21----%B5%E7%D3%B0----/ \t _blank 电影、计算机、电冰箱…，样样都是生活必须用品。一旦停电，日子不知怎么过。而我们生活中必不可少的电与我们物理中的电磁学息息相关。例如古代的电磁观察与应用，1936年，考古学家在巴格达附近挖出了一些铜罐，罐中铺了沥青，沥青上插着铁条。在大约同一地点，还发掘出了一些镀金物品。有研究者便认为这些铜罐就是巴比伦人发明的电池，而镀金物（如果是电镀）是这些东西确是电池之证据。而这些东西，其年代有早到公元前2000年以上的。如果这是真的，巴比伦人领先了近代电池(伏特，1793)与电镀(1800-35)，将近四千年。别的文明在电磁方面就没有这样可惊的成绩了。古希腊人发现了琥珀、毛皮等摩擦可以生电，至今英文electricity的字根，尚是希腊文的琥珀。但对他们说来，天上的雷电，仍然是宙斯大神的脱手武器。中国人很早就知道天然磁石会吸铁，带电物会吸小物体（东汉王充27-97「论衡」电磁力之记述：「顿牟拾介，慈石引针」），以及利用磁针导航，甚至对磁偏角有所记述(方以智，~1600)。「磁针导航」这技术，传到西方，促成了西方的「大探险时代」(15-16世纪。1492哥伦布发现美洲，1498达伽马绕过好望角到达印度，1519-22麦哲伦环绕世界一周，称为「三大航海」。他们都用磁针罗盘。)也引起了十八世纪以后的殖民主义。这些电磁的观察与应用，可以使我们感叹古人之智能，特别是巴比伦电池。但巴比伦电池即使是事实，对日后电磁学发展，却没有什么影响。摩擦生电与磁性现象却在停滞千余年之后，在十八世纪的西欧，成为电磁学发展的出发点。 二 关于电磁学的集中计算方法 1 电磁场数值方法的分类 电磁学问题的数值求解方法可分为时域和频域2大类。频域技术主要有矩量法、有限差分方法等，频域技术发展得比较早，也比较成熟。时域法主要有时域差分技术。时域法的引入是基于计算效率的考虑，某些问题在时域中讨论起来计算量要小。例如求解目标对冲激脉冲的早期响应时，频域法必须在很大的带宽内进行多次采样计算，然后做傅里叶反变换才能求得解答，计算精度受到采样点的影响。若有非线性部分随时间变化，采用时域法更加直接。另外还有一些高频方法，如gtd，utd和射线理论。 从求解方程的形式看，可以分为积分方程法（ie）和微分方程法（de）。ie和de相比，有如下特点：ie法的求解区域维数比de法少一维，误差限于求解区域的边界，故精度高；ie法适合求无限域问题，de法此时会遇到网格截断问题；ie法产生的矩阵是满的，阶数小，de法所产生的是稀疏矩阵，但阶数大；ie法难以处理非均匀、非线性和时变媒质问题，de法可直接用于这类问题〔1〕。 2 有限元方法 有限元方法是在20世纪40年代被提出，在50年代用于飞机设计。后来这种方法得到发展并被非常广泛地应用于结构分析问题中。目前，作为广泛应用于工程和数学问题的一种通用方法，有限元法已非常著名。 有限元法是以变分原理为基础的一种数值计算方法。其定解问题为： 应用变分原理，把所要求解的边值问题转化为相应的变分问题，利用对区域d的剖分、插值，离散化变分问题为普通多元函数的极值问题，进而得到一组多元的代数方程组，求解代数方程组就可以得到所求边值问题的数值解。一般要经过如下步骤： ①给出与待求边值问题相应的泛函及其变分问题。 ②剖分场域d，并选出相应的插值函数。 ③将变分问题离散化为一种多元函数的极值问题，得到如下一组代数方程组： 其中：kij为系数（刚度）矩阵；xi为离散点的插值。 ④选择合适的代数解法解式（2），即可得到待求边值问题的数值解xi（i＝1，2，…，n） 3 矩量法 很多电磁场问题的分析都归结为这样一个算子方程〔2〕： l（f）＝g（3）其中：l是线性算子，f是未知的场或其他响应，g是已知的源或激励。 在通常的情况下，这个方程是矢量方程（二维或三维的）。如果f能有方程解出，则是一个精确的解析解，大多数情况下，不能得到f的解析形式，只能通过数值方法进行预估。令f在l的定义域内被展开为某基函数系f1，f2，f3，…，fn的线性组合： 其中：an是展开系数，fn为展开函数或基函数。 对于精确解式（2）通畅是无限项之和，且形成一个基函数的完备集，对近似解，将式 （2）带入式（1），再应用算子l的线性，便可以得到： m＝1，2，3，… 此方程组可写成矩阵形式f，以解出f。矩量法就是