summaryrefslogtreecommitdiff
path: root/3872/CH12/EX12.4/Ex12_4.sce
diff options
context:
space:
mode:
authorprashantsinalkar2017-10-10 12:27:19 +0530
committerprashantsinalkar2017-10-10 12:27:19 +0530
commit7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 (patch)
treedbb9e3ddb5fc829e7c5c7e6be99b2c4ba356132c /3872/CH12/EX12.4/Ex12_4.sce
parentb1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (diff)
downloadScilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.gz
Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.bz2
Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.zip
initial commit / add all books
Diffstat (limited to '3872/CH12/EX12.4/Ex12_4.sce')
-rw-r--r--3872/CH12/EX12.4/Ex12_4.sce31
1 files changed, 31 insertions, 0 deletions
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))