三相电机设计程序清单.docx
- 文档编号:3287576
- 上传时间:2022-11-21
- 格式:DOCX
- 页数:19
- 大小:40.38KB
三相电机设计程序清单.docx
《三相电机设计程序清单.docx》由会员分享,可在线阅读,更多相关《三相电机设计程序清单.docx(19页珍藏版)》请在冰豆网上搜索。
三相电机设计程序清单
附录A程序清单
('主要参数')
formatshortg
Pn=input('Pn=');Un=input('Un=');m1=input('m1=');f=input('f=');p=input('p=');a=input('Y接法时a为1,Δ接法时a为2,a=')
switcha
case1
Unfai=Un/sqrt(3)%Y接法
case2
Unfai=Un%三角接法
end
Ikw=Pn/(m1*Unfai)
nnpie=0.88%效率
q1=input('定子每极每相槽数q1=')
Z1=2*m1*p*q1
cosfaipie=0.84%暂时规定(功率因数)
Zp1=Z1/(2*p)
ifp==1
Z2=16
elseifp==2
Z2=26
elseifp==3
Z2=44
elseZ2=56
end
Zp2=Z2/(2*p)
ifp==1
K=0.92+0.00866*log(Pn/1000)
elseK=0.931+0.0108*log(Pn/1000)-0.013*p
end
Ppie=K*Pn/(nnpie*cosfaipie);%计算功率
ifp==1
A=[16.5+8*log10(Pn/1000)]*1000;Bair=0.58;lemd=0.7;K2=0.57;Bair=0.57;%K2为Di1/D
elseifp==2
A=[19.8+5.95*log10(Pn/1000)]*1000,Bair=0.67;lemd=1.4;K2=0.64;Bair=0.67;
elseA=[18+6.5*log10(Pn/1000)]*1000;Bair=0.73;lemd=2.3;K2=0.675;Bair=0.65;
end
ap=0.68;Knm=1.1;Kdp1=0.96;npie=60*f/p-50;%暂时规定
V=6.1*Ppie/(ap*Knm*Kdp1*A*Bair*npie);
Di1=[2*p*V/(lemd*pi)]^(1/3);
D1=Di1/K2;
D=[0.120.1450.1670.210.2450.280.3270.3680.423];
HH=[0.090.010.1120.1320.160.180.20.2250.250.28];
i=1;
whileD1>D(i)
i=i+1;
end
D1=D(i-1)
Di1=D1*K2;
Di1=fix((Di1+0.005)*100)/100
HH1=HH(i-1);
lef=V/[Di1^2];
li=fix(lef*100)/100
air=0.3*[0.4+7*(Di1*lef)^(1/2)]*0.001%气隙
lef=lef+2*air
D2=Di1-2*air
Di2=D2/2.8
tao=pi*Di1/(2*p)%tao为极距
t1=pi*Di1/Z1%定子齿距
t2=pi*D2/Z2%转子齿距
bsk=fix(t2*2000)/2000;
Nfai1pie=fix(nnpie*cosfaipie*pi*Di1*A/(m1*Ikw));
a1=1;%规定
Ns1pie=m1*a1*Nfai1pie/Z1;
Ns1=fix(Ns1pie)
Nfai1=Ns1*Z1/(m1*a1)%每项串联导体数
N1=Nfai1/2%每项串联匝数
Jpie1=5.0;%初选?
?
?
?
?
?
?
?
?
?
?
?
?
?
?
?
?
?
?
?
?
?
?
I1pie=Ikw/(nnpie*cosfaipie)
J1=5.0;
S=I1pie/(a1*J1);
%导线并绕根数和每根截面积的乘积
d1pie=[0.500.530.560.600.630.670.690.710.750.770.800.830.850.900.930.951.001.061.121.181.251.301.351.401.451.501.561.601.701.801.902.002.122.242.362.50];
Spie=[0.19640.22060.24630.28270.31170.35260.37390.35590.44180.46570.52270.5410.56750.63620.67930.70880.78540.88250.98521.09361.22721.32731.43141.53941.65131.76721.91142.01062.26982.54472.83533.14163.52993.94084.37444.9087];
ifS>0.6;
i=1;
while(S/2)>Spie(i)
i=i+1;
end
Ao1pie=Spie(i)*10^(-6);
S=2*Spie(i)
Ni1=2;
elsei=1;
whileS>Spie(i)
i=i+1;
end
Ao1pie=Spie(i)*10^(-6);
S=Spie(i)%平方mm
Ni1=1
end
d1=d1pie(i)*0.001%mm
d=d1+0.00008
h01=0.0008;Kfe=0.95;Fs=1.46;
Bj1pie=1.25;rfappie=0.68;Bairpie=0.67;
hj1pie=tao*rfappie*Bairpie/(2*Kfe*Bj1pie);
K2=2;Sf=0.81;
whileSf>0.78
b01=Ni1*d+0.0014%m
rfapie=pi/6;K1=1.1;
%求方程求r21
b11=pi*[Di1*(1-1/(Kfe*K2))+2*h01-b01/(3^(1/2))]/(Z1-pi/(3^(1/2)))
bi1=pi*Di1/(Kfe*Z1*K2)
A=[2*Z1,-2*pi;2,3];
B=[pi*(Di1+2*h01)-Z1*bi1,1.5*(D1-Di1-2*0.0008)-0.75*K1*Z1*bi1/(p*Fs)]';
X=A\B;
r21=X
(1);H=X
(2);
h11=(b11-b01)/(2*3^(1/2))
h21=H-h11
As=(2*r21+b11)*(H-0.002)/2+pi*r21^
(2)/2;
%dta查表
ifHH1<0.113
dta=0.00025;
elseifHH1<0.17
dta=0.0003;
elsedta=0.00035
end
h=0.002;
Ai=dta*(2*h21+pi*r21);
Aef=As-Ai;
Sf=Ni1*Ns1*d1^2/Aef%槽满率
K2=K2+0.025;
end
K2=K2-0.025
Kp1=1;
rfa=p*2*pi/Z1;
Kd1=sin(q1*rfa/2)/(q1*sin(rfa/2));
Kdp1=Kd1*Kp1
Nfai1pei=fix(Nfai1*Kdp1+0.5)
%查表求Ki
cosfai=[0.650.700.750.800.850.900.95];
Ki=[0.740.770.820.860.900.950.985];
n1=1;
p1=polyfit(cosfai,Ki,n1);
Ki=p1
(1)*cosfaipie+p1
(2);
%查表求结束
I2pie=Ki*I1pie*3*Nfai1*Kdp1/Z2%%%预计转子导条电流
Jbpie=3.5;%初取
Abpie=I2pie/Jbpie;
Bi2pie=1.3;afappie=0.68;%初选
bi2pie=t2*Bairpie/(Kfe*Bi2pie)
Bj2pie=1.25%初选
hj2pie=tao*afappie*Bairpie/(2*Kfe*Bj2pie)
b02=0.001,h02=0.0005,
b12=D2/35,H2=D2/6.5,r22=b12/2,h12=(b12-0.001)/(2*3^(0.5)),h22=H2-h12
bi2=pi*[D2-4*(h02+h12+h22)/3]/Z2-b12
Ab=(b02+b12)*h12/2+b12*h22+pi*r22^
(2)/2%导条截面积
Irpie=I2pie*Z2/(2*pi*p)
Jrpie=0.6*Jbpie;
Arpie=Irpie/Jrpie*10^(-6)
Dr=D2-2*(h12+h22+r22);
('磁路计算')
Kepie=0.927;ppp=1;
Ke=0.927;%%不符要循环
whileppp>0.005
Kepie=Ke-(Ke-Kepie)/8
E1=Kepie*Unfai
%感应电机的Knm=f(Ks)曲线
%end
v=0.02;
Kspie=1.20
whilev>0.01
x=[1.0:
0.5:
3.0];
y1=[1.1091.081.061.0441.035];
n3=3;
p3=polyfit(x,y1,n3)
Knm=p3
(1)*Kspie^(3)+p3
(2)*Kspie^
(2)+p3(3)*Kspie+p3(4)
fai=E1/(4*Knm*Kdp1*f*N1)
Ai1=Kfe*li*bi1*Zp1
Ai2=Kfe*li*bi2*Zp2
hs1=h21+h11+h01+r21
hj1pie=(D1-Di1)/2-hs1+r21/3%定子轭部计算高度
hs2=r22+H2+h02
hj2pie=(D2-Di2)/2-hs2+r22/3%转子轭部计算高度
Aj1=Kfe*li*hj1pie
Aj2=Kfe*li*hj2pie
Aair=tao*lef%一极下空气隙截面积
%感应电机的α’p=f(Ks)曲线程序
x=[1.0:
0.5:
3.0];
y1=[0.640.7430.7850.8250.84];
n3=3;
p3=polyfit(x,y1,n3)
%end
rfappie=p3
(1)*Kspie^(3)+p3
(2)*Kspie^
(2)+p3(3)*Kspie+p3(4);
Fs=1/rfappie
Bair=Fs*fai/Aair%气隙磁密
Bi1=Fs*fai/Ai1%定子齿部磁密
Bi2=Fs*fai/Ai2%气隙磁密
%%%查表Hi1HI2
x=[0.40:
0.01:
1.89];
y=[1.381.401.421.441.461.481.501.521.541.561.581.601.621.641.661.691.711.741.761.781.811.841.861.891.911.941.972.002.032.062.102.132.162.202.242.282.322.362.402.452.502.552.602.652.702.762.812.872.932.993.063.133.193.263.333.413.493.573.653.743.833.924.014.114.224.334.445.564.674.804.935.075.215.365.525.685.846.006.166.336.526.726.947.167.387.627.868.108.368.628.909.209.509.8010.110.510.911.311.712.112.613.113.614.214.815.516.317.118.119.120.121.222.423.725.026.728.530.432.635.137.840.743.746.850.053.456.860.464.067.872.076.480.885.490.295.0100105110116122128134140146152158165172180];
n5=5;
p5=polyfit(x,y,n5);
('定子齿部磁场强度')
Hi1=p5
(1)*Bi1^5+p5
(2)*Bi1^4+p5(3)*Bi1^3+p5(4)*Bi1^2+p5(5)*Bi1+p5(6)
Hi2=p5
(1)*Bi2^5+p5
(2)*Bi2^4+p5(3)*Bi2^3+p5(4)*Bi2^2+p5(5)*Bi2+p5(6)
Kair1=t1*(4.4*air+0.75*b01)/[t1*(4.4*air+0.75*b01)-b01^
(2)]%气隙系数
Kair2=t2*(4.4*air+0.75*b02)/[t2*(4.4*air+0.75*b02)-b02^
(2)]%气隙系数
Kair=Kair1*Kair2
airef=Kair*air%有效气隙长度
Li1=h11+h21+r21/3%齿部磁路计算长度
Li2=h12+h22+r22/3%齿部磁路计算长度
Lj1pie=pi*(D1-hj1pie)/(4*p)
Lj2pie=pi*(Di2+hj2pie)/(4*p)
miu0=0.4*pi*10^(-6);
Fair=Kair*Bair*air/miu0%气隙磁压降
Fi1=100*Hi1*Li1%齿部磁压降
Fi2=100*Hi2*Li2%齿部磁压降
Ks=(Fair+Fi1+Fi2)/Fair
v=abs((Kspie-Ks)/Ks)
Kspie=Ks-(Ks-Kspie)/3;
end
Bj1=fai/(2*Aj1)%定子轭部磁密
Bj2=fai/(2*Aj2)%转子轭部磁密
('定子轭部磁场强度')%%查表Hj1Hj2
Hj1=p5
(1)*Bj1^5+p5
(2)*Bj1^4+p5(3)*Bj1^3+p5(4)*Bj1^2+p5(5)*Bj1+p5(6)
Hj2=p5
(1)*Bj2^5+p5
(2)*Bj2^4+p5(3)*Bj2^3+p5(4)*Bj2^2+p5(5)*Bj2+p5(6)
kk=hj1pie/tao
x=[1.2401.3951.550];%%图附1-3b
y=[7.4592-2.6729-0.2486];
n2=2;
p2=polyfit(x,y,n2);
ca=p2
(1)*Bj1^
(2)+p2
(2)*Bj1+p2(3)
x=[1.2401.3951.550];
y=[-3.56180.14450.9417];
n2=2;
p2=polyfit(x,y,n2);
cb=p2
(1)*Bj1^
(2)+p2
(2)*Bj1+p2(3)
x=[1.2401.3951.550];
y=[1.20330.84980.3081];
n2=2;
p2=polyfit(x,y,n2);
cc=p2
(1)*Bj1^
(2)+p2
(2)*Bj1+p2(3)
x=[1.2401.3951.550];
y=[0.44590.36160.3038];
n2=2;
p2=polyfit(x,y,n2);
cd=p2
(1)*Bj1^
(2)+p2
(2)*Bj1+p2(3)
Cj1=ca*kk^(3)+cb*kk^
(2)+cc*kk+cd;
Fj1=100*Cj1*Hj1*Lj1pie%轭部磁压降
kk2=hj2pie/tao
x=[1.2401.3951.550];%%图附1-3b
y=[-5.6566-3.8539-6.1538];
n2=2;
p2=polyfit(x,y,n2);
ca=p2
(1)*Bj1^
(2)+p2
(2)*Bj1+p2(3)
x=[1.2401.3951.550];
y=[4.50823.50583.4545];
n2=2;
p2=polyfit(x,y,n2);
cb=p2
(1)*Bj1^
(2)+p2
(2)*Bj1+p2(3)
x=[1.2401.3951.550];
y=[-0.7190-0.5929-0.6794];
n2=2;
p2=polyfit(x,y,n2);
cc=p2
(1)*Bj1^
(2)+p2
(2)*Bj1+p2(3)
x=[1.2401.3951.550];
y=[0.45010.35960.3001];
n2=2;
p2=polyfit(x,y,n2);
cd=p2
(1)*Bj1^
(2)+p2
(2)*Bj1+p2(3)
Cj2=ca*kk2^(3)+cb*kk2^
(2)+cc*kk2+cd;
Fj2=100*Cj2*Hj1*Lj2pie%轭部磁压降
F0=Fair+Fi1+Fi2+Fj1+Fj2%每极磁势
Im=2*p*F0/(0.9*m1*N1*Kdp1)%磁化电流
Imbiao=Im/Ikw%磁化电流标么值
Xms=4*f*miu0*m1*(N1*Kdp1)^
(2)*lef*tao/(pi*Ks*p*airef)%励磁电抗
Xmsbiao=Xms*Ikw/Unfai
('参数计算')
ifq1==2
('单层链式饶组'),bta=(Z1/(2*p)-1)/(Z1/(2*p))
vv=1;
elseifq1==3
('单层交叉式饶组'),bta=[(q1-1)*(Z1/(2*p))+Z1/(2*p)-2]/(3*Z1/(2*p))
vv=2;
elseifq1==4
('同心式饶组'),bta=1
vv=2;
elseifq1==5
('单层交叉式饶组'),bta=[(q1-1)*(Z1/(2*p))+Z1/(2*p)-2]/(3*Z1/(2*p))
elseifq1==6
('同心式饶组'),bta=1
vv=2;
elseifq1==7
('单层交叉式饶组'),bta=[(q1-1)*(Z1/(2*p))+Z1/(2*p)-2]/(3*Z1/(2*p))
vv=2;
elseifq1==8
('同心式饶组'),bta=1
vv=2;
end
taov=pi*[Di1+2*(h01+h11)+h21+r21]*bta/(2*p)
dd1=0.015
lb=li+2*dd1%直线部分
ifp==1
Ko=1.16
elseifp==4
Ko=1.25
elseKo=1.2
end
lo=lb+Ko*taov%半匝长
le=2*d1+Ko*taov
Cx=4*pi*f*miu0*(N1*Kdp1)^
(2)*lef*Pn/(m1*p*Unfai^
(2))
Kv1=1;Kl1=1;
lemdv1=h01/b01+2*h11/(b01+b11);
g=h21/(2*r21);gg=b11/(2*r21);
x=[5432.52.221.91.81.71.61.51.41.31.21.110.90.80.70.60.50.4];
y1=[0.05420.07860.09440.1270.15000.15560.168750.18570.19440.2060.206250.22080.25570.26870.2850.3360.35830.41250.4500.5400.63330.75];
n3=3;
p3=polyfit(x,y1,n3)
aa=p3
(1)*gg^(3)+p3
(2)*gg^
(2)+p3(3)*gg+p3(4)%%系数阿a
x=[5432.52.221.91.81.71.61.51.41.31.21.110.90.80.70.60.50.4];
y1=[0.01580.03560.0430.040.06440.061250.05430.06560.0700.073750.07920.09670.11830.11500.1300.1240.16170.15750.2000.1800.19670.230];
n3=3;
pp3=polyfit(x,y1,n3)
bb=pp3
(1)*gg^(3)+pp3
(2)*gg^
(2)+pp3(3)*gg+pp3(4)%%系数阿b
lemdl1=aa*g+bb
lemds1=Kv1*lemdv1+Kl1*lemdl1
Xs1biao=2*m1*p*li*lemds1*Cx/(Z1*Kdp1^
(2)*lef)
%%%查表4-10ss
ifq1==2
ss=0.2157*bta^3-0.4790*bta^2+0.3541*bta-0.0778;
elseifq1==3
ss=0.2157*bta^3-0.4790*bta^2+0.3541*bta-0.0778;
elseifq1==4
ss=0.1745*bta^3-0.3849*bta^2+0.2809*bta-0.0622;
elseifq1==5
ss=0.1703*bta^3-0.3726*bta^2+0.2676*bta-0.0591;
end
Xair1biao=m1*tao*ss*Cx/(Kdp1^
(2)*pi^
(2)*airef*Ks)
ifvv==1
Xe1biao=0.2*(le*Cx/(lef*Kdp1^
(2)))
elseXe1biao=0.47*(le-0.64*tao)*Cx/(lef*Kdp1^
(2))
end
X61biao=Xs1biao+Xair1biao+Xe1biao%定子漏抗标么值
lemdv2=h02/b02
g=h22/(2*r22);gg=b12/(2*r22);
aa=p3
(1)*gg^(3)+p3
(2)*gg^
(2)+p3(3)*gg+p3(4)%%系数阿a
bb=pp3
(1)*gg^(3)+pp3
(2)*gg^
(2)+pp3(3)*gg+pp3(4)%%系数阿b
lemdl=aa*g+bb
lemdl2=2*h12/(b02+b12)+lemdl
lemds2=lemdv2+lemdl2
Xs2biao=2*m1*p*li*lemds2*Cx/(Z2*lef)
ff=Z2/(2*p);
%%%查表RR4-11
RR=-6.6662*1/(ff^4)+8.3358/(ff^3)-0.8571/(ff^2)+0.0540/ff+0.0022;
Xair2biao=m1*tao*RR*Cx/(pi^
(2)*airef*Ks)
Xe2biao=0.757*Dr*Cx/(lef*2*p)
Xskbiao=0.5*(bsk/t2)^
(2)*Xair2biao
X62bia
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 三相 电机 设计 程序 清单