Substation Grounding Grid: Touch, Step, GPR and Conductor Size
Calculate substation grid resistance, GPR, mesh and step voltages, tolerable limits and copper conductor size. Compare grid spacings using a uniform-soil IEEE 80 method.
A substation grounding grid must control the voltage a person can bridge during a fault as well as carry fault current without unacceptable heating. This calculator supports preliminary design and independent hand checks for substation designers, electrical engineers and students. It estimates grounding grid resistance, ground potential rise (GPR), corner-mesh touch voltage, one-metre step voltage and the minimum annealed-copper conductor area for a specified thermal duty.
Start with square mesh spacing and the number of meshes along each side. These define a rectangular footprint: ten meshes at 7 m give a 70 m side and eleven conductor runs. Enter burial depth, the actual conductor outside diameter and any bonded perimeter rods. The outside diameter controls the voltage geometry; copper metal area controls heating. Those dimensions differ for a stranded conductor, so both are entered separately.
Use measured soil resistivity and the wet resistivity and thickness of the surface covering. This foot-contact-only model requires 0.08–0.15 m of covering and a covering thickness no greater than one quarter of burial depth; otherwise use layered-soil analysis. Select bare soil where a protective covering cannot be assured. The model applies the covering to both touch and step contact locations, including outside the perimeter. Choose the 50 kg or 70 kg body-weight criterion and an assessed shock duration. The 70 kg option requires an appropriate restricted-access design basis; it is not a way to relax a failed public-exposure assessment.
The fault-study current feeding remote earth is multiplied by the current division and decrement factors. This earth current determines mesh voltage, step voltage and GPR. Enter the conductor’s own worst equivalent rms fault current separately. Local transformer contributions, metallic returns and individual conductor paths can make its heating duty different from the earth-injected current. Include growth, dc offset and the relevant clearing sequence in the current and duration inputs.
Compare calculated mesh voltage with tolerable touch voltage, and calculated step voltage with tolerable step voltage. The copper result is a thermal minimum before allowances for mechanical strength, corrosion and connection requirements. A PASS requires the model coverage screens and these three numerical comparisons. Low grid resistance alone does not establish acceptable touch and step voltages.
The spacing plots keep the present footprint, soil, rods and fault duty fixed while varying conductor spacing. They interpolate the conductor counts continuously to show trends. Green denotes supported voltage utilization at or below one; red denotes an exceeded limit. Invalid samples are gaps, shown as gray spans in the utilization plot, with an omitted-sample count. Convert a promising spacing into whole mesh counts and check that actual layout again. Changing the spacing input itself also changes the footprint because the mesh counts stay fixed.
Sverak is the primary resistance estimate. Laurent–Niemann and an approximate Schwarz calculation provide comparisons, not a validated uncertainty band. The tool does not model layered soil, unequal meshes, transferred potentials, fences, lightning, very large-grid inductance or detailed surface potential fields. It follows the IEEE 80-2013 method basis; check against the standard in force. Final design requires a qualified engineer, and installation work belongs to qualified persons.
How it is calculated
Geometry and earth current
The grid has equally spaced square meshes inside a rectangle, with perimeter conductors included: , and . Optional rods are bonded at corners and distributed around the perimeter. The earth current is ; the conductor heating current is a separate fault-study input.
Resistance and GPR
Sverak is used for the headline resistance and GPR:
Laurent–Niemann adds the finite-length term to the equivalent surface-plate term . Open background: Balakrishnan (1990), printed pp. 7–10; published resistance checks: Silva et al. (2017), PDF p. 3.
The Schwarz comparison combines grid, rod-bed and mutual resistances using IEEE 80-2013, §14.3. Its shallow-grid coefficients are the explicit approximation in Perng et al. (2018), PDF p. 5: and , where is the longer-to-shorter side ratio. They do not replace the edition’s depth-dependent coefficient assessment. The rod logarithm follows IEEE Eq. (60), with rod radius defined explicitly; it differs from Perng Eq. (9). No rods returns the horizontal-grid resistance without evaluating a singular coupled expression.
Touch and step criteria
The geometry factors follow IEEE 80-2013, §16.5, checked against the licensed edition; the open development of the revised irregularity factor and geometry count is in Balakrishnan, printed pp. 29–32. For a rectangle, . The displayed steps expose , , , and the effective lengths. with perimeter rods, otherwise .
The surface factor uses a 0.09 m empirical length and reduces to one for bare soil. With , the body circuit uses 1000 Ω and contact coefficients 1.5/m for touch and 6/m for step. The body coefficients are 0.116 and 0.157 A·√s for the two weight criteria (Perng, PDF p. 3). Retain full precision; published IEEE examples round intermediate factors.
Copper thermal area
An adiabatic energy balance integrates temperature-dependent resistivity:
The result is converted explicitly from SI area to mm². Annealed-copper constants follow the IEEE 80-2013 thermal model; the displayed steps state each constant and its reference temperature. The independent KME manual, PDF p. 37, provides a rounded copper short-circuit-factor check.
Coverage of the spacing study
The selected layout and every sweep point are screened against the model envelope before a voltage comparison is accepted. The screen uses the benchmark coverage discussed in IEEE 80-2013, §16.7, together with this tool's narrower spacing range: area 6.25–10,000 m², one to forty meshes per side for the continuous trend, valid burial/diameter geometry, perimeter-rod node capacity and a strictly positive mesh factor. Invalid voltage samples are gaps, never zero or negative voltage estimates. Gray spans and the omitted-sample count explain these gaps. Passing this coverage screen is not a numerical accuracy guarantee.
The foot-contact-only covering model is restricted here to 0.08–0.15 m thickness, consistent with the usual surfacing band discussed in IEEE 80-2013, §§7.4 and 12.5. This tool also requires the covering to be no thicker than one quarter of burial depth. That additional ratio is a conservative scope screen adopted here, not an IEEE acceptance criterion or a quantified error bound. Outside either screen, use layered-soil analysis rather than crediting the covering only in the foot-contact circuit.
Assumptions
- Uniform, linear-resistivity soil at 50/60 Hz; negligible voltage drop along the bonded grid. Only the screened thin covering modifies foot contact; its effect on bulk grid resistance and surface potentials is neglected.
- Equal square meshes; a closed rectangular perimeter; rods of equal length bonded to perimeter nodes, including corners. Rod symbols in the plan are illustrative and are not an installation layout.
- The same wet covering exists at every assessed foot location, including at least the one-metre step outside the corner. Choose bare soil if that covering is not assured.
- Body weight is a choice of published criteria, not a continuous person-specific safety prediction. Footwear resistance is not credited.
- Fault currents include the assessed system growth. The user supplies the division/decrement factors and the equivalent rms conductor duty; this tool does not solve the network fault or protection study.
- Annealed copper with constant volumetric heat capacity and linear resistivity versus temperature. No heat loss, current-sharing credit, corrosion allowance or joint derating is included.
Where it applies
The method basis is IEEE Std 80-2013, published in 2015 with Cor 1-2015; check against the standard in force. The licensed edition was consulted for equation and worked-example verification; its text, tables and drawings are not reproduced.
Tool limits: 2.5–15 m spacing; 2–40 whole meshes per side for the selected layout; grid area no greater than 10,000 m²; 0.26–2.49 m burial; conductor diameter below one quarter of burial depth; rectangular aspect ratio at most 8:1; shock and thermal durations 0.03–3 s. Every spacing sample must separately satisfy the coverage screen, including a positive mesh factor and enough perimeter nodes for the fixed rods. The trend can include fractional mesh counts between one and forty, following the geometric envelope of IEEE 80-2013 §16.7. Its benchmark coverage is a boundary for this tool, not proof that every combination inside the envelope has the same accuracy.
Credited surface covering must be 0.08–0.15 m thick, at least as resistive as the underlying soil, and no thicker than one quarter of grid burial depth. The thickness band follows the usual surfacing discussed in §§7.4 and 12.5; the depth-ratio limit is an additional conservative screen adopted by this tool, not a universal IEEE thin-layer criterion. Merely putting the conductor below the covering is insufficient. Use layered-soil analysis outside these bounds; selecting bare soil is appropriate only if no covering is being credited.
For the approximate Schwarz coefficients this tool additionally requires and physically consistent self/mutual resistances. These extra restrictions are scope limits of this implementation, not general acceptance limits in IEEE 80.
Layered soil, irregular or unequal grids, complex rod locations, large interconnected stations and appreciable conductor inductance need detailed modelling. Check fences, accessible metal, metal-to-metal touch, transferred potentials, cable screens, remote bonds, protection and reclosing separately. Verify wet/seasonal soil and covering properties by site data.
The chart is a continuous analytical spacing study at fixed footprint and fixed rods, not an optimization or a buildable layout. Intermediate line counts can be fractional. Invalid samples are omitted; gray spans and the omission count identify portions outside the model. Green/red regions assess only supported mesh and step voltage samples. The dotted line marks the selected spacing. The native chart converts its metre/volt axis labels with the displayed units; schematic labels remain explicitly SI.
The PASS card requires the geometry and covering screens as well as voltage and copper thermal comparisons. It makes no claim about omitted design duties. Copper final temperature must respect the weakest conductor, connection and surrounding material.
A PASS here does not approve a substation design.
References
- IEEE Std 80-2013, IEEE Guide for Safety in AC Substation Grounding (2015 publication, incorporating Cor 1-2015). DOI: https://doi.org/10.1109/IEEESTD.2015.7109078. Method basis and licensed equation/worked-example verification; check against the standard in force.
- Silva et al. (2017), Calculation of Grounding Grids Parameter at Arbitrary Geometry, Transactions on Environment and Electrical Engineering 2(2). Open PDF, especially p. 3: https://distantreader.org/stacks/journals/teee/teee-95.pdf
- Perng, Kuo and Lu (2018), Grounding System Cost Analysis Using Optimization Algorithms, Energies 11, 3484. DOI: https://doi.org/10.3390/en11123484. Open equations for body current, Schwarz coefficients and Sverak resistance; PDF pp. 3–5.
- Balakrishnan (1990), Calculation of important design parameters for grounding systems in substations, Montana State University thesis. Open repository: https://scholarworks.montana.edu/server/api/core/bitstreams/81e41043-1c23-4936-8918-f92b708b5e80/content. Printed pp. 7–10 and 29–32.
- KME Italy, Mineral Insulated Cables MICO and their Terminations, instruction manual rev. 04, PDF p. 37 (printed p. 29). Independent rounded copper short-circuit factor: https://www.kme.com/fileadmin/DOWNLOADCENTER/COPPER%20DIVISON/6._Mineral_Insulated_Cables/Ueberarbeitet_MICO_ohne_Register_10.2022/KME_Italy_SpA_-_Instruction_Manual_-_Mineral_Insulated_Cables_MCO_and_their_terminations_rev_04.pdf
Inputs
| Input | Symbol | Unit | Default |
|---|---|---|---|
| Square mesh spacing Both directions use this spacing. Grid dimensions equal spacing × mesh counts. | D | m | 7 |
| Meshes along x 2–40 square meshes per side in this tool; total grid area must not exceed 10,000 m². | N_x | 10 | |
| Meshes along y 2–40 square meshes per side in this tool; total grid area must not exceed 10,000 m². | N_y | 10 | |
| Conductor burial depth | h | m | 0.5 |
| Installed conductor outside diameter Use the actual outside diameter for the voltage geometry, including stranding. | d | mm | 14 |
| Perimeter ground rods Use zero for no rods. Otherwise at least four, including the corners, bonded to the grid and distributed around its perimeter. | N_r | 20 | |
| Length of each rod | L_r | m | 7.5 |
| Rod diameter | d_r | mm | 16 |
| Uniform soil resistivity Use measured, seasonally conservative soil data. This is a single uniform-soil model. | rho | ohm*m | 100 |
| Surface covering | surface | High-resistivity surface layer | |
| Wet surface-layer resistivity | rho_s | ohm*m | 3000 |
| Surface-layer thickness Foot-contact-only model: 0.08–0.15 m and at most one quarter of burial depth. This extra depth-ratio screen is a tool limit; otherwise use layered-soil analysis. | h_s | m | 0.1 |
| Symmetrical fault current feeding remote earth Use the fault study contribution relevant to earth return. Local circulating transformer contributions may differ. | I_f | kA | 20 |
| Current division factor Fraction of the above fault current entering earth through this grid, after metallic return paths. | S_f | 0.3 | |
| Grid-current decrement factor From the fault study; includes the rms effect of dc offset over the fault duration. Growth allowance must already be included in the fault currents. | D_f | 1.1 | |
| Shock duration Use the assessed exposure/clearing time, including backup protection and reclosing where relevant. | t_s | s | 0.5 |
| Body-weight criterion Use 50 kg for public exposure unless the governing design basis establishes otherwise. No interpolation by actual weight. | body | 50 kg (110 lb) | |
| Conductor design current, equivalent rms Worst current through the conductor being sized, including dc-offset and growth allowances. Do not automatically use the split earth current. | I_c | kA | 22 |
| Conductor fault duration | t_c | s | 0.5 |
| Initial copper temperature | T_a | °C | 40 |
| Allowable final copper temperature Use the lowest allowable value for the conductor, joints and adjacent materials. This is an annealed-copper adiabatic calculation. | T_m | °C | 250 |
| Installed copper metal area Actual metal area from the conductor specification; separate from outside diameter for stranded conductors. | A_inst | mm^2 | 120 |
Results
With the default inputs above, the calc gives:
| Result | Value |
|---|---|
| Calculated voltage and copper criteria | PASS |
| Grid resistance — Sverak | 0.688 Ω |
| Ground potential rise — Sverak | 4,542 V |
| Maximum current entering earth | 6,600 A |
| Calculated mesh touch voltage | 602.8 V |
| Tolerable touch voltage | 680.8 V |
| Calculated one-metre step voltage | 474.9 V |
| Tolerable step voltage | 2,231.1 V |
| Minimum copper area — thermal duty only | 91.8 mm² |
| Installed copper current capacity — adiabatic | 28.75 kA |
Checked against sources
Each case below was calculated by this calc and agrees with the value from its source (within the tolerance where one is shown).
Published 20 m square, four meshes
Source: Silva et al. (2017), Table I, PDF p. 3: 20 m × 20 m, four square meshes, 100 Ω·m, 0.5 m burial. Published Laurent and Sverak entries are used. The paper does not state its Schwarz coefficients; its Schwarz result is therefore not used as an exact target for the explicit Perng coefficient approximation used here.
Inputs: Square mesh spacing = 10 m; Meshes along x = 2; Meshes along y = 2; Installed conductor outside diameter = 10 mm; Perimeter ground rods = 0; Installed copper metal area = 70 mm²; the others as above.
| Result | Source | This calc |
|---|---|---|
| Grid resistance — Laurent–Niemann | 3.049 Ω | 3.049 Ω |
| Grid resistance — Sverak | 2.957 Ω | 2.957 Ω |
Published 20 m square, sixteen meshes
Source: Silva et al. (2017), Table I, PDF p. 3: 20 m × 20 m with sixteen square meshes. Independent published resistance values; same soil and burial as the four-mesh case.
Inputs: Square mesh spacing = 5 m; Meshes along x = 4; Meshes along y = 4; Installed conductor outside diameter = 10 mm; Perimeter ground rods = 0; Installed copper metal area = 70 mm²; the others as above.
| Result | Source | This calc |
|---|---|---|
| Grid resistance — Laurent–Niemann | 2.716 Ω | 2.716 Ω |
| Grid resistance — Sverak | 2.624 Ω | 2.624 Ω |
IEEE 80-2013 worked example B.1: no rods
Source: IEEE Std 80-2013, Annex B.1, DOI 10.1109/IEEESTD.2015.7109078. Worked-example results checked against the licensed edition. The printed example rounds Rg, Km and Cs before evaluating voltages. Full precision here differs by 0.0044 Ω, 0.49 V mesh and 2.35 V tolerable touch; absolute tolerances of 0.005 Ω, 0.6 V and 2.5 V cover that rounding.
Inputs: Installed conductor outside diameter = 10 mm; Perimeter ground rods = 0; Uniform soil resistivity = 400 Ω·m; Wet surface-layer resistivity = 2,500 Ω·m; Surface-layer thickness = 0.102 m; Symmetrical fault current feeding remote earth = 3.18 kA; Current division factor = 0.6; Grid-current decrement factor = 1; Body-weight criterion = 70 kg (154 lb), restricted access; Installed copper metal area = 70 mm²; the others as above.
| Result | Source | This calc |
|---|---|---|
| Grid resistance — Sverak | 2.78 ± 0.005 Ω | 2.776 Ω |
| Calculated mesh touch voltage | 1,002 V | 1,002 V |
| Tolerable touch voltage | 838.2 ± 2.5 V | 840.5 V |
IEEE 80-2013 worked example B.2: perimeter rods
Source: IEEE Std 80-2013, Annex B.2, DOI 10.1109/IEEESTD.2015.7109078. Worked example with twenty 7.5 m perimeter rods, checked against the licensed edition. The printed example uses rounded geometry factors; full precision differs by 1.66 V mesh and 0.22 V step. Absolute tolerances are 0.005 Ω, 1.8 V mesh and 0.3 V step.
Inputs: Installed conductor outside diameter = 10 mm; Uniform soil resistivity = 400 Ω·m; Wet surface-layer resistivity = 2,500 Ω·m; Surface-layer thickness = 0.102 m; Symmetrical fault current feeding remote earth = 3.18 kA; Current division factor = 0.6; Grid-current decrement factor = 1; Body-weight criterion = 70 kg (154 lb), restricted access; Installed copper metal area = 70 mm²; the others as above.
| Result | Source | This calc |
|---|---|---|
| Grid resistance — Sverak | 2.75 ± 0.005 Ω | 2.753 Ω |
| Calculated mesh touch voltage | 747.4 ± 1.8 V | 749.1 V |
| Calculated one-metre step voltage | 548.9 V | 549.1 V |
Published 50 kg body-current coefficient at one second
Source: Perng et al. (2018), Eq. (3), PDF p. 3: body-current coefficient 0.116 A·√s; at one second this gives 0.116 A.
Inputs: Shock duration = 1 s; Installed copper metal area = 70 mm²; the others as above.
| Result | Source | This calc |
|---|---|---|
| Tolerable rms body current | 0.116 A | 0.116 A |
Independent copper short-circuit factor
Source: KME Italy, Mineral Insulated Cables and their Terminations, rev. 04, PDF p. 37 (printed p. 29): independently published copper factor 115 A·√s/mm² for a 70 °C to 160 °C excursion. The IEEE 80 annealed-copper constants retained here give 114.37; the 0.8 absolute tolerance allows the approximately 0.6% difference in rounded material/design conventions. This cross-check validates the thermal integral and SI-to-mm² conversion, not cable insulation selection.
Inputs: Initial copper temperature = 70 °C; Allowable final copper temperature = 160 °C; Conductor design current, equivalent rms = 10 kA; Conductor fault duration = 1 s; the others as above.
| Result | Source | This calc |
|---|---|---|
| Copper short-circuit factor | 115 ± 0.8 A·s⁰.5/mm² | 114.4 A·s⁰.5/mm² |
Related Resources
Related courses
IEC Electrical Schematics and Industrial Control Panels
Basic Electrical Controls & Diagram Reading

Responses