Power System Stability Classical Model
Power System Stability Classical Model
Power System Stability Classical Model
Assumptions
1 Exciter dynamics are neglected and the filed current is assumed to be
constant so that the induced voltage is always constant.
2 Damper winding dynamics are neglected.
3 The mechanical input power is assumed to be constant during the
period of study.
4 Rotor is assumed to be of cylindrical type so that no saliency is
present.
E δ V∞ = 1∠0◦
EV sin δ
Pe =
X
Where X = Xg + XTr + XTL in p.u.
d 2 θm
J = Ta = Tm − Te N-m
dt 2
where
J = the total moment of inertia of the rotor masses in kg-m2
θm = the angular displacement of the rotor with respect to a stationary
axis in mechanical radians (rad)
t = time in seconds (s)
Tm = the mechanical or shaft torque supplied by the prime mover in N-m
Te = the net electrical or electromagnetic torque in N-m
Ta = the net accelerating torque in N-m
θm = ωsm t + δm
where
ωsm is the synchronous speed of the machine in mechanical radians per
second
δm is the angular displacement of the rotor in mechanical radians from the
synchronously rotating reference axis.
d 2 δm
Jωm = Pa = Pm − Pe
dt 2
Siva (IIT P) EE549 5 / 61
where
Pm = shaft power input in MW
Pe = electrical power output in MW
Pa = accelerating power in MW
2H d 2 δ
= Pm − Pe per unit
ωs dt 2
2H dω
= Pm − Pe per unit
ωs dt
dδ
= ω − ωs
dt
Pmax
Pm
δ0 δmax δ
2H d 2 δ0 2H d 2 ∆δ
+ = Pm − Pmax sin δ0 − Pmax cos δ0 ∆δ
ωs dt 2 ωs dt 2
At the stable equilibrium point (δ0 ),
2H d 2 δ0
= Pm − Pmax sin δ0 = 0
ωs dt 2
Therefore,
2H d 2 ∆δ
= −Pmax cos δ0 ∆δ
ωs dt 2
Let Ps = Pmax cos δ0 . Ps is called the synchronizing power or torque. In
per unit system, torque and power are equal.
2H d 2 ∆δ
+ Ps ∆δ = 0
ωs dt 2
2H d 2 ∆δ d∆δ
2
+D + Ps ∆δ = 0
ωs dt dt
If Ps is positive
there will be two complex conjugate roots with negative real part.
hence the system will be stable.
If Ps is negative,
there will be two real roots with one positive.
hence the system will be unstable.
and
∆δ
∆δ(t) = p e −ζωn t sin(ωd t + θ)
1 − ζ2
p
where ωd = ωn 1 − ζ 2 and θ = cos−1 (ζ).
The motion of rotor relative to the operating point is
∆δ
δ = δ0 + p e −ζωn t sin(ωd t + θ)
1 − ζ2
ωn ∆δ −ζωn t
ω = ω0 − p e sin ωd t
1 − ζ2
0.2 X = 0.3
Ed0 V∞ = 1∠0◦
Assume the per unit damping power coefficient D = 0.2. Consider a small
disturbance of ∆δ = 10◦ . For example, the breakers open and quickly close.
Determine the motion of rotor angle and the generator frequency.
0.3 + 0.3
X = 0.3 + 0.2 + = 0.65
2
The direct axis transient internal voltage is
Ed0 V∞ 1.17 × 1
Pmax = = = 1.8
X1 0.65
δ0 = 26.38◦
26.38◦
t (s)
Figure: Swing Curve
f (Hz)
50
t (s)
Figure: Frequency
dδ d 2 δ ωs (Pm − Pe ) dδ
2 2
=
dt dt H dt
d dδ 2 ωs (Pm − Pe ) dδ
=
dt dt H dt
On integration,
2 δm
ωs (Pm − Pe )
Z
dδ
= dδ
dt δ0 H
where
Pmax
Pm1
Pm0
δ0 δ1 δm δ
Z δ1 Z δm
(Pm1 − Pe )dδ = (Pe − Pm1 )dδ
δ0 δ1
Area(A1 ) = Area(A2 )
E δ V∞ = 1∠0◦
XTL
X1 = Xd0 + XTr +
2
Pe1 = Pmax1 sin δ
where
EV
Pmax1 =
X1
+ +
E δ − − 1∠0◦
Pmax1
Pmax3
A2
Pm
A1
Pmax2
δ0 δcr δmax δ
δ
δm
δcr
δ0
tcr t
Figure: Swing Curve
E δ V∞ = 1∠0◦
1 Before Fault :
XTL
X1 = Xd0 + XTr +
2
Pe1 = Pmax1 sin δ
where
EV
Pmax1 =
X1
Pmax1
Pmax3
A2
Pm
A1
δ0 δcr δmax δ
We can find the critical clearing time for this case as follows:
d 2δ ωs
2
= (Pm − Pe )
dt 2H
Since Pe = 0 during fault,
d 2δ ωs
2
= Pm
dt 2H
Integrating this, Z t
dδ ωs ωs
= Pm = Pm t
dt 0 2H 2H
On further integration,
ωs
δ= Pm t 2 + δ0
4H
Siva (IIT P) EE549 39 / 61
If δ = δcr ,
ωs 2
δcr = Pm tcr + δ0
4H
s
4H(δcr − δ0 )
tcr =
ωs P m
Let
dx
= f(x, t)
dt
where x is the state vector and f(x, t) is a vector of non linear functions.
1 Explicit Methods
1 Euler Method
2 Modified Euler Method
3 Runge-Kutta Methods
2 Implicit Methods
1 Trapezoidal Rule
x1
∆x ∆t
x0
t0 t1 t
dx
x1 = x0 + ∆x = x0 + ∆t
dt x=x0
This has to be repeated till the time reaches the final simulation time.
Since it considers only the firs derivative of x, it is referred to as a
first order method.
∆t has to be small to achieve accuracy.
Since it uses only the first order information, it may introduce error.
This process has to be repeated until the desired accuracy or the final
simulation time.
Siva (IIT P) EE549 45 / 61
Runge-Kutta (R-K) Methods
Euler and the modified Euler method require smaller time steps.
R-K methods approximate the Taylor series solution. However they do
not need derivatives higher than the first.
R-K methods use the effectiveness of higher derivatives by several
evaluations of the first derivative.
They are classified based on the number of evaluations.
The value of x at t = t0 + ∆t is
k1 + k2
x1 = x0 + ∆x = x0 +
2
where
k1 = f (x0 , t0 )∆t
k2 = f (x0 + k1 , t0 + ∆t)∆t
In general,
k1 + k2
xn+1 = xn +
2
where
k1 = f (xn , tn )∆t
k2 = f (xn + k1 , tn + ∆t)∆t
Explicit methods calculate x at any time step from the knowledge of the
values of x at previous time steps.
f (x, t)
f (x1 , t1 )
f (x0 , t0 )
t0 ∆t t1 t
0.2 X = 0.3
Ed0 V∞ = 1∠0◦
A solid three phase fault occurs at point F and the fault is cleared by opening
the faulted line.
1 Determine the critical clearing angle and the critical fault clearing
0.3 + 0.3
X1 = 0.3 + 0.2 + = 0.65
2
The direct axis transient internal voltage is
Ed0 V∞ 1.17 × 1
Pmax1 = = = 1.8
X1 0.65
δ0 = 26.38◦
X2 = ∞
Pmax2 = 0
The reactance after the fault is
Pmax3 = 1.4625
Pm
δmax = π − sin−1 ( ) = 2.3887
Pmax3
δcr ≈ 54◦
1 The dynamics of rotor field winding and the damper winding on the
generator are totally neglected in the classical model. However, they
can effect the stability of a system significantly.
2 In the classical model, the internal voltage behind the transient reac-
tance was assumed to be constant. This is not true since the rotor
filed current is controlled through an exciter and automatic voltage
regulator (AVR). Their dynamics have to be included.
3 In the classical model, Pm is assumed to be constant. But Pm depends
on speed governor and turbine dynamics. Their dynamics need to be
considered.
4 Dynamic loads like induction motors, synchronous motors, power elec-
tronic devices do affect the stability. In the classical model, they were
not taken into consideration.