From 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 Mon Sep 17 00:00:00 2001
From: prashantsinalkar
Date: Tue, 10 Oct 2017 12:27:19 +0530
Subject: initial commit / add all books

---
 3472/CH17/EX17.28/Example17_28.sce | 28 ++++++++++++++++++++++++++++
 1 file changed, 28 insertions(+)
 create mode 100644 3472/CH17/EX17.28/Example17_28.sce

(limited to '3472/CH17/EX17.28/Example17_28.sce')

diff --git a/3472/CH17/EX17.28/Example17_28.sce b/3472/CH17/EX17.28/Example17_28.sce
new file mode 100644
index 000000000..e318b6bb1
--- /dev/null
+++ b/3472/CH17/EX17.28/Example17_28.sce
@@ -0,0 +1,28 @@
+// A Texbook on POWER SYSTEM ENGINEERING
+// A.Chakrabarti, M.L.Soni, P.V.Gupta, U.S.Bhatnagar
+// DHANPAT RAI & Co.
+// SECOND EDITION 
+
+// PART II : TRANSMISSION AND DISTRIBUTION
+// CHAPTER 10: POWER SYSTEM STABILITY
+
+// EXAMPLE : 10.28 :
+// Page number 310
+clear ; clc ; close ; // Clear the work space and console
+
+// Given data
+sin_delta_0 = 0.45     // Supplying percent of peak power capacity before fault
+x = 4.0                // Reactance under fault increased
+gamma_2 = 0.7          // Peak power delivered after fault clearance
+
+// Calculations
+delta_0 = asin(sin_delta_0)                         // δ_0(radians)
+delta_0_degree = delta_0*180/%pi                    // δ_0(°)
+gamma_1 = 1.0/x                                     // γ_1
+delta_m = %pi-asin(sin_delta_0/(gamma_2))           // δ_m(radians)
+delta_m_degree = delta_m*180/%pi                    // δ_m(°)
+delta_c = acosd(1.0/(gamma_2-gamma_1)*((delta_m-delta_0)*sin(delta_0)+(gamma_2*cos(delta_m)-gamma_1*cos(delta_0)))) // Clearing angle(°)
+
+// Results
+disp("PART II - EXAMPLE : 10.28 : SOLUTION :-")
+printf("\nCritical clearing angle, δ_c = %.f° ", delta_c)
-- 
cgit