电磁场数值计算.pdfVIP

电磁场数值计算中数学建模方法研究 李 伟, 严登俊，朱长江 （河海大学电气工程学院，江苏 南京，210098 ） 摘 要：总结了几种电磁场数值计算中数学建模方法—— 面体的四个顶点接近共面，体积接近零，尽管四 网格剖分方法，分析了他们的优缺点，并对网格剖分方法 个表面的三角形形状良好，也必须要消除薄元。 的未来发展方向和发展趋势作了一定的探讨。 关键词：有限元；网格；自动生成 0 引言 随着计算机技术的发展，人们对于科学和工 程计算的规模、复杂程度、效率、精度等各方面 a 图：生成网格的质量差 有了更高的要求，以期能更真实、准确地反映要 求对象。在有效的解决科学和工程计算问题中， 数学建模中网格剖分技术被广泛应用到实际问题 中，从而作为数学建模的先决条件的网格剖分技 术就显得尤为重要。本文就对电磁场数值计算中 的几种网格剖分方法进行了比较，并对网格剖分 [1] 技术的未来发展趋势作了一些分析 。 b 图：生成网格的质量好 1 网格剖分的基本方法 如图1 所示的a 图是不符合Delaunay 三角划分， 1.1 Delaunay 三 角 划 分 法 （Delaunay 其中 ∠BAD大于90º, ∠ADB角度很小,所以这样 Triangulation Method ） 划分容易产生薄元，b 图是符合Delaunay 三角划 Delaunay 三角划分(简称 DT)也是结点连元 分。常用的薄元处理方法有两种：一是结点抖动 法，是通用的全自动网格生成方法之一。在平面 法，将一个薄元的一个结点沿着某个方向移动， [2] 域的Delaunay 三角划分中，DT 法有两大优点 ： 以增加薄元的“厚度”，提高单元质量，但抖动效 首先，遵循最大-最小准则，即所有的三角形单中 果在实际应用中有许多困难，而且在某些情况下 最小角之和均大于非 Delaunay 划分所形成的三 会失效；二是薄元分解法，波源一般都相当于比 角形最小角之和，此时结点分布中不存在四点及 较薄的面片，可以在面片的形心处插入一个新的 四点以上共圆，从而避免了生成小内角的薄元。 结点，该结点可以同此薄元临近的4 个单元形成 另外，DT 法还有空外接圆特性，即Delaunay 三 新的单元，来代替薄元。 角划分中任意三角形的外接圆( 四面体为外接球) 1.2 映射法 内不包括其他结点，构造外接圆的基本过程：先 映射法出现于20 世纪70 年代，是最早采用 构造一个大外接圆，将所有结点都包含进去，然 的网格生成方法。根据映射函数的不同，主要分 后找出已有三角形中哪些三角形的外接圆内包含 [4] 为等参映射、超限映射和一致映射 。根据形体 新加入的结点，删除这些三角形内离心结点距离 边界的参数方程，利用适当的映射函数，将待分 最近的一条边，将新结点与周围的老结