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 atA 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.0102030405060strength (N/mm²)mean 30characteristic 23.2design 15.5one specimen in 155,818 falls below the design value
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=mexp ⁣(12σln21.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 meanThe 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.00.050.10.150.20.250.30.3500.20.40.60.81coefficient of variationcharacteristic ÷ meanas tested: 71%at 0.10: 84%worth a 19% rise in mean strength
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 tailThe 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.510152025302345678specimens testedfactor on the sample's scatterscatter estimatedscatter knownthe population itselfk = 3.37 at three specimens
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 computableThe 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.11010010000.60.70.80.91elements in the chain (logarithmic)median ÷ one element's median61% at 250a chain of 250 keeps 61% of one link's median
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(1F)n1 - (1-F)^n. Its median is the 10.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 elseThe 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.1032100316100031620.40.61.62.54.0size (mm, logarithmic)nominal strength (N/mm², logarithmic)the energetic lawthe statistical oneD₀ = 90 mmslope -0.05 at the small end and -0.49 at the large
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×1066.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 specimenNominal 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.10321003161000316201234size (mm, logarithmic)nominal strength (N/mm²)the specimen: 2.80the structure: 1.20the plastic limitfracture mechanicsthe test overestimates by 2.33× · D₀ = 80 mm
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.

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.

Two materials pulled until they stopTwo stress-strain curves — mild steel, high-strength steel — plotted to a strain of 2.0%. One of them has a plateau, so the stress at which yielding starts is something the specimen does rather than something anyone chooses. The 0.2% offset construction is drawn on the high-strength steel: a line of slope E from a strain of 0.002, cutting the curve at 460 N/mm².00.5%1%2%2%0100200300400500600strainstress, N/mm²the 0.2% proof stress: 460 N/mm²mild steelhigh-strength steel
Fig. 7 Two materials pulled until they stop. One has a plateau, so its yield is an event the specimen performs; the other has none, and its “yield stress” is where a construction line cuts the curve.

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. So for those materials the population being fractiled is a population of results of a construction, not of an event.

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.

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.

A short timber beam is a shear problem, and a steel one never isUtilisation of the bending and shear checks on one beam, against span-to-depth. The two cross where the ratio equals f_m ÷ f_v exactly — no load, no width and no span survives the cancellation — which for this timber is 7.5 and for steel is 1.73. So a timber beam shallower than about six times its depth is governed by shear parallel to the grain, and a steel beam would have to be shorter than twice its own depth before the same thing happened, which is not a beam. The third check is bearing across the grain, which does not move with the span at all: on the beam drawn it is at 0.48 and it is the one that governs.5101520253000.20.40.60.8span ÷ depthutilisationbendingshearbearing across the grainthey cross at 7.5f_m ÷ f_v = 7.5 for timber and 1.73 for steel
Fig. 8 A material with several strengths that are not exchangeable. Bending, shear and bearing on this timber beam are three checks with three scatters, and which of them arrives first is a geometric question, and the one that governs does not move with the span at all.

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.

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 17.2, 15.2 at 275, 355 N/mm², against quoted limits of 12.9, 11.4; A web, in bending (buckling coefficient 4) derives to 52.5, 46.2 at 275, 355 N/mm², against quoted limits of 38.8, 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.4317.2 at 27515.2 at 355quoted: 14ε1.33× the quoted limit, at every gradeweb, in bendingk = 452.5 at 27546.2 at 355quoted: 42ε1.35× the quoted limit, at every grade0102030405060width ÷ thickness
Fig. 9 What a fixed knockdown looks like from outside. The derived limit exceeds the quoted one by the same factor at every grade, which is the signature of an allowance applied to a theory rather than a theory of the allowance.

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.

The crack length at which the strength stops matteringFailure stress against crack length for a toughness of 70 MPa√m, with one steel grade drawn. The falling curve is fracture — Kc divided by Y times the root of pi a — and it does not know what the yield stress is. The horizontal lines are the grades. At 355 N/mm² the two cross at a crack 9.7 mm long. Below that the section yields and the crack is irrelevant; above it the crack decides and the 355 is irrelevant.102030400100200300400500crack length, mmstress at failure, N/mm²355 N/mm² crosses at 9.7 mmfracture: the crack decides
Fig. 10 The other thing a single strength hides. Below a crack of 9.7 mm the section yields and the crack is irrelevant; above it, the crack decides and the yield stress is irrelevant.
One criterion inside the other, touching at six pointsThe two yield criteria in principal stress space with the third principal stress zero, both normalised by the yield stress. Von Mises is the ellipse — σ₁² − σ₁σ₂ + σ₂² = f_y², which is a circle seen at an angle — and Tresca is the hexagon inscribed in it, touching at the six points where one principal stress is zero or the two are equal. Everywhere else Tresca is the smaller, by up to 15.5 per cent, and the widest gap is at pure shear, where σ₁ = −σ₂ and the two answers are 159 and 138 N/mm². The ratio there is exactly 2/√3, computed rather than quoted, and it is the whole reason a web is checked against f_y over root three.-1-0.50.51-1-0.50.51σ₁ ÷ f_yσ₂ ÷ f_yvon MisesTrescapure shearthe widest gap is 2/√3 = 1.1547, at pure shear
Fig. 11 And the fact that one number describes a surface. Two criteria fitted to the same uniaxial strength differ by up to 15.5% elsewhere, with the widest gap at pure shear.

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.

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.

Characteristic strengthCoefficient of variationDuctilityFractileFractureLognormalPartial factorProof stressReliabilityScatterSection classificationSize effectTestingWeakest linkYield criterion