Scribd
Upload
41 views3 pages

Admittance Matrix Formulation

Bus Admittance Matrix calculate the amount of current present in the bus, This Paper contains a code that asks for a text file input. Text file input contains line data. This line data is compose of from bus to bus, resistance in per unit, line charge admittance, and transformer tap changer. Then the code processes the input to generate a bus admittance matrix Y-bus.

Uploaded by

Juls Bacton
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
41 views3 pages

Admittance Matrix Formulation

Bus Admittance Matrix calculate the amount of current present in the bus, This Paper contains a code that asks for a text file input. Text file input contains line data. This line data is compose of from bus to bus, resistance in per unit, line charge admittance, and transformer tap changer. Then the code processes the input to generate a bus admittance matrix Y-bus.

Uploaded by

Juls Bacton
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 3

Admittance Matrix Formulation

Julius John L. Bacton Nathannaell Reiben M. Terio


DEET DEET
Mindanao State University -Iligan Institute of Technology Mindanao State University -Iligan Institute of Technology
Iligan City, Philippines Iligan City, Philippines
[email protected] [email protected]

Abstract—This paper contains a code that asks for a text frombus = Ldata(:,1);
file input. Text file input contains line data. This line data tobus = Ldata(:,2);
is compose of from bus, to bus, resistance in per unit, line r = Ldata(:,3);
charging admittance, and transformer tap changer. Then the x = Ldata(:,4);
code processes the input to generate a bus admittance matrix Bover2 = Ldata(:,5);
Y-bus. The output will then display in the command window. tap = Ldata(:,6);

Index Terms—Admittance, Matrix, Reactance y = zeros(a,1);


z = zeros(a,1);
I. I NTRODUCTION b = 1i * Bover2;
In a power system, the nodal admittance of the various buses for n = 1:1:a
is presented by Bus Admittance Matrix. Each bus is connected z(n) = r(n) + 1i*x(n);
to the various other buses with the help of the transmission y(n) = 1/z(n);
lines and admittance matrix is used to examine the data that end
is required the study of the load flow in a power system. The nbus = max(max(frombus),max(tobus));
admittance and the topology of the network are also explained Y = zeros(nbus,nbus);
with help of the Bus Admittance Matrix. Bus Admittance
Matrix can also calculate the amount of current present in the %off diagonal element
bus. Coming to the formation of the Bus Admittance Matrix for k = 1:a
Y(frombus(k),tobus(k)) =
can be understood by taking a deep insight to a power system Y(frombus(k),tobus(k)) - y(k)/tap(k);
and can be seen that the power is injected into a bus from Y(tobus(k),frombus(k)) =
generators and the load are being tapped from it. It is possible Y(frombus(k),tobus(k));
to have buses having either only generator or loads. There are end
also some buses which have both the generators and loads,
% diagonal element
and there also exists some other may which have only static for m = 1:nbus
capacitors for reactive power compensation. for n = 1:a
if frombus(n) == m
II. F ORMULATION Y(m,m) = Y(m,m) + y(n)/(tap(n)ˆ2) + b(n);
To start, the program will open and read an existing *.txt elseif tobus(n) == m
Y(m,m) = Y(m,m) + y(n) + b(n);
containing the line data to be used - the matlab funtion dlmread end
will be used. The program will start reading from the 3rd row end
to the last, bypassing the titles of the table header per column. end
The read data will be saved on the Ldata matrix. The numerical fprintf(fileOutput,’Bus Admittance
values per column will be assigned to a specific data handler Matrix\n\n’);
[d,e] = size(Y);
(frombus, tobus,r,x,...). Using the data, the Ybus matrix will
be solved using the Ybus formulation algorithm.
clear;
The MATLAB program will now save and display the
clc; calculated Ybus matrix on an external *.txt file created by
%diary on the user inline. Instead of using the diary function, the group
%diary output.txt created their own inline program function to display the result.
filename = ’Output.txt’;
fileOutput = fopen(filename,’wt+’); %% OUTPUT TO TXT FILE
for o1 = 1:d
Ldata = dlmread(’table3.txt’,’’,3,0); for o2 = 1:e
[a,c] = size(Ldata); if signIm(Y(o1,o2)) < 0
fprintf(fileOutput,’%10.5f -%9.5fj
’,real(Y(o1,o2)),abs(imag(Y(o1,o2))));
elseif signIm(Y(o1,o2)) == 0
fprintf(fileOutput,’%10.5f ’,real(Y(o1,o2)));
else
fprintf(fileOutput,’%10.5f +%9.5fj
’,real(Y(o1,o2)),abs(imag(Y(o1,o2))));
end
end
fprintf(fileOutput,’\n’);
end

fprintf(fileOutput,’\nby:\n’);
fprintf(fileOutput,’Julius Bacton\n’);
fprintf(fileOutput,’Nathannaell Terio’);
winopen(filename);
%disp(Y);
%diary off

III. R ESULTS AND D ISCUSSIONS


Fig. 3. 16-Bus Input

The very aim of this activity is to write a program using that


can formulate an Admittance Matrix. The input will be taken
from an external txt file and an output will also be saved on a
txt file. The program is verified by solving a sample problem as
shown in Fig.1. Comparing the results shown in Fig.2. which
is the output solved manually and Fig.6. which is solved using
the program - we can see an almost identical results.

Fig. 4. Cropped 16-Bus Output

Fig. 1. Sample Problem

Fig. 2. Theoretical Output Fig. 5. 3-Bus Input


Fig. 6. 3-Bus Output

IV. C ONCLUSION
Thus the admittance matrix for thegiven example were
determined and verified using Matlab. The group were able to
develop a program that solves for the bus admittance matrix
of the problem.
R EFERENCES
[1] https://circuitglobe.com/bus-admittance-matrix.html
[2] https://academia.edu/24224576/Formation-of-Bus-Admittance-Matrix-
using-MATLAB-Software

You might also like