中小型异步电动机的MATLAB计算程序

仅供个人参考不得用于商业用途中小型异步电动机的MATLAB计算程序%%%%%%%%%%%%%%%%%%%%%%%%%第一部分额定数据和主要尺寸%%%%%%%%%%%%%%%%%%%%%%%%%myflag1=1;%myflag1myflag1=1是三角形接法,myflag1=0是星形接法myflag2=0;%myflag2myflag2=1是双层槽绝缘占面积,myflag2=0是单层槽绝缘占面积myflag3=1;%myflag3myflag3=1是无径向通风道的铁心长度,myflag3=0是定转子径向通风道不交错，其它是通风道交错myflag4=1;%myflag4myflag4=1是无径向通风道的净铁心长度,myflag4=0是有径向通风道的净铁心长度myflag5=2;%myflag5myflag5=1是双层线圈,myflag5=2是单层线圈myflag6=0;%myflag6myflag6=1是平底槽,myflag6=0是圆底槽myflag7=0;%myflag7myflag7=1是平底槽,myflag7=0是圆底槽myflag8=1;%myflag8myflag8=1是圆底槽,myflag8=0是半开口平底槽,其它为开口平底槽myflag9=1;%myflag9myflag9=1是半开口槽和半闭开口槽,myflag9=0是开口槽myflag10=1;%myflag10myflag10=1是单层线圈,myflag10=0是双层线圈myflag11=1;%myflag11myflag11=1是无径向通风道,myflag11=0是有径向通风道myflag12=3;%myflag12myflag12=1是双层叠绕组,myflag12=2单层同心式,myflag12=3单层同心式(分组的）、交叉式,myflag12=4单层链式myflag13=1;%myflag13myfalg13=1是无径向通风道,myflag13=0是径向通风道仅供个人参考不得用于商业用途myflag14=4;%myflag14myflag14=1是二级防护式,myflag14=2是四级及以上防护式,myflag14=3是二级封闭型自扇冷式,myflag14=4是四级及以上封闭型自扇冷式myflag15=1;%myflag15myflag15=1是半闭口槽,myflag15=0是开口槽PN=11000;Un=380;%额定电压f=50;%频率cos_phi=0.84;eta_1=0.88;m1=3;p=2;%极对数q1=3;%每极每相槽数%myflag1myflag1=1是三角形接法，myflag1=0是星形接法。switchmyflag1case1UN_phi=Un;case0UN_phi=Un/sqrt(3);endIkw=PN/(m1*UN_phi);%功电流Z1=2*m1*p*q1;%定子槽数Z2=26;%转子槽数Zp1=Z1/(2*p);Zp2=Z2/(2*p);%定转子每极槽数KB_2=0.0108*log10(PN)-0.013*p+0.931;%满载电势标幺值alpha_p_1=0.68;Knm_1=1.10;Kdp1_1=0.96;A_1=26000;%由参考文献［电机设计］图10-2B_delta_1=0.67;n_1=1450;p_1=KB_2*PN/(eta_1*cos_phi);V=6.1*1*p_1/(alpha_p_1*Knm_1*Kdp1_1*A_1*B_delta_1*n_1);Lambda=1.4;%由参考文献［电机设计］表10-2Dt1_1=(2*p*V/(Lambda*pi))^(1/3);Dt1__D=0.64;%由参考文献［电机设计］表10-3Dt1__D表示Dt1/DD1_1=Dt1_1/(Dt1__D);%定子冲片外径D1=D1_1;Dt1=D1*(Dt1__D);lef=V/(Dt1^2);%铁心有效长度仅供个人参考不得用于商业用途lt=lef-0.006;delta=0.3*(0.4+7*sqrt(Dt1*lef))*(1/10^3);D2=Dt1-2*delta;%转子外径Dt2=0.06;%转子内径由转轴直径决定tou=pi*Dt1/(2*p);%极距t1=pi*Dt1/Z1;%定子齿距t2=pi*D2/Z2;%转子齿距bsk=0.0145;%转子斜槽宽一般取一个定子齿距t1alpha_1=1;%并联支路alpha_1为1N_phi1_1=(eta_1*cos_phi*pi*Dt1*A_1)/(m1*Ikw);Ns1=(m1*alpha_1*N_phi1_1)/Z1;N_phi1=Ns1*Z1/(m1*alpha_1);J1_1=5.0;%定子电密I1_1=Ikw/eta_1*cos_phi;%定子电流初步值NA=I1_1/(alpha_1*J1_1);%每相串联导体数用NA表示Nt1_1*Ac1_1KI=0.089;%查参考文献［电机设计］表10-10alpha=p*2*pi/Z1;Kd1=sin(q1*alpha/2)/(q1*sin(alpha/2));%分布系数Kp1=1;%短距系数Kdp1=Kd1*Kp1;%绕组系数I2_1=KI*I1_1*3*N_phi1_1*Kdp1/Z2;%转子导条电流IR_1=I2_1*Z2/(2*pi*p);%端环电流JR_1=2.1;AR_1=IR_1/JR_1;%端环面积h=0.002;%槽锲b01=0.0038;b11=0.0077;bt1=0.00745;bt2=0.01135;h01=0.0008;h11=0.0011;h21=0.0152;hs1=0.0211;hs2=0.0289;r21=0.0051;r22=0.0024;hs_1=0.0152;%定子槽型b02=0.001;b12=0.0048;h02=0.0005;h12=0.0011;h22=0.0249;%转子槽型As=(2*r21+b11)*(hs_1-h)/2+pi*r21*r21;%槽面积Delta_t=0.0003;仅供个人参考不得用于商业用途%myflag2myflag2=1是双层槽绝缘占面积myflag2=0是单层槽绝缘占面积switchmyflag2case1At=Delta_t*(2*hs_1+pi*r21+2*r21+b11);case0At=Delta_t*(2*hs_1+pi*r21);endAef=As-At;%槽有效面积Nt1=2;%并绕根数d=1.38*10^(-3);%绝缘后直径Sf=Nt1*Ns1*d^2/Aef;%×100%%槽满率%myflag3myflag3=1是无径向通风道的铁心长度，myflag3=0是定转子径向通风道不交错，其它是通风道交错。switchmyflag3case1lef=lt+2*delta;case0lef=lt-nv1*bv1_1;otherwiselef=lt-(nv1*bv1_1+nv2*bv2_1);end%lt铁心长度kfe=0.95;%myflag4myflag4=1是无径向通风道的净铁心长度myflag4=0是有径向通风道的净铁心长度switchmyflag4case1lfe=kfe*lt;case0lfe=kfe*(lt-nv*bv);endN_phi1_1=eta_1*cos_phi*pi*Dt1*A_1/(m1*Ikw);%每相串联导体数alpha_1=1;%并联支路Ns1_1=m1*alpha_1*N_phi1_1/Z1;%每槽导体数%myflag5myflag5=1是双层线圈myflag5=2是单层线圈switchmyflag5case1Ns1=2*Ns1_1;case2Ns1=Ns1_1;end%每线圈匝数N2=N_phi1*Kdp1;%每相有效串联导体数N1=N_phi1/2;%每相串联匝数仅供个人参考不得用于商业用途%%%%%%%%%%%%%%%%%%%%%%%%第二部分磁路计算%%%%%%%%%%%%%%%%%%%%%%%%%%KB_1=0.927;%KB_1=1-epsilon_l==0.927;%E1=(1-epsilon_l)*UN_phi==0.927*UN_phi;%E1为满载相电势E1=0.927*UN_phi;Ks_1=1.15;Ks=1.148;Knm=1.0975;Phi=E1/(4*Knm*Kdp1*f*N1);%每极磁通Kfe=0.95;At1=Kfe*lt*bt1*Zp1;%定子每极下齿部截面积At2=Kfe*lt*bt2*Zp2;%转子每极下齿部截面积Bt1_1=1.4;%初取值%myflag6myflag6=1是平底槽myflag6=0是圆底槽switchmyflag6case1hf1_1=(D1-Dt1)/2-hs1;case0hf1_1=(D1-Dt1)/2-hs1+r21/3;end%定子轭部计算高度Af1=Kfe*lt*hf1_1;%定子轭部截面积%myflag7myflag7=1是平底槽myflag7=0是圆底槽switchmyflag7case1hf2_2=(D2-Dt2)/2-hs2;case0hf2_2=(D2-Dt2)/2-hs2+r22/3;end%转子轭部计算高度Af2=Kfe*lt*hf2_2;%转子轭部截面积tau=0.1335;A_delta=tau*lef;%空气隙截面积alpha_p1=0.67;%计算极弧系数Fs=1/alpha_p1;%波幅系数B_s=Fs*Phi/A_delta;%气隙磁密Bt1=Fs*Phi/At1;%转子齿磁密Bt2=Fs*Phi/At2;%定子齿磁密Bf1=Phi/2*Af1;%转子轭磁密Bf2=Phi/2*Af2;%定子轭磁密B_delta=Fs*Phi/A_delta;%空气隙磁密Ht1=15.74;Ht2=9.5;%磁场强度%myflag8myflag8=1是圆底槽myflag8=0是半开口平底槽其它为开口平仅供个人参考不得用于商业用途底槽switchmyflag8case1Lt1=(h11+h21)+r21/3;Lt2=(h12+h22)+r22/3;case0Lt1=h11+h21;Lt2=h12+h22;otherwiseLt1=Hs1;Lt2=h12+h22;End%Lt1定子,Lt2转子齿部磁路计算长度Lf1_1=pi*(D1-hf1_1)/(2*p*2);%定子轭部磁路计算长度Lf2_1=pi*(Dt1-hf2_2)/(2*p*2);%转子轭部磁路计算长度%myflag9myflag9=1是半开口槽和半闭开口槽myflag9=0是开口槽switchmyflag9case1K_delta=t1*(4.4*delta+0.75*b01)/(t1*(4.4*delta+0.75*b01)-b01^2);case0K_delta=t1*(5*delta+b01)/(t1*(5*delta+b01)-b01^2);enddelta_ef=K_delta*delta;%有效气隙长度Ft1=Ht1*Lt1;%定子齿部磁压降Ft2=Ht2*Lt2;%转子齿部磁压降Cf1=0.557;Cf2=0.445;Hf1=8.23*10^2;Hf2=7.29*10^2;Ff1=Cf1*Hf1*Lf1_1;%定子轭部磁压降Ff2=Cf2*Hf2*Lf2_1;%转子轭部磁压降mu_0=1.25*10^(-6);F_delta=K_delta*delta*B_delta/mu_0;%空气隙磁压降K_s=(F_delta+Ft1+Ft2)/F_delta;%饱和系数Ks_2=Ks-(Ks-Ks_1)/3;Ks_1=Ks_2;F0=F_delta+Ft1+Ft2+Ff1+Ff2;%总磁压降Im=2*p*F0/(0.9*m1*N1*Kdp1);%满载磁化电流Im_=Im/Ikw;%满载磁化电流标幺值Xm_=1/Im_;%励磁电抗标幺值%%%%%%%%%%%%%%%%%%%%%%%%%第三部分参数计算%%%%%%%%%%