Magnetic Circuit Matlab Script: Difference between revisions

From Class Wiki
Jump to navigation Jump to search
No edit summary
(Error correction)
 
(3 intermediate revisions by the same user not shown)
Line 1: Line 1:
<font face="Courier New">
<font face="Courier New">

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

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

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

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

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

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

%&nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp _/&nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &nbsp &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
% John Hawkins


Line 26: Line 9:


% January 10, 2009
% January 10, 2009
<br /> <br />

% Problem Statement Constraints
% Problem Statement Constraints


Line 33: Line 16:
Bfg=.1;
Bfg=.1;


% Assumed relative permeability
<br />% Assumed relative permeability


mu0=4e-7*pi;
mu0=4e-7*pi;
Line 41: Line 24:
mu=mur*mu0
mu=mur*mu0


% Geometry
<br />% Geometry


lfg=1e-2
lfg=1e-2
Line 51: Line 34:
ldc=48e-2
ldc=48e-2


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


Afg=4e-2*8e-2
Afg=4e-2*8e-2
Line 63: Line 46:
Adabc=2e-2*4e-2
Adabc=2e-2*4e-2


% Air gap
<br />% Air gap


Rfg=lfg/(mu*Afg)
Rfg=lfg/(mu*Afg)
Line 71: Line 54:
Ffg=Rfg*Ifg
Ffg=Rfg*Ifg


% Right arms
<br />% Right arms


Idef=Ifg
Idef=Ifg
Line 85: Line 68:
Fghc=Fdef;
Fghc=Fdef;


% Center Column
<br />% Center Column


Fdc=Fdef+Ffg+Fghc
Fdc=Fdef+Ffg+Fghc
Line 93: Line 76:
Idc=Fdc/Rdc
Idc=Fdc/Rdc


% Left arm
<br />% Left arm


Idabc=Idc-Idef
Idabc=Idc-Idef
Line 101: Line 84:
Fdabc=Rdabc*Idabc
Fdabc=Rdabc*Idabc


% Results
<br />% Results


Ftotal = Fdabc+Fdc+Fdef+Ffg+Fghc
Ftotal = Fdabc+Fdc+Fdef+Ffg+Fghc

Latest revision as of 10:48, 18 January 2010

% 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=45e-2+48e-2+45e-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