Equilibrium

The ground is a mechanism

Bearing capacity is met as a formula with three terms and a table of coefficients, and that presentation hides what it is. Underneath is a plastic collapse mechanism — a rigid wedge, a fan of radial shear on a logarithmic spiral, and a passive wedge that has to be pushed up and out of the way — and every coefficient in the table is a property of that one drawing.

Assumes Two ways of being wrong, After the first yield, which is not the end and The load that depends on what carries it.

The bearing capacity of a footing is

qult=cNc+qNq+12γBNγq_{ult} = cN_c + qN_q + \tfrac{1}{2}\gamma BN_\gamma

and it is usually met as a formula with three terms and a table of coefficients. Read that way it is arbitrary: three products, three symbols nobody can derive, and a table to look them up in.

Read as a mechanism it is none of those things. Every coefficient is a property of one drawing, every angle in that drawing is a function of the friction angle alone, and the reason the three terms have different behaviours is that they are three different kinds of quantity.

A bearing capacity is a mechanism, and here it is. Prandtl's collapse mechanism under a 3.0 m footing in a soil of 32° friction. A rigid wedge is driven down with the footing at 61° to the horizontal; a fan of radial shear turns the stress through exactly ninety degrees on a logarithmic spiral whose growth rate is tanφ; and a passive wedge at 29° has to be pushed up and out of the way. Nothing here is empirical — every angle is a function of φ alone — and the mechanism reaches 15.9 m from the centre, which is 10.6 times the footing's half width. That is why two footings closer together than about four widths do not have separate bearing capacities.
Fig. 1 Prandtl’s collapse mechanism under a footing. A rigid wedge driven down at 45 + φ/2, a fan of radial shear turning the stress through exactly ninety degrees on a logarithmic spiral, and a passive wedge at 45 − φ/2 that has to be lifted out of the way.

Which free body produced the number

Three of them, in sequence, and the sequence is the mechanism.

Directly under the footing is a wedge of soil that goes down with it, rigid, bounded by two planes at 45+ϕ/245 + \phi/2 to the horizontal. It is at rest relative to the footing and it acts as part of it.

Outward from each edge of that wedge is a fan of radial shear. The soil here is at failure everywhere, and the direction of the major principal stress rotates through exactly ninety degrees between the two wedges. Because the soil obeys a friction law, the boundary of the fan is a logarithmic spiral r=r0eθtanϕr = r_0e^{\theta\tan\phi} — the curve whose tangent makes a constant angle ϕ\phi with its own radius, which is precisely the condition for the resultant on it to pass through the spiral’s centre.

Beyond the fan is a passive wedge, bounded by planes at 45ϕ/245 - \phi/2, which has to be pushed up and out to the surface for the footing to go down.

The whole thing is nine tenths geometry. The angles come from Mohr’s circle at failure and nothing else; the spiral’s growth rate is tanϕ\tan\phi; and the ninety degrees is the difference between the two wedge angles, which is (45+ϕ/2)((45ϕ/2))=90°(45 + \phi/2) - (-(45 - \phi/2)) = 90° at every friction angle. The fan is a quarter turn whatever the soil.

Where the coefficients come from

Follow the stress round the spiral and the ratio of the passive side to the active side comes out as

Nq=eπtanϕtan2 ⁣(45+ϕ2)N_q = e^{\pi\tan\phi}\tan^2\!\left(45 + \frac{\phi}{2}\right)

The exponential is the fan. The squared tangent is the two Rankine wedges either end of it. At ϕ=30°\phi = 30° it is 18.4, at 35° it is 33.3, at 40° it is 64.2 — an exponential in a variable engineers quote to the nearest degree, which is the first practical consequence of the derivation and not a small one.

Nc=(Nq1)/tanϕN_c = (N_q - 1)/\tan\phi is the same result written for a cohesion instead of a surcharge, and its limit as ϕ0\phi \to 0 is exactly 2+π2 + \pi.

Three coefficients, one exponential. The three bearing capacity coefficients against the soil's friction angle, on a logarithmic scale because they are exponentials: N_q is e^(π tanφ)·tan²(45 + φ/2) and the other two are written from it. At φ = 0 they are 5.14, 1 and 0 — the undrained case, where a footing carries 5.14 times the soil's shear strength and nothing else matters. At 32° they are 35.5, 23.2 and 30.2, and by 40° they have risen by another order of magnitude. Nothing on this plot is a fit except N_γ, which has no closed form: Vesic's expression gives 30.2 here where Meyerhof's gives 22.0 and Hansen's 20.8, a spread of tens of per cent inside a number quoted to three figures.
Fig. 2 The three coefficients against the friction angle, on a logarithmic scale because they are exponentials. Nothing on this plot is a fit except N_γ, and the three published expressions for it differ by tens of per cent.

NγN_\gamma is the odd one out and it is worth saying so plainly: it has no closed form. The mechanism above assumes weightless soil; putting the soil’s own weight inside the failing region makes the problem statically indeterminate in a way the spiral construction cannot handle, and every expression in use — Vesic’s 2(Nq+1)tanϕ2(N_q+1)\tan\phi, Meyerhof’s, Hansen’s — is a fit to a numerical solution. At ϕ=32°\phi = 32° they give 30.2, 22.0 and 26.9. A number quoted to three figures inside a formula whose authors disagree by a third.

Three terms, three kinds of quantity

The most useful thing the mechanism explains is why the three terms behave so differently, and the answer is dimensional.

The first two terms are stresses resisting stresses. A cohesion is a stress; a surcharge is a stress. Neither knows anything about how large the footing is, so neither term contains BB.

The third is a stress resisting a weight, and a weight inside the mechanism grows with the mechanism, which grows with BB. So the third term is proportional to the footing width — and the consequence is one of the strangest facts in the subject.

The only term with a width in it, and what it does. Ultimate bearing pressure against footing width, with the three terms drawn separately. The cohesion and surcharge terms are flat — 0 and 0 kN/m² here — because they are stresses resisting stresses and know nothing about the size of the footing. The weight term is ½γBN_γ and rises in proportion to B, so a wider footing carries a higher pressure as well as a larger load. Doubling the width of the footing drawn multiplies its allowable pressure by 2.00. And the intercept is the fact worth carrying: with no cohesion and no surcharge the line passes through the origin, so a footing of vanishing width on clean sand carries no pressure at all. Sand resting on sand holds nothing up.
Fig. 3 Ultimate pressure against footing width, with the three terms drawn separately. The flat part is cohesion and surcharge; the rising part is the weight term, and on clean sand with no surcharge the line passes through the origin.

A footing on clean dry sand at the surface has a bearing capacity of zero at zero width. Not a small capacity: none. Set a grain of sand on a bed of sand and it carries nothing, because there is no cohesion to hold the mechanism together and no depth of soil above the failure surface to weigh it down. The capacity is entirely bought with the size of the failure region, which is bought with the size of the footing.

That is why a small pad on sand is founded at a depth, and why a raft on sand is very much more efficient per unit area than a set of pads. Doubling the width of a footing on sand doubles the pressure it can carry as well as doubling the area — so it carries four times the load, not two.

It is also why the capacity is almost never what governs. A footing on sand wide enough for its capacity to be large is a footing whose settlement is the design criterion, and settlement is a stiffness problem rather than a strength one.

The building does not care how far it went down; it cares how much it tilted. Five footings on soil that is 45% as stiff under one of them, carrying 34 kN/m. They settle between 8 and 8 mm, and the number that matters is neither of those: it is the angular distortion between neighbours, 0.00 per thousand, or one in 8184909136 — against a limit of one in 500 for cracking in finishes, which this passes. A building that went down half a metre uniformly would be undamaged and would need a new front step; this one has moved a twentieth as far and has cracked.
Fig. 4 The other calculation the same ground needs. Bearing capacity is a collapse load and settlement is a serviceability displacement, and on almost every real foundation the second decides the size — and a soft patch under one support redistributes the frame’s moments long before anything approaches a mechanism.

The exponential, and what it does to a design

There is a consequence of the derivation that is worth a section on its own, because it decides how much site investigation a project should buy.

NqN_q is an exponential in tanϕ\tan\phi. Between 30° and 35° it rises from 18.4 to 33.3 — an 81% increase in capacity for a 5° change in a quantity that is measured indirectly, from a penetration test correlated to a friction angle by an empirical relation with a scatter of several degrees of its own.

So a foundation designed on ϕ=35°\phi = 35° where the soil is at 32° has about half the capacity assumed, and there is no partial factor anywhere in the calculation with a value of two. The response the profession has settled on is to apply the factor to tanϕ\tan\phi rather than to the capacity — which is the only place it does any good, because it is the only place where a small change is amplified.

That is a general property worth carrying beyond soil. Whenever a design quantity appears inside an exponential, the factor of safety belongs on the exponent and not on the answer, and a sensitivity study on the answer will systematically understate the risk. It is the same reasoning the load that never came near failing anything applies to a fatigue calculation, where the life goes as the cube of a stress range and the factor belongs on the range.

The bounds, which close

The undrained case is worth its own section, because it is one of a handful of problems in plasticity where the two bound theorems meet.

Set ϕ=0\phi = 0. The soil has a shear strength cuc_u and no friction; the spiral becomes a circle; the two wedges are at 45°.

Two ways of being wrong sets out the theorems: a mechanism gives an upper bound, because a mechanism that satisfies compatibility overestimates the collapse load, and a stress field in equilibrium that nowhere breaks the yield condition gives a lower bound.

The two theorems meet at 2 + π. Undrained bearing capacity factors from four calculations, for the one case in which the bounds close. A discontinuous stress field that nowhere breaks the yield condition gives N_c = 4 and is a lower bound: the footing certainly carries that. A circular slip through the footing edge is a mechanism and gives 2π = 6.28; optimising the circle brings it to 5.52; and Prandtl's wedge-fan-wedge mechanism gives 5.1416 = 2 + π. The mechanism cannot be improved and neither can the stress field that matches it, so the collapse load is not estimated here but known — which is why the undrained case is quoted to four figures and the drained one is quoted three different ways by three authors.
Fig. 5 Four calculations of the same undrained quantity, with this soil’s own drained factor beside them. A stress field gives 4, Prandtl’s mechanism gives 2 + π, and two circular slips give 5.52 and 2π — the mechanisms above and the stress field below, meeting at the exact answer, while the drained value sits six times higher and is bounded by nothing.

Four calculations, all of them classic:

A two-zone stress field, with a vertical stress under the footing and a horizontal one outside it, satisfies equilibrium everywhere and gives Nc=4N_c = 4. The footing certainly carries that.

A circular slip through the footing edge, centred on the edge, is a mechanism and gives Nc=2π=6.28N_c = 2\pi = 6.28.

Optimising the circle brings the mechanism down to 5.52.

Prandtl’s mechanism gives 2+π=5.14162 + \pi = 5.1416, and a stress field that matches it can be constructed. The mechanism cannot be improved and neither can the stress field, so the bounds have closed and the collapse load is not estimated but known.

That is why the undrained bearing capacity factor is quoted to four figures and the drained one is quoted three different ways by three authors. It is not that one soil is better understood; it is that one problem has been solved and the other has not.

There is a second reason the bounds are worth knowing about, and it is practical rather than aesthetic. Bounds tell a designer which way an approximation errs, and that is often more useful than the approximation itself. A circular slip is easy to compute, easy to check by hand, and easy to extend to a slope, a layered soil or an eccentric load — and because it is a mechanism it is always unsafe, so an answer from it is an upper bound and has to be treated as one. A stress field is harder to construct and always safe. Knowing which of the two a method is stops it being a number of unknown reliability and makes it a number with a direction.

The point the mechanism turns about is the same reasoning applied to a rigid-body collapse, and the same discipline: a mechanism gives a load, the load is an upper bound, and the way to improve it is to find a mechanism that gives a lower one.

Where the mechanism reaches, and what that means

The drawing has a width as well as angles. For ϕ=32°\phi = 32° and a 3 m footing the toe of the mechanism is about 13 metres from the centre — more than four times the footing width, on each side.

The angle of repose is where the demand crosses the coefficient. The friction a block on a plane demands of its contact, tan α, against the slope angle, with the coefficient μ = 0.62 drawn across it. The two cross at 31.8°, which is arctan μ and is the angle of repose: shallower than that and the demanded reaction is inside the cone, steeper and no reaction the contact can supply is. The whole calculation was run twice, at 100 and at 800, and both return 31.8° — the weight cancels out of tan α = F/N before the comparison is made, so it appears on neither axis and cannot move the crossing.
Fig. 6 Where the angles come from. A friction angle is the slope a granular material stands at, and it is the same angle that sets the wedge inclinations in the mechanism above.

Two consequences follow directly.

It also fixes the depth of ground that has to be investigated. A borehole taken to one footing width below founding level has not been taken through the mechanism, whose deepest point on the drawing is about a width below the underside and whose passive wedges reach four widths sideways. The free body is a choice decides, here, how deep to drill.

Two footings closer together than about four widths do not have separate bearing capacities. Their mechanisms overlap, and the honest calculation is one mechanism for both — which is generally better than two separate ones, because the passive wedges interfere, but which no formula covers.

A footing beside an excavation has lost most of it. The passive wedge on the open side has been dug away, and the loss is not proportional to how much was removed — the mechanism reorganises into a one-sided one whose capacity can be a third of the two-sided value. That is the arithmetic behind the rule that a trench must not be dug beside an existing foundation without checking, and it is the same reasoning as the load that depends on what carries it, where the pressure a wall receives depends on how far it has let the soil move.

And a footing near a slope has lost part of its mechanism. The passive wedge on the downhill side has to be pushed up into ground that is not there, so the capacity falls sharply — by half or more within a few widths of a crest, and by an amount that depends on the slope rather than on the soil.

Where the model stops

The soil was treated as rigid-plastic. It is not: it deforms a great deal before it reaches the mechanism, and the settlement that accompanies the last twenty per cent of the load is usually intolerable long before collapse. The mechanism is a real limit and is very rarely the design criterion.

The load was assumed vertical and central. An inclined or eccentric load changes the mechanism completely: the wedge is no longer symmetric, part of the footing lifts off, and the capacity is computed on an effective area centred on the resultant — which is the middle third applied to soil rather than to masonry.

One soil, one wall, and a factor of nine. The pressure on a 6 m wall retaining dry soil at 19 kN/m³ with a friction angle of 32°, in the three states Rankine's theory allows, drawn to one scale and with no surcharge so that the three thrusts stand in the ratio of the three coefficients exactly. Active is 0.307 and 105.1 kN/m; at rest is 0.470 and 160.8 kN/m; passive is 3.255 and 1113.1 kN/m. That is a spread of 10.59 from end to end, decided entirely by which way the wall moved and by how far — and the active and passive pair are exact reciprocals, Ka·Kp = 1.000. All three resultants sit at the same third point, 2.00 m above the base, because all three profiles are the same triangle scaled.
Fig. 7 The states the same soil can be in. Active, at-rest and passive are three different pressures from one material, decided by which way the wall has moved — and the passive wedge in the bearing mechanism is the third of them.

Water was ignored. A submerged soil has an effective unit weight of about half its bulk one, so the width term halves; and a rapidly loaded clay is undrained, which is a completely different calculation on the same ground.

The footing was assumed rigid and the pressure under it uniform. Neither is true: a rigid footing on clay has a pressure distribution peaked at its edges and a flexible one on sand has it peaked at the centre, and the difference decides the bending moment the footing itself has to carry. The beam that sits on the ground is that calculation, and it has nothing to do with the collapse mechanism above.

And the ground was assumed uniform. A stiff crust over soft clay, or a sand layer over peat, produces a mechanism that punches through rather than rotating — a different drawing with different arithmetic and no coefficients at all.

Why the mechanism is drawn at all

It is fair to ask why any of this matters, given that the mechanism is almost never the design criterion and the formula is available whether or not anybody understands it.

The answer is that the mechanism is what makes the formula extensible, and every real foundation is outside the case the formula covers. A footing beside a slope, a footing on a layered profile, a pair of footings close together, a footing carrying an inclined load, a footing on the crest of an excavation, a raft with a soft patch under one corner — none of them is in the table, all of them are ordinary, and each is a small modification to a drawing that anybody who has seen the drawing can make.

The formula, on its own, offers a set of correction factors: for shape, for depth, for inclination, for base tilt, for ground slope. Each is a number multiplying a term, each was derived from a modified mechanism by somebody else, and used without the picture they are a second table on top of the first. Drawing as calculation is the alternative, and here it is not a stylistic preference — it is the difference between a method that covers one case and a method that covers the cases that occur.

The generalisation

The idea to carry away is that a tabulated coefficient is usually a drawing that somebody has already made.

Three numbers looked up from a table are three numbers. The same three understood as properties of a wedge, a spiral and a passive block are a picture, and a picture can be interrogated: what happens near a slope, what happens with two footings, what happens if the soil is layered, what happens if the load is inclined. Every one of those questions is unanswerable from the table and straightforward from the drawing.

The second is about the bound theorems, and it is a rare piece of good news. Two ways of being wrong presents them as a pair of one-sided answers, which is how they usually behave. Here they close, and when bounds close the result stops being an engineering estimate and becomes a theorem — which is why 2+π2 + \pi appears in a geotechnical formula at all, and why it is the only number in this essay that will still be exactly right in a hundred years.

Named alongside this one

Essays reaching for the same objects. Nobody chose these; they are what the concept index makes visible.

The objects this essay names

Each one links to every other essay that touches it.

Bearing capacityCohesionCollapse mechanismEarth pressureFailure surfaceFree bodyFriction angleLower bound theoremPassive pressurePlastic mechanismSettlementSize effectSlip surfaceSurchargeUpper bound theorem