blob: c1c78023e184208f362e0ce91acda84ae1a4a419 (
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
25
26
27
28
29
30
31
32
33
34
35
36
|
//Book - Power system: Analysisi & Design 5th Edition
//Authors - J. Duncan Glover, Mulukutla S. Sarma, and Thomas J.Overbye
//Chapter-3 ;Example 3.2
//Scilab Version - 6.0.0; OS - Windows
clc;
clear;
Srated=20 //rated input power in kVA
E1rated=480 //Rated voltage across winding 1 in Volts
E2rated=120 //Rated voltage across winding 2 in Volts
F=60 //frequency in Hertz
Sl=15 //Load power in kVA
pf = 0.8 //power factor lagging
E2=118 //Load voltage in Volts
V1s=35 //Short circuit voltage in Volts
P1=300 //Short circuit power in Watts
I2=12 //Open circuit Winding 2 current in Amps
P2=200 //Open circuit power in Watts
I1rated=(Srated*1000)/E1rated //Rated current for winding 1
Req1=P1/(I1rated)^2 //Equivalent resistance of winding 1 in Ohms
Zeqm=abs(V1s/I1rated) //Magnitude of equivalent impedance of winding 1 in Ohms
Xeq1=sqrt(Zeqm^2-Req1^2) //Equivalent reactance of winding 1 in Ohms
Zeq1=Req1+%i*Xeq1 //Equivalent impedance of winding 1 in Ohms
V1o=E1rated //Since winding 1 open circuit voltage is equal to winding 1 rated volgage
Gc=P2/V1o^2
Ymm=abs((E2rated/E1rated)*I2/V1o)
Bm=sqrt(Ymm^2-Gc^2)
Ym=Gc-%i*Bm //Shunt admittance in Siemens
printf('The rated current for winding 1 is %.4f Ampere\n',I1rated);
printf('The Equivalent resistance of winding 1 is %.4f Ohms\n',Req1);
printf('The Equivalent reactance of winding 1 is %.4f Ohms\n',Xeq1);
printf('The magnitude of Equivalent impedance of winding 1 is %.4f Ohms and its angle is %.4f degree\n',abs(Zeq1),atand(imag(Zeq1),real(Zeq1)));
printf('The magnitude of Shunt admittance is %.4f Siemens and its angle is %.4f degree \n',abs(Ym),atand(imag(Ym),real(Ym)));
|
