无网格方法及在电磁场数值计算中应用.pdfVIP

无网格方法及其在电磁场数值计算中的应用 尹华杰 （华南理工大学广州510641 ） 摘要：本文之目的在于向电磁场数值计算界推介无网格法这一新的数值计算方法。文中着重针 对人们在第一次遇到无网格法时，可能提出的疑问，进行说明与阐述，并分析讨论了无网格法 在电磁场数值计算中的一些潜在应用可能。 关键词：无网格法，有限元法，电磁场数值计算，移动最小二乘法 1 简介 作为一种新的数值仿真技术，无网格法 （Meshless method ，简称为 MLM ）主要用于代替 和弥补有限元法 （FEM ）的不足。在网格需要反复生成的场合，如运动体仿真、区域变形问题、 裂纹生成仿真，以及存在奇点的场合，MLM 特别有用。然而 MLM 方法的边界条件以及交界 条件的难于处理，使得其研究和应用主要局限在机械、材料领域。 在电磁场数值计算 （EM ）领域，也存在利用MLM 的优点以克服FEM 缺点的许多理由。 从已发表的文献看，尽管MLM 在EM 领域的研究、应用尚十分有限，但它已在国外引起 EM 领域研究者的兴趣。经过仔细收集、分析关于MLM 的文献，我们总结出了7 个方面的问题， 在以下加以阐述，以期对EM 领域的研究者理解MLM 方法有所助益。 本文的第2 部分，着重回答EM 领域为何需要MLM 方法的问题。第3 部分回答什么是MLM 方法。第4 部分介绍MLM 方法的种类。第5 部分回答当前谁在研究MLM 方法的问题。第6 部分回答在EM 领域中，谁在应用MLM 方法的问题。第7 部分分析MLM 方法存在的问题及 局限。第8 部分回答EM 计算领域的研究者们面对MLM 方法能够做些什么的问题。为方便起 见，第9 部分对MLM 方法的数学理论作简单介绍。 2 为何需要MLM 方法？ 众所周知，在偏微分方程的数值计算中，FEM 是十分成熟而强大的工具。FEM 以网格单 元为基础，尽管网格自动剖分软件通过商业途径已很容易获得，但对普通研究者，它还过于昂 贵。而在技术上，FEM 则存在以下问题：（1）当求解问题区域变形时，FEM 要求区域反复剖 分，这将成为计算机的沉重负担；（2 ）FEM 仅仅提供 C0 连续的近似解，当使用标量位或矢量 位求解电磁场时，场解沿单元边界不连续；（3 ）网格单元形状对解的收敛性影响严重，在一些 逆向问题中 （如电磁设备的形状设计），为使求解收敛，即便区域变形是比例性质的，也必须 反复剖分。为了克服FEM 的这些缺陷，人们开发了MLM 方法。 3 什么是MLM 方法？ 根据E. Oñate(1996)的定义，任何方法，如果其近似解能够严格依节点 （而非单元）构建， [1] 就是MLM 方法 。MLM 法原理比较抽象，根据笔者的心得，用以下分步骤的操作来解释MLM 方法，有助于理解 （图1）： 本文受国家自然科学基金资助(批准号) ，为作者在美国伊利诺伊大学 （UIUC ）电磁计算中心作 访问学者时所作小组报告的翻译稿、略有增改。 1/9 依节点构建近似解 获得积分方程 获得代数方程组 边(交)界条件处理 获得形状函数 (伽辽金法等) (高斯积分法) (拉格朗日乘子法等) (MLS 方法等) 图1 MLM 方法图解 1°在求解区域中布一系列节点，记为x (I=1~N , N 为节点总数) 。节点可均布，也可以随机 I 分布。 2 °为每节点 x设置影响域 （或称支持域）。影响域以 x为中心，通常为圆形或矩形，大 I I 小可相同，也可不同。 3 °选定