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

1、中小型异步电动机的MATLAB计算程序%第一部分 额定数据和主要尺寸%myflag1 = 1; %myflag1 myflag1=1是三角形接法,myflag1=0是星形接法myflag2 = 0; %myflag2 myflag2=1是双层槽绝缘占面积,myflag2=0是单层槽绝缘占面积myflag3 = 1; %myflag3 myflag3=1是无径向通风道的铁心长度,myflag3=0是定转子径向通风道不交错，其它是通风道交错myflag4 = 1; %myflag4 myflag4=1是无径向通风道的净铁心长度,myflag4=0是有径向通风道的净铁心长度myflag5 = 2;

2、%myflag5 myflag5=1是双层线圈,myflag5=2是单层线圈myflag6 = 0; %myflag6 myflag6=1是平底槽,myflag6=0是圆底槽myflag7 = 0; %myflag7 myflag7=1是平底槽,myflag7=0是圆底槽myflag8 = 1; %myflag8 myflag8=1是圆底槽,myflag8=0是半开口平底槽,其它为开口平底槽myflag9 = 1; %myflag9 myflag9=1是半开口槽和半闭开口槽,myflag9=0是开口槽myflag10 = 1; %myflag10 myflag10=1是单层线圈,myflag1

3、0=0是双层线圈myflag11 = 1; %myflag11 myflag11=1是无径向通风道,myflag11=0是有径向通风道myflag12 = 3; %myflag12 myflag12=1是双层叠绕组,myflag12=2单层同心式,myflag12=3单层同心式(分组的）、交叉式,myflag12=4 单层链式myflag13 = 1; %myflag13 myfalg13=1是无径向通风道,myflag13=0是径向通风道myflag14 = 4; %myflag14 myflag14=1是二级防护式,myflag14=2是四级及以上防护式,myflag14=3是二级封闭型自

4、扇冷式,myflag14=4是四级及以上封闭型自扇冷式myflag15 = 1; %myflag15 myflag15=1是半闭口槽,myflag15=0 是开口槽PN = 11000;Un = 380; %额定电压f = 50; %频率cos_phi = 0.84;eta_1 = 0.88;m1 = 3;p = 2; %极对数q1 = 3; %每极每相槽数%myflag1 myflag1=1是三角形接法，myflag1=0是星形接法。switch myflag1 case 1 UN_phi = Un; case 0 UN_phi = Un/sqrt(3);endIkw = PN/(m1*UN

5、_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*

6、B_delta_1*n_1);Lambda=1.4; %由参考文献电机设计表10-2Dt1_1 = (2*p*V/(Lambda*pi)(1/3);Dt1_D = 0.64; %由参考文献电机设计表10-3 Dt1_D表示Dt1/DD1_1 = Dt1_1/(Dt1_D); %定子冲片外径D1 = D1_1;Dt1 = D1*(Dt1_D);lef = V/(Dt12); %铁心有效长度lt = lef - 0.006;delta = 0.3*(0.4+7*sqrt(Dt1*lef)*(1/103);D2=Dt1-2*delta; %转子外径Dt2 = 0.06; %转子内径由转轴直径决定to

7、u=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); %每

8、相串联导体数 用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

9、 = 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;%myflag2 myflag2=1是双层槽绝缘占面积myflag2=0是单层槽绝缘占面积swit

10、ch myflag2 case 1 At = Delta_t*(2*hs_1+pi*r21+2*r21+b11); case 0 At = Delta_t*(2*hs_1+pi*r21);endAef = As - At; %槽有效面积Nt1 = 2; %并绕根数d =1.38*10(-3); %绝缘后直径Sf = Nt1*Ns1*d2/Aef; % ×100% %槽满率%myflag3 myflag3=1是无径向通风道的铁心长度，myflag3=0是定转子径向通风道不交错，其它是通风道交错。switch myflag3 case 1 lef = lt + 2*delta; case

11、 0 lef = lt-nv1*bv1_1; otherwise lef = lt-(nv1*bv1_1+nv2*bv2_1);end %lt铁心长度kfe=0.95;%myflag4 myflag4=1是无径向通风道的净铁心长度 myflag4=0是有径向通风道的净铁心长度switch myflag4 case 1 lfe = kfe*lt; case 0 lfe = 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_ph

12、i1_1/Z1; %每槽导体数%myflag5 myflag5=1是双层线圈myflag5=2是单层线圈switch myflag5 case 1 Ns1 = 2*Ns1_1; case 2 Ns1 = 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;K

13、s = 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; %初取值%myflag6 myflag6=1是平底槽 myflag6=0是圆底槽switch myflag6 case 1 hf1_1 = (D1-Dt1)/2-hs1; case 0 hf1_1 = (D1-Dt1)/2-hs1+r21/3;end %定子轭部计算高度Af1 =Kfe*lt*hf1_1; %定子轭

14、部截面积%myflag7 myflag7=1是平底槽 myflag7=0是圆底槽switch myflag7 case 1 hf2_2 = (D2-Dt2)/2-hs2; case 0 hf2_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; %转子齿磁密B

15、t2 = Fs*Phi/At2; %定子齿磁密Bf1 = Phi/2*Af1; %转子轭磁密Bf2 = Phi/2*Af2; %定子轭磁密B_delta = Fs*Phi/A_delta; %空气隙磁密Ht1 =15.74;Ht2 =9.5; %磁场强度%myflag8 myflag8=1是圆底槽 myflag8=0是半开口平底槽 其它为开口平底槽switch myflag8 case 1 Lt1 = (h11+h21)+r21/3; Lt2 = (h12+h22)+r22/3; case 0 Lt1 = h11+h21 ; Lt2 = h12+h22; otherwise Lt1 = Hs1

16、; Lt2 = h12+h22;End %Lt1定子,Lt2转子齿部磁路计算长度Lf1_1 = pi*(D1-hf1_1)/(2*p*2); %定子轭部磁路计算长度Lf2_1 = pi*(Dt1-hf2_2)/(2*p*2); %转子轭部磁路计算长度%myflag9 myflag9=1是半开口槽和半闭开口槽 myflag9=0是开口槽switch myflag9 case 1 K_delta=t1*(4.4*delta+0.75*b01)/(t1*(4.4*delta+0.75*b01)-b012); case 0 K_delta = t1*(5*delta+b01)/(t1*(5*delta

17、+b01)-b012);enddelta_ef = K_delta*delta; %有效气隙长度Ft1 = Ht1*Lt1; %定子齿部磁压降Ft2 = Ht2*Lt2; %转子齿部磁压降Cf1 =0.557; Cf2=0.445;Hf1 = 8.23*102; Hf2 = 7.29*102;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_d

18、elta; %饱和系数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_; %励磁电抗标幺值%第三部分 参数计算%d1 = 0.015; %为线圈直线部分伸出铁心长度lb = lt +2*d1; %为直线部分长Kc = 1.2;beta_ = 0.852;tau_v = pi*(Dt1 + 2*(h01 + h11) + h21 + r21)/(2*p

19、)*beta_; %单层同心式或交叉式线圈beta取平均值Cs = tau_v/2*cos(alpha);%cos(alpha) = sqrt(1 - sin(alpha)2);sin_alpha = (b11 + 2*r21)/(b11 + 2*r21 + 2*b11);%myflag10 myflag10=1是单层线圈 myflag10=0 是双层线圈switch myflag10 case 1 lc = lb + Kc * tau_v; case 0 lc = lb + 2 * Cs;endfd = Cs*sin(alpha); %双层线圈端部轴向投影长lb = 2*d1 + Kc*ta

20、u_v; %单层线圈端部平均长Cx = 4*pi*f*mu_0*(N1*Kdp1)2*lef*PN/(m1*p*(UN_phi)2); % = 0.263*(N1*Kdp12)*lef*Pn/p*(U_Nphi)2*(1/103); %漏抗系数Ku1 =1.0; KL1 =1.0;lambda_u1 = h01/b01+2*h11/(b01+b11);lambda_L1 = 0.0765;lambda_s1 = Ku1*lambda_u1 + KL1*lambda_L1; %lambda_s1为定子槽比漏磁导 Ku1,KL1查参考文献电机设计图4-9，lambda_u1,lambda_L1查附

21、录% myflag11=1是无径向通风道 myflag11=0是有径向通风道switch myflag11 case 1 lt1_1 = lt; %无径向通风道时 case 0 lt1_1 = lt - nv1*bv1_1; end %有径向通风道时 bv1_1是通风道损失宽度 ，从图参考文献电机设计附1-9查出Xs1_ = 2*m1*p*lt1_1*lambda_s1/(Z1*Kdp12*lef)*Cx; %定子槽漏抗Sigma_s = 0.0129;X_delta1_ = m1*tau*Sigma_s/(pi2*delta_ef*Kdp12*Ks)*Cx; %定子谐波漏抗Sigma_s从参

22、考文献电机设计图4-10或附录八中查出tau_y = 0.1292;%myflag12 myflag12=1是双层叠绕组,myflag12=2单层同心式,myflag12=3单层同心式(分组的）、交叉式,myflag12=4 单层链式switch myflag12 case 1 XE1_ = 1.2*(d1 +0.5*fd)/lef*Cx; case 2 XE1_ = 0.67*(lb - 0.64*tau_y)/(lef*Kdp12)*Cx; case 3 XE1_ = 0.47*(lb-0.64*tau_y)/(lef*Kdp12)*Cx; case 4 XE1_ = 0.2*lb/(le

23、f*Kdp12)*Cx;end %双层叠绕组 XE1表示定子端部漏抗X_sigma1_ = Xs1_+X_delta1_+XE1_ ; %定子漏抗lambda_u2 = h02/b02;lambda_L = 1.806; lambda_L2 = 2*h12/(b02+b12)+lambda_L;lambda_s2 = lambda_u2+lambda_L2; %转子槽比漏抗 lambda_u2 lambda_L2查参考文献电机设计附录四%myflag13 myfalt13=1是无径向通风道 myflag13=0 是有径向通风道switch myflag13 case 1 lt2_1 =lt;

24、case 0 lt2_1 = lt-nv2*bv2_1;end %bv2从参考文献电机设计附录19中查出Xs2_ = 2*m1*p*lt2_1*lambda_s2/(Z2*lef)*Cx; %转子槽漏抗Sigma_R = 0.0195;X_delta2_ = m1*tau*Sigma_R/(pi2*delta_ef*Ks)*Cx; %转子谐波漏抗 Sigma_R从参考文献电机设计图4-11或附录九查出DR = 0.1285;XB2_ = 0.757/lef*(lb - lt)/1.13 + DR/2*p)*Cx;%转子端部漏抗 见参考文献电机设计图附1-5Xsk_ = 0.5*(bsk/t2)

25、2*X_delta2_; %转子斜槽漏抗X_sigma2_ = Xs2_ + X_delta2_ + XB2_ + Xsk_; %转子漏抗X_sigma_ = X_sigma1_ + X_sigma2_; %总漏抗rho_0 = 0.0217*(1/106);Ac1_1 = 1.327*10(-6);R1 = rho_0*(2*N1*lc/(Nt1*Ac1_1*alpha_1); %rho为铝 铜的密度R1_ = R1*Ikw/UN_phi; %定子相电阻标幺值C =1.1;rho_1 = 8.9*103;Gw = C*lc*Ns1*Z1*Ac1_1*Nt1*rho_1; %定子导线重量 C为

26、考虑导线绝缘和引线重量的系数 rho_1为导线密度KFe = 0.95;rho_F_1 = 7.8*103;GFe = KFe*lt*(D1 + delta)2*rho_F_1; %GFe为硅钢片重量delta = 5*(1/103); %delta为冲剪余量rho_F_1 = 7.8*103;KB = 1.04;rho = 0.0434*(1/106);lB = 0.155;AB = 131.75*10(-6);AR = 400*10(-6);RB_1 = rho*KB*lB/AB*(4*m1*(N1*Kdp1)2)/Z2; %导条电阻折算值 KB是叠片不整齐造成导条电阻增加的系数 rho为

27、电阻率RR_1 = rho*(DR*Z2)/(2*pi*p2*AR)*(4*m1*(N1*Kdp1)2/Z2); %端环电阻折算值RB_ = RB_1*Ikw/UN_phi; %导条电阻标幺值RR_ = RR_1*Ikw/UN_phi; %端环电阻标幺值R2_ = RB_ + RR_; %转子电阻标幺值%第四部分 工作性能计算%eta = 1;while (eta-eta_1)/eta) > (0.5/100) I1p_ = 1/eta_1; %满载时定子电流有功分量标幺值 Xms_ = 2.404; sigma_1 = 1+X_sigma_/Xms_;Ix_ = sigma_1*X_s

28、igma_*I1p_2*(1+(sigma_1*X_sigma_*I1p_)2); %满载时定子电流无功分量标幺值 I1Q_ = Im_+Ix_; %满载时定子电流无功分量标幺值KB = 1-(I1p_*R1_+I1Q_*X_sigma_); %满载电势标幺值KB表示1-epsilon_LEP = 1-Im_*X_sigma_; %空载电势标幺值 EP表示1-epsilon_0 Bt10 = EP/KB*Bt1; %空载时定子齿磁密 Bt20 = EP/KB*Bt2; %空载时转子齿磁密 Bf10 = EP/KB*Bf1; %空载时定子轭磁密 Bf20 = EP/KB*Bf2; %空载时转子轭

29、磁密 B_delta0 = EP/KB*B_delta; %空载时气隙磁密 Ht10 = 22.52*102; Ht20 = 11.86*102; Hf10 = 10*102; Hf20 = 8.79*102; Ft10 = Ht10*Lt1; %空载定子齿磁压降 Ft20 = Ht20*Lt2; %空载转子齿磁压降 Ff10 = Cf1*Hf10*Lf1_1; %空载定子轭磁压降 Ff20 = Cf2*Hf20*Lf2_1; %空载转子轭磁压降F_delta_0 = K_delta*delta*B_delta0/mu_0; %空载气隙磁压降 F00 = F_delta_0 + Ft10 +

30、Ft20 + Ff10 + Ff20; %空载总磁压降 Im0 = 2*p*F00/(0.9*m1*N1*Kdp1); %空载磁化电流 I1_ = sqrt(I1p_2+I1Q_2); %定子电流标幺值 I1 = I1_*Ikw; %定子电流实际值 J1 = I1/(alpha_1*Nt1*Ac1_1); %定子电流密度 A1 = m1*N_phi1*I1/(pi*Dt1); %线负荷 I2_ = sqrt(I1p_2+Ix_2); %转子电流标幺值 I2 = I2_*Ikw*m1*N_phi1*Kdp1/Z2; %转子电流实际值 IR = I2*Z2/(2*pi*p); %端环电流实际值 J

31、B = I2/AB; %转子电流导条电密 JR = IR/AR; %转子电流端环电密 Pcu1_ = I1_2*R1_; Pcu1 = Pcu1_*PN; %定子电气损耗 PAl2_ = I2_2*R2_; PAl2 = PAl2_*PN; %转子电气损耗 Ps_ = 0.02; Ps = Ps_*PN; %附加损耗%myflag14 myflag14=1是二级防护式,myflag14=2是四级及以上防护式,%myflag14=3是二级封闭型自扇冷式,myflag14=4是四级及以上封闭型自扇冷式 switch myflag14 case 1 Pfw=5.5*(3/p)2*(D2)3*103;

32、 case 2 Pfw=6.5*(3/p)2*(D2)3*103; case 3 Pfw=13*(1-D1)*(3/p)2*(D2)3*103; case 4 Pfw=(3/p)2*(D1)4*104; endPfw_ = Pfw/PN ; rho_Fe_1 = 7.8*103; At1 = 9873*10(-6); Lt1 = 16.9*10(-3); Af1 = 3770*10(-6); Lf1 = 92*10(-3); Gt = 2*p*At1*Lt1*rho_Fe_1; %定子齿重量 Gf = 2*p*Af1*Lf1*rho_Fe_1; %定子轭重量 Phet =5.388; Phef

33、 =3.765; %根据Bt10和Bf10查参考文献电机设计附录六 PFet = Phet*Gt; %定子齿损耗 PFef = Phef*Gf; %定子轭损耗 %myflag15 myflag15=1是半闭口槽myflag15=0 是开口槽 switch myflag15 case 1 k1 = 2.5; k2 = 2; case 0 k1 = 3.0; k2 = 2.5; endPFe = k1*Phet+k2*Phef; %定子铁耗 PFe_ = PFe/PN; %铁耗标幺值 Sigma_p_ = Pcu1_+PAl2_+Ps_+Pfw_+PFe_; %总损耗标幺值 PN1_ = 1+Si

34、gma_p_; %输入功率 eta = 1-Sigma_p_/PN1_; %效率 eta_1 = eta+1/100;endcos_phi = I1p_/I1_; %功率因数PFetr = (1-1/2)*163;PFetf = (1-1/2.5)*70.2;PFer_=(PFetr+PFetf)/PN;SN = PAl2_/(1+PAl2_+PFer_+Ps_+Pfw_); %转差率n_N = 60*f/p*(1-SN); %转速Tm_ = (1-SN)/(2*(R1_+sqrt(R1_2+X_sigma_2); %最大转矩倍数%第五部分 起动性能计算%Ist = 100.13;Ist_1

35、= 3.4*Tm_*Ikw; %起动电流假定值while (Ist-Ist_1)/Ist > (1/100)Fst = Ist_1*Ns1/alpha_1*0.707*(Ku1+Kd12*Kp1*Z1/Z2)*sqrt(EP);%起动时定子转子槽磁势平均值 beta_0 = 0.64+2.5*sqrt(delta/(t1+t2);BL = mu_0*Fst/2*delta*beta_0; %空气隙时漏磁场的虚拟磁密Ks = 0.46; %1-Ks = 0.54; %漏抗饱和系数由参考文献电机设计图10-18查Cs1 = (t1-b01)*(1-Ks); %齿顶漏磁饱和引起的定子齿顶宽度的减少Cs2= (t2-b02)*(1-Ks