变压器铁芯磁路的计算

1、课程设计的目的与作用1.1、设计目的1、学习电机的工作原理及电机设计的相关方法，利用电机设计仿真软件AnsoftRMxprt2、参数设计法和利用MATLAB软件编程的传统设计方法完成典型电机产品设计；3、完成电机主要尺寸的选择和确定、基本性能设计、磁路计算、参数设计、起动计算等；4、培养学生文献检索的能力，特别是如何利用Internet检索需要的文献资料。5、培养学生综合分析问题、发现问题和解决问题的能力。6、培养学生运用知识的能力和工程设计的能力。7、提高学生课程设计报告撰写水平。1.2、设计作用课程设计是培养和锻炼在校学生综合应用所学理论知识解决实际问题能力、进行工程实训的重要教学环节，它具有动手、动脑，理论联系实际的特点，是培养在校工科大学生理论联系实际、敢于动手、善于动手和独立自主解决设计实践中遇到的各种问题能力的一种较好方法。《电机学》是电气工程及自动化专业的一门专业基础课，具有应用性、实践性较强的特点，忽视了实践环节，学生不能很好的理解所学内容。通过设计，使学生系统、深入了解各种电机的工作原理和抽象出来的数学模型，对这门课程的认识和理解提高到一个新的水平。通过设计实践，培养学生查阅专业资料、工具书或参考书，掌握现代设计手段和软件工具，并能以仿真程序及仿真结果表达其设计思想的能力。通过设计，不但要培养和提高学生学习和应用专业知识的能力，而且要在实践过程中锻炼培养正确的设计思想，培养良好的设计习惯，牢固树立事实求是和严肃认真的科学工作态度。电机学课程设计是电机学课程学习的最后一个环节，通过设计不仅可以使学生更牢固的掌握所学知识，同时也可以为后续课程的学习打下扎实的理论基础。2、设计任务及所用MATLAB软件、MAXWELL软件环境介绍2.1、MATLAB软件环境作为和Mathematica、Maple并列的三大数学软件。其强项就是其强大的矩阵计算以及仿真能力。要知道Matlab的由来就是Matrix+Laboratory=Matlab，所以这个软件在国内也被称作《矩阵实验室》。每次MathWorks发布Matlab的同时也会发布仿真工具Simulink。在欧美很多大公司在将产品投入实际使用之前都会进行仿真试验，他们所主要使用的仿真软件就是Simulink。Matlab提供了自己的编译器：全面兼容C++以及Fortran两大语言。所以Matlab是工程师，科研工作者手上最好的语言，最好的工具和环境。Matlab已经成为广大科研人员的最值得信赖的助手和朋友！目前MATLAB产品族可以用来进行：－数值分析－数值和符号计算－工程与科学绘图－控制系统的设计与方针－数字图像处理－数字信号处理－通讯系统设计与仿真－财务与金融工程...Simulink是基于MATLAB的框图设计环境，可以用来对各种动态系统进行建模、分析和仿真，它的建模范围广泛，可以针对任何能够用数学来描述的系统进行建模，例如航空航天动力学系统、卫星控制制导系统、通讯系统、船舶及汽车等等，其中了包括连续、离散，条件执行，事件驱动，单速率、多速率和混杂系统等等。Simulink提供了利用鼠标拖放的方法建立系统框图模型的图形界面，而且Simulink还提供了丰富的功能块以及不同的专业模块集合，利用Simulink几乎可以做到不书写一行代码完成整个动态系统的建模工作2.2、MAXELL软件环境MAXWELL是主要建立在MAXWELL方程基础上的，有限元分析软件。MAXWELL2D：工业应用中的电磁软件，如传感器、调节器、电动机、变压器，以及其他工业控制系统比以往任何时候都使用得更加广泛。由于设计者对性能与体积设计封装的希望，因而先进而便于使用的数字场仿真技术的需求也显著的增长。在工程人员所关心的实用性及数字化功能方面，MAXWELL的产品遥遥领先其他的一流公司。MAXWELL2D包括交流/直流磁场、静电场以及瞬态电磁场、温度场分析，参数化分级；以及优化功能。此外MAXWELL2D还可以产生高精度的等效电路模型以供Ansoft的SIMPLORER模块和其他电路分析工具调用。MAXWELL3D向导式的用户界面、精度驱动的自适应剖分技术和强大的后处理器时的MAXWELL3D成为业界最佳的高性能三维电磁技术软件。可以分析涡流、位移电流、集肤效应和邻近效应具有不可忽视作用的系统，得到电机、母线、变压器、线圈等电磁部件的整体特性。功率损耗、线圈损耗、某一频率下的阻抗（R和L）、力、转矩、电感、储能等参数可以自己计算。同时也可以给出整个相位的磁力线、B和H分布线、能量密度、温度分布等图形结果。RMxprt是基于磁路法的旋转电机专业设计软件，能加快电机的设计和优化过程。通过RMxprt用户能快速地对成百上千种设计方案进行评估，并可对预选方案进行优化设计。RMxprt是Ansoft电机及其驱动系统综合设计流程的理想起点，它既可以生成系统的模型也可以生成物理模型，通过他们用户可以结合电力电子技术和控制电路对初步设计加以优化。Rmxprt能分析的电机有：三项感应电动机，单相感应电动机，三相同步发电机和电动机，永磁直流无刷发电机，永磁同步电动机和发电机，永磁直流电动机，开关磁阻电动机，通用电动机，普通直流发电机和电动机，抓极交流发电机。3、设计题目3.1、变压器铁芯磁路的计算铁心磁路如图所示，磁路尺寸为：1=3mm、2=2mm、w=125mm、h=150mm、N1=2N2=100匝，铁芯宽度l=50mm、铁心厚度d=50mm。铁芯宽度和厚度均匀，忽略铁芯磁场边缘效应。1、假设铁芯的磁导率为无穷大，（1）若I2=0，Φ1=6mWb，求I1（2）若I1=10A，I2=20A，求Φ1和Φ2。2、若铁芯材料为DR510钢片（磁化曲线见教材），若I2=0，Φ1=6mWb，求I1和Φ1。3.2、同步发电机的特性分析3.2.1、4、设计过程4.1、变压器铁芯磁路的计算4.1.1、原理描述由磁路基尔霍夫第一定律=0，可得213，由于铁芯宽度和厚度均匀，截面积相等，则321BBB由磁路基尔霍夫第二定律mRHlNiF，得:lhHBlhwHINlhHBlhwHIN30222222301111113232μμ磁场关系:HBμ则方程为：iFeiHBBBBlhHBlhwHINlhHBlhwHINμμμ213302222223011111132324.1.2、参数计算1、假设铁芯的磁导率为无穷大，（1）若I2=0，Φ1=6mWb，求I1（2）若I1=10A，I2=20A，求Φ1和Φ。铁芯的磁导率为无穷大，则铁芯磁路H=0，方程化简为0222201111μμBINBIN（1）将I2=0，Φ1=6mWb带入上式方程，得I1=57.2958A（2）将I1=10A，I2=20A带入上式方程，得B1=0.418879T，B2=-0.628319T，则Φ1=B1ld=1.04720mWb，Φ2=B2ld=-1.57080mWb2、若铁芯材料为DR510钢片（磁化曲线见图），若I2=0，Φ1=6mWb，求I1和Φ1。Φ1=6mWb,则B1=Φ1/(ld)=2.4T采用编程法求解，由方程iFeiHBBBBlhHBlhwHINμμ21330222222320可看作未知数为B2的一元非线性方程，通过二分法求解出B2，再计算出B3，然后由方程lhHBlhwHIN3011111132μ计算出I1。为提高计算精度，可使用拉格朗日插值法或曲线拟合法，对表中不能直接得到的数据进行近似计算。本题使用拉格朗日插值法。拉格朗日插值法：假设任意两个不同的xj都互不相同，那么应用拉格朗日插值公式所得到的拉格朗日插值多项式为：其中每个为拉格朗日基本多项式（或称插值基函数），其表达式为：拉格朗日基本多项式的特点是在上取值为1，在其它的点上取值为0。图表1B-H磁化曲线图表2B-H磁化曲线4.1.3、机械特性图、MATLAB代码代码如下：#includeiostream#includecstdlib#includecmath#includeiomanip#includestdlib.husingnamespacestd;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;doublebh[2][150]={{0.4,doubleLagrange(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-delta2)+b2*delta2/u0-b2h(b3)*(h+l))/(-n2);coutendl计算过程信息：二分法计算B2次数:i_diedai_root1endl;coutendl磁密计算结果:endlsetprecision(16)B1=b1endlB2=b2endlB3=b3endl;coutendl磁通计算结果:endlsetprecision(16)Φ1=b1*l*dendlΦ2=b2*l*dendlΦ3=b3*l*d;coutendlI1=i1;coutendl结果验算:setprecision(16)Φ1=(b2+b3)*l*dI2=i2;system(pause);return0;}doubleb2h(doubleb)//计算磁场强度H,输入磁感应强度B，输出磁场强度H{intt;doubleh_t;t=abs((int)(100*b))-40;if(t134){if(b=0)returnlinearfit2(b);elsereturn-linearfit2(-b);}if(t15){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=bh[0];y=bh[1];z=((int)(100*xx))-45;if(xx=0.45)z=0;n=10+z;if(xx1.84)n=6+z;for(j=z,sum=0;jn;j++){for(i=z,pai=1;in;i++){if(i==j)continue;pai=pai*(xx-*(x+i))/(*(x+j)-*(x+i));}sum=sum+*(y+j)*pai;}returnsum;}doubleli