blob: 7ccaed87a40c73bfc41796155532b1a980c12139 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
|
//Book - Power system: Analysisi & Design 5th Edition
//Authors - J. Duncan Glover, Mulukutla S. Sarma, and Thomas J.Overbye
//Chapter-11 ;Example 11.3
//Scilab Version - 6.0.0; OS - Windows
clc;
clear;
P=1.0 //Infinite bus received real power in per unit
Vbus=1.0 //Infinite bus voltage in per unit
pf=0.95 //Lagging power factor
Xdt=0.30
XTR=0.10
X12=0.20
X13=0.10
X23=0.20
Xeq=Xdt+XTR+(X12*(X13+X23))/(X12+(X13+X23)); //The equialent reactance between the machine internal voltage and infinite bus in per unit
theta=acos(pf);
I=(P/(Vbus*pf))*exp(-%i*theta); //Current into the infinite bus in per unit
Ei=Vbus+(%i*Xeq)*I; //The machine internal voltage in per unit
printf('The magnitude of he machine internal voltage in per unit is %.4f pu and its angle is %.4f degrees',abs(Ei),atand(imag(Ei),real(Ei)));
|
