initial commit / add all books

14 files changed, 703 insertions, 0 deletions

diff --git a/1820/CH9/EX9.1/Example9_1.sce b/1820/CH9/EX9.1/Example9_1.sce
new file mode 100755
index 000000000..8c0e92a2b
--- /dev/null
+++ b/1820/CH9/EX9.1/Example9_1.sce
@@ -0,0 +1,42 @@
+// ELECTRIC POWER TRANSMISSION SYSTEM ENGINEERING ANALYSIS AND DESIGN

+// TURAN GONEN

+// CRC PRESS

+// SECOND EDITION

+

+// CHAPTER : 9 : SYMMETRICAL COMPONENTS AND FAULT ANALYSIS

+

+// EXAMPLE : 9.1 :

+clear ; clc ; close ; // Clear the work space and console

+

+// GIVEN DATA

+V_a = 7.3 * exp(%i*12.5*%pi/180) ; // Phase voltage in V

+V_b = 0.4 * exp(%i*(-100)*%pi/180) ; // Phase voltage in V

+V_c = 4.4 * exp(%i*154*%pi/180) ; // Phase voltage in V

+a = 1 * exp(%i*120*%pi/180) ; // operator 'a' by application of symmetrical components theory to 3-Φ system . Refer section 9.3 for details

+

+// CALCULATIONS

+V_a0 = (1/3) * (V_a + V_b + V_c) ; // Analysis equ in V

+V_a1 = (1/3) * (V_a + a*V_b + a^2*V_c) ;

+V_a2 = (1/3) * (V_a + a^2*V_b + a*V_c) ;

+V_b0 = V_a0 ;

+V_b1 = a^2 * V_a1 ;

+V_b2 = a * V_a2 ;

+V_c0 = V_a0 ;

+V_c1 = a * V_a1 ;

+V_c2 = a^2 * V_a2 ;

+

+// DISPLAY RESULTS

+disp("EXAMPLE : 9.1 : SOLUTION :-") ;

+printf("\n The symmetrical components for the phase voltages V_a , V_b & V_c are\n") ;

+printf("\n V_a0 = %.2f<%.1f V \n",abs(V_a0),atand( imag(V_a0),real(V_a0) )) ;

+printf("\n V_a1 = %.2f<%.1f V \n",abs(V_a1),atand( imag(V_a1),real(V_a1) )) ;

+printf("\n V_a2 = %.2f<%.1f V \n",abs(V_a2),atand( imag(V_a2),real(V_a2) )) ;

+printf("\n V_b0 = %.2f<%.1f V \n",abs(V_b0),atand( imag(V_b0),real(V_b0) )) ;

+printf("\n V_b1 = %.2f<%.1f V \n",abs(V_b1),atand( imag(V_b1),real(V_b1) )) ;

+printf("\n V_b2 = %.2f<%.1f V \n",abs(V_b2),atand( imag(V_b2),real(V_b2) )) ;

+printf("\n V_c0 = %.2f<%.1f V \n",abs(V_c0),atand( imag(V_c0),real(V_c0) )) ;

+printf("\n V_c1 = %.2f<%.1f V \n",abs(V_c1),atand( imag(V_c1),real(V_c1) )) ;

+printf("\n V_c2 = %.2f<%.1f V \n",abs(V_c2),atand( imag(V_c2),real(V_c2) )) ;

+

+printf("\n NOTE : V_b1 = 3.97<-99.5 V & V_c2 = 2.52<-139.7 V result obtained is same as textbook answer V_b1 = 3.97<260.5 V & V_c2 = 2.52<220.3 V \n") ;

+printf("\n Changes is due to a^2 = 1<240 = 1<-120  where 1 is the magnitude & <240 is the angle in degree \n") ;

diff --git a/1820/CH9/EX9.10/Example9_10.sce b/1820/CH9/EX9.10/Example9_10.sce
new file mode 100755
index 000000000..5782c9c3e
--- /dev/null
+++ b/1820/CH9/EX9.10/Example9_10.sce
@@ -0,0 +1,84 @@
+// ELECTRIC POWER TRANSMISSION SYSTEM ENGINEERING ANALYSIS AND DESIGN

+// TURAN GONEN

+// CRC PRESS

+// SECOND EDITION

+

+// CHAPTER : 9 : SYMMETRICAL COMPONENTS AND FAULT ANALYSIS

+

+// EXAMPLE : 9.10 :

+clear ; clc ; close ; // Clear the work space and console

+

+// GIVEN DATA

+Z_0 = 0.2619 * %i ;

+Z_1 = 0.25 * %i ;

+Z_2 = 0.25 * %i ;

+v = 1 * exp(%i*0*%pi/180) ;

+a = 1 * exp(%i*120*%pi/180) ; // By symmetrical components theory to 3-Φ system

+A = [1 1 1; 1 a^2 a ;1 a a^2] ;

+

+// CALCULATIONS

+// For case (b)

+I_a0 = v/(Z_0 + Z_1 + Z_2) ; // Sequence currents at fault point F in pu A

+I_a1 = I_a0 ;

+I_a2 = I_a0 ;

+

+// For case (c)

+I_a1g1 = (1/2) * I_a1 ; // Sequence current at terminals of generator G1 in pu A

+I_a2g1 = (1/2) * I_a2 ;

+I_a0g1 = 0.5/(0.55 + 0.5)*I_a0 ; // By current division in pu A

+

+// For case (d)

+I_f = [A] * [I_a0g1 ; I_a1g1 ; I_a2g1] ; // Phase current at terminal of generator G1 in pu A

+

+// For case (e)

+V_a = [0 ; v ; 0] - [Z_0 0 0 ; 0 Z_1 0 ; 0 0 Z_2]*[I_a0g1 ; I_a1g1 ; I_a2g1] ; // Sequence voltage in pu V

+

+// For case (f)

+V_f = [A]*[V_a] ; // Phase voltage at terminal of generator G1 in pu V

+

+// For case (g)

+I_a1g2 = (1/2) * I_a1 ; // By symmetry for Generator G2

+I_a2g2 = (1/2) * I_a2 ;

+I_a0g2 = 0 ; // By inspection

+// V_a1(HV) leads V_a1(LV) by 30 degree & V_a2(HV) lags V_a2(LV) by 30 degree

+I_a0G2 = I_a0g2 ;

+I_a1G2 = abs(I_a1g2)*exp(%i * (atand( imag(I_a1g2),real(I_a1g2) ) - 30) * %pi/180) ; // (-90-30) = (-120)

+I_a2G2 = abs(I_a2g2)*exp(%i *(atand( imag(I_a2g2),real(I_a2g2) ) + 30) * %pi/180)  ; // (-90+30) = (-60)

+

+I_f2 = [A] * [I_a0G2 ; I_a1G2 ; I_a2G2] ; // Phase current at terminal of generator G2 in pu A

+

+ // Sequence voltage at terminal of generator G2 in pu V

+V_a0G2 = 0 ;

+V_a1G2 = abs(V_a(2,1))*exp(%i * (atand( imag(V_a(2,1)),real(V_a(2,1)) ) - 30) * %pi/180) ; // (0-30) = (-30)

+V_a2G2 = abs(V_a(3,1))*exp(%i * (atand( imag(V_a(3,1)),real(V_a(3,1)) ) + 30) * %pi/180) ; // (180+30)=(210)=(-150)

+

+V_f2 = A * [V_a0G2 ; V_a1G2 ; V_a2G2] ; // Phase voltage at terminal of generator G2 in pu V

+

+// DISPLAY RESULTS

+disp("EXAMPLE : 9.10 : SOLUTION :-") ;

+printf("\n (b) The sequence current at fault point F , I_a0 = I_a1 = I_a2 = %.4f<%.f pu A \n",abs(I_a0),atand(imag(I_a0),real(I_a0) )) ;

+printf("\n (c) Sequence currents at the terminals of generator G1 , \n") ;

+printf("\n     I_a0,G_1 = %.4f<%.f pu A ",abs(I_a0g1),atand( imag(I_a0g1),real(I_a0g1) )) ;

+printf("\n     I_a1,G_1 = %.4f<%.f pu A ",abs(I_a1g1),atand( imag(I_a1g1),real(I_a1g1) )) ;

+printf("\n     I_a2,G_1 = %.4f<%.f pu A ",abs(I_a2g1),atand( imag(I_a2g1),real(I_a2g1) )) ;

+printf("\n \n (d) Phase currents at terminal of generator G1 are , [I_af ; I_bf ; I_cf] = pu A \n") ;

+printf("\n       %.4f<%.f ",abs(I_f),atand(imag(I_f),real(I_f) )) ;

+printf("\n \n (e) Sequence voltages at the terminals of generator G1 , [V_a0 ; V_a1 ; V_a2 ] = pu V \n") ;

+printf("\n       %.4f<%.1f ",abs(V_a),atand(imag(V_a),real(V_a) )) ; 

+printf("\n \n (f) Phase voltages at terminal of generator G1 are , [V_af ; V_bf ; V_cf] = pu V \n") ;

+printf("\n       %.4f<%.1f ",abs(V_f),atand(imag(V_f),real(V_f) )) ; 

+printf("\n \n (g) Sequence currents at the terminals of generator G2 , \n") ;

+printf("\n     I_a0,G_2 = %.f<%.f pu A ",abs(I_a0G2),atand( imag(I_a0G2),real(I_a0G2) )) ;

+printf("\n     I_a1,G_2 = %.4f<%.f pu A",abs(I_a1G2),atand( imag(I_a1G2),real(I_a1G2) )) ;

+printf("\n     I_a2,G_2 = %.4f<%.f pu A",abs(I_a2G2),atand( imag(I_a2G2),real(I_a2G2) )) ;

+printf("\n \n     Phase currents at terminal of generator G2 are , [I_af ; I_bf ; I_cf] = pu A \n") ;

+printf("\n       %.4f<%.f ",abs(I_f2),atand(imag(I_f2),real(I_f2) )) ;

+printf("\n \n     Sequence voltages at the terminals of generator G2 , [V_a0 ; V_a1 ; V_a2 ] = pu V\n") ;

+printf("\n        %.f<%.f ",abs(V_a0G2),atand( imag(V_a0G2),real(V_a0G2) )) ;

+printf("\n        %.4f<%.f ",abs(V_a1G2),atand( imag(V_a1G2),real(V_a1G2) )) ;

+printf("\n        %.4f<%.f ",abs(V_a2G2),atand( imag(V_a2G2),real(V_a2G2) )) ;

+printf("\n \n     Phase voltages at terminal of generator G2 are , [V_af ; V_bf ; V_cf] = pu V \n") ;

+printf("\n        %.4f<%.1f ",abs(V_f2),atand(imag(V_f2),real(V_f2) )) ; 

+

+printf("\n \n NOTE : ERROR : Calclation mistake in textbook case(f) ") ;

+printf("\n In case (g) V_a2 = 0.1641<-150 is same as textbook answer V_a2 = 0.1641<210 , i.e (360-150)=210 \n") ;

diff --git a/1820/CH9/EX9.11/Example9_11.sce b/1820/CH9/EX9.11/Example9_11.sce
new file mode 100755
index 000000000..31f34cbbb
--- /dev/null
+++ b/1820/CH9/EX9.11/Example9_11.sce
@@ -0,0 +1,58 @@
+// ELECTRIC POWER TRANSMISSION SYSTEM ENGINEERING ANALYSIS AND DESIGN

+// TURAN GONEN

+// CRC PRESS

+// SECOND EDITION

+

+// CHAPTER : 9 : SYMMETRICAL COMPONENTS AND FAULT ANALYSIS

+

+// EXAMPLE : 9.11 :

+clear ; clc ; close ; // Clear the work space and console

+

+// GIVEN DATA

+kv = 230 ; // Line voltage in kV from Exa 9.9

+Z_0 = 0.56*%i ; // Zero-sequence impedance in pu

+Z_1 = 0.2618*%i ; // Zero-sequence impedance in pu

+Z_2 = 0.3619*%i ; // Zero-sequence impedance in pu

+z_f = 5 ; // Fault impedance in Ω

+v = 1*exp(%i*0*%pi/180) ; // 

+a = 1*exp(%i*120*%pi/180) ; // By symmetrical components theory to 3-Φ system

+A = [1 1 1; 1 a^2 a ;1 a a^2] ;

+

+// CALCULATIONS

+// For case (b)

+I_a0 = 0 ; // Sequence current in A

+Z_B = kv^2/200 ; // Base impedance of 230 kV line

+Z_f = z_f/Z_B ; // fault impedance in pu

+I_a1 = v/(Z_1 + Z_2 + Z_f) ; // Sequence current in pu A

+I_a2 = - I_a1 ; // Sequence current in pu A

+I_f = [A] * [I_a0 ; I_a1 ; I_a2] ; // Phase current in pu A

+

+// For case (c)

+V_a = [0 ; v ; 0]-[Z_0 0 0 ; 0 Z_1 0 ; 0 0 Z_2]*[I_a0 ; I_a1 ; I_a2] ; // Sequence voltages in pu V

+V_f = A*V_a ; // Phase voltages in pu V

+

+// For case (d)

+V_abf = V_f(1,1) - V_f(2,1) ; // Line-to-line voltages at fault points in pu V

+V_bcf = V_f(2,1) - V_f(3,1) ; // Line-to-line voltages at fault points in pu V

+V_caf = V_f(3,1) - V_f(1,1) ; // Line-to-line voltages at fault points in pu V

+

+

+

+// DISPLAY RESULTS

+disp("EXAMPLE : 9.11 :SOLUTION :-") ;

+printf("\n (b) Sequence currents are , \n") ;

+printf("\n  I_a0 = %.f pu A ",I_a0) ;

+printf("\n  I_a1 = %.4f<%.2f pu A ",abs(I_a1),atand( imag(I_a1),real(I_a1) )) ;

+printf("\n  I_a2 = %.4f<%.2f pu A ",abs(I_a2),atand( imag(I_a2),real(I_a2) )) ;

+printf("\n \n Phase currents are , [I_af ; I_bf ; I_cf] =  pu A \n") ;

+printf("\n       %.4f<%.1f ",abs(I_f),atand(imag(I_f),real(I_f) )) ;

+printf("\n \n (c) Sequence voltages are , [V_a0 ; V_a1 ; V_a2] =  pu V \n") ;

+printf("\n       %.4f<%.1f ",abs(V_a),atand(imag(V_a),real(V_a) )) ;

+printf("\n \n  Phase voltages are , [V_af ; V_bf ; V_cf] =  pu V \n") ;

+printf("\n       %.4f<%.1f ",abs(V_f),atand(imag(V_f),real(V_f) )) ;

+printf("\n \n (d) Line-to-line voltages at the fault points are \n") ;

+printf("\n     V_abf = %.4f<%.1f pu V \n",abs(V_abf),atand( imag(V_abf),real(V_abf) )) ;

+printf("\n     V_bcf = %.4f<%.1f pu V \n",abs(V_bcf),atand( imag(V_bcf),real(V_bcf) )) ;

+printf("\n     V_caf = %.4f<%.1f pu V \n",abs(V_caf),atand( imag(V_caf),real(V_caf) )) ;

+

+printf("\n \n NOTE : ERROR : Minor calclation mistake in textbook ") ;

diff --git a/1820/CH9/EX9.12/Example9_12.sce b/1820/CH9/EX9.12/Example9_12.sce
new file mode 100755
index 000000000..b887a9f22
--- /dev/null
+++ b/1820/CH9/EX9.12/Example9_12.sce
@@ -0,0 +1,56 @@
+// ELECTRIC POWER TRANSMISSION SYSTEM ENGINEERING ANALYSIS AND DESIGN

+// TURAN GONEN

+// CRC PRESS

+// SECOND EDITION

+

+// CHAPTER : 9 : SYMMETRICAL COMPONENTS AND FAULT ANALYSIS

+

+// EXAMPLE : 9.12 :

+clear ; clc ; close ; // Clear the work space and console

+

+// GIVEN DATA

+z_f = 5 ; // Fault-impedance in Ω

+z_g = 10 ; // Ground-impedance in Ω

+kv = 230 ; // Line voltage in kV from Exa 9.9

+Z_0 = 0.56*%i ; // Zero impedance in pu Ω

+Z_1 = 0.2618*%i ; // Positive sequence Impedance in pu Ω

+Z_2 = 0.3619*%i ; // Negative sequence Impedance in pu Ω

+v = 1*exp(%i*0*180/%pi) ;

+a = 1*exp(%i*120*%pi/180) ; // By symmetrical components theory to 3-Φ system

+A = [1 1 1; 1 a^2 a ;1 a a^2] ;

+

+// CALCULATIONS

+// For case (b)

+Z_B = kv^2/200 ; // Base impedance of 230 kV line

+Z_f = z_f/Z_B ; // fault impedance in pu Ω

+Z_g = z_g/Z_B ;

+I_a1 = v/( (Z_1 + Z_f) + ( (Z_2 + Z_f)*(Z_0 + Z_f + 3*Z_g)/((Z_2 + Z_f)+(Z_0 + Z_f + 3*Z_g)) )) ; // Sequence current in pu A

+I_a2 = -[(Z_0 + Z_f + 3*Z_g)/( (Z_2 + Z_f )+(Z_0 + Z_f + 3*Z_g) )]*I_a1 ; // Sequence current in pu A

+I_a0 = -[(Z_2 + Z_f)/( (Z_2 + Z_f)+(Z_0 + Z_f + 3*Z_g) )]*I_a1 ; // Sequence current in pu A

+I_f = A*[I_a0 ; I_a1 ; I_a2] ; // Phase currents in pu A

+

+// For case (c)

+V = [0 ; v ; 0] - [Z_0 0 0 ; 0 Z_1 0 ; 0 0 Z_2]*[I_a0 ; I_a1 ; I_a2] ; // Sequence Voltages in pu V

+V_f = A*[V] ; // Phase voltages in pu V

+

+// For case (d)

+V_abf = V_f(1,1) - V_f(2,1) ; // Line-to-line voltages at fault points a & b

+V_bcf = V_f(2,1) - V_f(3,1) ; // Line-to-line voltages at fault points b & c

+V_caf = V_f(3,1) - V_f(1,1) ; // Line-to-line voltages at fault points c & a

+

+// DISPLAY RESULTS

+disp("EXAMPLE : 9.12 : SOLUTION :-") ;

+printf("\n (b) Sequence currents are , \n") ;

+printf("\n   I_a0 = %.4f<%.2f pu A ",abs(I_a0),atand( imag(I_a0),real(I_a0) )) ;

+printf("\n   I_a1 = %.4f<%.2f pu A ",abs(I_a1),atand( imag(I_a1),real(I_a1) )) ;

+printf("\n   I_a2 = %.4f<%.2f pu A ",abs(I_a2),atand( imag(I_a2),real(I_a2) )) ;

+printf("\n \n   Phase currents are , [I_af ; I_bf ; I_cf] = pu A \n ") ;

+printf("\n       %.4f<%.1f ",abs(I_f),atand(imag(I_f),real(I_f) )) ;

+printf("\n \n (c) Sequence voltages , [V_a0 ; V_a1 ; V_a2] = pu V \n ") ;

+printf("\n       %.4f<%.1f ",abs(V),atand(imag(V),real(V) )) ;

+printf("\n \n  Phase voltages , [V_af ; V_bf ; V_cf] = pu V \n ") ;

+printf("\n       %.4f<%.1f ",abs(V_f),atand(imag(V_f),real(V_f) )) ;

+printf("\n \n (d) Line-to-line voltages at the fault points are , \n") ;

+printf("\n   V_abf = %.4f<%.1f pu V \n",abs(V_abf),atand( imag(V_abf),real(V_abf) )) ;

+printf("\n   V_bcf = %.4f<%.1f pu V \n",abs(V_bcf),atand( imag(V_bcf),real(V_bcf) )) ;

+printf("\n   V_caf = %.4f<%.1f pu V \n",abs(V_caf),atand( imag(V_caf),real(V_caf) )) ;

diff --git a/1820/CH9/EX9.13/Example9_13.sce b/1820/CH9/EX9.13/Example9_13.sce
new file mode 100755
index 000000000..2d314373a
--- /dev/null
+++ b/1820/CH9/EX9.13/Example9_13.sce
@@ -0,0 +1,56 @@
+// ELECTRIC POWER TRANSMISSION SYSTEM ENGINEERING ANALYSIS AND DESIGN

+// TURAN GONEN

+// CRC PRESS

+// SECOND EDITION

+

+// CHAPTER : 9 : SYMMETRICAL COMPONENTS AND FAULT ANALYSIS

+

+// EXAMPLE : 9.13 :

+clear ; clc ; close ; // Clear the work space and console

+

+// GIVEN DATA

+z_f = 5 ; // Fault-impedance in Ω

+Z_0 = 0.56*%i ; // Zero impedance in pu Ω

+Z_1 = 0.2618*%i ; // Positive sequence Impedance in pu Ω

+Z_2 = 0.3619*%i ; // Negative sequence Impedance in pu Ω

+kv = 230 ; // Line voltage in kV from Exa 9.9

+a = 1 * exp(%i*120*%pi/180) ; // By symmetrical components theory to 3-Φ system

+A = [1 1 1; 1 a^2 a ;1 a a^2] ;

+

+// CALCULATIONS

+// For case (b)

+Z_B = kv^2/200 ; // Base impedance of 230 kV line

+Z_f = z_f/Z_B ; // fault impedance in pu Ω

+v = 1*exp(%i*0*%pi/180) ;

+I_a0 = 0 ; // Sequence current in pu A

+I_a1 = v/(Z_1 + Z_f) ; // Sequence current in pu A

+I_a2 = 0 ; // Sequence current in pu A

+I_f = A*[I_a0 ; I_a1 ; I_a2] ; // Phase-current in pu A

+

+// For case (c)

+V = [0 ; v ; 0] - [Z_0 0 0 ; 0 Z_1 0 ; 0 0 Z_2]*[I_a0 ; I_a1 ; I_a2] ; // Sequence Voltages in pu V

+V_f = A*[V] ; // Phase voltages in pu V

+

+// For case (d)

+V_abf = V_f(1,1) - V_f(2,1) ; // Line-to-line voltages at fault points a & b

+V_bcf = V_f(2,1) - V_f(3,1) ; // Line-to-line voltages at fault points b & c

+V_caf = V_f(3,1) - V_f(1,1) ; // Line-to-line voltages at fault points c & a

+

+// DISPLAY RESULTS

+disp("EXAMPLE : 9.13 : SOLUTION :-") ;

+printf("\n (b) Sequence currents are , \n") ;

+printf("\n     I_a0 = %.1f pu A ",I_a0) ;

+printf("\n     I_a1 = %.4f<%.1f pu A ",abs(I_a1),atand( imag(I_a1),real(I_a1) )) ;

+printf("\n     I_a2 = %.1f pu A ",I_a2) ;

+printf("\n \n   Phase currents are , [I_af ; I_bf ; I_cf] = pu A \n ") ;

+printf("\n       %.4f<%.1f ",abs(I_f),atand(imag(I_f),real(I_f) )) ;

+printf("\n \n (c) Sequence voltages , [V_a0 ; V_a1 ; V_a2] = pu V \n ") ;

+printf("\n       %.4f<%.1f ",abs(V),atand(imag(V),real(V) )) ;

+printf("\n \n    Phase voltages , [V_af ; V_bf ; V_cf] = pu V \n ") ;

+printf("\n       %.4f<%.1f ",abs(V_f),atand(imag(V_f),real(V_f) )) ;

+printf("\n \n (d) Line-to-line voltages at the fault points are , \n") ;

+printf("\n    V_abf = %.4f<%.1f pu V \n",abs(V_abf),atand( imag(V_abf),real(V_abf) )) ;

+printf("\n    V_bcf = %.4f<%.1f pu V \n",abs(V_bcf),atand( imag(V_bcf),real(V_bcf) )) ;

+printf("\n    V_caf = %.4f<%.1f pu V \n",abs(V_caf),atand( imag(V_caf),real(V_caf) )) ;

+

+printf("\n \n NOTE : ERROR : Calclation mistake in textbook case(d) ") ;

diff --git a/1820/CH9/EX9.14/Example9_14.sce b/1820/CH9/EX9.14/Example9_14.sce
new file mode 100755
index 000000000..6fe7597fc
--- /dev/null
+++ b/1820/CH9/EX9.14/Example9_14.sce
@@ -0,0 +1,40 @@
+// ELECTRIC POWER TRANSMISSION SYSTEM ENGINEERING ANALYSIS AND DESIGN

+// TURAN GONEN

+// CRC PRESS

+// SECOND EDITION

+

+// CHAPTER : 9 : SYMMETRICAL COMPONENTS AND FAULT ANALYSIS

+

+// EXAMPLE : 9.14 :

+clear ; clc ; close ; // Clear the work space and console

+

+// GIVEN DATA

+VG_1 = 1*exp(%i*0*%pi/180) ;

+VG_2 = 1*exp(%i*0*%pi/180) ;

+

+// CALCULATIONS

+// For case (a)

+I_1 = 1*exp(%i*0*%pi/180) ;

+I_2 = 1*exp(%i*0*%pi/180) ;

+V_1 = 0.4522*exp(%i*90*%pi/180) ;

+V_2 = 0.4782*exp(%i*90*%pi/180) ;

+Y_11 = I_1/V_1 ; // When V_2 = 0 

+Y_21 = (-0.1087)*Y_11 ; // When V_2 = 0 

+Y_22 = I_2/V_2 ; // When V_1 = 0 

+Y_12 = Y_21 ;

+Y = [Y_11 Y_12 ; Y_21 Y_22] ; // Admittance matrix associated with positive-sequence n/w

+

+// For case (b)

+I_S1_12 = 2.0193*exp(%i*90*%pi/180) ; // Short-ckt F & F' to neutral & by superposition theorem

+I_S1_10 = 0.2884*exp(%i*90*%pi/180) ; // Short-ckt F & F' to neutral & by superposition theorem

+I_S2_12 = 0.4326*exp(%i*90*%pi/180) ;

+I_S2_10 = 1.4904*exp(%i*90*%pi/180) ;

+I_S1 = I_S1_12 + I_S1_10 ;

+I_S2 = I_S2_12 + I_S2_10 ;

+

+// DISPLAY RESULTS

+disp("EXAMPLE : 9.14 :SOLUTION :-") ;

+printf("\n (a) Admittance matrix associated with positive-sequence network , Y = \n") ; disp(Y) ;

+printf("\n (b) Source currents Two-port Thevenin equivalent positive sequence network are , \n") ;

+printf("\n     I_S1 = %.4f<%.f pu ",abs(I_S1),atand( imag(I_S1),real(I_S1) )) ;

+printf("\n     I_S2 = %.4f<%.f pu \n",abs(I_S2),atand( imag(I_S2),real(I_S2) )) ;

diff --git a/1820/CH9/EX9.15/Example9_15.sce b/1820/CH9/EX9.15/Example9_15.sce
new file mode 100755
index 000000000..e4ccd682a
--- /dev/null
+++ b/1820/CH9/EX9.15/Example9_15.sce
@@ -0,0 +1,37 @@
+

+// ELECTRIC POWER TRANSMISSION SYSTEM ENGINEERING ANALYSIS AND DESIGN

+// TURAN GONEN

+// CRC PRESS

+// SECOND EDITION

+

+// CHAPTER : 9 : SYMMETRICAL COMPONENTS AND FAULT ANALYSIS

+

+// EXAMPLE : 9.15 :

+clear ; clc ; close ; // Clear the work space and console

+

+// GIVEN DATA

+Y_11 = -2.2115*%i ;

+Y_12 = 0.2404*%i ;

+Y_21 = 0.2404*%i ;

+Y_22 = -2.0912*%i ;

+Y = [Y_11 Y_12 ; Y_21 Y_22] ;

+I_S1 = 2.3077*%i ;

+I_S2 = 1.9230*%i ;

+

+I_a1 = poly(0,'I_a1') ;

+I_a2 = poly(0,'I_a2') ;

+a = Y_12*I_S2 - Y_22*I_S1 ;

+b = (Y_12+Y_22)*I_a1 ;

+c = Y_12*I_S1 - Y_11*I_S2 ;

+d = (Y_12 + Y_11)*I_a1 ;

+V1 = (1/det(Y))*[(a-b) ; (c+d)] ; // Gives the uncoupled positive sequence N/W

+A = (Y_12+Y_22)*I_a2 ;

+B = (Y_12 + Y_11)*I_a2 ;

+V2 = (1/det(Y))*[A ; B] ; // Gives the uncoupled negative sequence N/W

+

+// DISPLAY RESULTS

+disp("EXAMPLE : 9.15 : SOLUTION :-") ;

+printf("\n (a) [V_a1 ; V_a11] = ") ; disp(V1) ;

+printf("\n     Values of Uncoupled positive-sequence network \n") ;

+printf("\n (b) [V_a2 ; V_a22] = ") ; disp(V2) ;

+printf("\n     Values of Uncoupled negative-sequence network \n") ;

diff --git a/1820/CH9/EX9.16/Example9_16.sce b/1820/CH9/EX9.16/Example9_16.sce
new file mode 100755
index 000000000..a30c5b2d6
--- /dev/null
+++ b/1820/CH9/EX9.16/Example9_16.sce
@@ -0,0 +1,58 @@
+// ELECTRIC POWER TRANSMISSION SYSTEM ENGINEERING ANALYSIS AND DESIGN

+// TURAN GONEN

+// CRC PRESS

+// SECOND EDITION

+

+// CHAPTER : 9 : SYMMETRICAL COMPONENTS AND FAULT ANALYSIS

+

+// EXAMPLE : 9.16 :

+clear ; clc ; close ; // Clear the work space and console

+

+// GIVEN DATA

+H_aa = 81.5 ;

+D_aa = 1.658 ;

+f = 60 ; // Freq in Hz

+I = 20 ;

+kV = 69 ; // Line voltage in kV

+MVA = 25 ; // Transformer T1 rating in MVA

+

+// CALCULATIONS

+// For case (a)

+C_0 = 29.842*10^-9/(log(H_aa/D_aa)) ; // Capacitance in F/mi

+b_0 = 2*%pi*f*C_0 ; // Susceptance in S/mi

+B_0 = b_0*I ; // For total system

+X_C0 = (1/B_0) ; // Total zero-sequence reactance in Ω

+TC_0 = B_0/(2*%pi*f) ; // Total zero-sequence capacitance in F

+

+// For case (c)

+X_1 = 0.05 ; // Leakage reactance of transformer T1 in pu

+X_0 = X_1 ;

+X_2 = X_1 ;

+Z_B = kV^2/MVA ;

+X_01 = X_0*Z_B ; // Leakage reactance in Ω

+V_F = 69*10^3/sqrt(3) ;

+I_a0PC = V_F/(17310.8915*%i) ; // Zero-sequence current flowing through PC in A

+I_PC = 3*abs(I_a0PC) ; // Continuous-current rating of the PC in A

+

+// For case (d)

+X_PC = (17310.8915 - X_01)/3 ; // Required reactance value for PC in Ω

+

+// For case (e)

+L_PC = X_PC/(2*%pi*f) ; // Inductance in H

+

+// For case (f)

+S_PC = (I_PC^2)*X_PC ; // Rating in VA

+S_PC1 = S_PC*10^-3 ; // Continuous kVA rating in kVA

+

+// For case (g)

+V_PC = I_PC * X_PC ; // continuous-voltage rating for PC in V

+

+// DISPLAY RESULTS

+disp("EXAMPLE : 9.16 :SOLUTION :-") ;

+printf("\n (a) Total zero-sequence susceptance per phase of system at 60 Hz , ΣX_C0 = %.4f Ω \n",X_C0) ;

+printf("\n     Total zero-sequence capacitance per phase of system at 60 Hz , ΣC_0 = %.4e F \n",TC_0) ;

+printf("\n (c) Continuous-current rating of the PC , I_PC = 3I_a0PC = %.4f A \n",abs(I_PC)) ;

+printf("\n (d) Required reactance value for the PC , X_PC = %.4f Ω \n",X_PC) ;

+printf("\n (e) Inductance value of the PC , L_PC = %.4f H \n",L_PC) ;

+printf("\n (f) Continuous kVA rating for the PC , S_PC = %.2f kVA \n",S_PC1) ;

+printf("\n (g) Continuous-voltage rating for PC , V_PC = %.2f V \n",V_PC) ;

diff --git a/1820/CH9/EX9.2/Example9_2.sce b/1820/CH9/EX9.2/Example9_2.sce
new file mode 100755
index 000000000..148df0c1d
--- /dev/null
+++ b/1820/CH9/EX9.2/Example9_2.sce
@@ -0,0 +1,39 @@
+// ELECTRIC POWER TRANSMISSION SYSTEM ENGINEERING ANALYSIS AND DESIGN

+// TURAN GONEN

+// CRC PRESS

+// SECOND EDITION

+

+// CHAPTER : 9 : SYMMETRICAL COMPONENTS AND FAULT ANALYSIS

+

+// EXAMPLE : 9.2 :

+clear ; clc ; close ; // Clear the work space and console

+

+// GIVEN DATA

+V_abc = [0 ; 50 ; -50] ; // Phase voltages of a 3-Φ system in V

+I_abc = [-5 ; 5*%i ; -5] ; // Phase current of a 3-Φ system in A

+

+// CALCULATIONS

+// For case (a)

+S_3ph = (V_abc)' * conj(I_abc) ; // 3-Φ complex power in VA

+

+// For case (b)

+a = 1*exp(%i*120*%pi/180) ; // By symmetrical components theory to 3-Φ system

+A = [1 1 1; 1 a^2 a ;1 a a^2] ;

+V_012 = inv(A) * (V_abc) ; // Sequence voltage matrices in V

+I_012 = inv(A) * (I_abc) ; // Sequence current matrices in A

+

+// For case (c)

+S_3ph1 = 3 * ([V_012(1,1) V_012(2,1) V_012(3,1)]) * (conj(I_012)) ; // Three-phase complex power in VA . Refer equ 9.34(a)

+

+// DISPLAY RESULTS

+disp("EXAMPLE : 9.2 : SOLUTION :-") ;

+printf("\n (a) Three-phase complex power using equ 9.30 , S_3-Φ = %.4f<%.f VA \n",abs(S_3ph) , atand(imag(S_3ph),real(S_3ph) )) ;

+printf("\n (b) Sequence Voltage matrices , [V_012] =  V \n") ;

+printf("\n     %.f<%.f ",abs(V_012(1,1)),atand( imag(V_012(1,1)),real(V_012(1,1)) )) ;

+printf("\n     %.4f<%.f ",abs(V_012(2,1)),atand( imag(V_012(2,1)),real(V_012(2,1)) )) ;

+printf("\n     %.4f<%.f ",abs(V_012(3,1)),atand( imag(V_012(3,1)),real(V_012(3,1)) )) ;

+printf("\n \n    Sequence current matrices , [I_012] = A \n") ;

+printf("\n     %.4f<%.1f ",abs(I_012(1,1)),atand( imag(I_012(1,1)),real(I_012(1,1)) )) ;

+printf("\n     %.4f<%.f ",abs(I_012(2,1)),atand( imag(I_012(2,1)),real(I_012(2,1)) )) ;

+printf("\n     %.4f<%.f ",abs(I_012(3,1)),atand( imag(I_012(3,1)),real(I_012(3,1)) )) ;

+printf("\n \n (c) Three-phase complex power using equ 9.34 , S_3-Φ = %.4f<%.f VA \n",abs(S_3ph1) , atand(imag(S_3ph1),real(S_3ph1) )) ;

diff --git a/1820/CH9/EX9.3/Example9_3.sce b/1820/CH9/EX9.3/Example9_3.sce
new file mode 100755
index 000000000..bd12b9d7e
--- /dev/null
+++ b/1820/CH9/EX9.3/Example9_3.sce
@@ -0,0 +1,44 @@
+// ELECTRIC POWER TRANSMISSION SYSTEM ENGINEERING ANALYSIS AND DESIGN

+// TURAN GONEN

+// CRC PRESS

+// SECOND EDITION

+

+// CHAPTER : 9 : SYMMETRICAL COMPONENTS AND FAULT ANALYSIS

+

+// EXAMPLE : 9.3 :

+clear ; clc ; close ; // Clear the work space and console

+

+// GIVEN DATA

+l = 40 ; // line length in miles

+// Conductor parameter from Table A.3

+r_a = 0.206 ; // Ohms per conductor per mile in Ω/mi

+r_b = r_a ; // r_a = r_b = r_c in Ω/mi

+D_s = 0.0311 ; // GMR in ft where D_s = D_sa = D_sb = D_sc

+D_ab = sqrt(2^2 + 8^2) ; // GMR in ft

+D_bc = sqrt(3^2 + 13^2) ; // GMR in ft

+D_ac = sqrt(5^2 + 11^2) ; // GMR in ft

+D_e = 2788.5 ; // GMR in ft since earth resistivity is zero

+r_e = 0.09528 ; // At 60 Hz in Ω/mi

+

+// CALCULATIONS

+// For case (a)

+Z_aa =[(r_a + r_e) + %i * 0.1213*log(D_e/D_s)]*l ; // Self impedance of line conductor in Ω

+Z_bb = Z_aa ;

+Z_cc = Z_bb ;

+Z_ab = [r_e + %i * 0.1213*log(D_e/D_ab)]*l ; // Mutual impedance in Ω

+Z_ba = Z_ab ;

+Z_bc = [r_e + %i * 0.1213*log(D_e/D_bc)]*l ;

+Z_cb = Z_bc ;

+Z_ac = [r_e + %i * 0.1213*log(D_e/D_ac)]*l ;

+Z_ca = Z_ac ;

+Z_abc = [Z_aa Z_ab Z_ac ; Z_ba Z_bb Z_bc ; Z_ca Z_cb Z_cc] ; // Line impedance matrix

+

+// For case (b)

+a = 1*exp(%i*120*%pi/180) ; // By symmetrical components theory to 3-Φ system

+A = [1 1 1; 1 a^2 a ;1 a a^2] ;

+Z_012 = inv(A) * Z_abc*A ; // Sequence impedance matrix

+

+// DISPLAY RESULTS

+disp("EXAMPLE : 9.3 : SOLUTION :-") ;

+printf("\n (a) Line impedance matrix , [Z_abc] = \n") ; disp(Z_abc) ;

+printf("\n (b) Sequence impedance matrix of line , [Z_012] = \n") ; disp(Z_012) ;

diff --git a/1820/CH9/EX9.4/Example9_4.sce b/1820/CH9/EX9.4/Example9_4.sce
new file mode 100755
index 000000000..8565dbd34
--- /dev/null
+++ b/1820/CH9/EX9.4/Example9_4.sce
@@ -0,0 +1,36 @@
+// ELECTRIC POWER TRANSMISSION SYSTEM ENGINEERING ANALYSIS AND DESIGN

+// TURAN GONEN

+// CRC PRESS

+// SECOND EDITION

+

+// CHAPTER : 9 : SYMMETRICAL COMPONENTS AND FAULT ANALYSIS

+

+// EXAMPLE : 9.4 :

+clear ; clc ; close ; // Clear the work space and console

+

+// GIVEN DATA

+l = 40 ; // line length in miles

+// Conductor parameter from Table A.3

+r_a = 0.206 ; // Ohms per conductor per mile in Ω/mi

+r_b = r_a ; // r_a = r_b = r_c in Ω/mi

+D_s = 0.0311 ; // GMR in ft where D_s = D_sa = D_sb = D_sc

+D_ab = sqrt(2^2 + 8^2) ; // GMR in ft

+D_bc = sqrt(3^2 + 13^2) ; // GMR in ft

+D_ac = sqrt(5^2 + 11^2) ; // GMR in ft

+D_e = 2788.5 ; // GMR in ft since earth resistivity is zero

+r_e = 0.09528 ; // At 60 Hz in Ω/mi

+

+// CALCULATIONS

+// For case (a)

+Z_s =[(r_a + r_e) + %i*0.1213*log(D_e/D_s)]*l ; // Self impedance of line conductor in Ω . From equ 9.49

+D_eq = (D_ab * D_bc * D_ac)^(1/3) ; // Equ GMR

+Z_m = [r_e + %i*0.1213*log(D_e/D_eq)]*l ; // From equ 9.50

+Z_abc = [Z_s Z_m Z_m ; Z_m Z_s Z_m ; Z_m Z_m Z_s] ; // Line impedance matrix

+

+// For case (b)

+Z_012 = [(Z_s+2*Z_m) 0 0 ; 0 (Z_s-Z_m) 0 ; 0 0 (Z_s-Z_m)] ; // Sequence impedance matrix . From equ 9.54

+

+// DISPLAY RESULTS

+disp("EXAMPLE : 9.4 : SOLUTION :-") ;

+printf("\n (a) Line impedance matrix when line is completely transposed , [Z_abc] = \n") ; disp(Z_abc) ;

+printf("\n (b) Sequence impedance matrix when line is completely transposed , [Z_012] = \n") ; disp(Z_012) ;

diff --git a/1820/CH9/EX9.5/Example9_5.sce b/1820/CH9/EX9.5/Example9_5.sce
new file mode 100755
index 000000000..15e6dc2c2
--- /dev/null
+++ b/1820/CH9/EX9.5/Example9_5.sce
@@ -0,0 +1,45 @@
+// ELECTRIC POWER TRANSMISSION SYSTEM ENGINEERING ANALYSIS AND DESIGN

+// TURAN GONEN

+// CRC PRESS

+// SECOND EDITION

+

+// CHAPTER : 9 : SYMMETRICAL COMPONENTS AND FAULT ANALYSIS

+

+// EXAMPLE : 9.5 :

+clear ; clc ; close ; // Clear the work space and console

+

+// GIVEN DATA

+Z_012 = [(19.6736 + 109.05044*%i) (0.5351182 + 0.4692097*%i) (- 0.5351182 + 0.4692097*%i) ; (- 0.5351182 + 0.4692097*%i) (8.24 + 28.471684*%i) (- 1.0702365 - 0.9384195*%i) ; (0.5351182 + 0.4692097*%i) (1.0702365 - 0.9384195*%i) (8.24 + 28.471684*%i)] ; // Line impedance matrix . result of exa 9.3

+Y_012 = inv(Z_012) ; // Sequence admittance of line

+

+// CALCULATIONS

+// For case (a)

+Y_01 = Y_012(1,2) ;

+Y_11 = Y_012(2,2) ;

+m_0 = Y_01/Y_11 ; // Per-unit unbalance for zero-sequence in pu from equ 9.67b

+m_0_per = m_0 * 100 ; // Per-unit unbalance for zero-sequence in percentage

+

+// For case (b)

+Z_01 = Z_012(1,2) ;

+Z_00 = Z_012(1,1) ;

+m_01 = -(Z_01/Z_00) ; // Per-unit unbalance for zero-sequence in pu from equ 9.67b

+m_01_per = m_01 * 100 ; // Per-unit unbalance for zero-sequence in percentage

+

+// For case (c)

+Y_21 = Y_012(3,2) ;

+Y_11 = Y_012(2,2) ;

+m_2 = (Y_21/Y_11) ; // Per-unit unbalance for zero-sequence in pu from equ 9.67b

+m_2_per = m_2 * 100 ; // Per-unit unbalance for zero-sequence in percentage

+

+// For case (d)

+Z_21 = Z_012(3,2) ;

+Z_22 = Z_012(3,3) ;

+m_21 = -(Z_21/Z_22) ; // Per-unit unbalance for zero-sequence in pu from equ 9.67b

+m_21_per = m_21 * 100 ; // Per-unit unbalance for zero-sequence in percentage

+

+// DISPLAY RESULTS

+disp("EXAMPLE : 9.5 : SOLUTION :-") ;

+printf("\n (a) Per-unit electromagnetic unbalance for zero-sequence , m_0 = %.2f<%.1f percent pu \n",abs(m_0_per),atand( imag(m_0_per),real(m_0_per) )) ;

+printf("\n (b) Approximate value of Per-unit electromagnetic unbalance for negative-sequence , m_0 = %.2f<%.1f percent pu \n",abs(m_01_per),atand( imag(m_01_per),real(m_01_per) )) ;

+printf("\n (c) Per-unit electromagnetic unbalance for negative-sequence , m_2 = %.2f<%.1f percent pu \n",abs(m_2_per),atand( imag(m_2_per),real(m_2_per) )) ;

+printf("\n (d) Approximate value of Per-unit electromagnetic unbalance for negative-sequence , m_2 = %.2f<%.1f percent pu \n",abs(m_21_per),atand( imag(m_21_per),real(m_21_per) )) ;

diff --git a/1820/CH9/EX9.6/Example9_6.sce b/1820/CH9/EX9.6/Example9_6.sce
new file mode 100755
index 000000000..846a349c3
--- /dev/null
+++ b/1820/CH9/EX9.6/Example9_6.sce
@@ -0,0 +1,56 @@
+

+// ELECTRIC POWER TRANSMISSION SYSTEM ENGINEERING ANALYSIS AND DESIGN

+// TURAN GONEN

+// CRC PRESS

+// SECOND EDITION

+

+// CHAPTER : 9 : SYMMETRICAL COMPONENTS AND FAULT ANALYSIS

+

+// EXAMPLE : 9.6 :

+clear ; clc ; close ; // Clear the work space and console

+

+// GIVEN DATA

+kv = 115 ; // Line voltage in kV

+

+// For case (a)

+h_11 = 90 ; // GMD b/w ground wires & their images

+r_a = 0.037667 ; // Radius in metre

+p_aa = 11.185 * log(h_11/r_a) ; // unit is F^(-1)m

+p_bb = p_aa ;

+p_cc = p_aa ;

+l_12 = sqrt(22 + (45 + 37)^2) ;

+D_12 = sqrt(2^2 + 8^2) ; // GMR in ft

+p_ab = 11.185*log(l_12/D_12) ; // unit is F^(-1)m

+p_ba = p_ab ;

+D_13 = sqrt(3^2 + 13^2) ; // GMR in ft

+l_13 = 94.08721051 ;

+p_ac = 11.185 * log(l_13/D_13) ; // unit is F^(-1)m

+p_ca = p_ac ;

+l_23 = 70.72279912 ;

+D_23 = sqrt(5^2 + 11^2) ; // GMR in ft

+p_bc = 11.185 * log(l_23/D_23) ; // unit is F^(-1)m

+p_cb = p_bc ;

+P_abc = [p_aa p_ab p_ac ; p_ba p_bb p_bc ; p_ca p_cb p_cc] ; // Matrix of potential coefficients

+

+// For case (b)

+C_abc = inv(P_abc) ; // Matrix of maxwells coefficients

+

+// For case (c)

+a = 1*exp(%i*120*%pi/180) ; // By symmetrical components theory to 3-Φ system

+A = [1 1 1; 1 a^2 a ;1 a a^2] ;

+C_012 = inv(A) * C_abc * A ; // Matrix of sequence capacitances

+

+// For case (d)

+C_01 = C_012(1,2) ;

+C_11 = C_012(2,2) ;

+C_21 = C_012(3,2) ;

+d_0 = C_01/C_11 ; // Zero-sequence electrostatic unbalances . Refer equ 9.115

+d_2 = -C_21/C_11 ; // Negative-sequence electrostatic unbalances . Refer equ 9.116

+

+// DISPLAY RESULTS

+disp("EXAMPLE : 9.6 : SOLUTION :-") ;

+printf("\n (a) Matrix of potential coefficients , [P_abc] = \n") ; disp(P_abc) ;

+printf("\n (b) Matrix of maxwells coefficients , [C_abc] = \n") ; disp(C_abc) ;

+printf("\n (c) Matrix of sequence capacitances , [C_012] = \n") ; disp(C_012) ;

+printf("\n (d) Zero-sequence electrostatic unbalances , d_0 = %.4f<%.1f \n",abs(d_0),atand( imag(d_0),real(d_0) )) ;

+printf("\n     Negative-sequence electrostatic unbalances , d_2 = %.4f<%.1f \n",abs(d_2),atand( imag(d_2),real(d_2) )) ;

diff --git a/1820/CH9/EX9.9/Example9_9.sce b/1820/CH9/EX9.9/Example9_9.sce
new file mode 100755
index 000000000..b712f6b58
--- /dev/null
+++ b/1820/CH9/EX9.9/Example9_9.sce
@@ -0,0 +1,52 @@
+// ELECTRIC POWER TRANSMISSION SYSTEM ENGINEERING ANALYSIS AND DESIGN

+// TURAN GONEN

+// CRC PRESS

+// SECOND EDITION

+

+// CHAPTER : 9 : SYMMETRICAL COMPONENTS AND FAULT ANALYSIS

+

+// EXAMPLE : 9.9 :

+clear ; clc ; close ; // Clear the work space and console

+

+// GIVEN DATA

+kv = 230 ; // Line voltage in kV

+Z_0 = 0.56 * %i ; // impedance in Ω

+Z_1 = 0.2618 * %i ; // Impedance in Ω

+Z_2 = 0.3619 * %i ; // Impedance in Ω

+z_f = 5 + 0*%i ; // fault impedance in Ω

+v = 1 * exp(%i*0*%pi/180) ;

+

+// CALCULATIONS

+// For case (a)

+Z_B = kv^2/200 ; // Imedance base on 230 kV line

+Z_f = z_f/Z_B ; // fault impedance in pu Ω

+I_a0 = v/(Z_0 + Z_1 + Z_2 + 3*Z_f) ; // Sequence currents in pu A

+I_a1 = I_a0 ;

+I_a2 = I_a0 ;

+a = 1 * exp(%i*120*%pi/180) ; // By symmetrical components theory to 3-Φ system

+A = [1 1 1; 1 a^2 a ;1 a a^2] ;

+I_f = A * [I_a0 ; I_a1 ; I_a2] ; // Phase currents in pu A

+

+// For case (b)

+V_a = [0 ; v ; 0] - [Z_0 0 0 ; 0 Z_1 0 ; 0 0 Z_2]*[I_a0 ; I_a1 ; I_a2] ; // Sequence voltage in pu V

+V_f = A*V_a ; // Phase voltage in pu V

+

+// For case (c)

+V_abf = V_f(1,1) - V_f(2,1) ; // Line-to-line voltages at fault points in pu V

+V_bcf = V_f(2,1) - V_f(3,1) ; // Line-to-line voltages at fault points in pu V

+V_caf = V_f(3,1) - V_f(1,1) ; // Line-to-line voltages at fault points in pu V

+

+// DISPLAY RESULTS

+disp("EXAMPLE : 9.9 : SOLUTION :-") ;

+printf("\n (b) Sequence currents , I_a0 = I_a1 = I_a2 = %.4f<%.1f pu A \n",abs(I_a0),atand(imag(I_a0),real(I_a0) )) ;

+printf("\n  Phase currents in pu A , [I_af ; I_bf ; I_cf] = pu A \n") ;

+printf("\n    %.4f<%.1f ",abs(I_f),atand(imag(I_f),real(I_f) )) ;

+printf("\n \n (c) Sequence voltages are , [V_a0 ; V_a1 ; V_a2 ] = pu V \n") ;

+printf("\n    %.4f<%.1f ",abs(V_a),atand(imag(V_a),real(V_a) )) ;

+printf("\n \n    Phase voltages are , [V_af ; V_bf ; V_cf ] = pu V \n") ;

+printf("\n    %.4f<%.1f ",abs(V_f),atand(imag(V_f),real(V_f) )) ;

+printf("\n \n (d) Line-to-line voltages at fault points are , V_abf = %.4f<%.1f pu V \n",abs(V_abf),atand(imag(V_abf),real(V_abf) )) ;

+printf("\n     Line-to-line voltages at fault points are , V_abf = %.4f<%.1f pu V \n",abs(V_bcf),atand(imag(V_bcf),real(V_bcf) )) ;

+printf("\n     Line-to-line voltages at fault points are , V_caf = %.4f<%.1f pu V \n",abs(V_caf),atand(imag(V_caf),real(V_caf) )) ;

+

+printf("\n NOTE : ERROR : Calclation mistake in textbook from case(c) onwards \n") ;