Equilibrium

Two ways of being wrong

Plastic analysis has two theorems and they point in opposite directions. Any equilibrium field that nowhere exceeds the plastic moment gives a load at or below the collapse load; any mechanism gives one at or above it. Only one of those errors is safe.

Assumes After the first yield, which is not the end, One support too many, and what it costs to know and Everything adds to nothing, and that is the whole of statics.

Almost every argument in this collection about a structure’s strength has quietly used one of two theorems, and they are worth separating, because they are not two ways of getting the same answer. They are two ways of getting a different answer, in opposite directions, and only one of the two errors is on the safe side.

The static, or lower-bound, theorem. If a distribution of internal forces can be found that is in equilibrium with the applied load and nowhere exceeds the strength of any member, the structure will carry that load. It does not matter whether the distribution is the real one.

The kinematic, or upper-bound, theorem. If a mechanism can be found, the load computed from its work equation is at or above the collapse load. It does not matter whether it is the mechanism that actually forms.

The two theorems close on the answer from opposite sides. A pinned-base portal frame under 200 kN at mid-span and 80 kN at the eaves, with a plastic moment of 200 kNm. The three horizontal lines are mechanisms — work equations, each an upper bound: the beam mechanism gives 1.000, sway 1.250 and the combined one 0.714. The curve is the lower bound: for each value of the one redundant, the largest load factor whose moment field stays inside the plastic moment everywhere. Its peak is 0.714, and it touches the lowest mechanism exactly — so the collapse load is known rather than bracketed. Using the beam mechanism instead would have claimed 1.40 times the real capacity, and the sway one 1.75 times, both of them on the wrong side.
Fig. 1 A pinned-base portal frame under 200 kN at mid-span and 80 at the eaves, with a plastic moment of 200 kNm. Three horizontal lines are mechanisms, each an upper bound; the curve is the best load factor each value of the one redundant admits, and its peak is the lower bound. They meet at 0.714.

Which free body produced the number

The frame has two pinned bases, so it is once redundant. Release the horizontal reaction at the right base and the primary structure is determinate: the moments at the two eaves and at mid-span follow from statics, and the redundant adds a moment of XhX h to each of them with a sign that depends on which side of the frame the section is.

Every value of XX therefore gives a complete equilibrium field for a given load factor. The lower-bound theorem says the frame will carry any load factor for which some XX keeps all three moments inside MpM_p — so the best lower bound is found by asking, for each XX, how large the load factor can be, and taking the largest.

That is a one-parameter search with a closed answer: each section imposes an interval on XX, and the field exists while the three intervals overlap. The largest load factor for which they do is 0.71429.

The three mechanisms, and the two that are wrong

The upper bounds come from work equations, and each is a line of arithmetic.

Beam mechanism — hinges at the two eaves and mid-span. Internal work 4Mpθ4M_p\theta, external λWLθ/2\lambda W L\theta/2, so λ=8Mp/WL=1.000\lambda = 8M_p/WL = 1.000.

Sway mechanism — hinges at the two eaves only, the bases being pinned. Internal 2Mpθ2M_p\theta, external λHhθ\lambda H h\theta, so λ=2Mp/Hh=1.250\lambda = 2M_p/Hh = 1.250.

Combined — the two added, with the hinge at one eaves cancelling. Internal 4Mpθ4M_p\theta, external λ(WL/2+Hh)θ\lambda(WL/2 + Hh)\theta, so λ=4Mp/(WL/2+Hh)=0.714\lambda = 4M_p/(WL/2 + Hh) = 0.714.

The lowest of the three is the answer, and the other two are not approximations to it — they are claims that the frame is 40 and 75 per cent stronger than it is, produced by work equations with no error in them.

That is the whole danger of the kinematic method. It is easy, it is quick, its arithmetic is elementary, and every mistake it makes is unsafe. A mechanism that has been overlooked does not produce a warning; it produces a higher number.

The minimisation is not confined to frames with several mechanisms to choose between. It is there on the simplest structure the theorem applies to, which is a single member yielding one section at a time.

The collapse mechanism of a propped cantilever. A collapse mechanism, with the hinge position found by searching rather than quoted. Every position gives an upper bound on the collapse load; the lowest is 36.43, at a hinge 58.6 per cent along, which is a coefficient of 11.657 times Mp over the square of the span.
Fig. 2 The same procedure on one member. A propped cantilever of 8 m with a plastic moment of 200 kNm has its span hinge at a position nobody can quote: every position gives an admissible mechanism and an upper bound, and the lowest of them is 7.29 kN/m, at a hinge 58.6 per cent along. That is a coefficient of 11.657 times Mp over the square of the span, and the two neighbouring positions give larger loads that are just as arithmetically correct.

Why the two meet, and what it means when they do

At λ=0.714\lambda = 0.714 the equilibrium field has ∣M∣=Mp|M| = M_p at mid-span and at one eaves, and 0.143 MpM_p at the other. Those two sections are exactly the two the combined mechanism puts hinges at.

That is not a coincidence; it is the reason the theorems are worth having. When a load factor is both an upper bound and a lower bound it is the collapse load exactly, and the collapse load has then been proved rather than computed — by two arguments that share no arithmetic and no assumptions beyond ductility.

It also says what a plastic analysis is for. There is no differential equation to solve, no stiffness to know, no sequence of yielding to follow. The answer is settled by a statement about equilibrium and a statement about work, and neither of them mentions how the structure got there.

Which way each method is wrong, and which one designers use

Put simply: the lower bound is the design method and the upper bound is the check.

A designer choosing reinforcement, sizing a bolt group, or laying out a strut-and-tie model is constructing an equilibrium field and satisfying it everywhere — which is a lower-bound calculation whose answer is safe whether or not it is the real distribution. That is why this collection has so many of them.

A thrust line inside the masonry is a lower bound: any line that stays inside proves the arch stands. A strut-and-tie model is a lower bound: any admissible truss proves the region works. Redistributing moments in a continuous beam is a lower bound: any moment diagram in equilibrium with the load may be designed for. Three quite different design methods, all of them the same theorem.

A line of thrust, and the masonry it has to stay inside. An arch ring of 10% of the span in thickness, rising 30% of the span, under its own weight as a uniform load. Any horizontal thrust between 4.29 and 5.85 puts a line of compression entirely inside the masonry, so the arch stands — and which of them it actually takes is not decided by statics. The two extremes are drawn: the minimum-thrust line, which rides high at the crown and low at the haunches, and the maximum-thrust line, which does the opposite.
Fig. 3 The oldest lower-bound calculation there is. A thrust line that fits inside the stonework is an equilibrium field that nowhere exceeds the material’s capacity, so the arch carries the load — and the theorem is what licenses the conclusion, since nobody claims the real line is the one drawn.

The upper bound is what an assessor uses. Given a structure that already exists, the fastest way to find out what it will carry is to look for the mechanism it will form — and the discipline that goes with the method is to look for all of them, because the one that was missed is the one that governs.

The price of admission is ductility

Both theorems assume the same thing, and it is not stated in either.

They assume a section that reaches its capacity and holds it while the rest of the structure catches up. The lower-bound theorem needs it because the equilibrium field it constructs is not the elastic one — getting from the elastic distribution to the assumed one requires some sections to yield and go on yielding while others load up. The upper-bound theorem needs it because a mechanism requires every hinge to be turning at once at full moment.

A structure that cannot do that is entitled to neither theorem. Three cases matter:

A brittle failure mode anywhere voids both. A shear failure, a bolt in tension, a weld, an anchorage — each reaches a capacity and then loses it, so the redistribution the theorems assume never happens. That is why plastic design of steelwork carries a requirement that shear and connections be stronger than the members they join: not for their own sake, but to keep the mechanism in the ductile parts.

A slender section voids both. A section that buckles locally before it reaches its plastic moment cannot form a hinge that holds, which is what the whole classification system exists to identify.

And a limited rotation capacity limits how far the redistribution may be taken, which is why codes cap moment redistribution at a percentage rather than allowing any equilibrium field at all.

Where the class limits come from. The width-to-thickness ratio at which two kinds of plate reaches its own elastic critical stress at the yield stress, for three steel grades. A flange outstand (buckling coefficient 0.43) derives to 17.2, 15.2, 13.3 at 275, 355, 460 N/mm², against quoted limits of 12.9, 11.4, 10.0; A web, in bending (buckling coefficient 4) derives to 52.5, 46.2, 40.6 at 275, 355, 460 N/mm², against quoted limits of 38.8, 34.2, 30.0. The derived number is the larger every time, and by the same factor at every grade — flange outstand 1.33, web, in bending 1.35 — because both the derivation and the quoted limit go as one over the root of the yield stress. A constant ratio is what a fixed knockdown looks like: the derivation is for a perfect plate and the quoted limit is for a rolled one, carrying residual stress and not quite flat.
Fig. 4 The classification that decides whether a section is allowed into the argument. A hinge that sheds moment as it rotates is not a hinge in the sense either theorem needs, and a section is admitted or not on that basis alone.

Counting the mechanisms, which is where the method actually fails

The kinematic method’s difficulty is not arithmetic. It is completeness.

A frame with nn critical sections and rr redundancies has n−rn - r independent mechanisms, and every other mechanism is a combination of those. The portal here has three critical sections and one redundancy, so there are two independent mechanisms — beam and sway — and the third is their sum.

That count is the discipline the method needs. Enumerate the independent mechanisms, combine them in every way that cancels a hinge, and minimise over the lot. Miss one and the answer is too high, with nothing to say so.

And the combinations are where the answer usually is. On this frame the two independent mechanisms give 1.00 and 1.25 and their combination gives 0.714 — the combination is 29 per cent below the lower of the two it is made of, because cancelling a hinge removes internal work while the external work of both is retained. A designer who checks the obvious mechanisms and stops has checked the two that are not the answer.

The count itself is the same subtraction that decides whether statics can answer a structure at all, read once for an analysis and once for a collapse. What it cannot say is which of the mechanisms it counts will govern, and that is decided by the load rather than by the frame.

The two theorems close on the answer from opposite sides. A pinned-base portal frame under 60 kN at mid-span and 80 kN at the eaves, with a plastic moment of 200 kNm. The three horizontal lines are mechanisms — work equations, each an upper bound: the beam mechanism gives 3.333, sway 1.250 and the combined one 1.429. The curve is the lower bound: for each value of the one redundant, the largest load factor whose moment field stays inside the plastic moment everywhere. Its peak is 1.250, and it touches the lowest mechanism exactly — so the collapse load is known rather than bracketed. Using the beam mechanism instead would have claimed 2.67 times the real capacity, and the sway one 2.67 times, both of them on the wrong side.
Fig. 5 The identical frame with its vertical load cut from 200 kN to 60 and everything else untouched. The order has changed completely: the beam mechanism is now 3.333, the sway 1.250, and the combination 1.429 — so the combination is no longer the answer and the bracket closes on the sway mechanism at 1.250. Either independent mechanism claims 2.67 times the real capacity here, and which of the three governs is a property of the load ratio, not of the structure.

Masonry is the interesting case

An arch of dry stone has no tensile strength at all, no ductility in the usual sense, and the lower-bound theorem applies to it perfectly.

The reason is that its limit state is not a material strength. The stones do not crush — the stresses in a masonry arch are a small fraction of the stone’s strength — so what limits the structure is geometry: whether a line of thrust in equilibrium with the load can be found that stays inside the masonry.

That is a lower-bound statement with a purely geometric admissibility condition, and it is exactly Heyman’s insight: masonry needs no ductility because it is never asked for any. Its collapse is a mechanism of rigid blocks rotating about hinges at the surfaces, which makes the upper-bound theorem available too, and the two bracket the answer in the ordinary way.

Which is a good place to see what the theorems really need. Not ductility as such, but the ability to reach an assumed distribution without anything failing on the way. Stone gets there by having enormous reserve in the only mode it can fail in; steel gets there by yielding; and concrete gets there by having reinforcement that yields, which is the same argument with the ductility supplied by an inclusion.

What the theorems are buying in every one of those cases is a moment held constant over a rotation, and the length of that plateau is what decides how much redistribution a real structure will deliver. It is the same property a semi-rigid joint has some of and a nominally pinned one has none of, measured on the member rather than on the connection.

What a plastic analysis buys, and what it does not

It is worth being precise about the size of the prize, because the theorems are often presented as though they were mainly an economy.

For a simply supported beam, plastic analysis buys the shape factor and nothing else — 13 per cent for an I-section, 50 for a rectangle. There is one critical section, one hinge turns it into a mechanism, and no redistribution is available.

For a continuous or framed structure it buys the shape factor and the redistribution, and the second is usually the larger. A fixed-ended beam under a uniform load has an elastic moment ratio of 2:1 between support and span; at collapse the two are equal, and the collapse load is 16Mp/L216M_p/L^2 against the 12Mp/L212M_p/L^2 at which the first hinge forms — a third more, before the shape factor is counted at all.

Which locates the economy exactly. It comes from redundancy, and a structure with none gets none of it. That is the same trade one support too many is about, priced in strength rather than in the effort of an analysis.

The two cases are worth drawing side by side, because the whole of the economy is the difference between them.

The collapse mechanism of a simple span. A collapse mechanism, with the hinge position found by searching rather than quoted. Every position gives an upper bound on the collapse load; the lowest is 25.00, at a hinge 50.0 per cent along, which is a coefficient of 8.000 times Mp over the square of the span.
Fig. 6 The determinate case, with nothing to redistribute. An 8 m simple span of 200 kNm plastic moment collapses at 25.00 kN/m, at a hinge exactly half way along — a coefficient of 8.000 times Mp over the square of the span, which is the elastic collapse of the same beam with the shape factor and nothing else added to it.

Now fix both ends of the same member and search again. Nothing about the section, the span or the material has changed; only the number of places a hinge has to form before the beam is a mechanism.

The collapse mechanism of a fixed-ended beam. A collapse mechanism, with the hinge position found by searching rather than quoted. Every position gives an upper bound on the collapse load; the lowest is 50.00, at a hinge 50.0 per cent along, which is a coefficient of 16.000 times Mp over the square of the span.
Fig. 7 The same beam, the same plastic moment, built in at both ends. The collapse load is 50.00 kN/m at a coefficient of 16.000 — exactly twice the simple span’s 8.000, because three hinges have to form instead of one. The first of those hinges arrives at 12 times Mp over the span squared, so a third of the capacity drawn here is redistribution that statics alone would never have found.

Why the theorems arrived so late

Both theorems are twentieth-century results, which is startling for statements this simple about structures people had been building for millennia.

The reason is that they are about materials rather than about statics. A lower-bound theorem is worthless unless a section can hold its capacity while its neighbours load up, and until steel was made ductile and reinforced concrete was made to yield in its steel rather than crush in its concrete, no structural material reliably could. Masonry is the exception that proves it: its limit state is geometric rather than material, so the theorem was available all along and was used by anyone who drew a thrust line, several centuries before anybody proved it.

That also explains the order in which the two were adopted. The kinematic method came first, in the 1930s and 1940s, because a work equation is a calculation a designer can do; the static method’s dominance came later, with the strut-and-tie and yield-line methods, once it was clear that the safe side of the answer was worth the extra construction.

The point the rafter turns about, which is off the frame. A pitched portal of 8 m span and 4.0 m to the eaves, with a 1.0 m rise, collapsing. Each rigid part of the mechanism rotates about some point: the left column about its base hinge, the right about its own. The rafter between them does neither, and its centre is found by one rule — two bodies joined at a hinge share that hinge, so the second body's centre lies on the line through the first body's centre and the hinge, extended. Two hinges give two lines and they cross at (8.0, 10.0) metres, which is 5.0 m above the ridge and outside any drawing of the frame itself. From there the whole collapse is two ratios of lengths and no trigonometry: the load factor is 0.804. Flatten the roof and the centre descends; make the two lines parallel and it goes to infinity, which is the statement that the rafter translates instead of turning.
Fig. 8 Why the work equation was the practical method: it is a drawing. The same frame given a metre of roof pitch has its rafter turning about a point at (8.0, 10.0) metres, five metres above the ridge, and from that one intersection the collapse load reads off as two ratios of lengths — a load factor of 0.804 against the flat frame’s 0.714. No equilibrium field is constructed anywhere, and the point the mechanism turns about is found with a straight edge.

Any admissible field is safe, and they are not all equally good

The static theorem hands a designer a freedom that looks total. Any equilibrium field within capacity proves the structure carries the load, so pick one — the elastic distribution, a uniform one, a convenient one, whatever makes the reinforcement easy to draw — and the theorem certifies the result. It is worth being clear about what that licence does not cover, because the freedom is real and unlimited use of it produces bad structures that are nonetheless safe.

The cost of choosing a field far from the elastic one is paid entirely in how much the structure has to move to get there. A frame loaded from zero follows its elastic distribution first. Reaching the assumed field means yielding the sections the assumption overloaded and holding them there while the others catch up, and the further the assumed field is from the elastic one, the more rotation the first hinges are asked for. The theorem says the structure will arrive; it says nothing about the state it is in when it does.

Three consequences follow, and none of them is a strength check.

Rotation demand. A section is only entitled to the theorem while it can still deliver the rotation being asked. The percentage caps codes put on redistribution are a blunt way of bounding a quantity — required hinge rotation — that the theorem itself never mentions.

Cracking at service. The assumed field is a collapse state, and the structure spends its life near its elastic one. A beam designed for a support moment 40 per cent below the elastic value has reinforcement sized for a moment the support genuinely sees in service, and it will crack accordingly.

Deflection. Every early hinge is a loss of stiffness, and a structure that has redistributed heavily is softer than the analysis that certified it assumed.

So the working rule is the opposite of what the theorem’s generality suggests: stay near the elastic field unless there is a reason to leave it, and treat every departure as a withdrawal against a ductility account. The theorem’s freedom is best used to move a moment out of an awkward place — a support congested with reinforcement, a joint that cannot be detailed — rather than to chase an economy, because the economy is collected at collapse and the price is paid every day the structure is in use.

Where the model stops

Everything above is first-order. A frame that sways significantly under load has second-order moments the mechanism’s work equation does not contain, and the whole analysis becomes an upper bound on something that is itself falling. Plastic design of sway frames therefore carries an explicit limit on slenderness, and beyond it the collapse load is a stability question rather than a strength one.

The plastic moment is treated as a constant. It is reduced by axial force and by shear, and in a frame’s columns the axial force is not small — so a real analysis iterates: find the collapse mechanism, compute the axial forces, reduce MpM_p, and repeat.

Only three critical sections were considered. That is right for a prismatic frame with point loads, and a member under a distributed load has its hinge at a position that has to be searched for, which is what the single-span figure above is doing.

And nothing here says how much deflection the collapse mechanism needs. A load factor of 0.714 is reached at a deflection that may be far outside anything serviceable, and a structure can be plastically adequate and unusable — which is the standing reason strength and stiffness are separate checks.

What the pictures cannot show

The bracket in the first figure closes to a point, which makes the answer look inevitable. It closed because a one-redundant frame with three critical sections has just enough freedom for the two bounds to meet; with more redundants and more sections the search is a linear program, and a hand calculation that stops early leaves a genuine gap between the two bounds rather than a point.

Nor can any figure show what a collapse mechanism looks like at the moment it forms, because it does not form at a moment. Hinges appear one at a time as the load rises, the structure softens each time one does, and the “mechanism” is the arrangement present at the end of a sequence that took the whole loading history to produce.

The assumption the figure rests on

The load is assumed to be applied proportionally — every component growing together, in a fixed ratio, characterised by a single factor λ\lambda.

Real loads do not arrive that way, and the consequence is more interesting than it sounds. Under varying repeated loads a frame can fail at a lower factor than either theorem gives, by ratcheting or alternating plasticity — accumulating deflection cycle by cycle at a load that would be perfectly safe if applied once and held.

So the collapse load computed here is a bound for one loading path, and the shakedown load, which is the honest limit for a structure under variable loading, is a different and smaller number obtained from a different pair of theorems. The two pairs are close relatives — one is a lower bound on a stress field and the other is a lower bound on a residual stress field — and knowing which one a design is relying on is the difference between a structure that stands and one that walks.

What this makes readable

Essays that name this one as a prerequisite.

Named alongside this one

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

What links here

The 8 essays that link to this one and share the most of its objects, of 21 that link here.

The objects this essay names

Each one links to every other essay that touches it.

Collapse mechanismDuctilityEquilibriumLimit analysisLoad factorLower-bound theoremMoment redistributionPlastic hingePlastic momentRedundancySafetyStrut-and-tieThrust lineUpper bound theoremVirtual work