变压器铁芯磁路的计算

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

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

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

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

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

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

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

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

9、.2、同步发电机的特性分析3.2.1、4、设计过程 4.1、变压器铁芯磁路的计算4.1.1、原理描述由磁路基尔霍夫第一定律=0,可得312,由于铁芯宽度和厚度均匀，截面积相等，则B1 B2 B3由磁路基尔霍夫第二定律NiHIRm，得:N1| 1 H1 2w3I1B1 口 0H3 hN2|2 H2 2W3l22 B2H 3 h l口0磁场关系：B卩H则方程为:Ni|i Hi 2w hN2|2 H2 2w13l 1 B1 H3 h l卩02h 3l 2B2H 3 h l口0B3 B1 B2Bi 口 FeHi4.1.2、参数计算1、假设铁芯的磁导率为无穷大,(1)若 12=0, 1=6mWb求 11

10、 (2)若 I1=10A, I2=2OA,求1和。铁芯的磁导率为无穷大，则铁芯磁路 H=0,方程化简为1N111 B1 卩02N2I 2 B2 口0(1)将|2=0,1=6mW带入上式方程，得I 1=57.2958A(2)将 I1=1OA, I2=2OA 带入上式方程，得B=O.418879T, B=-O.628319T，则1=Bld=1.04720mWb 2=Bld=-1.57080mWb2、若铁芯材料为DR510钢片（磁化曲线见图），若12=0,1=6mWb求11和1。 i=6mWI则 Bi=O i/(ld)=2.4T采用编程法求解，由方程N2I2 0 H2 2w h 3122B2H3 h

11、 lu0B3BiB2Bi卩 FeHi可看作未知数为E2的一元非线性方程,通过二分法求解出B2,再计算出B3,然后由方程1N1I1 H1 2w h 3l 1 B1 H3h lu0计算出11。为提高计算精度，可使用拉格朗日插值法或曲线拟合法，对表中不能直接得到的数据 进行近似计算。本题使用拉格朗日插值法。拉格朗日插值法：假设任意两个不同的xj都互不相同，那么应用拉格朗日插值公式所得到的拉格朗日插值多项式为：j=0（広-珂）其中每个&（工）为拉格朗日基本多项式（或称插值基函数），其表达式为：fl f富一网 _（H 国0）（忑 Nj_i）（常一叼+）'1=0,Xj 毎 （巧©

12、0）（叼弓-1）（巧叼+1）（叼业）拉格朗日基本多项式与（工）的特点是在巧上取值为1,在其它的点坯'¥ j上取值为0。图表1 B-H磁化曲线4.1.3、机械特性图、MATLA代码代码如下:#in elude <iostream>#in elude <cstdlib>#in elude <emath>#i nclude <ioma nip> #i nclude <stdlib.h> using n ames pace std;doubledelta1=0.003,delta2=0.002,w=0.125,h=0.15,

13、 n1=100, n2=50,l=0.05,d=0.05; double pi=3.1415926535898;double u0=4* pi* po w(10,-7);double i1,i2,f10,f20;int i_diedai_root1=0,i_diedai_root2=0,i_ddjs=0;double bh2150=0.4, double Lagra nge(double xx);double b2h(double b);double f2(double b2,double b1);double root2(double x1,double x2,double b1);doub

14、le lin earfit2(double x);double lin earfit1(double x);int mai n()double b1,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);分法计算 B2次coutvvendlvv"计算过程信息："<<"

15、1 =数:"vvi_diedai_root1vve ndl;coutvve ndlvv"磁密计算结果:"<<e ndlwset precisio n(16)vv"B仁"<<b1<<e ndlvv"B2= "<<b2<<e ndlvv"B3= "<<b3<<e ndl ;coutvvendlvv"磁通计算结果:"<<endlwsetprecision(16)vv""<

16、<b1*l*d<<e ndlvv" 2= "vvb2*l*dvve ndlvv" 3= "vvb3*l*d ;coutvve ndlvv"l仁"vvi1;coutvvendlvv"结 果验算： "vvsetprecision(16)vv" 1="vv(b2+b3)*l*dvv"I2= "vvi2 ;system(" pause");return 0;double b2h(double b) / 计算磁场强度H,输入磁感应强度B,输出磁场强

17、度 H int t;double h_t;t=abs(i nt)(100*b)-40;if(t>134)if(b>=0)return lin earfit2(b);else return -li nearfit2(-b);if(tv15)if(b>=0)return lin earfit1(b); else return -li nearfitl(-b);if(b>=0) return Lagra nge(b);else return -Lagra nge(-b); double Lagra nge(double xx) /拉格朗日插值，输入自变量，输出插值函数值int

18、 i=0,j=0 ,n=150,z=0;double p ai,sum;double *x,*y;x=bh0;y=bh1;z=(i nt)(100*xx)-45; if(xxv=0.45) z=0; n=10+z;if(xx>1.84) n=6+z;for(j=z,sum=0;j vn ;j+)for(i=z ,p ai=1;i vn ;i+)if(i=j) continue;p ai=pai*(xx-*(x+i)/(*(x+j)-*(x+i);sum=sum+*(y+j)* pai;return sum;double lin eafit2(double x)return (56450*

19、x-89188);double lin earfit1(double x)return (202.21*x+57.015);double f2(double b2,double b1)double b3,h2,h3;b3=b1-b2;h2=b2h(b2);h3=b2h(b3);return h2*(2*w+h+3*l-delta2)+b2*delta2/u0-h3*(h+l)+n2*i2;double root2(double x1,double x2,double b1) double x,y,y1;y仁f2(x1,b1);docoutvv'!'i_diedai_root1+

20、; 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(y)>=2* pow(10,-13); return x;4.2、5、结果分析5.1、变压器铁芯磁路的计算计算结果（设定收敛精度10-13）:磁由计箕结黑：E1- 2,4E2- 0_7217486249482798E3= i.G792G137&e&172磁通1+算结果：计算过程信息，二分法计算B2次数洱九C>2= e.0018043715623707eg A

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