cckt1.m

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

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%  cckt1.m - Analysis of C-shaped magnetic circuit.
%            Calculates coil current for specified
%            value of flux & plots Flux vs. mmf.
%            (Type 1 problem)
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear;
% Units of input dimensions - 'MKS' or 'FPS' (British)
dim='MKS';
if dim=='MKS'; k=1; elseif dim=='FPS'; k=39.37; else; end
w=0.15; h=0.10;      % Window width & height
d=0.05; del=0.002;   % Core depth & air gap length
l1=0.05; l2=0.05; l3=0.05; l4=0.05; % Core widths
phi=0.0040;          % Flux for analysis point - Wb
N=500; SF=0.95;      % Coil turns, stacking factor
w=w/k; h=h/k; d=d/k; del=del/k;
l=[ l1 l2 l3 l4]/k;  % All dimensions now in MKS units
lm=[h+(l(2)+l(4))/2 w+(l(1)+l(3))/2 h+(l(2)+l(4))/2-del ...
      w+(l(1)+l(3))/2];
Pg=4e-7*pi*(l(3)*d/del+0.52*(l(3)+d)+0.308*del); F1=0;
% Specified point analysis
for i=1:4
   F1=F1+hm27(phi/l(i)/d/SF)*lm(i);
end
F1=F1+phi/Pg; I1=F1/N;
disp([blanks(4) 'C-SHAPED MAGNETIC CIRCUIT, No Leakage']);
disp(' '); disp([blanks(8) 'COIL CURRENT = ', num2str(I1)]);
% Flux vs. mmf plot
phimax=2.0*d*min(l)*SF;
phip=linspace(0,phimax,100); npts=length(phip);
for j=1:npts
   F(j)=0;
   for i=1:4;
      F(j)=F(j)+hm27(phip(j)/l(i)/d/SF)*lm(i);
   end
   F(j)=F(j)+phip(j)/Pg;
end
plot(F,phip,F1,phi,'o'); grid
title('C-shaped magnetic circuit');
xlabel('Coil mmf, A-t'); ylabel('Flux, Wb');
legend('Flux vs mmf', 'Specified pt.', 4);