Magnetic Circuit Matlab Script

%&nbsp &nbsp _-_,,&nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp ,,&nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp _-_-&nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp ,,&nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp

%&nbsp (&nbsp &nbsp //&nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp ||&nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp /,&nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp _&nbsp &nbsp &nbsp ;&nbsp &nbsp &nbsp &nbsp &nbsp &nbsp ||&nbsp &nbsp &nbsp &nbsp '&nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp

%&nbsp &nbsp &nbsp _||&nbsp &nbsp &nbsp /'\\&nbsp ||/\\&nbsp \\/\\&nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp ||&nbsp __&nbsp &nbsp &nbsp &nbsp <&nbsp \,&nbsp \\/\/\&nbsp ||/\&nbsp \\&nbsp \\/\\&nbsp _-_,&nbsp

%&nbsp &nbsp &nbsp _||&nbsp &nbsp ||&nbsp ||&nbsp ||&nbsp ||&nbsp ||&nbsp ||&nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp ~||-&nbsp &nbsp -&nbsp &nbsp &nbsp /-||&nbsp ||&nbsp |&nbsp |&nbsp ||_<&nbsp ||&nbsp ||&nbsp ||&nbsp ||_.&nbsp &nbsp

%&nbsp &nbsp &nbsp &nbsp ||&nbsp &nbsp ||&nbsp ||&nbsp ||&nbsp ||&nbsp ||&nbsp ||&nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp ||===||&nbsp ((&nbsp ||&nbsp ||&nbsp |&nbsp |&nbsp ||&nbsp |&nbsp ||&nbsp ||&nbsp ||&nbsp ~&nbsp ||&nbsp

%&nbsp -__-,&nbsp &nbsp \\,/&nbsp &nbsp \\&nbsp |/&nbsp \\&nbsp \\&nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp (&nbsp \_,&nbsp |&nbsp &nbsp &nbsp \/\\&nbsp \\/\\/&nbsp \\,\&nbsp \\&nbsp \\&nbsp \\&nbsp ,-_-&nbsp &nbsp

%&nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp _/&nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp `&nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp

%

% John Hawkins

% EMEC

% Magnetic Circuit Calculations

% Professor Rob Frohne

% January 10, 2009

% Problem Statement Constraints

N=1600;

Bfg=.1;

% Assumed relative permeability

mu0=4e-7*pi;

mur=500;

mu=mur*mu0

% Geometry

lfg=1e-2

ldef=31e-2+24.5e-2

lghc=ldef

ldc=48e-2

ldabc=46e-2+48e-2+46e-2

Afg=4e-2*8e-2

Adef=Afg

Aghc=Afg

Adc=8e-2*12e-2

Adabc=2e-2*4e-2

% Air gap

Rfg=lfg/(mu*Afg)

Ifg=Bfg*Afg

Ffg=Rfg*Ifg

% Right arms

Idef=Ifg

Rdef=ldef/(mu*Adef)

Fdef=Rdef*Idef

Ighc=Idef;

Rghc=Rdef;

Fghc=Fdef;

% Center Column

Fdc=Fdef+Ffg+Fghc

Rdc=ldc/(mu*Adc)

Idc=Fdc/Rdc

% Left arm

Idabc=Idc-Idef

Rdabc=ldabc/(mu*Adabc)

Fdabc=Rdabc*Idabc

% Results

Ftotal = Fdabc+Fdc+Fdef+Ffg+Fghc

i=Ftotal/N