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

《用有限元法分析异步电动机》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、用有限元法分析异步电动机第11、引言电气工把分布参数系统问题归结为集中参数系统问题求解，是屡见不鲜的。有限元数值解法，就是这种相同的思考方法。但是，有限元问世后，在电工行业中却远不及力学领域等应用得充分。有限元法的特点是适用于求解各种形式（几何上、物理上）复杂的问题，精度高，通用性强，对问题的处理既彻底又系统，适用于采用电子计算机方式。它本是线性问题的解法，但通过迭代法（如牛顿-拉裴森迭代法）也能巧妙地解决非线性问题。用来求解电磁场问题远较电工行业中常用的图解法、电解槽法等优越。因此，随着科学技术的发展，各类超高压、大容量、高精度的电工产品的

2、研制，其磁场等的精确计算直接关系到该产品的优异性能和技术经济指标。有限元法这一有力工具，在电工行业中，在电磁研究领域里，越来越有广泛的应用。用有限元法对异步电动机性能进行分析的方法，在实际电机产品中得到验证。在磁场中向量场函数B的旋度：(1)通常用向量位A来描述，即(2)则磁场的基本方程(3)式中B——磁感应强度μ——介质常数J——电流密度从电机中截出横断面，在这里的磁场可视为平行平面场。设选定坐标使B和A在Z轴方向的分量为零，向量磁位A、电流密度向量J只有Z轴方向的分量，令磁阻率，于是方程（3）化为二维形式。(4)电机定、转子为硅钢片迭成，

3、由于铁磁物质受饱和影响，磁化特性非线性,如图1示，磁阻是磁场强度H的函数，所以方程（4）是一个偏微分方程，直接求解这样的非线性偏微分方程是比较困难的，使用数值计算法——有限元法，可得到较高精度的解。图1DG41(厚度0.35mm)硅钢片磁化曲线2、有限元法有限元法的基本原理是以变分原理和剖分插值为基础的一种数值计算方法，把所要求的电磁场问题即偏微分方程的边值问题化为与之等价的变分问题即所谓泛函数的极值问题。从而得到一个高阶非线性方程组。最后求解方程组，即得待求的电磁场问题的近似解。为了用有限元法解方程（3），首先必须确定一个恰当的泛函。根据微

4、积分中函数的极值原理，可以证明（证明从略）泛函。（5）与式（4）等价，即上述偏微分方程的定解问题等价为条件变分问题。条件变分问题与其等价的偏微分方程边值问题比较，在边界问题上要简便得多，且在这样的基础上才能进行剖分插值。根据电机的结构，转子位置不同时磁场的分布情况不一样，我们尽量利用结构磁场的对称性，缩小其求解区域。将确定的求解区域，剖分成有限个三角形单元，剖分时，不同介质的交换面必须是三角形单元的边，关键部分剖分密度大一些，避免出现太尖太钝三角元，以保证计算的精度。2.1剖分插值将问题区域剖分成若干三角形单元，其顶点为i,j,m,假设位函数

5、用式(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）在非线性情况下的离散，基本上按照线

6、性问题所采用的同样方法进行。假定问题的求解域D离散成一组互不重叠的有限单元，并着重考虑一个单元，在单个单元内，磁位A可用式（7）表示，借助于这一代换，ARGIN-TOP:0px;MARGIN-BOTTOM:0px;LINE-HEIGHT:17pt"align=center>(8)式中E0——求解区域剖分的总单元数L0——求解区域剖分的总节点数即式中(9)将式（9）代入式（8）且将代入，则方程（8）变成短阵形式：SA=J（10）式中A——节点磁位值向量J——电流密度列向量，各项为S——单元系数矩阵，包含下列各项磁阻率υ是A值的函数，而且还与

7、磁场相关。由于铁磁材料的饱和，故方程（8）是非线性方程组，对于空气介质中的单元，则按其材料及所在频率下的磁化曲线计算磁阻率。为了适用电子计算机计算，进行分段插值，实测的磁化曲线，如图1示。从曲线拐弯处开始取曲线上25个节点Bi和对应的函数Hi(i=1,2,…,25),由B和H构成拉格朗日多项式。从而得2.3求解非线性方程组的牛顿迭代法设A为所求的准确解，A(k)表示一个不准确但却是适当接近A的估计值。A(k)=A-ΔA(k)(11)将ARGIN-TOP:0px;MARGIN-BOTTOM:0px;LINE-HEIGHT:17pt"align=

8、center>（12）若略去泰勒级数中超过第二项的各项则方程（12）就提供了计算A(k)对A的偏差的方法ΔA=-P-1V(13)式中，P点牛顿迭代的雅可比矩阵，其元