UNIVERSITIMALAYSIAPERLIS
FACULTY OF ELECTRICAL ENGINEERING
TECHNOLOGY
EMJ32004 POWER SYSTEMENGINEERING
LABORATORY MODULE2
TRANSMISSION LINEMODEL
Name: MOHAMOUD ADEN OSMAN
Matrix No: 191090039-1
Programme: RK 23
Group: G1B
EMJ32004 Power System Engineering Laboratory Module 2
TASK 1: PRE – LAB CALCULATION (40 Marks)
(Please use your own paper if the space is not enough.)
1. A 345 kV, three phase transmission line is 130 km long. The resistance per phase is 0.036
Ω/km and the inductance per phase is 0.8 mH/km. The shunt capacitance is 0.0112 μF/km.
The receiving end load is 270 MVA with 0.8 power factor lagging at 325 kV. Use the
medium line model to find the voltage and power at the sending end and the voltage
regulation/
2
Faculty of Electrical Engineering Technology, Universiti Malaysia Perlis (UniMAP)
2. A 345 kV, three phase transmission line is 130 km long. The series impedance is z = 0.036
+ j0.3 Ω per phase per km, and the shunt admittance is y = j4.22 x 10-6 siemens per phase
per km. The sending end voltage is 345 kV and the sending end current is 400 A at 0.95
power factor lagging. Use the medium line model to find the voltage, current and power at
the receiving end and voltage regulation.
1
Faculty of Electrical Engineering Technology, Universiti Malaysia Perlis (UniMAP)
TASK 2: WRITING MATLAB COMMAND
A 345 kV, three phase transmission line is 130 km long. The resistance per phase is 0.036 Ω/km
and the inductance per phase is 0.8 mH/km. The shunt capacitance is 0.0112 μF/km. The
receiving end load is 270 MVA with 0.8 power factor lagging at 325 kV. Use the medium line
model to find the voltage and power at the sending end and the voltage regulation.
1. Type the following command at M-File window and save it in your own folder
(e.g.: exp2.m)
r = 0.036; g = 0; f = 60;
L = 0.8; % milli-Henry
C = 0.0112; % micro-Farad
Length = 130; VR3ph = 325;
VR = VR3ph/sqrt(3) + j*0; % kV(receiving end phase voltage)
[Z, Y, ABCD] = rlc2abcd(r, L, C, g, f, Length);
AR = acos(0.8);
SR = 270*(cos(AR) + j*sin(AR)); % MVA (receiving end power)
IR = conj(SR)/(3*conj(VR)); % KA (receiving end power)
VsIs = ABCD* [VR; IR]; % column vector [Vs; Is]
Vs = VsIs(1);
Vs3ph = sqrt(3)*abs(Vs); % KA (sending end L-L voltage)
Is = VsIs(2); Ism = 1000*abs(Is); % A (sending end current)
pfs= cos(angle(Vs)- angle(Is)); % A (sending end power factor)
Ss = 3*Vs*conj(Is); % MVA (sending end power)
REG = (Vs3ph/abs(ABCD(1,1)) - VR3ph)/VR3ph *100;
fprintf(' Is = %g A', Ism), fprintf(' pf = %g\n', pfs)
fprintf(' Vs = %g L-L kV\n', Vs3ph)
fprintf(' Ps = %g MW', real(Ss)),
fprintf(' Qs = %g Mvar\n', imag(Ss))
fprintf(' Percent voltage Reg. = %g\n', REG)
1
Faculty of Electrical Engineering Technology, Universiti Malaysia Perlis (UniMAP)
RESULTS TASK 2:
Table 1: Comparison of pre-lab calculation and simulation results for Task 2.
Parameters Pre-Lab Calculation Simulation results
Z 4.68 + j39.26 Ω 4.68 + j 39.2071 ohms
Y 0+j0.549m S 0 + j 0.000548899 Siemens
A 0.989+j1.285m 0.98924+ j 0.0012844
B 4.68 + j39.26 4.68 + j 39.207
C -0.352u + j0.548m -3.5251e-07 + j 0.00054595
D 0.989+j1.285m 0.98924 + j 0.0012844
Current at the sending end (Is) 420.840 A 421.132 A
Voltage at the sending end (Vs) 348.790k V 345.002 L-L k V
Real power at the sending end (Ps) 218.768M W 218.851 M W
Reactive power at the sending end 123.961 Mvar 124.23 Mvar
(Qs)
Power Factor (PF) 0.8 0.869657
% Voltage regulation (%VR) 7.32% 7.3091%
1
Faculty of Electrical Engineering Technology, Universiti Malaysia Perlis (UniMAP)
TASK 3:
Write the MATLAB command for the problem below. Use the program function
zy2abcd.m to execute your program. Verify your answer by comparing it with your pre-
lab calculation.
A 345 kV, three phase transmission line is 130 km long. The series impedance isz
= 0.036 + j0.3 Ω per phase per km, and the shunt admittance is y = j4.22 x 10-6 siemens per
phase per km. The sending end voltage is 345 kV and the sending end current is 400 A at 0.95
power factor lagging. Use the medium line model to find the voltage, current and power at the
receiving end and voltage regulation.
Show your complete programming and the result in the provided section.
1
Faculty of Electrical Engineering Technology, Universiti Malaysia Perlis (UniMAP)
EMJ32004 Power System Engineering Laboratory Module 2
RESULTS TASK 3:
z = 0.036 + j*0.3;
y = j * 4.22/1000000;
Length = 130;
Vs3ph = 345;
Ism = 0.4;
As = -acos (0.95);
Vs = Vs3ph/sqrt (3) + j*0;
Is = Ism * (cos(As) + j*sin (As));
[Z,Y,ABCD] = zy2abcd (z,y,Length);
VrIr = inv (ABCD)*[Vs;Is];
Vr = VrIr(1);
Vr3ph = sqrt(3)*abs(Vr);
Ir = VrIr(2);
Irm = 1000*abs(Ir);
pfr = cos (angle(Vr)- angle(Ir));
Sr = 3* Vr* conj (Ir);
REG = (Vs3ph /abs(ABCD(1,1)) - Vr3ph)/Vr3ph*100;
fprintf('Ir = %gA',Irm), fprintf ('pf = %g\n', pfr)
fprintf('Vr = %g L-L kV\n', Vr3ph)
fprintf(' Pr = %g MW', real(Sr)),
fprintf(' Qr = %g Mvar\n', imag(Sr))
fprintf(' Percent voltage Reg. = %g\n', REG)
1
Faculty of Electrical Engineering Technology, Universiti Malaysia Perlis (UniMAP)
EMJ32004 Power System Engineering Laboratory Module 2
Table 2: Comparison of pre-lab calculation and simulation results for Task 3.
Parameters Pre-Lab Calculation Simulation results
Z 4.68+j39 Ω 4.68 + j 39 Ω
Y 0+j0.549m S 0 + j 0.0005486 S
A 0.989+j1.268m 0.9893 + j 0.0012837
B 4.67+j39 4.68 + j 39
C -0.352u+j0.546m -3.5213e-07 + j 0.00054567
D 0.989+j1.268m 0.9893 + j 0.0012837
Current at the receiving end (Ir) 442.198 A 441.832A
Voltage at the receiving end (Vr) 330.115 kV 330.68 L-L kV
Real power at the receiving end (Pr) 224.904 MW 224.592 MW
Reactive power at the receiving end 116.978 Mvar 116.612 Mvar
(Qr)
Power Factor (PF) 0.95 0.887501
% Voltage regulation (%VR) 5.67% 5.45863%
1
Faculty of Electrical Engineering Technology, Universiti Malaysia Perlis (UniMAP)