电磁场数值计算之5-西安交通大学电气工程学院.docVIP

PAGE PAGE 17 三维有限元分析 §5-1 有限元分析 所有问题都是三维的，本节用非常直接的方法将二维有限元方法直接推广到三维 5.1.1 边值问题与变分公式 一般情况下的标量场边值问题： （5-1） （5-2） （5-3） 在媒质分界面上 （5-4） （5-5） 在静电场中，—电位，—介电常数，电荷密度。恒定磁场中，—标量磁位，—磁导率，。（矢量场另外讲） 等价变分问题 式中 （5-6） 与二维场相似，式（5-4）、（5-5）是自动满足的，式（5-2）是强加边界条件。 5.1.2 有限元分析 一、区域离散 以四面体单元为例， ① 一组体单元整数编码，e =1，2，3，； ② 一组四面体顶点处所有节点编码 ③ 用4*M整型数组将单元编码和节点编码联系起来(举例说明)； ④ 用一组整型数组第一类边界上节点的全局编码； ⑤ 用整型数组表示第三类边界（含第二类边界）面上三角形单元及与它们相关节点的关系，用，是面上的三角形单元总数； ⑥ 节点坐标数组；材料参数数组等。 二、建立基函数 四面体单元e内任一点位函数表示为 （5-7） 将4个顶点坐标代入上式，得到4个顶点上的未知变量表达式，解四阶方程，得到 单元体积： 代入式（5-7），整理后，可以得到标量函数用基函数和节点函数的展开式 （5-8） 其中， （注：与不相等,……） 可以证明， ，即当场点位于四面体单元的第j个节点的对面上时（含面上的3个节点），，类似二维有限元，即面上任一点的值由该面的3个顶点决定，这就保证了面两侧单元解的连续性。 三、建立有限元方程 离散后，泛函方程可以写为 （5-9） 是离散单元综述。下面用里兹法推导 1、单元方程 取一阶变分，先令 写成矩阵形式 （5-10） 式中 2、组合方程组 组合个单元， （5-11） 即 3、第三类边界条件 在或的情况下 设由个三角形面组成，上式可以写成 三角形单元内的未知函数可近似成 ，代入上式，然后对求微分，得到 写成矩阵形式为 式中 因此，式（5-11）修正为 最终 （5-12） 4、第一类边界条件 与二维问题的处理方法相同，在第一类边界节点上赋值。令 当时 以及 当时 下面针对不同的电磁场问题，分别讨论。 §5-2 静电场问题 上面建立的公式能够容易地应用于所有静电场问题，静电场通常用下面边值问题 泊松方程 （5-13） 边界条件 （5-14） 等价变分问题 当研究区域是开域或无界区域时，可以引入虚构球面来截断区域。如果虚构面离源足够远，则该面上的电位具有渐进式 （5-15） 式中，A与无关，因此 （5-16） 可作为近似边界条件，移项后，是第三类齐次边界条件。按照§5-1介绍的方法得到后，可以计算出线性四面体单元内的电场 （5-17） 它在一个单元内是常矢量，单元越小越精确。 （若用线性矩形单元时，可以发现在每一单元内，只有在p方向是常数，而在另外两个方向上均按线性变化。这是因为单元插值函数是 （5-18） 具有二次和三次项。，计算精度比四面体四节点单元高。 §5-3 恒定磁场问题 恒定磁场的边值问题为 （5-19） 常用边界条件