用有限元法分析异步电动机.pdfVIP

用有限元法分析异步电动机 newmaker 【摘要】阐述了有限元方法和牛顿迭代算法，建立了电机有限元磁场模型， 并对异步电动 机进行性能分析，在实际电机上得到验证。 关键词：有限元法；异步电动机；牛顿迭代法 1、引 言 电气工作者把分布参数系统问题归结为集中参数系统问题求解，是屡见不鲜的。有限元数 值解法，就是这种相同的思考方法。但是，有限元问世后，在电工行业中却远不及力学领 域等应用得充分。 有限元法的特点是适用于求解各种形式（几何上、物理上）复杂的问题，精度高，通用性 强，对问题的处理既彻底又系统，适用于采用电子计算机方式。它本是线性问题的解法， 但通过迭代法（如牛顿-拉裴森迭代法）也能巧妙地解决非线性问题。用来求解电磁场问题 远较电工行业中常用的图解法、电解槽法等优越。因此，随着科学技术的发展，各类超高 压、大容量、高精度的电工产品的研制，其磁场等的精确计算直接关系到该产品的优异性 能和技术经济指标。有限元法这一有力工具，在电工行业中，在电磁研究领域里，越来越 有广泛的应用。 用有限元法对异步电动机性能进行分析的方法，在实际电机产品中得到验证。在磁场中向 量场函数B的旋度： (1) 通常用向量位A来描述，即 (2) 则磁场的基本方程 (3) 式中 B——磁感应强度 μ——介质常数 J——电流密度 从电机中截出横断面，在这里的磁场可视为平行平面场。 设选定坐标使B和A在Z轴方向的分量为零，向量磁位A、电流密度向量J只有Z轴方向 的分量，令磁阻率 ，于是方程（3）化为二维形式。 (4) 电机定、转子为硅钢片迭成，由于铁磁物质受饱和影响，磁化特性非线性,如图1示，磁阻 是磁场强度H的函数，所以方程（4）是一个偏微分方程，直接求解这样的非线性偏微分 方程是比较困难的，使用数值计算法——有限元法，可得到较高精度的解。 图1 DG41(厚度0.35mm)硅钢片磁化曲线 2、有限元法 有限元法的基本原理是以变分原理和剖分插值为基础的一种数值计算方法，把所要求的电 磁场问题即偏微分方程的边值问题化为与之等价的变分问题即所谓泛函数的极值问题。从 而得到一个高阶非线性方程组。最后求解方程组，即得待求的电磁场问题的近似解。 为了用有限元法解方程（3），首先必须确定一个恰当的泛函。根据微积分中函数的极值原 理，可以证明（证明从略）泛函。 （5 ） 与式（4 ）等价，即上述偏微分方程的定解问题等价为条件变分问题。条件变分问题与其等 价的偏微分方程边值问题比较，在边界问题上要简便得多，且在这样的基础上才能进行剖 分插值。 根据电机的结构，转子位置不同时磁场的分布情况不一样，我们尽量利用结构磁场的对称 性，缩小其求解区域。将确定的求解区域，剖分成有限个三角形单元，剖分时，不同介质 的交换面必须是三角形单元的边，关键部分剖分密度大一些，避免出现太尖太钝三角元， 以保证计算的精度。 2.1 剖分插值 将问题区域剖分成若干三角形单元，其顶点为i,j,m, 假设位函数用式(6)来近似： A(x,y)=a+bx+cy (6) 式中的系数a,b,c 可由3 个联立的独立方程来确定。 假定位函数值为顶点值Ai, Aj,Am, 将3 个顶点的位值Ai, Aj,Am 及其坐标（x1,y1 ）,（x2,y2 ）， （x3,y3 ）代入式（6），求解联立方程，得系数a,b,c，并将结果代入式（6 ），就得到 将x,y 和逆系数矩阵的元素组合为新的位置函数，则可写成 (7) 式中 它仅是位置的线性函数，而Δ代表三角形的面积。可通过下标的循环置换而获得，可以证 明式（7）是三角形三个顶点的内插函数。 2.2 有限元的求极小值 泛函式（5）在非线性情况下的离散，基本上按照线性问题所采用的同样方法进行。假定问 题的求解域D离散成