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.1/Ex12_1.JPG | Bin 0 -> 16531 bytes 3872/CH12/EX12.1/Ex12_1.sce | 30 +++++++++++++++++++++++++++++ 3872/CH12/EX12.2/Ex12_2.JPG | Bin 0 -> 21849 bytes 3872/CH12/EX12.2/Ex12_2.sce | 35 +++++++++++++++++++++++++++++++++ 3872/CH12/EX12.3/Ex12_3.JPG | Bin 0 -> 14455 bytes 3872/CH12/EX12.3/Ex12_3.sce | 19 ++++++++++++++++++ 3872/CH12/EX12.4/Ex12_4.JPG | Bin 0 -> 34953 bytes 3872/CH12/EX12.4/Ex12_4.sce | 31 +++++++++++++++++++++++++++++ 3872/CH12/EX12.5/Ex12_5.JPG | Bin 0 -> 31060 bytes 3872/CH12/EX12.5/Ex12_5.sce | 40 ++++++++++++++++++++++++++++++++++++++ 3872/CH12/EX12.6/Ex12_6.JPG | Bin 0 -> 44869 bytes 3872/CH12/EX12.6/Ex12_6.sce | 24 +++++++++++++++++++++++ 3872/CH12/EX12.7/Ex12_7.JPG | Bin 0 -> 40777 bytes 3872/CH12/EX12.7/Ex12_7.sce | 46 ++++++++++++++++++++++++++++++++++++++++++++ 3872/CH12/EX12.9/Ex12_9.JPG | Bin 0 -> 21453 bytes 3872/CH12/EX12.9/Ex12_9.sce | 30 +++++++++++++++++++++++++++++ 16 files changed, 255 insertions(+) create mode 100644 3872/CH12/EX12.1/Ex12_1.JPG create mode 100644 3872/CH12/EX12.1/Ex12_1.sce create mode 100644 3872/CH12/EX12.2/Ex12_2.JPG create mode 100644 3872/CH12/EX12.2/Ex12_2.sce create mode 100644 3872/CH12/EX12.3/Ex12_3.JPG create mode 100644 3872/CH12/EX12.3/Ex12_3.sce create mode 100644 3872/CH12/EX12.4/Ex12_4.JPG create mode 100644 3872/CH12/EX12.4/Ex12_4.sce create mode 100644 3872/CH12/EX12.5/Ex12_5.JPG create mode 100644 3872/CH12/EX12.5/Ex12_5.sce create mode 100644 3872/CH12/EX12.6/Ex12_6.JPG create mode 100644 3872/CH12/EX12.6/Ex12_6.sce create mode 100644 3872/CH12/EX12.7/Ex12_7.JPG create mode 100644 3872/CH12/EX12.7/Ex12_7.sce create mode 100644 3872/CH12/EX12.9/Ex12_9.JPG create mode 100644 3872/CH12/EX12.9/Ex12_9.sce (limited to '3872/CH12') diff --git a/3872/CH12/EX12.1/Ex12_1.JPG b/3872/CH12/EX12.1/Ex12_1.JPG new file mode 100644 index 000000000..407efe386 Binary files /dev/null and b/3872/CH12/EX12.1/Ex12_1.JPG differ diff --git a/3872/CH12/EX12.1/Ex12_1.sce b/3872/CH12/EX12.1/Ex12_1.sce new file mode 100644 index 000000000..2240ec1c0 --- /dev/null +++ b/3872/CH12/EX12.1/Ex12_1.sce @@ -0,0 +1,30 @@ +//Book - Power System: Analysis & Design 5th Edition +//Authors - J. Duncan Glover, Mulukutla S. Sarma, and Thomas J. Overbye +//Chapter - 12 ; Example 12.1 +//Scilab Version - 6.0.0 ; OS - Windows + +clc; +clear; + +Tr=0; +Ka=100; +Ta=0.05; +Vrmax=5; +Vrmin=-5; +Ke=1; +Te=0.26; +Kf=0.01; +Tf=1; + +Efd=2.9135; //Value taken from Example 11.10 +Vt=1.0946; //Value taken from Example 11.10 + +Vr=Ke*Efd; //Initial value of Vr +Vf=0; //Initial value of vf +Vref=(Vr/Ka)+Vt+Vf; //Initial value of Vref + +printf('The initial value of Vr is %.4f\n',Vr) +printf('The initial value of Vf is %.4f\n',Vf) +printf('The initial value of Vref is %.4f\n',Vref) + +//Section 'b' of this problem cannot be simulated using current version of Scilab diff --git a/3872/CH12/EX12.2/Ex12_2.JPG b/3872/CH12/EX12.2/Ex12_2.JPG new file mode 100644 index 000000000..b1e885d20 Binary files /dev/null and b/3872/CH12/EX12.2/Ex12_2.JPG differ diff --git a/3872/CH12/EX12.2/Ex12_2.sce b/3872/CH12/EX12.2/Ex12_2.sce new file mode 100644 index 000000000..edf79be4b --- /dev/null +++ b/3872/CH12/EX12.2/Ex12_2.sce @@ -0,0 +1,35 @@ +//Book - Power System: Analysis & Design 5th Edition +//Authors - J. Duncan Glover, Mulukutla S. Sarma, and Thomas J. Overbye +//Chapter - 12 ; Example 12.2 +//Scilab Version - 6.0.0 ; OS - Windows + +clc; +clear; + +KQi=0.4; +KVi=40; +XIqmax=1.45; +XIqmin=0.5; +Vmax=1.1; +Vmin=0.9; + +Tr=0; +Ka=100; +Ta=0.05; +Vrmax=5; +Vrmin=-5; +Ke=1; +Te=0.26; +Kf=0.01; +Tf=1; +vt=0.5; +Vf=0; //Initial value of vf +Vref=1.0239; //Initial value of Vref from Example 11.12 +Isorq=-XIqmax/0.8; // Reactive component of Isorc +Qcmd=0.22; //Obtained from Example 11.12 +Qnet=(vt)*abs(Isorq)-(vt)^2/0.8; // Net reactive power injection in pu + +printf('The initial value of reference voltage is %.4f pu\n',Vref) +printf('The initial value of reactive power Qcmd = %.4f pu = %.4f Mvar\n',Qcmd,Qcmd*100) +printf('The maxximum net reactive power Qnet = %.4f pu = %.4f Mvar\n',Qnet,Qnet*100) + diff --git a/3872/CH12/EX12.3/Ex12_3.JPG b/3872/CH12/EX12.3/Ex12_3.JPG new file mode 100644 index 000000000..310cac9cb Binary files /dev/null and b/3872/CH12/EX12.3/Ex12_3.JPG differ diff --git a/3872/CH12/EX12.3/Ex12_3.sce b/3872/CH12/EX12.3/Ex12_3.sce new file mode 100644 index 000000000..7b71197ab --- /dev/null +++ b/3872/CH12/EX12.3/Ex12_3.sce @@ -0,0 +1,19 @@ +//Book - Power System: Analysis & Design 5th Edition +//Authors - J. Duncan Glover, Mulukutla S. Sarma, and Thomas J. Overbye +//Chapter - 12 ; Example 12.3 +//Scilab Version - 6.0.0 ; OS - Windows + +clc; +clear; + +M=500; //MVA rating of the generator +f=60; //requency in Hertz +R=0.05; //Regulation constant in pu +delF=0.01; //Increase in frequency in Hertz + +delFpu=delF/f; //Frequency increase in pu +delPref=0; //Since fixed reference power setting is assumed +delPmpu=delPref-(1/R)*delFpu //Change in mechanical power in pu +delPm=delPmpu*M; //Actual value of mechanical power in MW + +printf('The turbine mechanical power output decreases by %.3f MW.',abs(delPm)) 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)) diff --git a/3872/CH12/EX12.5/Ex12_5.JPG b/3872/CH12/EX12.5/Ex12_5.JPG new file mode 100644 index 000000000..bdd340193 Binary files /dev/null and b/3872/CH12/EX12.5/Ex12_5.JPG differ diff --git a/3872/CH12/EX12.5/Ex12_5.sce b/3872/CH12/EX12.5/Ex12_5.sce new file mode 100644 index 000000000..054733001 --- /dev/null +++ b/3872/CH12/EX12.5/Ex12_5.sce @@ -0,0 +1,40 @@ +//Book - Power System: Analysis & Design 5th Edition +//Authors - J. Duncan Glover, Mulukutla S. Sarma, and Thomas J. Overbye +//Chapter - 12 ; Example 12.5 +//Scilab Version - 6.0.0 ; OS - Windows + +clc; +clear; + +f=60; //Frequency of the system in Hertz +G1=2000; //Total generation of area 1 in MW +G2=4000; //Total generation of area 2 in MW +beta1=700; //Area frequency response characteristic of area 1 in MW/Hz +beta2=1400; //Area frequency response characteristic of area 2 in MW/Hz +delPD1=100; //Load increment of area 1 in MW +delPD2=0; + +//WITHOUT LFC +delF=(delPD1+delPD2)/(-(beta1+beta2)); //Frequency Change in Hertz +delPm1=-beta1*delF; //Change in power of area 1 +delPm2=-beta2*delF; //Change in power of area 2 +delPtie1=-delPm2 //Tie line power flow from area 1 to 2 +delPtie2=delPm2 //Tie line power flow from area 2 to 1 + +disp('RESULTS WITHOUT LFC') +printf('\nThe steady state frequency error is %.4f Hz',delF) +printf('\nThe tie-line power flow from area 1 is %.4f MW',delPtie1) +printf('\nThe tie-line power flow from area 2 is %.4f MW\n',delPtie2) + +//WITH LFC + +delFl=0/(beta1+beta2); //Frequency Change in Hertz (as ACE1+ACE2=0) +delPm1=-beta1*delFl; //Change in power of area 1 +delPm2=-beta2*delFl; //Change in power of area 2 +delPtie1=-delPm2 //Tie line power flow from area 1 to 2 +delPtie2=delPm2 //Tie line power flow from area 2 to 1 + +disp('RESULTS WITH LFC') +printf('\nThe steady state frequency error is %.4f Hz',delFl) +printf('\nThe tie-line power flow from area 1 is %.4f MW',delPtie1) +printf('\nThe tie-line power flow from area 2 is %.4f MW',delPtie2) diff --git a/3872/CH12/EX12.6/Ex12_6.JPG b/3872/CH12/EX12.6/Ex12_6.JPG new file mode 100644 index 000000000..7956950a9 Binary files /dev/null and b/3872/CH12/EX12.6/Ex12_6.JPG differ diff --git a/3872/CH12/EX12.6/Ex12_6.sce b/3872/CH12/EX12.6/Ex12_6.sce new file mode 100644 index 000000000..19fe7085a --- /dev/null +++ b/3872/CH12/EX12.6/Ex12_6.sce @@ -0,0 +1,24 @@ +//Book - Power System: Analysis & Design 5th Edition +//Authors - J. Duncan Glover, Mulukutla S. Sarma, and Thomas J. Overbye +//Chapter - 12 ; Example 12.6 +//Scilab Version - 6.0.0 ; OS - Windows + +clc; +clear; + +C1=[8e-3 10 0] //Coefficients of cost equation for unit 1 +C2=[9e-3 8 0] //Coefficients of cost equation for unit 2 + +dC1=[2*C1(1) C1(2)] //Coefficients of incremental cost equation for unit 1 +dC2=[2*C2(1) C2(2)] //Coefficients of incremental cost equation for unit 2 +result=[]; +for PT=500:100:1500 + P1=(dC2(1)*PT+(dC2(2)-dC1(2)))/(dC2(1)+dC1(1)); + P2=PT-P1; + dC1value=dC1(1)*P1+dC1(2); + dC2value=dC2(1)*P2+dC2(2); + CT=C1(1)*P1^2+C1(2)*P1+C1(3)+C2(1)*P2^2+C2(2)*P2+C2(3); + result=[result;PT P1 P2 dC1value CT] +end + +disp(result,' PT(MW) P1(MW) P2(MW) dC1=dC2($/MWhr) CT($/hr)'); diff --git a/3872/CH12/EX12.7/Ex12_7.JPG b/3872/CH12/EX12.7/Ex12_7.JPG new file mode 100644 index 000000000..1ecd121a2 Binary files /dev/null and b/3872/CH12/EX12.7/Ex12_7.JPG differ diff --git a/3872/CH12/EX12.7/Ex12_7.sce b/3872/CH12/EX12.7/Ex12_7.sce new file mode 100644 index 000000000..6bf46aba1 --- /dev/null +++ b/3872/CH12/EX12.7/Ex12_7.sce @@ -0,0 +1,46 @@ +//Book - Power System: Analysis & Design 5th Edition +//Authors - J. Duncan Glover, Mulukutla S. Sarma, and Thomas J. Overbye +//Chapter - 12 ; Example 12.7 +//Scilab Version - 6.0.0 ; OS - Windows + +clc; +clear; + +C1=[8e-3 10 0] //Coefficients of cost equation for unit 1 +C2=[9e-3 8 0] //Coefficients of cost equation for unit 2 + +dC1=[2*C1(1) C1(2)] //Coefficients of incremental cost equation for unit 1 +dC2=[2*C2(1) C2(2)] //Coefficients of incremental cost equation for unit 2 + +P1lim=[100 600]; //Lower and upper generation limit for unit 1 +P2lim=[400 1000]; //Lower and upper generation limit for unit 2 + +result=[]; +for PT=[500 600 700 725 800 900 1000 1100 1200 1244 1300 1400 1500] + P1=(dC2(1)*PT+(dC2(2)-dC1(2)))/(dC2(1)+dC1(1)); + P2=PT-P1; + dC1value=dC1(1)*P1+dC1(2); + dC2value=dC2(1)*P2+dC2(2); + + if P1P1lim(2) //Checking for limits of P1 + if P1P2lim(2) //Checking for limits of P2 + if P2