电磁场数值分析.docx
- 文档编号:1867658
- 上传时间:2022-10-24
- 格式:DOCX
- 页数:13
- 大小:247.95KB
电磁场数值分析.docx
《电磁场数值分析.docx》由会员分享,可在线阅读,更多相关《电磁场数值分析.docx(13页珍藏版)》请在冰豆网上搜索。
电磁场数值分析
电磁场数值分析
电和磁现象在自然界普遍存在,两者相互依存形成一个不看分割的整体。
电能产生磁,磁能生电。
很早以前人们就注意到电现象和磁现象,但是两者之间的这种相互联系在很长的一段时间内都没有被人们认识。
直到奥斯特首先发现了通电直导线周围存在磁场这一现象人们才开始把电和磁放在一起来研究。
然而这个
时候人们依然没有办法揭示电和磁中间的秘密,只是停留在实验研究阶段,没有
形成科学的理论。
1831年法拉第发现了电磁感应定律,从此电和磁的计算可以量化了,人类历史也开启了一个新的时代一电气时代。
由于法拉第的杰出工作,电和磁不再是不可触摸的了,人们已经掌握了运用它的钥匙。
在法拉第之后,另一位杰出的科学家麦克斯韦则更进一步,建立了麦克斯韦方程组,电和磁的理论已经到了相当完美的程度。
现代电机,不管结构多么复杂,都是基于法拉第电磁感应定律和麦克斯韦方程组的原理来运行的,其电和磁的相关量都可以利用这两个定律来进行精确地分析,在设计电机时,我们也是基于这两个定律对电机的电磁过程来进行精确的设计,从而设计出理想的电机。
学会电磁场分析,主要是基于麦克斯韦方程组的相关计算,对电机的学习非常重要。
它为我们今后的学习打下基础。
在学习过程中,主要要把握以下几个度之间的关系:
梯度、旋度、散度,这三者的变换正体现了电和磁之间的转换。
一基本原理
电磁场的内在规律由电磁场基本方程组一麦克斯韦(Maxwell)方程组表达。
这些方程是由麦克斯韦对大量实验结果及基本概念进行了数学加工和推广归纳而成的。
麦克斯韦方程组是分析和计算电磁场问题的出发点,它既可写成微分形
式,又可写成积分形式。
微分形式的麦克斯韦方程组为
'H=J—ct
(1)
„cB
e=
(2)
ct
、B=0
(3)
(4)
式中,E为电场强度(V/m);B为磁感应强度(T);D为电位移矢量(C/m2);
H为磁场强度(A/m);J为电流密度(A/m2);'为电荷密度(C/m2)。
通常可将式
(1)称为麦克斯韦第一方程,将式
(2)称为麦克斯韦第二方程。
在麦克斯韦方程组中,有关场量之间的关系可表示为
D=E(5)
B-H(6)
J=注(7)
式中,;为介电常数(电容率);"为磁导率;二为电导率。
对于各向异性
媒质,这些参数是张量;对于各向同性媒质,它们是标量。
只有在线性且各向同性媒质的情况下,才是常数。
在SI单位制中,对应于自由空间的介电常数
-12IJ.它-7
0=8.85410F/m,磁导率o=4二10H/m。
积分形式麦克斯韦方程组为:
仆dl=i=JJdS+DdS(8)
•S£t'S
;EdlsBdS(9)
|:
tS
、sBdS=O(10)
;DdSdV(11)
SV
二、电磁场数值计算方法
1有限差分法
有限差分法是利用网格剖分将定解区域离散化为网格离散节点的集合,然
后,以差分原理为基础,以各离散点上函数的差商来近似替代该点的偏导数,把
要求解的边值问题转化为一组相应的差分方程问题,解出各离散点上的待求函数
值,即为所求定解问题的离散解,若再应用插值方法,便可从离散解得到定解问题在整个场域上的近似解。
1964年,Winslow利用向量位,采用有限差分离散,求解了二维非线性磁场问题。
优点:
网格剖分容易,数据准备省时,编制程序方便。
缺点:
对不规则的边界,如曲线边界,处理不方便。
当区域的边界线和内部媒介分界线形状比较复杂,以及场域的分布变化较大时,差分法的网格剖分缺少灵活性,给使用带来极大的不便。
有限差分法主要适用于边界形状规则的第一类边界,第二类齐次边界;静态场,时变场;线性场,非线性场等。
2有限元法
有限元法是根据变分原理和离散化而取得近似解的一种方法。
它首先从偏微分方程边值问题出发,找出一个能量泛函的积分式,并令其在满足第一类边界条件的前提下取极值,即构成条件变分问题。
然后,利用剖分插值,将变分问题离散化为普通多元函数的极值问题,解之即得待求边值问题的数值解。
3无网格Galerkin法
无网格Galerkin法与有限元法相似之处:
两者都是将边值问题等价为一个条件变分问题,然后由条件变分问题通过数值积分离散为代数方程组。
不同之处:
有限元法是对逐个有限单元进行数值积分,形成单元矩阵,然后将其叠加到单元节点所对应的方程中;而无网格Galerkin法是在积分单元上进行数值积分,然后将每个高斯点上的积分值叠加到该高斯点所支撑的若干节点所对应的方程中。
优点:
只需节点,不需单元,适合处理复杂边界问题,场函数的近似解连续可导,计算精度高,收敛速度快。
4小波分析算法目前,小波及其小波分析在电磁场工程问题中的应用已成为计算电磁学和工程界广泛关注的一个新的研究方向。
人们利用小波函数特有的消失矩、紧支集、正则性等性质,解决电磁场数值计算中的一些特殊问题。
小波分析在电磁场工程问题中的应用有两个方面:
小波变换与小波展开。
而小波变换又引申出小波包变换。
无论是小波变换还是小波包变换,都是为了使矩阵变得稀疏,便于方程求解。
小波展开的思想使用小波变换作为基函数,对电场或磁场作时域或空间域上的展开,并引申出一系列应用。
三、用ANSYS对永磁同步电机进行有限元分析
1永磁同步电动机空载及负载磁场的有限元析
样机结构尺寸如下:
极对数:
p=3;定子内、外径:
Di1=18018mm,D1=260mm;转子内、外径:
Di2=60mm,D2=180mm;定、转子槽数:
Q1=36,Q2=42;定子槽型选用通用圆底槽,转子槽选用圆顶圆底槽;转子采用气隙隔磁结构,永磁体内置、切间放置。
选取永磁同步电动机的一对磁极范围和定子外侧表面、电动机轴外表面作为计算区域,在Ansys中建模如图1所示。
图1样机的计算区域
1.1永磁同步电动机空载气隙磁场分析
对图1所示电动机的一对极区域进行求解,不加载荷,即定子电流全部赋零
仅永磁体作用。
通过后处理显示出磁场分析,即为样机的空载磁场,如图2所示
图2空载磁场磁力线分布图(一对极区域)
在Ansys后处理中,改变Ansys坐标为极坐标,以电动机轴心为圆心,经过最小气隙中心处取一圆弧,在这一圆弧上取出极坐标下的磁通密度Bx,即为气
隙径向磁通密度Br,把Br作为电动机空载气隙磁场的磁通密度数值。
以上述圆弧由电动机对称轴处到左右的长度为横坐标,Br为纵坐标,得到电动机空载磁场气隙磁密的空间分布曲线如图3所示。
图3空载气隙磁密分布曲线
以电动机轴心为圆心,经过定子内圆取一圆弧,在这一圆弧上取出磁位Az。
以上述圆弧由电动机对称轴处到左右的长度为横坐标,Az为纵坐标绘图,得到如图4所示的定子内圆节点磁位分布曲线。
图4定子内圆节点磁位分布曲线
1.2永磁同步电动机负载气隙磁场分析
永磁同步电动机负载时定子电流不为零,若选取A相轴线与时轴一致,并令
t=0时刻电动机的A相绕组轴线与转子q轴轴线重合,贝U三相电流瞬时值为:
iA=J2ICOSG=ImCOSOt
iB-21cos:
-120=Imcos:
-120
(1)
ic--21cosixT20=lmcosi很T20
由三相电流产生的电枢反应磁势与q轴有一固定角180-,如图5所示
图5电枢反应电势
由图5可以看出,如果假设:
=0,即iA=Im,iB二ic=Tm/2,此时的电枢反应磁势只有q轴分量。
在前面空载磁场分析基础上,再在图1所示电动机模型上加载如上三相定子电流密度,在Ansys中分析求解后,得到样机的交轴电枢反应磁场图如图6所示。
图7为交轴电枢反应气隙磁密分布图。
图6交轴电枢反应磁场图
图7交轴电枢反应气隙磁密图
如果:
=90,即iA=0,iB「3lm/2,ic~3lm/2,此时的电枢反应磁势
只有d轴分量。
在Ansys中分析求解后,得到样机的直轴电枢反应磁场图如图8
所示。
图9为直轴电枢反应气隙磁密分布图。
图8直轴电枢反应磁场图
DT5T
图9直轴电枢反应气隙磁密图
2永磁同步电动机空载反电动势计算
空载反电动势E0是永磁同步电动机的重要参数。
空载反电动势由电动机中永磁体产生的空载气隙磁密基波磁通在电枢绕组中感应产生。
由于永磁同步电动
机的励磁不能调节,无法像电励磁同步电动机通过调节励磁改变功率因数以达到改善电网功率因数的目的,因而必须合理选取空载反电动势。
合理设计E0,
可降低定子电流,提高电动机效率,降低永磁材料用量。
计算空载反电动势的方法很多,如磁链微分法、磁密分解法、磁位分解法等。
根据上述样机Ansys分析可知,经过有限元计算,可以求出样机气隙磁密Br和定子内圆节点矢量磁位Az的分布,分别如图3和图4所示。
由此可以得到计算永磁同步电动机空载反电动势的两种较为简便的方法。
2.1磁密分解法
在Ansys后处理中,不仅可以显示曲线图形,而且可以列出构成曲线的各点坐标值,根据这些坐标值,运用Ansys软件绘制气隙磁密Br波形,如图10波形a所示。
再对波形作博立叶分解,得到气隙磁密基波和各次谐波的幅值,取1、3、5次谐波,得出如式⑵所示的傅立叶级数展开式:
3x'i
—0.0302sin0.1233cos
lL丿
式中,L为半周期的圆弧长(一对极)。
画出气隙磁密基波波形如图10波形
图10气隙磁密波形图及基波波形
b所示。
计算结果如下:
气隙磁密基波幅值:
Bm二0.6057(T)
每极下空载基波气隙磁通幅值为:
2
mBmLef=0.0034(Wb)⑶
n
式中,Lef为电枢计算长度;Lef=9018mm;为极距,=2nx9012/2p=
94146mm。
从而得到空载励磁反电势为:
E。
二■2fNKdPm"0152fNKdp⑷
2.2磁位分解法
11
Z土汀Azdz二Az2.dz-Az1dz二Az2-Az1-Az⑸
00
对于二维平面场,通过单位轴向长度面积内的磁通量恰好等于这两点的标
量磁位之差。
对定子内圆上节点磁位Az作傅立叶分解,得到基波波形圆如图11波形b
所示。
傅立叶级数展开式为
其中,L为半周期的圆弧长(一对极)。
0,025
0B2
0.015
0.01
0.005
0
-Q.D05
-0.01
-0,015
-0.02
-0.025
图11定子内圆节点磁位分布曲线及其基波波形图
由于两点磁矢量位差的绝对值就是单位长度内两点之间的磁通量,所以上
式基波余弦项系数代表电动机单位长度每极q轴基波磁通的1/2,基波正弦项系
数代表电动机单位长度每极d轴基波磁通的1/2,则每极下空载基波气隙磁通幅值为:
m-2Lef-0.001320.0178^0.00324(Wb)(7)
可求得空载励磁反电动势为:
E。
二..2:
fNKdp爲=0.0144fNKdp(8)
通过对磁密分解法和磁位分解法计算得到的空载基波气隙磁通的计算结果分析对比,可见采用这两种方法得到的计算结果是相当的。
四、结语
由电机电磁场的数值计算,使人们对电机内部电磁场的分布和电磁性能了解得更透彻。
由于对电机深入研究的迫切需要,电磁场数值计算理论不断进步,
许多过去认识不清或无法分析的问题得到新的了解,许多概念更加深化,计算方法发
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 电磁场 数值 分析
![提示](https://static.bdocx.com/images/bang_tan.gif)