From 7bc77cb1ed33745c720952c92b3b2747c5cbf2df Mon Sep 17 00:00:00 2001 From: prashantsinalkar Date: Sat, 3 Feb 2018 11:01:52 +0530 Subject: Added new code --- 3871/CH13/EX13.1/Ex13_1.sce | 21 +++++++++++++++++++++ 3871/CH13/EX13.10/Ex13_10.sce | 30 ++++++++++++++++++++++++++++++ 3871/CH13/EX13.11/Ex13_11.sce | 19 +++++++++++++++++++ 3871/CH13/EX13.12/Ex13_12.sce | 22 ++++++++++++++++++++++ 3871/CH13/EX13.13/Ex13_13.sce | 26 ++++++++++++++++++++++++++ 3871/CH13/EX13.14/Ex13_14.sce | 28 ++++++++++++++++++++++++++++ 3871/CH13/EX13.15/Ex13_15.sce | 17 +++++++++++++++++ 3871/CH13/EX13.16/Ex13_16.sce | 23 +++++++++++++++++++++++ 3871/CH13/EX13.17/Ex13_17.sce | 26 ++++++++++++++++++++++++++ 3871/CH13/EX13.2/Ex13_2.sce | 25 +++++++++++++++++++++++++ 3871/CH13/EX13.4/Ex13_4.sce | 20 ++++++++++++++++++++ 3871/CH13/EX13.5/Ex13_5.sce | 22 ++++++++++++++++++++++ 3871/CH13/EX13.6/Ex13_6.sce | 15 +++++++++++++++ 3871/CH13/EX13.7/Ex13_7.sce | 24 ++++++++++++++++++++++++ 3871/CH13/EX13.8/Ex13_8.sce | 33 +++++++++++++++++++++++++++++++++ 3871/CH13/EX13.9/Ex13_9.sce | 25 +++++++++++++++++++++++++ 16 files changed, 376 insertions(+) create mode 100644 3871/CH13/EX13.1/Ex13_1.sce create mode 100644 3871/CH13/EX13.10/Ex13_10.sce create mode 100644 3871/CH13/EX13.11/Ex13_11.sce create mode 100644 3871/CH13/EX13.12/Ex13_12.sce create mode 100644 3871/CH13/EX13.13/Ex13_13.sce create mode 100644 3871/CH13/EX13.14/Ex13_14.sce create mode 100644 3871/CH13/EX13.15/Ex13_15.sce create mode 100644 3871/CH13/EX13.16/Ex13_16.sce create mode 100644 3871/CH13/EX13.17/Ex13_17.sce create mode 100644 3871/CH13/EX13.2/Ex13_2.sce create mode 100644 3871/CH13/EX13.4/Ex13_4.sce create mode 100644 3871/CH13/EX13.5/Ex13_5.sce create mode 100644 3871/CH13/EX13.6/Ex13_6.sce create mode 100644 3871/CH13/EX13.7/Ex13_7.sce create mode 100644 3871/CH13/EX13.8/Ex13_8.sce create mode 100644 3871/CH13/EX13.9/Ex13_9.sce (limited to '3871/CH13') 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); -- cgit