Magnetic Circuit Matlab Script

% _-_,, ,, _-_- ,,

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

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

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

% || || || || || || || ||===|| (( || || | | || | || || || ~ ||

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

% _/ `

%

% 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