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

---
 3872/CH12/EX12.4/Ex12_4.JPG | Bin 0 -> 34953 bytes
 3872/CH12/EX12.4/Ex12_4.sce |  31 +++++++++++++++++++++++++++++++
 2 files changed, 31 insertions(+)
 create mode 100644 3872/CH12/EX12.4/Ex12_4.JPG
 create mode 100644 3872/CH12/EX12.4/Ex12_4.sce

(limited to '3872/CH12/EX12.4')

diff --git a/3872/CH12/EX12.4/Ex12_4.JPG b/3872/CH12/EX12.4/Ex12_4.JPG
new file mode 100644
index 000000000..d4245e76a
Binary files /dev/null and b/3872/CH12/EX12.4/Ex12_4.JPG differ
diff --git a/3872/CH12/EX12.4/Ex12_4.sce b/3872/CH12/EX12.4/Ex12_4.sce
new file mode 100644
index 000000000..602114e9e
--- /dev/null
+++ b/3872/CH12/EX12.4/Ex12_4.sce
@@ -0,0 +1,31 @@
+//Book - Power System: Analysis & Design 5th Edition
+//Authors - J. Duncan Glover, Mulukutla S. Sarma, and Thomas J. Overbye
+//Chapter - 12 ; Example 12.4
+//Scilab Version - 6.0.0 ; OS - Windows
+
+clc;
+clear;
+funcprot(0);
+f=60;                          //Frequency in Hertz
+G=[1000 750 500];              //Rating of unit 1,2 &3 respectively in MVA
+R=0.05;                        //Regulation constant of each unit in pu
+delP=200;                      //Load increment in MW
+SBnew=1000;                    //New MVA base of the entire system
+
+Rnew=R*(SBnew./G);             //Regulation of each generators with common base
+beta=sum(1 ./Rnew);             //area frequency response characteristic, beta
+
+printf('The area frequency response characteristic beta is %.2f per unit\n',beta)
+
+delPpu=delP/SBnew;             //Load increment in per unit
+delFpu=(-1/beta)*delPpu        //Frequency drop in per unit
+delF=delFpu*f;                 //Frequency drop in Hertz
+
+printf('The steady-state frequency drop is %.4f Hz\n',abs(delF))
+
+delPm=delFpu*(-1 ./Rnew);
+delPmact=SBnew*delPm;
+
+printf('The increase in turbine mechanical power output of unit1=%.4f pu = %.4f MW\n',delPm(1),delPmact(1))
+printf('The increase in turbine mechanical power output of unit2=%.4f pu = %.4f MW\n',delPm(2),delPmact(2))
+printf('The increase in turbine mechanical power output of unit3=%.4f pu = %.4f MW',delPm(3),delPmact(3))
-- 
cgit