1、中小型异步电动机的matlab计算程序中小型异步电动机的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是有径向通风道
2、的净铁心长度myflag5 = 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 myfl
3、ag10=1是单层线圈,myflag10=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是四级及以上防护
4、式,myflag14=3是二级封闭型自扇冷式,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)
5、;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_
6、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-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
7、.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
8、/(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.007
9、7;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;%myflag2 myflag2=1是双层槽绝缘占面积myf
10、lag2=0是单层槽绝缘占面积switch 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 +
11、2*delta; case 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
12、*alpha_1*N_phi1_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;
13、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; %初取值%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*l
14、t*hf1_1; %定子轭部截面积%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
15、/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; %磁场强度%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; otherw
16、ise Lt1 = 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); %转子轭部磁路计算长度%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)
17、/(t1*(5*delta+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 + F
18、t1 + 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_; %励磁电抗标幺值%第三部分 参数计算%d1 = 0.015; %为线圈直线部分伸出铁心长度lb = lt +2*d1; %为直线部分长Kc = 1.2;beta_ = 0.852;tau_v = pi*(Dt1 + 2*(h01 + h11) + h2
19、1 + r21)/(2*p)*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 =
20、 2*d1 + Kc*tau_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_u
21、1,lambda_L1查附录% 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; %定子
22、谐波漏抗Sigma_s从参考文献电机设计图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_
23、 = 0.2*lb/(lef*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
24、1 lt2_1 =lt; 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_ =
25、 0.5*(bsk/t2)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_
26、1; %定子导线重量 C为考虑导线绝缘和引线重量的系数 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是叠片不整齐造成导
27、条电阻增加的系数 rho为电阻率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_ = sig
28、ma_1*X_sigma_*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 + 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/(p
