From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001
From: priyanka
Date: Wed, 24 Jun 2015 15:03:17 +0530
Subject: initial commit / add all books

---
 1646/CH12/EX12.1/Ch12Ex1.sce   | 12 ++++++++++++
 1646/CH12/EX12.10/Ch12Ex10.sce | 11 +++++++++++
 1646/CH12/EX12.2/Ch12Ex2.sce   | 17 +++++++++++++++++
 1646/CH12/EX12.3/Ch12Ex3.sce   | 20 ++++++++++++++++++++
 1646/CH12/EX12.4/Ch12Ex4.sce   | 16 ++++++++++++++++
 1646/CH12/EX12.5/Ch12Ex5.sce   | 19 +++++++++++++++++++
 1646/CH12/EX12.6/Ch12Ex6.sce   | 13 +++++++++++++
 1646/CH12/EX12.7/Ch12Ex7.sce   | 26 ++++++++++++++++++++++++++
 1646/CH12/EX12.8/Ch12Ex8.sce   | 13 +++++++++++++
 1646/CH12/EX12.9/Ch12Ex9.sce   | 11 +++++++++++
 10 files changed, 158 insertions(+)
 create mode 100755 1646/CH12/EX12.1/Ch12Ex1.sce
 create mode 100755 1646/CH12/EX12.10/Ch12Ex10.sce
 create mode 100755 1646/CH12/EX12.2/Ch12Ex2.sce
 create mode 100755 1646/CH12/EX12.3/Ch12Ex3.sce
 create mode 100755 1646/CH12/EX12.4/Ch12Ex4.sce
 create mode 100755 1646/CH12/EX12.5/Ch12Ex5.sce
 create mode 100755 1646/CH12/EX12.6/Ch12Ex6.sce
 create mode 100755 1646/CH12/EX12.7/Ch12Ex7.sce
 create mode 100755 1646/CH12/EX12.8/Ch12Ex8.sce
 create mode 100755 1646/CH12/EX12.9/Ch12Ex9.sce

(limited to '1646/CH12')

diff --git a/1646/CH12/EX12.1/Ch12Ex1.sce b/1646/CH12/EX12.1/Ch12Ex1.sce
new file mode 100755
index 000000000..65b5e9fed
--- /dev/null
+++ b/1646/CH12/EX12.1/Ch12Ex1.sce
@@ -0,0 +1,12 @@
+// Scilab Code Ex12.1 Page-603 (2011)
+clc; clear;
+H = 5e+3;....// Coercivity of a bar magnet, A/m
+L = 0.1;....// Length of the solenoid, m
+N = 50;....// Turns in solenoid
+n = 500;....// Turns/m
+// Using the relation
+I = H/n;....// where I is the current through the solenoid
+printf("\nThe current through the solenoid is = %2d A", I);
+
+// Result
+// The current through the solenoid is = 10 A
diff --git a/1646/CH12/EX12.10/Ch12Ex10.sce b/1646/CH12/EX12.10/Ch12Ex10.sce
new file mode 100755
index 000000000..0c36bfa60
--- /dev/null
+++ b/1646/CH12/EX12.10/Ch12Ex10.sce
@@ -0,0 +1,11 @@
+// Scilab Code Ex12.10 : Page-608 (2011)
+clc; clear;
+e = 1.6e-19;....// Charge on anlectron, C
+m = 9.1e-31;....// Mass of the electron, kg
+r = 5.1e-11;....// Radius of the electronic orbit, m
+B = 2.0;....// Applied magnetic field, weber per metre-square
+delta_pm = e^2*r^2*B/(4*m);
+printf("\nThe change in the magnetic dipole moment of the electron = %3.1e A-metre square", delta_pm);
+
+// Result
+// The change in the magnetic dipole moment of the electron = 3.7e-29 A-metre square 
diff --git a/1646/CH12/EX12.2/Ch12Ex2.sce b/1646/CH12/EX12.2/Ch12Ex2.sce
new file mode 100755
index 000000000..da83f3011
--- /dev/null
+++ b/1646/CH12/EX12.2/Ch12Ex2.sce
@@ -0,0 +1,17 @@
+// Scilab Code Ex12.2 : Page-603 (2011)
+clc; clear;
+n = 500;....// Number of turns wound per metre on the solenoid
+i = 0.5;....// Current through the solenoid, A
+V = 1e-03;....// Volume of iron rod, per metre cube
+mu_r = 1200;    // Relative permeability of the iron
+H = n*i;    // Magnetic intensity inside solenoid, ampere-turn per metre
+// As B = mu_o * (H + I) => I = B/mu_o -  H 
+// But B = mu_o * mu_r * H and solving for I 
+I = (mu_r - 1) * H;
+printf("\nThe Intensity of magnetisation inside the solenoid, I = %5.3e A/m", I);
+M = I * V;    // Magnetic moment of the rod, ampere metre square 
+printf("\nThe magnetic moment of the rod, M = %3d ampere metre square", M)
+
+//Result
+// The Intensity of magnetisation inside the solenoid, I = 2.998e+005 A/m
+// The magnetic moment of the rod, M = 299 ampere metre square 
diff --git a/1646/CH12/EX12.3/Ch12Ex3.sce b/1646/CH12/EX12.3/Ch12Ex3.sce
new file mode 100755
index 000000000..e3c39005b
--- /dev/null
+++ b/1646/CH12/EX12.3/Ch12Ex3.sce
@@ -0,0 +1,20 @@
+// Scilab Code Ex12.3 : Page-604 (2011)
+clc; clear;
+n = 300;....// Number of turns wound per metre on the solenoid
+i = 0.5;....// Current through the solenoid, A
+V = 1e-03;....// Volume of iron rod, per metre cube
+mu_r = 100;    // Relative permeability of the iron
+H = n*i;    // Magnetic intensity inside solenoid, ampere-turn per metre
+// As, I = (B-mu_o* H)/mu_o
+//But, B= mu * H = mu_r * mu_o * H and I = (mu_r-1)* H
+I = (mu_r-1)*n*i; 
+printf("\nThe Intensity of magnetisation inside the solenoid, I = %5.3e A/m", I);
+l = 0.2;....//length of the rod,m
+r = 5e-3;....//radius of the rod,m
+V = 1.57e-5;....//V=%pi*r^2*l where the volume of the rod having radius r and length,m
+ M = I * V ;    // Magnetic moment of the rod, ampere metre square
+printf("\nThe magnetic moment of the rod, M = %5.3f ampere metre square",M)
+
+//Result
+// The Intensity of magnetisation inside the solenoid, I = 1.485e+004 A/m
+// The magnetic moment of the rod, M = 0.233 ampere metre square 
diff --git a/1646/CH12/EX12.4/Ch12Ex4.sce b/1646/CH12/EX12.4/Ch12Ex4.sce
new file mode 100755
index 000000000..12db96406
--- /dev/null
+++ b/1646/CH12/EX12.4/Ch12Ex4.sce
@@ -0,0 +1,16 @@
+// Scilab Code Ex12.4 : Page-605 (2011)
+clc; clear;
+B = 0.0044;....// Magnetic flux density, weber/meter square
+mu_o = 4*%pi*1e-07;....// Relative permeability of the material, henery/m 
+I = 3300;....// Magnetization of a magnetic material, A/m
+//B = mu_o*(I+H), solving for H  
+H = (B/mu_o)- I;....// Magnetizing force ,A/m
+printf("\nThe magnetic intensity,H = %3d A/m",H);
+// Relation between intensity of magnetization and relative permeability 
+mu_r = (I/H)+1;....//substitute the value of I and H 
+printf("\nThe relative permeability, mu_r = %5.2f",mu_r);
+
+//Result
+// The magnetic intensity,H = 201 A/m
+// The relative permeability, mu_r = 17.38 
+
diff --git a/1646/CH12/EX12.5/Ch12Ex5.sce b/1646/CH12/EX12.5/Ch12Ex5.sce
new file mode 100755
index 000000000..5b3acb1c6
--- /dev/null
+++ b/1646/CH12/EX12.5/Ch12Ex5.sce
@@ -0,0 +1,19 @@
+// Scilab Code Ex12.5 : Page-605 (2011)
+clc; clear;
+mu_o = 4*%pi*1e-07;....// Magnetic permeability of the free space, henery/m
+mu_r = 600;
+mu = mu_o*mu_r;     // Magnetic permeability of the medium, henery/m
+n = 500;...// Turns in a wire
+i = 0.3;....// Current flows through a ring,amp
+r = 12e-02/2;....// Mean radius of a ring, m
+B = mu_o*mu_r*n*i/(2*%pi*r);
+printf("\nThe magnetic flux density = %2.1f weber/meter-square", B);
+H = B/mu;   // Magnetic intensity, ampere-turns/m
+printf("\nThe magnetic intensity = %5.1f ampere-turns/m", H);
+// As B = mu_o*(I + H) => mu_o*I = B - mu_o*H
+printf("\nThe percentage magnetic flux density due to electronic loop currents = %5.2f percent", (B - mu_o*H)/B*100);
+
+// Result
+// The magnetic flux density = 0.3 weber/meter-square
+// The magnetic intensity = 397.9 ampere-turns/m
+// The percentage magnetic flux density due to electronic loop currents = 99.83 percent 
diff --git a/1646/CH12/EX12.6/Ch12Ex6.sce b/1646/CH12/EX12.6/Ch12Ex6.sce
new file mode 100755
index 000000000..2cdffd680
--- /dev/null
+++ b/1646/CH12/EX12.6/Ch12Ex6.sce
@@ -0,0 +1,13 @@
+// Scilab Code Ex12.6 : Page-606 (2011)
+clc; clear;
+M_i = 4.5;....// Intial value of total dipole moment of the sample
+H_i = 0.84;....// External magnetic field, tesla
+T_i = 4.2;....// Cooling temerature of the sample, K
+H_f = 0.98;....// External magnetic field, tesla
+T_f = 2.8;....// Cooling temerature of the sample, K
+// According to the curie's law, Mf/Mi = (Hf/Hi)*(Ti/Tf)
+M_f = M_i*H_f/H_i*T_i/T_f;
+printf("\nThe total dipole moment of the sample = %5.3f joule/tesla",M_f);
+
+// Result
+// The total dipole moment of the sample = 7.875 joule/tesla 
diff --git a/1646/CH12/EX12.7/Ch12Ex7.sce b/1646/CH12/EX12.7/Ch12Ex7.sce
new file mode 100755
index 000000000..653aec2f4
--- /dev/null
+++ b/1646/CH12/EX12.7/Ch12Ex7.sce
@@ -0,0 +1,26 @@
+// Scilab Code Ex12.7 : Page-606 (2011)
+clc; clear;
+mu_o = 4*%pi*1e-07;....// Magnetic permeability of free space, henry/m
+n = 1e+29;....// Number density of atoms of iron, per metre cube 
+p_m = 1.8e-23;....// Magnetic moment of each atom, ampere-metre square  
+k_B = 1.38e-23;....// Boltzmann constant, J/K
+B = 0.1;    // Magnetic flux density, weber per metre square
+T = 300;....// Absolute room temperature, K
+l = 10e-02; // Length of the iron bar, m
+a = 1e-04;  // Area of cross-section of the iron bar, metre square
+V = l*a;    // Voluem of the iron bar, metre cube
+chi = n*p_m^2*mu_o/(3*k_B*T);
+printf("\nThe paramagnetic susceptibility of a material = %5.3e", chi);
+pm_mean = p_m^2*B/(3*k_B*T); // Mean dipole moment of an iron atom, ampere metre-square
+P_m = n*pm_mean;  // Dipole moment of the bar, ampere metre-square
+I = n*p_m;  // Magnetization of the bar in one domain, ampere/metre
+M = I*V;    // Magnetic moment of the bar, ampere metre-square
+printf("\nThe dipole moment of the bar = %5.3e ampere metre-square", P_m);
+printf("\nThe magnetization of the bar in one domain = %3.1e ampere/metre", I);
+printf("\nThe magnetic moment of the bar = %2d ampere metre-square", M);
+
+// Result
+// The paramagnetic susceptibility of a material = 3.278e-03
+// The dipole moment of the bar = 2.609e+02 ampere metre-square
+// The magnetization of the bar in one domain = 1.8e+06 ampere/metre
+// The magnetic moment of the bar = 18 ampere metre-square 
diff --git a/1646/CH12/EX12.8/Ch12Ex8.sce b/1646/CH12/EX12.8/Ch12Ex8.sce
new file mode 100755
index 000000000..c2ecbf37b
--- /dev/null
+++ b/1646/CH12/EX12.8/Ch12Ex8.sce
@@ -0,0 +1,13 @@
+// Scilab Code Ex12.8 :Page-607 (2011)
+clc; clear;
+A = 500;....// Area of the B-H loop, joule per metre cube
+n = 50;....// Total number of cycles, Hz
+m = 9;....// Mass of the core, kg
+d = 7.5e+3;....// Density of the core, kg/metre cube
+t = 3600;....// Time during which the energy loss takes place, s
+V = m/d;....// Volume of the core, metre cube
+E = n*V*A*t;....// Hystersis loss of energy per hour, joule
+printf("\nThe hystersis loss per hour = %5.2eJ", E);
+
+// Result
+// The hystersis loss per hour = 1.08e+005J 
diff --git a/1646/CH12/EX12.9/Ch12Ex9.sce b/1646/CH12/EX12.9/Ch12Ex9.sce
new file mode 100755
index 000000000..644613b67
--- /dev/null
+++ b/1646/CH12/EX12.9/Ch12Ex9.sce
@@ -0,0 +1,11 @@
+// Scilab Code Ex12.9 : Page-607 (2011)
+clc; clear;
+n = 50;....// Total number of cycles per sec, Hz
+V = 1e-03;....// Volume of the specimen, metre cube
+t = 1;....// Time during which the loss occurs, s
+A = 0.25e+03;....// Area of B-H loop, joule per metre cube
+E = n*V*A*t;    // Energy loss due to hysteresis, J/s
+printf("\nThe hystersis loss = %4.1f J/s", E); 
+
+// Result
+// The hystersis loss = 12.5 J/s 
-- 
cgit