Structural Engineering & Statics
where the load goes
Series

Strain energy — the series

3 essays on one idea, from the one that introduces it to the one that assumes the rest.
  1. 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.

    part 1 · deflection
  2. 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.

    part 2 · deflection
  3. 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.

    part 3 · deflection

All series