The closed‑form solution for stresses and displacements in a homogeneous, isotropic, linear elastic semi‑infinite half‑space produced by a concentrated normal point load applied at the surface; used to estimate the spatial distribution of elastic stress increments and surface/subsurface settlement under localized loads. Assumptions: linear elasticity, isotropy, homogeneity, semi‑infinite domain, small strains, static loading, and a point (or approximated concentrated) load at the surface.

Principle

Under the stated assumptions the vertical stress and displacement fields are determined uniquely by the point load and decay with distance; stress components scale with the load and fall off with depth and radial distance according to the analytic Green’s‑function for a half‑space (allowing distributed loads to be obtained by superposition/integration).

Demonstration

Illustrative scenario: A point load P is applied at the ground surface. The vertical stress at depth z and radial distance r is given (analytically) by σ_z(r,z)=3P z^3/[2π (r^2+z^2)^(5/2)]. An engineer integrates this kernel over a rectangular footing footprint to estimate stress increments beneath the foundation and then uses elastic moduli to estimate elastic settlement.

Misapplication

Applying the Boussinesq Solution directly to layered, anisotropic, highly non‑linear, or plastically yielding soils (e.g., active failure) assumes linear elastic response and homogeneity; the semantic error is treating an elastic Green’s function as valid where inelastic or boundary effects dominate.

Consequence

When correctly applied it provides a first‑order elastic estimate of stress changes and influence zones beneath surface loads and supports superposition for distributed loads; misapplication can under‑ or overestimate stresses and settlements, leading to unsafe designs or unnecessary conservatism.

Reversal

The solution no longer applies when: soils behave nonlinearly or plastically under the load; the medium has significant layering or anisotropy at the scale of interest; the load is dynamic or cyclic; or the domain is finite or constrained by nearby rigid bodies—then numerical or layered elastic/plastic models are required.

Boundary

Clearly within: a small surface load on a thick, uniform, elastic granular medium where strains remain small. Boundary case: a shallow layered profile with slightly different stiffnesses—elastic superposition may give rough guidance but accuracy is limited. Clearly outside: deeply layered, saturated undrained clay exhibiting plastic consolidation, or pile foundations transferring loads to discrete depths.

Semantic Tension

Simplicity and analytic clarity (elastic Green’s function) versus the need to represent real soil nonlinearity, layering, and finite geometry; choosing Boussinesq approximations trades model tractability for potential inaccuracy in complex ground conditions.

Synthesis

Boussinesq gives an exact elastic kernel for point‑load response that is invaluable for influence‑zone reasoning and superposition, but it is a mechanical idealization: use it to quantify elastic influence and to seed more realistic layered or nonlinear analyses rather than as a definitive predictor when inelastic or geometric effects matter.