Energy and Power Engineering
Vol.10 No.06(2018), Article ID:85654,18 pages

Transient Stability Simulation of 33 kV Power Grid

Nader Barsoum, Nadira Zulkeffley, Myjessie Songkin

Department of Electrical and Electronics Engineering, Faculty of Engineering, Universiti Malaysia Sabah (UMS), Kota Kinabalu Sabah, Malaysia

Copyright © 2018 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

Received: June 2, 2018; Accepted: June 25, 2018; Published: June 28, 2018


An example of Sandakan power grid problem is presented in this paper. Sandakan is a suburb in east coast of Sabah state of Malaysia. Stability problem occurs due to the increase in load demand, lack of generation sources and inadequate supply. The tripping disturbances occur frequently in the network which is contributing to voltage instability. In this paper dynamics stability of 33 kV power grid as related to Sandakan network is analyzed and simulated. The analysis is completed by modelling the network data in Power System Simulation for Engineering (PSS/E) software and simulate the transient stability of generator, exciter and governor during a three phase fault occurs on a far and close distance from a bus, and determine the critical clearing time as well as swing curve of rotor angle. The output values of electrical power, machine speed, rotor angle and bus voltage are observed.


Transient Stability, PSS/E, Load Flow, Generator, Exciter, Governor, Critical Clearing Time

1. Introduction

Generally, a power system under normal operating conditions may face a contingency such as transmission element outages, generator outages, loss of transformer, and sudden change in the load or faults [1] . Transient is an event occurs when a power system subjected to large disturbances under dynamic stability [2] . Disturbances include loss of synchronism, loss of generation, loss of load in transformers or faults on transmission element and lines. Transient stability is one of major analysis in the power system in order to ensure the system stability to withstand a major disturbance and to ensure that the transmission system is operated safely, steady state and contingency analysis must be performed [7] .

The round rotor generator model (GENROU) represents as solid rotor generator at sub-transient level is used to produce machine rotor angle for transient stability. Rotor angle has the ability of interconnecting the synchronous machines with power system to remain in synchronism. Stability is related to generator electromagnetic torque and mechanical torque which cause the rotor to accelerate or decelerate [1] . Voltage stability is related to change in the load. This stability is the ability of a power system to maintain steady voltages at all busses from a given initial condition after being subjected to a disturbance [6] . Voltage instability increases by load demand or change in system condition which cause the incontrollable drop in voltage. With abnormal low voltage it is lead to voltage collapse also contributes to blackout of the grid system [2] . The problems were reported with power flow and contingency in terms of blackout after main grid supply outages with overload and high fault current on distribution system [4] [6] . Critical clearing time is known as maximum time duration that a fault may occurs in a power system without loss of stability. There are three type of fault condition which pre-fault system conditions, fault structure (type and location) and post fault conditions [3] . The three-phase fault is the most serious kind of fault and its critical clearing time can reflect the transient stability of power system. Critical clearing time (CCT) can be obtained by trial and error method [4] .

Synchronous generator is the source of electrical energy where in the generator the mechanical energy usually transformed into electrical energy. This transformation is provided by excitation of synchronous generator and is regulated by excitation system. An IEEE Type 2 Excitation System (IEEET2) and turbine governor such as Gas Turbine (GAST) and Turbine IEEE Type 1 Speed Governing Model (IEEEG1) are used.

Generator excitation is defined as generator output voltage and output reactive power. It means that the excitation is actually output energy of generator regulation and this can impact the stability of the power system. The use of an excitation system is for maintaining the output voltage, control the shaft’s speed and enhancing the generator performance.

In this paper, PSS/E will be used for characterizes the power system transmission network and generation performance for both load flow analysis and transient analysis [5] . All the sources referred from the books, articles, research papers and journals.

2. Modelling the Network in PSS/E

The following information needed for modelling the network. This includes bus data, branch data, load data, generator data and transformer data. These data are saved in data.sav file in PSS/E and are given in the Appendix. All data and parameters are taken from Sandakan power grid system.

The power network consists of 26 buses, 8 generators and 22 loads. The highest bus base is 33 kV and the lowest is 6.6 kV. Figure 1 represents the existing network.

3. Transient Stability Analysis

In order to achieve stability of power system, load flow study is an important tool that gives a numerical solution. In power flow analysis a per-unit system is used for voltage magnitudes and angles, real and reactive powers.

In conducting transient analysis, there are three machine model must be taken into account such as generator, exciter and governor. These are as follows:

- Round Rotor Generator Model (GENROU)

- Exciter IEEE Type 2 Excitation System (IEEET2)

- Turbine Governor GAST

- Turbine IEEE Type 1 Speed Governing Model (IEEEG1)

Tables 1-4 show the parameters of the above models

4. Load Flow Result

The diagram Figure 1 is the operation of load flow for 26 buses divided into 1 slack bus which is BN_11 as a swing bus. The transient stability analysis requires the solution of a system of coupled non-linear differential equation. Load flow simulation on PSS/E using Newton Raphson gives the following values in Tables 5-7.

Figure 1. Sandakan power grid model in PSS/E.

Table 1. Parameters of genrou genrator model.

Table 2. Parameters of IEEET2 exciter model.

Table 3. Parameters of GAST turbine model.

Table 4. Parameters of IEEEG1 governor model.

Table 5. Swing bus summary.

Table 6. Voltage performance under normal conditions (Pre-Disturbance).

Table 7. Voltage performance under contingency conditions (Post-Fault).

Table 5 shows the swing bus BN_11 holds the highest real power generation, Pgen is 27.8 MW and reactive power generation, Qgen is 7.4 MVar. This swing bus is a special generator bus serving as the reference bus. The steady-state supply voltage limits applicable for the pre-disturbance and post-fault state defined in Table 6 and Table 7. These variations of voltages are normally applied for pre and post fault

5. Bus Voltage before and after Load Flow Simulation

Table 8 shows all the voltages before load flow analysis are under normal condition which are not exceed than 105% (overvoltage) and not below 95% (under-voltage). Thus, the grid does not have any critical busses.

These parameters are found from Sandakan power grid of 26 buses. Load flow increases the voltage bus as shown.

6. Transient Stability Result

The transient stability analysis approach by applying a fault on a bus and run the simulation from time = 1 until time = breaker open. Based on the Sandakan grid’s simulation in PSS/E, it can be seen from Figure 2 that all the machines in the system are in good initial condition.

The data from channel output file in PSS/E of Sandakan grid such as swing curve of rotor angle, impulse response of shaft speed, electrical power and bus voltage are plotted in Figures 3-9.

6.1. Generator

In the normal state, the machines in a power system network operate at equilibrium corresponding to the mechanical power input, Pm being equal to the electrical power output, Pe. When a fault occurs in the system at time = 1.0 seconds, the mechanical power input become greater than the electrical power output (Pm > Pe), then the speed of the machines increase as it will accelerate the rotor.

Table 8. Voltage performance before and after load flow.

Figure 2. Initial condition of machines.

Figure 3. Swing rotor angle vs time for generators.

Figure 4. Impulse response of speed vs time for machines.

Figure 5. Active power vs time for machines.

Figure 6. Terminal voltage vs time for exciters.

Figure 7. Step response of speed vs time for governors.

Figure 8. Bus near and far from generator.

Figure 9. Voltage vs Time when faulted at bus BS_33.

With the fault cleared at time = 1.18 seconds, the electrical power becomes low but the rotor still running above synchronous speed, hence the angle and the electrical power continue to increase. When the electrical power is greater than mechanical power (Pe > Pm), it causing the rotor to decelerate toward synchronous speed until the angle reaches its critical value. When the system reached the critical value, the rotor angle will continue to oscillate back and forth at its natural frequency until it becomes stable.

In running the simulation the time of fault clear is from breaker open until 10 seconds.

6.2. Exciter

The results of the excitation system response in Figure 6, steady-state value for LD1_11, LD2_11, LD3_11 and LD4_11 are 0.994996 pu (terminal voltage). Steady-state value for GN_11 is 0.662346 pu (terminal voltage). Steady-state value for BN_11 is 0.455770 pu (terminal voltage). Steady-state value for KB_6.6 and SB_6.6 are 0.401253343 pu (terminal voltage).

Based on the results in Figure 6 indicate that GN_11, BN_11, KB_6.6 and SB_6.6 exciters perform slower than LD1_11, LD2_11, LD3_11 and LD4_11 exciters which are this exciter meet the 1 pu the excitation system voltage requirement. These bus conditions are shown in Table 9.

6.3. Governor

Figure 7 shows the variation of speed with time for governor IEEEG1 and GAST types. All final frequencies were determined by the droop, R of the responding governors, Table 10. The frequency drops depends upon the generator inertia values. The least frequency deviation occurs with high inertia and fast governors. Governor condition are given in Table 11.

6.4. Critical Clearing Time (Result)

To determine the critical clearing time when the fault occurs at bus that close and far from generator, refer to Figure 8. It shows that bus BS_33 is near to generator BN_11, while for bus SA_33 is 14 kilometers far from generator GN_11.

Figure 9 and Figure 10 have been obtained by the technique of trial and error method in order to determining the critical clearing time of the system where the fault duration was increased gradually using the step time of Δt = 0.01 seconds until the system appears to be unstable by observing machine’s rotor angle as a reference point.

Figure 9 shows the variation of voltage with time for a three phase fault applied on bus BS_33 (bus near to generator). Since the three phase fault applied at time = 1.0 seconds, then the fault is clear at time = 1.18 seconds. Thus, the critical clearing time is 0.18 seconds and the system becomes stable.

Figure 10 shows the variation of voltage with time for a three phase fault applied on bus SA_33 (bus far from generator). Since the three phase fault applied at time = 1.0 seconds, then the fault is clear at time = 1.23 seconds. Thus, the critical clearing time is 0.23 seconds and the system becomes stable.

It can be seen that transient stability is greatly affected by the location of a fault from bus to generator. Table 12 shows the critical clearing times in seconds determined for all the twenty-six buses on the Sandakan Power Grid.

Transient stability analysis is run starting with a clearing time of 0.01 seconds. If the system is proved a stable condition, another analysis run is made by increasing the clearing time higher than first run. If the second run is still in stable condition, then more runs are made until the system becomes unstable. If the run showed an unstable system, then the clearing time of previous run gives the desired result.

Figure 10. Voltage vs Time when faulted at bus SA_33.

Table 9. Exciter conditions.

Table 10. Steady-state of governor.

Table 11. Governor conditions.

Table 12. Critical clearing time with different location.

7. Conclusion

This paper presents a modeling and simulating case data of 26 buses, 8 generators and 22 loads 33 kV power grid that in service mode. Through load flow analysis, the voltage performance under different conditions can be determined. The desired analysis of the transient stability of the system based on output values such as machine rotor angle, electrical power, machine speed and bus voltage were found to be stable after fault is cleared. The rotor angle, power, speed and the voltage in the grid system is back to its normal condition where there is no generator set that will out of phase and when the fault is cleared. The theory of the critical clearing time (CCT) when the fault occurs close and far from the generator was proved in this paper by using trial and error method. Thus, as the distance between the bus and the generator increases, the critical clearing time also increases.

Cite this paper

Barsoum, N., Zulkeffley, N. and Songkin, M. (2018) Transient Stability Simulation of 33 kV Power Grid. Energy and Power Engineering, 10, 301-318.


  1. 1. Kundur, P. (1994). Power System Stability and Control. In: Lauby, N.J., Ed., Methods of Improving Stability, McGraw-Hill Inc., p. 1104.

  2. 2. Kundur, P., et al. (2004) Definition and Classification of Power System Stability. IEEE Transactions on Power Systems, 19, 1387-1401.

  3. 3. Hooda, P.R. (2012) Transient Stability Analysis of Power System Using Matlab. International Journal of Engineering Sciences and Research Technology (IJESRT).

  4. 4. Ren, L.-Y., Fang, T., Yan, J.-F., Yu, Z.-H., Su, F. and Wu, T. (2012) Online Application and Fast Solving Method for Critical Clearing Time of Three-Phase Short Circuit in Power System. International Journal of Smart Grid and Clean Energy, 2, 93-99.

  5. 5. Mohamad, A.M., Hashim, N., Hamzah, N., Ismail, N.F.N. and Latip, M.F.A. (2011) Transient Stability Analysis on Sarawak’s Grid using Power System Simulator for Engineering (PSS/E). IEEE Symposium on Industrial Electronics and Application (IEA), Langkawi, 25-28 September 2011, 521-526.

  6. 6. Barsoum, N., Asok, C.B., Kwong, D.T.S. and Kit, C.G.T. (2017) Power Analysis for a Limited Bus Grid System with Distribution Generators. Global Journal of Technology and Optimization, OMICS Online, 8, 1-10.

  7. 7. Barsoum, N. and Lee, C.Z. (2013) Simulation of Power Flow and Protection of a Limited Bus Grid System with Injected Solar Power. Journal of Energy and Power Engineering EPE, Scientific Research on Optimization, 5, 59-69.


Table I. Buses data used in PSS/E.

Table II. Transformer branch data used in PSS/E.

Table III. Machines data used in PSS/E.

Table IV. Load data used in PSS/E.

Table V. Branch/Distribution line data used in PSS/E.