Materials

The strength no specimen had

A material property is written into a calculation as a number, and a material does not have one. It has a population of strengths with a mean and a spread, and the number used is a low fractile of that population — a value that need not have been measured, that most of the material exceeds, and whose distance below the mean is decided entirely by the scatter.

Assumes The bigger one is the weaker one, The flaw that sets the strength and The stress at which nothing in particular happens.

Every calculation in this collection has divided something by a strength. The strength was a number, and the material does not have a number.

It has a population. Cast thirty cubes from one mix, cure them together and crush them, and thirty different answers come back, spread over a range of perhaps a third of their mean. The number used in design is one particular point on the distribution those thirty results describe, chosen by a rule, and the choice is more consequential than almost anything else in a set of material properties.

The characteristic strength, which nothing was measured at. A lognormal population of strengths with a mean of 30 N/mm² and a coefficient of variation of 0.15. The characteristic value is the 5% fractile — 23.2 N/mm², which is 77% of the mean, and which need not be the strength of any specimen that was tested. Dividing it by 1.50 gives 15.5, and the shaded sliver below that is the fraction of the population that would fail to reach it: 6.4e-6, or one in 155,818. A factor applied to a fractile is not covering the scatter, because the scatter has already been spent getting to the fractile.
Fig. 1 A population of strengths with a mean of 30 N/mm² and a coefficient of variation of 0.15. The characteristic value is the 5% fractile and the design value is that divided by a partial factor.

Which free body produced the number

There is no free body. This is the one essay in the collection whose quantity comes from a distribution rather than from an equilibrium, and it is worth naming that as a change of subject rather than sliding past it.

What replaces the free body is a rule for choosing a point on a distribution. The rule almost everywhere is the 5% fractile: the value 95% of the population exceeds. Fit a lognormal to a mean mm and a coefficient of variation VV — lognormal because a strength cannot be negative and a normal fitted to 20% scatter puts real probability below zero — and the fractile has a closed form:

xk=m exp⁡ ⁣(−12σln⁡2−1.645 σln⁡),σln⁡=ln⁡(1+V2)x_k = m\,\exp\!\left(-\tfrac{1}{2}\sigma_{\ln}^2 - 1.645\,\sigma_{\ln}\right), \qquad \sigma_{\ln} = \sqrt{\ln(1+V^2)}

At V=0.15V = 0.15 that is 77% of the mean. At V=0.30V = 0.30 it is 60%. The whole of the gap between a material’s average behaviour and the number a designer is allowed to use is a function of one quantity, and that quantity is the spread.

Scatter beats mean, and it is not close

The consequence is the most useful thing in this essay.

Cutting the scatter beats raising the mean. The characteristic strength as a fraction of the mean, against the coefficient of variation, for a lognormal population at the 5% fractile. At a scatter of 0.20 the characteristic value is 71% of the mean; at 0.10 it is 84%. That improvement is worth the same as raising the mean strength by 19% and leaving the scatter alone — which is usually the more expensive of the two, and always the one that gets proposed. Nothing about the material's best specimens has changed in either move; the whole of the difference is in the tail.
Fig. 2 The characteristic strength as a fraction of the mean, against the coefficient of variation. At 0.20 it is 71% of the mean and at 0.10 it is 84%.

Halving the coefficient of variation from 0.20 to 0.10 takes the characteristic value from 71% of the mean to 84% — a gain of 19%, which is exactly what raising the mean strength by 19% and leaving the scatter alone would have bought.

Those two routes to the same answer have wildly different costs. Raising a mean concrete strength by 19% means more cement, a lower water-cement ratio, and a mix that is harder to place. Halving the scatter means better batching control, better curing, and more consistent testing — a process improvement rather than a material one, usually cheaper, and one that improves everything else at the same time.

It is also the route almost nobody proposes, because the specification names a strength and the scatter is nobody’s line item. The same asymmetry appears wherever this argument does: the mean is a property somebody sells and the spread is a property somebody controls.

Three specimens cannot see the tail

The second consequence is about testing, and it catches people who reason from a small sample as though the sample were the population.

Three specimens cannot see the tail. The factor k applied to the sample's own scatter when a characteristic value is estimated from n specimens. With the scatter known in advance it is z·sqrt(1 + 1/n) and barely moves; with the scatter estimated from the same n results it is the Student t quantile instead, and it runs from 7.73 at two specimens to 1.73 at 30. At n = 3 the characteristic strength comes out at 18.1 N/mm² against 23.2 for a population known exactly — 22% lower, for a material that is identical. A small test programme does not report a worse estimate of the strength; it reports a worse strength.
Fig. 3 The factor applied to the sample’s own scatter when a characteristic value is estimated from n specimens. With the scatter known it barely moves; with it estimated from the same results it is a Student t quantile.

If the population’s scatter is known in advance, the estimator uses z1+1/nz\sqrt{1 + 1/n} — 1.90 at three specimens against 1.645 at infinity, so a small sample costs little. If the scatter has to be estimated from the same nn results, the normal quantile is replaced by a tt quantile and the factor runs from 7.73 at two specimens to 1.73 at thirty.

At three specimens the characteristic strength of this material comes out at 18.1 N/mm² against 23.2 for a population known exactly — 22% lower, for a material that is identical. That is not a wider confidence interval; it is a lower number, entering a design as a strength.

That is why an approval by testing is expensive in a way that surprises people: the cost is not the tests, it is the number of them needed before the statistics stop taking the answer away.

It has a practical corollary that comes up on every project where something unusual is tested: a three-specimen programme on an excellent material can report a worse characteristic strength than a thirty-specimen programme on a mediocre one. The rule is not being unfair. It is refusing to believe a tail it has not seen.

The chain, which is the same argument about geometry

The third consequence connects this essay to one about size.

The chain is weaker than its links, and by how much is computable. The median strength of the weakest of n identical elements, as a fraction of one element's median, for a population with a coefficient of variation of 0.18. Nothing about the material changes along this axis: the distribution of the minimum of n draws is 1 − (1 − F)ⁿ, and its median is the fractile 1 − 0.5^(1/n) of a single draw. A chain of 250 is 39% weaker than one link, and the fall is slow because it is logarithmic — which is the same statement as the statistical size effect, arrived at without mentioning size at all. The scatter is the whole of the mechanism: at zero scatter the line would be flat.
Fig. 4 The median strength of the weakest of n identical elements, as a fraction of one element’s median. A chain of 250 is 39% weaker than one link, and the fall is logarithmic.

A structure made of many similar elements fails at the weakest, and the distribution of the minimum of nn draws is 1−(1−F)n1 - (1-F)^n. Its median is the 1−0.51/n1 - 0.5^{1/n} fractile of one draw, which for 250 elements at a coefficient of variation of 0.18 is 61% of one element’s median.

Nothing about the material changes along that axis. The scatter is the whole of the mechanism: at zero scatter the line is flat, and a chain of any length is exactly as strong as one link.

Two explanations that agree about the direction and nothing else. The same size effect under two theories. Weibull's is statistical — a larger specimen holds more flaws and fails at the worst one — and gives a straight line on log axes with slope −2/8, which has no size in it anywhere and so predicts a strength that falls forever. Bažant's is energetic: a crack releases energy in proportion to a volume and consumes it in proportion to an area, so there is a size at which the two balance and the curve bends from a plateau onto the −½ slope of fracture mechanics. They differ by up to 175% across this range, and the difference is not a detail: only one of them says where the transition is.
Fig. 5 Two explanations of the size effect. Weibull’s is this statistical argument and gives a straight line on log axes that predicts a strength falling forever; the energetic one bends from a plateau onto the −½ slope of fracture mechanics.

That is the bigger one being the weaker one arrived at without mentioning size at all — and it is why the two explanations of the size effect are different theories rather than variants. The statistical one is this argument; the energetic one is about crack growth. They agree about the direction and disagree by up to 175% about everything else.

What the partial factor is not doing

The last consequence is about the number that comes after the fractile.

Divide the characteristic value by a partial factor of 1.5 and the result sits at a fractile of 6.4×10−66.4\times10^{-6} — one specimen in 155,818. No test programme in the world has established anything about that part of a distribution, and no lognormal fitted to thirty results is a description of it.

So the factor is not covering the scatter, because the scatter was already spent getting to the 5% fractile. What it covers is a different list: the difference between a specimen and the material in the structure, the difference between a laboratory test and a real load, the geometric tolerances, and the possibility that the model used is wrong. Those are model uncertainties, and they are not distributed like a material property at all.

The clearest evidence for that reading is the variety of the factors. Concrete carries 1.5 and steel 1.0 in the same code, on materials whose scatters are 0.15 and 0.05 — a ratio of three where the scatters differ by three, which looks like agreement until the fractiles are computed. The characteristic values already differ by the scatter; applying a further factor in the same ratio applies the correction twice. What actually separates them is that concrete in a structure is not the concrete that was tested, and steel in a structure is.

A strength that is a property of the specimen. Nominal strength against size for geometrically similar specimens of one material. On the left the specimen is too small for a crack to run and the strength is a plateau — a plastic limit, and the regime laboratory specimens sit in. On the right a crack releases more energy than it consumes as soon as it starts and the strength falls as the inverse square root of size, which is the regime real structures sit in. The turn happens at D₀ = 80 mm. A 100 mm specimen reads 2.80 N/mm² and a 900 mm member of the same material carries 1.20: the test overestimates the structure by a factor of 2.33.
Fig. 6 One item on that list, made explicit. A 100 mm specimen reads 2.80 N/mm² and a 900 mm member of the same material carries 1.20 — the test overestimates the structure by a factor of 2.33.

The number the whole system is actually aiming at

The 5% fractile and the partial factor look like two independent conventions. They are not: both are calibrated against a single target, and knowing what the target is resolves the puzzle the previous section left open.

The target is a reliability index β\beta, defined so that the probability of the resistance falling below the load over a fifty-year reference period is Φ(−β)\Phi(-\beta). The usual value is 3.8, which is a probability of about 7×10−57\times10^{-5} — one structure in fourteen thousand, over its life, for the combination of everything going wrong at once.

That target is then split between the two sides. A sensitivity factor of about 0.8 is assigned to the resistance and 0.7 to the load, so the resistance’s design value is placed at 0.8β=3.040.8\beta = 3.04 standard deviations below its mean, which is a fractile of 1.2×10−31.2\times10^{-3} — one in eight hundred.

Compare that with the 6.4×10−66.4\times10^{-6} the material’s own distribution gave. They differ by a factor of about two hundred, and the difference is not an error in either. The reliability target is on the resistance of a member, and a member’s resistance is more variable than its material: it carries the geometric tolerances, the difference between the specimen and the in-situ material, and the error in the model that converts a strength into a capacity. A resistance with a coefficient of variation of 0.20 sitting at its one-in-eight-hundred fractile is at a much less extreme point of a much wider distribution — and lands at a similar absolute value.

A scatter that belongs to the site, not to the material. The as-built effective depth of a 200 mm slab, and the capacity that follows from it. The drawing says 169 mm. The distribution is centred 4 mm below it, because reinforcement in a slab is walked on and pushed down far more often than it is lifted, and it has a standard deviation of 8 mm. Two losses follow and they are separate: the bias takes 2.5% off the capacity before any scatter is considered, and the 5% fractile takes 10.8%. Neither appears in any material partial factor, because neither is about the material — the concrete and the steel are exactly as strong as they were tested to be, and the member is 10.8% weaker than the drawing says because a dimension nobody can measure after the pour came out different.
Fig. 7 One item on that wider list, drawn as its own distribution. The drawing puts the effective depth of a 200 mm slab at 169 mm; the as-built population is centred 4 mm below that, with a standard deviation of 8 mm, because reinforcement is walked on far more often than it is lifted. The bias alone costs 2.5% of the capacity and the 5% fractile costs 10.8%.

That figure is the clearest available demonstration of what the material factor is not for. The concrete and the steel in that slab are exactly as strong as they were tested to be — every specimen behaved, every cube reached its mean — and the member is 10.8% weaker than the drawing says because a dimension that nobody can measure after the pour came out different. The scatter is real, it is large, and it belongs to the site rather than to the material. It has its own bias too, which a symmetric material factor does not resemble at all: a strength distribution is centred on what the material does, and this one is centred 4 mm below what the drawing asked for, because the error has a direction.

So the partial factor is doing two jobs at once, and the arithmetic separates them cleanly. Part of it moves along the material’s distribution, and that part is small. The larger part is buying coverage of a different and wider distribution that no material test has anything to say about. Which is why a factor of 1.0 on steel and 1.5 on concrete is not a statement about their scatters — those were already spent — but about how much of the member’s resistance the specimen failed to describe.

And the target is not one number either

The 3.8 is itself a choice, and it varies with what happens if the structure fails.

Codes carry three consequence classes: a lower target of about 3.3 where failure has small consequences, 3.8 for ordinary buildings, and 4.3 where it would be severe. In probability those are one in two thousand, one in fourteen thousand and one in a hundred thousand over fifty years — a spread of fifty.

The interesting part is where the adjustment is applied. It is not applied to the material. The characteristic strength stays at the 5% fractile and the partial factor stays where it was; what changes is a multiplier on the loads, of 0.9, 1.0 and 1.1.

That is a deliberate choice and it is the right one. A material’s distribution is a fact about the material and does not know what the structure is for; the consequence of failure is a fact about the structure and nothing to do with the concrete. Putting the adjustment on the load side keeps each factor a statement about the thing it is derived from — which is the same discipline this whole essay is about, applied to the factors rather than to the strengths.

Where the number comes from for other properties

Not every material property is a strength, and the fractile rule is applied differently to each.

The modulus is used at its mean, because stiffness matters for load-sharing between parts and a conservative estimate for one member is unconservative for its neighbour.

The tensile strength of concrete — the property a flaw sets and the one with the widest scatter of any in common use — is quoted at three fractiles at once — a lower one for anchorage, a mean for deflection, an upper one for capacity design — because which end of the distribution is unsafe depends on the check.

The yield stress of steel is quoted at a low fractile for strength and at a high one for capacity design, for the same reason: a beam that is stronger than expected forces its connection to be stronger still.

That last case is worth pausing on because it compounds the problem. The stress at which nothing happens is a definition rather than a measurement for every material without a plateau — a line of slope EE from a strain of 0.002, cutting the curve wherever it does. Mild steel has a plateau, so its yield is an event the specimen performs and a testing machine records; a high-strength steel has none, and its yield stress is wherever a construction line happens to cut a smooth curve. So for those materials the population being fractiled is a population of results of a construction, not of an event, and part of its scatter is the scatter of the construction.

Which materials have which spread

The coefficient of variation is the quantity everything in this essay turns on, and its value is a property of the material and its manufacture rather than of the science.

Steel is the tightest: about 5% on yield stress for a rolled section, because it is made under control in a plant with continuous testing and a specification that the producer has to hit. A characteristic value 92% of the mean is the reward.

Concrete is around 10 to 20% depending on the batching, which puts its characteristic value between 71 and 84% of the mean — the whole range of the figure above, decided by how the material is made rather than by what it is.

Timber is 20 to 30% before grading and much less after it, which is the entire purpose of a grading system: it does not improve any piece of timber, it sorts the population into narrower ones. A grade is a scatter-reduction device, and a material with a direction needs it more than most because its properties depend on where the tree’s defects happen to be.

Masonry is the widest in ordinary construction, at 25% or more, because the unit, the mortar and the workmanship all vary and the third is not measured at all.

The characteristic strength, which nothing was measured at. A lognormal population of strengths with a mean of 30 N/mm² and a coefficient of variation of 0.25. The characteristic value is the 5% fractile — 19.4 N/mm², which is 65% of the mean, and which need not be the strength of any specimen that was tested. Dividing it by 1.50 gives 12.9, and the shaded sliver below that is the fraction of the population that would fail to reach it: 5.0e-4, or one in 2,007. A factor applied to a fractile is not covering the scatter, because the scatter has already been spent getting to the fractile.
Fig. 8 The first figure’s population again, with nothing changed but the scatter: a coefficient of variation of 0.25 rather than 0.15. The characteristic value falls from 77% of the mean to 65%, and the design value that the same partial factor of 1.5 produces now sits at one specimen in 2,007 rather than one in 155,818.

That pair is the whole essay in two pictures of the same population. The rule did not change, the factor did not change and the mean did not change; the scatter went from 0.15 to 0.25 and the design value moved from a fractile no test programme could reach to one that a run of two thousand units would be expected to contain. A masonry wall and a concrete column carrying the same nominal factor are therefore not being treated alike, and nothing in either calculation says so.

Read against the 19%-for-halving rule, that list explains a great deal about how the four materials are used. Steel is worked close to its mean and its members are slender; masonry is worked far below its mean, which is one reason a masonry wall’s calculated capacity looks so unlike what a wall obviously does.

It also explains why a stronger steel does not change the modulus but does change the scatter. Higher-strength grades are made to tighter specifications, so their characteristic values sit closer to their means, and part of the gain from specifying a stronger steel is not strength at all.

Where the scatter comes from

It is worth separating the sources, because they respond to different interventions.

The material itself. Concrete’s aggregate, cement and water vary batch to batch; timber contains knots; steel’s chemistry varies within a specification band.

The specimen. Its size, its curing, its surface, its alignment in the machine.

The test. Its rate of loading, its temperature, the operator. Loading rate alone is worth several per cent on a yield stress and much more on concrete — a steel that is stronger in a millisecond is the same steel, measured differently.

Of these the second and third are usually the largest for well-made materials, which is why standardising a test does more for a characteristic value than improving a material does. It is also why a material tested one way and used another has a characteristic strength that describes neither.

There is a further complication for any material with several strengths that are not exchangeable. A timber beam is checked in bending, in shear and in bearing, and those are three populations with three means and three scatters, sorted by three different grading rules. Which of them arrives first is then a geometric question rather than a material one, and the check that governs does not move with the span at all — so a single “characteristic strength of the timber” is not a quantity that exists, and the phrase has to be finished with the word for what is being resisted before it means anything.

Where the model stops

The distribution is fitted, not observed. A lognormal through a mean and a coefficient of variation is a two-parameter description of thirty numbers, and the 5% fractile lies at the edge of the data. The fifth percentile of the fit and the fifth percentile of the sample are different things.

Strengths are not independent. The chain argument assumed the elements’ strengths were independent draws. Adjacent concrete in one pour, or adjacent bolts from one batch, are correlated, and a positive correlation makes the chain stronger than the formula says.

Nothing here is about consequences. A fractile is a probability of a value being exceeded; a design decision is about the cost of exceeding it. Two materials with the same fractile in a component whose failure is ductile and in one whose failure is sudden should not carry the same factor, which is why the property that appears in none of the equations is decisive here without ever being computed, and the treatment of that is a subject the fractile does not enter.

And a fractile is not the only way a code hides an allowance. Several limits that look derived turn out to exceed the theory they are quoted against by the same factor at every grade of steel, which is the signature of a knockdown applied to a theory rather than a theory of the knockdown. Those allowances behave like a partial factor and are not called one, so they do not appear in any reliability calibration and cannot be reasoned about the way the material factor can.

What the picture cannot show

A distribution drawn on a page has smooth tails, and the tails are the part nobody has measured. Everything to the left of about the tenth percentile in the first figure is an extrapolation of a functional form, and the design value sits several times further out than that.

Nor does the picture show which population is being described. The thirty cubes were made from one mix on one day by one gang. The material in the structure is a year of mixes, several gangs, and weather — a wider population with a larger scatter and a lower fractile, described by no test that was performed.

And it does not show that the quantity being fractiled may stop being the governing one altogether. A single tensile strength assumes the section yields, and it does — below a crack of about ten millimetres in a plate of ordinary toughness. Above that length the crack decides, the yield stress is irrelevant to the answer, and a distribution of yield stresses is a distribution of a number that no longer enters the calculation. The same is true across a stress state rather than along a crack length: one uniaxial strength describes a whole yield surface, and two criteria fitted to the identical uniaxial number disagree by up to about a sixth of it elsewhere, with the widest gap in pure shear. A characteristic value is a summary of a population of measurements of one loading, and the structure is not obliged to load the material that way.

The generalisation

The habit worth carrying is that every material property in a calculation is a summary of a distribution, and which summary depends on which way the answer is unsafe.

A strength is taken low because a weak specimen is what fails. A modulus is taken at its mean because both directions are unsafe somewhere. A yield stress for capacity design is taken high because the thing being protected is the rest of the structure. Three different summaries of three distributions, each chosen by asking what would be bad.

That question generalises well past materials, and it is the same one that decides whether a load is factored up or down, whether a friction coefficient is taken at its lower or upper bound, and whether a stiffness should be over- or under-estimated in a system with two paths. In every case the arithmetic is the easy half. The judgement is knowing which tail of which distribution the structure is standing on.

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 20 that link here.

The objects this essay names

Each one links to every other essay that touches it.

Characteristic strengthCoefficient of variationDuctilityFractileFractureLognormalPartial factorProof stressReliabilityScatterSection classificationSize effectTestingWeakest-linkYield criterion