Stability Slide 2
Stability Slide 2
-PAKISTAN CENTER
FOR ADVANCED STUDIES
IN ENERGY (USPCAS-E)
Table 1. Rotor Dynamics of Synchronous Machine (Synchronous Generator and Synchronous Motor)
Table 1. Rotor Dynamics of Synchronous Machine (Synchronous Generator and Synchronous Motor)
9
3. Rotor Dynamics and Swing Equation.
Contd.
According to the def. of rotor dynamics, “The acceleration torque is the product
of inertia and angular acceleration” , Swing equation can-be derived via
differential equation as:
Symbol Meaning
𝐽 Moment of inertia in Kg-m3
𝜃𝑚 Angular displacement of rotor
wrt. a stationary axis in
mechanical radians
(It is cal. Wrt stationary axis)
𝑇𝑚 Mechanical or shaft torque
supplied by prime mover
Note: 𝑇𝑚 and 𝑇𝑒 have positive (minus retarding torque due to
signs for synchronous generator rotational losses) in N-m
operation and negative signs for 𝑇𝑒 Net elect. or Electromagnetic
synchronous motor operation torque in N-m
𝑇𝑎 Net accelerating torque
10
3. Rotor Dynamics and Swing Equation.
Contd.
It is appropriate to measure the rotor angular position wrt reference axis
which rotates at synchronous speed.
𝜃𝑚 = 𝜔𝑠𝑚 t + 𝛿𝑚 (2)
𝝎𝒔𝒎 : It is synchronous speed of machine in mechanical rad/sec.
𝜹𝒎 : Angular displacement of rotor in mechanical radians from sync rotating ref.
𝜽𝒎 : Changes continuously under SS condition b/c its measured wrt stationary axis.
𝜽𝒎 : Can-be measure wrt a synchronous rotation axis
11
3. Rotor Dynamics and Swing Equation.
Contd.
Second derivative of eqn. (3) wrt t:
𝑑2 𝑑2
𝜃 = 𝛿𝑚 (4)
𝑑𝑡 2 𝑚 𝑑𝑡 2
𝑑 2 𝜃𝑚 𝑑2
𝐽 2 = 𝐽 2 𝛿𝑚 = 𝑇𝑎 = 𝑇𝑚 − 𝑇𝑒 (5)
𝑑𝑡 𝑑𝑡
Since, all terms are torque in N-m and we are comfortable with power terms.
Lets donate 𝜔𝑚 in rad/s in eqn. (6) and multiply with eqn. (5), we get eqn. (7):
𝑑
𝜔𝑚 = 𝜃𝑚 (6)
𝑑𝑡
𝑑2
𝐽𝜔𝑚 2 𝛿𝑚 = 𝜔𝑚 𝑇𝑎 = 𝜔𝑚 𝑇𝑚 − 𝜔𝑚 𝑇𝑒 (7)
𝑑𝑡
Kg-m2 Rad/s Rad/s N-m Rad/s N-m Electrical Power
12
Shaft power i/p to machine less rotational loss Crossing the air gap
3. Rotor Dynamics and Swing Equation.
Contd.
Eqn. (7) in watts can-be written as eqn. (8):
𝑑2
𝐽𝜔𝑚 2 𝛿𝑚 = 𝑃𝑚 − 𝑃𝑒 = 𝑃𝑎 (8)
𝑑𝑡
Practically 𝝎𝒎 ≅ 𝝎𝒔𝒎 (Sync. Speed) and has a very small difference. If difference
is large, it is when machine losses synchronism, as swing equation in eqn. (9).
𝑑2
𝑀 2 𝛿𝑚 = 𝑃𝑚 − 𝑃𝑒 (9)
𝑑𝑡
Note/Justification of inertia constant/angular momentum M: M is normally
not constant and with assumption 𝜔𝑚 ≅ 𝜔𝑠𝑚 ; thus M becomes constant.
Further, error is negligible and justify the coefficient of acceleration as constant.
M varies over wide range depending on size and type of machine.
Hydro and Turbo generators have widely different values of M.
13
4. Inertial Constant “H”
According to the def. of rotor dynamics, “The acceleration torque is the product
of inertia and angular acceleration” , Swing equation can-be written as follows:
𝐾.𝐸 𝑖𝑛 𝐾𝐽/𝑀𝐽 @ 𝑆𝑦𝑛𝑐 𝑆𝑝𝑒𝑒𝑑 Symbol Meaning
𝐻=
𝑀𝑎𝑐𝑖𝑛𝑒 𝑅𝑎𝑡𝑖𝑛𝑔 𝑖𝑛 𝐾𝑉𝐴/𝑀𝑉𝐴 𝐻 Inertia constant of rotor and turbine
𝜔𝑠 Synchronous speed of machine
2H d 2
Pm Pe 𝑃𝑚 Mechanical power applied
• Hydro and Turbo generators of same MVA have widely different values of Kinetic
energy and momentum and to express them in common way, we use H.
• Frequency of stator electrical quantities “depends on/is synchronized with” rotor
mechanical speed
• In steady state, rotor field and revolving stator field have the same speed –
synchronous speed 𝜔𝑠
• But there is an angular separation 𝛿 between them depending on mechanical power
input 𝑃𝑚 and electrical power output 𝑃𝑒 of the generator 14
4. Inertial Constant “H”
Now we will establish a relationship between M and H:
1 1
𝐾.𝐸 𝑖𝑛 𝐾𝐽/𝑀𝐽 @ 𝑆𝑦𝑛𝑐 𝑆𝑝𝑒𝑒𝑑 ∗𝐽𝑤𝑠𝑚 2 2∗𝑀𝑤𝑠𝑚
𝐻= = 2
= (10)
𝑀𝑎𝑐𝑖𝑛𝑒 𝑅𝑎𝑡𝑖𝑛𝑔 𝑖𝑛 𝐾𝑉𝐴/𝑀𝑉𝐴 𝑆(𝑀𝑎𝑐) 𝑆(𝑀𝑎𝑐)
2𝐻
Substituting M=J𝑤𝑠𝑚 = S(𝑀𝑎𝑐ℎ) in eqn. (9), we get:
𝑤𝑠𝑚
2𝐻 𝑑 2 𝑃𝑚 − 𝑃𝑒
𝑤𝑠𝑚 𝑑𝑡 2
𝛿𝑚 = S(𝑀𝑎𝑐)
= 𝑃𝑚 (p.u) − 𝑃𝑒 (p.u) (11)
𝟐𝑯 𝒅𝟐
• Where, 𝒘 is speed in mech rad/s and elect rad/s; 𝒅𝒕𝟐 𝜹𝒎 is change in angle in mech
𝒔𝒎
rad and elect rad; Pm and Pe are mech and electrical powers in per units,
respectively. 𝜹𝒎 and 𝒘𝒔𝒎 is consistent, thus By dropping subscript terms, we get
swing equation of electrical machine, that is second order differential equation.
2𝐻 𝑑 2 𝛿
𝑤𝑠 𝑑𝑡 2
= 𝑃𝑚 − 𝑃𝑒 (12) 15
4. Inertial Constant “H”
For system with electrical frequency f in Hz, eqn. 12 can be written as (13):
𝐻 𝑑2 𝛿
𝜋𝑓 𝑑𝑡 2
= 𝑃𝑚 − 𝑃𝑒 (13)
𝐻 𝑑2 𝛿
180∗𝑓 𝑑𝑡 2
= 𝑃𝑚 − 𝑃𝑒 (14)
• Where, eqn. (12-14) are basic swing equations around which stability of power system
depends.
• Electrical power (Pe) is o/p of the synchronous machine and is non-linear function of 𝜹.
Where, we have 𝜹 as function of time.
• Swing equation is second order non-linear differential equation.
• Suppose, we have several machines, we have to write down swing equation of each
machine.
• The solution of swing equation is called swing curve. 16
5. Swing Curve
• Frequency of stator electrical quantities “depends on/is synchronized
with” rotor mechanical speed
• In steady state, rotor field and revolving stator field have the same
speed – synchronous speed 𝜔𝑠
• But there is an angular separation 𝛿 between them depending on
mechanical power input 𝑃𝑚 and electrical power output 𝑃𝑒 of the
generato
• Read page 19 of Prabha Kunder for further details
Symbol Meaning
𝐻 Inertia constant of rotor and turbine
2H d 2
Pm Pe
𝜔𝑠 Synchronous speed of machine
17
5. Swing Curve. Contd.
2H d 2
Pm Pe
s dt 2
• The swing equation has pivotal role in stability analysis of power system
• Here 𝛿 is the unknown, Solution of equation gives 𝛿 values with changing time
18
5. Swing Curve. Contd.
• We have to maintain equilibrium between the torque supplied by the prime over
and the electromagnetic torque developed by the synchronous generator.
• The Typical swing curves are as follows:
19
6. Power Vs. Angle relationship.
Two Machine System
• An important characteristic that has a strong bearing on power system
stability is the relationship between interchange power and angular
positions of rotors of the synchronous machines.
• This relationship between two or more synchronous machines is highly
non-linear. (And it will make our swing equation non-linear too!)
SS stability
Limit
∆𝑇𝑒 =Ts ∆𝛿 + j Td ∆𝜔
When rotor is oscillating (oscillations are sinusoidal in
nature), we can represent this torque as a phasor in plane
where, ∆𝛿 is in x-axis and ∆𝜔 is in y-axis.
EEE 802 Advanced Power System Stability
8. Small signal stability. Contd.
A. Non-oscillatory stability: A. Non-oscillatory stability:
Traits: -ve Ts and +ve Td. Traits: +ve Ts and +ve Td.
(With constant field voltage) (With constant field voltage)
(unstable) (stable if both Ts and Td are +ve)
B. Oscillatory stability:
Traits: +ve Ts and -ve Td.
(With excitation control)
(unstable)
30
8. Small signal stability. Contd.
Various type of oscillatory stability/instability.
We came across various modes of instability in large interconnected
power system.
1. Local Mode (Frequency of oscillations 0.7 - 2Hz): Mode in which
synchronous machines in power plant oscillates wrt. rest of system.
2. Inter-area Mode (Frequency of oscillations 0.2 - 1Hz): Mode in
which a group of machine oscillates with another group of
machine, connected by a tie-line. (Also, H is considered large)
3. Control Modes: The frequency of different modes depends upon
actual control phenomenon with no particular frequency specified.
4. Torsional Modes: The torsional modes is another phenomena
which occurs in the steam turbine and the different frequencies
which we come across are for typical mass system it comes out be
16.3, 24.1, 30.3 and 44 Hz.
EEE 802 Advanced Power System Stability
9. Transient stability
Definition:
Transient stability is the ability of power system to maintain
synchronism when subjected to severe transient disturbances. i.e.
1. Transmission system faults.
2. Sudden large load changes.
3. Loss of generating units.
4. Line tripping/switching.
Reasoning:
A. Whenever large disturbances occur then the resulting response involves
large excursion of rotor angles, the delta varies over a wide range and
therefore the power output is highly non-linear with respect to delta.
B. Thus, transient stability analysis requires the solution of non-linear swing
equations.
Reasoning:
A. For each operating condition the behavior of the system is going to
be different.
B. The severity of the disturbance, depends upon the type of fault,
location of fault and duration of fault.
C. The quantum of load which has come out of the system.
D. The amount of generation that has been taken out of the system.
Reasoning:
A. The system is initially operating in a bulk steady state operate
condition, the moment disturbance occurs there will be a change in
the topology of the system.
B. If the (post disturbance) system is stable, the new steady state
operating condition will be different from the previous (last) one.