 Open Access
 Total Downloads : 63
 Authors : Rumana Ali, Anser Pasha C A
 Paper ID : IJERTCONV7IS08023
 Volume & Issue : RTESIT – 2019 (VOLUME 7 – ISSUE 08)
 Published (First Online): 13062019
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
Modelling and Analysis of a Multi Machine Power System using MATLAB/Simulink
Rumana Ali Anser Pasha C A
Assistant Professor Assistant Professor
Department of Mechatronics Department of Information Science & Engg. Mangalore Institute of Technology and Engineering, Moodbidri Yenepoya Institute of Technology, Moodbidri.
Abstract Transient stability is the ability of the power system to maintain synchronism when subjected to sever transient disturbance [1]. This article involves the modelling and investigation of the synchronous machine through the angular position of the rotor with respect to time after a fault occurs in the system. The article presents analysis of the test system which includes multimachine power systems. The transient stability analysis of the multi machine system is analysed with the help of MATLAB R2013a.
 INTRODUCTION
The transient stability studies involve the determination of whether or not synchronism is maintained after the machine has been subjected to severe disturbance. This may be sudden application of load, loss of generation, loss of large load, or a fault on the system. Transient stability is the ability of the power system to maintain synchronism when subjected to sever transient disturbance [1]. The response to this type of disturbance involves large excursions of rotor angles and is influenced by nonlinear powerangle relationship. Stability depends on the initial operating state of the system and the severity of the disturbance. The system usually altered after the disturbance which may cause the system to operate in a different steadystate status from that prior the disturbance.
transient stability study of a multimachine power system was developed using Simulink. It is basically a transfer function and block diagram representation of the system equations. The major objectives are:
 To determine whether the system is stable or unstable with a critical clearing time for a threephase fault.
 To analyze the effect of fault location and critical clearing time on the system stability, a threephase fault is located at two different locations.
 Disturbances for preventing the serious power
disruptions.
 Increasing value of H
 Increasing voltage set point and
 Decreasing mechanical input power.
.
 SYSTEM MODELING
The complete system has been represented in terms of Simulink blocks in a single integral model. It is self explanatory with the mathematical model given below. One of the most important features of a model in Simulink is its tremendous interactive capacity. It makes the display of a signal at any point readily available; all one has to do is to add a Scope block or, alternatively, an output port. Giving a feedback signal is also as easy as drawing a line. A parameter within any block can be controlled from a MatLab command line or through an mfile program. This is particularly useful for a transient stability study as the power system configurations differ before, during and after fault.
 Changing Generator Inertia Constants, H
Inertia constant H is defined as the ratio of kinetic energy at rated speed to rated apparent power of the machine.
H=stored energy in megajoules
rating in MVA
(1)
Fig 1. WSCC 3machine, 9bus system
Power systems are designed to be stable for a selected set of contingencies. The contingencies usually considered are short circuits of different types: phasetoground, phaseto phaseto ground, or threephase [2]. Fig.1 depicts the standard 9Bus system considered for the analysis. A complete model for
Inertia constant H simply quantifies the kinetic energy of the rotor at the synchronous speed in terms of the number of seconds it would take for the generator to provide an equivalent amount of electrical energy when operating at a power output equal to its MVA rating [6].With the transmission line open, and with D = 0, the swing equation becomes:
(t)=fo P0 (2)
H M
The machine inertia constant H plays an important role in stability assessment. By increasing H the system stability can be improved. From the equation perspective, equation (2) indicates that enhancing H would reduce the rotor acceleration and therefore help in the stability of the machine.
 Generator Voltage Set points
2
2
IEaIIVI sin + V ( 1 1
) sin 2 (3)
XL+Xd
2 XL+Xq
XL+Xd
Increasing generator voltage set point enhances the system stability. Equation (3) shows that decreasing Ea would lower the generator real power output curve and therefore reduce the deceleration area in the equal area criterion diagram.
Fig 2: Angular positions of the individual generators
 Generator Real Power Output
Increasing generator real power output enhances the system transient stability. In steady state condition.
M
M
PG (o) = P0 (4)
M
M
M
M
As a result, setting up generator real power output in prefault period settle the value for P0 .From the equal area criterion it is noted that increasing P0 would decrease the deceleration area and reduce the stability. To improve the power system stability after disturbance, it is required to either decrease the fault clearing time, mechanical input power, or increase the generator inertia constants H, terminal voltage, and armature voltage.
 Changing Generator Inertia Constants, H
 SIMULATION RESULTS.
The system is simulated using power system tool box of MATLAB with a fault clearing time of 0.1s. The angular positions, relative angular positions and relative angular velocities of generators is depicted in Fig 24. Fig.5 shows the accelerating powers of the generators Table 1 shows the summary of the transient stability analysis of IEEE 14bus power system under three phase fault at selected buses and lines.
Fig.3: Relative angular positions of 21
and 31
Fig 4: Relative angular velocities of 12 and 13
Fig 5: Generator accelerating powers
Table 1: Summary of IEEE 14Bus Fault Analysis
Sl. No Faulte d Bus No Remove d Faulty Line Critical Clearin g Time Clear ing Time Remarks 1 11 611 0.33 0.3 Stable 2 11 611 0.33 0.33 Unstable Fig.6 shows the response of the machine for change in inertia constant from 5.148 to 3, it is observed in that the response improve as the inertia constant is increased from 5.148 to 7.
(a)
(b)
Fig 6: Response with Change in Inertia Constant H (a)Decrease in H (b) Increase in H
Fig. 7 shows the response with change in Generator Voltage Set points. It is observed that the system gives significant stability response with generator 1s voltage set points from 1.0 to 0.95.
(a)
(b)
Fig 7 : Responses with Change in Voltage Set Point (a)Decrease in Set Point (b) Increase in Set Point
Fig. 8 shows the responses with real power output. .It is noted that the system remains stable and has better performance when its subjected to higher real power output.
(a)
(b)
Fig 8: Responses with Change in Real Power output (a)Decrease in Real Power output (b) Increase in Real Power output
 CONCLUSIONS
nonlinear differential equations namely the swing equation. The modeling and analysis of a multimachine system greatly improves the stability of a system.
REFERENCES

 P. Kundur, Power System Stability and Control, EPRI Power System Engineering Series.(Mc GrawHill, New York, 1994).
 I. J. Nagrath and D. P. Kothari, Power System Engineering (Tata McGrawHill, New Delhi, 1994).
 W. Long et al., EMTP a powerful tool for analyzing power system transients, IEEE Computer. Appl.Power, 3 (July 1990), pp 3641.
 L. W. Nagel, SPICE 2 A computer program to simulate semiconductor circuits, University of California, Berkeley, Memo. ERLM520, 1975.
 LouisA Dessaint et al., Power system simulation tool based on Simulink, IEEE Trans. Industrial Electronics, 46 (6) (1999), pp.12521254.
 Simulink Users Guide (The Math works, Natick, MA, 1999).
 Hadi Sadat, Power System Analysis (McGrawHill, New York, 1999).
 Power System Block set Users Guide (The Math works, Natick, MA, 1998).
 A. Arapostathis, 5. Sastry, and P. Varaiya, Analysis of powerflow equation, Electrical Energy Power Syst., vol. 3, no. 3, pp.
 T. Athay, V. R. Sherket, R. Podmore, S. Virmani, and C. Puech, Transient energy stability analysis, in Proc. Conf. on System Engineering for Power: Emergency Operating State Control Section /V (Davos, Switzerland, 1979).
 T. Athay, R. Posmore, and S. Virmani, A practical method for direct analysis of transient stability, /E Trans. Power App. Syst., vol. PAS%, no. 2, pp. 573584, Mar./Apr. 1979.
 P. D. Aylett, The energyintegral criterion of transient stability limits of power systems, Proc. Inst. Elec. fng. (London),
 1. Bailleul and C. I. Byrnes, Geometric critical point analysis of lossless power system models, /E Trans. Circuits Syst., vol. CAS29, no. 11, pp. 724737, Nov. 1982.
 Berggren, B. and Andersson, G. (1993). On the nature of unstable equilibrium points in power systems. Power Systems, IEEE Transactions on, 8(2), pp738745.