Internal forces

How far a wrong load reaches

Every figure in this collection applies a load as a point, a line or a uniform pressure, and no real load is any of those. The licence is Saint-Venant's, it is usually quoted as a principle, and it is really a statement about a wavelength.

Assumes What a cut reveals, and why it was there all along, The free body is a choice, and choosing it well is the whole skill and The section that cannot stay flat.

Every load in this collection is a lie of some kind. A wheel is drawn as a point and is really a contact patch. A wall’s weight is drawn as a uniform line and is really a stack of bricks. A bearing reaction is drawn as an arrow and is really a pressure distribution over a plate, which is itself an idealisation of what happens between two surfaces that are neither flat.

The reason any of it works is Saint-Venant’s principle, and the reason the principle works is not a theorem about equilibrium. It is a statement about a length.

Three ways to apply the same force, and one depth to forget the differenceThree end loads on a member 400 mm deep, all with the same resultant and the same moment: a point load, the same force spread over a fifth of the depth, and the same force split in two. What is plotted is the *difference* between each of them and the beam-theory answer — the self-equilibrating remainder — as a fraction of the mean stress. The point load starts at 20 times it and is under a tenth of it by 0.77 depths; all three are under one per cent by about 1.18. That distance is the licence every figure in this collection is drawn under, and the exact strip eigenvalue agrees with it: 2.106 + 1.125i, whose real part puts one per cent at 1.09 depths and whose imaginary part means the remainder changes sign on the way out, which no statement of the principle mentions.00.511.5200.20.40.60.81distance from the loaded end (depths)self-equilibrating remainder ÷ mean stressa point loadspread over h/5split in twoone depthsame resultant, same moment · one per cent left by 1.18 depths
Fig. 1 Three end loads on a member 400 mm deep, with the same resultant and the same moment, distributed three different ways. What is plotted is the difference between each of them and the beam-theory answer — the self-equilibrating remainder. The point load starts at twenty times the mean stress and is under a tenth of it three quarters of a depth in.

Which free body produced the number

Two statically equivalent load systems differ by a system with no resultant force and no resultant moment. Subtract one from the other and what is left is self-equilibrating, and the whole principle is a claim about what such a system does.

So take the simplest self-equilibrating load there is: a normal traction on the boundary of a half-plane varying as cos(2πy/λ)\cos(2\pi y/\lambda). It pushes as much as it pulls, over every wavelength, so its resultant is zero by construction.

The Airy stress function

ϕ=pα2(1+αx)eαxcos(αy),α=2πλ\phi = \frac{p}{\alpha^2}(1 + \alpha x)\,e^{-\alpha x}\cos(\alpha y), \qquad \alpha = \frac{2\pi}{\lambda}

satisfies the biharmonic equation and both boundary conditions — the applied normal traction, and no shear on the loaded face. Differentiate it twice and

σxx=p(1+αx)eαxcos(αy)\sigma_{xx} = -p\,(1 + \alpha x)\,e^{-\alpha x}\cos(\alpha y)

The decay length is the load’s own wavelength, and nothing else appears. No modulus, no Poisson’s ratio, no thickness, no size of the loaded region except through λ\lambda itself.

The decay length is the load's own wavelength, and nothing elseThe stress left from a self-equilibrating end load that varies as a cosine of wavelength λ, against distance in depths. The Airy function (A + Bx)e^(−αx)cos(αy) satisfies both boundary conditions and gives a factor of (1 + αx)e^(−αx) with α = 2π/λ, so the curves are the same curve stretched: at one wavelength in there is 1.4 per cent left, whatever λ was. No modulus, no Poisson's ratio and no thickness appears anywhere. A self-equilibrating load across a depth must change sign at least twice, so its slowest component has a wavelength of about the depth — which is the whole of Saint-Venant's principle, with a number in it.00.511.5200.20.40.60.81distance from the load (depths)stress remaining ÷ stress appliedλ = 0.25 depthsλ = 0.50 depthsλ = 1.00 depthsλ = 2.00 depthsat one wavelength: 1.36% · at two: 0.005%
Fig. 2 The factor (1+αx)eαx(1 + \alpha x)e^{-\alpha x} for four wavelengths. They are the same curve stretched, which is the content of the result: at one wavelength in, 1.36 per cent is left, whatever λ\lambda was.

Which converts the principle into arithmetic

A self-equilibrating traction across a depth hh must change sign at least twice — once to have no resultant force and once more to have no resultant moment. Its slowest Fourier component therefore has a wavelength of about hh.

At x=hx = h that component is down to 1.4 per cent, and every faster component is down further. One depth is the answer, and it is an answer rather than a rule of thumb: it is (1+2π)e2π(1 + 2\pi)e^{-2\pi}.

That single number is the licence under which every figure on this site is drawn. A point load is drawn as a point because a distance of one depth from it, the difference between the point and whatever the real contact was is under two per cent — and every stress the essay is about is quoted somewhere else.

The exact answer, which oscillates

The half-plane is not a strip, and the strip has its own decay rates: the Papkovich–Fadle eigenvalues, the roots of sin2λ±2λ=0\sin 2\lambda \pm 2\lambda = 0 for a strip with traction-free long edges.

They are complex. Solving the symmetric equation by Newton’s method in the complex plane gives 2.1062+1.1254i2.1062 + 1.1254i, and the antisymmetric one gives 3.7488+1.3843i3.7488 + 1.3843i.

The real part is the decay, and it agrees with the half-plane argument closely: measured from the centre of a strip of depth hh, the first symmetric mode loses 99 per cent of itself in 1.09 depths against the half-plane’s 1.06. Two quite different routes, one number.

The imaginary part is the interesting one, because nothing in any statement of the principle mentions it. The residual stress changes sign as it dies away. It is not a decaying bump; it is a decaying oscillation, with a half-period of about 1.4 depths — so the stress a per cent past the end of the disturbed region is on the other side of the beam-theory answer, and a designer who has computed a peak near a discontinuity has not found the only place where the section is not doing what the model says.

It is worth saying plainly that the principle was a conjecture for more than a century. Saint-Venant stated it in 1855 as an observation about torsion and bending of prisms; Boussinesq gave it its general form thirty years later; and nobody proved anything until the middle of the twentieth century, when the honest statements turned out to be about energy rather than about stress. The theorem that exists says the strain energy stored beyond a distance falls exponentially, at a rate belonging to the body’s geometry — which is exactly what the eigenvalue above is, and which is why a counterexample has to be a body with a second length scale rather than a body with an unusual load.

Three times the stress, and it does not matter how big the hole isThe hoop stress around a circular hole in a wide plate pulled at 100 N/mm², from Kirsch's exact solution. At the sides of the hole it is 3.0 times the applied stress — 300 N/mm² — and the factor is the same for a hole of any radius, because the radius cancels. At the top and bottom of the hole it is -1.0 times the applied stress, which is compression in a plate that nothing is pushing. The disturbance dies quickly: the stress is within 5% of the applied value by 3.5 hole radii, which is Saint-Venant's principle with a number on it.pulled at 100 N/mm², left and right300-100 — compressionhoop stress, tinted3.0× at the edgewithin 5% by 3.5 radiithe applied stressdistance from the centre, in hole radii12345
Fig. 3 Where a stress does not decay at all. A hole raises the stress by a factor that depends on its shape and not on its size, and the elevation dies away over a distance of the hole’s own dimension — the same wavelength argument with the loaded region being the hole rather than the end.

Three statically equivalent loads, and how long they stay different

The principle is a statement about a difference, so the cleanest test of it is to build three loads with the same resultant and the same moment and watch them converge.

A point load at mid-depth starts at twenty times the mean stress. The same force spread over a fifth of the depth starts at five times. The same force split into two starts at nearly ten. All three are under a tenth of the mean by about three quarters of a depth, and under a hundredth of it by about 1.2.

That is the whole of it, and it is worth noticing how fast it is. The distinction between a point load and a patch — which decides everything about the bearing, the plate, the local crushing and the reinforcement — has ceased to be measurable within one depth of the place where it was decided.

The corollary is the one every essay on connections in this collection depends on: the region where it is not yet true is exactly the region a connection lives in. A joint is not a point precisely because it occupies the length over which Saint-Venant has not yet arrived, and a strut-and-tie model is the tool for that region for the same reason.

Where a section exists, and where it does notThe same beam divided into the regions the two theories own. Within about one depth of a support, a concentrated load, a corner or an opening, the strain is not linear across the section and every calculation on this site that begins by choosing one is inapplicable — those are the D-regions, marked here. What is left between them is the B-region, where beam theory is exact enough to have been trusted for two centuries. On a beam this deep the D-regions are most of it, which is the practical reason the strut-and-tie model exists at all: 38% of this span is a region a section cannot describe.DBDBDshaded: one depth either side of every discontinuity, where no section describes the strain
Fig. 4 The regions the principle divides a member into. Everything more than about a depth from a load or a support is a B-region, where plane sections is close enough to true; everything within one is a D-region, where it is not. The boundary is a depth, and this essay is where the depth comes from.

What plane sections is really assuming

Plane sections stay plane is usually presented as a kinematic assumption about bending. It is better read as a consequence of what has just been derived.

Any deviation from a plane section is a strain distribution with no net force and no net moment — because the force and the moment are exactly what the plane part of the distribution carries. So a deviation from plane sections is a self-equilibrating system, and self-equilibrating systems decay over their own wavelength.

Plane sections is therefore not an assumption about materials or about bending. It is Saint-Venant’s principle, restated for the interior of a member. Which explains why it holds so well in the middle of a span and fails so completely at a support: nothing has changed about the material, and everything has changed about how far away the disturbance is.

Where plane sections stop staying planeStrain across a cut face at four span-to-depth ratios, with the straight line the theory assumes drawn faintly behind. For a slender beam the two coincide; for a beam as deep as its span the real distribution is nothing like a straight line, and beam theory has no claim on it.span ÷ depth = 2off by 19%span ÷ depth = 4plane sections holdspan ÷ depth = 8plane sections holdspan ÷ depth = 16plane sections holdthe assumption is the theory — everything else is arithmetic on top of it
Fig. 5 How good the assumption is, against how slender the member is. The two arguments are the same argument: a deep member is one where a large fraction of its length is inside a depth of one end or the other, and a slender one is mostly interior.

The load that does not decay

There is one self-equilibrating load in this collection that travels metres rather than depths, and it is the exception that says what kind of statement the principle is.

A bimoment on a thin-walled open section is four flange forces — two pushing, two pulling — arranged so that there is no axial force, no shear and no bending moment about either axis. Every resultant is zero. It is exactly the kind of load system Saint-Venant’s principle is about.

It decays over EIw/GJ\sqrt{EI_w/GJ}, which for the 400 mm section here is 1.5 metres, or 3.8 depths. On a longer, thinner section — a channel purlin, a cold-formed zed — the same expression gives ten depths or more, and on a member only a few metres long the bimoment applied at one end has not decayed by the time it reaches the other.

The one self-equilibrating load that does travelHow far a self-equilibrating disturbance reaches, for two kinds of it, against the depth of the member. An ordinary end load is gone within about a depth — the bars on the left, all of them about 1.1 depths long. A bimoment on a thin-walled open section is the exception: it has no force and no moment, nothing a resultant can see, and it decays over √(EI_w/GJ), which for the section drawn is 1.52 m — 3.8 depths, and 16 times as far as the principle would suggest. Saint-Venant's is a statement about a length scale rather than a theorem about statics, and the length scale is not always the section's.200 mm300 mm400 mm600 mm900 mm0246distance to one per cent (m)an ordinary end loada bimoment on an open sectionsection depth
Fig. 6 Two kinds of self-equilibrating disturbance, against the depth of the member. An ordinary end load is gone within about a depth for every section drawn. The bimoment is not, and the difference is not a matter of degree: they are governed by different length scales, and only one of them is the section’s own.

The reason is visible in the expression. EIw/GJ\sqrt{EI_w/GJ} contains the ratio of two stiffnesses — the section’s resistance to warping against its resistance to twisting — and an open thin-walled section is enormously stiffer against the first than the second. Saint-Venant’s principle assumes a body with only one length scale in it. A thin-walled member has at least two, and the second is not a dimension of the cross-section at all.

A longitudinal stress with nothing on any diagram to predict itWarping normal stress at a flange tip along the member, against the Saint-Venant shear stress, which is usually the only stress a torsion calculation produces. The bimoment reaches 3.469 kN·m² at the held section itself — a stress resultant in units no diagram carries, being four flange forces whose force, moment and torque resultants all vanish — and the longitudinal stress it raises is 84.7 N/mm², against a peak Saint-Venant shear of 42.4 N/mm². The stress is computed two ways that share no algebra, B·ω/I_w through the sectorial coordinate ω = 19200 mm² and M_f/Z_f through flange bending, and the two agree to 1e-14 N/mm². It lands at the flange corner where the bending stress is already highest, and nothing in a bending calculation knows it is there.012345020406080distance along the member (m)stress (N/mm²)B = −3.469 kN·m² → 84.7 N/mm²warping normalstress at a tipSaint-Venant shear42.4 N/mm²
Fig. 7 The stress resultant with no resultant. A bimoment appears on no diagram, contributes to no equilibrium equation, and is the reason the principle that licenses this whole collection has an exception with a name.

The two-to-one spread, which is this principle wearing a hard hat

Every design office rule about load dispersion is Saint-Venant’s principle with the arithmetic replaced by a slope.

A load on a bearing plate is taken to spread through the plate at 45 degrees; a load on a slab is taken to spread at two horizontal to one vertical; a patch load on a girder web is taken to be effective over a length that grows with the flange thickness and the fillet. None of those is derived anywhere. All of them are saying the same thing: the width over which a load is effective grows with distance from where it was applied, at a rate of order one.

The dispersion angle is that rate. A spread of 45 degrees means a load is distributed over a width of about 2x2x at a distance xx, which is the same statement as “the disturbance has a wavelength of about xx and has therefore decayed” — and the reason the rules use a slope rather than an exponential is that a slope can be drawn on a detail and an exponential cannot.

What the slope loses is the fact that the convergence is quick. A designer who spreads a load at 45 degrees through 200 mm of concrete has, according to the arithmetic above, already lost more than nine tenths of the difference between a point load and a patch. The rule is conservative in the direction it usually needs to be, and its conservatism is not where it appears to be — it is not in the angle, it is in treating a continuous decay as though the load became uniform at some particular depth and not before.

The bearing is one length and the web is loaded over anotherA load applied over a stiff bearing of 100 mm on the flange of a girder with a 1200 × 6 mm web. The flange bends under it and the yield lines that form spread the load along the web over 703 mm — 7.0 times the bearing, and 86% of the yield resistance is that spread rather than the bearing. The effective length is not a decision anybody made: it is what the flange's own bending stiffness against the web's own strength works out to.s_s = 100l_y = 703flangewebspread 7.0 times the bearingyield resistance 1498 kN · elastic critical 215 kNresistance 284 kN at a slenderness of 2.64
Fig. 8 The dispersion rule where the code writes it down. A patch load’s effective length is a length the web chooses rather than one the detail imposes, and the choosing is the same convergence this essay is about.

The principle decides how long a test specimen has to be

Saint-Venant’s is usually met as a licence to simplify a calculation. Its most practical use is in the opposite direction: it says how much material an experiment needs before it is measuring anything.

A tensile coupon is gripped at both ends, and the grips apply a load distribution nothing like the uniform stress the test is about. So the gauge length has to sit far enough from both grips for the difference to have decayed — one specimen width is the usual rule, and it comes from the arithmetic above rather than from custom. A short coupon does not give a wrong number with a wide error bar; it gives a confident number about the grip.

The same reasoning fixes the length of a column test, the position of strain gauges on a beam specimen, and the distance from a weld at which a residual stress is quoted. In each of them the question is not what the load is but how far away the measurement has to be before the load’s shape stops mattering — and the answer is a section dimension, for the reason this page derives.

Two materials pulled until they stopTwo stress-strain curves — mild steel, aluminium alloy — plotted to a strain of 2.0%. One of them has a plateau, so the stress at which yielding starts is something the specimen does rather than something anyone chooses. No offset construction is drawn.00.5%1%2%2%050100150200250300350strainstress, N/mm²mild steelaluminium
Fig. 9 Where the numbers under every material property come from. Each is a measurement made at a gauge length chosen so that the grips have stopped being visible, which is Saint-Venant’s principle used as a design rule for an experiment.

Where the model stops

The half-plane solution is exact for a half-plane and an approximation for a strip. It gets the mechanism and the length scale right and does not know that the strip’s far face is free, which is what makes the exact eigenvalues complex.

The Fourier decomposition treats the end load as a boundary traction on an elastic half-space. A real end load is delivered through a plate, a weld or a bearing, whose own stiffness redistributes it before the member ever sees it — so the “point load” case here is the worst case rather than a case.

Everything is linear elastic. A local yield at the point of application redistributes the peak and shortens the decay, in the direction of safety, which is why the principle survives in structures whose materials do not obey any of the arithmetic above.

And the whole essay is about a plane problem. In three dimensions a self-equilibrating patch on a surface decays as the inverse cube of distance for a doublet and faster for higher-order systems, which is quicker than anything drawn here — so the plane result is the conservative one, which is a pleasant place for an approximation to land.

What the pictures cannot show

The curves above are drawn against a distance measured in depths, and the region they are most interesting in — the first quarter of a depth — is where the model is least reliable, because that is where the loading detail the model has abstracted away actually is.

Nor can they show the thing that makes the principle useful in practice, which is not that stresses die away but that nobody has to model the end. The value of Saint-Venant’s principle is a permission to stop: to draw a bearing as an arrow, compute the member, and know that the arrow’s dishonesty is confined to a region that will be detailed separately by somebody applying a different set of rules.

The assumption the figure rests on

The disturbance is assumed to be applied to a member that is otherwise uniform and continues far enough for it to decay in.

Both parts fail routinely. A beam with a bearing at one end and a splice a metre away has two disturbed regions overlapping, and their sum is not either of them. A short cantilever bracket is entirely disturbed region — there is nowhere in it that is a depth from both the load and the support — which is why brackets are designed by strut-and-tie rather than by section, and why the anchorage zone has an essay of its own.

The honest form of the principle is therefore conditional: stresses converge to the beam-theory answer at a distance of about one depth from a disturbance, provided there is a depth of member to converge in. A great many of the places where structures actually fail are places where there was not.

A truss cut through panel 2The truss severed through one panel, with everything to the right removed and the three cut member forces drawn on the exposed faces. Taking moments about the marked joint removes two of the three unknowns, so one equation gives the third: -56.47.30.0the cut — three members, three unknowns35.29-56.4727.79moments about here kill two of the threesolved without touching any of the other 18 members
Fig. 10 The cut this collection is built on. Everything a free body says about internal forces is exact and says nothing about how they are distributed across the cut. Saint-Venant’s principle is the bridge between the two, and one depth is the length of the bridge.

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.

Airy stress functionBimomentDecay lengthDisturbed regionEigenvalueElasticityEnd effectFourier seriesFree bodyPlane sectionsSaint venants principleSelf equilibratingStress concentrationSuperpositionWarping