教案有限元法计算电机磁场.doc

用有限元法计算直流电动机的磁场分布 引言：电机磁场存在磁饱和现象，是一非线性问题再加上电机结构复杂，无法用解析法求得场的精确解,而图解法既不准又费时,采用差分法则难于处理复杂的界面形状。而有限元法是以变分原理为基础的一种数值计算方法。把要求解的边值问题转化为相应的变分问题，利用对区域D的剖分、插值，将变分问题离散化为普通多元函数的极值问题。进一步得到一组多元的代数方程组，求解得到所求边值问题的数值解rot = (1) div= 0 (2) =v (3) 式中为磁场强度,为磁通密度, v为磁阻率,为电流密度。 为便于计算,引人辅助矢函数,使=rot， (4) 从矢量分析知div(rot )=0. 故式 ( 2 ) 必满足。 由式 ( 1 )( 3 ) 和 ( 4 ) 可 得rot(v rot )= (5) 在分析电机磁场时,我们假定: 1磁场分布与轴向位置无关,忽略绕组端部的影响； 2磁化曲线是单值的,不计磁滞； 3机座外部空间导磁率为0。 由假定1知电流密度 只有轴向分量 ,并与轴向位置无关。磁场被简化为二维平面场。矢量磁位也只有轴向分量,且与轴向位置无关.于是式 ( 5 ) 简化为 （v）+(v)=－J (6) 为简单起见,这里略去轴向分量、的下角标Z。用A、J分别示、的轴向分量。这时式( 4 )简化为 =i－j (7) 2有限单元 将求解域R分成许多小的三角形单元,并使介质分界线成为三角形单元的边,也就是说使每个三角形单元只有一种磁阻率和一种电流密度。假定在每个三角形单元中是线性地变化。对顶点为K ( ,)、l(，,)、m(,)的三角形单元(图 3 ), 图3 则有 = (13) 其中 即三角形的面积逆时针方向时,面积为正;顺时针方向时,面积为负.且、、为三角形三顶点 的矢量磁位的值。 在此线性变化的假定下,该单元中各点的磁通密度彼此大小相等。因为将式 (13 ) 代人式 ( 8 ),即得 显然与坐标 ( ) 无关。由于磁阻率v只是磁通密度的函数,因此每个三角形单b中各点有相同的磁阻率。 能量积分F对每个三角形顶点(不包括边界十L上的顶点)的矢量磁位值的偏导数皆为0,F即取极小值。 这就要求下式成立,即 k=1,2,3……n (14) n为区域 R内三角形顶点( 节点 )总数(边界 + L 上的顶点除外 )。将F分成各三角形单元的能量积分之和： F= N 为三角形单元总数。其中第j项为 = 因此式（14）可以写成 k=1,2,…..n (15) 如上所述,每个三角形单元中各点的磁通密度大小相等,电流密度相 同。所以对三角形有 = = (16) 其中、分别为三角形的磁阻率和电流密度。 显然，除非顶点K是三角形的一个顶点,否则设三角形 的三个顶点为K、l、m(图 3 ),则式(1 6)中右端两项分别为 = = (17) (18) 其中为三角形单元的的面积。、、分别为三角形的三个内角。 设顶点 k 的周围有L 十 1个三角形,依次为、、、…其顶点依次为l,m,o…..t,（图4），则式（15）简化为 (19) 将式 (1 6 ) 及式 (1 7 ) 和式 (1 8 )代人式 (1 9 )得： = （20） 对每个顶点都有这样一个方程,只需将相应的三角形单元编号和顶点 编号代人即得。这样共有n个代数方程,写成矩阵形式为 (21) 其中[A]是矢量磁位顶点值的列矢量,[Q]是相应顶点的电流源列矢量,[S] 被称为磁阻矩阵。由于磁阻率v是磁通密度的函数,因此磁阻矩阵[S]是[A]的函数。此方程为非线性代数方程。至此,一连续函数 的偏微分方程被离散为有限个未知数的代数方程