Transformer optimization calculation sheet

( MATLAB 5.3.0.620a (R11) – Student Edition )

Подпись: 0 1 2 1 2 2
Transformer optimization calculation sheet Подпись: %Base model: 18mOhm

% Experimental data used (Efficiency x power) % File: n48v. mat

%|-

– Experimental………

%

Pout

Efficiency[%]

%

60

76.27

%

120

86.36

%

238

92.3

%

297

93.52

%

355

94.31

%

449

95.08

%

561

95.59

%

673

95.84

%

783

95.95

%

945

96

%

1057

95.93

%

1167

96.01

%

1275

95.86

%

1359

95.72

%

1463

95.53

%

1565

95.33

%

1666

95.12

%

1917

94.7

%

2015

94.48

%

2112

94.24

%

2208

94

%

2283

94.01

%

2377

93.76

%

2469

93.49

%

2560

93.23

%

2631

93.02

%

2808

92.51

%

2963

92.07

%

3318

91.24

%

3475

90.73

for i=1:1:70, for j=1:1:50, Ax(50*(i-1)+j)=A(i); end; end;

%……………………

load – ascii n48v. mat NL = 18.6;

nc_schmid = 100*n48v(:,1)./(((n48v(:,1))./(0.01*n48v(:,2)))-NL);

polschmid=polyfit(n48v(:,1),nc_schmid,12); % EXPERIMENTAL DATA

FitSchmid = polyval(polschmid, n48v(:,1));

nmedschmid=0;

Area=0; for i=1:1:3500,

fitncschmid(i) = polyval(polschmid, i);

nmedschmid = nmedschmid + fitncschmid(i)*Ax(i);

Area = Area + Ax(i);

end;

nmedschmid = nmedschmid / Area %figure(1); hold off;

%plot(n48v(:,1),n48v(:,2),’r’);

%hold;

%plot(n48v(:,1),nc_schmid,’b’);

%plot( 10:1:3499,fitncschmid( 10:3499),’g’); %grid;

%

%——————————

load – ascii trafo. mat

polmag = polyfit(trafo(:,4),trafo(:,6),5); Fitmag = polyval(polmag,10:1:20); figure(1); hold off; plot(10:1:20,Fitmag,’g’); hold;

%plot(trafo(:,4),trafo(:,6),’*r’);

%axis([10 20 0.010 0.055]);

hold;

grid;

poltra = polyfit(trafo(:,1),trafo(:,4),5);

Fittra = polyval(poltra,35:1:54); figure(2); hold off; plot(35:1:54,Fittra,’g’); hold;

plot(trafo(:,1),trafo(:,4),’vr’);

hold;

grid;

NL = 18.6;

R = polyval(polmag, NL);

ntrafo186 = (P-R.*((P. A2)./(39.2A2)))./P;

figure(3); hold off;

plot(P,100*ntrafo186,’g’);

hold;

ninv186 = (P-(10.6)-R.*((P. A2)./(39.2A2)))./P; plot(P(100:3500),100*ninv186(100:3500),’r’); plot(n48v(:,1),n48v(:,2),’b’);

 

% Estimated transformer Req[Ohm] x PNL[W] % Fit Curve

 

% Inverter Characteristic PNL[W] x Vin[V] % Fit Curve

 

% Nominal no-load % ReqxV characteristic

% Rated Trafo Conversion Efficiency x power curve % Transformer conversion effic. x power curve (Model) % Transformer rated efficiency x power curve (Model)

% Inverter rated efficiency x power curve (Experimental)

 

%NL = 12.0; for j=1:1:100,

NL(j)=j/10+10; % 10W a 20W R(j) = polyval(polmag, NL(j)); ntrafox = (P-R(j).*((P. A2)./(39.2A2)))./P; nw(j)=0;

Area=0; for i=1:1:3500,

nw(j) = nw(j) + (fitncschmid(i)*ntrafox(i)/ntrafo186(i))*Ax(i); Area = Area + Ax(i); end;

nw(j) = nw(j) / Area; end;

figure(4); hold off; plot(NL, nw,’r’); axis([10 20 97.4 98.5]); hold;

%plot(NL(j),nw(j),’r’);

grid;

%…………………………………………..

EE=200;

Pw = 0.72.*NL + ((1-nw/100)./(nw/100)).*EE; [x, y]=min(Pw); nl1=NL(y)

figure(5); hold off;

plot(NL, Pw,’r’);

hold;

EE=350;

Pw = 0.72.*NL + ((1-nw/100)./(nw/100)).*EE; [x, y]=min(Pw); nl2=NL(y) plot(NL, Pw,’b’);

EE=500;

Pw = 0.72.*NL + ((1-nw/100)./(nw/100)).*EE; [x, y]=min(Pw); nl3=NL(y)

plot(NL, Pw,’g’);

grid;

%…………………………………………..

figure(6); hold off;

plot(n48v(:,1),n48v(:,2),’r’);

hold;

axis([0 3000 80 97]);

%plot( 10:1:3499,fitncschmid( 10:3499),’g’); grid;

R1 = polyval(polmag, nl1)

ntrafox = (P-R1*((P. A2)./(39.2A2)))./P;

ncw=0;

Area=0; for i=1:1:3500,

nct1(i) = (fitncschmid(i)*ntrafox(i)/ntrafo186(i)); nt1(i) = 100*(i/(nl1 + 100*i/nct1(i)));

ncw = ncw + nct1(i)*Ax(i);

Area = Area + Ax(i);

end;

nmed1 = ncw / Area

R2 = polyval(polmag, nl2)

ntrafox = (P-R2*((P. A2)./(39.2A2)))./P;

ncw=0;

Area=0; for i=1:1:3500,

nct2(i) = (fitncschmid(i)*ntrafox(i)/ntrafo186(i)); nt2(i) = 100*(i/(nl2 + 100*i/nct2(i)));

ncw = ncw + nct2(i)*Ax(i);

Area = Area + Ax(i); end;

nmed2 = ncw / Area

R3 = polyval(polmag, nl3)

ntrafox = (P-R3*((P. A2)./(39.2A2)))./P;

ncw=0;

Area=0; for i=1:1:3500,

nct3(i) = (fitncschmid(i)*ntrafox(i)/ntrafo186(i)); nt3(i) = 100*(i/(nl3 + 100*i/nct3(i))); ncw = ncw + nct3(i)*Ax(i);

Area = Area + Ax(i); end;

nmed3 = ncw / Area

%plot(10:3499,nct1(10:3499),’k’); plot(50:3499,nt1(50:3499),’k’); %plot( 10:3499,nct2( 10:3499),’k’); plot(50:3499,nt2(50:3499),’k’); %plot( 10:3499,nct3( 10:3499),’k’); plot(50:3499,nt3(50:3499),’k’);

%………………………

Подпись: NL1=18.6;Подпись: NL2=16; NL3=10; n1 = 0.984; n2 = 0.969; n3 = 0.963;

n4 = nmedl/100; NL4=nl1; n5 = nmed2/100; NL5=nl2; n6 = nmed3/100; NL6=nl3;

P1 = 0.72*NL1 + ((1-n1)/n1).*E; P2 = 0.72*NL2 + ((1-n2)/n2).*E;

P3 = 0.72*NL3 + ((1-n3)/n3).*E; P4 = 0.72*NL4 + ((1-n4)/n4).*E; P5 = 0.72*NL5 + ((1-n5)/n5).*E; P6 = 0.72*NL6 + ((1-n6)/n6).*E;

figure(7); hold off; plot(E, P1,’r’); axis([0 500 6 20]); hold;

plot(E, P2,’b’);

plot(E, P3,’g’);

plot(E, P4,’k’);

plot(E, P5,’k’);

plot(E, P6,’k’);

grid;

%———— END —

Leave a reply

You may use these HTML tags and attributes: <a href="" title=""> <abbr title=""> <acronym title=""> <b> <blockquote cite=""> <cite> <code> <del datetime=""> <em> <i> <q cite=""> <s> <strike> <strong>