13.1 Symmetrical Three-Phase Bolted Fault Calculations

13.1 Symmetrical Three-Phase Bolted Fault Calculations

Executive Overview: Symmetrical three-phase bolted faults represent the standard benchmark for short-circuit duty in electric power systems. Although three-phase faults account for less than $5%$ of all power system short circuits, they typically produce the highest mechanical and thermal stresses on substation buswork, switchgear, and protective equipment. On the NCEES PE Power examination, mastery of three-phase fault analysis requires calculating positive-sequence Thevenin impedances ($Z_{1,th}$), converting between Short-Circuit MVA ($SCMVA$) and physical fault currents, accounting for time-varying generator reactances ($X_d'', X_d', X_d$), evaluating $X/R$ ratios and DC offset decay, and integrating induction and synchronous motor backfeed contributions into circuit breaker interrupting and momentary ratings.

---

1. Balanced Short-Circuit Physics & Positive-Sequence Network

A three-phase bolted short circuit occurs when all three phase conductors ($a, b, c$) are solidly connected together with zero fault impedance ($Z_f = 0$). Because the network remains completely balanced and symmetrical during the fault:

1. Phase currents remain equal in magnitude and displaced by exactly $120^\circ$: $\mathbf{I}_a + \mathbf{I}_b + \mathbf{I}_c = 0$
2. Negative-sequence and zero-sequence currents are identically zero: $\mathbf{I}_{a2} = 0, \qquad \mathbf{I}_{a0} = 0$
3. The positive-sequence current equals the total fault current in phase $a$: $\mathbf{I}_{a1} = \mathbf{I}_a = \mathbf{I}_{f,3\phi}$

Consequently, three-phase fault calculations require solving only the single-phase positive-sequence Thevenin equivalent circuit viewed from the fault point.

Positive-Sequence Thevenin Equivalent Circuit for 3-Phase Bolted Fault:

         +---[ Z_1,th = R_1 + j*X_1 ]---+----o (+) Fault Bus
         |                              |
        (~) V_F = 1.0 /_ 0° pu          |  I_f,3ph = V_F / Z_1,th
         |                              v  (Bolted Fault: Z_f = 0)
         +------------------------------+----o (-) Neutral / Ref

The Fundamental Fault Current Equation

Using the prefault line-to-neutral voltage at the fault point (typically $\mathbf{V}_F = 1.0\angle 0^\circ\text{ pu}$ under nominal no-load or operating conditions):

$\mathbf{I}_{f,3\phi,pu} = \frac{\mathbf{V}_F}{\mathbf{Z}_{1,th}} = \frac{1.0\angle 0^\circ}{\mathbf{Z}_{1,th}} \quad [\text{pu}]$

To convert the per-unit fault current into physical amperes, multiply by the system base current ($I_{base}$) at the fault voltage level:

$I_{base} = \frac{S_{base,3\phi}}{\sqrt{3} \cdot V_{base,LL}} \quad \Longleftrightarrow \quad I_{f,3\phi} = |\mathbf{I}_{f,3\phi,pu}| \cdot I_{base} \quad [\text{A}]$

Short-Circuit Apparent Power ($SCMVA$)

Utilities and industrial facilities frequently specify fault availability in terms of Short-Circuit MVA ($SCMVA$) or Short-Circuit Capacity ($SCC$):

$SCMVA = \frac{S_{base,3\phi}}{|\mathbf{Z}_{1,th,pu}|} = \sqrt{3} \cdot V_{LL,kV} \cdot I_{sc,kA} \quad [\text{MVA}]$

Rearranging gives the per-unit Thevenin impedance directly from utility-provided $SCMVA$:

$|\mathbf{Z}_{1,th,pu}| = \frac{S_{base,3\phi}}{SCMVA}$

$\begin{array}{|l|l|l|} \hline \textbf{Parameter} & \textbf{Per-Unit Formulation} & \textbf{Physical Units Formulation} \\ \hline \text{Fault Current } (I_{sc}) & I_{pu} = \frac{V_F}{Z_{1,th,pu}} & I_{sc,kA} = \frac{SCMVA}{\sqrt{3} \cdot V_{LL,kV}} \\ \hline \text{Thevenin Impedance } (Z_{th}) & Z_{th,pu} = \frac{S_{base}}{SCMVA} & Z_{th,\Omega} = \frac{V_{LL,kV}^2}{SCMVA} \\ \hline \text{Short-Circuit Power } (SCMVA) & SCMVA = \frac{S_{base}}{|Z_{1,th,pu}|} & SCMVA = \sqrt{3} \cdot V_{LL,kV} \cdot I_{sc,kA} \\ \hline \end{array}$

---

2. Time-Varying Machine Reactances ($X_d'', X_d', X_d$)

When a short circuit occurs near synchronous generators or large motors, the sudden increase in armature current induces opposing currents in the rotor damper windings (amortisseur circuits) and field windings to preserve flux linkages (Constant Flux Linkage Theorem). As these induced rotor currents decay due to internal rotor resistance, the effective machine reactance increases over time, causing the AC symmetrical fault current to decay from a high initial value to a lower steady-state value.

Time Decay of Symmetrical Fault Current in a Synchronous Machine:

Armature | 
Current  |   Subtransient Region (0 to 2 cycles): X_d''
 |I_ac|   |   |---| 
          |   \   \   Transient Region (2 to 30 cycles): X_d'
          |    \   ` - - - - - - - 
          |     \                  ` - - - - - - - - - -> Steady-State: X_d
          +--------------------------------------------------> Time (cycles/s)
              0   2               30

Formulations for Machine Symmetrical Currents

$I_{sym}'' = \frac{E_g''}{X_d''} \quad [\text{Subtransient}], \qquad I_{sym}' = \frac{E_g'}{X_d'} \quad [\text{Transient}], \qquad I_{ss} = \frac{E_g}{X_d} \quad [\text{Steady-State}]$

where $E_g'', E_g', E_g$ represent the machine internal voltages behind the respective reactances prior to the fault.

---

3. DC Offset Transient & Asymmetrical Currents

Because electric power circuits contain both inductance ($L$) and resistance ($R$), current cannot change instantaneously ($\Delta i_L = 0$). When a fault occurs at voltage angle $\alpha$, a unidirectional DC transient offset current ($i_{dc}$) is injected to satisfy boundary conditions:

$i_{total}(t) = i_{ac}(t) + i_{dc}(t) = \sqrt{2} I_{sym} \sin(\omega t + \alpha - \theta) + i_{dc}(t)$

where $\theta = \arctan(X/R)$ is the system impedance angle. To maintain continuity at $t = 0$ (assuming zero prefault current):

$i_{total}(0) = 0 \implies i_{dc}(0) = -\sqrt{2} I_{sym} \sin(\alpha - \theta) = \sqrt{2} I_{sym} \sin(\theta - \alpha)$

DC Offset Decay and Time Constant

The DC component decays exponentially according to the system $X/R$ ratio:

$i_{dc}(t) = \sqrt{2} I_{sym} \sin(\theta - \alpha) \cdot e^{-t/\tau} \quad \text{where } \tau = \frac{L}{R} = \frac{X}{\omega R} = \frac{X}{2\pi f R} \quad [\text{seconds}]$

[!IMPORTANT] Maximum DC Offset Condition: Maximum DC offset occurs when the fault initiates at a voltage zero-crossing ($\alpha = 0^\circ$) in a purely inductive circuit ($\theta = 90^\circ$), yielding $\sin(90^\circ - 0^\circ) = 1.0$ and $i_{dc}(0) = \sqrt{2} I_{sym}$. Conversely, zero DC offset occurs when the fault initiates at voltage peak ($\alpha = 90^\circ$).

Waveform Components of Asymmetrical Fault Current:

 Current | 
   i(t)  |          Envelope of Total Asymmetrical Current
         |         /------------------------------------
         |        /     ^  Peak Current: I_peak ~ 2.6 * I_sym
         |       /     / \
         |      /     /   \   DC Offset Decay: i_dc(t)
  i_dc   | - - o - - / - - \ - - - - - - - - - - - - - - - -> 0
         |          /       \     /
         |         /         \   /
         |        /           v /
         |       /-------------------------------------
         +----------------------------------------------------> Time (t)
                 0    8.33 ms (1/2 cycle)

Total RMS Asymmetrical Current

The total effective (RMS) asymmetrical current combines AC symmetrical RMS and instantaneous DC offset:

$I_{asym,RMS}(t) = \sqrt{I_{sym,RMS}^2 + \left[ i_{dc}(t) \right]^2} = I_{sym} \sqrt{1 + 2 e^{-2t/\tau}}$

Peak Instantaneous Current ($I_{peak}$)

The maximum mechanical force on busbars is proportional to the square of the instantaneous peak current ($F \propto i^2$). The peak current occurs during the first half-cycle ($t = 1/2f = 1/120\text{ s} \approx 8.333\text{ ms}$):

$I_{peak} = \sqrt{2} I_{sym} \left( 1 + e^{-\frac{1}{2 f \tau}} \right) = \sqrt{2} I_{sym} \left( 1 + e^{-\frac{\pi}{X/R}} \right) = k_{peak} \cdot I_{sym}$

For typical high-voltage substations with $X/R \ge 30$, $k_{peak} \approx 2.6 - 2.7$.

---

4. Motor Backfeed Contributions

When a short circuit depresses bus voltage, connected electric motors temporarily act as generators driven by mechanical load inertia, feeding energy back into the fault.

Induction Motors

• Operating principle: The rotating magnetic field is maintained briefly by trapped rotor flux linkages.
• Reactance model: Modeled as a voltage source behind subtransient reactance $X_m'' \approx 15% - 25%$ (roughly equal to locked-rotor impedance $Z_{LR} = 1/LRA$).
• Decay rate: Because induction machines lack a DC excitation winding, rotor flux collapses rapidly. The backfeed contribution decays to zero within $1$ to $4$ cycles ($16 - 66\text{ ms}$).
• Sizing standard: Per IEEE C37.010 / C37.13, low-voltage induction motors ($< 50\text{ hp}$) are often lumped at $4.0 - 5.0\times$ rated full-load current for momentary duties and neglected for breaker interrupting duties $> 4\text{ cycles}$.

Synchronous Motors

• Operating principle: Synchronous motors have continuous DC field excitation supplied by an external exciter or brushless system.
• Reactance model: Modeled using three distinct reactances ($X_d'', X_d', X_d$), identical to synchronous generators.
• Decay rate: Sustains fault current for several seconds through subtransient and transient periods as long as field excitation is maintained.

---

5. Circuit Breaker Application: Momentary vs. Interrupting Ratings

Circuit breakers must satisfy two distinct operational rating criteria under IEEE C37 standards:

1. Closing and Latching (Momentary / First-Cycle) Rating:
   • Tests the breaker's mechanical ability to withstand magnetic repulsion forces during the first half-cycle.
   • Sourced from subtransient reactances ($X_d''$) of all generators, synchronous motors, and induction motors.
   • Rated in peak amperes ($I_{peak}$) or first-cycle RMS asymmetrical amperes ($I_{asym,1\text{-}cycle} = 1.6 \times I_{sym}''$).
2. Interrupting (Contact Parting) Rating:
   • Tests the breaker's ability to extinguish the electrical arc across contacts opening at $2, 3, 5,\text{ or } 8\text{ cycles}$.
   • Sourced from generator $X_d'$, synchronous motor $X_d'$, and adjusted induction motor contributions.
   • Symmetrical current is multiplied by an asymmetrical multiplying factor ($MF$) derived from the system $X/R$ ratio and breaker contact parting time.

$\begin{array}{|c|c|c|c|} \hline \textbf{Breaker Speed} & \textbf{Contact Parting Time} & \textbf{Standard } X/R \text{ Baseline} & \textbf{Multiplying Factor } (MF) \text{ for } X/R = 25 \\ \hline 8\text{ Cycle} & 4.0\text{ cycles } (66.7\text{ ms}) & X/R = 15 & 1.02 \\ \hline 5\text{ Cycle} & 3.0\text{ cycles } (50.0\text{ ms}) & X/R = 15 & 1.08 \\ \hline 3\text{ Cycle} & 2.0\text{ cycles } (33.3\text{ ms}) & X/R = 15 & 1.18 \\ \hline 2\text{ Cycle} & 1.5\text{ cycles } (25.0\text{ ms}) & X/R = 15 & 1.28 \\ \hline \end{array}$

---

6. Comprehensive Worked Calculation: Industrial Facility Fault Study

Problem Statement

An industrial distribution facility operates at $4.16\text{ kV}$ (line-to-line). The system is fed from a $13.8\text{ kV}$ utility grid through a step-down transformer and includes local generation and motor loads connected to the $4.16\text{ kV}$ switchgear bus:

• System Base: $S_{base} = 10\text{ MVA}$, $V_{base,LV} = 4.16\text{ kV}$, $V_{base,HV} = 13.8\text{ kV}$.
• Utility Source: $S_{sc,util} = 500\text{ MVA}$ at $13.8\text{ kV}$ with $X/R = 15$.
• Transformer $T_1$: $10\text{ MVA}$, $13.8\text{ kV} : 4.16\text{ kV}$, impedance $Z_T = 5.5% = 0.055\text{ pu}$ with $X/R = 10$.
• Local Synchronous Generator $G_1$: $5.0\text{ MVA}$, $4.16\text{ kV}$, subtransient reactance $X_d'' = 15% = 0.15\text{ pu}$ (on $5\text{ MVA}$ base), $X/R = 25$.
• Induction Motor Load $M_1$: Total rating $3.0\text{ MVA}$, $4.16\text{ kV}$, subtransient reactance $X_m'' = 20% = 0.20\text{ pu}$ (on $3\text{ MVA}$ base), $X/R = 8$.
• Prefault Bus Voltage: $\mathbf{V}_F = 1.0\angle 0^\circ\text{ pu}$.

Calculate:

1. Per-unit branch impedances on the $10\text{ MVA}$ system base.
2. The combined positive-sequence Thevenin impedance $\mathbf{Z}_{1,th}$ and effective $X/R$ ratio at the $4.16\text{ kV}$ bus.
3. Initial symmetrical subtransient fault current ($I_{sym}''$) in per-unit and physical kiloamperes ($kA$).
4. Total Short-Circuit MVA ($SCMVA$).
5. Maximum instantaneous peak asymmetrical current ($I_{peak}$) at $t = 8.33\text{ ms}$.
6. Breaker asymmetrical interrupting duty at $3\text{ cycles}$ ($50\text{ ms}$).

============================== STEP-BY-STEP SOLUTION ==============================

Step 1: Convert All Branch Impedances to Common 10 MVA Base
  Base Current at 4.16 kV:
  I_base,LV = S_base / (sqrt(3) * V_base,LV) = 10,000 kVA / (sqrt(3) * 4.16 kV) 
            = 1,387.86 A = 1.38786 kA

  a) Utility Grid Impedance:
     |Z_util,pu| = S_base / S_sc,util = 10 MVA / 500 MVA = 0.0200 pu
     With X/R = 15: theta_util = arctan(15) = 86.186 deg
     R_util = 0.0200 / sqrt(1 + 15^2) = 0.0200 / 15.0333 = 0.00133 pu
     X_util = 15 * 0.00133 = 0.01996 pu
     Z_util = 0.00133 + j0.01996 pu

  b) Step-Down Transformer T_1:
     |Z_T,pu| = 0.0550 pu (Already on 10 MVA base)
     With X/R = 10: theta_T = arctan(10) = 84.289 deg
     R_T = 0.0550 / sqrt(1 + 10^2) = 0.0550 / 10.0499 = 0.00547 pu
     X_T = 10 * 0.00547 = 0.05473 pu
     Z_T = 0.00547 + j0.05473 pu

     Combined Utility + Transformer Branch:
     Z_grid = Z_util + Z_T = (0.00133 + 0.00547) + j(0.01996 + 0.05473)
            = 0.00680 + j0.07469 pu = 0.07500 /_ 84.80 deg pu
     Y_grid = 1 / Z_grid = 1.209 - j13.279 pu

  c) Local Synchronous Generator G_1 (Base Change from 5 MVA to 10 MVA):
     X_G1,new = X_G1,old * (S_base,new / S_base,old) = 0.15 * (10 / 5) = 0.3000 pu
     R_G1 = X_G1 / (X/R) = 0.3000 / 25 = 0.0120 pu
     Z_G1 = 0.0120 + j0.3000 pu = 0.30024 /_ 87.71 deg pu
     Y_G1 = 1 / Z_G1 = 0.133 - j3.328 pu

  d) Induction Motor Load M_1 (Base Change from 3 MVA to 10 MVA):
     X_M1,new = X_M1,old * (S_base,new / S_base,old) = 0.20 * (10 / 3) = 0.6667 pu
     R_M1 = X_M1 / (X/R) = 0.6667 / 8 = 0.0833 pu
     Z_M1 = 0.0833 + j0.6667 pu = 0.6719 /_ 82.87 deg pu
     Y_M1 = 1 / Z_M1 = 0.185 - j1.477 pu

Step 2: Calculate Combined Thevenin Admittance and Impedance at Fault Bus
  Total Parallel Admittance:
  Y_th = Y_grid + Y_G1 + Y_M1
       = (1.209 + 0.133 + 0.185) - j(13.279 + 3.328 + 1.477)
       = 1.527 - j18.084 pu
  |Y_th| = sqrt(1.527^2 + 18.084^2) = 18.148 pu  (angle = -85.17 deg)

  Thevenin Equivalent Impedance:
  Z_1,th = 1 / Y_th = (1.527 + j18.084) / (18.148^2)
         = 0.004636 + j0.054907 pu = 0.05510 /_ 85.17 deg pu
  
  Effective System X/R Ratio:
  (X/R)_eff = 0.054907 / 0.004636 = 11.84

Step 3: Determine Symmetrical Subtransient Short-Circuit Current
  Per-unit Symmetrical Current:
  I_sym,pu'' = V_F / |Z_1,th| = 1.00 / 0.05510 = 18.148 pu
  
  Physical Symmetrical Current:
  I_sym'' = 18.148 * 1,387.86 A = 25,187 A = 25.19 kA

Step 4: Determine Short-Circuit MVA (SCMVA)
  SCMVA = S_base / |Z_1,th,pu| = 10 MVA / 0.05510 = 181.48 MVA
  Check: sqrt(3) * 4.16 kV * 25.187 kA = 181.48 MVA  (Exact match)

Step 5: Compute Peak Instantaneous Asymmetrical Current (t = 8.33 ms)
  Time constant tau = (X/R) / (2 * pi * f) = 11.84 / (2 * pi * 60) = 0.03141 s = 31.41 ms
  Peak multiplier at first half-cycle (t = 1/120 s = 8.333 ms):
  k_peak = sqrt(2) * [ 1 + exp(-8.333 ms / 31.41 ms) ]
         = 1.4142 * [ 1 + exp(-0.2653) ]
         = 1.4142 * [ 1 + 0.7670 ] = 1.4142 * 1.7670 = 2.499
  
  I_peak = k_peak * I_sym'' = 2.499 * 25.187 kA = 62.94 kA peak

Step 6: Compute Circuit Breaker Interrupting Duty at 3 Cycles (t = 50 ms)
  DC offset current at contact parting (t = 50 ms):
  i_dc(50 ms) = sqrt(2) * I_sym'' * exp(-50 ms / 31.41 ms)
              = 1.4142 * 25.187 kA * exp(-1.5918)
              = 35.620 kA * 0.2036 = 7.25 kA
  
  Total Asymmetrical RMS Current at 3 Cycles:
  I_asym,RMS(3 cycles) = sqrt( (I_sym'')^2 + (i_dc)^2 )
                       = sqrt( 25.187^2 + 7.25^2 ) = sqrt( 634.38 + 52.56 )
                       = sqrt( 686.94 ) = 26.21 kA RMS
  
  Asymmetrical Multiplying Factor: MF = 26.21 / 25.19 = 1.040
===================================================================================

---

7. Common Exam Traps & Strategic Pitfalls

• Omitting Base Power Conversions for Generator and Motor Reactances: Directly plugging in manufacturer nameplate reactances (given on the machine's individual MVA base) without converting to the system study base ($S_{base}$). Always apply $Z_{new} = Z_{old} \cdot \frac{S_{base,new}}{S_{base,old}} \cdot \left(\frac{V_{base,old}}{V_{base,new}}\right)^2$.
• Treating Induction Motors as Permanent Sources: Including induction motor contributions in steady-state ($X_d$) or long-time backup protection studies. Induction motors possess no internal field excitation; their backfeed vanishes completely within $1-4\text{ cycles}$.
• Confusing Peak Current ($I_{peak}$) with First-Cycle Asymmetrical RMS ($I_{asym}$): $I_{peak}$ is the instantaneous crest value (factor of $\sim 2.5 - 2.7 \times I_{sym}$), whereas first-cycle asymmetrical RMS is the heating/effective value (factor of $\sim 1.6 \times I_{sym}$).