Modelling of Load Flow Analysis in MatlabSimulink Software2
Modelling of Load Flow Analysis in MatlabSimulink Software2
Gauss-Seidel,
Newton Raphson
& Fast
Decoupled Load
Flow Method
Using
By Dr. J. A. Laghari MATLAB
Power System Modelling in MATLAB
Load Flow Analysis Methods
Load flow analysis is very important in the determination
of losses, and voltages at different buses in the power system.
Chapter 6,
Example 6.7.
Example
Figure shows the one line diagram of a simple power three bus power
system with generation at bus 1. The magnitude of voltage at bus 1 is
adjusted to 1.05 p.u. the scheduled loads at buses 2 and 3 are as marked on
the diagram. Line impedances are marked in per unit on a 100 MVA base
and the line charging susceptance are neglected.
R 0.02 p.u.
X L 0.04 p.u.
Actual Impedance is given by
V 2 base
Z base
Sbase
Base MVA is given as 100 MVA and Base kV is not given in the example. However,
to put the values in Matlab/Simulink, we need the base kV. The default base kV in
Matlab/Simulink is 25 kV. Therefore, we assume base kV as 25 kV. However, you may
choose any base kV value, it will have no effect in the Matlab/Simulink Results.
Solution:
Actual Impedance of the Transmission Line 1-2
V 2 base (25 103 ) 2 625 106
Z base 6.25
Sbase 100 10 6
100 10 6
R 0.02 6.25
The Matlab/Simulink ask the value of Inductance in
Henry instead of Inductive Reactance. We know that
X L 2 f L
Solution:
Actual Impedance of the Transmission Line 1-2
XL
L
2 f
Therefore, actual inductance will be
R 0.01 6.25
Therefore, actual inductance will be
R 0.0125 6.25
Therefore, actual inductance will be
1 1 1 0.02 j 0.04
y12 10 j 20
Z12 0.02 j 0.04 0.02 j 0.04 0.02 j 0.04
1 1 1 0.01 j 0.03
y13 10 j 30
Z13 0.01 j 0.03 0.01 j 0.03 0.01 j 0.03
1 1 1 0.0125 j 0.025
y23 16 j 32
Z 23 0.0125 j 0.025 0.0125 j 0.025 0.0125 j 0.025
Solution:
1. Calculation of Phase Voltages:
Y11 y12 y13 10 j 20 10 j 30 20 j 50
Y22 y12 y23 10 j 20 16 j 32 26 j 52
Y33 y13 y23 10 j 30 16 j 32 26 j 62
Y12 Y21 y12 10 j 20
Y13 Y31 y13 10 j 30
Y22 V2
Solution:
1. Calculation of Phase Voltages:
2.566 j1.102
1 ( 10 j 20 ) (1.05 j 0.0)
V2(1) 1.0 j 0.0
26 j52
( 16 j 32 ) (1.0 j 0.0)
Y33 V3
1.386 j 0.452
1 ( 10 j 30 ) (1.05 j 0.0)
V3(1) 1.0 j 0.0
26 j 62
( 16 j 32 ) ( 0.9825 j 0.0310 )
Solution:
1. Calculation of Phase Voltages:
V 3
(1)
1.0011 j 0.0353
For second iteration,
1 P2 jQ2 (1)
V 2
( 2)
(1)*
Y21 V1 Y23 V3
(0)
Y22 V2
2.566 j1.102
1 ( 10 j 20 ) (1.05 j 0.0)
V2( 2) 0.9825 j 0.031
26 j52
( 16 j 32 ) (1.0011 j 0.0353 )
Solution:
1. Calculation of Phase Voltages:
V 3
(1)
1.0011 j 0.0353
For second iteration,
1 P2 jQ2 (1)
V 2
( 2)
(1)*
Y21 V1 Y23 V3
(0)
Y22 V2
2.566 j1.102
1 ( 10 j 20 ) (1.05 j 0.0)
V2( 2) 0.9825 j 0.031
26 j52
( 16 j 32 ) (1.0011 j 0.0353 )
V 2
( 2)
0.9816 j 0.052
Solution:
1. Calculation of Phase Voltages:
1 P3 jQ3 ( 2)
V 3
( 2)
(1)*
Y31 V1 Y32 V2
(0)
Y33 V3
1.386 j 0.452
1 (10 j30) (1.05 j 0.0)
V3( 2) 1.0011 j 0.0353
26 j 62
( 16 j 32 ) ( 0.9816 j 0.052 )
V 2
(7)
0.98 j 0.06 0.98183 3.5035 p.u
V 3
(7)
1.0 j 0.05 1.00125 2.8624 p.u
Solution:
2. Slack Bus Real and Reactive Power:
With the knowledge of all bus voltages, the slack bus power is obtained from:
It converts impedances to admittances and obtains the bus admittance matrix.
Solution:
Gauss Seidel Power Flow Programs:
lfgauss MTLAB Function:
This program obtains the power flow solution by the Gauss-
Seidel method and requires the files named busdata and linedata.
Loads and generation are converted to per unit quantities on the base MVA selected.
Solution:
Gauss Seidel Power Flow Programs:
busout MTLAB Function:
This program produces the bus output result in a tabulated form.
The bus output result includes the voltage magnitude and angle, real and
reactive power of generators and loads and the shunt capacitor/reactor Mvar.
Total generation and total load are also included as outlined in the example case.
Solution:
Gauss Seidel Power Flow Programs:
lineflow MTLAB Function:
This program prepares the line output data.
It is designed to display the active and reactive power flow entering
the line terminals and line losses as well as the net power at each bus.
The total real and reactive losses in the system are also included.
The output of this portion is also shown in the sample case.
Solution:
Gauss Seidel Power Flow Programs:
Data Preparation:
In order to perform a power flow analysis by the Gauss-seidel method
in the MATLAB environment, the following variables must be defined:
Power system base MVA, Power mismatch accuracy,
Acceleration factor, and Maximum number of iterations.
The name (in lowercase letters) reserved for these variables are basemava,
accuracy, accel, and maxiter iteratively. Typical values are as follows:
basemava=100; accuracy = 0.001;
accel = 1.8; maxiter = 80
The initial step in the preparation of input file is the numbering of each bus.
Solution:
Gauss Seidel Power Flow Programs:
In addition, the following data files are required.
Bus Data File – busdata:
The format for the bus is shown to facilitate the required data for each case in a
single row. The information required must be included in a matrix called busdata.
Column 1 is the bus number.
Column 2 contains the bus code.
Column 3 and 4 are voltage magnitude in per unit and phase angles in degrees.
Column 5 and 6 are load MW and Mvar.
Column 7 through 10 are MW and Mvar, minimum Mvar and maximum Mvar of
generation in that order. The last column is the injected Mvar of shunt capactiors.
Solution:
Gauss Seidel Power Flow Programs:
Bus Data File – busdata:
The bus code entered in column 2 is used for identifying load, voltage controlled
and slack buses as outlined below:
1 = This code is used for slack bus. The only necessary information for this bus
is the voltage magnitude and its phase angle.
0 = This code is used for load buses. The loads are entered positive in MW & Mvars.
2 = This code is used for the voltage controlled buses. For this bus, voltage
magnitude, real power generation in megawatts, and the minimum and maximum
limits of the megavar demand must be specified.
Solution:
Gauss Seidel Power Flow Programs:
Line Data File – linedata:
Lines are identified by the node-pair method. The information required must be
incuded in a matrix called linedata.
The last column is for the transformer tap setting for lines, 1 must be entered in
this column.
Method # 2
What Changes are required for Newton Raphson
Load Flow Method
Solution:
Power Flow Programs:
To write the MATLAB Program for Newton Raphson Load
Flow Method, Only change the lfgauss with lfnewton Matlab
function already developed by Hadi Sadat in his book.
Method # 3
What Changes are required for Fast Decoupled
Load Flow Method
Solution:
Power Flow Programs:
To write the MATLAB Program for Fast Decoupled Load
Flow Method, Only change the lfnewton with decouple Matlab
function already developed by Hadi Sadat in his book.
Solution:
IEEE 9 Bus Test System Data:
Bus2 T2 Bus7 Bus8 Bus9 T3 Bus3
250 MVA 150 MVA
Gen 2 19 kV / 230 kV 13.8 kV / 230 kV
Gen 3
Bus5 Bus6
Bus4
T1
100 MVA
16.5 kV / 230 kV
Bus1
Gen1
Solution:
IEEE 14 Bus Test System Data:
Solution:
IEEE 30 Bus Test System Data:
Solution:
IEEE 69 Bus Test System Data:
THANK YOU VERY MUCH