讲义2012-有限元分析的后处理.doc

第7章 电磁场数值计算的后处理 各种数值计算模型，最终都归结为一个以节点或棱边物理量为未知数的离散化代数方程组。求解这一方程组，就得到所研究区域内各节点或棱边物理量的值。通常这些物理量是电位或磁位，或场矢量与电、磁位的组合。但是电位和磁位只是为简化计算而引入的辅助量，还不能直接供研究者和设计工程师使用。为了应用数值分析的结果进行电气设备的设计和优化，需要根据已经求出的节点或棱边物理量，进一步算出诸如磁感应强度、电流密度、电磁力与力矩、磁场能量、损耗等局部和总体电磁量。此外，由于上述计算结果所形成的大量数据难于检索分析，还需要将它们转化成图形形式。这些工作统称为数值分析的后处理。本章主要介绍节点有限元分析的后处理方法，包括局部与总体电磁量的计算和力与力矩计算，这些方法容易推广到其它数值分析方法的后处理过程中去。关于矢量场的图示，可参考文献[1]、[2] 和其它有关计算机图示学文献。 7.1 局部与总体电磁量的计算 下面以正弦稳态涡流场为例，由已知节点矢量磁位和标量电位，给出各种电磁量的离散化表达式。对于静态磁场，则与电导率有关的项应消失。 在2.3.2节中已经指出，对于任一单元e，矢量磁位的插值函数可用该单元各节点处的磁位近似表示为 (7-1) 上式亦即 ， ， (7-2) 类似地，单元中标量电位的插值函数为 (7-3) 7.1.1磁感应强度 由于 ，因而单元磁感强度的三个分量、、分别为 ， ， (7-4) 将式 (7-1) 代入 (7-4) ，有 (7-5) 7.1.2 电流密度 考虑到式 (1-11)和(1-20)，单元电流密度与和的关系为 (7-6) 其中为单元源电流密度。将式(7-1)、(7-2) 代入上式，可得单元电流密度的各分量为 (7-7) 7.1.3 电感、能量与涡流损耗 根据经典电磁理论，当媒质为线性时，磁场能量密度的瞬时值为 (7-8) 其中H为磁场强度的瞬时值。在正弦电磁场中，H可表示为 (7-9) 磁场强度瞬时值与复数值的关系为 (7-10) 式中表示取实部算符；为电角频率；为最大值复矢量，又称复振幅， (7-11) 式(7-10)是磁场强度复矢量的定义式，但不便于用来进行运算。对时间的依赖关系还可用更简单的缩写形式表示[3]： (7-12) 式中 为的复共轭矢量。式(7-12)方便地解决了场量的瞬时值与复数值的混合运算问题。将(7-12)代入(7-8)，有 (7-13) 上式的第二项展开后为随时间变化的二倍频率余弦函数，在一周期内的平均值为零；第一项不随时间变化，正好等于磁能密度在一周期内的平均值，记为， (7-14) 式中，符号表示复矢量各复数分量模平方的和，它等于瞬时值矢量的模平方平均值的两倍。根据式(7-14)，体积V内的磁场能量的时间平均值为 (7-15) 式中的下标R和I分别表示复矢量相应分量的实部和虚部，上标e表示单元值；为单元i的体积，E表示体积V内的单元总数，为单元i中的有效磁导率（参见8.3节）。 平均磁场能量得到之后，容易求出相应的电感参数。由于 (7-16) 因而 (7-17) 其中I表示源电流的正弦有效值。 由经典电磁理论可知，实心导体区域中涡流损耗密度的瞬时值为 (7-18) 式中J表示涡流电流密度的瞬时值，为电阻率。应用与 (7-12) 类似的表示法，式 (7-18) 可进一步写成 (7-19) 比较式 (7-19) 与 (7-