The School of Engineering and Physics: Experiment 1
The School of Engineering and Physics: Experiment 1
Date: 14/02/2020
GROUP MEMBERS:
Nischal Kavitesh Kumar– S11146567
0
AIM
This lab is aimed at teaching students on how to use:
INTRODUCTION
Matlab software is widely used in engineering as it has many features that can support
calculations, simulation and using scripting tools which is very important for this lab. Also
Matlab is also very useful when it comes to analyzing power systems, with proper codes, delta
connected or wye connected connections can both be implemented on Matlab. Other than just
the coding part of the Matlab, it also has additional features add-on it like Simulink which is not
a written programming but a graphical programming language which is needed in order to do
simulation, analyzing and modeling. The Simulink has many built in tools which is located in the
SimPowerSystems toolbox which has access to all the tools required to carry out simulation. [1]
RESULTS
Prelab
I)
wt = 0:0.05:2*pi; % describing range
Va = 1200*cos(wt); % declaring the voltageVa
Vb = 1200*cos(wt-(120*pi/180)); % declaring the voltage Vband phase
shift of 120 degrees
Vc = 1200*cos(wt+(120*pi/180));% declaring the voltage Vc and phase
shift of 120 degrees
Ia = 10*cos(wt-(36.87*pi/180)); % declaring the current Ia (I= V/Z)
1
Ib = 10*cos(wt-(156.87*pi/180)); % declaring the current Ia (I= V/Z)Ic
= 10*cos(wt+(83.13*pi/180)); % declaring the current Ia (I= V/Z)
Pa = Va.*Ia; % finding power Pa
Pb = Vb.*Ib; % finding power Pb
Pc = Vc.*Ic; % finding power Pc
Psum= Pa + Pb +Pc; %finding total power
plot(wt,Pa,'r',wt,Pb,'b',wt,Pc,'g',wt,Psum,'y') % graphing the power
sinusoids
ylabel('Instantaneous Powers') % y axis
P = mean(Pa+Pb+Pc); % finding the average of the 3 powers
legend('Pa','Pb','Pc','Psum') % labeling the legends
xlabel('wt') % x axes
title('Graph of Power Vs wt') % title
Ptotal=(3/2)*1200*10*cos(36.87*pi/180); % value of the total 3 phase
power
disp('Total Average Power (graphically) = ')
disp(P)
disp('Total Average Power (using the Formula 3*|V|*|I|*cos(phi) = ')
disp(Ptotal)
II)
2
t = 0:0.0001:0.0167;
V =200* cos(377*t);% declaring volatge value
P = 800+1000*cos(754*t-36.87*pi/180);% declaring power value
I = P./V; %declaring current
wt = 180/pi*377*t; % declaring the values of wt.
figure
subplot(2,2,1)
plot(wt,V),grid %graph of voltage
xlabel('wt') % x axes label
ylabel('Voltage') % y axes label
title('V(t)') % title label
subplot(2,2,2)
plot(wt,P),grid % graph of power
xlabel('wt') %x axes label
ylabel('Power') % y axes label
title('P(t)') %title label
subplot(2,2,3)
plot(wt,I),grid % graph of current
ylabel('current')%y axes label
title('I(t)') %title label
3
The circuit from the lab manual was made on Simulink and generated with
frequencies of firstly 50 Hz and then 60 Hz.
4
Output Graphs
Figure 4: The output current and voltage graph when the simulation is done at a frequency of 60Hz.
Figure 5: The output current and voltage graph when the simulation is done at a frequency of 50Hz.
5
Figure 6: Graph of Active power = 320W & Reactive power = 0Var.
DISCUSSION
In this experiment, we were introduced on how to use Simulink on Matlab and to analyze
power systems. For part 1 of the lab, we aimed to analyze power systems by using Matlab. At
first, a balanced 3 phase source was used on Matlab to plot the instantaneous power for each
phase with respect to wt. Figure 1 shows the graph of the 3 phase individual power against wt.
When all these instantaneous powers are added we get total 3 phase real power, which if
compared with the calculations done analytically are showing the same result. Secondly, a
single phase load is supplied with sinusoidal voltage (v(t) = 200cos (377t)) and by using Matlab,
the voltage, power and current were plotted over a range of 0 to 16.67. From observing the
current flow in figure 2 and comparing with analytical calculation, it is clear that the peak
amplitude, phase angle and the angular frequency are the same.
Adding on, the power system figure from the lab handout was implemented on Simulink by
using SimPowerSystem toolbox, specifically electrical tools. From looking at figure 3, it is clearly
shown that the respective currents, voltage and phasors angles for the 60Hz system. When
respective formulas were used we were able to find the total active and total reactive power
when supplied and consumed by the system. For the purpose of getting correct waveform we
changed simulation time from 2 seconds to 0.1 seconds. Moving on, as for the 50Hz system, the
same steps were done. Figure 4 shows the currents, voltage and phasors angles.
6
CONCLUSION
To conclude, it could be said that the two main aims of the lab was achieved successfully. The
power analysis was deeply understood by the use of Matlab where the important power
systems were found by using specific codes. Both the three phases and a single phase load were
found by Matlab. The second aim of the lab was also met as it once learned, the Simulink was
easy to begin with because the toolbox was to get the required power outcome which taught
us what to expect as output of the Simulink.
REFERENCE
[1]EE321 Power System Analysis: Lab 1 Manual - Simulation and Analysis of Power Systems
Using Matlab/ Simulink and PSIM
[2] Halvorsen.blog, 2020. [Online]. Available:
https://www.halvorsen.blog/documents/tutorials/resources/simulink/Introduction%20to
%20Simulink.pdf. [Accessed: 20- Feb- 2020].