電機(jī)設(shè)計(jì)matlab程序

電機(jī)設(shè)計(jì)matlab程序電機(jī)設(shè)計(jì)matlab程序電機(jī)設(shè)計(jì)matlab程序xxx公司電機(jī)設(shè)計(jì)matlab程序文件編號(hào)：文件日期：修訂次數(shù)：第1.0次更改批準(zhǔn)審核制定方案設(shè)計(jì)，管理制度%電機(jī)設(shè)計(jì)程序clearallformatshortem1=3;p=2;f=50%1.額定功率PN=*10^3;%2.額定電壓（單位V，三角形接法）UN=380;UN0=380;%3.功電流(單位A)IKW=PN/(m1*UN0)%4.效率eta按照技術(shù)條件的規(guī)定eta=eta=;%5.功率因數(shù)cos(phi)=,按照技術(shù)條件的規(guī)定cos(phi)=phi=acos;cos(phi);%6.極對(duì)數(shù)p=2p=2;%7.定轉(zhuǎn)子槽數(shù)：每極每相槽數(shù)取整數(shù)。參考類(lèi)似規(guī)格電機(jī)取q1=3,則Z1=2m1pq1，再查表10-8選Z2=32，并采用轉(zhuǎn)子斜槽。q1=3;Z1=2*m1*p*q1Z2=32;%8.定轉(zhuǎn)子每極槽數(shù)Zp1=Z1/(2*p)Zp2=Z2/(2*p)%9.確定電機(jī)的主要尺寸；一般可參考類(lèi)似電機(jī)的主要尺寸來(lái)確定Di1和lef.現(xiàn)按10-2中的KE1=*log(PN/1000)*p+P1=KE1*PN/(eta*cos(phi))alphap1=;KNm1=;Kdp1=;A1=25000;Bdelta1=;n1=1450;V=(alphap1*KNm1*Kdp1))*(1/(A1*Bdelta1))*(P1/n1)D1=;%鐵心的有效長(zhǎng)度Di1=;lef=V/((Di1)^2)%氣隙的確定%參考類(lèi)似產(chǎn)品或由經(jīng)驗(yàn)公式（10-10a），得lt=;delta=lef=lt+2*deltaD2=Di1-2*delta%轉(zhuǎn)子內(nèi)徑先按轉(zhuǎn)軸直徑?jīng)Q定（以后再校驗(yàn)轉(zhuǎn)子軛部磁密）Di2=;%11.極距tautau=pi*Di1/(2*p)%12.定子齒距t1t1=(pi*Di1/Z1)%轉(zhuǎn)子齒距t2t2=(pi*D2/Z2)bsk=;%15.設(shè)計(jì)定子繞組Nphi11=eta*cos(phi)*pi*Di1*A1/(m1*IKW)%取并聯(lián)支路a1=1,由式（10-15），可得每槽導(dǎo)體數(shù)a1=1;Ns1=47%16.每相串聯(lián)導(dǎo)體數(shù)Nphi1Nphi1=Ns1*Z1/(m1*a1)%每相串聯(lián)匝數(shù)N1N1=Nphi1/2%17.繞組線(xiàn)規(guī)設(shè)計(jì)%初選定子電密J11=5.0A/mm^2,由式（10-16），計(jì)算導(dǎo)線(xiàn)并繞根數(shù)和每根導(dǎo)線(xiàn)面積的乘積。J11=;%其中定子電流初步估計(jì)值I11=IKW/(eta*cos(phi))Nt1Ac11=I11/(a1*J11) %在附錄二中選用截面積相近的銅線(xiàn)，高強(qiáng)度漆包線(xiàn)Nt1=2;d1=;d=;Ac11=;Nt1Ac11=*2;%18.設(shè)定子槽形%因定子繞組為原導(dǎo)線(xiàn)散嵌，故采用梨形槽，齒部平行。初步取Bt11=,按式（10-18）估計(jì)定子齒寬Bt11=;KFe=;bt1=t1*Bdelta1/(KFe*Bt11)%初步取Bj11=,按式（10-19），估計(jì)定子軛部計(jì)算高度Bj11=;hj11=tau*alphap1*Bdelta1/(2*KFe*Bj11)%按齒寬和定子軛部計(jì)算高度的估算值作出定子槽形如圖10-24，槽口尺寸參考類(lèi)似產(chǎn)品決定，齒寬計(jì)算如下：b01=;b11=;h01=;h11=(b11-b01)/2*tan(pi/6);h21=r21=;bi11=pi*(Di1+2*h01+2*h21+2*h11)/Z1-2*r21;bi12=pi*(Di1+2*h01+2*h11)/Z1-b11;bt1=*(bi12+bi11)%齒部基本平行，齒寬%19.槽滿(mǎn)率%曹面積hs1=h11+h21;h=;As=((2*r21+b11)/2)*(hs1-h)+pi*(r21)^2/2%按附錄三，槽絕緣采用DMD復(fù)合絕緣，,槽楔為h=2mm層壓板，槽絕緣占面積Aideltai=;Ai=deltai*(2*hs1+pi*r21)%槽的有效面積Aef=As-Ai%槽滿(mǎn)率sf=Nt1*Ns1*d^2/Aef%20.繞組系數(shù)alpha1=(p*360/Z1)Kp1=1;alpha=alpha1*pi/180,q=q1;Kd1=sin(q*alpha/2)/(q*sin(alpha/2))Kdp1=Kd1*Kp1%每相有效串聯(lián)導(dǎo)體數(shù)NefNef=Nphi1*Kdp1%21.設(shè)計(jì)轉(zhuǎn)子槽形與轉(zhuǎn)子繞組%按式（10-39），預(yù)計(jì)轉(zhuǎn)子導(dǎo)條電流：Ki=;I21=Ki*I11*(3*Nphi1*Kdp1)/Z2%初步取轉(zhuǎn)子導(dǎo)條電密JB1=3.5A/mm^2,于是導(dǎo)條截面積JB1=;AB1=I21/JB1%初步取Bt21=,估算轉(zhuǎn)子齒寬Bt21=;bt21=(t2/10*Bdelta1)/(KFe*Bt21)%初步取Bj21=,估計(jì)轉(zhuǎn)子軛部計(jì)算高度Bj21=;hj21=tau*alphap1*Bdelta1/(2*KFe*Bj21)%為獲得比較好的起動(dòng)性能，采用平行槽，作槽形圖如圖10-25所示，取槽口尺寸b02=1mm,h02=0.5mmb02=;h02=;%齒壁不平行的槽形的寬度（按3-3）計(jì)算如下：b12=;h12=;h22=導(dǎo)條截面積（轉(zhuǎn)子槽面積）AB=(b02+b12)/2*h12+(b12+b22)/2*h22%按式（11-41）估計(jì)端環(huán)電流IR1=I21*Z2/(2*p*pi)%端環(huán)所需面積JR1=*JB1;AR11=IR1/JB1AR11=(a+b)/2*(h02+h12+h22)*10^6AR1=input('輸入端環(huán)面積AR1=')%其中端環(huán)電密JR1=,JB1=2.1A/mm^2.按照工藝要求有所需面積確定端環(huán)內(nèi)徑及厚度如圖10-25b,得端環(huán)面積AE=40*10^(-6)%(二)磁路計(jì)算%22KE1=;epsilonl1=;E1=(1-epsilonl1)*UN%23Ks(1)=;KNm(1)=;precision=1;whileprecision>forn=1:3KNm(n)=input('輸入波形系數(shù)KNm(n)=')alphap(n)=input('輸入極弧系數(shù)alphap(n)=')phi(n)=E1/(4*KNm(n)*Kdp1*f*N1)%24At1=KFe*lt*bt1*Zp1At2=KFe*lt*bt213*Zp2%25hs1=h21+h01+h11+r21;hs2=h02+h12+h22;hj11=(D1-Di1)/2-hs1+r21/3hj21=(D2-Di2)/2-hs2Aj1=KFe*lt*hj11Aj2=KFe*lt*hj21%26Adelta=tau*lef%27Fs(n)=input('輸入Fs(n)=')%28Bdelta(n)=Fs(n)*phi(n)/Adelta%29Bt1=Fs(n)*phi(n)/At1%30Bt2=Fs(n)*phi(n)/At2%31Ht1(n)=input('輸入Ht1(n)=');Ht2(n)=input('輸入Ht2(n)=');%32Kdelta1=t1*100**delta*100+*b01*100)/(t1*100**delta*100+*b01*100)-(b01*100)^2)Kdelta2=t2*100**delta*100+*b02*100)/(t2*100**delta*100+*b02*100)-(b02*100)^2)Kdelta=Kdelta1*Kdelta2deltaef=Kdelta*delta%33Lt1=(h11+h21)+1/3*r21Lt2=(h12+h22)%34Lj11=pi*(D1-hj11)/(4*p)Lj21=pi*(Di2+hj21)/(4*p)%35mu0=4*pi*10^(-7);Fdelta(n)=Kdelta*Bdelta(n)*delta/mu0%36Ft1(n)=Ht1(n)*Lt1*100Ft2(n)=Ht2(n)*Lt2*100%37Ks(2*n)=(Fdelta(n)+Ft1(n)+Ft2(n))/Fdelta(n)Ks(2*n+1)=Ks(2*n)-(Ks(2*n)-Ks(2*n-1))/3precision=abs((Ks(2*n-1)-Ks(2*n))/Ks(2*n))endendKs=Ks(6);alphap1=alphap(3);phi=phi(3);Bdelta=Bdelta(3);Fdelta=Fdelta(3);Ft1=Ft1(3);Ft2=Ft2(3);KNm=KNm(3);Fs=Fs(3);Ht1=Ht1(3);Ht2=Ht2(3);%迭代的出來(lái)的結(jié)論%38Bj1=*phi/Aj1%39Bj2=*phi/Aj2%40Hj1=input('輸入Hj1=')Hj2=input('輸入Hj2=')%41hj11/tauCj1=input('根據(jù)Bj1和hj11/tau輸入Cj1=')Fj1=Cj1*Hj1*Lj11*10^2hj21/tauCj2=input('根據(jù)Bj2和hj21/tau輸入Cj2=')Fj2=Cj2*Hj2*Lj21*10^2%42F0=Fdelta+Ft1+Ft2+Fj1+Fj2%43Im=2*p*F0/*m1*N1*Kdp1)%44Ims=Im/IKW%45Xms=(4*f*mu0*m1/pi)*((N1*Kdp1)^2/(Ks*p))*lef*(tau/deltaef)Xmss=Xms*IKW/UN%（三)參數(shù)計(jì)算%46beta=(2*8+1*7)/(3*9)tauy=pi*(Di1+2*(h01+h11)+h21+r21)*beta/(2*p)d1=;lB=lt+2*d1K0=;l0=lB+K0*tauy%47lE=2*d1+K0*tauy%48Cx=4*pi*f*mu0*(N1*Kdp1)^2*lef*PN/(m1*p*UN^2)%49KU1=;KL1=;lambdaU1=h01/b01+2*h11/(b01+b11)h212r21=h21/(2*r21),b112r21=b11/(2*r21)lambdaL1=input('輸入lambdaL1=')lambdas1=KU1*lambdaU1+KL1*lambdaL1%50Xs1s=2*m1*p/(Z1*Kdp1^2)*(lt/lef)*lambdas1*Cx%51sigmas=;Xdelta1s=(m1/(pi^2)*(tau/deltaef)*(sigmas/((Kdp1^2)*Ks)))*Cx%52XE1s=**tauy)*Cx/(lef*Kdp1^2)%53Xsigma1s=Xs1s+Xdelta1s+XE1s%54%55lambdaU2=h02/b02hb2=h22/b22,b1b2=b12/b22lambdaL=input('輸入lambdaL=');lambdaL2=2*h12/(b02+b12)+lambdaLlambdas2=lambdaU2+lambdaL2%56Xs2s=2*m1*p/Z2*(lt/lef)*lambdas2%57Z22p=Z2/(2*p)sigmaR=input('根據(jù)Z22p輸入sigmaR=')Xdelta2s=m1*tau*sigmaR/(pi^2*deltaef*Ks)%58DR=;XE2s=*Z2*DR/(2*p*lef*2*p))*(2*m1*p/Z2)%59Xsks=*(bsk/t2)^2*Xdelta2s%60Xsigma2s=(Xs2s+Xdelta2s+Xsks+XE2s)*Cx%61Xsigmas=Xsigma1s+Xsigma2s%62rhow=*10^(-6);R1=rhow*(2*N1*l0)/(Nt1*Ac11*a1)%63R1s=R1*IKW/UN%64C=;rhoCu=*10^3;GCu=C*l0*Ns1*Z1*Ac11*Nt1*rhoCurhoFe=*10^3;GFe=KFe*lt*(D1+delta)^2*rhoFe%65AR1代替ARKB=;rhow1=*10^(-6);R21=rhow1*(KB*lt/AB+Z2*DR/(2*pi*p^2*AR1*10^(-6)))*(4*m1*(N1*Kdp1)^2/Z2)RB1=rhow1*(KB*lt/AB)*(4*m1*(N1*Kdp1)^2/Z2)RBs=RB1*IKW/UNRR1=rhow1*Z2*DR/(2*pi*p^2*AR1*10^(-6))*(4*m1*(N1*Kdp1)^2/Z2)RRs=RR1*IKW/UNR2s=R21*IKW/UN%66eta(1)=;I1Ps(1)=1/eta(1);precisioneta=1;whileprecisioneta>form=1:2I1Ps(m)=1/eta(2*m-1)%67sigma1=1+Xsigma1s/XmssIXs(m)=sigma1*Xsigmas*(I1Ps(m))^2*(1+(sigma1*Xsigmas*I1Ps(m))^2)%68I1Qs(m)=Ims+IXs(m)%69epsilonL(m)=I1Ps(m)*R1s+I1Qs(m)*Xsigma1sKE=1-(I1Ps(m)*R1s+I1Qs(m)*Xsigma1s)%70epsilon0=Ims*Xsigma1s1-epsilon0%71Bt10(m)=(1-epsilon0)/(1-epsilonL(m))*Bt1Ht10(m)=input('輸入Ht10(m)=')Bt20(m)=(1-epsilon0)/(1-epsilonL(m))*Bt2Ht20(m)=input('輸入Ht20(m)=')%73Bj10(m)=(1-epsilon0)/(1-epsilonL(m))*Bj1Hj10(m)=input('輸入Hj10(m)=')%74Bj20(m)=(1-epsilon0)/(1-epsilonL(m))*Bj2Hj20(m)=input('輸入Hj20(m)=')%75Bdelta0(m)=(1-epsilon0)/(1-epsilonL(m))*Bdelta%76Ft10(m)=Ht10(m)*Lt1*100%77Ft20(m)=Ht20(m)*Lt2*100%78hj11/tau,Bj10Cj1(m)=input('輸入Cj1(m)=')Fj10(m)=Cj1(m)*Hj10(m)*Lj11*100%79hj21/tau,Bj20Cj2(m)=input('輸入Cj2(m)=')Fj20(m)=Cj2(m)*Hj20(m)*Lj21*100%80Fdelta0(m)=Kdelta*delta*Bdelta0(m)/mu0%81F00(m)=Fdelta0(m)+Fj20(m)+Fj10(m)+Ft20(m)+Ft10(m)%82Im0=2*p*F00(m)/*m1*N1*Kdp1)%(四)工作性能計(jì)算%83I1s(m)=(I1Ps(m)^2+I1Qs(m)^2)^I1(m)=I1s(m)*IKW%84J1(m)=(I1(m)/(a1*Nt1*Ac11))/10^6%85A1(m)=m1*Nphi1*I1(m)/(pi*Di1)%86I2s(m)=(I1Ps(m)^2+IXs(m)^2)^I2(m)=I2s(m)*IKW*m1*Nphi1*Kdp1/Z2IR(m)=I2(m)*Z2/(2*pi*p)%87JB(m)=I2(m)/(AB*10^6)JR(m)=IR(m)/(AR1)%88pCu1s(m)=I1s(m)^2*R1spCu1(m)=pCu1s(m)*PN%89pA12s(m)=I2s(m)^2*R2spA12(m)=pA12s(m)*PN%90pss=;ps=pss*PN%91pfw=(3/p)^2*(D1)^4*10^4pfws=pfw/PN%92Bj10phej(m)=input('輸入phej(m)=')k2=2;Gj=4*p*Aj1*Lj11*rhoFepFej(m)=k2*phej(m)*GjBt10phet(m)=input('輸入phet(m)=')Gt=2*p*At1*Lt1*rhoFek1=2;pFet(m)=k1*phet(m)*GtpFe(m)=pFej(m)+pFet(m)pFes(m)=pFe(m)/PN%93sigmaps(m)=pCu1s(m)+pA12s(m)+pss+pfws+pFes(m)%94PN1s(m)=1+sigmaps(m)%95eta(2*m)=1-sigmaps(m)/PN1s(m)eta(2*m+1)=eta(2*m)+(eta(2*m)-eta(2*m-1))/5precisioneta=abs((eta(2*m)-eta(2*m-1))/eta(2*m))endend%96I1s=I1s(2);I1Ps=I1Ps(2);phi=acos(I1Ps/I1s);cos(phi)%97pFers=((1-1/2)*pFej+(1-1/*pFet)*10^(-3)/PN*1000sN=pA12s/(1+pA12s+pFers+pss+pfws)%98nN=60*f/p*(1-sN)%99Tms=(1-sN)/(2*(R1s+(R1s^2+Xsigmas^2)^))%(五)啟動(dòng)性能%100K(1)=3;%101precisionist=1;K(1)=3;Ist(1)=K(1)*Tms*IKWwhileprecisionist>fork=1:3Fst(k)=Ist(2*k-1)*Ns1/a1**(KU1+Kd1^2*Kp1*Z1/Z2)*(1-epsilon0)^beta0=+*(delta/(t1+t2))^BL=mu0*Fst(k)/(2*delta*beta0)%102KZ(k)=input('輸入啟動(dòng)漏抗飽和KZ（k）=')%103cs1(k)=(t1-b01)*(1-KZ(k))%104cs2(k)=(t2-b02)*(1-KZ(k))%105deltalambdaU1(k)=(h01+*h11)/b01*(cs1(k)/(cs1(k)+*b01))lambdas1st(k)=KU1*(lambdaU1-deltalambdaU1(k))+KL1*lambdaL1%106Xs1sts(k)=lambdas1st(k)/lambdas1*Xs1s%107Xdelta1sts(k)=KZ(k)*Xd