Sections and stress

Deliberately the wrong shape

Concrete in compression follows a curve, and no design office has ever integrated it. Every code in the world replaces it with a rectangle of reduced depth and reduced intensity, and the answer is right to a fraction of a per cent — not because the shapes are similar, which they visibly are not, but because a bending calculation only ever asks a stress distribution two questions.

Assumes Plane sections stay plane, and what the assumption costs, Bending is a pair of forces, pushing and pulling and The section that yields from the outside in.

Plane sections fixes the strain across a bending section: linear, zero at the neutral axis, largest at the face. The material then decides what stress goes with each strain, and for concrete that relationship is a curve — rising to a peak at about 0.2% strain, flat or falling from there, crushing at 0.35%.

Read that curve off the linear strain profile and the compression zone of a beam carries a stress distribution shaped like this:

Wrong in shape, right in two integrals. The compression zone of a C30 section with its neutral axis 150 mm down, drawn twice. The curved outline is the real parabolic-rectangular stress distribution — the material's own law read off the linear strain profile plane sections supplies. The rectangle over it is what every design office uses instead: intensity η f_cd = 16.5 MPa over a depth λx = 125 mm. The two shapes are visibly different and give the same answer, because a bending calculation asks a stress distribution only two questions — how much compression there is, and where its resultant acts. Both are 619 kN at 62.4 mm from the face. The factors are α = 0.8095 and β = 0.4160, and λ = 2β follows from wanting the same centroid. A triangle and a full rectangle match neither integral and are nowhere near.
Fig. 1 The compression zone of a C30 section with its neutral axis 150 mm down, drawn twice. The curved outline is the real stress distribution; the rectangle over it is what every design office uses instead. The two are visibly different and give the same answer.

Nobody integrates it. Every code in the world replaces the curve with a rectangle — reduced intensity, reduced depth — and gets the same answer. The interesting question is not whether that works, since it has worked for a century, but why it works, because the reason is not that the shapes are similar. They are not.

The two questions, which is the whole argument

A bending calculation interrogates a stress distribution exactly twice:

C=∫σ dAyˉ=∫σy dA∫σ dAC = \int \sigma\, dA \qquad\qquad \bar{y} = \frac{\int \sigma y \, dA}{\int \sigma \, dA}

How much compression is there, and where does its resultant act. Nothing else about the distribution enters the moment. The lever arm is d−yˉd - \bar{y}, the moment is C(d−yˉ)C(d - \bar{y}), and two distributions agreeing on those two numbers give identical moments whatever else they differ in.

So the equivalent rectangle is not chosen to look like the curve. It is chosen to have the same two integrals, and once that is done the resemblance is beside the point.

The construction takes two lines. Write the curve’s resultant as a fraction of the box it sits in and its centroid as a fraction of the depth:

α=Cfcd b x,β=yˉx\alpha = \frac{C}{f_{cd}\, b\, x}, \qquad \beta = \frac{\bar{y}}{x}

For the standard parabolic-rectangular curve taken to 0.35% strain these are

α=0.8095,β=0.4160\alpha = 0.8095, \qquad \beta = 0.4160

A rectangle of intensity ηfcd\eta f_{cd} over a depth λx\lambda x has resultant ηλfcdbx\eta \lambda f_{cd} b x and centroid λx/2\lambda x / 2. Matching both gives

λ=2β=0.832,η=α/λ=0.973\lambda = 2\beta = 0.832, \qquad \eta = \alpha/\lambda = 0.973

and that is the derivation in its entirety. λ=2β\lambda = 2\beta is not a coincidence and not a fit; it is what “same centroid” means for a rectangle.

Which free body produced the number

The compression zone above the neutral axis, taken as a free body with the tension steel’s force below it. Horizontal equilibrium fixes the neutral axis depth, and the moment about the steel gives the capacity.

For the section drawn above — 300 by 550 effective, C30, 1,800 mm² of grade 500 steel — the steel pulls 783 kN, the neutral axis settles at 189.6 mm, and the moment comes out 368.7 kNm. Solved with the curve, integrated over two thousand strips: 189.56 mm and 368.73 kNm. Solved with the rectangle, in three lines of arithmetic: 189.56 mm and 368.73 kNm.

Not close. Identical, to the precision the numbers are printed at, because the two distributions were built to have the same two integrals and the two integrals are all that equilibrium and moment use.

Wrong in shape, right in two integrals. The compression zone of a C30 section with its neutral axis 189.6 mm down, drawn twice. The curved outline is the real parabolic-rectangular stress distribution — the material's own law read off the linear strain profile plane sections supplies. The rectangle over it is what every design office uses instead: intensity η f_cd = 16.5 MPa over a depth λx = 158 mm. The two shapes are visibly different and give the same answer, because a bending calculation asks a stress distribution only two questions — how much compression there is, and where its resultant acts. Both are 783 kN at 78.9 mm from the face. The factors are α = 0.8095 and β = 0.4160, and λ = 2β follows from wanting the same centroid. A triangle and a full rectangle match neither integral and are nowhere near.
Fig. 2 The same section redrawn at the neutral axis equilibrium actually put it at. The compression zone is now 189.6 mm deep, the rectangle over it is λx = 158 mm at η f_cd = 16.5 MPa, and both distributions give 783 kN acting 78.9 mm from the face. That resultant is the steel’s pull turned round, and the lever arm is 550 − 78.9 = 471.1 mm — which over that force is the 368.7 kNm above.

Two shapes that do not work, and by how much

The claim that this is about the integrals rather than about smoothness is testable, so it is worth testing. Take the same compression zone and put two other plausible distributions on it.

A triangle — the elastic distribution, which is what the section carried before it yielded. Resultant 0.5fcdbx0.5f_{cd}bx, centroid at x/3x/3 from the face. It gives 383 kN at a lever arm of 500 mm: 36.7% low.

A full rectangle — the crudest plastic idealisation, everything at fcdf_{cd}. Resultant fcdbxf_{cd}bx, centroid at x/2x/2. It gives 765 kN at 475 mm: 20.3% high.

Both are smooth, both are simple, both sit over exactly the same compression zone, and both are useless. The equivalent rectangle is not in a different class because it is a better-shaped approximation. It is in a different class because it is the only one of the three that was constructed from the curve’s own integrals.

Neither of those two shapes gets a figure of its own here, and the reason is the argument rather than economy: there is nothing to draw. Both are one line over the same compression zone, both are computed above in a sentence each, and the picture that settles the matter is the first one — where the curve and the rectangle sit over one another with the same resultant at the same height.

How the neutral axis is found, which is where the rectangle earns its keep

The exactness above is a statement about a given neutral axis. In a real calculation the neutral axis is not given: it is the unknown, and it is found from horizontal equilibrium.

With the curve, that means guessing xx, integrating the compression over the zone, comparing with the steel’s tension, and iterating. It is four lines of code and it was never four lines of anything else. With the rectangle the compression is ηfcdbλx\eta f_{cd} b \lambda x — linear in xx — so equilibrium is a single division:

x=Asfydηλfcdbx = \frac{A_s f_{yd}}{\eta \lambda f_{cd} b}

and the moment follows directly. That is the actual saving, and it is larger than the saving on the moment: the moment was one integral either way, while the neutral axis was a root-finding problem and is now not.

It is also why the block survived the arrival of computers. A closed form for the neutral axis is what lets a designer see how the answer moves when the steel changes, and a section calculation that has to be re-run to answer “what if” is a different tool from one that can be read.

Wrong in shape, right in two integrals. The compression zone of a C30 section with its neutral axis 273.9 mm down, drawn twice. The curved outline is the real parabolic-rectangular stress distribution — the material's own law read off the linear strain profile plane sections supplies. The rectangle over it is what every design office uses instead: intensity η f_cd = 16.5 MPa over a depth λx = 228 mm. The two shapes are visibly different and give the same answer, because a bending calculation asks a stress distribution only two questions — how much compression there is, and where its resultant acts. Both are 1131 kN at 113.9 mm from the face. The factors are α = 0.8095 and β = 0.4160, and λ = 2β follows from wanting the same centroid. A triangle and a full rectangle match neither integral and are nowhere near.
Fig. 3 The same section with 2,600 mm² of steel instead of 1,800 — the one division run again. The neutral axis goes to 273.9 mm, which is 189.6 × 2600/1800, and the block is λx = 228 mm carrying 1,131 kN at 113.9 mm from the face. The lever arm has shortened from 471.1 mm to 436.1 mm, so 44% more steel buys 34% more moment: 493 kNm against 369.

Where the numbers came from, which was not a derivation

The construction above is how the factors are justified now. It is not how they arrived.

Charles Whitney proposed the equivalent rectangle in 1937 as a fit: he had a body of test results on beams failing in compression, he wanted an arithmetic that reproduced them without an integration, and he found a pair of factors that did. The two-integral argument came afterwards, as the explanation of why the fit was so good.

That order of events is worth knowing for a general reason. A rule that is a fit and a rule that is a theorem behave identically inside their range and completely differently outside it — and the whole of the difference shows up at the edges. Whitney’s factors held for sixty years and then began to move, above C50, exactly where the tests he fitted them to had no specimens. The theorem-shaped account says why: the shape of the normalised curve had changed, and a factor that describes a shape cannot survive the shape changing.

What the factors are properties of, which is not the concrete

Compute α\alpha and β\beta at C20, C25, C30, C35, C40, C45 and C50 and they come out identical: 0.8095 and 0.4160 at every grade.

That is not a rounding. The design strength fcdf_{cd} multiplies the whole curve, and both integrals are normalised by it — one divided by fcdbxf_{cd}bx and one a ratio of two integrals — so a uniform scaling of the stress axis cancels out completely. The factors describe the shape of the normalised curve and nothing else.

What does move them is the shape itself:

what changes α\alpha β\beta
the standard curve, εcu=0.0035\varepsilon_{cu} = 0.0035, n=2n = 2 0.8095 0.4160
a cubic rather than a parabola, n=3n = 3 0.8571 0.4357
a lower crushing strain, εcu=0.0030\varepsilon_{cu} = 0.0030 0.7778 0.4048

Both of those are real variations. High-strength concretes above about C50 are less ductile, crush at a lower strain and have a less pronounced plateau — which is exactly why every code that goes above C50 starts varying λ\lambda and η\eta with the grade at that point, and not before. The variation is not the strength arriving in the formula; it is the strength changing the shape of the curve.

Wrong in shape, right in two integrals. The compression zone of a C20 section with its neutral axis 150 mm down, drawn twice. The curved outline is the real parabolic-rectangular stress distribution — the material's own law read off the linear strain profile plane sections supplies. The rectangle over it is what every design office uses instead: intensity η f_cd = 11.0 MPa over a depth λx = 125 mm. The two shapes are visibly different and give the same answer, because a bending calculation asks a stress distribution only two questions — how much compression there is, and where its resultant acts. Both are 413 kN at 62.4 mm from the face. The factors are α = 0.8095 and β = 0.4160, and λ = 2β follows from wanting the same centroid. A triangle and a full rectangle match neither integral and are nowhere near.
Fig. 4 The opening section’s compression zone at C20 instead of C30. The rectangle is thinner — η f_cd is 11.0 MPa rather than 16.5 — and the resultant falls from 619 kN to 413 kN. The two factors do not move at all: α = 0.8095 and β = 0.4160, and the resultant still acts 62.4 mm from the face.

The interesting comparison is the other end of the range rather than the middle of it, because a factor that survives a change of 50 per cent in the design strength is a factor that is not about strength.

Wrong in shape, right in two integrals. The compression zone of a C50 section with its neutral axis 150 mm down, drawn twice. The curved outline is the real parabolic-rectangular stress distribution — the material's own law read off the linear strain profile plane sections supplies. The rectangle over it is what every design office uses instead: intensity η f_cd = 27.6 MPa over a depth λx = 125 mm. The two shapes are visibly different and give the same answer, because a bending calculation asks a stress distribution only two questions — how much compression there is, and where its resultant acts. Both are 1032 kN at 62.4 mm from the face. The factors are α = 0.8095 and β = 0.4160, and λ = 2β follows from wanting the same centroid. A triangle and a full rectangle match neither integral and are nowhere near.
Fig. 5 The same 150 mm compression zone at C50. The resultant is 1,032 kN against C20’s 413 — two and a half times as much — and α is 0.8095 and β is 0.4160 to the last digit printed, exactly as at C20 and C30. The block depth λx is 125 mm in all three, and the resultant sits 62.4 mm from the face in all three.

Two and a half times the stress and not a digit of movement in either factor. Now change the shape of the curve rather than its height, by taking the concrete to a crushing strain of 0.30% instead of 0.35% — which is what a genuinely high-strength concrete does — and the same two integrals come back different.

Wrong in shape, right in two integrals. The compression zone of a C30 section with its neutral axis 150 mm down, drawn twice. The curved outline is the real parabolic-rectangular stress distribution — the material's own law read off the linear strain profile plane sections supplies. The rectangle over it is what every design office uses instead: intensity η f_cd = 16.3 MPa over a depth λx = 121 mm. The two shapes are visibly different and give the same answer, because a bending calculation asks a stress distribution only two questions — how much compression there is, and where its resultant acts. Both are 595 kN at 60.7 mm from the face. The factors are α = 0.7778 and β = 0.4048, and λ = 2β follows from wanting the same centroid. A triangle and a full rectangle match neither integral and are nowhere near.
Fig. 6 The C30 section again with the crushing strain at 0.0030 rather than 0.0035. The curve now spends less of the zone on its plateau, so α falls to 0.7778 and β to 0.4048, the block shortens from 125 mm to 121 mm, and the resultant drops from 619 kN to 595 kN on identical concrete at an identical neutral axis. The strength did not change; the shape did.

What the rounding costs, which is the only real error

The codes do not print 0.832 and 0.973. They print 0.8 and 1.0, and for a good reason — nobody wants a lever-arm calculation with four significant figures in it when the concrete strength is known to one.

That rounding is a second approximation on top of a construction that was exact, and it is the only place in this whole procedure where an error actually lives. On the section above it is 0.69%, and it is on the safe side. It stays under one per cent across the useful range of neutral axis depths, and it changes sign nowhere — because λ\lambda was rounded down and η\eta up, and the two errors partly cancel in the resultant while adding in the lever arm.

Which is worth saying plainly, because the received account has it backwards. The stress block is usually described as “an approximation that works well”. It is an identity that has been rounded, and the rounding is the approximation. If the printed factors were 0.832 and 0.973 there would be no error at all.

Where the model stops

Four places, and the first two are ordinary while the third is the one that catches people.

Above C50. The curve’s shape changes with the grade and the factors have to vary with it, as described above.

Under axial load. Everything here assumed the strain profile runs from zero at the neutral axis to εcu\varepsilon_{cu} at the face. A section carrying substantial axial compression may be entirely in compression, with the strain profile a trapezium rather than a triangle, and the factors for that are different ones — which is where the second branch of every column interaction diagram comes from.

Where the compression zone is not rectangular. The construction gives α\alpha and β\beta for a stress distribution, and turns them into a rectangle for a section of constant width. A tee beam whose neutral axis falls in the web, a circular column, a section with a chamfer — for all of these the equivalent rectangle has to be applied to the actual width at each depth, and quoting a single λx\lambda x without checking the width over that depth is a common and quiet error.

In the serviceability calculation, where it does not apply at all. The block is an ultimate-limit-state device. At working load the concrete is nowhere near crushing, the distribution really is close to linear, and the section to use is the cracked transformed one. Using a stress block to compute a deflection is using the wrong end of the material’s curve.

The neutral axis is wherever the first moment vanishes. A 300 by 611 section with 1800 mm² of steel at a depth of 550, carrying 150 kNm after the tension side has cracked. The neutral axis is no longer at mid-depth: it has risen to 182.0 mm from the top, which is where the first moment of the compression zone plus the transformed steel vanishes. The compression is 11.2 N/mm² at the top fibre and the steel carries 170 N/mm²; the resulting couple is 307 kN on a lever arm of 489 mm, which multiplies back to the 150 kNm applied. The uncracked section would have had 6363×10⁶ mm⁴ against the cracked 2431×10⁶ — a loss of 62% of the stiffness.
Fig. 7 The other section, for the other question. At working load the stress is linear and the neutral axis is where the first moment of the transformed section vanishes — a different calculation with a different answer, on the same beam.

The third of those deserves a worked number rather than a sentence, because it is the one that produces wrong numbers in practice rather than merely imprecise ones. A tee beam with a 150 mm flange and a neutral axis at 190 mm has a stress block that is 150 mm of full width and 40 mm of web — and a designer who applies λx=158\lambda x = 158 mm at the flange width has assumed 158 mm of a section that is only 150 mm wide there. The error is in the direction of over-estimating the capacity, and it grows as the neutral axis approaches the flange boundary from below.

That 158 mm is the λx drawn earlier in this essay, on the same C30 section at the same 189.6 mm neutral axis, and nothing in the construction knows what width it is being applied over. The block has to be split at the change of width and the two parts summed separately, and a calculation that does not split it is over-estimating the capacity by the area of flange that is not there.

What the picture cannot show

The figure draws a stress distribution as though the compression zone were a continuum with a stress at every point. It is not. Concrete at 0.35% strain in a real beam is a network of cracks and crushed patches, and the “stress” plotted is a smeared quantity over a region a good deal larger than the aggregate. The curve is a description of a specimen’s average behaviour, transferred to a region of a member on the assumption that the region behaves like a small specimen — which is precisely the assumption the size effect says is not safe in general.

It happens to be safe here, because the compression zone of a bending member is confined by the material around it and fails in a much more ductile way than a cylinder in a testing machine does. But the reason is a physical one about confinement rather than anything in the arithmetic, and it is why the same block cannot be used for a plain concrete member with nothing holding it together sideways. Confining the compression zone changes both the peak and the crushing strain, and so changes the factors — which is the one material intervention that moves them at all.

The calibration assumed a rectangle, and not every compression zone is one

The two factors were fixed by matching two integrals over a compression zone of constant width. Nothing in that match survives a zone whose width changes with depth, and several ordinary sections have one.

A circular column is the clearest. Its compression zone is a circular segment: widest at the neutral axis and narrowing to a point at the extreme fibre. A triangular section does the same, as does a T-beam bent the wrong way up — a cantilever, or the hogging region of a continuous beam, where the flange is in tension and the compression is in the narrowing web-plus-fillet.

In all of them the real stress distribution and the equivalent rectangle disagree, because the two integrals weight the width differently. The rectangle spreads a uniform stress over the top 80 per cent of the zone, which is the part that has least width; the real distribution is zero at the neutral axis and rises toward the face, so it also weights the narrow part — but not in the same proportion, and the resultant and its centroid come out slightly wrong.

The codes acknowledge it with a blunt correction: where the width decreases toward the compression face, reduce the design strength by ten per cent. That number is not a material factor, not a partial factor and not a statistical allowance. It is a correction to a geometric approximation, applied to a material strength because that is where there was a place to put it.

Which is worth knowing for what it says about reading a code. A clause that reduces fcdf_{cd} looks like a statement about concrete, and this one is a statement about the shape of an integral — the same kind of substitution the essay’s opening was about, arriving in the opposite direction and much less honestly labelled.

The thing the block quietly decides, which is not the moment

There is one output of a section calculation that the equivalent rectangle does not reproduce, and it matters more than any of the small errors above.

The rectangle gets the resultant and the centroid right. It does not get the strain distribution right, because it does not have one — it is a stress shape with no material behind it. So anything that depends on how far the concrete has actually been strained is outside what the block can say: the curvature at failure, the rotation capacity of a hinge, whether the section is ductile enough to redistribute, and how much warning it gives.

Those are exactly the properties a plastic analysis depends on, which is why the two calculations are usually run with different machinery. The moment comes from the block in three lines; the rotation capacity comes from a full integration of the real curve, and the two live in different parts of the same design.

It is also why a section that “just passes” on a block calculation can be a poor section. The block will report the same moment for a heavily reinforced section with its neutral axis at 0.6d0.6d as any other method would, and will say nothing whatever about the fact that such a section crushes before its steel yields.

What it costs to reach the plastic moment, for two shapes. Moment against curvature for two cross-sections of identical area (165000 mm²) and identical depth (550 mm), in concrete, 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 — times the curvature at first yield; The I-section has a shape factor of 1.54 and reaches 98% of its plastic moment at — 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.
Fig. 8 The same beam over-reinforced. The moment is higher and the curvature at failure is a fraction of what it was — a difference the stress block cannot report, because it has no strain in it.

The generalisation, which is worth more than the stress block

The habit is this: before approximating a distribution, find out what the calculation actually asks of it.

A moment asks two integrals. A shear check asks a different pair — a first moment over a second moment — which is why the shear stress distribution cannot be replaced by a rectangle even though the direct stress can. A deflection asks for the whole curvature field along the member, which is why it needs the real stiffness rather than an equivalent one. A fatigue check asks for the peak, which no smeared distribution contains at all. And a stability check asks whether the compression zone can reach any of these stresses without buckling, which is a question about the geometry of the region rather than about the stress in it.

Once the question is named, the approximation is usually obvious and often exact. The stress block is the best-known instance and it is a good one to carry, because it is the case where the approximation looks least like the truth and reproduces it best — which is a useful corrective to the instinct that a picture which looks right is a picture that computes right.

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

Every essay whose body links to this one.

The objects this essay names

Each one links to every other essay that touches it.

Bending stressCentroidConcreteCrushing strainDesign strengthEquilibriumLever armMoment-curvatureNeutral axisPlane sectionsReinforcementResultantSectionStress blockUltimate strength