Power Systems and Simulation Laboratory: Department of Electrical and Electronics Engineering

POWER SYSTEMS AND SIMULATION
LABORATORY
YEAR: III/IV SEM: II
STUDENT MANUAL
Department of
ELECTRICAL AND ELECTRONICS ENGINEERING
Prepared by
Mr. Anil Kumar Chinta
Assistant Professor
Department of Electrical and Electronics Engineering
DHANEKULA INSTITUE OF ENGINEERING AND

TECHNOLOGY
(Approved by AICTE, New Delhi, Affliated to JNTUK, KAKINADA)
GANGURU, VIJAYAWADA
INDEX
POWER SYSTEMS AND SIMULATION LABORATORY
Sno Description
1 JNTU KAKINADA SYLLABUS
2 LIST OF EXPERIMENTS
3 INSTITUTE VISION, MISSION, DEPARTMENT VISION, MISSION, PEO’S
4 PO’S, PSO’S
5 COURSE OBJECTIVES, COURSE OUTCOMES
6 MAPPINGS – CO VS PO, CO VS PSO
7 DO’S AND DON’TS
8 SEQUENCE IMPEDANCES OF 3 PHASE TRANSFORMER
9 SEQUENCE IMPEDANCES OF 3 PHASE ALTERNATOR BY FAULT ANALYSIS
10 SEQUENCE IMPEDANCES OF 3 PHASE ALTERNATOR BY DIRECT METHOD
11 ABCD PARAMETERS OF TRANSMISSION NETWORK
12 PERFORMANCE OF LONG TRANSMISSION LINE WITHOUT COMPENSATION
13 PERFORMANCE OF LONG TRANSMISSION LINE WITH SHUNT
COMPENSATION
14 DETERMINATION OF YBUS USING DIRECT INSPECTION METHOD
15 LOAD FLOW STUDIES BY USING GAUSS SEIDEL METHOD
16 LOAD FLOW STUDIES BY USING NEWTON RAPSHON METHOD
17 ECONOMIC LOAD DISPATCH WITH & WITHOUT LOSSES
18 LOAD FREQUENCY CONTROL OF A TWO AREA POWER SYSTEM WITHOUT &
WITH PI CONTROLLER
19 TRANSIENT STABILITY ANALYSIS
R-20 Syllabus for EEE-JNTUK w.e.f. 2020 –21
JAWAHARLAL NEHRU TECHNOLOGICAL UNIVERSITY KAKINADA

KAKINADA–533003, Andhra Pradesh, India
DEPARTMENT OF ELECTRICAL AND ELECTRONICS ENGINEERING
L T P C
III Year – II SEMESTER 0 0 3 1.5
POWER SYSTEMS AND SIMULATION LAB
Course Objectives:
To impart the practical knowledge of functioning of various power system components and
determination of various parameters and simulation of load flows, transient stability, LFC and Economic
dispatch.
Any of 5 experiments are to be conducted from each section:

Section I: Power Systems Lab:
1. Estimation of sequence impedances of 3-phase Transformer

2. Estimation of sequence impedances of 3-phase Alternator by Fault Analysis
3. Estimation of sequence impedances of 3-phase Alternator by Direct method
4. Estimation of ABCD parameters on transmission line model
5. Performance of long transmission line without compensation
6. Performance of long transmission line with shunt compensation
7. Analyze the Ferranti effect on long transmission line
Section II: Simulation Lab

8.Determination of Ybus using direct inspection method
9.Load flow solution of a power system network using Gauss-Seidel method
10. Load flow solution of a power system network using Newton Raphson method.
11. Formation of Zbus by building algorithm.
12. Economic load dispatch with & without losses
13. Load frequency control of a two area Power System without & with PI controller
14. Transient Stability analysis of single machine connected to an infinite bus (SMIB) using equal area
criterion.
Course Outcomes:
After the completion of the course the student should be able to:
 Estimate the sequence impedances of 3-phase Transformer and Alternators
 Evaluate the performance of transmission lines
 Analyse and simulate power flow methods in power systems
 Analyse and simulate the performance of PI controller for load frequency control.
 Analyse and simulate stability studies of power systems
LIST OF EXPERIMENTS DIET/7.5.1/FT L05
DHANEKULA INSTITUTE OF ENGINEERING & TECHNOLOGY

POWER SYSTEM AND SIMULATION LABORATORY
LIST OF EXPERIMENTS
Cycle
S No Name of the Experiment Connected CO
No
Experiments as per curriculum

1 Estimation of sequence impedances of 3-phase Transformer I R20C318.1
Estimation of sequence impedances of 3-phase Alternator by
2 I R20C318.1
Fault Analysis
3 Estimation of ABCD parameters on transmission line model I R20C318.1

4 Performance of long transmission line with shunt compensation I R20C318.1
5 Analyze the Ferranti effect on long transmission line I R20C318.1
6 Determination of Ybus using direct inspection method I R20C318.1
Load flow solution of a power system network using Gauss-
7 I R20C318.1
Seidel method
Load flow solution of a power system network using Newton
8 I R20C318.1
Raphson method
Load frequency control of a two area Power System without &
9 I R20C318.2
with PI controller
Transient Stability analysis of single machine connected to an
10 I R20C318.2
infinite bus (SMIB) using equal area criterion.
Power Systems Lab Student Manual
------------------------------------------------------------------------------------------------------------------------
INSTITUTION VISION
Pioneering Professional Education through Quality
INSTITUTION MISSION
 Providing Quality Education through state-of-art infrastructure, laboratories and committed
staff.
 Moulding Students as proficient, competent, and socially responsible engineering personnel

with ingenious intellect.
 Involving faculty members and students in research and development works for betterment of
society.
DEPARTMENT VISION
Emerge as Quality Human Resource Provider for Industry and Society in the field of
Electrical & Electronics Engineering
DEPARTMENT MISSION
 Providing Quality Education through State-of-art resources.
 To develop innovative, proficient Electrical engineers.
 Promoting Ethical and moral values among the students so as to make them responsible
professionals for the society.
PROGRAM EDUCATIONAL OBJECTIVES(PEO’s)
PEO1: Have strong foundation in Electrical Engineering along with Mathematics, Sciences and
allied Engineering subjects.
PEO2: Possess good problem solving, design skills, capability to use modern engineering
tools, ability to pursue higher education and research.
PEO3: Seek employment in various engineering or technological positions of their interest and
continue to achieve their aspirations through lifelong learning.
PEO4: Exhibit professional and ethical attitude, effective communication skills, Teamwork
and multidisciplinary approach.
---------------------------------------------------------------------------------------------------------------------------------------------
DHANEKULA INSTITUE OF ENGINEERING AND TECHNOLOGY :: DEPARTMENT OF ELECTRICAL
AND ELECTRONICS ENGINEERING
Page | I
------------------------------------------------------------------------------------------------------------------------
PROGRAM OUTCOMES
Engineering Graduates will be able to:
1. Engineering knowledge: Apply the knowledge of mathematics, science, engineering fundamentals, and an
engineering specialization to the solution of complex engineering problems.
2. Problem analysis: Identify, formulate, review research literature, and analyze complex engineering problems
reaching substantiated conclusions using first principles of mathematics, natural sciences, and engineering
sciences.
3. Design/development of solutions: Design solutions for complex engineering problems and design system
components or processes that meet the specified needs with appropriate consideration for the public health and
safety, and the cultural, societal, and environmental considerations.
4. Conduct investigations of complex problems: Use research-based knowledge and research methods including
design of experiments, analysis and interpretation of data, and synthesis of the information to provide valid
conclusions.
5. Modern tool usage: Create, select, and apply appropriate techniques, resources, and modern engineering and IT
tools including prediction and modeling to complex engineering activities with an understanding of the
limitations.
6. The engineer and society: Apply reasoning informed by the contextual knowledge to assess societal, health,
safety, legal and cultural issues and the consequent responsibilities relevant to the professional engineering
practice.
7. Environment and sustainability: Understand the impact of the professional engineering solutions in societal and
environmental contexts, and demonstrate the knowledge of, and need for sustainable development.
8. Ethics: Apply ethical principles and commit to professional ethics and responsibilities and norms of the
engineering practice.
9. Individual and team work: Function effectively as an individual, and as a member or leader in diverse teams,
and in multidisciplinary settings.
10. Communication: Communicate effectively on complex engineering activities with the engineering community
and with society at large, such as, being able to comprehend and write effective reports and design
documentation, make effective presentations, and give and receive clear instructions.
11. Project management and finance: Demonstrate knowledge and understanding of the engineering and
management principles and apply these to one’s own work, as a member and leader in a team, to manage
projects and in multidisciplinary environments.
12. Life-long learning: Recognize the need for, and have the preparation and ability to engage in independent and
life-long learning in the broadest context of technological change.
---------------------------------------------------------------------------------------------------------------------------------------------
Page | II
------------------------------------------------------------------------------------------------------------------------
PROGRAM SPECIFIC OUTCOMES (PSOs)
PSO 1: Ability to design solutions for identified problems by using latest engineering tools like
MATLAB, Simulink, PSPICE, plc etc.
PSO 2: Able to design and develop the Green Electrical systems.
---------------------------------------------------------------------------------------------------------------------------------------------
Page | III
------------------------------------------------------------------------------------------------------------------------
COURSE OBJECTIVES:
At the end of the Course/Subject, the students will:
S.No Course Objectives (Cobs)
Impart the practical knowledge of determination of various parameters of power
R20Cobj318.1
system components
Impart the practical knowledge of Load Frequency Control and Transient
R20Cobj318.2
Stability
COURSE OUTCOMES:
At the end of the Course/Subject, the students will be able to:
S. No Course Outcomes(COs)
Determine the parameters of various power system components occur in power
R20C318.1
system studies
Design the frequency control of a system With and Without Controllers and Transient
R20C318.2
stability
---------------------------------------------------------------------------------------------------------------------------------------------
Page | IV
------------------------------------------------------------------------------------------------------------------------
Course Outcomes vs POS Mapping:

Courses Out
P01 P02 P03 P04 P05 P06 P07 P08 P09 P10 P11 P12
Comes
R20C318.1 - - 2 - 3 - - 2 3 - - -
R20C318.2 - 2 3 - 3 - - 2 3 - - -
Total 0 2 5 0 6 0 0 4 6 0 0 0
Level of
- 2 3 - 3 - - 2 3 - - -
Mapping
Course Outcomes vs PSOs Mapping:
Course Out Comes PSO1 PSO2

R20C318.1 3 3
R20C318.2 3 3
Total 6 6
Level of Mapping 3 3
Justification of Mapping of Course Outcomes with Program Outcomes and Program Specific
Outcomes:
1. R20C318.1 maps highly with PO5, PO9, PSO1, PSO2 as it needs modern tools to determine the
parameters of various power system components occur in power system studies and improves
Individual and teamwork skills. R20C318.1 maps moderately with PO3, PO8 as it requires design of
solutions for determining the parameters of various power systems by applying ethical principles.
2.R20C318.2 maps highly with PO3, PO5, PO9, PSO1 and PSO2 as it requires design by using
modern tools and improves Individual and teamwork skills. R20C318.1 maps moderately with PO2,
PO8 as it requires problem analysis, ethical principles.
---------------------------------------------------------------------------------------------------------------------------------------------
Page | V
------------------------------------------------------------------------------------------------------------------------
DO’S AND DON’TS IN THE LAB

DO’S:-
1. Proper dress code has to be maintained while entering in to the Lab.
2. Students should carry observation notes and record completed in all aspects.
3. Correct specifications of the equipment have to be mentioned in the circuit diagram.
4. Student should be aware of operating the equipment.
5. Students should be at their concerned experiment table, unnecessary moment is restricted.
6. The readings must be shown to the Lecturer In-Charge for verification.
7. Students must ensure that all switches are in the OFF position, all the connections are removed.
8. All patch cords and stools should be placed at their original positions.
DON’Ts:-
1. Don’t come late to the Lab.
2. Don’t enter into the Lab with valuable ornaments like Golden rings, bracelets and bangles.
3. Don’t make or remove the connections with power ON.
4. Don’t switch ON the supply without verifying by the Staff Member.
5. Don’t switch OFF the machine with load.
6. Don’t leave the lab without the permission of the Lecturer In-Charge.
---------------------------------------------------------------------------------------------------------------------------------------------
Page | VI
------------------------------------------------------------------------------------------------------------------------
CIRCUIT DIAGRAM:
Positive or Negative Sequence Impedance:
Zero Sequence Impedance:
----------------------------------------------------------------------------------------------------------------------------- ----------------
Page No
------------------------------------------------------------------------------------------------------------------------
Experiment No: Date:

SEQUENCE IMPEDANCE OF 3 PHASE TRANSFORMER
AIM:
To determine the Positive, Negative and Zero sequence (sequence impedance) of given
three phase transformer
NAME PLATE DETAILS:

3 Phase Transformer
S.No
Specifications Ratings
APPARATUS:
Sl. No Name of Apparatus Range Type Quantity
THEORY:
Need to be written by Students by referring Text Books
PROCEDURE:
1. Connect the given three phase transformer as shown in the Circuit Diagram Fig1
2. After short-circuiting the low voltage side adjusts the voltage on high voltage side with the
help of the autotransformer such that the rated current flows in the windings.
3. Note down the voltage and current.
4. From these readings determine the transformer positive sequence Impedance which is also
equal to negative sequence impedance.
5. Determine the zero sequence impedance of the transformer by making connections as
shown in Circuit DiagramFig2.
6. Note down the voltage and current.
THEORITICAL CALCULATIONS:
Positive and Negative Sequence Impedance, Zero Sequence Impedance,
Z1 = Z2 = V/(3*I) Z0 = V/(3*I)
----------------------------------------------------------------------------------------------------------------------------- ----------------
Page No
------------------------------------------------------------------------------------------------------------------------
R1 = W / I2 R0 = W / (3I2)
X1 = X2 = (Z21 – R21) X0 = (Z20 – R20)
TABULAR FORM:
Positive Sequence Impedance Zero Sequence Impedance
S. No V(volts) I(Amps) Z1 S. V(volts) I(Amps) Z0
No
PRECAUTIONS:
1. Connections are to be made tightly and according to the circuit diagram.
2. There should be no parallax error while taking readings.
3. While making or breaking connections supply has to switched OFF.
RESULT:
----------------------------------------------------------------------------------------------------------------------------- ----------------
Page No
------------------------------------------------------------------------------------------------------------------------
CIRCUIT DIAGRAM:
L-G FAULT:
L-L FAULT:
LL-G FAULT
----------------------------------------------------------------------------------------------------------------------------- ----------------
Page No
------------------------------------------------------------------------------------------------------------------------

SEQUENCE IMPEDANCE OF 3 PHASE ALTERNATOR BY FAULT ANALYSIS
AIM:
To determine the fault currents of an unloaded synchronous generator for
(a) Line to ground fault (L-G Fault)
(b) Line to Line fault (L-L fault)
(c) Double Line to ground fault (LL-G Fault)
NAME PLATE DETAILS:

D.C. Shunt Motor Salient Pole Synchronous Machine
H.P rating KVA Rating
Voltage Voltage
Armature Armature
Current Current
Field Current Field Current
Speed Synchronous Speed
APPARATUS:
Theory:
PROCEDURE:
(a) L-G Fault:
1. Connect the circuit as per the circuit diagram for a line to ground fault on phase R.
2. Calculate the determinate value of the fault current from impedances (Positive, Negative
and Zero sequences).
3. Run the generator rated speed.
4. Increase the field current of excitation so that terminal voltage is constant value.
5. Close the switch to create the L-G fault on Phase R.
6. Note the current and voltage in the ammeter and voltmeter.
7. Open the switch and remove the L-G fault on phase R.
8. Reduce the excitation and open the field circuit switch and switch off the prime mover.
----------------------------------------------------------------------------------------------------------------------------- ----------------
Page No
------------------------------------------------------------------------------------------------------------------------
Note: This voltage must be such that it does not cause the rated current of the machine
to be exceeded.
(b) L-L Fault:

1. Connect the circuit as per the circuit diagram for a line to line fault on any two phases.
5. Close the switch to create the L-L fault on any two phases.
7. Open the switch and remove the L-L fault.
Note: This voltage must be such that it does not cause the rated current of the machine to be
exceeded.
(c) LL-G Fault:

1. Connect the circuit as per the circuit diagram for a line to line to ground fault.
5. Close the switch to create the LL-G fault.
7. Open the switch and remove the LL-G fault on phase R.
Note: This voltage must be such that it does not cause the rated current of the machine to be
exceeded.
----------------------------------------------------------------------------------------------------------------------------- ----------------
Page No
------------------------------------------------------------------------------------------------------------------------
Tabulation:
EB
EY
After Fault
ER
LG- Fault
If
EB
Before Fault
EY
ER
If
VExc
itation
S.No
----------------------------------------------------------------------------------------------------------------------------- ----------------
Page No
------------------------------------------------------------------------------------------------------------------------
EB
EY
After Fault
ER
LL- Fault
If
EB
Before Fault
EY
ER
If
VExc
itation
S.No
----------------------------------------------------------------------------------------------------------------------------- ----------------
Page No
------------------------------------------------------------------------------------------------------------------------
EB
EY
After Fault
ER
LLG- Fault
If
EB
Before Fault
EY
ER
If
VExc
itation
S.No
PRECAUTIONS:
RESULT:
----------------------------------------------------------------------------------------------------------------------------- ----------------
Page No
------------------------------------------------------------------------------------------------------------------------
Circuit Diagram:
Figure1: Circuit for open circuit test for calculating A and C values
Figure2: Circuit for short circuit test for calculating B and D values
----------------------------------------------------------------------------------------------------------------------------- ----------------
Page No
------------------------------------------------------------------------------------------------------------------------

ABCD PARAMETERS OF TRANSMISSION NETWORK
AIM:
To determine ABCD constants of 3-phase transmission line with Distributed Connection
APPARATUS:
THEORY:
PROCEDURE:
1. Make connections shown in figure-2 below. Now the receiving end terminals are kept open.
Turn the phase shift knob in Source-1 to zero.
2. Switch ON power to the module. Slowly turn the variac in Source-1 and apply a voltage of
about 100kV. This is indicated by the voltmeter at the sending end. Now observe the voltmeter
reading at the receiving end. It reads a value much higher than the input voltage. This effect is
known as Ferranti effect. This effect occurs because of high shunt capacitance of the
transmission line. This effect is predominant only when the line is not loaded.
3. Now record VS and VR. Magnitude of these two are read by the respective voltmeters and
their phase angles are read by respective phase angle meters. Since VSis taken as reference
phasor, its phase angle will be zero. If VR is lagging, as indicated by the lag LED at the
receiving end, then, its phase angle should be taken as negative. If VR is leading as indicated by
the lead LED at the receiving end, then, its phase angle should be taken as positive.
4. Also record IS magnitude as read by the ammeter and its phase angle as read by the phase
angle meter. Now calculate A and C. Turn variac to zero and switch OFF power supply.
5. Now make connections shown in figure-3. In this case, the receiving end terminals are
shorted and we need to measure the receiving end and sending end currents.
6. Switch ON the power supply. Slowly apply about 100kV at the sending by operating the
variac in Source-1. Now measure both IS and IR. Take care to measure phase angle properly. The
phase angle of VSis zero since it is taken as reference phasor. Now Calculate B and D.
----------------------------------------------------------------------------------------------------------------------------- ----------------
Page No
------------------------------------------------------------------------------------------------------------------------
For A and C: (figure-2)
For B and D: (figure-3)
Formulae: A = [VS/VR]IR=0; C = [IS/VR]IR=0;

B = [VS/IR]VR=0; D = [IS/IR]VR=0
Tabular form:
For A and C: For B and D:
S.No VS IS VR IR S.No VS IS VR IR
PRECAUTIONS:
RESULT:
----------------------------------------------------------------------------------------------------------------------------- ----------------
Page No
------------------------------------------------------------------------------------------------------------------------
Circuit Diagram:
Figure1: Power transfer capability of Transmission line without compensation
----------------------------------------------------------------------------------------------------------------------------- ----------------
Page No
------------------------------------------------------------------------------------------------------------------------

PERFORMANCE OF LONG TRANSMISSION LINE WITHOUT COMPENSATION
AIM:
To observe the Performance of long transmission line without compensation
APPARATUS:
THEORY:
PROCEDURE:
1. Make connections as shown.
2. Turn the phase shift control knobs to zero in both Source-1 and Source-2.
3. Switch ON Source-1 and slowly increase its output voltage by operating the variac until the
sending end voltage is 100 KV as indicated by the voltmeter on the sending end side.
4. Switch ON Source-2 and slowly increase its output voltage to 100 KV by operating its variac.
This is indicated by the voltmeter on the receiving end side. The phase angle meter on the
receiving side displays the load angle δ. This can be changed by operating the phase shift
control knob in Source-2.
5. Turn the phase shift control knob of Source-2 to 1 position and record δ shown by the phase
angle meter on the receiving side in Table-1. If the magnitude of receiving end voltage has
changed, set it 200 KV by operating the variac in Source-2. Record sending and receiving end
active powers in Table-1.
6. Repeat the previous step for each position of phase shift control knob of Source-2 up to 400.
7. The maximum power that can be transferred without compensation corresponds to δ = 300
which is the practical limit.
----------------------------------------------------------------------------------------------------------------------------- ----------------
Page No
------------------------------------------------------------------------------------------------------------------------
Hardware Connection diagram:
Tabular form:
Vs, kV VR, kV δ0 PS, MW PR, MW
PRECAUTIONS:
RESULT:
----------------------------------------------------------------------------------------------------------------------------- ----------------
Page No
------------------------------------------------------------------------------------------------------------------------
Circuit Diagram:
Figure1: Inductive Load connected at the receiving end
Figure2: Capacitive Load connected at the receiving end
----------------------------------------------------------------------------------------------------------------------------- ----------------
Page No
------------------------------------------------------------------------------------------------------------------------

PERFORMANCE OF LONG TRANSMISSION LINE WITH SHUNT COMPENSATION
AIM:
To observe the Performance of long transmission line with shunt compensation
APPARATUS:
THEORY:
PROCEDURE:
1. Make connections as shown in above figure but without connecting Source-2 connect
inductive load at the receiving end.
2. Switch ON Source-1 and slowly increase its output voltage by operating the variac until the
sending end voltage is 100 KV as indicated by the voltmeter on the sending end side. Switch on
the inductive load.
3. When the inductive load is switched on, the receiving end voltage decreases. Note down this
voltage.
4. Disconnect the load. connect a capacitor in parallel with the load and then switch on power to
the load keeping the sending end voltage same as it was earlier. Now, we can see that the
receiving end voltage improves.
----------------------------------------------------------------------------------------------------------------------------- ----------------
Page No
------------------------------------------------------------------------------------------------------------------------
Hardware Connection diagram:
Tabular form:
With Inductive load:
S.No L VS IS VR IR
With Capacitive load:

S.No C VS IS VR IR
PRECAUTIONS:
RESULT:
----------------------------------------------------------------------------------------------------------------------------- ----------------
Page No
------------------------------------------------------------------------------------------------------------------------
Circuit Diagram:
Figure1: Circuit for observing Ferranti Effect
----------------------------------------------------------------------------------------------------------------------------- ----------------
Page No
------------------------------------------------------------------------------------------------------------------------

ANALYZE THE FERRANTI EFFECT ON LONG TRANSMISSION LINE
AIM:
To analyze the Ferranti effect on long transmission line
APPARATUS:
THEORY:
PROCEDURE:
1. Make connections shown in figure-1 above. Now the receiving end terminals are kept open.
Turn the phase shift knob in Source-1 to zero.
2. Switch ON power to the module. Slowly turn the variac in Source-1 and apply a voltage of
about 100kV. This is indicated by the voltmeter at the sending end. Now observe the voltmeter
reading at the receiving end. It reads a value much higher than the input voltage. This effect is
known as Ferranti effect. This effect occurs because of high shunt capacitance of the
transmission line. This effect is predominant only when the line is not loaded.
Hardware Connections:
----------------------------------------------------------------------------------------------------------------------------- ----------------
Page No
------------------------------------------------------------------------------------------------------------------------
Tabular form:
S.No VS IS VR IR
PRECAUTIONS:
RESULT:
----------------------------------------------------------------------------------------------------------------------------- ----------------
Page No
------------------------------------------------------------------------------------------------------------------------

YBUS USING DIRECT INSPECTION METHOD
AIM:
To determine the Ybus by using Direct Inspection method.
APPARATUS REQUIRED: Personal Computer with MATLAB software
THEORY:
Problem:
Form the Ybus by using Direct Inspection method for the data shown below.
Element Bus code Self impedance in p.u
1 0–1 0.15
2 0–2 0.1
3 0–3 0.3
4 0–4 0.2
5 1–2 0.4
6 2–3 0.15
7 3–4 0.2
8 1–4 0.1
PROCEDURE:
1. Switch ON the Personal Computer and open the Start Menu and Open “Matlab Icon”.
2. Click on the New Script on the Main Window.
3. A new window opens, Type the program in the window.
4. After completion, save the file by clicking on Save.
5. Then click on the Run button, which executes the program.
6. Then note down the output displayed on the command window.
7. Then close the window.
PROGRAM:
% Ybus by Direct Inspection Method
clc
z01 = 0.15
z02 = 0.1
z03 = 0.3
----------------------------------------------------------------------------------------------------------------------------- ----------------
Page No
------------------------------------------------------------------------------------------------------------------------
z04 = 0.2
z12 = 0.4
z23 = 0.15
z34 = 0.2
z14 = 0.1
disp(‘Obtain the values of self admittances from self impedances’)

y01 = 1/z01
y02 = 1/z02
y03 = 1/z03
y04 = 1/z04
y12 = 1/z12
y23 = 1/z23
y34 = 1/z34
y14 = 1/z14
disp('Formulate Y-Bus Matrix')

Y11 = y01 + y12 + y14;
Y12 = -y12;
Y13 = 0;
Y14 = 0;
Y21 = Y12;
Y22 = y02+y12+y23;
Y23 = -y23;
Y24 = 0;
Y31 = Y13;
Y32 = Y23;
Y33 = y03+y23+y34;
Y34 = -y34;
Y41 = Y14;
Y42 = Y24;
Y43 = Y34;
Y44 = y04+y34+y14;
disp(‘Ybus matrix is’)

YBus = [Y11 Y12 Y13 Y14; Y21 Y22 Y23 Y24; Y31 Y32 Y33 Y34; Y41 Y42 Y43 Y44]
RESULT:
----------------------------------------------------------------------------------------------------------------------------- ----------------
Page No
------------------------------------------------------------------------------------------------------------------------
Experiment No: 9 Date:

LOAD FLOW STUDIES BY USING GAUSS SEIDEL METHOD
AIM:
To determine the power flow solution using Gauss Seidel method
THEORY:
Problem:
In the system shown in figure below, bus 1 is slack bus and remaining busses are load busses.
By taking a flat voltage profile, determine the bus voltage at the end of iterations.
Bus data:
Bus code P Q V Remarks
1 - - 1.04 Slack bus
2 -0.5 - 1.0 ∟00 PV bus
3 1.0 -0.55 - PQ bus
4 -0.3 0.1 - PQ bus
PROCEDURE:
----------------------------------------------------------------------------------------------------------------------------- ----------------
Page No
------------------------------------------------------------------------------------------------------------------------

PROGRAM:
%Gauss Seidal Method clc
clc
y12 = 0.873-6.55j;
y13 = 1-3j;
y23 = 0.666-2j;
y24 = 1-3j;
y34 = 2-6j;
p2 = 0.5;
p3 = -1.0;
p4 = 0.3;
q2 = -0.2
q3 = 0.55;
q4 = -0.1;
disp('STEP1: Formulate Y-Bus Matrix')

Y11 = y12+y13;
Y12 = -y12;
Y13 = -y13;
Y14 = 0;
Y21 = -y12;
Y22 = y12+y23+y24;
Y23 = -y23;
Y24 = -y24;
Y31 = -y13;
Y32 = -y23;
Y33 = y13+y23+y34;
Y34 = -y34;
Y41 = 0;
Y42 = -y24;
Y43 = -y34;
Y44 = y24+y34;
YBus = [Y11 Y12 Y13 Y14; Y21 Y22 Y23 Y24; Y31 Y32 Y33 Y34; Y41 Y42 Y43 Y44]
disp('STEP2: Initialize bus voltages')

v1 = 1.04
v2 = 1.04
v3 = 1.0
v4 = 1.0
disp('STEP3: Calculate Q2 value for PV Bus')
Q2 = v2*[(Y21*v1)+(Y22*v2)+(Y23*v3)+(Y24*v4)];
Q2cal = -imag(Q2)
----------------------------------------------------------------------------------------------------------------------------- ----------------
Page No
------------------------------------------------------------------------------------------------------------------------
disp('STEP4: Calculate V2 value')

V2 = (1/Y22)*[((p2-(Q2cal*j))/v2)-(Y21*v1)-(Y23*v3)-(Y24*v4)]
delta2 = angle(V2)
V2new = v2*(cos(delta2)+j*sin(delta2))
V3 = (1/Y33)*[((p3-(q3*j))/v3)-(Y31*v1)-(Y32*V2new)-(Y34*v4)]
V4 = (1/Y44)*[((p4-(q4*j))/v4)-(Y41*v1)-(Y42*V2new)-(Y43*V3)]
RESULT:
----------------------------------------------------------------------------------------------------------------------------- ----------------
Page No
------------------------------------------------------------------------------------------------------------------------

LOAD FLOW STUDIES BY USING NEWTON- RAPHSON METHOD
AIM:
To determine the power flow solution using Newton- Raphson method
THEORY:
Problem:
Consider the 2 bus system having series impedance of 0.0839 + j0.5183p.u and a total shunt
admittance of 0.0636jpu. The specified quantities at the buses are given below:
Bus PD QD PG QG Voltage Specified
1 2 1 - - V1=1.05
2 0 0 0.3 0.1 Unspecified
PROCEDURE:
PROGRAM:
%Newton Raphson Method
clc
z12 = 0.0839+0.5183j;
hlc = 0.0636j;
----------------------------------------------------------------------------------------------------------------------------- ----------------
Page No
------------------------------------------------------------------------------------------------------------------------
disp('STEP1: Formulate Y-Bus matrix')

Y11 = 1/z12 + hlc;
Y12 = -1/z12;
Y21 = -1/z12;
Y22 = 1/z12 + hlc;
yBus = [Y11 Y12;

Y21 Y22]
Y11m = abs(Y11);
Y12m = abs(Y12);
Y21m = abs(Y21);
Y22m = abs(Y22);
Y11a = angle(Y11);
Y12a = angle(Y12);
Y21a = angle(Y21);
Y22a = angle(Y22);
disp('STEP2: Initialize bus voltages')

v1 = 1.05
v2 = 1.0
d1 = 0
d2 = 0
X0 = [d2;
v2]
disp('STEP3: Calculate P2cal, Q2cal, dP2, dQ2')

P2cal = v2*((v1*Y21m*cos(Y12a+d2-d1))+(v2*Y22m*cos(Y22a+d2-d1)))
P2spec = -0.3
dP2 = P2spec-P2cal
Q2cal = -v2*((v1*Y21m*sin(Y12a+d1-d2))+(v2*Y22m*sin(Y22a+d2-d1)))
disp('Check for q limit')
disp('Q2cal < Q2min')
Q2spec = -0.1
dQ2 = Q2spec-Q2cal
disp('STEP4: Form Jacobian matrix')

J11 = (v2*v1*Y21m*sin(Y12a+d1-d2))+(v2*v2*Y22m*0);
J12 = (v1*Y21m*cos(Y12a+d1-d2))+(2*v2*Y22m*cos(Y22a));
J21 = (v2*v1*Y21m*cos(Y12a+d1-d2))-(v2*v2*Y22m*0);
J22 = (-v1*Y21m*sin(Y12a+d1-d2))-(2*v2*Y22m*sin(Y22a));
Jmatrix = [J11 J12;
J21 J22]
Jinv = inv(Jmatrix)
disp('STEP5: Compute dX')

pq = [dP2;
----------------------------------------------------------------------------------------------------------------------------- ----------------
Page No
------------------------------------------------------------------------------------------------------------------------
dQ2];
dX = Jinv*pq
delta2 = dX(1,1);
dV2 = dX(2,1);
Xnew = X0+dX
v2new = v2+dV2
Delta2new = d2+delta2
RESULT:
----------------------------------------------------------------------------------------------------------------------------- ----------------
Page No
------------------------------------------------------------------------------------------------------------------------

(a) LOAD FREQUENCY CONTROL WITHOUT CONTROL
AIM:
To become familiar with modelling and analysis of the frequency and tie-line flow
dynamics of a two area power system without control for Load Frequency Control (LFC).
THEORY:
PROBLEM:
A two area system connected by a tie line has the following parameters on a 1000MVA
common base. The unit are operated in parallel at the nominal frequency of 50Hz. The
synchronizing power coefficient is computed from the initial operating condition and is given to
be Ps = 0.2PU. A load change of 187.5MW occurs in Area1. Construct the Simulink block
diagram and obtain the frequency deviation response.
AREA 1 2
Speed regulation R1 = 0.05 R2 = 0.0625
Frequency sensitive load coefficient D1 = 0.6 D2 = 0.9
Inertia constant H1 = 5 H2 = 4
Base Power 1000MVA 1000MVA
Governor Time Constant Tg1 = 0.2sec Tg2 = 0.3sec
Turbine Time Constant Tt1 = 0.5sec Tt2 = 0.6sec
PROCEDURE:
1. Enter the command window of the MATLAB.
2. Create a new Model by selecting File - New – Model
3. Pick up the blocks from the simulink library browser and form a block diagram.
4. After forming the block diagram, save the block diagram.
5. Double click the scope and view the result.
----------------------------------------------------------------------------------------------------------------------------- ----------------
Page No
------------------------------------------------------------------------------------------------------------------------
Simulink Model for Single Area with LFC Without Control:
MODEL GRAPH:
Frequency plot of Area-1:
----------------------------------------------------------------------------------------------------------------------------- ----------------
Page No
------------------------------------------------------------------------------------------------------------------------
Tie Line Power:
Mechanical Input Power of Area-1:
RESULT:
----------------------------------------------------------------------------------------------------------------------------- ----------------
Page No
------------------------------------------------------------------------------------------------------------------------

(b) LOAD FREQUENCY CONTROL WITH CONTROL
AIM:
To become familiar with modelling and analysis of the frequency and tie-line flow
dynamics of a Two Area Power System with control for Load Frequency Control (LFC).
APPARATUS: Personal Computer with MATLAB software
THEORY :
PROBLEM:
A two area system connected by a tie line has the following parameters on a 1000MVA
common base. The unit are operated in parallel at the nominal frequency of 50Hz. The
synchronizing power coefficient is computed from the initial operating condition and is given to
be Ps = 0.2PU. A load change of 187.5MW occurs in Area1. Construct the Simulink block
diagram and obtain the frequency deviation response.
AREA 1 2
Speed regulation R1 = 0.05 R2 = 0.0625
Frequency sensitive load coefficient D1 = 0.6 D2 = 0.9
Inertia constant H1 = 5 H2 = 4
Base Power 1000MVA 1000MVA
Governor Time Constant Tg1 = 0.2sec Tg2 = 0.3sec
Turbine Time Constant Tt1 = 0.5sec Tt2 = 0.6sec
PROCEDURE:
4. After forming the block diagram, save the block diagram and click run.
----------------------------------------------------------------------------------------------------------------------------- ----------------
Page No
------------------------------------------------------------------------------------------------------------------------
Simulink Model for Two Area LFC With Control:
MODEL GRAPH:
----------------------------------------------------------------------------------------------------------------------------- ----------------
Page No
------------------------------------------------------------------------------------------------------------------------
Tie Line Power:
RESULT:
----------------------------------------------------------------------------------------------------------------------------- ----------------
Page No
------------------------------------------------------------------------------------------------------------------------

TRANSIENT STABILITY ANALYSIS
AIM:
To determine the transient stability of a single machine connected to infinite bus.
APPARATUS: Personal Computer with MATLAB Software

THEORY:
PROCEDURE:
4. After forming the block diagram, save the block diagram and click run.
PROBLEM:
A 20 MVA, 50 Hz generator delivers 18 MW over a double circuit line to an infinite
bus. The generator has KE of 2.52 MJ/MVA at rated speed. The generator transient reactance is
X’d=0.35 p.u. Each transmission circuit has R=0 and a reactance of 0.2 pu on a 20 MVA
base.|E’|=1.1 p u and infinite bus voltage V=1.0/_00.A three-phase short circuit occurs at the
mid point of one of the transmission lines. Plot swing curves with fault cleared by simultaneous
opening of breakers at both ends of the line at 2.5 cycles and 6.25 cycles after the occurrence of
fault. Also plot the swing curve over the period of 0.5s if the fault is sustained.
Note: Before running simulation, integrator 1 has to be initialized to pre fault value of δ, i.e., δ0.
This can be done by double-clicking on integrator 1 block and changing the initial value from 0
to δ0 (in radians). Also double click the switch block and change the threshold value from 0 to
the fault clearing time (in sec).
----------------------------------------------------------------------------------------------------------------------------- ----------------
Page No
------------------------------------------------------------------------------------------------------------------------
SIMULINK MODEL:
SUBSYSTEM DETAILS:
----------------------------------------------------------------------------------------------------------------------------- ----------------
Page No
------------------------------------------------------------------------------------------------------------------------
SUBSYSTEM1 DETAILS:
EXPECTED OUTPUT:
RESULT:
----------------------------------------------------------------------------------------------------------------------------- ----------------
Page No

Power Systems and Simulation Laboratory: Department of Electrical and Electronics Engineering

Uploaded by

Copyright:

Available Formats

Power Systems and Simulation Laboratory: Department of Electrical and Electronics Engineering

Uploaded by

Document Information

Original Description:

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Power Systems and Simulation Laboratory: Department of Electrical and Electronics Engineering

Uploaded by

Copyright:

Available Formats

POWER SYSTEMS AND SIMULATION

DHANEKULA INSTITUE OF ENGINEERING AND

JAWAHARLAL NEHRU TECHNOLOGICAL UNIVERSITY KAKINADA

Any of 5 experiments are to be conducted from each section:

1. Estimation of sequence impedances of 3-phase Transformer

Section II: Simulation Lab

DHANEKULA INSTITUTE OF ENGINEERING & TECHNOLOGY

Experiments as per curriculum

3 Estimation of ABCD parameters on transmission line model I R20C318.1

 Moulding Students as proficient, competent, and socially responsible engineering personnel

 To develop innovative, proficient Electrical engineers.

PROGRAM EDUCATIONAL OBJECTIVES(PEO’s)

PROGRAM SPECIFIC OUTCOMES (PSOs)

Course Outcomes vs POS Mapping:

Course Outcomes vs PSOs Mapping:

Course Out Comes PSO1 PSO2

DO’S AND DON’TS IN THE LAB

3. Correct specifications of the equipment have to be mentioned in the circuit diagram.

4. Student should be aware of operating the equipment.

5. Students should be at their concerned experiment table, unnecessary moment is restricted.

6. The readings must be shown to the Lecturer In-Charge for verification.

3. Don’t make or remove the connections with power ON.

4. Don’t switch ON the supply without verifying by the Staff Member.

5. Don’t switch OFF the machine with load.

Zero Sequence Impedance:

Experiment No: Date:

NAME PLATE DETAILS:

Experiment No: Date:

NAME PLATE DETAILS:

(b) L-L Fault:

(c) LL-G Fault:

Experiment No: Date:

For A and C: (figure-2)

For B and D: (figure-3)

Formulae: A = [VS/VR]IR=0; C = [IS/VR]IR=0;

Figure1: Power transfer capability of Transmission line without compensation

Experiment No: Date:

Hardware Connection diagram:

Figure1: Inductive Load connected at the receiving end

Figure2: Capacitive Load connected at the receiving end

Experiment No: Date:

Hardware Connection diagram:

With Capacitive load:

Figure1: Circuit for observing Ferranti Effect

Experiment No: Date:

Experiment No: Date:

APPARATUS REQUIRED: Personal Computer with MATLAB software

disp(‘Obtain the values of self admittances from self impedances’)

disp('Formulate Y-Bus Matrix')

disp(‘Ybus matrix is’)

Experiment No: 9 Date:

APPARATUS REQUIRED: Personal Computer with MATLAB software

5. Then click on the Run button, which executes the program.

disp('STEP1: Formulate Y-Bus Matrix')

disp('STEP2: Initialize bus voltages')

disp('STEP4: Calculate V2 value')

Experiment No: Date:

APPARATUS REQUIRED: Personal Computer with MATLAB software