6 files changed, 300 insertions, 0 deletions
diff --git a/1820/CH10/EX10.1/Example10_1.sce b/1820/CH10/EX10.1/Example10_1.sce
new file mode 100755
index 000000000..52fc514ff
--- /dev/null
+++ b/1820/CH10/EX10.1/Example10_1.sce
@@ -0,0 +1,25 @@
+// ELECTRIC POWER TRANSMISSION SYSTEM ENGINEERING ANALYSIS AND DESIGN
+// TURAN GONEN
+// CRC PRESS
+// SECOND EDITION
+
+// CHAPTER : 10 : PROTECTIVE EQUIPMENT AND TRANSMISSION SYSTEM PROTECTION
+
+// EXAMPLE : 10.1 :
+clear ; clc ; close ; // Clear the work space and console
+
+// GIVEN DATA
+X_d = 0.14*%i ; // Reactance of generator in pu
+E_g = 1*exp(%i*0*%pi/180) ;
+S_B = 25*10^3 ; // voltage in kVA
+V_BL_V = 13.8 ; // low voltage in kV
+
+// CALCULATIONS
+I_f = E_g/X_d ; // Subtransient fault current in pu
+I_BL_V = S_B/( sqrt(3)*V_BL_V) ; // Current base for low-voltage side
+I_f1 = abs(I_f)*I_BL_V ; // magnitude of fault current in A
+
+// DISPLAY RESULTS
+disp("EXAMPLE : 10.1 : SOLUTION :-") ;
+printf("\n Subtransient fault current for 3-Φ fault in per units = pu \n") ; disp(I_f) ;
+printf("\n Subtransient fault current for 3-Φ fault in ampere = %.f A \n",I_f1) ;
diff --git a/1820/CH10/EX10.2/Example10_2.sce b/1820/CH10/EX10.2/Example10_2.sce
new file mode 100755
index 000000000..3438dc066
--- /dev/null
+++ b/1820/CH10/EX10.2/Example10_2.sce
@@ -0,0 +1,42 @@
+// ELECTRIC POWER TRANSMISSION SYSTEM ENGINEERING ANALYSIS AND DESIGN
+// TURAN GONEN
+// CRC PRESS
+// SECOND EDITION
+
+// CHAPTER : 10 : PROTECTIVE EQUIPMENT AND TRANSMISSION SYSTEM PROTECTION
+
+// EXAMPLE : 10.2 :
+clear ; clc ; close ; // Clear the work space and console
+
+// GIVEN DATA
+// For case (a)
+I_f = 7.1428571 ; // Subtransient fault current in pu . Result of exa 10.1
+
+// For case (d)
+V_pf = 13800 ; // voltage in V
+zeta = 1.4 ;
+I_f1 = 7471 ; // magnitude of fault current in A
+
+// CALCULATIONS
+// For case (a)
+I_fdc_max = sqrt(2)*I_f ; // Max dc current in pu
+
+// For case (b)
+I_f_max = 2*I_fdc_max ; // Total max instantaneous current in pu
+
+// For case (c)
+I_momt = 1.6*I_f ; // Total rms momentary current
+
+// For case (d)
+S_int = sqrt(3)*(V_pf)*I_f1*zeta*10^-6 ; // Interrupting rating in MVA
+
+// For case (e)
+S_momt = sqrt(3)*(V_pf)*I_f1*1.6*10^-6 ; // Momentary duty of CB in MVA
+
+// DISPLAY RESULTS
+disp("EXAMPLE : 10.2 : SOLUTION :-") ;
+printf("\n (a) Maximum possible dc current component , I_fdc_max = %.1f pu \n",I_fdc_max) ;
+printf("\n (b) Total maximum instantaneous current , I_max = %.1f pu \n",I_f_max) ;
+printf("\n (c) Momentary current , I_momentary = %.2f pu \n",I_momt) ;
+printf("\n (d) Interrupting rating of a 2-cycle CB , S_interrupting = %.f MVA \n",S_int) ;
+printf("\n (e) Momentary duty of a 2-cycle CB , S_momentary = %.2f MVA \n",S_momt) ;
diff --git a/1820/CH10/EX10.4/Example10_4.sce b/1820/CH10/EX10.4/Example10_4.sce
new file mode 100755
index 000000000..2e166e257
--- /dev/null
+++ b/1820/CH10/EX10.4/Example10_4.sce
@@ -0,0 +1,41 @@
+// ELECTRIC POWER TRANSMISSION SYSTEM ENGINEERING ANALYSIS AND DESIGN
+// TURAN GONEN
+// CRC PRESS
+// SECOND EDITION
+
+// CHAPTER : 10 : PROTECTIVE EQUIPMENT AND TRANSMISSION SYSTEM PROTECTION
+
+// EXAMPLE : 10.4 :
+clear ; clc ; close ; // Clear the work space and console
+
+// GIVEN DATA
+z_l = 0.2 + %i * 0.7 ; // Line impedance in pu
+f_l = 0.7 ; // Fault point at a distance from A in pu
+f_m = 1.2 ; // magnitude of fault current in pu
+l = 10.3 ; // Line spacing in ft
+p = 100 ; // Power in MVA
+v = 138 ; // voltage in kV
+i = 418.4 ; // current in A
+z = 190.4 ; // Impedance in Ω
+
+// CALCULATIONS
+// For case (a)
+I = f_m * i ; // Current in arc in A
+R_arc = 8750 * l/(I^1.4) ; // Arc resistance in Ω
+R_arc1 = R_arc/z ; // Arc resistance in pu
+
+// For case (b)
+Z_L = z_l * f_l ;
+Z_r = Z_L + R_arc1 ; // Impedance seen by the relay in pu
+
+// For case (c)
+phi_1 = atand( imag(Z_L),real(Z_L) ) ; // Line impedance angle without arc resistance in degree
+phi_2 = atand( imag(Z_r),real(Z_r) ) ; // Line impedance angle with arc resistance in degree
+
+// DISPLAY RESULTS
+disp("EXAMPLE : 10.4 : SOLUTION :-") ;
+printf("\n (a) Value of arc resistance at fault point in Ω , R_arc = %.2f Ω \n",R_arc) ;
+printf("\n     Value of arc resistance at fault point in pu , R_arc = %.2f pu \n",R_arc1) ;
+printf("\n (b) Value of line impedance including the arc resistance , Z_L + R_arc = pu \n") ; disp(Z_r) ;
+printf("\n (c) Line impedance angle without arc resistance , Φ = %.2f degree \n",phi_1) ;
+printf("\n     Line impedance angle with arc resistance , Φ = %.2f degree \n",phi_2) ;
diff --git a/1820/CH10/EX10.5/Example10_5.sce b/1820/CH10/EX10.5/Example10_5.sce
new file mode 100755
index 000000000..a9361d6d5
--- /dev/null
+++ b/1820/CH10/EX10.5/Example10_5.sce
@@ -0,0 +1,104 @@
+// ELECTRIC POWER TRANSMISSION SYSTEM ENGINEERING ANALYSIS AND DESIGN
+// TURAN GONEN
+// CRC PRESS
+// SECOND EDITION
+
+// CHAPTER : 10 : PROTECTIVE EQUIPMENT AND TRANSMISSION SYSTEM PROTECTION
+
+// EXAMPLE : 10.5 :
+clear ; clc ; close ; // Clear the work space and console
+
+// CALCULATIONS
+// For case (a)
+// Coordinate Values taken here are only for reference . Refer exa 10.5
+
+T = 0:0.01:300 ;
+
+for i = 1:int(length(T)/1.1) ; 
+    po(i) = 4 ;
+end
+for i = int(length(T)/1.1):length(T)
+    po(i) = 5 ;
+end
+for i = 1:int(length(T)/1.1)
+    io(i) = 4 ;
+    end
+for i = int(length(T)/1.1):length(T)
+    io(i) = 3 ;
+end
+
+a= gca() ;
+subplot(2,1,1) ; // To plot 2 graph in same graphic window
+a.thickness = 2 ; // sets thickness of plot of points
+plot2d(T,po,3,'012','',[0 0 310 7]) ;
+plot2d(T,io,3,'012','',[0 0 310 7]) ;
+xtitle("Fig 10.5 (a)   Zones of protection for relay R_12") ;
+xset('thickness',2); // sets thickness of axes
+xstring(25,3.8,'[]') ;
+xstring(45,4.2,'(1)') ;
+plot(45,4,'+') ;
+xstring(60,3.8,'[]') ;
+xstring(60,4.2,'B_12') ;
+xstring(120,3.8,'[]') ;
+xstring(120,4.2,'B_21') ;
+xstring(140,4.2,'(2)') ;
+plot(140,4,'+') ;
+xstring(155,3.8,'[]') ;
+xstring(155,4.2,'B_23') ;
+xstring(220,3.8,'[]') ;
+xstring(220,4.2,'B_32') ;
+xstring(270,5.0,'(3)') ;
+xstring(285,2.8,'[]') ;
+xstring(285,3.2,'B_35') ;
+xstring(285,4.8,'[]') ;
+xstring(285,5.2,'B_34') ;
+xstring(85,3.4,'TL_12') ;
+xstring(180,3.4,'TL_23') ;
+xstring(60,3,'ZONE 1') ;
+xstring(100,2,'ZONE 2') ;
+xstring(190,1,'ZONE 3') ;
+
+// For case (b)
+
+for i = 1:int(length(T)/4) ;
+    vo(i) = 0.5;
+end
+for i = int(length(T)/4):length(T/1.7)
+    vo(i) = 2;
+end
+for i = int(length(T)/1.7):length(T)
+    vo(i) = 4
+end
+
+for i = int(length(T)/2.14):length(T/1.35) ; // plotting Voltage values
+    uo(i) = 0.5;
+end
+for i = int(length(T)/1.35):length(T)
+    uo(i) = 2;
+end
+
+a = gca() ;
+a.thickness = 2 ;
+subplot(2,1,2)
+plot2d(T,vo,2,'012','',[0 0 310 7]) ;
+plot2d(T,uo,2,'012','',[0 0 310 7]) ;
+ylabel("OPERATING TIME") ;
+xlabel("IMPEDANCE") ;
+xtitle("Fig 10.5 (b) Coordination of distance relays , Operating time v/s Impedance") ;
+xset('thickness',2); // sets thickness of axes
+xstring(0.1,0.3,'T_1') ;
+xstring(30,0.6,'R_12') ;
+xstring(58,1.3,'T_2') ;
+xstring(100,2.0,'R_12') ;
+xstring(160,3.0,'T_3') ;
+xstring(230,4.0,'R_12') ;
+xstring(160,0.6,'R_23') ;
+xstring(260,2.1,'R_23') ;
+
+// DISPLAY RESULTS
+disp("EXAMPLE : 10.5 : SOLUTION :-") ;
+printf("\n (a) The zone of protection for relay R_12 is shown in Fig 10.5 (a) \n") ;
+printf("\n  ZONE 1 lies b/w (1) & B_21 \n") ;
+printf("\n  ZONE 2 lies b/w (1) & TL_23 \n") ;
+printf("\n  ZONE 3 lies after (1) \n") ;
+printf("\n (b) The coordination of the distance relays R_12 & R_21 in terms of Operating time v/s Impedance is shown in Fig 10.5 (b)") ;
diff --git a/1820/CH10/EX10.6/Example10_6.sce b/1820/CH10/EX10.6/Example10_6.sce
new file mode 100755
index 000000000..07e13e031
--- /dev/null
+++ b/1820/CH10/EX10.6/Example10_6.sce
@@ -0,0 +1,55 @@
+// ELECTRIC POWER TRANSMISSION SYSTEM ENGINEERING ANALYSIS AND DESIGN
+// TURAN GONEN
+// CRC PRESS
+// SECOND EDITION
+
+// CHAPTER : 10 : PROTECTIVE EQUIPMENT AND TRANSMISSION SYSTEM PROTECTION
+
+// EXAMPLE : 10.6 :
+clear ; clc ; close ; // Clear the work space and console
+
+// GIVEN DATA
+kv = 230 * 10^3 ; // transmission system voltage in V
+VA = 100 * 10^6 ; // Maximum peak load supplied by TL_12 in VA
+ZTL_12 = 2 + %i * 20 ; // Positive-sequence impedances of line TL_12
+ZTL_23 = 2.5 + %i * 25 ; // Positive-sequence impedances of line TL_23
+pf = 0.9 ; // Lagging pf
+
+// CALCULATIONS 
+// For case (a)
+I_max = VA/(sqrt(3)*kv) ; // Maximum load current in A
+
+// For case (b)
+CT = 250/5 ; // CT ratio which gives about 5A in secondary winding under the maximum loading
+
+// For case (c)
+vr = 69 ; // selecting Secondary voltage of 69 V line to neutral
+VT = (kv/sqrt(3))/vr ; // Voltage ratio
+
+// For case (d)
+Z_r = CT/VT ; // impedance measured by relay . Z_r = (V/VT)/(I/CT)
+Z_TL_12 = Z_r * ZTL_12 ; // Impedance of lines TL_12 as seen by relay
+Z_TL_23 = Z_r * ZTL_23 ; // Impedance of lines TL_23 as seen by relay
+
+// For case (e)
+Z_load = vr * CT * (pf + %i*sind(acosd(pf)))/(I_max) ; // Load impedance based on secondary ohms
+
+// For case (f)
+Z_r1 = 0.80 * Z_TL_12 ; // Zone 1 setting of relay R_12
+
+// For case (g)
+Z_r2 = 1.20 * Z_TL_12 ; // Zone 2 setting of relay R_12
+
+// For case (h)
+Z_r3 = Z_TL_12 + 1.20*(Z_TL_23) ; // Zone 3 setting of relay R_12
+
+// DISPLAY RESULTS
+disp("EXAMPLE : 10.6 : SOLUTION :-") ;
+printf("\n (a) Maximum load current , I_max = %.2f A \n",I_max) ;
+printf("\n (b) CT ratio , CT = %.1f \n",CT) ;
+printf("\n (c) VT ratio , VT = %.1f \n",VT) ;
+printf("\n (d) Impedance measured by relay  = %.3f Z_line \n",Z_r) ;
+printf("\n (e) Load impedance based on secondary ohms , Z_load = Ω(secondary) \n") ; disp(Z_load) ;
+printf("\n (f) Zone 1 setting of relay R_12 , Z_r = Ω(secondary) \n") ; disp(Z_r1) ;
+printf("\n (g) Zone 2 setting of relay R_12 , Z_r = Ω(secondary) \n") ; disp(Z_r2) ;
+printf("\n (h) Zone 3 setting of relay R_12 , Z_r = Ω(secondary) \n") ; disp(Z_r3) ;
diff --git a/1820/CH10/EX10.7/Example10_7.sce b/1820/CH10/EX10.7/Example10_7.sce
new file mode 100755
index 000000000..13610f8ae
--- /dev/null
+++ b/1820/CH10/EX10.7/Example10_7.sce
@@ -0,0 +1,33 @@
+// ELECTRIC POWER TRANSMISSION SYSTEM ENGINEERING ANALYSIS AND DESIGN
+// TURAN GONEN
+// CRC PRESS
+// SECOND EDITION
+
+// CHAPTER : 10 : PROTECTIVE EQUIPMENT AND TRANSMISSION SYSTEM PROTECTION
+
+// EXAMPLE : 10.7 :
+clear ; clc ; close ; // Clear the work space and console
+
+// GIVEN DATA
+Z_r1 = 0.0415692 + %i*0.4156922 ; // Required zone 1 setting . From result of exa 10.6
+Z_r2 = 0.0623538 + %i*0.6235383 ; // Required zone 2 setting . From result of exa 10.6
+Z_r3 = 0.1299038 + %i*1.2990381 ; // Required zone 3 setting . From result of exa 10.6
+
+// CALCULATIONS 
+// For case (a)
+theta1 = atand(imag(Z_r1),real(Z_r1)) ;
+Z_1 = abs(Z_r1)/cosd(theta1 - 30) ; // Zone 1 setting of mho relay R_12
+
+// For case (b)
+theta2 = atand(imag(Z_r2),real(Z_r2)) ;
+Z_2 = abs(Z_r2)/cosd(theta2 - 30) ; // Zone 2 setting of mho relay R_12
+
+// For case (b)
+theta3 = atand(imag(Z_r3),real(Z_r3)) ;
+Z_3 = abs(Z_r3)/cosd(theta3 - 30) ; // Zone 3 setting of mho relay R_12
+
+// DISPLAY RESULTS
+disp("EXAMPLE : 10.7 : SOLUTION :-") ;
+printf("\n (a) Zone 1 setting of mho relay R_12 = %.4f Ω(secondary) \n",Z_1) ;
+printf("\n (b) Zone 2 setting of mho relay R_12 = %.4f Ω(secondary) \n",Z_2) ;
+printf("\n (c) Zone 3 setting of mho relay R_12 = %.4f Ω(secondary) \n",Z_3) ;