International Journal of Engineering Research and Development

International Journal of Engineering Research and Development
e-ISSN: 2278-067X, p-ISSN: 2278-800X, www.ijerd.com
Volume 10, Issue 8 (August 2014), PP.35-50
Improving Transient Stability of Multi-Machine AC/DC Systems
via Energy-Function Method
1M Shobha, 2K.Durga Rao, 3Tegala.Srinivasa Rao,
1P.G Student Scholar, 2Assistant Professor, 3 Associate Professor
1,2,3Department of Electrical & Electronics Engineering
1,2,3Avanti Institute of Engineering & Technology, Makavaripalem(P),Vishakhapatnam
(Dt), Andhra Pradesh, India.
Abstract:- In this paper, the direct method of stability analysis using energy functions is applied for multi-machine
AC/DC power systems. The system loads including the terminal characteristics of the DC link are
represented as constant current type loads, and their effects on the generators at the internal nodes are obtained
as additional bus power injections using the method of distribution factors, thus avoiding transfer conductance
terms. Using the centre of angle formulation, a modified form of the energy-function method is used for the
swing equations and the DC link dynamical equations to compute the critical clearing time for a given fault.
Numerical results of critical clearing time for a single and multi-machine system using the energy-function
method agree well with the step-by-step method.
Index Terms:- Direct Current Link, Energy Function, External Control Signal.
I. INTRODUCTION
Direct methods for analyzing power system stability have been applied successfully so far for pure AC
systems. The literature on this topic is vast and was summarized in a survey paper by Fouad in 1975. Since
1975, there has been a significant advance in this research area, which has helped to remove the conservative
nature of the results associated with this method in the past. Hence, the possibility of using this technique for
transient security assessment is now quite good.
Since DC links are now being introduced for economic and other reasons, there is a need to extend the
direct method of stability analysis to systems containing such links. It is well known that the quick response of
the DC link, as opposed to an AC line combined with an effective control scheme, can enhance transient
stability. The degree of transient stability for given fault is either the critical clearing time or the critical energy.
The application of the direct method of stability analysis to AC/DC systems is not a routine extension of the
method as applied to AC systems. Instead, it requires a different approach based on treating the post fault DC
link dynamics as a parameter variation in the swing equations. A simplified first-order model of the DC link
controller is proposed, which augments the usual swing equation for the machine. For the multi-machine
case, the method using distribution factors is proposed to reflect (at the internal nodes) the terminal
characteristics of the DC link and the system loads as additional power injections. This eliminates automatically
the problem of transfer conductances in the swing equations. The computational algorithm and results for a
multi-machine system are presented.
II. SYSTEM DESCRIPTION
35
A. Multi-Machine AC/DC System
A 3-machine 9-bus system whose single-line diagram is shown in Fig. 3.1 is considered. For details of
the AC system date refer to [16]. A DC link is added to the system across Buses 9 and 4. The following
parameters are chosen for the DC link:

Improving Transient Stability of Multi-Machine AC/DC Systems via Energy-Function Method
퐾푎 = 1.0 pu/rad per sec, 푇푑푐 = 0.1 sec, 푃푟푒푓 = 0.0
for both prefault and postfault conditions. 푃푟푒푓 is assumed to be zero for the sake of convenience. In general,
푃푟푒푓 will have different values in the prefault and postfault states, as in the single-machine case; and AC/DC
load flow calculations have to be performed for each condition. Maximum 푃푑푐 = −2.0 pu; minimum; 푃푑푐 =
−2.0 pu; 푞푟 = 0.5. The external control signal (ECS) is chosen to be the difference between the rotor speeds of
the generator nearest to the rectifier and inverter terminals, i.e., 푢 = 휔3 − 휔1.
The extension of the method to multi-machine AC/DC systems involves a new method of handling system loads
and DC link characteristics in the swing equations, as well as the use of the potential energy boundary surface
method [3, 13] for computing 푉푐푟 .
TABLE I Comparison of 풕풄풓 by Energy-Function method and Actual simulation
Postfault Loading 푡푐푟 by Energy- 푡푐푟 by Actual on DC Link
Function method simulation
Case 1: Line (5, 7) 0.181 0.18
Case 2: Line (7, 8) 0.182 0.183
Case 3: Line (4, 6) 0.56 0.57
36
B. Representation of the Effect of Loads
It is well known that the transfer conductances present in the internal bus description using the classical
model pose a problem in constructing a valid V-function, as well as in computing 푡푐푟 . These transfer
conductances are mostly due to the system loads being converted to constant impedances and subsequent
elimination of the load buses. In the method proposed here, which also applies to the DC link element, the
effect of loads is reflected at the internal buses in the form of additional bus power injections.
Consider a power system network consisting of n buses and m generators. The bus admittance
matrix 푌퐿퐿 for the transmission network alone, excluding the loads and DC link, is formulated and is thus
augmented with the network elements corresponding to direct axis reactances of the m machines. The resulting
augmented matrix 푌퐵푢푠 has (n + m) buses altogether, and is represented as
(1)
where 푌퐺퐺 , 푌퐺퐿 , 푌퐿퐺 , and 푌퐿퐿 are submatrices of dimensions (m x m), (m x n), (n x m), and (n x n) respectively.
The overall network representation is
(2)

−1 (5)
37
where
퐼퐺푡
푡 = 퐼퐿1, 퐼퐿2 … . 퐼퐿푛
= 퐼퐺1, 퐼퐺2, … 퐼퐺푚 , 퐼퐿
퐼퐺 and 퐼퐿 are the current injections at the internal nodes of the generators and the transmission network nodes,
respectively; 퐸퐺 and 푉퐿 are the associated voltages. 푌퐵푈푆 is computed for the faulted and postfault conditions by
properly taking the corresponding network changes into account. The method of distribution factors are
suggested in [15] is now used for reflecting loads at the internal buses. Eliminating 푉퐿 from Eq. (3.2), we get
−1퐼퐿 − 푌퐿퐿
푉퐿 = 푌퐿퐿
−1푌퐿퐺 퐸퐺 (3)
and
퐼퐺 = 푌′ 퐸퐺 + 퐷퐿 퐼퐿 (4)
where
−1푌퐿퐺
푌′ = 푌퐺퐺 − 푌퐺퐿 푌퐿퐿
and the distribution factor matrix for loads is given by
퐷퐿 = 푌퐺퐿 푌퐿퐿
Also, we have
(6)
where 푃퐿푗 and 푄퐿푗 are the active and reactive power components of load at the jth bus. The additional bus
power injections at the internal bus of the kth generator (k = 1, 2….m) due to the load at jth bus (j = 1, 2 ….n) is
obtained as follows
(7)
Where 푑푘푗 is the appropriate (k, j) element of 퐷퐿 . The following assumption is made regarding the load
characteristics: the complex ratio of voltages
퐸푘
푉퐿푗
is assumed to be constant, corresponding to the prefault values. This is a deviation from the conventional type
of representation of loads as constant impedances. Since only active power is of interest in the swing equation,
we get
Δ푃푘퐿푗 = 푎푘퐿푗 푃퐿푗 − 푏푘퐿푗 푄퐿푗 (8)
The effect of all the loads at the internal bus of the kth generator is then obtained as
(9)
C. Representation of the Effect of DC Link
The effect of DC link is represented in a manner similar to that of the loads. For simplicity, we assume only one
DC link to be present. The analysis however, easily extends to cases of more than one DC link. In the 푌퐵푈푆 of
Eq. (1), all buses except the internal buses of the generators and the bus pair corresponding to the rectifier and
inverter terminals of the DC link are eliminated. The reduced network may be represented as

(10)
−1 (12)
38
where
퐼퐷푡
= 퐼푟 , 퐼푖 , 푉퐷푡
= 푉푟 , 푉푖
′ , 푌퐺퐷 , 푌퐷퐺 , 푌퐷퐷 are submatrices of
Subscripts r and I refer to the rectifier and inverter sides, respectively, and 푌퐺퐺
dimensions (m x m), (m x 2), (2 x m), and (2 x 2), respectively. From Eq. (10), we get
−1퐼퐷 − 푌퐷퐷
푉퐷 = 푌퐷퐷
−1푌퐷퐺 퐸퐺 (11)
and
퐼퐺 = 푌′′ 퐸퐺 + 퐷퐷 퐼퐷
where
푌′′ = 푌퐺퐺
−1푌퐷퐺
′ − 푌퐺퐷 푌퐷퐷
and the distribution factor matrix for the DC link is given by
퐷퐷 = 푌퐺퐷 푌퐷퐷
Now, we represent the effect of the DC link currents 퐼퐷 as additional bus power injections at the internal buses
of the generators. We have
and
(13)
where
푃푟 = −푃푖 = 푃푑푐
and
푄푟 = 푄푖 = 푄푑푐
It is assumed here that the DC link is lossless and the power factors at the rectifier and inverter stations are
equal. 푃푑푐 and 푄푑푐 are the active and reactive power components of the DC link that depend upon the DC link
controller dynamics. The effect of the rectifier and inverter ends of the DC link as additional bus power
injections at the internal bus of the generator is given by
(14)
(15)
where 푑푘푟 and 푑푘푖 are the appropriate (k, 1) and (k, 2) elements of the matrix 퐷퐷 .
From Eqs. (14) and (15), we get
Δ푃푘푟 = 푎푘푟 푃푑푐 − 푏푘푟 푄푑푐 (16)

39
and
Δ푃푘푖 = −푎푘푖 푃푑푐 − 푏푘푖 푄푑푐
As in the case of the load model representation, here also the ratios
퐸푘
푉푟
and
퐸푘
푉푖
are assumed to be
constant, corresponding to their prefault values. Since a simple structure is assumed for the DC link controller,
the output of which is Pdc, let Qdc = qr Pdc, where qr is a constant.
From Eq. (16), we get the total bus power injections at the kth generator due to the DC link as
Δ푃푘퐷 = Δ푃푘푟 + Δ푃푘푖 = 푎푘푟 − 푎푘푖 − 푞푟(푏푘푟 + 푏푘푖 )}푃푑푐
= 푐푘퐷 푃푑푐 푘 = 1,2 … . 푚 (17)
Where ckD is the expression in brackets in Eq. (17).
The parameters, (k = 1, 2…m, j = 1, 2…n),
which reflect the effect of the loads and the DC link, are thus computed for both the faulted and postfault
condition. By kron reduction technique, the bus admittance matrix 푌퐵푈푆 is reduced to the internal nodes of the
generators for these two conditions.
C. Inclusion of DC link Dynamics
A structure similar to that described earlier is assumed for the DC link controller whose equations in terms of
푃푑푐 are 푃푑푐 = − 1
푇푑푐
푃푑푐 +
푃푟푒푓
푇푑푐
+
퐾푎
푇푑푐
푢 (18)
Where u is the external control signal (ECS) obtained from the AC system quantities, such as the difference in
rotor speed of adjacent generators. The DC link dynamics are incorporated into the transient stability analysis in
manner similar to the approach described earlier. Also, 푃푑푐 is constrained to vary with in the specified practical
limits. While the faulted system equations are integrated, Eq. (18) also is solved for 푃푑푐 . At the end each time
step, the additional bus power injections at the internal buses of the generators are calculated using Eq. (17).
The effect of DC link is thus represented as the term that modifies the power input of the generator.
III. System Equations
Under the usual assumptions [1] for the classical model, and following notation in [2], the system equations in
the centre-of-angle reference frame are
푀푘휔 푘 = 푃푘 − 푃푒푘 −
푀푘
푀푇
푃퐶푂퐴
휃푘 = 휔 푘 k = 1, 2….m
(19)
where
푃푘 = 푃푚푘 − Δ푃푘퐿 − Δ푃푘퐷 − 퐸푘 2퐺푘푘
푃푒푘 = 퐶푘푗 sin 휃푗푘 + 퐷푘푗 cos 휃푘푗 푚푗
=1
≠푘
퐶푘푗 = 퐸푘 퐸푗 퐵푘푗 ; 퐷푘푗 = 퐸푘 퐸푗 퐺푘푗
and
휃푘 = 훿푘 − 훿표
Where 훿표 is the centre of angle defined by
푚푘
푀푇 훿표 = 푀푘 훿푘
푚푘
=1 , 푀푇 = =1
푀 푘
The following equation for 휔표 = 훿표
is easily derived

40
푚푘
푀푇휔표 = 푃푘
푚푗
=1 − 2 퐷푘푗 cos 휃푘푗
=푘+1
푚 −1
푘=1
≜ 푃퐶푂퐴 (20)
In our formulation, since the system loads are not converted into constant impedances, the transfer conductance
terms are only due to the transmission lines and hence can be neglected, i.e., 퐷푘푗 = 0. If the angle is constant,
i.e., 푃푘 = 0, then 푃퐶푂퐴 = 0 [13]. The postfault SEP is obtained by solving the set of nonlinear equations
푃푘 = 푃푒푘 k = 1, 2…..m-1 (21)
Solution of these power flow equations is discussed extensively in the literature [13, 14].
A. Transient Energy Function
The transient energy function used is that given in [2], assuming the damping to be zero
= kinetic energy (KE) + rotor potential energy
(PE) + magnetic potential energy (PE)
= 푉푘 휔 + 푉푝 휃 (22)
where
푉푝 휃 = rotor PE + magnetic PE
B. Computing 푉푐푟
Following [3], 푉푐푟 is computed as the value of 푉푝 along the sustained fault trajectory at the instant 푉푝 = 0. This
happens to be a point on the so-called potential energy boundary surface (PEBS) [13]. An assumption is made
that the PEBS crossing of the faulted trajectory is a good approximation to the value of 푉푐푟 , which is the value of
V(x) at the controlling UEP [3, 13].
C. Computational Algorithm
The algorithm for calculating the critical clearing time based on the proposed method is as follows:
1. Load flow calculation is performed for the prefault AC/DC system.
2. For the faulted and postfault states, the following computations are performed by augmenting the
passive network with generator reactance.
a. The overall 푌퐵푈푆 is computed excluding the loads and the DC link.
b. The distribution factors due to the system loads and DC link characteristics are computed as explained
earlier.
c. The 푌퐵푈푆 above is reduced to the internal buses of generators by eliminating all other buses. In
doing so, the transmission line resistance is neglected.
3. The postfault SEP is computed by solving the nonlinear Eq. (21).
4. The faulted Eqs. (18) & (19) are numerically integrated to obtain values of 휃, 휔 , 푃푑푐 at 푡 = Δ푡. At the
end of the integration interval, the following computations are done.
a. 푃푑푐 obtained from Eq. (18) is used in updating the bus power injections for both faulted and
postfault states.
b. 푃푘 is accordingly modified in Eq. (19), and the new postfault SEP 휃푠 is computed by solving Eq.
(21).
c. Using the updated values of 휃푠 , the V-function in (22), as well as 푉푝 and 푉푝 are calculated.

d. The integration is now continued for the faulted Eqs. (18) and (19) and steps (a) - (c) are
repeated at 푡 = 2Δ푡. For the sustained fault trajectory, this is continued until 푉푝 changes sign from positive to
negative. The value of 푉푝 at this instant is an estimate of 푉푐푟 .
5. Using this value of 푉푐푟 , the integration of the faulted equations is carried out and 푡푐푟 is reached when
푉 휃, 휔 = 푉푐푟 . Steps 4(a) and (b) are incorporated during the integration.
IV. SIMULATION RESULTS AND DISCUSSIONS
Case 1: Line (5, 7), t = 0.18sec, Unstable without HVDC
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
40
30
20
10
0
-10
-20
-30
Figure. 4.2. Variation of Energy with time without HVDC
Case1: line (5, 7), t = 0.18sec, Stable with HVDC
8
7
6
5
4
3
2
1
Figure. 4.3. Variation of Energy with time with HVDC
41
-40
Time (secs)
Energy in (p.u.)
W
pe
W
ke
W
T
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
0
Time (secs)
Energy in (p.u.)
W
pe
W
ke
W
T

0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
42
Case1: line (5, 7), t = 0.18sec
5
0
-5
-10
-15
-20
-25
-30
Figure. 4.4. Variation of potential energy with time
Case1: line (5, 7), t = 0.18sec
40
35
30
25
20
15
10
5
Figure. 4.5. Variation of kinetic energy with time
Case1: line (5, 7), t = 0.18sec
Figure. 4.6. Variation of Total energy with time
-35
Time (secs)
Potential Energy in (p.u.)
No HVDC
HVDC
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
0
Time (secs)
Kinetic Energy in (p.u.)
No HVDC
HVDC
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
25
20
15
10
5
0
-5
-10
Time (secs)
Total Energy in (p.u.)
No HVDC
HVDC

0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
43
Case1: line (5, 7), t = 0.181sec of Machine-1
Figure. 4.7. Variation of Rotor angle with time
100
0
-100
-200
-300
-400
-500
2500
2000
1500
1000
500
0
-600
Time (secs)
Rotor Angle in (degrees)
No HVDC
HVDC
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
-500
Time (secs)
No HVDC
HVDC
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
100
0
-100
-200
-300
-400
-500
-600
Time (secs)
No HVDC
HVDC

0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
40
30
20
10
0
-10
-20
-30
-40
Time (secs)
44
Case2: line (7, 8), t = 0.182sec, Stable with HVDC
8
6
4
2
0
-2
-4
-6
Case2: line (7, 8), t = 0.182sec
5
0
-5
-10
-15
-20
-25
-30
Figure. 4.12. Variation of potential energy with time
Energy in (p.u.)
W
pe
W
ke
W
T
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
-8
Time (secs)
Energy in (p.u.)
W
pe
W
ke
W
T
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
-35
Time (secs)
No HVDC
HVDC

0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
45
Case2: line (7, 8), t = 0.182sec
40
35
30
25
20
15
10
5
Case2: line (7, 8), t = 0.182sec
25
20
15
10
5
0
-5
0
Time (secs)
No HVDC
HVDC
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
-10
Time (secs)
No HVDC
HVDC
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
100
0
-100
-200
-300
-400
-500
-600
Time (secs)
No HVDC
HVDC

No HVDC
HVDC
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
Time (secs)
No HVDC
HVDC
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
Time (secs)
W
pe
W
ke
W
T
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
Time (secs)
46
2500
2000
1500
1000
500
0
-500
100
0
-100
-200
-300
-400
-500
-600
-700
20
10
0
-10
-20
-30
-40
-50
Energy in (p.u.)

0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
10
5
0
-5
-10
-15
-20
Time (secs)
47
Case 3: Line (4, 6), t = 0.56sec
5
0
-5
-10
-15
-20
-25
-30
-35
-40
Figure. 4.20. Variation of Potential energy with time
Case 3: Line (4, 6), t = 0.56sec
18
16
14
12
10
8
6
4
2
Energy in (p.u.)
W
pe
W
ke
W
T
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
-45
Time (secs)
No HVDC
HVDC
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
0
Time (secs)
No HVDC
HVDC

0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
48
Case 3: Line (4, 6), t = 0.56sec
10
5
0
-5
-10
-15
-20
-25
Case 3: Line (4, 6), t = 0.57sec of Machine-1
100
0
-100
-200
-300
-400
-500
-30
Time (secs)
No HVDC
HVDC
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
-600
Time (secs)
No HVDC
HVDC
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
1400
1200
1000
800
600
400
200
0
-200
Time (secs)
No HVDC
HVDC

No HVDC
HVDC
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
Time (secs)
49
V. CONCLUSION
1400
1200
1000
800
600
400
200
0
-200
In this paper, a technique is proposed for applying to the direct method of stability analysis to multi-machine
AC/DC systems. A new method of handling transfer conductances is presented that is also useful in
representing the DC link characteristics in the swing equations. The centre-of-angle formulation is used. A 3-
machine, nine-bus system illustrates the validity of the method and the effects of the DC link in improving
transient stability. Three case studies have been done and the computational algorithm and simulation results are
presented.
VI. APPENDIX
Generator data: Base 100MVA
Gen 1: 16.5/230 kv
Gen 2: 18/230 kv
Gen 3: 13.8/230 kv
Ka = 1.0 pu/rad per sec
Tdc = 0.1 sec
Pref = 0.0
Maximum Pdc = 2.0 pu
Minimum Pdc = -2.0 pu
qr = 0.5
REFERENCES
[1]. A.A. Fouad, “Stability Theory—Criteria for Transient Stability.” Proceedings of the Engineering
Foundation Conference on System Engineering on Power, Henniker, New Hampshire, Publication No.
CONF-750867, pp. 421-450, 1975.
[2]. T. Athay, R. Podmore, and S. Virmani, “A Practical Method for Direct Analysis of Transient
Stability.” IEEE Transactions on Power Apparatus and System, vol. PAS-98, No. 2, pp. 573-584,
1979.
[3]. N. Kakimoto, Y. Ohsawa, and M. Hayashi, “Transient Stability Analysis of Electric Power System via
Lure’ Type Lyapunov Function, Parts I and II.” Transactions of IEEE of Japan, vol. 98, No. 5/6,
May/June 1978.
[4]. M. Ribbens Pavella, Lj T. Gruijc, T. Sabatel, and A. Bouffioux, “Direct Methods of Stability Analysis
of Large Scale Power Systems.” Proceedings of the IFAC Symposium on Computer Applications in
Large-Scale Power Systems, New Delhi, India, vol. II, pp. 168-175, August 1979.

[5]. A. A. Fouad, and S. E. Stanton, “Transient Stability of a Machine Power System, Parts I and II.”
Papers No. 81 WM 078-5 and 81 WM 079-3, IEEE Winter Power Meeting, Atlanta, Georgia, Feb. 1-6,
1981.
[6]. E. Uhlmann, “Stabilization of an AC link by a Parallel DC link.” Direct Current, vol. 12, No. 8, pp.
50
89-94, 1964.
[7]. H. A. Peterson, and P. C. Krause, “Damping of Power Swings in a Parallel AC and DC system.” IEEE
Transactions on Power Apparatus and Systems, vol. PAS-85, pp. 1231-1239, 1966.
[8]. W. K. Marshall, K. R. Padiyar, L. M. Denton, W. J. Smolinski, and E. F. Hill, “ A Simplified HVDC
Link representation for Stability Studies.” Paper C-74-434-7, IEEE Summer Power Meeting, July
1974.
[9]. M. A, Pai, M. A. Mohan, and J. G. Rao, “Power System Transient Stability Region Using Popov’s
Method.” IEEE Transactions on Power Apparatus and Systems, vol. PAS-89, No. 5, pp. 788-794,
May/June 1970.
[10]. J. L. Willems, and J. C. Willems, “The Application of Lyapunov Methods to the Computation of
Stability Regions of Multi-Machine Power Systems, IEEE Transactions on Power Apparatus and
Systems, vol. Pas-89, No. 5, pp. 795-801, May/June 1970.
[11]. J. L. Abatut and J. Fantin, “Etude de la Stabilite Transitoire dans un Systeme de Puissance Influence
des Parametres.” Automatic Control Theory and Application, vol. 2, pp. 14-18, 1973.
[12]. H. Sasaki, “An Approximate Incorporation of Field Flux Decay into Transient Stability Analysis of
Multi-Machine Power Systems by the Second Method of Lyapunov.” IEEE Transactions on Power
Apparatus and Systems, vol. PAS-98, No. 2, pp. 473-483, March/ April 1979.
[13]. T. Athay, V. R. Sherkey, R. Podmore, S. Virmani, and C. Puech, “Transient Energy Stability
Analysis.” Systems Engineering for Power: Emergency Operating State Control- -Section IV, U.S.
Department of Energy Publication No. CONF-790904-PL, October 1979.
[14]. A. H. El-Abiad, and K. Nagappan, “Transient Stability Region of Multi-Machine Power
Systems.” IEEE Transactions on Power Apparatus and Systems, vol. PAS-85, pp. 169-179, February
1966.
[15]. W. F. Tinney, and W. L. Powell, “The REI Approach to Power Network Equivalents.” 1977 PICA
Conference, Toronto, Canada, May 1977.
[16]. P. M. Anderson, and A. A. Fouad, Power System Control and Stability. Ames, Iowa State University
Press, 1977.
[17]. J. J. Dougherty, and T. Hillesland, “Power System Stability Considerations with Dynamically
Responsive DC Transmission Lines.” IEEE Transactions on Power Apparatus and Systems, vol. PAS-
95, pp. 536-541, March/April 1976.
[18]. J. Reeve, “Multiterminal HVDC Power Systems.” IEEE Transactions on Power Apparatus and
Systems, vol. PAS-99, No. 2, pp. 729-737, March/April 1980.

International Journal of Engineering Research and Development

More Related Content

What's hot (17)

Viewers also liked (9)

Similar to International Journal of Engineering Research and Development (20)

More from IJERD Editor (20)

Recently uploaded (20)

International Journal of Engineering Research and Development