diff options
Diffstat (limited to '3871/CH13')
| -rw-r--r-- | 3871/CH13/EX13.1/Ex13_1.sce | 21 | ||||
| -rw-r--r-- | 3871/CH13/EX13.10/Ex13_10.sce | 30 | ||||
| -rw-r--r-- | 3871/CH13/EX13.11/Ex13_11.sce | 19 | ||||
| -rw-r--r-- | 3871/CH13/EX13.12/Ex13_12.sce | 22 | ||||
| -rw-r--r-- | 3871/CH13/EX13.13/Ex13_13.sce | 26 | ||||
| -rw-r--r-- | 3871/CH13/EX13.14/Ex13_14.sce | 28 | ||||
| -rw-r--r-- | 3871/CH13/EX13.15/Ex13_15.sce | 17 | ||||
| -rw-r--r-- | 3871/CH13/EX13.16/Ex13_16.sce | 23 | ||||
| -rw-r--r-- | 3871/CH13/EX13.17/Ex13_17.sce | 26 | ||||
| -rw-r--r-- | 3871/CH13/EX13.2/Ex13_2.sce | 25 | ||||
| -rw-r--r-- | 3871/CH13/EX13.4/Ex13_4.sce | 20 | ||||
| -rw-r--r-- | 3871/CH13/EX13.5/Ex13_5.sce | 22 | ||||
| -rw-r--r-- | 3871/CH13/EX13.6/Ex13_6.sce | 15 | ||||
| -rw-r--r-- | 3871/CH13/EX13.7/Ex13_7.sce | 24 | ||||
| -rw-r--r-- | 3871/CH13/EX13.8/Ex13_8.sce | 33 | ||||
| -rw-r--r-- | 3871/CH13/EX13.9/Ex13_9.sce | 25 |
16 files changed, 376 insertions, 0 deletions
diff --git a/3871/CH13/EX13.1/Ex13_1.sce b/3871/CH13/EX13.1/Ex13_1.sce new file mode 100644 index 000000000..854c1566e --- /dev/null +++ b/3871/CH13/EX13.1/Ex13_1.sce @@ -0,0 +1,21 @@ +//=====================================================================================
+//Chapter 13 example 1
+
+clc;clear all;
+
+//variable declaration
+l = 0.6; //length of solenoid in m
+N = 600; //number of turns
+I = 2; //current passing through solenoid in A
+ur = 1; //air coiled solenoid
+r = 0.025; //radius in m
+
+//calculations
+H = (N*I)/(l); //magnetic field at the centre in AT/metre
+u0 = 4*(%pi)*(10^-7); //flux
+a = ((%pi)/(4))*(r^2); //area
+phi = ur*u0*H*a; //flux passing through thesecondary coil
+
+//calculations
+mprintf("magnetic field = %3.2f AT/metre",H);
+mprintf("flux = %3.2e Wb",phi);
diff --git a/3871/CH13/EX13.10/Ex13_10.sce b/3871/CH13/EX13.10/Ex13_10.sce new file mode 100644 index 000000000..21bae3cdd --- /dev/null +++ b/3871/CH13/EX13.10/Ex13_10.sce @@ -0,0 +1,30 @@ +//=====================================================================================
+//Chapter 13 example 10
+
+clc;clear all;
+
+//variable declaration
+d = 0.3; //diameter in m
+a = 4*10^-4; //cross sectional area of iron ring in m**2
+N = 80; //number of turns on magnetising coil
+Ns = 30; //number of turns on secondary coil
+F = 0.1*10^-3; //flux meter constant in Wb-turn
+D = 46; //deflection factor
+I = 2; //current in A
+
+//calculations
+//phi = (N*I*u0*ur*a)/l
+// phi = x*l
+//lat x = (N*I*u0*a)/l
+l = d*%pi;
+u0 = 4*(%pi)*10^-7;
+x =(N*I*u0*a)/(l);
+// total change in Wb-turns y = 2*phi*Ns = 2*x*ur*Ns
+y = 2*x*Ns;
+df = F*D; //change in flux measuredby the flux meter in wb-turns
+ur =df/y; //relative permeability
+
+//result
+mprintf("relative permeabitlity = %3.0d",ur);
+mprintf("\n Note:textbook answer represents the approximate value")
+
diff --git a/3871/CH13/EX13.11/Ex13_11.sce b/3871/CH13/EX13.11/Ex13_11.sce new file mode 100644 index 000000000..65e751633 --- /dev/null +++ b/3871/CH13/EX13.11/Ex13_11.sce @@ -0,0 +1,19 @@ +//=====================================================================================
+//Chapter 13 example 13
+
+clc;clear all;
+
+//variable declaration
+Q = 1000; //Charge passed through galvanometer in uC
+theta1 = 60;
+d = 10; //defelction in mm
+r = 1000; //m=circular scale
+
+//calculations
+theta2 = %pi/(3) //conversion of degrees to radians
+K = Q/(theta2); //galvanometer constant in uC/radian
+theta = d/(2*r); //angle turned through by reflected ray for aswing of 10 mm
+Q1 = K*theta; //charge for a swing of 0.005 radian in uC
+
+//result
+mprintf("charge for a swing of 0.005 radian = %3.2f uC",Q1);
diff --git a/3871/CH13/EX13.12/Ex13_12.sce b/3871/CH13/EX13.12/Ex13_12.sce new file mode 100644 index 000000000..8c5da438c --- /dev/null +++ b/3871/CH13/EX13.12/Ex13_12.sce @@ -0,0 +1,22 @@ +//=====================================================================================
+//Chapter 13 example 12
+
+clc;clear all;
+
+//variable declaration
+f = 50; //number of reversals
+m = 1; //mass
+d = 7.8*10**3; //density
+A = 4800; //area of the loop m^3
+x = 200; //in AT/m
+x1 = 10; // 1 unit in mm
+y1 = 10; // 1 unit in mm
+y = 0.1; //in T
+
+//claculations
+V = m/d; //volume of magnetic material in m^3
+l = A*(x/x1)*(y/y1);
+l1 = l*V*f; //hysteresis loss in watts per kg at 50 Hz
+
+//result
+mprintf("hysteresis loss at 50 Hz = %3.3f watts per kg ",l1);
diff --git a/3871/CH13/EX13.13/Ex13_13.sce b/3871/CH13/EX13.13/Ex13_13.sce new file mode 100644 index 000000000..24a319699 --- /dev/null +++ b/3871/CH13/EX13.13/Ex13_13.sce @@ -0,0 +1,26 @@ +clc;
+clear all;
+
+//variable declaration
+f = 60; //supply frequency in Hz
+Pi = 360; //iron loss in W
+f = 60;
+//Pe =6*Ph;
+//Pi = Pe+Ph
+//360= (6*Ph)+Ph
+Ph = Pi/7; //hysteresis loss in W
+Pe = Pi-Ph; //eddy current loss in W
+ //Ph1 = (f1/f)*Ph
+Ph1 = (1/f)*Pe; //hysteresis loss in watts
+//Ph1 =Ph1*f1
+//Pe1 = ((f1/f)^2)*Pe
+Pe1 =((1/f)^2)*Pe; //eddy current loss
+//Pe1 = Pe1*Pe
+Pi1 =Ph1+Pe1;
+Pi1 = 2*Pi;
+//720 = 0.857*f1+(0.0857*f1^2)
+f1 =86.8
+
+
+//result
+mprintf("new supplyfrequency = %3.2f HZ",f1);
diff --git a/3871/CH13/EX13.14/Ex13_14.sce b/3871/CH13/EX13.14/Ex13_14.sce new file mode 100644 index 000000000..71df4d8d5 --- /dev/null +++ b/3871/CH13/EX13.14/Ex13_14.sce @@ -0,0 +1,28 @@ +//=====================================================================================
+//Chapter 13 example 14
+
+clc;
+clear all;
+
+//variable declaration
+//Ph = A*f
+//Pe = B*f^2
+//Pi = Ph+Pe
+Pi = 17.2; //power in W
+f = 50; //frequency in Hz
+Pi1 = 28.9; //iron loss in W
+
+m = 13; //weight in kg
+
+//calculations
+//17.2 = 40*A+((40)^2)*B
+//28.9 = 60*A+((60)^2)*B
+A = 0.326667
+B = 0.002583
+Ph = (A*f)/m //hysteresis loss per kg in W
+Pe = (B*(f^2))/m //eddy current loss per kg in W
+
+//result
+mprintf("hysteresis loss per kg = %3.2f W",Ph);
+mprintf("\neddy current loss per kg = %3.3f W",Pe);
+
diff --git a/3871/CH13/EX13.15/Ex13_15.sce b/3871/CH13/EX13.15/Ex13_15.sce new file mode 100644 index 000000000..a8dd818cd --- /dev/null +++ b/3871/CH13/EX13.15/Ex13_15.sce @@ -0,0 +1,17 @@ +//=====================================================================================
+//Chapter 13 example 15
+clc;
+clear all;
+
+//variable declaration
+A = 0.5;
+B = 0.01;
+f = 50;
+n = 10;
+
+//calculations
+Pe = B*(f^2); //eddy current loss at 50 Hz in W
+Pe1 = Pe/n; //eddy current loss per kg at 50 Hz in watts
+
+//result
+mprintf("eddy current loss per kg at 50 Hz = %3.2f watts",Pe1);
diff --git a/3871/CH13/EX13.16/Ex13_16.sce b/3871/CH13/EX13.16/Ex13_16.sce new file mode 100644 index 000000000..eb6052971 --- /dev/null +++ b/3871/CH13/EX13.16/Ex13_16.sce @@ -0,0 +1,23 @@ +//=====================================================================================
+//Chapter 13 example 16
+
+clc;
+clear all;
+
+//variable declaration
+x = 0.8; //Kf2/Kf1
+y =1.2;
+
+//Pe2/Pe1 = (Kf2/Kf1)^2
+p = x^2;
+//Pe2 = p*Pe1; //
+//p1 = (Pe1-Pe2)/Pe1;
+p1 = (1-p)*100; //percentage change in hysteresis current loss
+p2 = y^2;
+p12 = (y-1)*100; //percentage change in hysteresis current loss
+p3 =(p2-1)*100; //percentage change in eddy current loss in %
+
+//result
+mprintf("percentage change in hysteresis current loss = %3.3f percentage decrease",p1);
+mprintf("\npercentage change in hysteresis current loss = %3.3f percentage increase",p12);
+mprintf("\npercentage change in eddy current loss in = %3.2f percentage increase",p3);
diff --git a/3871/CH13/EX13.17/Ex13_17.sce b/3871/CH13/EX13.17/Ex13_17.sce new file mode 100644 index 000000000..693a5f32f --- /dev/null +++ b/3871/CH13/EX13.17/Ex13_17.sce @@ -0,0 +1,26 @@ +//=====================================================================================
+//Chapter 13 example 17
+
+clc;
+clear all;
+
+//variable declaration
+w = 0.03; //width of plates in m
+n = 51; //number of plates
+t = 0.000489; //thickness in m^3
+f = 50; //frequency in Hz
+Bmax = 1;
+N = 600;
+P1 = 3; //copper loss in watts
+m = 11; //weight in kg
+
+//calculations
+A = w*n*t; //mean area of plates in m^3
+E = 4.44*f*Bmax*A*N; //induced voltage in V
+//from graph corresponding to voltage of 100 volts
+P2 = 30.5; //total losses in watts
+P = P2-P1; //iron loss in watts
+PL = P/m; //loss per kg in watts
+
+//result
+mprintf("iron loss per kg = %3.2f watts",PL);
diff --git a/3871/CH13/EX13.2/Ex13_2.sce b/3871/CH13/EX13.2/Ex13_2.sce new file mode 100644 index 000000000..5238cb299 --- /dev/null +++ b/3871/CH13/EX13.2/Ex13_2.sce @@ -0,0 +1,25 @@ +//=====================================================================================
+//Chapter 13 example 2
+clc;clear all;
+
+//variable declaration
+m = 100; //number of turns
+n = 1000; //turns per m
+theta1 = 10; //first throw in mm
+theta2 = 9.5; //second throw in mm
+I =10; //current in A
+R = 500; // resistance in Ω
+A = 0.002; //area in m**2
+
+//calculations
+//Q = (8*(math.pi)*(10**-7)*N*Ns*I*A)/(l*R) //in columbs
+//Q = (8*(math.pi)*(10**-7)*n*l*m*I*A)/(l*R)
+//Q =(8*(math.pi)*(10**-7)*n*m*I*A)/(R)
+lamda = log(theta1/(theta2));
+theta = theta1*(1+(lamda/(2)));
+K =(8*(%pi)*(10^-7)*n*m*I*A)/(R*theta); //galvanometer constant in C/mm
+
+//result
+
+mprintf("galvanometer constant = %3.2e C/mm",K);
+
diff --git a/3871/CH13/EX13.4/Ex13_4.sce b/3871/CH13/EX13.4/Ex13_4.sce new file mode 100644 index 000000000..8ec1e4a27 --- /dev/null +++ b/3871/CH13/EX13.4/Ex13_4.sce @@ -0,0 +1,20 @@ +//=====================================================================================
+//Chapter 13 example 4
+
+
+clc;clear all;
+
+//variable declaration
+T0 = 4; //time of swing in seconds
+Ig =0.001; //current in A
+lamda = 0;
+theta = 50; //steady deflection in scale divisions
+theta1 = 220; //maximum throw in scale division
+V =100; //potential of the condenser in V
+
+//calculations
+Q = (T0/(2*%pi))*(Ig/theta)*(1+(lamda/2))*theta1; //quantity of electricity discharged in uC
+C = Q/(V); //capacity of the condenser in F
+
+//result
+mprintf("capacity of the condenser = %3.2d uF",(C*10^6));
diff --git a/3871/CH13/EX13.5/Ex13_5.sce b/3871/CH13/EX13.5/Ex13_5.sce new file mode 100644 index 000000000..2c0bf2960 --- /dev/null +++ b/3871/CH13/EX13.5/Ex13_5.sce @@ -0,0 +1,22 @@ +//=====================================================================================
+//Chapter 13 example 5
+clc;clear all;
+
+//variable declaration
+N = 1; //number of turns on search coil
+Rc = 0.025; //resistance of search coil in Ω
+Nw = 1.5*10^-4; //number of wb-turns required for deflection of 1 division
+M = 120000; //reluctane of magnetic circuit
+MMF = 8000; //magnetic circuit is excited in ampere-turn
+f = 1.5*10^-4; //fluxmeter without shunt (K/N = phi/theta)
+theta =120;
+
+//calculations
+phi = (MMF/(M)); //flux produced in WB
+//phi = ((Rs+Rc)/Rs)*((K*theta)/N)
+Rs = (Rc*f*theta)/(phi-(f*theta)); //resistance of shunt in Ω
+
+//result
+mprintf("resistance of shunt = %3.2e Ω",Rs);
+
+
diff --git a/3871/CH13/EX13.6/Ex13_6.sce b/3871/CH13/EX13.6/Ex13_6.sce new file mode 100644 index 000000000..91a4c46bc --- /dev/null +++ b/3871/CH13/EX13.6/Ex13_6.sce @@ -0,0 +1,15 @@ +//=====================================================================================
+//Chapter 13 example 6
+
+clc;clear all;
+
+//variable declaration
+Rc = 1; //resistance in Ω
+N = 5; //multiplying factor
+
+//calculations
+//N = (Rs+Rc)/Rs
+Rs = Rc/(N-1); //shunt resistance in Ω
+
+//result
+mprintf("shunt resistance = %3.2f Ω",Rs);
diff --git a/3871/CH13/EX13.7/Ex13_7.sce b/3871/CH13/EX13.7/Ex13_7.sce new file mode 100644 index 000000000..96ca3a745 --- /dev/null +++ b/3871/CH13/EX13.7/Ex13_7.sce @@ -0,0 +1,24 @@ +//=====================================================================================
+//Chapter 13 example 7
+
+clc;clear all;
+
+//variable declaration
+R1 = 180; //resistance in Ω
+R2 = 20; //resistance in Ω
+A = 0.005; //area in m^2
+Ns = 1000; //number of turns on search coil
+G1 = 100*10^-6; //galvanometer constant C
+G2 = 100; //galvanometer throw
+
+//calculations
+Rs = R1+R2; //total resistance of secondary circuit in Ω
+Q = G1*G2; //charge passed through ballistic galvanometer in C
+//Q = i*t = (E/Rs)*t = ((2*phi*Ns)/(t*Rs))*t = (2*phi*Ns)/Rs
+phi = (Q*Rs)/(2*Ns); //flux in Wb
+B = phi/(A); //flux density in Wb/m^2
+
+//result
+
+mprintf("flux density = %3.2f Wb/m^2",B);
+
diff --git a/3871/CH13/EX13.8/Ex13_8.sce b/3871/CH13/EX13.8/Ex13_8.sce new file mode 100644 index 000000000..bba9c53ce --- /dev/null +++ b/3871/CH13/EX13.8/Ex13_8.sce @@ -0,0 +1,33 @@ +//=====================================================================================
+//Chapter 13 example 8
+
+clc;clear all;
+
+//variable declaration
+d = 0.1; //diameter in m
+a = 33.5*10^-6; //cross sectional area of iron ring in m^2
+Ns = 220; //number of turns on secondary coil
+Nm = 320; //number of turns on magnetising winding
+I = 10; //current in A
+B = 2.5*10^-3; //flux in Wb
+n = 102; //reading of scale
+g = 272; //galvanometer throw
+
+
+//calculations
+l = (%pi)*d; //mean length of iron ring in m
+H = (Nm*I)/(l); //magnetising force with 10 A current
+K = B/(n);
+//2*phi*Ns = K*g
+phi = (K*g)/(2*Ns); //flux in Wb
+B1 = phi/(a); //flux density in Wb/m**2
+u0 = 4*%pi*10^-7;
+//B = u0*ur*H
+x = u0*H;
+//B = x*ur
+//ur = B/x
+ur = B1/x; //relative permeability
+
+//result
+mprintf("relative permeability = %3.1f",ur);
+
diff --git a/3871/CH13/EX13.9/Ex13_9.sce b/3871/CH13/EX13.9/Ex13_9.sce new file mode 100644 index 000000000..b7989bc34 --- /dev/null +++ b/3871/CH13/EX13.9/Ex13_9.sce @@ -0,0 +1,25 @@ +//=====================================================================================
+//Chapter 13 example 9
+clc;clear all;
+
+//variable declaration
+R = 2000; //resistance in Ω
+l = 1; //mean length of iron ring in m
+A = 350*10^-6; //area in m**2
+Ns = 200; //number of turns on secondary coil
+G1 = 1*10^-6; //galvanometer constant C
+G2 = 100; //galvanometer throw
+N =100;
+
+//calculations
+u0 = 4*(%pi)*10^-7;
+H = (N*I)/(l); //magnetising force with 10 A current
+Q = G1*G2; //charge passed through ballistic galvanometer in C
+//Q = i*t = (E/Rs)*t = ((2*phi*Ns)/(t*Rs))*t = (2*phi*Ns)/Rs
+phi = (Q*R)/(2*Ns); //flux in Wb
+B = phi/(A); //flux density in Wb/m**2
+ur = (B/(u0*H)); //relative permeability
+
+//result
+mprintf("flux density = %3.3f Wb/m**2",B);
+mprintf("\nrelative permeability = %3.0f",ur);
|