Activity 1 Admittance Matrix Formulation Lab

Download as pdf or txt
Download as pdf or txt
You are on page 1of 4

POWER SYSTEM ANALYSIS II

Activity No. 1
Admittance Matrix Formulation
I.

DELIVERABLES
1. Softcopy of source code and required output text files (.txt).
2. Hardcopy of documentation. Submit in A4 paper, IEEE double-column format.
3. Place all files in an encrypted compressed folder (RAR) and upload it at our class e-group with the
following filename format: SURNAME_GROUP#_ACT# (e.g. CRUZ_GROUP1_ACT1.rar)

II.

OBJECTIVES
1. To be able to understand common data format used in solving load flow problems.
2. To be able to use MATLAB or other programming tool/language to read and parse data.
3. To be able to develop a program that will formulate Bus Admittance Matrix.

III.

DISCUSSION

Formation of Bus Admittance Matrix

Most power system networks are analyzed by first forming the admittance matrix. The admittance
matrix is based upon Kirchhoff's current law (KCL), and it is easily formed and very sparse.
Consider the three-bus network shown in Figure that has five branch impedances and one current
source.

ZA

ZC

ZD

ZB

ZE

3
I3

Figure 1. Three-Bus Network


Applying KCL at the three independent nodes yields the following equations for the bus voltages (with
respect to ground):
At bus 1,

V1 V1 V2

0 ,
ZE
ZA

At bus 2,

V2 V2 V1 V2 V3

0 ,
ZB
ZA
ZC

At bus 3,

V3 V3 V2

I3 .
ZD
ZC

Collecting terms and writing the equations in matrix form yields

1
1

ZE Z A
1

ZA

1
ZA
1
1
1

Z A Z B ZC
1

ZC

V1 0
1

V2 0 ,
ZC

1
1 V3 I 3

Z C Z D
0

or in matrix form,

YV I ,
where Y is the admittance matrix, V is a vector of bus voltages (with respect to ground), and I is a vector
of current injections.
Voltage sources, if present, can be converted to current sources using the usual network rules. If a bus
has a zero-impedance voltage source attached to it, then the bus voltage is already known, and the
dimension of the problem is reduced by one.
A simple observation of the structure of the above admittance matrix leads to the following rule for
building Y:
1. The diagonal terms of Y contain the sum of all branch admittances connected directly to the
corresponding bus.

2. The off-diagonal elements of Y contain the negative sum of all branch admittances connected
directly between the corresponding busses.

Modeling of Tap changers

Figure 1. Tap changer connection


This is valid for real values of a. Taking the y-bus formation, break the diagonal elements into two
components:
the off-diagonal element represent the impedance across the two buses.
the remainder form the shunt element.

Figure 2. Equivalent model for Tap changers


2

Sample input data:


Table 1. Bus data
|Bus | Type | Vsp | theta | PGi |
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16

1
2
2
3
3
2
3
2
2
3
3
2
3
3
3
3

1.040
1.010
0.985
1.000
1.000
0.980
1.000
1.005
0.980
1.000
1.000
1.015
1.000
1.000
1.000
1.000

0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0

0.0
3.0
41.0
0.0
13.0
75.0
0.0
150.0
121.0
5.0
0.0
377.0
18.0
10.5
22.0
43.0

QGi | PLi |

QLi |

Qmin |

0.0
0.0
88.0
0.0
21.0 40.0
0.0
0.0
4.0
0.0
2.0
0.0
0.0
0.0
22.0 450.0
26.0
0.0
2.0
0.0
0.0
0.0
24.0 310.0
2.3
0.0
5.3
0.0
5.0
0.0
3.0
0.0

0.0
-0.8
-1.0
0.0
0.0
0.8
0.0
62.1
2.2
0.0
0.0
128.5
0.0
0.0
0.0
0.0

0.0
50.0
60.0
0.0
0.0
25.0
0.0
200.0
9.0
0.0
0.0
155.0
0.0
0.0
0.0
0.0

Qmax |
0.0
-17.0
-10.0
0.0
0.0
-8.0
0.0
-140.0
-3.0
0.0
0.0
-150.0
0.0
0.0
0.0
0.0

(succeeding data are intentionally deleted)

Table 2. Line data


|
|

From | To
Bus | Bus
1
2
3
4
4
6
6
8
9
9
9
9
13
13
1
1

2
3
4
5
6
7
8
9
10
11
12
13
14
15
15
16

|
|

R
pu
0.0083
0.0298
0.0112
0.0625
0.0430
0.0200
0.0339
0.0099
0.0369
0.0258
0.0648
0.0481
0.0132
0.0269
0.0178
0.0454

|
|

X
pu
0.0280
0.0850
0.0366
0.1320
0.1480
0.1020
0.1730
0.0505
0.1679
0.0848
0.2950
0.1580
0.0434
0.0869
0.0910
0.2060

|
|

B/2
pu
0.0645
0.0409
0.0190
0.0129
0.0174
0.0138
0.0235
0.0274
0.0220
0.0109
0.0386
0.0203
0.0055
0.0115
0.0494
0.0273

| X'mer |
| TAP (a) |
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1

(succeeding data are intentionally deleted)

IV.

PROBLEM SPECIFICATION
1. The input text files should follow the format shown below.
2. The high-level programming language Matlab can be used but not limited to other programming
language such as C/C++, Java, Visual Basic/C++, etc..
3. The documentation should include but not limited to the following sections.
i. Abstract brief summary of the problem and your program
ii. Implementation include a short description of each program module/part
iii. Testing emphasize all assumptions you made in your program.
iv. Results and Discussion
v. Conclusions answer the problem objectives
vi. On the first page/cover page of your documentation, write the boxed statement below (or any
statement of the same effect) and sign.
3

We certify upon our honor that we solved this activity independently, without asking help
from anybody. We may have, however, consulted published references and reading
materials, and we hereby acknowledge (write proper bibliography) them accordingly. It
took us _______ hours to finish this activity.

___________________
Sign over printed name

V.

___________________
Sign over printed name

DISCUSSION
1. What programming language/tool did you use in this activity and why?
2. Where you able to meet all objectives of this activity? State all the challenges you have encountered
in this activity and discuss how you come up with the solutions.

You might also like