电机学-研究学习报告-三相感应电动机机械特性曲线的仿真分析 精品.docx

电机学研究性学习三相感应电动机机械特性曲线的仿真分析一MATLAB的介绍：MATLAB是美国MathWorks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB和Simulink两大部分。MATLAB是matrix&laboratory两个词的组合，意为矩阵工厂（矩阵实验室）。是由美国公司MathWorks发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。MATLAB的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来解算问题要比用C，FORTRAN等语言完成相同的事情简捷得多，并且MATLAB也吸收了像Maple等软件的优点，使MATLAB成为一个强大的数学软件。在新的版本中也加入了对C，FORTRAN，C+，JAVA的支持。可以直接调用，用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用，此外许多的MATLAB爱好者都编写了一些经典的程序，用户可以直接进行下载就可以用。二训练载体研究目的本次研究性训练的目的是通过利用MATLAB软件或C语言编程，仿真计算一台实际的三相感应电动机稳态运行时的机械特性曲线，并分析影响感应电动机机械特性的因素。在三相感应电机的稳态运行中，当负载变化时，其电磁转矩与电机转速之间的曲线称为感应电动机的机械特性，机械特性曲线的不同形状代表电机的性能差异，三相感应电动机机械特性的公式为：可见，当上述公式中的各参数发生变化时，将直接影响电机的机械特性曲线形状的改变。本次专题中需要应用MATLAB软件编程或C语言编程，MATLAB软件具有强大的绘图功能，能将计算的结果进行直观的图像显示，能方便对计算结果正确与否进行直观的判断。通过本专题可以使学生把所学的MATLAB软件应用到实际电机的计算、绘图中，用直观的图形变化来认识抽象的电机理论。三训练载体研究内容1 利用三相感应电动机的机械特性公式，进行MATLAB 编程，计算并绘制当电机转速变化时，电机的机械特性曲线。2 改变电源电压，计算且绘制当电机转速变化时，电机的机械特性曲线。并比较特性曲线的变化特点，有何结论？3 改变电源频率，计算且绘制当电机转速变化时，电机的机械特性曲线。并比较特性曲线的变化特点，有何结论？4 改变转子电阻，计算且绘制当电机转速变化时，电机的机械特性曲线。并比较特性曲线的变化特点，有何结论？5 改变定子电阻，计算且绘制当电机转速变化时，电机的机械特性曲线。并比较特性曲线的变化特点，有何结论？6 改变定子电抗，计算且绘制当电机转速变化时，电机的机械特性曲线。并比较特性曲线的变化特点，有何结论？7 改变转子电抗，计算且绘制当电机转速变化时，电机的机械特性曲线。并比较特性曲线的变化特点，有何结论？8 对影响三相感应电动机机械特性曲线的变化因素及其曲线的变化特点进行总结。谈谈本次研究性训练的收获。四、训练载体题目要求要求学生参考课本第五章“Induction Machine”之中三相感应电动机机械特性公式及其曲线的形状等，通过MATLAB编程计算，绘制出电机参数改变时，电机的具体机械特性曲线，理解和掌握影响电机机械特性曲线变化的因素。通过本次训练，能够让学生对三相感应电动机机械特性曲线的特点加深理解，加深对电机后续内容的理解和掌握。五、MATLAB仿真数据及仿真结果截图（一）绘制如下参数的三相四极异步电动机各运行特性：Pn=7.5kW，U1n=380V，f1n=50Hz，Nn=1455r/min，R1=1.42，X1=5.41，R2=1.25，X2 =7.75，Xm=90.4，Rm=10.6 ，T0=94.1Nm （Y接法）三相异步电动机的运行特性:1.转速（N）与P2的关系:（1）MATLAB文本：clcclearsymsU1nNphPolesNnR1R2pX10X20pRmXm0X1X2pT0.Ns0Z1Z2ZmImTem1iI1I2cudE1I2sn; U1n=380/sqrt(3);Nph=3;Poles=4;Fe0=50; Nn=1455;R1=1.42;R2p=1.25;X10=5.41; X20p=7.75;Rm=10.6;Xm0=90.4; Ns0=120*Fe0/Poles; Z1=R1+j*X10; Zm=Rm+j*Xm0;Z2=R2p+j*X20p;c=1+X10/Xm0;Im=U1n/(Z1+Zm);Os=4*pi*Fe0/Poles; fori=1:1:7500 Pr1=iI2=(Pr1/(Z1+c*Z2)(0.5);d=angle(-Im+I2);I1=abs(-Im+I2);I2s=abs(I2);E1=abs(-Im*Zm)Pe=Nph*E1*I2s*cos(d)s=Nph*(I2s)(2)*R2p/Pe; n=Ns0-s*Ns0;plot(Pr1,n); holdon; end（2）MATLAB仿真结果截图：2.电流（I1）与P2的关系： （1）MATLAB文本：clcclearsymsU1nNphPolesNnR1R2pX10X20pRmXm0X1X2pT0.Ns0ZeqZeq1Zeq2Z1Z2ZmI1niF1I1c; U1n=380/sqrt(3);Nph=3;Poles=4;Fe0=50; Nn=1455;R1=1.42;R2p=1.25;X10=5.41; X20p=7.75;Rm=10.6;Xm0=90.4; Ns0=120*Fe0/Poles; Sn=(Ns0-Nn)/Ns0; Z1=R1+j*X10; Zm=Rm+j*Xm0;Z2=R2p+j*X20p;c=1+X10/Xm0;Zeq=abs(Z1+Zm);Zeq1=abs(Z1+c*Z2);Im=U1n/Zeq;Os=4*pi*Fe0/Poles;fori=1:1:7500 Pr1=iI2=(Pr1/Zeq1)(0.5);I1=-Im+I2;plot(Pr1,I1); holdon; end（2）MATLAB仿真结果截图：3.效率（）与P2的关系：（1）MATLAB文本：clcclearsymsU1nNphPolesNnR1R2pX10X20pRmXm0X1X2pT0.Ns0XmZ1Z2ZmiF1I1cud; U1n=380/sqrt(3);Nph=3;Poles=4;Fe0=50; Nn=1455;R1=1.42;R2p=1.25;X10=5.41; X20p=7.75;Rm=10.6;Xm0=90.4; Ns0=120*Fe0/Poles; Sn=(Ns0-Nn)/Ns0; Z1=R1+j*X10; Zm=Rm+j*Xm0;Z2=R2p+j*X20p;c=1+X10/Xm0;Im=U1n/(Z1+Zm);Os=4*pi*Fe0/Poles;fori=1:1:7500 Pr1=iI2=(Pr1/(Z1+c*Z2)(0.5);d=angle(-Im+I2);I1=abs(-Im+I2);P1=3*U1n*I1*cos(d)u=Pr1/P1plot(Pr1,u); holdon; end（2）MATLAB仿真结果截图：4.功率因数(cos)与P2的关系：（1）MATLAB文本：clcclearsymsU1nNphPolesNnR1R2pX10X20pRmXm0X1X2pT0.Ns0XmZ1Z2ZmiF1I1ccos; U1n=380/sqrt(3);Nph=3;Poles=4;Fe0=50; Nn=1455;R1=1.42;R2p=1.25;X10=5.41; X20p=7.75;Rm=10.6;Xm0=90.4; Ns0=120*Fe0/Poles; Sn=(Ns0-Nn)/Ns0; Z1=R1+j*X10; Zm=Rm+j*Xm0;Z2=R2p+j*X20p;c=1+X10/Xm0;Im=U1n/(Z1+Zm);Os=4*pi*Fe0/Poles;fori=1:1:7500 Pr1=iI2=(Pr1/(Z1+c*Z2)(0.5);cos=angle(Im+I2);plot(Pr1,cos); holdon; end（2）MATLAB仿真结果截图：（3）MATLAB文本：clcclearsymsU1nNphPolesNnR1R2pX10X20pRmXm0X1X2pT0.Ns0XmZ1Z2ZmiF1I1ccos; U1n=380/sqrt(3);Nph=3;Poles=4;Fe0=50; Nn=1455;R1=1.42;R2p=1.25;X10=5.41; X20p=7.75;Rm=10.6;Xm0=90.4; Ns0=120*Fe0/Poles; Sn=(Ns0-Nn)/Ns0; Z1=R1+j*X10; Zm=Rm+j*Xm0;Z2=R2p+j*X20p;c=1+X10/Xm0;Im=U1n/(Z1+Zm);Os=4*pi*Fe0/Poles;fori=1:1:7500 Pr1=iI2=(Pr1/(Z1+c*Z2)(0.5);cos=angle(-Im+I2);plot(Pr1,cos); holdon; end（4）MATLAB仿真结果截图：5.电磁转矩（Te）与P2的关系：（1）MATLAB文本：clcclearsymsU1nNphPolesNnR1R2pX10X20pRmXm0X1X2pT0.Ns0XmZ1Z2ZmI1niF1I1cudE1I2s; U1n=380/sqrt(3);Nph=3;Poles=4;Fe0=50; Nn=1455;R1=1.42;R2p=1.25;X10=5.41; X20p=7.75;Rm=10.6;Xm0=90.4; Ns0=120*Fe0/Poles; Sn=(Ns0-Nn)/Ns0; Z1=R1+j*X10; Zm=Rm+j*Xm0;Z2=R2p+j*X20p;c=1+X10/Xm0;Im=U1n/(Z1+Zm);Os=4*pi*Fe0/Poles; fori=1:1:7500 Pr1=iI2=(Pr1/(Z1+c*Z2)(0.5);d=angle(Im+I2);I1=abs(Im+I2);I2s=abs(I2);E1=abs(-Im*Zm)Pe=Nph*E1*I2s*cos(d)Te=Pe/Os; plot(Pr1,Te); holdon; end（2）MATLAB仿真结果截图：如下参数的三相四极异步电动机机械特性的绘制： U1n=380V，f1n=50Hz，Nn=1455r/min，R1=1.42，X1=5.41，R2=1.25，X2 =7.75，Xm=90.4，Rm=10.6（Y接法）(二)改变电源电压，绘制当电机转速变化时电机的机械特性曲线。（1）MATLAB文本：clcclearsymsU1nNphPolesFe0Fe1Fe2Fe3Fe4NnR1R2pX10X20pRmXm0X1X2p.Ns0XmZeq1Z1I1nSnE1NsmSNr1Tem1iF1; U1n=380/sqrt(3);Nph=3;Poles=4;Fe0=50; Nn=1455;R1=1.42;R2p=1.25;X10=5.41; X20p=7.75;Rm=10.6;Xm0=90.4; Ns0=120*Fe0/Poles; Sn=(Ns0-Nn)/Ns0; Z1=R1+j*X10; Zeq1=(Rm+j*Xm0)*(R2p/Sn+X20p)/(Rm+j*Xm0)+(R2p/Sn+j*X20p); I1n=abs(Zeq1/(Zeq1+Z1)*U1n); Ue1=380;Ue2=400;Ue3=420;Ue4=440; for m=1:4 if m=1U1=Ue1/sqrt(3); elseif m=2 U1=Ue2/sqrt(3); elseif m=3 U1=Ue3/sqrt(3); elseif m=4 U1=Ue4/sqrt(3); endfori=1:2000 S=i/2000; Nr1=Ns0*(1-S); Tem1=Nph*Poles/(4*pi)*(U1/Fe0)2*(Fe0*R2p/S/(R1+(R2p/S)2+(X10+X20p)2); plot(Tem1,Nr1); holdon; endend（2）MATLAB仿真结果截图：（三）采用恒U1=f1控制，通过MATLAB编程，绘出下列不同供电频率(f1=50，60，70，80Hz）下三相异步电动机的机械特性：（1）MATLAB文本：clcclearsymsU1nNphPolesFe0Fe1Fe2Fe3Fe4NnR1R2pX10X20pRmXm0X1X2p.Ns0XmZeq1Z1I1nSnE1NsmSNr1Tem1iF1; U1n=380/sqrt(3);Nph=3;Poles=4;Fe0=50; Nn=1455;R1=1.42;R2p=1.25;X10=5.41; X20p=7.75;Rm=10.6;Xm0=90.4; Ns0=120*Fe0/Poles; Sn=(Ns0-Nn)/Ns0; Z1=R1+j*X10; Zeq1=(Rm+j*Xm0)*(R2p/Sn+X20p)/(Rm+j*Xm0)+(R2p/Sn+j*X20p); I1n=abs(Zeq1/(Zeq1+Z1)*U1n); Fe1=50;Fe2=60;Fe3=70;Fe4=80; for m=1:4 if m=1F1=Fe1; elseif m=2 F1=Fe2; elseif m=3 F1=Fe3; elseif m=4 F1=Fe4; endNs=120*F1/Poles; X1=X10*(F1/Fe0); X2p=X20p*(F1/Fe0); Xm=Xm0*(F1/Fe0); U1=U1n; fori=1:2000 S=i/2000; Nr1=Ns*(1-S); Tem1=Nph*Poles/(4*pi)*(U1/Fe0)2*(Fe0*R2p/S/(R1+(R2p/S)2+(X10+X20p)2); plot(Tem1,Nr1); holdon; endend（2）MATLAB仿真结果截图：（四）改变转子电阻，绘制当电机转速变化时电机的机械特性曲线。（1）MATLAB文本：clcclearsymsU1nNphPolesFe0Fe1Fe2Fe3Fe4NnR1R2pX10X20pRmXm0X1X2p.Ns0XmZeq1Z1I1nSnE1NsmSNr1Tem1iF1; U1n=380/sqrt(3);Nph=3;Poles=4;Fe0=50; Nn=1455;R1=1.42;R2p=1.25;X10=5.41; X20p=7.75;Rm=10.6;Xm0=90.4; Ns0=120*Fe0/Poles; Sn=(Ns0-Nn)/Ns0; Z1=R1+j*X10; Zeq1=(Rm+j*Xm0)*(R2p/Sn+X20p)/(Rm+j*Xm0)+(R2p/Sn+j*X20p); I1n=abs(Zeq1/(Zeq1+Z1)*U1n); Re1=1.25;Re2=1.5;Re3=1.75;Re4=2; for m=1:4 if m=1R2p=Re1; elseif m=2 R2p=Re2; elseif m=3 R2p=Re3; elseif m=4 R2p=Re4; endU1=U1nfori=1:2000 S=i/2000; Nr1=Ns0*(1-S); Tem1=Nph*Poles/(4*pi)*(U1/Fe0)2*(Fe0*R2p/S/(R1+(R2p/S)2+(X10+X20p)2); plot(Tem1,Nr1); holdon; endend（2）计算结果截图：（3）MATLAB仿真结果截图：（五）改变定子电阻，绘制当电机转速变化时电机的机械特性曲线。（1）MATLAB文本：clcclearsymsU1nNphPolesFe0Fe1Fe2Fe3Fe4NnR1R2pX10X20pRmXm0X1X2p.Ns0XmZeq1Z1I1nSnE1NsmSNr1Tem1iF1; U1n=380/sqrt(3);Nph=3;Poles=4;Fe0=50; Nn=1455;R1=1.42;R2p=1.25;X10=5.41; X20p=7.75;Rm=10.6;Xm0=90.4; Ns0=120*Fe0/Poles; Sn=(Ns0-Nn)/Ns0; Z1=R1+j*X10; Zeq1=(Rm+j*Xm0)*(R2p/Sn+X20p)/(Rm+j*Xm0)+(R2p/Sn+j*X20p); I1n=abs(Zeq1/(Zeq1+Z1)*U1n); Re1=1.5;Re2=2;Re3=2.5;Re4=3; for m=1:4 if m=1R1=Re1; elseif m=2 R1=Re2; elseif m=3 R1=Re3; elseif m=4 R1=Re4; endU1=U1nfori=1:2000 S=i/2000; Nr1=Ns0*(1-S); Tem1=Nph*Poles/(4*pi)*(U1/Fe0)2*(Fe0*R2p/S/(R1+(R2p/S)2+(X10+X20p)2); plot(Tem1,Nr1); holdon; endend（2）MATLAB仿真结果截图：（六）改变定子电抗，计算且绘制当电机转速变化时，电机的机械特性曲线。（1）MATLAB文本：clcclearsymsU1nNphPolesFe0Fe1Fe2Fe3Fe4NnR1R2pX10X20pRmXm0X1X2p.Ns0XmZeq1Z1I1nSnE1NsmSNr1Tem1iF1; U1n=380/sqrt(3);Nph=3;Poles=4;Fe0=50; Nn=1455;R1=1.42;R2p=1.25;X10=5.41; X20p=7.75;Rm=10.6;Xm0=90.4; Ns0=120*Fe0/Poles; Sn=(Ns0-Nn)/Ns0; Z1=R1+j*X10; Zeq1=(Rm+j*Xm0)*(R2p/Sn+X20p)/(Rm+j*Xm0)+(R2p/Sn+j*X20p); I1n=abs(Zeq1/(Zeq1+Z1)*U1n); Xe1=5.41;Xe2=6.41;Xe3=7.41;Xe4=8.41; for m=1:4 if m=1X10=Xe1; elseif m=2 X10=Xe2; elseif m=3 X10=Xe3; elseif m=4 X10=Xe4; endU1=U1nfori=1:2000 S=i/2000; Nr1=Ns0*(1-S); Tem1=Nph*Poles/(4*pi)*(U1/Fe0)2*(Fe0*R2p/S/(R1+(R2p/S)2+(X10+X20p)2); plot(Tem1,Nr1); holdon; endend（2）MATLAB仿真结果截图：（七）改变转子电抗，绘制当电机转速变化时电机的机械特性曲线。（1）MATLAB文本：clcclearsymsU1nNphPolesFe0Fe1Fe2Fe3Fe4Nn