Subsynchronous Resonance (SSR) Instability in Series-Compensated Generator Systems

Electromechanical Fundamentals and the Origin of Subsynchronous Resonance

Subsynchronous resonance (SSR) instability represents one of the most complex and destructive phenomena in very high-voltage alternating current (HVAC) electrical power systems. It typically occurs in synchronous generation plants connected to ultra-long-distance transmission grids that employ series capacitive compensation. The fundamental purpose of series compensation is to reduce the total inductive reactance of the transmission line, thereby enhancing active power transfer capability, optimizing the phase angle, and improving the transient stability of the interconnected system. However, the introduction of capacitor banks in series with the grid inductances creates series LC resonant circuits whose natural frequencies lie below the fundamental system frequency (50 Hz or 60 Hz).

From the perspective of distributed and lumped circuit theory, the natural electrical oscillation frequency of a series-compensated transmission system is defined by the following fundamental expression:

fer = f_0 \sqrt{ \frac{X_C}{X_L} } = f_0 \sqrt{ kse }

Where $f_0$ represents the fundamental system frequency (50 Hz or 60 Hz), $X_C$ is the total capacitive reactance of the series banks, $X_L$ is the total equivalent inductive reactance of the system viewed from the generator terminals, and $kse$ is the series compensation degree expressed in per-unit or percentage. Given that compensation degrees in long lines commonly range between 25% and 70%, the resonant electrical frequency $fer$ turns out to be strictly lower than the fundamental frequency. Consequently, any transient or harmonic disturbance in the grid excites this natural electrical frequency, giving rise to subsynchronous currents and voltages.

When these subsynchronous currents circulate through the stator windings of the synchronous generator, they generate a rotating magnetic field in the air gap of the machine that rotates at an angular velocity lower than the synchronous speed of the rotor. This subsynchronous rotating field interacts with the direct-current magnetic field of the rotor, producing oscillatory electromagnetic torques on the generator shaft whose pulsation frequency $f_m$ is given by the difference between the fundamental frequency and the resonant electrical frequency:

f_m = f_0 - fer = f_0 \left( 1 - \sqrt{kse} \right)

If the frequency of this oscillatory torque coincides with or closely approximates one of the natural torsional vibration frequencies of the generating unit's shaft train (which comprises the high-pressure, intermediate-pressure, and low-pressure turbines, as well as the generator rotor and the exciter), severe mechanical-electrical coupling occurs. This mechanism is formally known as Torsional Interaction (TI). The machine acts as a mechanical energy amplifier through negative or positive feedback of the torsional damping, potentially leading to extreme torsional fatigue and catastrophic shaft fracture within seconds.

Classification and Physical Mechanisms of SSR

The global phenomenon of subsynchronous resonance is not a monolithic event; rather, it is formally subdivided into three distinct yet interrelated physical mechanisms according to the IEEE SSR Task Force:

Rotor Induction Generator Effect (RIGE)

This phenomenon is strictly electrical in nature and can occur even with an ideally rigid rotor (devoid of torsional degrees of freedom). When subsynchronous currents of frequency $fer$ flow through the stator, they establish a magnetic field that rotates at a speed lower than that of the rotor. Viewed from the rotor reference frame spinning at synchronous speed \omega_0, this field appears at a subsynchronous slip frequency $fslip = f_0 - fer$ . This relative frequency induces currents in the damper circuits and in the massive solid-steel rotor body.

The equivalent impedance viewed from the stator at the subsynchronous frequency includes the reflected effective rotor resistance. Due to the negative slip for frequencies lower than the fundamental ( $fer < f_0$ ), the apparent rotor resistance viewed by the stator at that subsynchronous frequency can become negative. A total negative stator and rotor resistance at frequency $fer$ means that the circuit absorbs energy from the grid and amplifies it, or that the circuit negatively damps the electrical oscillations, provoking electrical self-excitation if the net circuit resistance Rtotal(fer) < 0.

Torsional Interaction (TI)

Torsional Interaction occurs when the electrical system and the mechanical shaft train couple via the air-gap magnetic field. Suppose the shaft train possesses a natural torsional frequency $fmech$ . If a torsional oscillation occurs in the shaft at said frequency, the angular variations of the rotor modulate the stator currents, generating sidebands in the electrical spectrum at frequencies f_0 \pm fmech. If one of these sidebands coincides with the complementary natural frequency of the series electrical circuit ( $fer = f_0 - fmech$ ), a closed feedback loop is established.

Mathematically, this coupling is evaluated by calculating the torsional electrical damping coefficient $D_e$ for each torsional mode $i$ . The logarithmic decrement or total damping $Dtotal$ of a torsional mode is the sum of the inherent mechanical damping $D_m$ and the induced electrical damping $D_e$ :

D_{total, i} = D_{m, i} + D_{e, i}

If the grid resistance and the characteristics of the LC circuits cause $D_{e, i}$ to be negative and its magnitude exceeds the positive mechanical damping ( $|D_{e, i}| > D_{m, i}$ ), the total system damping becomes negative ( $D_{total, i} < 0$ ). This leads to an exponential instability where torsional oscillations grow uncontrollably until the mechanical destruction of the shaft couplings.

Torque Amplification (TA)

Unlike RIGE and TI, which are dynamic and small-signal instability phenomena, Torque Amplification is a large-magnitude transient phenomenon. It occurs when severe short circuits take place in the vicinity of series-capacitor-compensated lines, followed by the opening and closing operations of power circuit breakers. The sudden discharge and recharge of the series capacitor banks interact with transient harmonics and fault currents, injecting a severe electromagnetic torque impulse into the generator shaft.

If the frequency of this impulse coincides with the natural modes of the shaft, the shaft train experiences mechanical stresses that can vastly exceed the nominal design torque, accumulating severe plastic fatigue damage in a single transient event.

Mathematical Modeling and Small-Signal Analysis in Multimode Systems

To perform a rigorous analysis of subsynchronous resonance stability in advanced industrial and transmission environments, linearized state-space modeling around a steady-state operating point is employed. The complete system is decomposed into three coupled subsystems: the series-compensated transmission electrical network, the synchronous generator with its Park's equations (d-q axes), and the multi-mass mechanical model of the shaft train.

The mechanical model of the shaft train with $N$ rotating masses is represented by a system of second-order ordinary differential equations:

J_i \frac{d^2 \theta_i}{dt^2} + Dmi \frac{d \theta_i}{dt} + \sum_{j=1}^{N} Kij (\theta_i - \theta_j) = Tmi - Tei, \quad i = 1, 2, \dots, N

Where $J_i$ is the moment of inertia of mass $i$ , \theta_i is the rotor angular position, $Dmi$ is the self-damping coefficient, $Kij$ is the torsional shaft stiffness between masses $i$ and $j$ , $Tmi$ is the mechanical torque applied by the turbine, and $Tei$ is the electromagnetic torque induced on the generator mass.

The electrical system comprises Park's transformation equations for the stator and rotor, coupled with the differential equations of the series R-L-C transmission circuits:

\frac{d i_d}{dt} = \frac{1}{L_d} \left( -R_a i_d + \omega \psi_q - v_d \right)

\frac{d i_q}{dt} = \frac{1}{L_q} \left( -R_a i_q - \omega \psi_d - v_q \right)

The compensated transmission line equations are expressed in network-synchronized axis coordinates via differential equations incorporating load dynamics and series capacitors:

\frac{d vcd}{dt} = \omega_0 i_d + \frac{1}{Cse} i_d

\frac{d vcq}{dt} = -\omega_0 i_q + \frac{1}{Cse} i_q

By linearizing this set of nonlinear equations around a stable operating point ( $X_0, U_0$ ), the closed-loop system state matrix is obtained:

\Delta \dot{X} = A \Delta X + B \Delta U

\Delta Y = C \Delta X + D \Delta U

The small-signal stability analysis is performed by calculating the eigenvalues \lambda_i = \sigma_i \pm j \omega_i of the state matrix $A$ . The real part \sigma_i of each eigenvalue determines the stability of the corresponding mode:

If \sigma_i < 0, the mode is asymptotically stable (positive damping).
If \sigma_i > 0, the mode is unstable (negative damping, exponential growth of oscillation).
If \sigma_i = 0, the mode lies on the marginal stability boundary.

For torsional electromechanical modes, the real part of the eigenvalue is directly related to the total damping coefficient: \sigma_i = -\frac{D_{total, i}}{2 J_i}. Therefore, an eigenvalue with a positive real part directly indicates the occurrence of Torsional Interaction (TI) instability.

Comparative Matrix of International Standards and Evaluation Criteria

SSR evaluation and mitigation are governed by strict international standards issued by the IEEE and the IEC. The following comparative matrix details critical parameters, normative limits, and operational consequences associated with SSR in high-voltage networks.

Standard / Code	Evaluation Parameter	Critical Limit / Acceptance Criterion	Associated Fault Condition	Operational and Dielectric Consequence
IEEE Std 1133	Torsional Electrical Damping ( $D_e$ )	Dtotal \ge 0 under any dispatch condition ( $N-0, N-1$ )	Torsional Interaction (TI) Instability	Accumulated torsional fatigue, deformation, and catastrophic shaft fracture.
IEEE Std 693	Coupled Seismic and Torsional Response	Maximum mechanical stresses below 80% of yield strength	Combined torque transient amplification	Structural damage to turbogenerator supports and substation foundations.
IEC 60383 / IEC 60871	Overvoltages in Series Capacitor Banks	Vpeak \le 1.8 p.u. under steady-state subsynchronous regime	Prolonged subsynchronous resonance and ferroresonance	Dielectric puncture of capacitor elements, triggering of metal-oxide varistors (MOV).
IEEE Std 1547 / C50.13	Subsynchronous Currents in Stator Windings	Isub \le 0.05 p.u. continuous; Isub \le 0.15 p.u. for 10s	Rotor Induction Generator Effect (RIGE) and localized heating	Critical overheating in rotor wedges, retaining rings, and stator core.

Forensic Failure Analysis of Critical Equipment Under SSR

The presence of subsynchronous currents and torques generates severe stresses that transcend the synchronous generator, profoundly affecting all associated substation infrastructure, transformation, and switching systems. Forensic failure analysis reveals degradation mechanisms in the following components:

Power Transformers and Autotransformers

Subsynchronous current components ( $fer$ ) superimposed on the fundamental frequency generate temporal asymmetries in the magnetic flux of the transformer core. Given that magnetic flux is the time integral of voltage, the presence of low-frequency magnetomotive components elevates the peak flux density \hat{B}, driving the core deep into the nonlinear magnetic saturation region.

Periodic core saturation causes:

Massive generation of higher-order harmonics (especially 2nd, 3rd, and 5th harmonics) that distort the voltage waveform.
Drastic increase in excitation currents and core losses due to hysteresis and eddy currents (Foucault currents).
Severe mechanical vibrations in the transformer tanks and windings due to nonlinear magnetostriction and Lorentz forces weighted by asymmetric currents, resulting in the loosening of mechanical clamping and premature failure of solid insulation (oil-impregnated paper).

Circuit Breakers and Switching Systems

Circuit breakers located on lines compensated with series capacitors experience extremely severe voltage recovery conditions during SSR events. Transient recovery voltages (TRV) are modulated by the natural frequency of the series circuit $fer$ , presenting very high rates of rise of transient voltage ( $du/dt$ ) and peak voltage values that can exceed the breaker's nominal interrupting capacity.

Additionally, when a breaker operates to clear a fault in the presence of subsynchronous currents, the probability of re-strike or re-conduction of the electric arc increases considerably because the current passes through zero at abnormal rates or contains superimposed low-frequency bidirectional components. This can cause the catastrophic destruction of the SF6 arc-extinction chambers and explosive failure of the breaker pole.

High-Voltage Cables and Measurement Systems

High-voltage insulated cables interconnected with the substation suffer an increase in dielectric losses and induced current losses in metallic shields (copper or lead shielding). Subsynchronous currents induce low-frequency alternating magnetic fields that traverse the shields, generating high circulating currents that overheat the polymeric insulation (XLPE). This accelerates the formation of water and electrical treeing and drastically reduces the estimated service life of the cable.

On the other hand, current transformers (CTs) and potential transformers (PTs) suffer significant ratio and phase angle errors due to magnetic saturation induced by subsynchronous current components, which disables or severely degrades the performance of numerical protection relays based on fundamental-frequency phasors.

Advanced Mitigation Strategies and Engineering Design

To guarantee the operational viability of transmission systems with massive series compensation, modern engineering implements a set of advanced technological countermeasures, divided into grid-level mitigation, generation-equipment-level mitigation, and high-speed power electronics systems.

Subsynchronous Resonance Blocking Filters (SSBF)

SSBF filters are tuned circuits installed in series with the generator stator windings or directly on the transmission line. They consist of an arrangement of reactors and capacitors connected in parallel, tuned precisely to present an extremely high impedance (a virtual open circuit) at the critical resonant electrical frequency $fer$ .

The impedance of the SSBF filter at the tuning frequency is expressed as:

ZSSBF(j\omegaer) = \frac{R + j\omegaer L}{1 - \omegaer^2 L C + j\omegaer R C} \approx \infty