Structural Engineering & Statics
where the load goes
Concept

Strain energy — where it appears

The work stored in a structure as it deforms, recoverable while it stays elastic, and the quantity every energy method is written in. Differentiating it with respect to a load gives the deflection at that load, which is Castigliano's theorem and the same sum the unit-load method writes another way.

Named by 12 essays across 4 fields — each of them below, with the objects they name alongside it.

20 kN applied at once and held, on a structure of 0.500 s period. Displacement against time for a single-degree-of-freedom structure of natural period 0.500 s and 2.0% damping, under 20 kN applied at once and held. The static deflection under the same peak force is 12.67 mm and the peak response is 24.56 mm — a factor of 1.94.

Twice the deflection, for the same load

A weight placed gently on a beam deflects it by one amount. The same weight let go from rest, a millimetre above the same beam, deflects it by twice as much — and the factor of two is exact, for every structure ever built.

dynamics · Dynamic amplification What a drop height is worth, as a factor on the answer for a weight placed slowly. The peak displacement as a multiple of the static deflection, against the height a weight is dropped from divided by the deflection that weight causes when it is placed. The curve is 1 + √(1 + 2h/δ), which is conservation of energy and nothing else: the weight does work over the height it falls PLUS the distance the structure then gives, and the structure stores work only over the second. At a ratio of 40 the factor is 10.00, and at zero it is exactly 2 — the marked point, where a dropped weight becomes a placed one.

The weight that was dropped

A half-tonne load lowered onto a beam produces 5 kN. The same load dropped one metre onto the same beam produces 160 kN — and onto a beam ten times softer, 54 kN. The stiff structure is the one that suffers, which is the opposite of nearly every other rule about structures.

dynamics · Impact factor Four guesses at one buckling mode. A pin-ended column, with four assumed shapes and the load each of them gives. The reference is a ten-term Ritz expansion solved as an eigenvalue problem, at 9.8696 EI/L² — which is π², as it must be. a half sine gives 9.870, its own sag shape gives 9.882, a mid-span sag gives 10.000, a parabola gives 12.000. Every one of them is high and none of them is low, because an assumed shape is a constraint on the column and a constraint can only stiffen it.

Guessing the shape, and getting the load anyway

A column's buckling load can be had from a shape that is wrong everywhere, because the energy criterion is stationary at the true mode. The error in the load is the square of the error in the shape, and it is always high.

stability · Stability energy A derivative taken with a ruler, and the step that makes it worst. Castigliano's theorem says the deflection is ∂U/∂P, and the derivative here is taken numerically — two solves at ±dQ and a central difference. Against the unit-load answer of 1.720635e-2 it agrees to 1.6e-13, which for a linear structure it must: ∂N/∂P is exactly the force a unit load produces, so the two expressions are the same sum written twice. The error against step size is the classic pair of straight lines — truncation falling as the step shrinks, round-off rising as the difference of two nearly equal energies loses its digits — meeting near dQ = 1.2e+1. For a linear structure the truncation term is exactly zero, so what is drawn here is round-off alone.

The deflection that is a derivative

A structure's strain energy is one number. Differentiate it with respect to a load and out comes the displacement under that load — and the trick that makes it a method rather than an identity is that the load does not have to be there.

deflection · Strain energy Half the beam does nearly all of the deflecting. The virtual-work integrand M·m/EI along the member, normalised to its own peak, with the running share of the answer beside it. The integrand is a density: it says how much of the deflection each millimetre of the beam produced. For this case the half nearest the root supplies 87.5 per cent of it, and the rest of the member supplies the remainder. Stiffening the busy 50 per cent by 1.5 times takes the deflection down by 29.2 per cent; the same material spent on the quiet end takes it down by 4.2 — a factor of 7.0 for the same steel. The map of what is contributing is not the map of where the moment is largest, and the second is the one that gets drawn.

Where a deflection comes from

The unit-load method gives a deflection as an integral, and this collection has treated that integral as a number to evaluate. It is not a number. It is a density, and it says which millimetres of the member produced the answer — which is not the same map as where the moment is largest.

deflection · Deflection distribution A section has two areas and the tables give one of them. Peak shear stress divided by the mean, for four sections of exactly the same gross area and depth. The mean is V/A and is the number a first calculation uses; the peak is what the material actually sees, and the ratio between them is a property of shape alone. A rectangle's is 1.5 — the parabola's peak over its average — and it is one of the few numbers in this subject that is exactly derivable and universally ignored. An I-section's is near 1.98, and the reason is on the second bar: 97% of the shear is inside a web that is 56% of the area. So the flanges carry the moment and almost none of the shear, and the web carries the shear and almost none of the moment — which is why a shear check on an I-section uses the web area and a moment check uses the whole section, and why the two checks are about two different pieces of steel.

The section has two areas

A shear force divided by the area of the section is not the shear stress anywhere in it. A rectangle's peak is exactly one and a half times that number and an I-section's web carries nearly all of the shear over a fifth of the area, which is why a moment check and a shear check on the same member are checks on two different pieces of steel.

sections · Shear area Maxwell's reciprocal theorem. A load at one point and the deflection it causes at another, against the same load moved to the second point and the deflection read at the first. Both integrals return 307.5006, and neither calculation was told about the other. The two deflected shapes are entirely different; the two readings are identical.

An influence line is a deflected shape

Finding where a load has to stand to be worst means solving the structure once for every position it could stand in. Reciprocity says the answer is a single deflected shape — release the quantity being asked about, move it by a unit, and the shape the structure takes is the influence line.

deflection · Reciprocity Every member's share of the movement, and they are not the members expected. A Pratt truss of eight panels at a depth of 1, carrying 15 kN at each top node, with the movement of the bottom chord at mid-span attributed member by member. The unit-load sum δ = ΣF·f·L/EA gives 2019.41 at EA = 1: 49.4% from six top chords, 30.8% from eight bottom chords, 16.8% from eight diagonals, 3.0% from seven verticals. The single worst member is a top chord at mid-span at 11.9% of the whole. Each member is drawn at the width of its own share. The same deflection from a stiffness solution that shares none of this arithmetic is 2019.41, a relative residual of 1.4e-14.

The member that is not worth stiffening

A truss's deflection is a sum of one term per member, and a term is zero whenever either force in its product is. A vertical carrying the whole of a panel load can contribute nothing at all to the movement — which a total can never show and a per-member sum shows nothing else.

deflection · Truss deflection The coefficient is an envelope, and its scallops are whole half-waves. The plate buckling coefficient against aspect ratio α = a/b. Each faint branch is one half-wave count m: k = (m/α + α/m)², a curve whose own minimum is exactly 4 at α = m. The plate takes whichever branch is lowest, so the answer is the bold envelope — four touching 4 at α = 1, 2, 3 and 4, with cusps between them at α = √(m(m+1)) = 1.41, 2.45, 3.46, 4.47, where the plate is indifferent between m and m+1 half-waves. The first cusp reaches k = 4.50 and every later one is lower — 4.17, 4.08, 4.05. Past α = 1 the envelope never exceeds 4.49, which is why the length of a plate drops out of a formula that is otherwise entirely geometry.

Four was never a fact about plates

The coefficient every plate calculation starts from is quoted as 4, derived nowhere and remembered by everyone. It is the minimum of a quantity that has nothing to do with plates in it — and what the plate supplies is not the four but the restriction that produces the scallops around it.

stability · Plate buckling Three routes to one deflection, and the one that is wrong. Two 2500 mm bars of 250 MPa proof stress meeting at a loaded apex, the drop of the apex against the load. The line is geometry: each bar's extension from its own stress, divided by the sine of its slope. The dots are the derivative of the total complementary energy with respect to the load, and lie on it. The dashed line is the derivative of the strain energy — Castigliano's theorem applied to a material that is not linear — which leaves the truth by ten per cent at 99 kN and at 170 kN gives 308 mm for a deflection of 46.0.

The other area under the curve

Castigliano's theorem says a deflection is the derivative of the strain energy with respect to the load, and it is true only while the material is linear. Past that, the right energy is the area on the other side of the stress–strain curve. On two aluminium bars at their proof stress the strain energy gives a deflection four times too large, and on a redundant truss minimising it picks a set of forces in perfect equilibrium that no deformed shape can produce.

deflection · Strain energy The same mode's curvature, which can. The first mode's drift in each storey — the difference between the floors above and below it, which for a shear building plays the part a beam's curvature plays — intact, solid, and with storey 3 10 per cent less stiff, dashed, each scaled to the roof. The damaged storey's drift rises by 9.6 per cent; the storeys above it change by at most 1.8 and below it by at most 1.5. The largest change is in storey 3: the drift finds the damage and says where it is, which the frequency's fall of 0.91 per cent cannot.

What a mode shape notices that a frequency does not

Take a tenth of the stiffness out of one storey of a ten-storey building and its first frequency falls by a per cent at most, and by almost nothing if the storey is near the top. The mode shape looks unchanged, its match to the old one 0.99997. But the mode's storey drift — its curvature — rises by about a tenth in the damaged storey and hardly anywhere else, wherever that storey is. The instrument that sees the damage is the one almost nobody installs.

dynamics · Mode shapes The strain-energy route halves the sag. The midpoint deflection of a cable 10.0 m long, pretensioned to 60 kN, with an axial stiffness EA of 64 MN, loaded across its span at midpoint, against the load, found three ways. Solid: the deflection itself, from the cable's geometry; the dots, the derivative of the complementary energy with respect to the load, lie on it. Dashed: the derivative of the strain energy. Dotted: the pretension alone, as if the cable's tension did not grow. At 10.0 kN the cable sags 212 mm; the strain-energy route gives 107 mm, 0.51 of it, and the pretension alone 417 mm. At 20.0 kN: 294, 128 and 833 mm. The strain-energy route understates the deflection of a member that stiffens by the same mechanism that made it overstate one that softens.

The sag the strain energy halves

The derivative of the strain energy overstates the deflection of a member that softens. Turn the curve the other way up — a cable that grows stiffer as it sags, a hanger that takes up its play — and the same derivative understates it: by half for a pretensioned cable carrying a modest load, by two thirds for one with no pretension, and by the whole of the slack for a member that has any. The error is one ratio in both directions, and the least-work shortcut built on it hands the load to the member that stiffens.

deflection · Strain energy

Named alongside it

The objects these essays reach for when they reach for this one.

Virtual workFree bodyDeflectionStiffnessCastiglianoComplementary energyConservation of energyEigenvalueReciprocityRedundantUnit loadBuckled mode shape

All concepts