Xcos examples from textbooks and for blocksHEADmaster

3 files changed, 93 insertions, 0 deletions

diff --git a/Working_Examples/293/CH15/EX15.1/eg15_1.sce b/Working_Examples/293/CH15/EX15.1/eg15_1.sce
new file mode 100755
index 0000000..2417766
--- /dev/null
+++ b/Working_Examples/293/CH15/EX15.1/eg15_1.sce
@@ -0,0 +1,42 @@
+//a

+phi = 6*10^-4; //given magnetic flux (in Wb)

+A = 0.001; // cross sectional area (in meter square)

+B = phi/A ; //

+Ha = 10; //magnetic field intensity of material a needed to establish   the given magnetic flux 

+Hb = 77; // magnetic field intensity of material b

+Hc = 270; // magnetic field intensity of material c

+La = 0.3; //arc length of material a (in meters)

+Lb = 0.2; //arc length of material b (in meters) 

+Lc = 0.1; //arc length of material c (in meters)

+

+F = Ha*La + Hb*Lb + Hc*Lc; //magnetomotive force

+disp("a")

+disp(F, "magnetomotive force needed to establish a flux of 6*10^-4(in At) = ")

+

+//b

+N = 100; //no. of turns 

+I = F/N; //current in amps

+disp("b")

+disp(I,"current that must be made to flow through the coil(in amps) = ")

+

+//c

+MU0 = 4*%pi*10^-7; 

+MUa = B/Ha; //permeability of material a

+MUb = B/Hb; //permeability of material b

+MUc = B/Hc; //permeability of material c

+

+MUra = MUa/MU0; //relative permeability of material a

+MUrb = MUb/MU0; //relative permeability of material b

+MUrc = MUc/MU0; //relative permeability of material c

+

+Ra = Ha*La/phi; //reluctance of material a 

+Rb = Hb*Lb/phi; //reluctance of material b

+Rc = Hc*Lc/phi; //reluctance of material c

+

+disp("c")

+disp(MUra,"relative permeability of material a = ")

+disp(MUrb,"relative permeability of material b = ")

+disp(MUrc,"relative permeability of material c = ")

+disp(Ra,"reluctance of material a = ")

+disp(Rb,"reluctance of material b = ")

+disp(Rc,"reluctance of material c = ")
\ No newline at end of file
diff --git a/Working_Examples/293/CH15/EX15.3/eg15_3.sce b/Working_Examples/293/CH15/EX15.3/eg15_3.sce
new file mode 100755
index 0000000..db0e188
--- /dev/null
+++ b/Working_Examples/293/CH15/EX15.3/eg15_3.sce
@@ -0,0 +1,39 @@
+mu0 = 4*%pi*10^-7;

+A = 0.0025; //cross sectional area of the coil

+//dimensions of the coil (in meters)

+Lg = 0.002; //air gap length (in meters)

+Lbd = 0.025; 

+Lde = 0.1;

+Lef = 0.025;

+Lfk = 0.2;

+Lbc = 0.175;

+Lcab = 0.5;

+

+Lbghc = 2*(Lbd + Lde + Lef + (Lfk/2)) - Lg;//length of the ferromagnetic material involved here

+

+phig = 4*10^-4; //air gap flux (in Wb)

+Bg = phig/A ; //air gap flux density (in tesla)

+Hg = Bg/mu0 ; //feild intensity of the air gap 

+mmfg = Hg*Lg ; //mmf produced in the air gap (in At)

+

+Bbc = 1.38 ; //flux density corresponding to cast steel

+

+Hbghc = 125; //field intensity corresponding to flux density of 0.16T in the steel

+mmfbghc = Hbghc*Lbghc ; // mmf corresponding to bghc

+

+mmfbc = mmfg + mmfbghc ; //mmf across path bc

+Hbc = mmfbc/Lbc;

+phibc = Bbc*A ; //flux produced in bc  

+

+phicab = phig + phibc; //total fiux existing in leg cab 

+Bcab = phicab/0.00375; //flux density

+Hcab = 690; 

+mmfcab = Hcab*Lcab; //mmf in leg cab

+

+mmf = mmfbc + mmfcab ; //mmf produced by the coil

+

+disp(mmf,"mmf produced by the coil(in At) = ")

+

+

+

+

diff --git a/Working_Examples/293/CH15/EX15.5/eg15_5.sce b/Working_Examples/293/CH15/EX15.5/eg15_5.sce
new file mode 100755
index 0000000..151ccb4
--- /dev/null
+++ b/Working_Examples/293/CH15/EX15.5/eg15_5.sce
@@ -0,0 +1,12 @@
+//b

+mu0 = 4*%pi*10^-7 ;

+//plunger magnet dimensions (in meters)

+x = 0.025; 

+h = 0.05;

+a = 0.025;

+g = 0.00125; 

+

+mmf = 1414; //(in At)

+

+F = %pi*a*mu0*(mmf^2)*(h^2)*(1/(x + h)^2)/g; //magnitude of the force

+disp(F, "magnitude of the force (in Newtons) = ")
\ No newline at end of file