Program Matlab
Pn=76000 %puterea nominala,la ax,in W
p=2
f1n=50
cosfi1n=0.95 % capacitiv,cadranul 2!!
%fi1n=acos(cosfi1n)% fi1n >0,pt cosfin capacitiv, vezi tema proiect,
atentie!!!
fi1n=-acos(cosfi1n)% fi1n <0,pt cosfin inductiv!!Pt cei cu cosfi1n inductiv!!
sinfi1n=sin(fi1n)
rann=.92 % randamentul nominal
efn=cosfi1n*rann
Mm=1.8 % cuplul maxim,raportat la cuplul nominal, data de catalog (Mm=1.8
[Link])
Sn=Pn/(rann*cosfi1n)% puterea aparenta ,nominala,la bornele inf.
trifazate(stator)
sumapn=(1-rann)*cosfi1n %suma pierderilor,in u.r.,raportate la Sn
prefn=sumapn/3 %pierderi de referinta pt pj1n=pjen=pfen+pmecn
'Date calculate'
mm=Mm*efn % cuplul maxim ,raportat la cuplul de baza (Mb=Sn/(2*pi*n))
%ien=mm % curent de excitatie,in u.r., nominal,din triunghiul simplificat al
curentilor
ienc=Mm*(1+efn)/2 %curent de excitatie nominal,corectat >ien
imnc=sqrt(ienc^2-cosfi1n^2)-sinfi1n%curent de magnetizare,din triunghiul
curentilor
xmnc=1/imnc % estimare reactanta de magnetizare
tetanc=atan(cosfi1n/(imnc+sinfi1n))%unghi intern nominal,in rad,din triunghiul
curentilor
tetagnc=tetanc*180/pi %unghi intern in grade
% observatie: se va construi,la scara,DIAGRAMA CURENTILOR(cadranul 1,pt
cosfin inductiv, respectiv in cadranul 2 pt. cosfin capacitiv)
'parametri schemei echivalente'
pj1nc=sumapn/3 % in ipoteza pierderilor uniform distribuite:pj1=pje=pmec+pfe
ponc=sumapn-pj1nc % ponc=pmec+pfe1+pjen =const
r1n=pj1nc %rezistenta,in u.r., [Link] stator
xs1n=0.1*xmnc %se estimeaza reactanta de scapari a inf. stator
z1n=r1n+j*xs1n;
u1z=1; %origine de faza,diagrama curentilor (motor alimentat de la o sursa
de tensiune)
i1zc=exp(j*fi1n);%Atentie!! fi1n >0, in regim capacitiv,si <0 in [Link]
umz=u1z-z1n*i1zc;
umc=abs(umz)
xmn=umc/imnc %valoare, corectata, reactanta magnetizare
rwn=umc^2/ponc %rezistenta echivalenta pierderilor constante
zmn=j*xmn*rwn/(rwn+j*xmn);
'caracteristici electromecanice din schema in T(ie=const.)'
tetag=[1:tetagnc/20:1.4*tetagnc];%pt. caracteristici electromecanice ,gol-
nominal
%tetag=[2:tetagnc/10:180];%se activeaza [Link] unghiulara
extinsa,0.180 grade
teta=tetag.*pi/180;
iez=ienc*exp(-j*(pi/2+teta));%vezi diagrama curentilor!
i1z=(u1z-zmn.*iez)./(z1n+zmn);%vezi schema echivalenta paralel
umz=u1z-z1n.*i1z;um=abs(umz);
u1=abs(u1z);
i1=abs(i1z);
fi1=angle(i1z);
alfa=i1.*cos(fi1)/cosfi1n;% alfa =p1/p1n,pt i1=1 si cosfi1=cosfi1n alfa=1
p1=(u1.*i1).*cos(fi1); %puterea activa,electrica, absorbita de la retea
q1=(u1.*i1).*sin(fi1); %puterea reactiva schimbata cu reteaua
p2=p1-r1n.*i1.^2-um.^2/rwn;%puterea mecanica, disponibila la ax
ran=p2./p1; %randamentul motorului
Me=p2./efn; %cuplul, la ax,raportat la cuplul nominal, in u.r.
%plot(alfa,i1) %,fig.1, pt alfa =1 (sarcina nominala), i1nx=.... fata de
i1n=1
%plot(alfa,cos(fi1)) %fig.2, pt alfa =1,cosfi1nx=.....,fata de cosfi1n=...
%plot(alfa,ran) % fig.3, pt alfa =1 ,rannx=....,fata de rann=
%plot(alfa,ran,'-',alfa,i1,'+',alfa,cos(fi1),'*')
Mmi=(Mm.*tetag)./tetag; %cuplul maxim impus
plot(tetag,Me,'-',tetag,Mmi,'x',tetag,i1,'+',tetag,cos(fi1),'.')
grid
Rezultate Matlab
Pn = 76000
p =2
f1n = 50
cosfi1n = 0.9500
fi1n = -0.3176
sinfi1n =-0.3122
rann = 0.9200
efn =0.8740
Mm =1.8000
Sn =8.6957e+004
sumapn = 0.0760
prefn = 0.0253
Date calculate
mm =1.5732
ienc =1.6866
imnc = 1.7058
xmnc = 0.5862
tetanc = 0.5983
tetagnc =34.2818
parametri schemei echivalente
pj1nc = 0.0253
ponc =0.0507
r1n = 0.0253
xs1n =0.0586
umc = 0.9588
xmn = 0.5621
rwn =18.1448
caracteristici electromecanice din schema in T(ie=const.)
Grafice