变压器铁芯磁路的计算要点

1、1、课程设计的目的与作用1.1、 设计目的1、学习电机的工作原理及电机设计的相关方法，利用电机设计仿真软件AnsoftRMxprt2、参数设计法和利用MATLA歆件编程的传统设计方法完成典型电机产品设计；3、完成电机主要尺寸的选择和确定、基本性能设计、磁路计算、参数设计、起动计算等；4、培养学生文献检索的能力，特别是如何利用Internet检索需要的文献资料。5、培养学生综合分析问题、发现问题和解决问题的能力。6、培养学生运用知识的能力和工程设计的能力。7、提高学生课程设计报告撰写水平。1.2、 设计作用课程设计是培养和锻炼在校学生综合应用所学理论知识解决实际问题能力、进行工程实训的重要教学环

2、节，它具有动手、动脑，理论联系实际的特点，是培养在校工科大学生理论联系实际、敢于动手、善于动手和独立自主解决设计实践中遇到的各种问题能力的一种较好方法。电机学是电气工程及自动化专业的一门专业基础课，具有应用性、实践性较强的特点，忽视了实践环节，学生不能很好的理解所学内容。通过设计，使学生系统、深入了解各种电机的工作原理和抽象出来的数学模型，对这门课程的认识和理解提高到一个新的水平。通过设计实践，培养学生查阅专业资料、工具书或参考书，掌握现代设计手段和软件工具，并能以仿真程序及仿真结果表达其设计思想的能力。通过设计，不但要培养和提高学生学习和应用专业知识的能力，而且要在实践过程中锻炼培养正确的设

3、计思想，培养良好的设计习惯，牢固树立事实求是和严肃认真的科学工作态度。电机学课程设计是电机学课程学习的最后一个环节，通过设计不仅可以使学生更牢固的掌握所学知识，同时也可以为后续课程的学习打下扎实的理论基础。2、设计任务及所用MATLA欧件、MAXWE附件环境介绍2.1、 MATLA颤件环境作为和Mathematica、Maple并列的三大数学软件。其强项就是其强大的矩阵计算以及仿真能力。要知道Matlab的由来就是Matrix+Laboratory=Matlab,所以这个软件在国内也被称作矩阵实验室。每次MathWorks发布Matlab的同时也会发布仿真工具Simulinko在欧美很多大公司

4、在将产品投入实际使用之前都会进行仿真试验，他们所主要使用的仿真软件就是Simulink0Matlab提供了自己的编译器：全面兼容C+©及Fortran两大语言。所以Matlab是工程师，科研工作者手上最好的语言，最好的工具和环境。Matlab已经成为广大科研人员的最值得信赖的助手和朋友！目前MATLAB产品族可以用来进行：- 数值分析- 数值和符号计算- 工程与科学绘图- 控制系统的设计与方针- 数字图像处理- 数字信号处理- 通讯系统设计与仿真- 财务与金融工程.Simulink是基于MATLAB的框图设计环境，可以用来对各种动态系统进行建模、分析和仿真，它的建模范围广泛，可以针对

5、任何能够用数学来描述的系统进行建模，例如航空航天动力学系统、卫星控制制导系统、通讯系统、船舶及汽车等等，其中了包括连续、离散，条件执行，事件驱动，单速率、多速率和混杂系统等等。Simulink提供了利用鼠标拖放的方法建立系统框图模型的图形界面，而且Simulink还提供了丰富的功能块以及不同的专业模块集合，利用Simulink几乎可以做到不书写一行代码完成整个动态系统的建模工作2.2、 MAXELL件环境MAXWEL是主要建立在MAXWELL程基础上的，有限元分析软件。MAXWELL2D工业应用中的电磁软件，如传感器、调节器、电动机、变压器，以及其他工业控制系统比以往任何时候都使用得更加广泛。

6、由于设计者对性能与体积设计封装的希望，因而先进而便于使用的数字场仿真技术的需求也显著的增长。在工程人员所关心的实用性及数字化功能方面，MAXWELL产品遥遥领先其他的一流公司。MAXWELL2包括交流/直流磁场、静电场以及瞬态电磁场、温度场分析，参数化分级；以及优化功能。此外MAXWELL2D可以产生高精度的等效电路模型以供Ansoft的SIMPLORE模块和其他电路分析工具调用。MAXWELL3晌导式的用户界面、精度驱动的自适应剖分技术和强大的后处理器时的MAXWEL3D成为业界最佳的高性能三维电磁技术软件。可以分析涡流、位移电流、集肤效应和邻近效应具有不可忽视作用的系统，得到电机、母线、变

7、压器、线圈等电磁部件的整体特性。功率损耗、线圈损耗、某一频率下的阻抗（R和L）、力、转矩、电感、储能等参数可以自己计算。同时也可以给出整个相位的磁力线、B和H分布线、能量密度、温度分布等图形结果。RMxprt是基于磁路法的旋转电机专业设计软件，能加快电机的设计和优化过程。通过RMxprt用户能快速地对成百上千种设计方案进行评估，并可对预选方案进行优化设计。RMxprt是Ansoft电机及其驱动系统综合设计流程的理想起点，它既可以生成系统的模型也可以生成物理模型，通过他们用户可以结合电力电子技术和控制电路对初步设计加以优化。Rmxprt能分析的电机有：三项感应电动机，单相感应电动机，三相同步发电

8、机和电动机，永磁直流无刷发电机，永磁同步电动机和发电机，永磁直流电动机，开关磁阻电动机，通用电动机，普通直流发电机和电动机，抓极交流发电机。3、设计题目3.1、 变压器铁芯磁路的计算铁心磁路如图所示，磁路尺寸为：61=3mm$2=2mmw=125mmh=150mmN=2N=100匝，铁芯宽度l=50mm铁心厚度d=50mm铁芯宽度和厚度均匀，忽略铁芯磁场边缘效应1、 假设铁芯的磁导率为无穷大，（1）若12=0,i=6mWb,求Ii（2）若Ii=10A,I2=20A,求1和2。2、 若铁芯材料为DR510钢片（磁化曲线见教材），若12=0,1=6mWb,求I1和1。3.2、 同步发电机的特性分析

9、3.2.1、4、设计过程4.1、 变压器铁芯磁路的计算4.1.1、 原理描述由磁路基尔霍夫第一定律Z*=0,可得则=曲-则，由于铁芯宽度和厚度均匀，截面积相等，则B1=B2-B3由磁路基尔霍夫第二定律工F=£Ni=£Hl=£4Rm,得:Ni|i=H1，2w+h+3l-6i+Bi+H3(h+l)J0°一N2|2=H22wh3l-、.2B2-H3hl°磁场关系：BH则方程为：Ni|i=Hi，2w+h+3-6i+Bi+H3(h+l)"°r162<_N2|2=H2，bw+h+3l62】+B2H3(h+l)°°

10、;B3=Bi-B2Bi=NFeHi4.1.2、 参数计算i、假设铁芯的磁导率为无穷大，(i)若I2=°,i=6mWb求Ii(2)若Ii=i°A,I2=2°A,求i和。铁芯的磁导率为无穷大，则铁芯磁路H=°,方程化简为c6iNi11=Bi°0°62N2I2=B2L.(I。(i)将12=。，i=6mW®入上式方程，得Ii=57.2958A(2)将Ii=i°A,12=2°A带入上式方程，得B=°.4i8879T,B=-°.6283i9T,则i=Bld=i.0472°mWh2=Bld

11、=-i.57°8°mWb2、若铁芯材料为DR5i硼片(磁化曲线见图)，若I2=。，i=6mWb求Ii和i。i=6mWb!UB=i/(ld)=2.4T采用编程法求解，由方程N2|2=0=H2.bw+h+3l&2+B2«H3（h+l）0°/B3=B1-B2Bi=FFeHi可看作未知数为B的一元非线性方程，通过二分法求解出B2,再计算出国然后由方程Ni|i=Hi2wh31-、11BiH3hl（10计算出I1。为提高计算精度，可使用拉格朗日插值法或曲线拟合法，对表中不能直接得到的数据进行近似计算。本题使用拉格朗日插值法。拉格朗日插值法:假设任意两个不同的

12、Xj都互不相同，那么应用拉格朗日插值公式所得到的拉格朗日插值多项式为：k£3）:二£%匕（0j=0其中每个（f）为拉格朗日基本多项式（或称插值基函数），其表达式为：,f.yr£一与（工一10）（一叼-1）（1（£1矗）彳印,*当有一与（叼一Q）（叼叼-1）（叼叼+1）（有一了的）拉格朗日基本多项式2K1）的特点是在叼上取值为1,在其它的点踢，*于上取值为00xW图表1B-H磁化曲线图表2B-H磁化曲线4.1.3、 机械特性图、MATLA跃码代码如下：#include<iostream>#include<cstdlib>#inclu

13、de<cmath>#include<iomanip>#include<stdlib.h>usingnamespacestd;doubledelta1=0.003,delta2=0.002,w=0.125,h=0.15,n1=100,n2=50,l=0.05,d=0.05;doublepi=3.1415926535898;doubleu0=4*pi*pow(10,-7);doublei1,i2,f10,f20;inti_diedai_root1=0,i_diedai_root2=0,i_ddjs=0;doublebh2150=0.4,doubleLagrang

14、e(doublexx);doubleb2h(doubleb);doublef2(doubleb2,doubleb1);doubleroot2(doublex1,doublex2,doubleb1);doublelinearfit2(doublex);doublelinearfit1(doublex);intmain()doubleb1,b2,b3,h1,h2,h3,f1,f2;b1=2.4;b2=root2(0,2.4,b1);b3=b1-b2;i1=(h1*(2*w+h+3*l-delta1)+b1*delta1/u0+h3*(h+l)/n1;i2=(b2h(b2)*(2*w+h+3*l-d

15、elta2)+b2*delta2/u0-b2h(b3)*(h+l)/(-n2);cout<<endl<<"计算过程信息："<<"二分法计算B2次数:"<<i_diedai_root1<<endl;cout<<endl<<"磁密计算结果："<<endl<<setprecision(16)<<"B1="<<b1<<endl<<"B2="<&

16、lt;b2<<endl<<"B3="<<b3<<endl;cout<<endl<<"磁通计算结果:"<<endl<<setprecision(16)<<"1="<<b1*l*d<<endl<<"2="<<b2*l*d<<endl<<"3="<<b3*l*d;cout<<endl<<&

17、quot;I1="<<i1;cout<<endl<<"结果验算："<<setprecision(16)<<"1="<<(b2+b3)*l*d<<”I2="<<i2;system("pause");return0;doubleb2h(doubleb)/计算磁场强度H,输入磁感应强度B,输出磁场强度Hintt;doubleh_t;t=abs(int)(100*b)-40;if(t>134)if(b>=0)retu

18、rnlinearfit2(b);elsereturn-linearfit2(-b);if(t<15)if(b>=0)returnlinearfit1(b);elsereturn-linearfit1(-b);if(b>=0)returnLagrange(b);elsereturn-Lagrange(-b);doubleLagrange(doublexx)/拉格朗日插值，输入自变量，输出插值函数值inti=0,j=0,n=150,z=0;doublepai,sum;double*x,*y;x=bh0;y=bh1;z=(int)(100*xx)-45;if(xx<=0.45

19、)z=0;n=10+z;if(xx>1.84)n=6+z;for(j=z,sum=0;j<n;j+)(for(i=z,pai=1;i<n;i+)(if(i=j)continue;pai=pai*(xx-*(x+i)/(*(x+j)-*(x+i);sum=sum+*(y+j)*pai;returnsum;doublelinearfit2(doublex)(return(56450*x-89188);doublelinearfit1(doublex)(return(202.21*x+57.015);doublef2(doubleb2,doubleb1)(doubleb3,h2,h

20、3;b3=b1-b2;h2=b2h(b2);h3=b2h(b3);returnh2*(2*w+h+3*l-delta2)+b2*delta2/u0-h3*(h+l)+n2*i2;doubleroot2(doublex1,doublex2,doubleb1)(doublex,y,y1;y1=f2(x1,b1);do(cout<<'!'i_diedai_root1+;if(i_diedai_root1>=1000000)break;x=(x1+x2)/2;y=f2(x,b1);if(y*y1>0)(x1=x;yi=y;)elsex2=x;while(fabs

21、(y)>=2*pow(10,-13);returnx;4.2、5、结果分析5.1、 变压器铁芯磁路的计算计算结果（设定收敛精度10A-13）:二二7sersAdmi川市木口八Dhsiktop、电机篌铁芯2-利亚binDeb<j,exe»?!?!?!?!?!计算过程信息二分法计算股次数W4磁密计算结果：Bl=2,4B2-0,721?4862494B2798B3-1.67825137505172磁通计算结果：力1=0.006中2=0,30180437156237073-3,fl(1419562B437629301(11-57.29577951309219结果验算二力1=0.0

22、的12=5.1超8嘱46722646妊-015请按任意键继续.计算结果：B1=2.4T,B2=0.721749T,B3=1.67825T1=0.006Wb2=0.0018044Wb3=0.00419563WbI1=57.2958A5.2、6、设计总结此次课设内容是铁心磁路计算，只牵涉到电机学第一章的内容，所以难度并不是很大。但即便如此，既然是选择这门课程设计，就必须与组员合作，认真努力地完成，不能图完成任务，而且实际上也碰到了一些问题。我负责的是第二问，编程的时候所得出的磁感应强度超过了书上所给表格的对应范围，当时还特意问了熊老师，老师说这个需要自己解决，并建议了我用线性插值法。回来的时候，和张丰伟一起讨论了一下，便将磁化曲线进行了拟合，还运用了拉格朗日插值法（这个是从电路实验上看到的，记得当时有提过，线性插值法精度好像不够），但是一开始n设置的过大，导致了无法收敛，在网上查了原因，解释是当插值点比较多的时候，拉格朗日插值多项