Mechanics

The one number the tolerances cannot touch

A redundant structure's load sharing depends on stiffnesses and manufacturing errors that nobody knows. Its collapse load does not depend on either. Once members yield they hold a known force instead of a force proportional to a displacement, the compatibility equations that needed the unknowns drop out, and the load at which the structure becomes a mechanism follows from a work balance with no stiffness in it at all.

Assumes: The table statics cannot settle · A state no load could reach

Three earlier arguments have been about a quantity nobody can compute. Four legs and three equations leave the load sharing undetermined; the stiffnesses that would settle it are known to a factor of two at best, and a leg that is stiff draws load to itself for no better reason; a manufacturing error of a tenth of a millimetre swamps them; and a temperature difference of one kelvin does more than the load.

Buildings nevertheless stand, and they are designed by people who are not troubled by any of this. The reason is that the quantity which decides whether a structure stands is a different quantity, and it does not depend on any of the unknowns.

Four tolerances, four elastic answers, one collapse load. The force in each of a table's four legs against the load on it, for four different manufacturing errors, with the legs made of a material that yields at 25 kilonewtons. While everything is elastic the short diagonal pair takes more than its share by a fixed amount that depends on the error and not at all on the load, so the load at which the first leg reaches its capacity runs from 25 to 100 kilonewtons — a spread of more than a factor of three. Past that point the yielded legs hold a constant force and the others take the rest, and the difference the tolerance made is erased. Every one of the four cases collapses at 100 kilonewtons, which is four times one leg's capacity, checked here to a part in a million across the four. The quantity nobody could compute and the quantity that decides whether the structure stands are not the same quantity.
Fig. 1 The force in each of a table’s four legs against the load on it, for four different manufacturing errors, with legs that yield at twenty-five kilonewtons. The load at which the first leg reaches its capacity runs from twenty-five kilonewtons to a hundred — a factor of four, set entirely by a tolerance. All four collapse at a hundred kilonewtons, which is four times one leg’s capacity, checked to a part in a million.

Why yielding erases the difficulty

The reason a redundant problem has no answer is a counting one: there are more unknown forces than equilibrium equations, and the missing equations have to come from compatibility — from the requirement that the members’ deformations fit together — which brings the stiffnesses in with them.

A member that has yielded does not obey that requirement. Once a leg is squashing plastically its force is its yield force, whatever its displacement does; the relation between force and deformation, which was the source of the extra equations and of all the trouble in them, has been replaced by a constant.

So each member that yields removes one unknown from the problem rather than adding one. Yield enough of them and the count comes back into balance: the structure becomes determinate, and then over-determinate, and at the moment it can no longer satisfy equilibrium at all it has become a mechanism and it collapses.

At that moment every member that is going to yield has yielded, so every force in the structure is a known capacity rather than an unknown share. The collapse load is a sum of capacities, and a sum of capacities has no stiffness in it, no tolerance, no temperature and no assembly sequence. It is the one quantity in the whole subject that the idealisation does not throw away.

The leg that is short, and the wobble that follows. The four reactions against how much shorter one leg is than the other three, with the load held fixed and the legs treated as springs of the same stiffness. At the shortfall of 0.02 drawn on the right-hand end of the chart the redistribution is -0.125, 0.125, -0.125, 0.125 of the load, and it is the same whatever the load is doing — computed at two quite different load positions and agreeing to 1.4e-16, which is what makes it a property of the table rather than of what is on it. The short leg's own reaction falls and the diagonally opposite one falls with it, while the other two rise; the pattern is exactly the null direction the rigid problem could not fix, which is the point. Past a shortfall of about 0.024 in these units a reaction goes negative, and since a leg cannot pull, the top lifts off and rocks between two three-legged states. That is a wobble, and it is why a three-legged stool has none: with three supports there is no null direction for a manufacturing error to excite. A stiffer leg makes the same shortfall worse in proportion, which is why a heavy rigid table on a stone floor rocks and a light one on a carpet does not.
Fig. 2 The elastic answer, for comparison: the four reactions against how much shorter one leg is, at a fixed load. The redistribution runs along the direction statics could not fix, is proportional to the shortfall, and is independent of the load entirely — so the gap between the two pairs at the start of the hero figure is this, and nothing in a drawing or a specification pins it down.

The contrast the two figures make is the whole point. The quantity that is unknowable and the quantity that decides safety are not the same quantity, and a great deal of nineteenth-century anxiety about indeterminate structures was anxiety about the wrong one.

Following a beam through it

The table collapses when its legs run out, which is simple enough to hide the interesting stage. A beam shows the redistribution happening.

The moment that stops growing, and the one that catches up. The bending moment at the fixed end of a propped cantilever and at the section under its central load, both against the load, for a beam whose section yields fully at a moment of 100 kilonewton-metres. Elastically the fixed end carries the larger moment and reaches that limit first, at 88.9 kilonewtons. It then stops growing — the section is yielded through and can carry no more — while the beam goes on taking load as though it were simply supported with a constant moment applied at one end, and the moment under the load rises faster than before. The second hinge forms at 100.0 kilonewtons, twelve and a half per cent above the first, and the beam becomes a mechanism. That collapse load is checked here against a virtual-work calculation on the mechanism alone, which knows nothing about the elastic stage and agrees to nine figures: the reserve depends on how the moments were distributed and the collapse load does not.
Fig. 3 The bending moment at the fixed end of a propped cantilever and at the section under its central load, both against the load. Elastically the fixed end carries the larger moment and reaches the section’s plastic capacity first, at eighty-nine kilonewtons. It then stops growing while the moment under the load rises faster than before, and the second hinge forms at a hundred — twelve and a half per cent later. The collapse load is checked against a work balance that knows nothing about the elastic stage.

Elastically a propped cantilever with a central load carries three-sixteenths of WLWL at its fixed end and five thirty-seconds under the load, a ratio of six to five. So the fixed end reaches the section’s plastic moment first, and a design that stopped there would call that the failure load.

What happens instead is that the fixed end yields through its whole depth and can carry no more moment — but it can still rotate, and while it does it holds the plastic moment constant. The beam from that point on behaves as a simply supported one with a fixed moment applied at one end, and every further increment of load goes into the section under the load. That section catches up, reaches the plastic moment in its turn, and with two hinges the beam has a mechanism and deflects without limit.

Twelve and a half per cent is the reserve, and its size is a measure of how unequal the elastic moments were. A beam whose elastic moments are nearly equal has almost no reserve; one whose moments are very unequal has a great deal. The reserve is therefore a property of the elastic distribution — which is to say of the stiffnesses — and the collapse load is not, which the figure checks by computing it a second way.

That second way is a virtual-work calculation on the mechanism alone. Give the midspan a small downward movement, work out how far each hinge rotates from geometry, multiply each rotation by the plastic moment to get the work dissipated, and set it equal to the work the load does. The answer agrees with the incremental calculation to nine figures, and it never mentioned a modulus, a second moment of area or an elastic moment distribution.

The one thing it costs

The redistribution above requires the fixed end to go on rotating after it has yielded, at constant moment, while the rest of the beam catches up. How far it has to rotate is calculable, and whether the section can rotate that far is a question about the material and the shape rather than about the structure.

A steel section can, comfortably: mild steel strains by twenty times its yield strain before it begins to harden, and by a hundred and fifty times before it breaks, so a compact steel section will rotate several times as far as any redistribution requires. That is the whole reason the method is a steel method.

A material that cannot does not get the reserve. Cast iron, unreinforced concrete in tension, glass, a high-strength bolt loaded into its own plastic range and a slender steel section that buckles locally before it yields fully all reach the first hinge and then break rather than rotate. For those, the first-yield load is the collapse load, and every unknown catalogued earlier comes straight back — which is why a prestressed concrete section is designed to stay uncracked rather than to redistribute, why the design rules for brittle materials are so much more conservative than the rules for ductile ones, and why they stay conservative no matter how well the loads are known.

Ductility is not a virtue of a material; it is what buys the indeterminacy back. A structure of a brittle material has the same free direction, the same unknowable load sharing and none of the relief, so its weakest member decides everything and its weakest member is decided by a tolerance.

Where else a mechanism decides the answer

The step the beam takes is general enough to be worth recognising elsewhere, because in each case a hard problem is replaced by a count of ways the thing can come apart.

A granular pack is rigid above the isostatic coordination number and flows below it, and the count of contacts against freedoms is exactly the count performed here — with the difference that a pack rearranges continuously and a frame does it in a handful of discrete steps. A soil slope fails along a surface found by searching over candidate slip circles for the least factor of safety, which is the third figure’s scan with a different parameter. And whether a pushed block slides or topples is the same question asked of a determinate body: two mechanisms, each with its own load, and the lower one happens.

The pattern is that strength is a property of the cheapest way out, and the cheapest way out is a geometric object rather than a stress. That is why the collapse calculation is indifferent to everything the elastic one needed, and why it has to be searched for rather than solved.

Two theorems, and only one of them is safe

The virtual-work calculation above needed a mechanism to be assumed. There is generally more than one, and which is assumed matters in a way that is worth being exact about.

Every guess is too high, so the answer is the lowest one. The collapse load a work balance returns for a propped cantilever under a uniform load, against where the sagging hinge is assumed to be, in units of the plastic moment over the span squared. Each point on the curve is a legitimate calculation: assume a mechanism, equate the work the load does to the work dissipated at the hinges, and solve. Every one of them is an overestimate, because a structure given a mechanism it does not want has to be pushed harder to use it. The least value is 11.657 times Mₚ/L², at a hinge 58.6 per cent of the way along — found here by scanning four hundred thousand positions and agreeing with the closed form to six figures. Guessing midspan, which is where a reader would put it, gives 12.00 and is 2.9 per cent unsafe if trusted.
Fig. 4 The collapse load a work balance returns for a propped cantilever under a uniform load, against where the sagging hinge is assumed to be. Every point is a legitimate calculation and every one is an overestimate. The least is 11.657 times the plastic moment over the span squared, at a hinge fifty-nine per cent of the way along, found by scanning four hundred thousand positions. Guessing midspan gives 12.00 and is nearly three per cent unsafe.

Assume a mechanism and equate work, and the load that comes out is the load at which the structure would collapse if it were forced to collapse that way. A structure given a mechanism it does not want has to be pushed harder to use it, so any assumed mechanism returns a load at least as large as the true one. That is the upper-bound theorem, and the true collapse load is the least value over all mechanisms.

The complementary statement runs the other way. Find any distribution of internal forces that balances the applied load and nowhere exceeds any member’s capacity, and the structure will carry that load. That is the lower-bound theorem, and it is safe: guessing badly gives a load that is too small.

The pair are what make plastic design usable, and the asymmetry between them is the thing to carry away. A mechanism is easy to imagine and gives an unsafe answer; an equilibrium stress field is harder to construct and gives a safe one. Practical design finds several mechanisms, takes the least, and then checks that a corresponding equilibrium field exists — at which point the upper and lower bounds coincide and the answer is exact.

The scan makes the shape of the risk visible. The curve is flat near its minimum, so a hinge position ten per cent wrong costs under one per cent in the answer, which is what makes the method workable by hand and why the classic textbook mechanisms give nearly the right answer. It is also what makes it treacherous: the error is always in the unsafe direction, and flat means there is no warning near the minimum to say whether it has been reached. The mechanism nobody thought of is the one that is not on the curve.

Where the method came from, and what it was resisted for

Plastic analysis was worked out in the 1930s and 1940s and met considerable opposition, which is instructive because the objection was reasonable.

Elastic design had an intelligible discipline: compute the stresses by resolving every force into components that suit the member, keep them below a fraction of the yield stress, and let nothing anywhere in the structure yield. Plastic design says that parts of the structure will yield in service, deliberately, and that this is fine. To an engineer trained on the first, the second sounds like a licence to accept damage.

What settled it was partly J. F. Baker’s wartime work on air-raid shelters, where the problem was to protect people with a minimum of steel and the elastic method’s answer was unaffordable. The shelters were designed by plastic analysis, they were permitted to deform substantially, and the criterion was that the occupants survive rather than that the steel stay elastic. That is exactly the right use of the method — it answers the question “at what load does this become a mechanism” and not the question “at what load does it stop looking new”.

The modern arrangement keeps both. Strength is checked against collapse, using this method; deflections, vibration and cracking are checked against the elastic behaviour at working load, using the other. Two calculations answering two different questions, and the reason there are two is precisely that the plastic one throws away the information the serviceability check needs.

What the collapse load will not tell anybody

It is worth being clear about how much is lost along with the stiffnesses, because the method’s own advantage is the source of its limits.

What the legs give is what decides. The four reactions again, this time with each leg treated as a spring rather than as a rigid point. The top is still rigid, so it settles into some plane, and a plane has three numbers in it — a height and two slopes. Three unknowns and three equations: the problem closes, and the answer is 0.256, 0.124, 0.426, 0.194 of the load. It lies on the family the rigid calculation produced, at the parameter 0.0850, to 9.7e-17 — so elasticity has not overruled statics, it has supplied the one piece of information statics was missing. The chart sweeps the load along a diagonal and shows the four reactions following it, with the near pair rising and the far pair falling. The stiffnesses used are 1, 3.5, 1, 0.3 relative to one another; making all four equal makes the answer the symmetric one, and the interesting cases are the unequal ones, where a stiff leg takes more than its geometric share simply for being stiff. The idealisation to unpack, and it is the one this whole figure is about, is that a real floor is not four independent springs either.
Fig. 5 What is thrown away: four legs of unequal stiffness under one rigid top, where a stiff leg draws load to itself simply for being stiff. Nothing about that arrangement survives into a collapse calculation, which is what makes the collapse calculation reliable — and means it can say nothing about how much anything moves, how the load is shared today, or which member will fatigue first.

It gives no deflections. At collapse the deflection is undefined; below collapse it depends on the elastic problem with all its unknowns. A structure can be perfectly safe against collapse and unusable because it sags, and the sag has to be computed by the method this one replaced — which needs every stiffness the collapse calculation was so pleased to be rid of, and inherits every uncertainty in them. Whether a contact grips at all is a related difficulty one level down: a quantity that has to be known for the serviceability answer and does not appear in the strength one.

It gives no fatigue life. Fatigue depends on the stress range at working load, which is an elastic quantity, and a member that carries a small share elastically and a full share at collapse is exactly the member a plastic calculation is indifferent to and a fatigue calculation is not.

It assumes the loads are applied once and monotonically. A structure loaded, unloaded and loaded again in a different pattern can accumulate plastic deformation cycle after cycle — incremental collapse, or ratcheting — at loads well below the single-application collapse load. Where a thermal load cycles, which is a self-loading structure’s subject, this is the governing case rather than an exotic one.

It assumes the capacity of a member is a number. A leg’s yield force, a section’s plastic moment: both are treated here as properties known in advance, and both are statistical quantities with their own spread in a real batch of material. The specified value is a lower percentile rather than a mean, which is one of the reasons measured capacities exceed calculated ones.

It has a third load hidden between its two. Under loads that are applied and removed repeatedly rather than once, there is a level — above first yield and below collapse — at which the structure yields during the first few cycles, arrives at a locked-in residual stress state, and from then on responds purely elastically for ever. That is the shakedown load, and the residual state it settles into is exactly the self-equilibrating state a temperature change produces, found by the structure itself rather than put there by a designer. Above the shakedown load and below the collapse load a structure either accumulates deformation every cycle or alternates plastically in one place until it cracks, and neither is visible in a monotonic calculation. Where the repeated load is a temperature, which is the commonest case, shakedown rather than collapse is what governs.

And it says nothing about stability. A mechanism forms at a computable load; a column buckles at a load computed by an entirely different argument, and for a slender member the second can be far lower. Plastic design applies to sections stocky enough to reach their plastic moment before anything buckles, and the section classification rules that decide which those are are most of the length of a steel code.

Three supports, and exactly one way to share the load. A rigid top on three supports, seen from above, with a load applied off centre. Three unknown reactions and three equations — the vertical forces balance, and the moments balance about two horizontal axes — so there is one answer and it is drawn: 0.430, 0.070, 0.500 of the load. Those three numbers are the barycentric coordinates of the load point in the triangle of support, which is why they sum to exactly one and why each is the area of the opposite sub-triangle divided by the whole. Nothing about the material enters: a top of steel and a top of cheese on legs of steel and legs of cheese give the same three numbers, provided both stay rigid enough to keep the geometry. Every reaction is positive, so the load lies inside the triangle of support and the top stands. Move it outside and one of these numbers goes negative, which is a leg pulling downwards — the tipping condition, and the only way this calculation can warn about it.
Fig. 6 And the determinate case, which has none of this on either side: three supports, three equations, one answer at every load. There is no redistribution available because there is nothing to redistribute to, so the first member to reach its capacity is the last — the reserve is exactly zero, and the reliability of the collapse calculation costs nothing because there was never anything unknown about the elastic one.

What perfect plasticity is standing in for

Yield is taken as perfectly plastic. Real steel strain-hardens, so a yielded section goes on gaining moment slowly rather than holding constant, which makes the true collapse load a little higher than computed — a conservative error, and the reason it is accepted without much comment.

The plastic moment is a section property computed for bending alone. A section also carrying axial force or shear reaches its bending capacity at a lower moment, and the interaction is nonlinear. In a frame where the columns carry substantial axial load, using the pure bending capacity is not conservative.

The mechanism is assumed to form at fixed geometry. A structure deflecting substantially before it collapses carries its load on a deformed shape, and the extra moments that produces can bring the collapse load down by ten or twenty per cent in a tall frame. That effect is second-order in the deflection and first-order in the consequences.

And the four-legged table’s collapse load is four times a leg’s capacity only for a central load. Move the load off centre and the collapse mechanism is a rotation about a line through two legs, the work balance is a different one, and the answer is lower — which is the upper-bound scan of the third figure applied to a different structure, and is the version of the tipping condition that a yielding structure has.

The settlement that tells the three stages apart

The hero figure draws four leg forces and cannot show the displacement, which is the thing that distinguishes the stages. Before first yield the top settles by micrometres; between first yield and collapse it settles by millimetres as the yielded legs squash; at collapse it moves without limit. A plot of force against load makes the second stage look like a mild change of slope and it is a change of kind. The same omission hides the warning a real structure gives: a steel frame approaching collapse sags visibly and audibly for some time first, which is a large part of why ductile design is trusted, and none of it is in the figures.

The moment figure draws two moments and hides the rest of the beam. A bending moment diagram is a function along the span, and what happens between first hinge and collapse is that its whole shape changes — the point of maximum sagging moment moves, the point of contraflexure moves, and the diagram at collapse is not a scaled version of the elastic one. Two numbers from it are the two that matter and they are not the picture.

And the mechanism scan draws one family of mechanisms, parameterised by one number. The real space of mechanisms for a frame is not a line: it is a set of combinations of beam mechanisms, sway mechanisms and joint mechanisms, combined in ways that cancel hinges, and the number of them grows quickly with the number of members. Drawing the search as a curve with a visible minimum is honest for this beam and flattering to the method in general.

Still open: how much of a real structure’s reserve is this reserve

The reserve computed here — twelve and a half per cent for the propped cantilever, a factor of up to four for the table — is the redistribution reserve, and it is the part of a structure’s margin that can be calculated.

It is not the whole margin and is usually not most of it. A real structure also has strain hardening, a yield stress typically ten to twenty per cent above the specified minimum, member sizes rounded up to what is manufactured, composite action with floor slabs that nobody counted, and restraint from cladding and partitions that no model contains. Load tests on old buildings routinely find capacities two or three times the calculated one, and the excess is mostly in that list rather than in the redistribution.

Which means the honest description of a structural safety factor is that it covers a set of uncertainties of very different characters: a variability in the loads that is statistical and estimable, a variability in the materials that is controlled and small, and a set of modelling omissions that are systematically in the safe direction and are not quantified at all. Attempts to make design fully probabilistic founder on the third, because a margin nobody has measured cannot be given a distribution.

The habit worth carrying away is the one the first two figures make side by side. When a quantity cannot be computed, ask whether it is the quantity the decision depends on. The elastic load sharing in a redundant structure genuinely cannot be known, and three earlier arguments are about why. The thing anybody actually wants to know — at what load does it fall down — turns out not to need it, and the reason is that the mechanism which destroys the structure is also the mechanism that destroys the ignorance.

Part 6 of 6

This essay is one argument about Free-body. The others:

The objects named here

The third axis, after the field and the reading path: the things themselves, and every essay that touches each one.

Constraint countingDuctilityPlasticityReaction forceSafety factorStatic equilibriumStatically indeterminateStiffnessVirtual-workYield