Sections and stress

What is left after the first fibre yields

The elastic section modulus stops at the moment the outermost fibre reaches yield. Nothing else in the section has, so it goes on taking load — and how much more it takes turns out to be a property of the shape alone, with no dimension, no stress and no material anywhere in the answer.

Assumes Bending is a pair of forces, pushing and pulling, Plane sections stay plane, and what the assumption costs and After the first yield, which is not the end.

Elastic bending theory produces one number for a section and stops: the section modulus ZZ, and with it the moment My=fyZM_y = f_y Z at which the outermost fibre reaches yield. Everything on this site that computes a bending stress computes it against that.

The section, however, has not stopped. One fibre out of a few thousand has reached its limit, and the rest are carrying somewhere between nothing and fyf_y. Put more moment on and the yielded zone eats inwards from the top and the bottom, the elastic core between them shrinks, and the section keeps taking load until every fibre is at ±fy\pm f_y and there is nothing left to recruit. The moment then is the plastic moment Mp=fySM_p = f_y S, and the ratio

ν=SZ\nu = \frac{S}{Z}

is the shape factor. The remarkable thing about it is not that it exceeds one. It is what it contains.

What is left after the first fibre yields, which is a property of shapeThe shape factor — plastic modulus over elastic — for six sections, computed by finding each one's equal-area axis and summing ±f_y over it. The numbers contain no dimension, no stress and no material: a rectangle is exactly 3/2 whatever its size, a diamond exactly 2, a circle 16/3π. The spread is the argument. A i section keeps only 13 per cent in reserve past first yield, because nearly all its material is already at the extreme fibre and there is nothing further in to recruit; a diamond keeps 100 per cent, because most of its material is near the middle and doing very little elastically. So the section shapes that are best at elastic bending are the ones with the least left afterwards, which is exactly backwards from the way the reserve is usually described.1.129i section1.194box1.316tube→ 4/π1.500rectangle3/21.698circle16/3π2.000diamond211.21.41.61.822.2plastic modulus ÷ elastic modulusthe closed forms, where there is one, underneath
Fig. 1 The shape factor of six sections, each computed by finding its equal-area axis and summing ±f_y over it. No dimension appears in any of the numbers, and no material.

Two neutral axes, not one

The first surprise is not the ratio at all, it is the axis.

The elastic neutral axis is the centroid. That follows from the stress being proportional to distance: for the axial force on the section to vanish, ydA\int y \, dA must vanish about the axis, which is the definition of the centroid.

The plastic neutral axis is not. When every fibre carries ±fy\pm f_y and nothing in between, the axial force is fy(AtensionAcompression)f_y (A_{tension} - A_{compression}), and for that to vanish the two areas must be equal. The plastic axis is therefore the equal-area axis, and the first moment of area has nothing to do with it.

For a section symmetric about the axis of bending the two coincide, which is why the distinction almost never comes up: rectangles, circles, I-sections about their major axis, boxes. For anything else they do not. A tee of 200 mm flange and 400 mm depth has its centroid well down in the web and its equal-area axis up inside the flange — sixty-nine millimetres apart, a sixth of the depth.

A tee at 85% of its plastic momentThe same tee drawn three ways: the shape, the strain across its depth, and the stress that strain produces in mild steel. The strain diagram is a straight line, because plane sections stay plane whatever the material is doing. The stress diagram is not: 22% of the area has yielded, working inward from both faces, and the neutral axis sits at 158.8 mm against a centroid at 149.0 mm. The compression resultant is 243.7 kN and the tension resultant 243.7 kN, on a lever arm of 128.8 mm, which multiplies back to the 31.4 kNm the section is carrying.neutral axisteestrainalways a straight linestressthe material's own curve, sidewaysC = 243.7 kN · T = 243.7 kN · lever arm 129 mm · M = 31.4 kNm22% of the area has yielded — 0 mm from the top, 91 mm from the bottom · Mp = 36.9 kNm · shape factor 1.78
Fig. 2 The stress distribution in a partly yielded section. The elastic core is a wedge between two blocks of constant stress, and the axis between them is not where it was at first yield.

So a yielding asymmetric section migrates its own neutral axis, continuously, from the centroid at first yield to the equal-area axis at full plasticity. No elastic calculation contains that motion, and neither does the plastic one — each of them describes an end of it. What happens in between is a family of stress distributions with a shrinking elastic wedge that is not centred on either axis.

Four closed forms with nothing in them

For simple shapes the ratio comes out exactly, and it is worth writing the numbers down because they are so bare.

A rectangle of breadth bb and depth dd: Z=bd2/6Z = bd^2/6 and S=bd2/4S = bd^2/4, since the plastic stress block is two rectangles of area bd/2bd/2 with their centroids d/4d/4 either side. So ν=3/2\nu = 3/2, exactly, whatever bb and dd are.

A circle of diameter DD: Z=πD3/32Z = \pi D^3/32 and S=D3/6S = D^3/6, giving ν=16/3π=1.6977\nu = 16/3\pi = 1.6977.

A diamond — a rhombus on its point: Z=bd2/24Z = bd^2/24 and S=bd2/12S = bd^2/12, so ν=2\nu = 2 exactly.

A thin circular tube: ν4/π=1.2732\nu \to 4/\pi = 1.2732 as the wall gets thin.

Not one of those contains a length, a stress, a modulus or a material. That is the sense in which the shape factor is a property of shape: it is a pure number attached to a form, in the way that a rectangle’s ratio of area to bounding box is, and it survives every scaling of the section that keeps the form.

The ordering, which is backwards

Read the six values on the figure above and something uncomfortable falls out. The sections that are best at elastic bending have the least left afterwards.

A rolled I-section — the shape this whole subject exists to justify, the one that gets the material as far from the middle as it can — has a shape factor of about 1.13. A diamond, which is the worst conceivable arrangement of a given area for bending, has 2.00. A rectangle, halfway between, has 1.50.

The reason is exactly the reason the I-section is good. Its material is already at the extreme fibre, so when the extreme fibre yields there is very little further in that was under-stressed and can be recruited. The diamond’s material is nearly all near the middle, doing almost nothing elastically, and the plastic stress block puts every bit of it to work at full stress.

Efficiency and reserve are the same quantity read twice. A section that wastes nothing elastically has nothing to recover; a section that wastes most of itself has a great deal. That symmetry is worth carrying, because it means the phrase “plastic design is worth thirteen per cent” is a statement about I-sections specifically and not about plastic design.

Every strip counts by the square of its distanceA rectangular section divided into equal strips, with each strip's contribution to the second moment of area drawn beside it. The strips are identical in size; only their distance from the neutral axis differs.neutral axiscontribution of each striptotal I = 1066.67 × 10⁶the outer strips do almost all of the work
Fig. 3 Where the second moment of area comes from — the material far out, squared. The same argument read the other way says where the plastic reserve does not come from.

The reserve is approached and not reached

The second uncomfortable fact is that MpM_p is an asymptote.

For a rectangle the moment–curvature relation past first yield is exact and short. With the elastic core reaching to ±ye\pm y_e and the rest yielded,

MMy=3212(κyκ)2\frac{M}{M_y} = \frac{3}{2} - \frac{1}{2}\left(\frac{\kappa_y}{\kappa}\right)^2

which reaches 3/23/2 only as κ\kappa \to \infty. At three times the yield curvature the section is at 96.3 per cent of MpM_p; at ten times, 99.5.

What it costs to reach the plastic moment, for two shapesMoment against curvature for two cross-sections of identical area (3000 mm²) and identical depth (200 mm), in mild steel, each divided by its own first-yield moment and its own first-yield curvature. The rectangle has a shape factor of 1.50 and reaches 98% of its plastic moment at 4.3 times the curvature at first yield; The I-section has a shape factor of 1.09 and reaches 98% of its plastic moment at 1.2 times the curvature at first yield. The dashed lines are the rigid-plastic moments, computed from the equal-area axis rather than read off the curves, and no curve reaches its own.02468101200.511.5curvature ÷ curvature at first yieldmoment ÷ moment at first yieldrectangle: 1.50× the yield moment, at 4.3× the yield curvatureI-section: 1.09× the yield moment, at 1.2× the yield curvature
Fig. 4 Moment against curvature for two sections. The approach to the plastic moment is asymptotic in both, and the section with the larger reserve is the one that takes longer to collect it.

Two consequences follow, and the second is the one that decides real design.

The first is that the plastic moment is never actually developed. Every plastic hinge in every collapse mechanism on this site is carrying a few per cent less than MpM_p, and the calculation is conservative by that amount — which is a small comfort next to how large the modelling assumptions are.

The second is that the reserve has to be paid for in curvature, and the section has to survive supplying it. A stocky section can rotate a long way past first yield with its compression flange still flat; a slender one buckles locally somewhere on the way and stops. That is the entire content of section classification: a class 1 section can reach MpM_p and hold it long enough for a mechanism to form, a class 2 can reach it and not hold it, a class 3 can reach MyM_y only, and a class 4 cannot even reach that.

Where the class limits come fromThe width-to-thickness ratio at which two kinds of plate reaches its own elastic critical stress at the yield stress, for two steel grades. A flange outstand (buckling coefficient 0.43) derives to 18.6, 15.2 at 235, 355 N/mm², against quoted limits of 14.0, 11.4; A web, in bending (buckling coefficient 4) derives to 56.8, 46.2 at 235, 355 N/mm², against quoted limits of 42.0, 34.2. 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.flange outstandk = 0.4318.6 at 23515.2 at 355quoted: 14ε1.33× the quoted limit, at every gradeweb, in bendingk = 456.8 at 23546.2 at 355quoted: 42ε1.35× the quoted limit, at every grade0102030405060width ÷ thickness
Fig. 5 The classification limits, which are a statement about how much curvature a section can supply before its compression parts stop being flat. The shape factor says how much reserve there is; this says whether it can be collected.

So the shape factor and the classification are two halves of one question, and they pull in opposite directions on the same variable. Making a section more efficient — thinner flanges further out — lowers the shape factor and pushes the section towards a higher class. The reserve gets smaller and harder to reach at the same time.

Which free body produced the number

The plastic modulus deserves the same treatment every other quantity here gets: a free body, and an equation that could have come out wrong.

Cut the section on its axis of bending and take the piece above. Every fibre in it carries fyf_y in compression, so the total force on it is fyAcf_y A_c, and it acts through the centroid of that area. The piece below carries fyAtf_y A_t in tension through its own centroid. Axial equilibrium demands Ac=AtA_c = A_t, which fixes the axis; moment equilibrium about any point then gives

Mp=fy(Acyˉc+Atyˉt)=fyA2(yˉc+yˉt)M_p = f_y (A_c \bar{y}_c + A_t \bar{y}_t) = f_y \cdot \frac{A}{2}(\bar{y}_c + \bar{y}_t)

so SS is half the area times the distance between the two half-centroids, which is a lever arm in the ordinary sense. For a rectangle that is (bd/2)(d/2)=bd2/4(bd/2)(d/2) = bd^2/4 in one line, and it is the same construction a stress block in a concrete section uses.

The check that matters is that the same number can be reached from the other end. Sweeping the curvature of a fibre model to infinity and reading off the moment gives SS by a completely different route — an integration rather than a pair of centroids — and the two agree to the seventh figure. Neither derivation contains the other.

What the picture cannot show

Everything above is a section under pure bending, and three things are missing from it.

Axial force. The equal-area axis is where it is because the net force is zero. Put an axial load on the section and the axis moves off it, the two areas stop being equal, and the plastic moment falls — which is the whole of the N–M interaction diagram. For a section carrying half its squash load, the plastic moment can be down by a quarter.

Two ways to fail, and the curve between themThe exact plastic interaction between axial force and moment for one section of identical area, both normalised by their own squash load and their own plastic moment. The I-section stands 6.0% of its plastic moment outside the straight line at an axial ratio of 0.12. Every section here is symmetric about its centroid, so the equal-area axis and the centroid coincide and it makes no difference which the moments are taken about. The straight line is the rule that says the two capacities share out in proportion, and everything between it and a curve is capacity that rule gives away.00.20.40.60.8100.20.40.60.81moment ÷ plastic momentaxial force ÷ squash loadI-section: 6.0% of Mp outside the linethe straight-line rule
Fig. 6 The plastic moment as a function of the axial force on the section. The shape factor quoted anywhere is the value at zero axial load, which is a corner of this curve.

Shear. The web is carrying shear as well, and where the shear is a large fraction of the plastic shear capacity the material in the web is using part of its yield criterion up on it. The standard treatment is to reduce the web’s yield stress, which reduces SSand the reduction only bites above about half the shear capacity.

Which axis. Every number here is about one bending axis. A section bent about a diagonal has a different equal-area axis, a different plastic modulus and a different shape factor, and for an angle or a channel that axis is nowhere near either principal direction. The shape factor is not one number attached to a section — it is one number attached to a section and a direction.

The forty years it took to be allowed

The arithmetic in this essay is not hard and it was not new when it became legal. Every term in it was available to a competent engineer in 1900, and plastic design entered British practice in the 1950s.

Gábor Kazinczy tested fixed-ended beams in Budapest in 1914 and reported what the moment–curvature curve above predicts: they carried roughly twice the load at which the first section reached yield, and they did so by forming hinges in a sequence. Hermann Maier-Leibnitz did the same on continuous beams in the 1920s and made the further point that the collapse load did not depend on the support settlements — which is the statement that a plastic calculation does not need the compatibility a continuous elastic analysis spends most of its effort on.

What took the remaining decades was not evidence but the theorems. Until the upper- and lower-bound theorems were stated and proved, a plastic calculation was a plausible estimate with no way of telling which side of the answer it fell on, and no code can be written on that. Once the two bounds existed, a mechanism became a demonstrable overestimate and an equilibrium field a demonstrable underestimate, and the collapse load became a quantity that could be bracketed rather than guessed.

John Baker’s work at Cambridge supplied the third thing needed, which was a reason to care: wartime air-raid shelters designed elastically were heavy, and designed plastically they were not. The Morrison shelter is a plastic design, and the argument that got the method accepted was made out of steel that had to be saved rather than out of theorems.

Where plane sections stop staying planeStrain across a cut face at three 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 = 4plane sections holdspan ÷ depth = 8plane sections holdspan ÷ depth = 16plane sections holdthe assumption is the theory — everything else is arithmetic on top of it
Fig. 7 The assumption every line of this essay rests on, and the one that survives yielding unchanged: sections stay plane. What stops being true when the material yields is the proportionality between strain and stress, not the linearity of the strain.

Where it changes what gets built

The place plastic reserve is worth most is a continuous member, and it is not because of the shape factor.

A continuous beam gets its economy from redistributing moment: the support section yields first, forms a hinge, and passes further load to midspan until a second hinge turns the span into a mechanism. That is worth far more than 13 per cent — a uniformly loaded fixed-ended beam has a collapse load twice its first-yield load — and it needs the shape factor only to be non-zero.

The collapse mechanism of a fixed-ended beamA 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 10.00, at a hinge 50.0 per cent along, which is a coefficient of 16.000 times Mp over the square of the span.sagging hinge at 4.00hinge at the fixed endand herelowest upper bound: 10.00every hinge position gives an upper bound on the collapse loadassumed position of the sagging hingecoefficient 16.00 Mp ÷ L²
Fig. 8 The mechanism a continuous member collapses by. The shape factor is worth a few per cent of one hinge’s capacity; the redistribution between hinges is worth a great deal more.

There is a third place it matters and it is a detail rather than a member. A bolt group, a base plate and a beam-to-column endplate are all sections in the sense used here — a set of areas resisting a moment — and each of them has an elastic distribution and a plastic one with a ratio between them. The endplate’s is often the largest shape factor in the whole structure, because the plate’s material is spread evenly rather than concentrated at the extremes, and it is the reason a connection designed elastically and tested to failure so consistently outperforms its calculation.

A simply supported beam gets nothing but the shape factor, because there is only one hinge and it is the mechanism. For a rolled I-section that is 13 per cent, on a member whose deflection is usually the thing that governs anyway. That is why plastic design changed continuous frames and left simply supported floor beams almost exactly as they were.

Bending is a push and a pullA section carrying a bending moment, with the stress at every height computed as the moment times the distance from the neutral axis divided by the second moment of area. It is compression above and tension below, and zero exactly at the neutral axis.neutral axiscompressiontensionI = 29.97 × 10⁶Z = 299.7 × 10³peak stress 200.2σ = M y ÷ I, at every height
Fig. 9 The stress block partway between the two limits, which is the state a real section is actually in at the moment a hinge is said to have formed.

The generalisation

Two sentences are worth separating before the general one, because they are the practical residue of everything above.

The shape factor is a property of a form and a direction, and it is a pure number. It does not scale, it does not depend on the steel, and it can be written down for any section anybody can draw by finding one axis and taking two centroids. Where a section is symmetric the axis is the centroid; where it is not, the axis moves as the section yields and neither end of that motion is where an elastic calculation put it.

And the reserve is a trade rather than a gift. It is collected in curvature, the curvature has to be survivable, and the survivability is a separate calculation with its own limits — which is why the two most efficient shapes in the collection have the smallest reserves and the hardest time reaching them.

The habit this essay is really about is the difference between a limit set by a point and a limit set by a body.

ZZ is a point limit: it is decided entirely by the single most stressed fibre, and every other fibre in the section is irrelevant to it. SS is a body limit: every fibre appears, weighted by area and lever arm, and no single one of them decides anything. The same pair recurs everywhere on this site — a bolt group’s elastic and plastic distributions, a weld group’s, a frame’s first-yield load and its collapse load, and a soil bearing pressure’s peak against its resultant.

In every one of those pairs, the ratio between the two answers is a pure number attached to the geometry, and reaching the second requires the material to redistribute. The shape factor is the simplest member of the family and the only one where the number can be written down exactly for a shape a reader can draw.

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.

Bound theoremsDuctilityEqual area axisLever armMoment curvatureNeutral axisPlane sectionsPlastic hingePlastic modulusSecond momentSection classificationSection modulusShape factorStress blockYield stress