资源描述:
《基于MATLAB的异步电动机仿真.doc》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、课程报告COURSEREPORT课程名称:——课程号:——授课教师:——学号:——姓名:西木小卒所属:上大桂树林基于MATLAB的异步电动机仿真1.电机仿真模型一台三相六机鼠笼式异步电动机,定子绕组为Y型连接,额定电压为UN=380V,额定转速nN=975r/min,电源频率f1=50Hz,定子电阻R1=2.08Ω,定子漏电抗X1=3.12Ω,转子电阻折合值=1.53Ω,转子漏电抗折合值=4.25Ω。要求:绘制以上参数电动机的固有机械特性曲线、定子串电阻人为特性曲线、电子串电抗人为特性曲线、转子串电阻人为特性曲线、降电压人为特性曲线;给出仿真源代码。2
2、.仿真代码实现clcclearm1=3;%%电机相数U1=220;%%额定定子相电压n1=[-1000,1000];%%两个相的同步转速,+为规定正相,-为反相p=3;%%电机极对数f1=50;%%输入电流频率r1=2.08;%%定子侧电阻r2=1.53;%%转子侧电阻折合值w1=2*pi*f1/p;%%电机同步角速度x1=3.12;%%定子侧电抗x2=4.25;%%转子侧折合电抗s=-1:0.001:1;%%设定转差率,是画图的关键自变量ns=[-1,1];%用来标定转矩的方向,规定ns=1为正,-1为反%%绘制定子串电阻前的固有机械特性曲线figu
3、re(1);fori=1:length(n1);n=n1(i)*(1-s);%%计算转速T=ns(i)*(m1*p*U1^2*r2)./s./(w1.*((r1+r2./s).^2+(x1+x2)^2));%%计算转矩plot(T,n,'--');ifi==1text(max(T),800,strcat('r1=',num2str(r1),'Omega'),'FontSize',9,'Color','black');%%曲线标注位置设置,每条曲线的标注横轴上对齐其Tm以便于区分,纵坐标无严格限制endholdon;endxlabel('电磁转矩{i
4、tT}/(Ncdotm)');ylabel('转速{itn}/(r/min)');title('定子串电阻的机械特性曲线')gridon;holdon;%%绘制定子串电阻人为机械特性曲线r1p=r1;str_y=800;fori=1:length(n1);%%绘制正反相序的波形图,n=n1(i).*(1-s);forcoef=1:3%%三条人为曲线r1p=r1p+1;%%串电阻的步进值T1=ns(i)*(m1*p*U1^2*r2)./s./(w1.*((r1p+r2./s).^2+(x1+x2)^2));plot(T1,n,'k-');holdon
5、;ifi==1str=strcat('r1=',num2str(int16(r1p)),'Omega');str_y=str_y-150;text(max(T1),str_y,strcat('r1=',num2str(r1p),'Omega'),'FontSize',9,'Color','black');endendend%%绘制降电压前的固有机械特性曲线figure(2);fori=1:length(n1);n=n1(i)*(1-s);T=ns(i)*(m1*p*U1^2*r2)./s./(w1.*((r1+r2./s).^2+(x1+x2)^2
6、));plot(T,n,'--');ifi==1text(max(T),800,strcat('U1=',num2str(int16(U1)),'V'),'FontSize',9,'Color','black');endholdon;endxlabel('电磁转矩{itT}/(Ncdotm)');ylabel('转速{itn}/(r/min)');title('降低定子电压的机械特性曲线')gridon;holdon;%%绘制降电压人为机械特性曲线str_y=800;fori=1:length(n1);n=n1(i).*(1-s);forcoef
7、=0.75:-0.25:0.25;U1p=U1*coef;T1=ns(i)*(m1*p*U1p^2*r2)./s./(w1.*((r1+r2./s).^2+(x1+x2)^2));plot(T1,n,'k-');ifi==1str=strcat('U1=',num2str(int16(U1p)),'V');str_y=str_y-150;text(max(T1),str_y,str,'FontSize',9,'Color','black');endholdon;endend%%绘制定子串电抗前的固有机械特性曲线figure(3);fori=1:leng
8、th(n1);n=n1(i)*(1-s);T=ns(i)*(m1*p*U1^2*r2)./s.