7.3 Symmetrical Components & Fortescue Sequence Transformation Matrices

7.3 Symmetrical Components & Fortescue Sequence Transformation Matrices

Executive Overview: Symmetrical components represent the most powerful mathematical tool in power systems engineering for analyzing unsymmetrical faults (Single Line-to-Ground, Line-to-Line, and Double Line-to-Ground) and unbalanced load distributions. Formulated by Charles LeGeyt Fortescue in 1918, the method decouples a set of three coupled, unbalanced phase vectors into three independent, symmetrical sets of balanced sequence phasors. The PE Power exam rigorously tests the complex $\alpha$ operator, sequence transformation matrices, zero-sequence neutral current relationships, and the physical impacts of sequence currents on rotating machinery.

---

1. Fortescue's Theorem & The Three Symmetrical Sets

Fortescue's Theorem states that any set of $N$ unbalanced, coupled phasors can be resolved into $N$ symmetrical sets of balanced phasors. For a three-phase power system ($N = 3$), any arbitrary unbalanced voltage or current phasors $(\mathbf{V}_a, \mathbf{V}_b, \mathbf{V}_c)$ can be decomposed into three balanced components:

POSITIVE SEQUENCE (1):        NEGATIVE SEQUENCE (2):        ZERO SEQUENCE (0):
(Same sequence as system)     (Reversed phase sequence)     (Identical in-phase vectors)

          Va1                           Va2                       Va0 = Vb0 = Vc0
           ^                             ^\                              ^ ^ ^
           |                             |                              | | |
           |                             |                              | | |
      /    |    \                   /    |    \                         | | |
     /     |     \                 /     |     \                        | | |
    v      |      v               v      |      v                       | | |
   Vc1     |     Vb1             Vb2     |     Vc2                      +-+-+
(Displaced by 120° lag)       (Displaced by 120° lead)        (Magnitude & Phase Equal)

Definitions of the Sequence Components

1. Positive-Sequence Components ($\mathbf{V}_{a1}, \mathbf{V}_{b1}, \mathbf{V}_{c1}$):
   • Three phasors of equal magnitude, displaced by $120^\circ$, having the same phase sequence ($abc$) as the original power system:
   $$\mathbf{V}_{a1} = V_1 \angle \theta_1, \quad \mathbf{V}_{b1} = \mathbf{V}_{a1} \angle -120^\circ = \alpha^2 \mathbf{V}_{a1}, \quad \mathbf{V}_{c1} = \mathbf{V}_{a1} \angle +120^\circ = \alpha \mathbf{V}_{a1}$$
2. Negative-Sequence Components ($\mathbf{V}_{a2}, \mathbf{V}_{b2}, \mathbf{V}_{c2}$):
   • Three phasors of equal magnitude, displaced by $120^\circ$, having the reversed phase sequence ($acb$) relative to the original power system:
   $$\mathbf{V}_{a2} = V_2 \angle \theta_2, \quad \mathbf{V}_{b2} = \mathbf{V}_{a2} \angle +120^\circ = \alpha \mathbf{V}_{a2}, \quad \mathbf{V}_{c2} = \mathbf{V}_{a2} \angle -120^\circ = \alpha^2 \mathbf{V}_{a2}$$
3. Zero-Sequence Components ($\mathbf{V}_{a0}, \mathbf{V}_{b0}, \mathbf{V}_{c0}$):
   • Three phasors of equal magnitude and identical phase angle (zero relative phase displacement):
   $$\mathbf{V}_{a0} = \mathbf{V}_{b0} = \mathbf{V}_{c0} = V_0 \angle \theta_0$$

---

2. The Complex Phase Rotation Operator ($\alpha$)

The complex operator $\alpha$ (often designated as $a$ in older literature) rotates any complex phasor by $120^\circ$ counter-clockwise without altering its magnitude:

$$\alpha = 1\angle 120^\circ = e^{j 2\pi/3} = -\frac{1}{2} + j\frac{\sqrt{3}}{2} \approx -0.5000 + j0.8660$$

+---------------------------------------------------------------------------------------------------+
| ALGEBRAIC IDENTITIES OF THE COMPLEX ALPHA OPERATOR (HIGH-YIELD MEMORIZATION)                     |
+---------------------------------------+-----------------------------------------------------------+
| Fundamental Power Operations          | Trigonometric / Rectangular Equivalence                   |
+---------------------------------------+-----------------------------------------------------------+
| $\alpha = 1\angle 120^\circ$          | $-\frac{1}{2} + j\frac{\sqrt{3}}{2} = -0.500 + j0.866$    |
| $\alpha^2 = 1\angle 240^\circ = 1\angle -120^\circ$ | $-\frac{1}{2} - j\frac{\sqrt{3}}{2} = -0.500 - j0.866$ |
| $\alpha^3 = 1\angle 360^\circ = 1\angle 0^\circ$    | $1.000 + j0.000 = 1$                                |
| $\alpha^4 = \alpha, \quad \alpha^5 = \alpha^2$      | Periodic modulo 3                                         |
+---------------------------------------+-----------------------------------------------------------+
| Summation & Difference Identities     | Rectangular & Polar Values                                |
+---------------------------------------+-----------------------------------------------------------+
| $1 + \alpha + \alpha^2 = 0$           | Fundamental zero-sum vector closure                       |
| $1 - \alpha = \sqrt{3}\angle -30^\circ$ | $\frac{3}{2} - j\frac{\sqrt{3}}{2} = 1.500 - j0.866$      |
| $1 - \alpha^2 = \sqrt{3}\angle +30^\circ$ | $\frac{3}{2} + j\frac{\sqrt{3}}{2} = 1.500 + j0.866$   |
| $\alpha - \alpha^2 = j\sqrt{3} = \sqrt{3}\angle +90^\circ$ | $0.000 + j1.732$                            |
| $\alpha^2 - \alpha = -j\sqrt{3} = \sqrt{3}\angle -90^\circ$ | $0.000 - j1.732$                           |
+---------------------------------------+-----------------------------------------------------------+

---

3. Transformation Matrices ($\mathbf{A}$ and $\mathbf{A}^{-1}$)

Synthesis Equation (Sequence $\to$ Phase Transformation)

Expressing each physical phase phasor as the superposition of its symmetrical components:

$$\mathbf{V}_a = \mathbf{V}_{a0} + \mathbf{V}_{a1} + \mathbf{V}_{a2}$$ $$\mathbf{V}_b = \mathbf{V}_{b0} + \mathbf{V}_{b1} + \mathbf{V}_{b2} = \mathbf{V}_{a0} + \alpha^2 \mathbf{V}_{a1} + \alpha \mathbf{V}_{a2}$$ $$\mathbf{V}_c = \mathbf{V}_{c0} + \mathbf{V}_{c1} + \mathbf{V}_{c2} = \mathbf{V}_{a0} + \alpha \mathbf{V}_{a1} + \alpha^2 \mathbf{V}_{a2}$$

In compact matrix notation, defining the Transformation Matrix $\mathbf{A}$:

$$\begin{bmatrix} \mathbf{V}_a \\ \mathbf{V}_b \\ \mathbf{V}_c \end{bmatrix} = \begin{bmatrix} 1 & 1 & 1 \\ 1 & \alpha^2 & \alpha \\ 1 & \alpha & \alpha^2 \end{bmatrix} \begin{bmatrix} \mathbf{V}_{a0} \\ \mathbf{V}_{a1} \\ \mathbf{V}_{a2} \end{bmatrix} \implies \mathbf{V}_{abc} = \mathbf{A} \mathbf{V}_{012}$$

Analysis Equation (Phase $\to$ Sequence Transformation)

Inverting matrix $\mathbf{A}$ yields $\mathbf{A}^{-1}$:

$$\mathbf{A}^{-1} = \frac{1}{3} \begin{bmatrix} 1 & 1 & 1 \\ 1 & \alpha & \alpha^2 \\ 1 & \alpha^2 & \alpha \end{bmatrix}$$ $$\begin{bmatrix} \mathbf{V}_{a0} \\ \mathbf{V}_{a1} \\ \mathbf{V}_{a2} \end{bmatrix} = \frac{1}{3} \begin{bmatrix} 1 & 1 & 1 \\ 1 & \alpha & \alpha^2 \\ 1 & \alpha^2 & \alpha \end{bmatrix} \begin{bmatrix} \mathbf{V}_a \\ \mathbf{V}_b \\ \mathbf{V}_c \end{bmatrix} \implies \mathbf{V}_{012} = \mathbf{A}^{-1} \mathbf{V}_{abc}$$

Individual Component Expansion Formulas

$$\mathbf{I}_{a0} = \frac{1}{3}(\mathbf{I}_a + \mathbf{I}_b + \mathbf{I}_c)$$ $$\mathbf{I}_{a1} = \frac{1}{3}(\mathbf{I}_a + \alpha \mathbf{I}_b + \alpha^2 \mathbf{I}_c)$$ $$\mathbf{I}_{a2} = \frac{1}{3}(\mathbf{I}_a + \alpha^2 \mathbf{I}_b + \alpha \mathbf{I}_c)$$

Exam Trap Alert: Pay meticulous attention to row 2 and row 3 of matrices $\mathbf{A}$ and $\mathbf{A}^{-1}$:

• In $\mathbf{A}$ (Sequence to Phase): Row 2 (Phase $b$) is $[1, \alpha^2, \alpha]$ because positive sequence lags by $120^\circ$ ($\alpha^2$).
• In $\mathbf{A}^{-1}$ (Phase to Sequence): Row 2 (Sequence 1) is $\frac{1}{3}[1, \alpha, \alpha^2]$ because computing positive sequence rotates phase $b$ forward by $+120^\circ$ ($\alpha$).

4. Physical Significance of Sequence Quantities in Power Apparatus

Symmetrical components are not merely mathematical abstractions; each sequence component corresponds to distinct physical electromagnetic phenomena in power systems equipment.

1. Zero-Sequence ($0$) & Ground Return Paths

• The zero-sequence current represents the co-phasal component flowing in all three phase conductors simultaneously.
• Kirchhoff's Current Law at a grounded Wye neutral junction demonstrates: $$\mathbf{I}_N = \mathbf{I}_a + \mathbf{I}_b + \mathbf{I}_c = 3 \mathbf{I}_{a0}$$
• Ground Path Prerequisite: Zero-sequence current CANNOT flow unless there is a complete physical neutral or earth return path back to a grounded system source.
• In 3-wire Delta systems or ungrounded Wye systems, $\mathbf{I}_{a0} = 0$ identically under all steady-state and fault conditions.

2. Negative-Sequence ($2$) & Rotor Surface Thermal Damage

• When negative-sequence current ($I_2$) flows in the stator windings of a synchronous generator or induction motor, it produces a stator magnetic field rotating at synchronous speed $\omega_s$ in the opposite direction to the mechanical rotor rotation (which spins at $+\omega_s$).
• The relative velocity between the negative-sequence stator field and the rotor is: $$\omega_{relative} = \omega_s - (-\omega_s) = 2 \omega_s \quad (120\text{ Hz in a } 60\text{ Hz system})$$
• This $120\text{ Hz}$ flux cuts the rotor iron forging, slot wedges, and retaining rings, inducing massive double-frequency eddy currents.
• Because of the extreme skin effect at $120\text{ Hz}$, these currents concentrate entirely within the outer $2\text{ to }5\text{ mm}$ of the rotor surface, causing localized melting and catastrophic mechanical failure within seconds.
• Synchronous generators are assigned strict thermal limits per ANSI C50.13: $I_2^2 t = K$ (where $K$ ranges from 5 to 30 depending on generator cooling), protected by Device 46 (Negative-Sequence Overcurrent Relay).

3. Symmetrical Component Complex Power Invariance

Total three-phase complex power can be calculated directly in the sequence domain without transforming back to phase quantities:

$$\mathbf{S}_{3\phi} = \mathbf{V}_{abc}^T \mathbf{I}_{abc}^* = 3 \left( \mathbf{V}_{a0}\mathbf{I}_{a0}^* + \mathbf{V}_{a1}\mathbf{I}_{a1}^* + \mathbf{V}_{a2}\mathbf{I}_{a2}^* \right)$$

---

5. Comprehensive Step-by-Step Worked Symmetrical Component Decomposition

Problem Statement

An unsymmetrical fault on a 3-phase, 4-wire grounded system results in the following unbalanced line currents (RMS):

• $\mathbf{I}_a = 120.0\angle 0^\circ\text{ A}$
• $\mathbf{I}_b = 100.0\angle -120^\circ\text{ A}$
• $\mathbf{I}_c = 60.0\angle +90^\circ\text{ A}$

Calculate:

1. The zero-sequence current phasor $\mathbf{I}_{a0}$ and the total neutral return current $\mathbf{I}_N$.
2. The positive-sequence current phasor $\mathbf{I}_{a1}$.
3. The negative-sequence current phasor $\mathbf{I}_{a2}$.
4. Synthesize phase $b$ positive and negative sequence components ($\mathbf{I}_{b1}, \mathbf{I}_{b2}$).

=========================================================================================
CALCULATION WORKFLOW & SOLUTION:
=========================================================================================

Step 1: Convert Phase Currents to Rectangular Form
  I_a = 120.00 /_ 0°   = 120.00 + j0.00 A
  I_b = 100.00 /_ -120° = 100.00 * (-0.5000 - j0.866025) = -50.00 - j86.60 A
  I_c = 60.00 /_ +90°  = 0.00 + j60.00 A

Step 2: Compute Zero-Sequence Current (I_a0) and Neutral Current (I_N)
  Sum the phase currents:
    I_sum = I_a + I_b + I_c
          = (120.00 - 50.00 + 0.00) + j(0.00 - 86.6025 + 60.00)
          = 70.00 - j26.6025 A

  Divide by 3 for zero-sequence current:
    I_a0 = I_sum / 3 = (70.00 - j26.6025) / 3
         = 23.333 - j8.868 A

  Convert I_a0 to polar form:
    |I_a0| = sqrt(23.333^2 + (-8.868)^2) = sqrt(544.43 + 78.64) = sqrt(623.07) = 24.96 A
    theta_0 = arctan(-8.868 / 23.333) = -20.81°
    -> I_a0 = 24.96 /_ -20.81° A

  Compute Neutral Return Current I_N:
    I_N = 3 * I_a0 = I_sum = 70.00 - j26.60 A
    |I_N| = 3 * 24.96 A = 74.88 A
    -> I_N = 74.88 /_ -20.81° A

Step 3: Compute Positive-Sequence Current (I_a1)
  Formula: I_a1 = (1/3) * [ I_a + alpha * I_b + alpha^2 * I_c ]

  Evaluate each term:
    Term 1: I_a = 120.00 /_ 0° = 120.00 + j0.00 A
    Term 2: alpha * I_b = (1 /_ 120°) * (100.00 /_ -120°) = 100.00 /_ 0° = 100.00 + j0.00 A
    Term 3: alpha^2 * I_c = (1 /_ 240°) * (60.00 /_ 90°) = 60.00 /_ 330° = 60.00 /_ -30° A
            = 60.00 * (cos(-30°) + j sin(-30°))
            = 60.00 * (0.866025 - j0.5000)
            = 51.962 - j30.00 A

  Sum terms inside bracket:
    I_a1_sum = (120.00 + 100.00 + 51.962) + j(0.00 + 0.00 - 30.00)
             = 271.962 - j30.00 A

  Divide by 3:
    I_a1 = (271.962 - j30.00) / 3
         = 90.654 - j10.00 A

  Convert I_a1 to polar form:
    |I_a1| = sqrt(90.654^2 + (-10.00)^2) = sqrt(8218.15 + 100.0) = sqrt(8318.15) = 91.20 A
    theta_1 = arctan(-10.00 / 90.654) = -6.29°
    -> I_a1 = 91.20 /_ -6.29° A

Step 4: Compute Negative-Sequence Current (I_a2)
  Formula: I_a2 = (1/3) * [ I_a + alpha^2 * I_b + alpha * I_c ]

  Evaluate each term:
    Term 1: I_a = 120.00 /_ 0° = 120.00 + j0.00 A
    Term 2: alpha^2 * I_b = (1 /_ 240°) * (100.00 /_ -120°) = 100.00 /_ 120° A
            = 100.00 * (-0.5000 + j0.866025) = -50.00 + j86.603 A
    Term 3: alpha * I_c = (1 /_ 120°) * (60.00 /_ 90°) = 60.00 /_ 210° = 60.00 /_ -150° A
            = 60.00 * (-0.866025 - j0.5000) = -51.962 - j30.00 A

  Sum terms inside bracket:
    I_a2_sum = (120.00 - 50.00 - 51.962) + j(0.00 + 86.603 - 30.00)
             = 18.038 + j56.603 A

  Divide by 3:
    I_a2 = (18.038 + j56.603) / 3
         = 6.013 + j18.868 A

  Convert I_a2 to polar form:
    |I_a2| = sqrt(6.013^2 + 18.868^2) = sqrt(36.156 + 355.989) = sqrt(392.145) = 19.80 A
    theta_2 = arctan(18.868 / 6.013) = +72.33°
    -> I_a2 = 19.80 /_ +72.33° A

Step 5: Verify Synthesis of Phase a Current
  I_a = I_a0 + I_a1 + I_a2
      = (23.333 - j8.868) + (90.654 - j10.000) + (6.013 + j18.868)
      = (23.333 + 90.654 + 6.013) + j(-8.868 - 10.000 + 18.868)
      = 120.00 + j0.00 A = 120.00 /_ 0° A  (EXACT MATCH - FULLY VERIFIED)

Step 6: Phase b Sequence Components
  I_b1 = alpha^2 * I_a1 = (1 /_ 240°) * (91.20 /_ -6.29°) = 91.20 /_ 233.71° = 91.20 /_ -126.29° A
  I_b2 = alpha * I_a2   = (1 /_ 120°) * (19.80 /_ 72.33°) = 19.80 /_ 192.33° = 19.80 /_ -167.67° A
=========================================================================================

---

6. Common Exam Traps & Strategic Pitfalls

• The $\frac{1}{3}$ Factor Omission in Sequence Decompositions: Forgetting to divide by 3 when computing $\mathbf{I}_{a0}, \mathbf{I}_{a1}, \mathbf{I}_{a2}$ from phase currents. The synthesis matrix $\mathbf{A}$ has no $1/3$, but the analysis matrix $\mathbf{A}^{-1}$ always carries the $1/3$ factor.
• Neutral Current vs Zero-Sequence Current Equivalence: Equating neutral current to zero-sequence current ($I_N = I_{a0}$). The neutral current is the sum of all three phases, making it three times the zero sequence current: $\mathbf{I}_N = 3\mathbf{I}_{a0}$.
• Assuming Zero Sequence Exists in 3-Wire Lines: Attempting to calculate a nonzero zero-sequence current for a 3-phase, 3-wire system without ground return. By definition, in any 3-wire ungrounded system, $I_a + I_b + I_c = 0$, so $\mathbf{I}_{a0} = 0$ identically.
• Negative Sequence Frequency Error: Believing negative sequence currents in the stator produce a negative frequency on the rotor. The physical frequency of rotor currents induced by negative sequence stator currents is $2f = 120\text{ Hz}$ ($2\times$ system frequency), not $-60\text{ Hz}$ or $0\text{ Hz}$.