The role of shallow foundations—such as isolated, strip, and raft footings—is fundamental to structural engineering, as they transfer superstructural loads directly to the upper soil strata. The safety of any structure depends not only on the structural capacity of its reinforced concrete elements, but critically on the ability of the underlying soil to support these applied stresses without undergoing shear failure. In geotechnical engineering, evaluating bearing capacity represents the most vital limit state check to prevent general or local shear failure mechanisms. An inaccurate bearing capacity evaluation can lead to excessive or differential settlements, structural tilting, and ultimately, catastrophic structural failure.
Developed in 1943 by Karl Terzaghi, the founding father of modern soil mechanics, the classical bearing capacity theory remains the cornerstone of geotechnical design. Decades later, it continues to serve as the analytical foundation for shallow foundation design within modern design codes, including TBDY 2018 (Turkish Building Earthquake Code) and Eurocode 7 (EN 1997). By applying the superposition principle, Terzaghi integrated the contributions of soil cohesion, surcharge (embedment depth), and soil unit weight into a practical and reliable framework. Although contemporary codes have introduced additional correction factors and modern safety philosophies, the underlying physical mechanism still relies on the failure wedge model defined by Terzaghi.
In practical engineering and routine design, errors rarely stem from the mathematical structure of Terzaghi’s equation itself, but rather from the misinterpretation of input parameters and boundary conditions. Critical errors frequently arise from neglecting the impact of the ground water table on unit weight, failing to account for effective foundation dimensions ($B’, L’$) under eccentric loading, misinterpreting net versus gross stress considerations regarding backfill and basement conditions, or improperly averaging parameters in multi-layered soil profiles. Consequently, mastering geotechnical design requires far more than blindly applying formulas; it demands a thorough understanding of the underlying engineering mechanics and regulatory guidelines.
2. Regulatory Nuances: Backfill & Basement Conditions (Eurocode7 Guidelines)
Evaluating bearing capacity for shallow foundations requires careful attention to how backfill and excavation affect soil stresses, particularly when aligning design practices with Eurocode 7 standards. A common pitfall in geotechnical analysis is applying a uniform assumption to excavation weight and backfill, whereas post-construction operations completely alter the governing equilibrium equations and limit state verifications.
When isolated or continuous footings are constructed and subsequently backfilled, the soil weight removed during excavation is effectively reapplied to the system. Under Eurocode 7 principles, the design check accounts for the total applied bearing pressure, expressed as $q_a + q < q_t$, where $q_a$ represents the applied net foundation pressure at the base, $q$ is the overburden pressure ($\gamma \cdot D_f$), and $q_t$ is the characteristic design bearing resistance of the soil. In this scenario, excavation weight is not deducted from the bearing capacity; instead, the weight of the placed backfill resting on the footing must be included in the total base pressure calculation and checked against the total resistance.
Conversely, for basement structures where no permanent soil backfill is placed over the raft or foundation base, the soil excavated during construction provides permanent stress relief at the foundation level. In accordance with Eurocode 7 verification approaches, safety is evaluated directly using the net bearing capacity framework, formulated as $q_a < q_{tnet}$. Here, the net characteristic design resistance ($q_{tnet} = \frac{q_k – q}{\gamma_{Rv}}$) isolates the net structural action, subtracting the weight of the excavated soil from both the capacity and the foundation pressure balance. The net structural stress applied to the soil is thus strictly bounded by the net bearing capacity of the ground.
Software solutions like SETAF2018 streamline these regulatory nuances by automatically executing code-compliant checks without requiring manual intervention for boundary conditions. When assessing static ($1.4G + 1.6Q$) and dynamic ($G + Q + E$) load combinations under rigid foundation assumptions, the software dynamically identifies whether the structure includes a basement or utilizes isolated, continuous, or raft footings. It automatically computes net foundation pressure ($q_{net}$), accounts for surcharge effects, and verifies design resistance ($q_t$ or $q_{tnet}$) in full compliance with Eurocode 7 guidelines. This automated approach eliminates duplicate accounting of soil weights and prevents critical errors in limit state verification.
3. Core Formulation & Parameter Integration
One of the primary advantages of Terzaghi’s classical bearing capacity theory is its ability to combine soil strength parameters and foundation geometry within a single unified framework. While basic analytical calculations often concentrate solely on bearing capacity factors ($N_c, N_q, N_\gamma$), neglecting foundation shape factors ($s_c, s_q, s_\gamma$) can lead to substantial errors. Because failure surface development varies between continuous, rectangular, and circular geometry due to three-dimensional stress distribution effects, shape corrections directly dictate the ultimate bearing capacity. For continuous (strip) footings, shape factors are taken as unity since plane strain conditions prevail, whereas rectangular and circular footings benefit from three-dimensional load dispersion, which increases bearing resistance through specific geometry multipliers. Consequently, two foundations situated on identical soil profiles can yield significantly different bearing capacities based purely on their plan geometry.
The classical bearing capacity factors themselves depend strictly on the internal friction angle ($\phi$) of the soil. In granular soils, as $\phi$ increases, the shear resistance grows, causing $N_q$ and $N_\gamma$ to exhibit a steep, non-linear logarithmic increase. Because small variations in friction angles dramatically affect bearing capacity—especially at higher $\phi$ values—accurate determination of shear strength parameters from field or laboratory testing is essential. Conversely, in undrained cohesive soils under total stress analysis where $\phi = 0$, the formulation simplifies considerably: $N_q = 1$, $N_\gamma = 0$, and $N_c = 5.14$. In this undrained scenario, the third term of the equation vanishes completely, leaving ultimate bearing capacity dependent almost entirely on undrained shear strength ($c_u$) and surcharge overburden pressure. Selecting the correct drainage condition is therefore far more decisive than factor selection alone.
SETAF2018 automates the selection and computation of both bearing capacity factors and geometry correction factors based on user-defined soil parameters and foundation shapes. Whether analyzing strip, rectangular, or circular footings, the software applies the exact formulations required for drained or undrained conditions. This automated parameter integration eliminates manual table lookups, linear interpolation errors, and incorrect scenario selections, ensuring consistent and reliable results across complex foundation designs.
4. Impact of Ground Water Table (GWT) on Unit Weight ($\gamma$)
The ground water table (GWT) represents one of the most frequent sources of calculation error in Terzaghi’s bearing capacity equation. Water table location alters not only effective stresses within the soil mass, but also directly dictates the unit weight ($\gamma$) parameters applied across the bearing capacity terms. Incorrect unit weight selection—particularly within the soil self-weight term ($\frac{1}{2} \gamma B N_\gamma$)—can lead to dangerous overestimations or overly conservative underestimations of foundation capacity.
Geotechnical analysis evaluates the impact of the water table across three critical scenarios defined relative to foundation depth ($D_f$) and width ($B$):
- Case A ($z \ge B$ below base): When the water table is located at a depth greater than or equal to the foundation width $B$ below the foundation base, the shear failure zone develops entirely within unsaturated or moist soil. The natural moist unit weight ($\gamma$) is used throughout without modification.
- Case B ($0 < z \le D_f$ or $0 < z \le B$ below base): When the water table lies within the influence zone below the foundation base, the reduced effective unit weight must be accounted for within the failure wedge. As the water table approaches the foundation level, buoyant forces reduce the soil’s effective weight, requiring an averaged effective unit weight ($\gamma’$) based on the submerged portion.
- Case C (GWT within embedment depth $D_f$): When the water table rises above the foundation base level into the embedment zone, both the overburden surcharge term ($q = \gamma D_f$) and the foundation base soil weight are affected. Submerged soil zones must adopt the buoyant unit weight ($\gamma’ = \gamma_{sat} – \gamma_w$), directly lowering effective surcharge stress and significantly reducing ultimate bearing capacity.
In multi-layered soil profiles, determining the effective unit weight requires rigorous weighted-average calculations across different soil strata and groundwater boundaries. Manual calculations are time-consuming and prone to human error, particularly when assessing seasonal water table fluctuations.
SETAF2018 addresses this complexity through direct borehole log integration. By processing groundwater levels defined in borehole profiles, the software dynamically computes effective vertical stresses and adjusts surcharge and unit weight parameters across every impacted soil layer. This ensures that effective unit weight corrections and overburden stress calculations are applied accurately without requiring manual multi-layer averaging.
5. Eccentric Loading & Effective Foundation Dimensions ($B’, L’$)
Terzaghi’s classical bearing capacity theory assumes that structural loads are uniformly distributed and applied concentrically through the centroid of the foundation base. In practical structural design, however, columns and shear walls transfer overturning bending moments ($M_x, M_y$) alongside vertical loads ($P$). This eccentricity shifts the resultant load away from the foundation center, generating a non-uniform contact pressure distribution across the soil interface. Because soil possesses negligible tensile strength, large eccentricities can cause gapping or loss of contact along the footing edges. To maintain equilibrium under compressed zones, bearing capacity calculations must disregard uncompressed areas and rely exclusively on an effective reduced foundation area.
To account for eccentric loading, the physical dimensions of the foundation ($B \times L$) are reduced to effective dimensions ($B’ \times L’$) centered directly under the load resultant. The eccentricity in each principal direction is calculated as:
$$e_x = \frac{M_y}{P} \quad \text{and} \quad e_y = \frac{M_x}{P}$$
Correspondingly, the effective width and length are defined as:
$$B’ = B – 2e_x \quad \text{and} \quad L’ = L – 2e_y$$
In Terzaghi’s formulation, substituting $B’$ for $B$ in the soil self-weight term—and utilizing the effective area $A’ = B’ \cdot L’$ to evaluate average contact pressure—directly captures the capacity-reducing effects of load eccentricity. As eccentricity increases, the effective contact area shrinks, raising soil contact stresses while simultaneously lowering allowable bearing resistance.
While Terzaghi provides a fundamental framework for eccentric loads, advanced methodologies such as Meyerhof and Vesic extend this concept by introducing explicit correction factors for load inclination, base tilt, and ground slopes. SETAF2018 automates this entire pipeline by evaluating vertical load and moment components ($P, M_x, M_y$) to compute bi-axial eccentricities ($e_x, e_y$) and effective dimensions ($B’, L’$) instantly. Furthermore, the software seamlessly supports advanced Meyerhof and Vesic formulations, automatically adjusting shape, depth, and load inclination factors to eliminate manual computation errors in complex loading environments.
6. Engineering Software Solution: SETAF2018 Overview
Modern geotechnical design requires balancing complex analytical equations—such as groundwater adjustments, multi-layer stress distribution, and bi-axial eccentricity—with strict regulatory compliance. Traditional workflows often force engineering firms to rely on manual spreadsheets prone to human error or fragment their budget across costly foreign software suites where individual modules for footing design, slope stability, and liquefaction are sold separately at exorbitant prices.
SETAF2018 addresses these challenges by delivering an all-in-one, integrated geotechnical engineering suite. Fully aligned with both TBDY 2018 (Turkish Building Earthquake Code) and Eurocode 7 standards, the software automates the end-to-end bearing capacity workflow. From dynamic borehole integration and automatic groundwater corrections to multi-combination limit state checks ($1.4G + 1.6Q$ and $G + Q + E$), SETAF2018 eliminates repetitive manual calculations and ensures reliable, code-compliant designs within a single unified workspace.
Back-office efficiency and regulatory confidence should not come with a prohibitive price tag. With an intuitive user interface that takes only hours to master, comprehensive video training libraries, and dedicated technical support, SETAF2018 provides a complete engineering package built for modern design offices.
7. Frequently Asked Questions (FAQ)
Which soil types are suitable for Terzaghi’s bearing capacity analysis?
Terzaghi’s bearing capacity theory is primarily formulated for homogeneous, uniform soil layers beneath shallow foundations under general shear failure conditions. While ideal for dense sands, stiff clays, and well-compacted gravels, its basic framework can also be adapted to soft clays or loose sands by applying local shear failure correction factors (modifying cohesion $c’ = \frac{2}{3}c$ and friction angle $\tan \phi’ = \frac{2}{3}\tan \phi$). For highly non-homogeneous or multi-layered profiles, modified parameter averaging techniques or numerical limit equilibrium approaches are recommended.
How do Terzaghi, Meyerhof, and Vesic methods differ?
While Terzaghi established the classical three-term superposition model ($c$, $q$, $\gamma$), later theories refined his assumptions to handle more complex field conditions:
- Terzaghi (1943): Assumes vertical concentric loads on a continuous shallow footing with a horizontal base and ground surface. Shape factors ($s_c, s_q, s_\gamma$) are applied for square or circular geometries.
- Meyerhof (1963): Extends the failure surface above the base level into the overburden zone and introduces explicit depth factors ($d_i$) and load inclination factors ($i_i$) to account for inclined and eccentric loads.
- Vesic (1973): Incorporates advanced parameters for ground slope, tilting foundation bases, and soil compressibility (rigidity index), offering a more comprehensive failure wedge geometry particularly suitable for deep or inclined shallow foundations.
Does a rising water table always decrease bearing capacity?
Yes, in almost all scenarios. A rising water table increases pore water pressure and reduces effective stress within the soil matrix. When groundwater rises above the base level or enters the failure zone beneath the foundation (within depth $B$), the submerged unit weight ($\gamma’ = \gamma_{sat} – \gamma_w$) replaces the moist unit weight ($\gamma$). Because $\gamma’$ is roughly half of $\gamma$, the self-weight term ($\frac{1}{2} \gamma B N_\gamma$) and surcharge overburden term ($q = \gamma D_f$) are significantly diminished, reducing the overall bearing capacity by up to 50%.
Why is load eccentricity critical in footing design?
Overturning moments ($M_x, M_y$) caused by wind, seismic forces, or eccentric column placement cause non-uniform contact stress distributions at the foundation base. Since soil cannot sustain tension, large eccentricities reduce the compression zone, leading to edge gapping or localized bearing failure. Calculating reduced effective dimensions ($B’ = B – 2e_x$ and $L’ = L – 2e_y$) ensures that bearing capacity checks are evaluated strictly over the compressed soil area, preventing unsafe overestimations of contact resistance.
How should multi-layered soil profiles be evaluated?
In multi-layered soil profiles, applying parameters from a single layer across the entire depth introduces severe calculation errors. Evaluating layered soils requires determining the influence depth of the failure surface beneath the footing (typically $1.5B$ to $2.0B$). If a strong soil layer overlies a weak sublayer, punching shear failure into the softer layer must be checked. Conversely, for gradually changing strata, weighted average shear strength parameters ($c_{avg}, \phi_{avg}$) and effective unit weights ($\gamma’_{avg}$) must be computed across the active failure zone, or automated multi-layer analysis algorithms should be utilized.
