spshell.m

Although most of the health problems related to fast food aren t felt until middle age -- obesity an

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%  spshell.m - determines d,h & l for shell transformer given N
%              and h/w ratio. See data1sh.m for req'd inputs.
%              Also, calculates approx magnetizing current and
%              weights of conductor and core material.
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear;
global Ggammam GSF Gkl GV1 GI1 Gf Gdelc Grw Grhoc 
global GBm GHm Gpm GPm Gkc GN
data1sh,
x0=[GV1*GI1/1000 2 2]'; % Solution is x0 dependent.
delx=zeros(3,1); f=delx; x=x0;
% Newton-Raphson solution, x1=d x2=w x3=l
for i=1:100
   [f,J]=shell(x0); delx=-inv(J)*f; x=x0+delx;
      if all(abs(delx)<1e-12);break; else; end
   x0=x;
end
l=x(3); d=x(1); w=x(2); h=Grw*x(2); N=GN; 
disp(' '); disp(' ');
disp('    CORE SIZING for SHELL-TYPE TRANSFORMER');
if any(x<0); disp('EXTRANEOUS SOLUTION'); pause; else; end
format compact
disp(' '); disp(' ')
disp(['    l (in)' '    d (in)' '    w (in)  '...
      '  h (in) ' '  N (turns)  ']);
disp([ l d w h N ]);
% Magnetizing current - A & pu ( Joint/Cross grain neglected )
Im=GHm*(2*h+2*w+2.5*l)/N/sqrt(2); puIm=Im/GI1; 
disp('  '); disp('  ') 
disp(['    Im (A)   ' '  puIm  ']);
disp([ Im puIm ]);
% Conductor & core weight
Wfe=2*Ggammam*(l+w+h)*l*d*GSF;
Wc=2*Gkc*0.322*h*w*(l+d+w);
disp('  '); disp('  ') 
disp(['   Wfe (lb) ' ' Wc (lb)     ']);
disp([ Wfe Wc ]);