对转双转子永磁同步电机电磁设计

本文格式为Word版，下载可任意编辑——对转双转子永磁同步电机电磁设计%%|已知数据的输入|%clear;

%|内外转子的公共数据|

pi=3.1415926;%圆周率

PN=5.5;%内外转子等效电机的额定功率(KW)m=3;%相数p=3;%极对数

UNl=190;%额定线电压(V)f=50;%频率(HZ)etaN=93.5;%效率(%)cosfaiN=0.90;%功率因数UN=190;%额定相电压V

%注：如无特别说明，本程序默认长度单位为厘米（cm）%内外转子、定子、气隙的相关尺寸gi=0.12;go=0.585;%内外气隙长度D1=23.05;Di1=12.55;%定子内外径Dm=18.3;%定子中心线

D2i=11.52;D2o=28.30;%内外转子的外径Di2i=8.62;Di2o=24.9;%内外转子的内径

L1=19.8;L2=19.8;La=19.8;Lb=19.8;%定转子铁心长Qi=36;Qo=36;Q=36;%定子内外槽数

hi01=0.105;bi01=0.3;bi1=0.67;ri1=0.395;hi12=1.212;his1=0.1782;bi2=0.79;hi2=1.0338;%定子内槽的槽型尺寸ho01=0.175;bo01=0.42;bo1=1.44;ro1=1.154;

ho12=0.782;hos1=0.316;bo2=1.36;ho2=0.466;%定子外槽槽型尺寸h=0.2;%槽契厚度Ci=0.035;%槽绝缘层厚度

%绕组、导线、端部的相关尺寸

Nsi=20;Nso=20;Ns=20;%每槽导体数q1=2;q2=2;q=2;%内外转子的级相组a=1;%并联支路数

Nt1=2;d11=1.25;%并绕根数-线径hd1=0.012;%导线双边绝缘厚度y=6;%节距

do=2;di=2;d=2;%绕组直线部分伸出长LE1=5.25;%端部长度

%永磁体的相关尺寸%永磁体牌号N33SH

hMi=0.395;hMo=0.34;%内外的磁化方向长度Br=1.0741;%剩磁密度Hc=876;%矫顽力

Ami=114.85;Amo=218.87;%内外转子提供每极磁通的截面积bMi=5.80;bMo=11.05;%内外永磁体宽度LMi=19.8;LMo=19.8;%内外永磁体的轴向长度zata0i=1.05;zata0o=1.15;%空载漏磁系数假定值bm01i=0.560;bm01o=0.540;%永磁体空载工作点

%相关的一些固定系数mur=1.046;mu0=4*pi*1e-7;%相对磁导率RouCu=0.217;%铜电阻率1e-3*Ω*mm^2/cmrouCu=8.9;%铜线密度g/cm^3rouFe=7.8;%硅钢片密度g/cm^3deita=0.5;%剪切余量

roum=7.4;%永磁体密度g/cm^3KFe=0.93;%叠加系数cata=pi/18;%转矩角

%磁化曲线插值表格1查磁化曲线

Bcz=0.40:0.01:1.89;

H=[1.381.401.421.441.461.481.501.521.541.56...

1.581.601.621.641.661.691.711.741.761.78...1.811.841.861.891.911.941.972.002.032.06...2.102.132.162.202.242.282.322.362.402.45...2.502.552.602.652.702.762.812.872.932.99...3.063.133.193.263.333.413.493.573.653.74...3.833.924.014.114.224.334.444.564.674.80...4.935.075.215.365.525.685.846.006.166.33...6.526.726.947.167.387.627.868.108.368.62...8.909.209.509.8010.110.510.911.311.712.1...12.613.113.614.214.815.516.317.118.119.1...20.121.222.423.725.026.728.530.432.635.1...37.840.743.746.850.053.456.860.464.067.8...72.076.480.885.490.295.0100105110116...122128134140146152158165172180];

%轭部磁路校正曲线表格Bjz=[0.30.801.091.241.3951.551.711.86];

Xj=[00.0250.050.0750.100.1250.150.1750.200.2250.25];

%定子

Cj1=[0.6360.6360.6360.6370.6370.6370.6370.6370.6370.6370.6370.610.610.6120.6170.6190.6230.6270.6300.6340.6350.6360.550.5510.5530.5670.5720.5790.5910.6110.6340.6350.6260.450.4510.4530.4580.4620.4750.4850.5040.5230.5350.550.360.3650.370.3750.380.390.3980.4040.4130.4220.4380.300.300.3020.3030.3050.310.3130.3170.3250.330.340.270.270.270.270.270.2730.2750.280.2940.3050.3200.230.230.2320.2340.2350.2390.240.2520.2670.280.291];%转子

Cj2=[0.6360.6360.6360.6370.6370.6370.6370.6370.6370.6370.6370.610.610.6120.6170.6190.6230.6270.6300.6340.6350.6360.550.5510.5530.5670.5720.5790.5910.6110.6340.6350.6260.450.4510.4530.4580.4620.4750.4850.5040.5230.5350.550.360.3650.370.3750.380.390.3980.4040.4130.4220.438

0.300.300.3020.3030.3050.310.3150.3270.340.3580.3730.270.270.270.270.2710.2750.280.2950.3060.3260.360.230.230.230.230.2320.2370.2420.2550.2720.2910.31];

%%|程序计算部分|%

IN=PN*1e5/(m*UN*etaN*cosfaiN);%额定电流nN=60*f/p;%额定转速

TN=9.549*PN*1000/nN;%额定和转矩

%定、转子冲片相关计算tao1i=pi*Di1/2/p;

tao1o=pi*Di2o/2/p;%内外极距Lefi=La+2*gi;

Lefo=La+2*go;%电枢计算长度

mFe=rouFe*Lb*KFe*(D2o+deita)^2*1.0*1e-3;%硅钢片质量mmi=2*p*Ami*hMi*roum*0.001;

mmo=2*p*Amo*hMo*roum*0.001;%内外永磁体的质量t1i=pi*Di1/Qi;

t1o=pi*D1/Qo;%内外定子齿距bt11i=pi*(Di1+2*(hi01+hi12))/Qi-2*ri1;

bt11o=pi*(D1+2*(ho01+ho12))/Qo-1.36;%内外定子齿宽的计算bt12i=pi*(Di1+2*(hi01+his1))/Qi-bi1;bt12o=pi*(D1+2*(ho01+hos1))/Qo-bo1;ifbt12iaii=54.6/60;aio=52.3/60;%内外极弧系数Kfi=4*sin(aii*pi/2)/pi;

Kfo=4*sin(aio*pi/2)/pi;%内外气隙磁密波形系数Kfaii=8*sin(aii*pi/2)/aii/pi^2;

Kfaio=8*sin(aio*pi/2)/aio/pi^2;%内外气隙磁通波形系数Kg1i=t1i*(4.4*gi+0.75*bi01)/(t1i*(4.4*gi+0.75*bi01)-bi01^2);Kg1o=t1o*(4.4*go+0.75*bo01)/(t1o*(4.4*go+0.75*bo01)-bo01^2);Kg2i=1;Kg2o=1;

Kgi=Kg1i*Kg2i;Kgo=Kg1o*Kg2o;%内外气隙系数bm0i=0;bm0o=0;

whileabs((bm01i-bm0i)/bm01i)>0.01|abs((bm01o-bm0o)/bm01o)>0.001faig0i=bm01i*Br*Ami*1*1e-4/zata0i;

faig0o=bm01o*Br*Amo*1e-4/zata0o;%内外空载主磁通Bgi=10000*faig0i/aii/tao1i/Lefi;

Bgo=10000*faig0o/aio/tao1o/Lefo;%内外气隙磁密Fgi=2*Bgi*Kgi*gi*0.01/mu0;Fgqi=Fgi;

Fgo=2*Bgo*Kgo*go*0.01/mu0;Fgqo=Fgo;%内外交直轴气隙磁位差Bt1i=Bgi*t1i*Lefi/bt1i/KFe/L1;

Bt1o=Bgo*t1o*Lefo/bt1o/KFe/L1;%内外定子齿磁密Ht1i=interp1(Bcz,H,Bt1i,'spline');Ht1o=interp1(Bcz,H,Bt1o,'spline');%插值

Ft1i=2*ht1i*Ht1i;Ft1o=2*ht1o*Ht1o;%内外定子齿磁位差Bj1=(faig0i+faig0o)*1e4/2/L1/KFe/hj1;Bj1i=Bj1;Bj1o=Bj1;%定子轭部磁密bzj1=hj1/tao1o;%定子bzj2i=hj2i/tao1i;%内转子bzj2o=hj2o/tao1o;%外转子Hj1=interp1(Bcz,H,Bj1,'spline');C1=interp2(Xj,Bjz,Cj1,bzj1,Bj1);Fj1=2*C1*Hj1*Lj1;%定子轭部磁位差Bj2i=faig0i*1e4/2/L2/KFe/hj2i;

Bj2o=faig0o*1e4/2/L2/KFe/hj2o;%内外转子轭部磁位差Hj2i=interp1(Bcz,H,Bj2i,'spline');Hj2o=interp1(Bcz,H,Bj2o,'spline');%插值

C2i=interp2(Xj,Bjz,Cj2,bzj2i,Bj2i);C2o=interp2(Xj,Bjz,Cj2,bzj2o,Bj2o);Fj2i=2*Hj2i*C2i*Lj2i;

Fj2o=2*Hj2o*C2o*Lj2o;%内外转子轭部磁位差sumFi=Fgi+Ft1i+Fj1+Fj2i;

sumFo=Fgo+Ft1o+Fj1+Fj2o;%每对极总磁位差Ksti=(Fgqi+Ft1i)/Fgqi;

Ksto=(Fgqo+Ft1o)/Fgqo;%磁路齿饱和系数Agi=faig0i/sumFi;

Ago=faig0o/sumFo;%主磁导lemtgi=200*Agi*hMi/mur/mu0/Ami;

lemtgo=200*Ago*hMo/mur/mu0/Amo;%主磁导标幺值lemtni=zata0i*lemtgi;

lemtno=zata0o*lemtgo;%外磁路总磁导标幺值lemtzi=(zata0i-1)*lemtgi;

lemtzo=(zata0o-1)*lemtgo;%漏磁导标幺值bm0i=lemtni/(lemtni+1);

bm0o=lemtno/(lemtno+1);%永磁体空载工作点bm01i=(bm01i+bm0i)/2;

bm01o=(bm01o+bm0o)/2;%程序循环变量end

Bg1i=Kfi*faig0i*1e4/aii/tao1i/Lefi;

Bg1o=Kfo*faig0o*1e4/aio/tao1o/Lefo;%气隙磁密基波幅值E0i=4.44*f*Kdp*N*faig0i*Kfaii;

E0o=4.44*f*Kdp*N*faig0o*Kfaio;%空载反电动势

%参数计算R1=RouCu*2e-3*Lavi*N/(pi*a*Nt1*(d11^2/4));%定子直流电阻Cxi=4*pi*f*mu0*Lefi*(Kdp*N)^2*0.01/p;

Cxo=4*pi*f*mu0*Lefo*(Kdp*N)^2*0.01/p;%漏抗系数

KU1i=1;KU1o=1;KL1i=1;KL1o=1;%内外槽的上下部节距漏抗系数lemtU1i=hi01/bi01+2*his1/(bi01+bi1);lemtU1o=ho01/bo01+2*hos1/(bo01+bo1);rfai=bi1/bi2;rfao=bo1/bo2;

bt2i=hi2/bi2;bt2o=ho2/bo2;

Kr1i=1/3-(1-rfai)/4*(1/4+1/3/(1-rfai)+1/2/(1-rfai)^2+1/(1-rfai)^3+log(rfai)/(1-rfai)^4);

Kr1o=1/3-(1-rfao)/4*(1/4+1/3/(1-rfao)+1/2/(1-rfao)^2+1/(1-rfao)^3+log(rfao)/(1-rfao)^4);

Kr2i=(2*pi^2-9*pi)/1536/bt2i^2+pi/16/bt2i-pi/8/(1-rfai)/bt2i-(pi^2/64/(1-rfai)/bt2i^2+pi/8/(1-rfai)^2/bt2i)*log(rfai);

Kr2o=(2*pi^2-9*pi)/1536/bt2o^2+pi/16/bt2o-pi/8/(1-rfao)/bt2o-(pi^2/64/(1-rfao)/bt2o^2+pi/8/(1-rfao)^2/bt2o)*log(rfao);

lemtL1i=bt2i*(Kr1i+Kr2i)/(pi/8/bt2i+(1+rfai)/2)^2;lemtL1o=bt2o*(Kr1o+Kr2o)/(pi/8/bt2o+(1+rfao)/2)^2;lemts1i=KU1i*lemtU1i+KL1i*lemtL1i;

lemts1o=KU1o*lemtU1o+KL1o*lemtL1o;%定子槽比漏磁导Xs1i=2*p*m*L1*lemts1i*Cxi/Lefi/Kdpi^2/Qi;

Xs1o=2*p*m*L1*lemts1o*Cxo/Lefo/Kdpo^2/Qo;%定子槽漏抗sums=0.0156;%查表得

Xd1i=m*tao1i*sums*Cxi/pi^2/Kgi/gi/Kdpi^2/Ksti;

Xd1o=m*tao1o*sums*Cxo/pi^2/Kgo/gi/Kdpo^2/Ksto;%定子谐波漏抗XE1i=0.2*(LEi/Lefi*Kdpi^2)*Cxi;

XE1o=0.2*(LEo/Lefo*Kdpo^2)*Cxo;%定子端部漏抗X1i=Xs1i+Xd1i+XE1i;

X1o=Xs1o+Xd1i+XE1o;%定子漏抗

Kadi=1/Kfi;Kado=1/Kfo;%直轴电枢折算系数Kaqi=Kadi;Kaqo=Kado;%交轴电枢磁动势折算系数Idi=IN/2;Ido=IN/2;

Fadi=0.45*m*Kadi*Kdpi*N*Idi/p;Fado=0.45*m*Kado*Kdpo*N*Ido/p;fa1i=2*Fadi/zata0i/hMi/Hc/10;fa1o=2*Fado/zata0o/hMo/Hc/10;bmNi=lemtni*(1-fa1i)/(lemtni+1);bmNo=lemtno*(1-fa1o)/(lemtno+1);

faigNi=(bmNi-(1-bmNi)*lemtzi)*Ami*Br*1e-4;faigNo=(bmNo-(1-bmNo)*lemtzo)*Amo*Br*1e-4;

Edi=4.44*f*Kdpi*N*faigNi*Kfaii;Edo=4.44*f*Kdpo*N*faigNo*Kfaio;

Xadi=abs(E0i-Edi)/Idi;Xado=abs(E0o-Edo)/Ido;%直轴电枢反应电抗Xdi=Xadi+X1i;Xdo=Xado+X1o;%直轴同步电抗

%工作特性计算pfw=150;%机械损耗

catai=pi/9;catao=catai;%转矩角Xqi=Xdi;Xqo=Xdo;

Xaqi=Xadi;Xaqo=Xado;%交轴同步电抗

P1i=m*(E0i*UN*(Xqi*sin(catai)-R1*cos(catai))+R1*UN^2)/(Xdi^2+R1^2);P1o=m*(E0o*UN*(Xqo*sin(catao)-R1*cos(catao))+R1*UN^2)/(Xdo^2+R1^2;P1=P1i+P1o;%输入功率

Idi=(R1*UN*sin(catai)+Xqi*(E0i-UN*cos(catai)))/(Xdi^2+R1^2);

Ido=(R1*UN*sin(catao)+Xqo*(E0o-UN*cos(catao)))/(Xdo^2+R1^2);%直轴电流

Iqi=(Xdi*UN*sin(catai)-R1*(E0i-UN*cos(catai)))/(Xdi^2+R1^2);

Iqo=(Xdo*UN*sin(catao)-R1*(E0o-UN*cos(catao)))/(Xdo^2+R1^2);%交轴电流cosfaii=cos(catai-atan(Idi/Iqi));

cosfaio=cos(catao-atan(Ido/Iqo));%功率因数I1i=sqrt(Idi^2+Iqi^2);I1o=sqrt(Ido^2+Iqo^2);I1=(I1i+I1o)/2;%定子电流pCu=m*I1^2*R1;%定子电阻损耗Egi=sqrt((E0i-Idi*Xadi)^2+(Iqi*Xaqi)^2);Ego=sqrt((E0o-Ido*Xado)^2+(Iqo*Xaqo)^2);faigi=Egi/4.44/f/Kdpi/N/Kfaii;

faigo=Ego/4.44/f/Kdpo/N/Kfaio;%负载气隙磁通Bgdi=faigi*1e4/aii/tao1i/Lefi;

Bgdo=faigo*1e4/aio/tao1o/Lefo;%负载气隙磁密Bt1di=Bgdi*t1i*Lefi/bt1i/KFe/L2;

Bt1do=Bgdo*t1o*Lefo/bt1o/KFe/L2;%负载定子齿磁密Bj1d=(faigi+faigo)*1e4/2/L1/KFe/hj1;%负载定子轭磁密pt1di=interp1(Bcz,H,Bt1di,'spline');pt1do=interp1(Bcz,H,Bt1do,'spline');pj1d=interp1(Bcz,H,Bj1d,'spline');k1i=2.5;k1o=k1i;k2=2;

pFe=(k1i*pt1di*Vt1i+k1o*pt1do*Vt1o+k2*pj1d*Vj1)*1e-3;%铁耗psNxing=0.015;%经验数据

ps=(I1/IN)^2*psNxing*PN*1e3;%杂散损耗sump=pCu+pFe+pfw+ps;%总损耗P2=P1-sump;%输出功率eta=(P2/P1)*100;%效率

%一系列递增的转矩角对应的工作特性循环catai1=pi/18:pi/36:2*pi/3;catao1=catai1;Xqi=Xdi;Xqo=Xdo;Xaqi=Xadi;Xaqo=Xado;

P1i1=m*(E0i*UN*(Xqi*sin(catai1)-R1*cos(catai1))+R1*UN^2)/(Xdi^2+R1^2);P1o1=m*(E0o*UN*(Xqo*sin(catao1)-R1*cos(catao1))+R1*UN^2)/(Xdo^2+R1^2);

P11=P1i1+P1o1;

Idi1=(R1*UN*sin(catai1)+Xqi*(E0i-UN*cos(catai1)))/(Xdi^2+R1^2);Ido1=(R1*UN*sin(catao1)+Xqo*(E0o-UN*cos(catao1)))/(Xdo^2+R1^2);Iqi1=(Xdi*UN*sin(catai1)-R1*(E0i-UN*cos(catai1)))/(Xdi^2+R1^2);Iqo1=(Xdo*UN*sin(catao1)-R1*(E0o-UN*cos(catao1)))/(Xdo^2+R1^2);cosfaii1=cos(catai1-atan(Idi1/Iqi1));cosfaio1=cos(catao-atan(Ido1/Iqo1));I1i1=sqrt(Idi1.^2+Iqi1.^2);I1o1=sqrt(Ido1.^2+Iqo1.^2);I11=(I1i1+I1o1)/2;pCu1=m*I11.^2*R1;

Egi1=sqrt((E0i-Idi1*Xadi).^2+(Iqi1*Xaqi).^2);Ego1=sqrt((E0o-Ido1*Xado).^2+(Iqo1*Xaqo).^2);faigi1=Egi1/4.44/f/Kdpi/N/Kfaii;faigo1=Ego1/4.44/f/Kdpo/N/Kfaio;Bgdi1=faigi1*1e4/aii/tao1i/Lefi;Bgdo1=faigo1*1e4/aio/tao1o/Lefo;Bt1di1=Bgdi1*t1i*Lefi/bt1i/KFe/L2;Bt1do1=Bgdo1*t1o*Lefo/bt1o/KFe/L2;Bj1d