1 files changed, 56 insertions, 0 deletions
diff --git a/1319/CH1/EX1.41/1_41.sce b/1319/CH1/EX1.41/1_41.sce
new file mode 100644
index 000000000..95aa107a5
--- /dev/null
+++ b/1319/CH1/EX1.41/1_41.sce
@@ -0,0 +1,56 @@
+// To determine the parameters of a toroid
+
+clc;
+clear;
+
+//Arc length of different materials
+lni=0.3; // nickel iron alloy
+lss=0.2;//Silicon Steel
+lcs=0.1;//Cast steel
+
+A=1*(10^-3); // Area of cross section of Toroid
+
+N=100; // Hundred Turns
+
+phi=6*(10^-4); // Flux
+B=0.6; // Flux densities
+m0=4*%pi*(10^-7);
+
+//Field intensities
+Hni=10; // Ni alloy
+Hss=77;// Si Steel
+Hcs=270; // CAst steel
+
+mmf=(Hni*lni)+(Hss*lss)+(Hcs*lcs); // Total mmf
+
+B=phi/A;
+
+I=mmf/N;
+
+//P = Permeability, RP = Relative permeablity
+
+//For Nickel Alloy
+mni=B/Hni; // P
+mrni=mni/m0; // RP
+
+//For Si Steel
+mss=B/Hss; //P
+mrss=mss/m0; // RP
+
+//For Cast Steel
+mcs=B/Hcs; //P
+mrcs=mcs/m0; // RP
+
+deff('x=rel(y,z)','x=(y*z)/phi'); // Fucntion to find out Reluctance
+
+//Reluctances
+Rni=rel(Hni,lni);
+Rss=rel(Hss,lss);
+Rcs=rel(Hcs,lcs);
+
+printf('i) The mmf required to establish a magnetic flux of 6 *(10^-4) Wb = %g AT\n',mmf)
+printf('ii) The Current through the coil = %g A\n',I)
+printf('iii) Relative permeability and reluctance of each ferro magnetic material\n')
+printf('    Nickel Iron alloy is %g and %g ohms \n',mrni,Rni)
+printf('    Silicon Steel is %g and %g ohms \n',mrss,Rss)
+printf('    Cast Steel is %g and %g ohms \n',mrcs,Rcs)