initial commit / add all books

15 files changed, 284 insertions, 0 deletions

diff --git a/1544/CH4/EX4.1/Ch04Ex1.sce b/1544/CH4/EX4.1/Ch04Ex1.sce
new file mode 100755
index 000000000..ba2b63a3c
--- /dev/null
+++ b/1544/CH4/EX4.1/Ch04Ex1.sce
@@ -0,0 +1,9 @@
+// Scilab code Ex4.1: Pg 116 (2008)

+clc; clear;

+A = 6e-04;                  // Cross-sectional area of pole face, metre-square

+phi = 30e-06;              // Flux, Wb

+B = phi/A;                  // Flux density, T

+printf("\nThe flux desity at the pole face = %2d mT", B/1e-03);

+

+// Result

+// The flux desity at the pole face = 50 mT

diff --git a/1544/CH4/EX4.10/Ch04Ex10.sce b/1544/CH4/EX4.10/Ch04Ex10.sce
new file mode 100755
index 000000000..2140b3f3c
--- /dev/null
+++ b/1544/CH4/EX4.10/Ch04Ex10.sce
@@ -0,0 +1,24 @@
+// Scilab code Ex4.10: Pg 125 (2008)

+clc; clear;

+N = 1000;                                     // Number of turns on coil

+r = 0.1;                                      // Mean radius of toroid, m

+phi = 0.1775e-03;                             // Flux density(value from graph), Wb

+A = %pi*1e-04;                               // Csa of toroid, m^2

+H = 88;                                       // Magnetic field strength(value from graph), At/m

+B = phi/A;                                    // Flux density, T

+

+// Part (a)

+l = 2*%pi*r;                                  // Effective length of toroid, m

+// Since H = (N*I)/l, solving for I

+I = (H*l)/N ;                                 // Electric current in coil, A

+printf("\nCoil current = %4.1f mA", I/1e-03);

+

+// Part (b)

+mew_o = 4*(%pi)*1e-07;                        // Pemeability for free space

+// Since B = mew_o * mew_r * H, solving for mew_r

+mew_r = B/(mew_o*H);                         //Relative permeability of toroid

+printf("\nThe relative permeability of toroid = %4d",mew_r);

+

+// Result

+// Coil current = 55.3 mA

+// The relative permeability of toroid = 5109

diff --git a/1544/CH4/EX4.11/Ch04Ex11.sce b/1544/CH4/EX4.11/Ch04Ex11.sce
new file mode 100755
index 000000000..e9c833f28
--- /dev/null
+++ b/1544/CH4/EX4.11/Ch04Ex11.sce
@@ -0,0 +1,17 @@
+// Scilab code Ex4.11: Pg 125-126 (2008)

+clc; clear;

+mew_o = 4*(%pi)*1e-07;          // Pemeability for free space

+l = 0.15;                       // Mean length, m

+N = 2500;                      // Number of turns

+I = 0.3;                       // Electric current, A

+// Since magnetic field strength is defined as the mmf per metre length of the magnetic circuit, therefore, we have

+H = (N*I)/l;                   // Magnetic field strength, At/m

+B = 0.75;                      // Flux density( value taken from graph ), T

+// Since B = ( mew_r * mew_o * H ), solving for mew_r

+mew_r = B/(mew_o * H);               // Relative permeability

+printf("\nThe flux desity of given toroid = %3.2f T ", B);

+printf("\nThe relative permeability of given toroid = %5.1f", mew_r);

+

+// Result

+// The flux desity of given toroid = 0.75 T

+// The relative permeability of given toroid = 119.4

diff --git a/1544/CH4/EX4.12/Ch04Ex12.sce b/1544/CH4/EX4.12/Ch04Ex12.sce
new file mode 100755
index 000000000..da206322b
--- /dev/null
+++ b/1544/CH4/EX4.12/Ch04Ex12.sce
@@ -0,0 +1,25 @@
+// Scilab code Ex4.12: Pg 126-127 (2008)

+clc; clear;

+mew_o = 4*(%pi)*1e-07;              // Permeability for free space

+l = 0.1875;                         // Mean length, m

+A = 8e-05;                          // Cross- sectional area of of coil, metre-square

+N = 750;                            // Number of turns

+phi = 112e-06;                      // Flux, Wb

+l_gap = 0.5e-03;                    // Average length of the magnetic circuit,m

+B = phi/A;                          // Flux density, Wb

+H = 2000;                           // Magnetic field strength( value taken from graph ), At/m

+F_Fe = H*l;                         // The m.m.f in the iron part of the circuit, At

+// Since F = I*N, solving for I

+I = F_Fe/N;                         // Coil current under normal conditions, A

+// Since B = mew_o * H_gap, solving for H_gap

+H_gap = B/mew_o;                    // Magnetic field strength, At/m

+// Since H_gap = F_gap/l_gap, solving for F_gap

+F_gap = H_gap * l_gap;              // The mmf in the air part of the circuit, At

+F = F_Fe + F_gap;                  // Total circuit mmf, At

+I_new = F/N;                       // Current required to maintain the flux at its original value, A

+printf("\nThe coil current required to produce a flux of %3d micro-weber in the toroid = %3.1f A ", phi/1e-06, I);

+printf("\nCurrent required to maintain the flux at its original value = %5.3f A", I_new);

+

+// Result

+// The coil current required to produce a flux of 112 micro-weber in the toroid = 0.5 A 

+// Current required to maintain the flux at its original value = 1.243 A

diff --git a/1544/CH4/EX4.13/Ch04Ex13.sce b/1544/CH4/EX4.13/Ch04Ex13.sce
new file mode 100755
index 000000000..583fbebcd
--- /dev/null
+++ b/1544/CH4/EX4.13/Ch04Ex13.sce
@@ -0,0 +1,21 @@
+// Scilab code Ex4.13: Pg 127-128 (2008)

+clc; clear;

+l_A = 0.25;                  // Mean length of cicuit A, m

+l_B = 0.15;                  // Mean length of cicuit A, m

+A_A = 11.5e-04;              // Cross-sectional area of circuit A, metre-square

+A_B = 12e-04;                // Cross-sectional area of circuit B, metre-square

+phi = 1.5e-03;               // Flux, Wb

+N = 1000;                    // Number of turns

+B_A = phi/A_A;              // Flux density linked with circuit A, T

+B_B = phi/A_B;              // Flux density linked with circuit B, T

+H_A = 1470;                 // Magnetic field strength of cicuit A( value taken from graph ), At/m

+H_B = 845;                  // Magnetic field strength of cicuit B( value taken from graph ), At/m

+// Snce H = F/l, solving for F 

+F_A = H_A * l_A;                // Magnetic field strength of circuit A, At/m

+F_B = H_B * l_B;                // Magnetic field strength of circuit B, At/m

+F = F_A + F_B;                   // Total circuit m.m.f, At/m

+I = F/N;                          // Coil current, A

+printf("\Coil current in the magnetic circuit = %5.3f A", I);

+

+// Result

+//  Coil current in the magnetic circuit = 0.494 A 

diff --git a/1544/CH4/EX4.14/Ch04Ex14.sce b/1544/CH4/EX4.14/Ch04Ex14.sce
new file mode 100755
index 000000000..2408835ab
--- /dev/null
+++ b/1544/CH4/EX4.14/Ch04Ex14.sce
@@ -0,0 +1,19 @@
+// Scilab code Ex4.14: Pg 129-130 (2008)

+clc; clear;

+A = 8e-04;                                // Cross-sectional area, metre-square

+d = 24e-02;                               // Mean diametre of iron ring, m 

+phi = 1.2e-03;                            // Flux, Wb

+mew_r = 1200;                             // Relative permeability

+mew_o = 4*(%pi)*1e-07;                   // Pemeability for free space

+mew_air = 1;                             // Pemeability for air

+l_gap = 3e-03;                           // Mean length, m

+l_Fe = (%pi) * d;                        // Mean length of iron circuit, m

+S_Fe = l_Fe/(mew_r * mew_o *A);          // Reluctance of iron circuit, At/Wb

+S_gap = l_gap/(mew_air * mew_o *A);        // Reluctance of gap, At/Wb

+S = S_Fe + S_gap;                        // Total circuit reluctance, At/Wb

+// Since phi = F/S, solving for F

+F = phi*S;                               // Magnetomotive force, At

+printf("\nThe required mmf = %5.1f At", F);

+

+// Result

+// The required mmf = 4331 At

diff --git a/1544/CH4/EX4.15/Ch04Ex15.sce b/1544/CH4/EX4.15/Ch04Ex15.sce
new file mode 100755
index 000000000..e60426be0
--- /dev/null
+++ b/1544/CH4/EX4.15/Ch04Ex15.sce
@@ -0,0 +1,32 @@
+// Scilab code Ex4.15: Pg 130-131 (2008)

+clc; clear;

+N = 500;                                       // Number of turns on first section's coil

+phi = 2e-03;                                   // Flux produced by first section, Wb

+l_1 = 85e-02;                                  // Length of first section, m

+l_2 = 65e-02;                                  // Length of second section, m

+l_3 = 0.1e-02;                                 // Length of third section, m

+A_1 = 10e-04;                                  // Csa of first section, m^2

+A_2 = 15e-04;                                  // Csa of second section, m^2

+A_3 = 12.5e-04;                                // Csa of second section, m^2

+mew_o = 4*(%pi)*1e-07;                         // Pemeability for free space

+mew_r1 = 600;                                  // Relative permeability of first section

+mew_r2 = 950;                                  // Relative permeability of second section

+mew_r3 = 1;                                    // Relative permeability of third section

+

+// Part (a)

+S_1 = l_1/(mew_r1 * mew_o * A_1);              // Reluctance of first section, At/Wb

+S_2 = l_2/(mew_r2 * mew_o * A_2);              // Reluctance of first section, At/Wb

+S_3 = l_3/(mew_r3 * mew_o * A_3);              // Reluctance of first section, At/Wb

+S = S_1 + S_2 + S_3;                           // Total reluctance of the circuit, At/Wb

+printf("\nTotal reluctance of the circuit = %4.2fe+06 At/Wb", S*1e-06);

+

+// Part (b)

+// Since phi = F/S, solving for F

+F = phi*S;                                    // Magnetomotive force, At

+// Since F = N*I, solving for I

+I = F/N;                                      // Electric current in first section, A

+printf("\nElectric current in first section = %4.2f A", I);

+

+// Result

+// Total reluctance of the circuit = 2.13e+06 At/Wb

+// Electric current in first section = 8.51 A

diff --git a/1544/CH4/EX4.2/Ch04Ex2.sce b/1544/CH4/EX4.2/Ch04Ex2.sce
new file mode 100755
index 000000000..72a96f8c0
--- /dev/null
+++ b/1544/CH4/EX4.2/Ch04Ex2.sce
@@ -0,0 +1,10 @@
+// Scilab code Ex4.2: Pg 116 (2008)

+clc; clear;

+A = 45e-06;                    // Cross sectional area of pole face, metre-square

+B = 0.6;                      // Flux density, T

+// Using formula B = phi/A, solving for phi

+phi = B*A;                   // Flux, Wb

+printf("\nThe flux produced by pole face = %2d micro-wWb", phi/1e-06);

+

+// Result 

+//The flux produced by pole face = 27 micro-Wb 

diff --git a/1544/CH4/EX4.3/Ch04Ex3.sce b/1544/CH4/EX4.3/Ch04Ex3.sce
new file mode 100755
index 000000000..cda33f323
--- /dev/null
+++ b/1544/CH4/EX4.3/Ch04Ex3.sce
@@ -0,0 +1,14 @@
+// Scilab code Ex4.3: Pg 117 (2008)

+clc; clear;

+N = 1500;                      // Number of turns in a coil

+A = 5e-04;                     // Cross- sectional area of of coil, metre-square

+phi = 0.2e-03;                 // Flux, Wb

+I = 0.75;                      // Coil-current, A

+// Since m.m.f is the product of the current and the number of turns, therefore, we have

+F = N*I;                      // Magnetomotive force, At

+B = phi/A;                    // Flux density, T

+printf("\The m.m.f and flux density produced are %4d At and %3.1f T respectively", F, B);

+

+// Result

+// The m.m.f and flux density produced are 1125 At and 0.4 T respectively

+

diff --git a/1544/CH4/EX4.4/Ch04Ex4.sce b/1544/CH4/EX4.4/Ch04Ex4.sce
new file mode 100755
index 000000000..354134ca6
--- /dev/null
+++ b/1544/CH4/EX4.4/Ch04Ex4.sce
@@ -0,0 +1,10 @@
+// Scilab code Ex4.4:Pg 117 (2008)

+clc; clear;

+N = 600;                      // Number of turns in a coil

+F = 1500;                      // Magnetomotive force, At

+// Since magnetomotive force,F = N*I, solving for I

+I = F/N;                      // Excitation-current, A

+printf("\nThe excitation current required = %3.1f A", I);

+

+// Result

+// The excitation current required = 2.5 A

diff --git a/1544/CH4/EX4.5/Ch04Ex5.sce b/1544/CH4/EX4.5/Ch04Ex5.sce
new file mode 100755
index 000000000..839d55d14
--- /dev/null
+++ b/1544/CH4/EX4.5/Ch04Ex5.sce
@@ -0,0 +1,12 @@
+// Scilab code Ex4.5: Pg 118 (2008)

+clc; clear;

+I = 0.4;                  // Current, A

+N = 550;                  // Number of turns in a coil

+d = 8e-02;                // Diameter, m

+l = (%pi*d);              // Average length of the magnetic circuit, m

+// Since magnetic field strength is defined as the mmf per metre length of the magnetic circuit, therefore, we have

+H = (N*I)/l;                // Magnetic field strength, At/m

+printf("\nThe magnetic field strength inside the toroid = %6.2f At/m", H);

+

+// Result

+// The magnetic field strength inside the toroid = 875.35 At/m 

diff --git a/1544/CH4/EX4.6/Ch04Ex6.sce b/1544/CH4/EX4.6/Ch04Ex6.sce
new file mode 100755
index 000000000..b4301c306
--- /dev/null
+++ b/1544/CH4/EX4.6/Ch04Ex6.sce
@@ -0,0 +1,16 @@
+// Scilab code Ex4.6: Pg 119-120 (2008)

+clc; clear;

+A = 15e-04;                   // Cross-sectional area of core, metre-square

+mew_r1 = 65;                  // Relative permeability of core

+phi_1 = 2e-04;                // Flux, Wb

+mew_r2 = 800;                 // Changed relative permeability of core

+B_1 = phi_1/A;                // Flux density, T

+mew_r = mew_r2/mew_r1;        // Relative permeability of core

+// Since cross-sectional area of core A remains constant, therefore, we have mew_r = B_1/B_2 , solving for B_2

+B_2 = mew_r*B_1;              //  New flux density, T

+// Since B_2 = phi_2/A, solving for phi_2

+phi_2 = B_2*A;                // New flux, Wb

+printf("\nThe new flux and flux density are %5.3f mWb and %5.3f T respectively", phi_2/1e-03, B_2);

+

+// Result

+// The new flux and flux density are 2.462 mWb and 1.641 T respectively

diff --git a/1544/CH4/EX4.7/Ch04Ex7.sce b/1544/CH4/EX4.7/Ch04Ex7.sce
new file mode 100755
index 000000000..0a401b831
--- /dev/null
+++ b/1544/CH4/EX4.7/Ch04Ex7.sce
@@ -0,0 +1,29 @@
+// Scilab code Ex4.7: Pg 120 (2008)

+clc; clear;

+r = 0.04;                           // Mean radius of torod, m

+A = 3e-04;                          // Csa of toroid, m^2

+mew_o = 4*(%pi)*1e-07;              // Permeability of free space

+mew_r = 150;                        // Relative permeability of toroid

+N = 900;                            // Number of turns on coil

+I = 1.5;                            // Coil current, A

+l = 2*(%pi)*r;                      // Effective length of toroid, m

+

+// Part (a)

+// Since m.m.f is the product of the current and the number of turns, therefore, we have

+F = N*I;                            // Magnetomotive force, At

+printf("\nThe m.m.f of toroid = %4d At", F);

+

+// Part (b)

+// Since magnetic field strength is defined as the mmf per metre length of the magnetic circuit, therefore, we have

+H = F/l;                               // Magnetic field strength, At/m

+printf("\nThe magntic field strength = %6.1f At/m", H);

+

+// Part (c)

+B = (mew_r*mew_o*H);                   // Flux density, T

+phi = B*A;                             // Flux, Wb

+printf("\nThe flux and flux density are %6.2f micro-weber and %6.4f T respectively", phi/1e-06, B)

+

+// Result

+// The m.m.f of toroid = 1350 At

+// The magntic field strength = 5371.5 At/m

+// The flux and flux density are 303.75 micro-weber and 1.0125 T respectively

diff --git a/1544/CH4/EX4.8/Ch04Ex8.sce b/1544/CH4/EX4.8/Ch04Ex8.sce
new file mode 100755
index 000000000..af2a4f265
--- /dev/null
+++ b/1544/CH4/EX4.8/Ch04Ex8.sce
@@ -0,0 +1,20 @@
+// Scilab code Ex4.8: Pg 120-121 (2008)

+clc; clear;

+r = 3e-02;                            // Radius of toroid, m

+A = 4.5e-04;                          // Cross-sectional area of toroid, metre-square

+N = 500;                              // Number of turns 

+phi = 250e-06;                        // Flux, Wb

+mew_o = 4*(%pi)*(1e-07);              // Permeability of free space

+mew_r = 300;                          // Relative permeability

+l = 2*(%pi)*r;                        // Effective length, m

+B = phi/A;                            // Flux density, T

+// Since B = (mew_r)*(mew_o)*H, solving for H

+H = B /((mew_r)*(mew_o));              // Magnetic field strength, At/m

+// Since H = F/l, solving for F

+F = H*l;                               // Magnetomotive force, At

+// Since mmf,F = N*I, solving for I

+I = F/N;                               // Electric current, A

+printf("\nThe value of current needs to be passed through the coil = %4.2f A", I);

+

+// Result

+// The value of current needs to be passed through the coil = 0.56 A

diff --git a/1544/CH4/EX4.9/Ch04Ex9.sce b/1544/CH4/EX4.9/Ch04Ex9.sce
new file mode 100755
index 000000000..1e65c2533
--- /dev/null
+++ b/1544/CH4/EX4.9/Ch04Ex9.sce
@@ -0,0 +1,26 @@
+// Scilab code Ex4.9: Pg 121-122 (2008)

+clc; clear; 

+// Part (a)

+I = 0.2;                    // Electric current, A

+l = 5e-02;                  // Effective length, m

+A = 7e-04;                  // Cross-sectional area, metre-square

+d = 0.5e-03;                // Diametre, m

+mew_r = 1;                  //Relative permeability for wood

+mew_o = 4*(%pi)*1e-07;      // Pemeability for free space

+N = l/d;                    // Number of turns 

+// Since mmf is the product of the current and the number of turns, therefore, we have

+F = N*I;                     // Magnetomotive force, At

+// Part (b)

+// Since magnetic field strength is defined as the mmf per metre length of the magnetic circuit, therefore, we have

+H = F/l;                               // Magnetic field strength, At/m

+B = ( mew_r * mew_o * H );             //  Flux density, T

+// Part (c)

+phi = B * A;                           // Flux, Wb

+printf("\nThe mmf produced = %2d At", F);

+printf("\nThe flux density produced = %3d micro-tesla", B/1e-06);

+printf("\nThe flux produced = %5.3f micro-weber", phi/1e-06);

+

+// Result

+// The mmf produced = 20 At

+// The flux density produced = 502 micro-tesla

+// The flux produced = 0.352 micro-weber