Structural Engineering & Statics
where the load goes
Concept

Flexibility — where it appears

Displacement per unit of force, the reciprocal of stiffness, and the quantity a compatibility equation is naturally written in. It is what adds along a load path, which is why a compatibility equation is written in flexibilities and a stiffness assembly is written in stiffnesses.

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

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 93.3335, and neither calculation was told about the other. The two deflected shapes are entirely different; the two readings are identical.

The theorem that swaps the question round

Push here and measure there; push there and measure here. The two readings are identical, for every elastic structure, whatever its shape — and that fact turns an influence line into something a model can be asked for directly.

deflection · Reciprocity 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 The tie is a redundancy, so its stiffness decides the thrust. Thrust and rib bending for a 60 m tied arch of 0.15 rise ratio, against the stiffness of its tie. Cut the tie and the structure is a curved simply supported beam, so the tie force is the one redundant and the force method gives it: with a rigid tie the answer is 2227 kN, within 1.0 per cent of the funicular wL²/8f, and the shortfall is the arch's own axial shortening. A real tie stretches 68 mm and returns 2166 kN — 2.7 per cent of the flexibility is the tie — and whatever thrust the arch does not get, it carries as bending: 760 kNm at 30 m. A tenth of the tie stiffness is not a tenth of the problem; it is a different structure.

The thrust that never reaches the ground

Every arch on this site has ended at the same sentence — the foundation is where an arch is really decided. A tie changes the sentence without changing the arithmetic: the horizontal force is still there, still the same size, and it now closes on itself through a bar at deck level.

structures · Tied arch Two different structures released, and one bending moment diagram. The bending moment in a continuous beam of 8, 10, 8 m under 12 kN/m, solved twice by the force method with different redundants. The first release puts a hinge over each interior support, so the released structure is a row of simple spans and the redundants are moments. The second removes each interior support, so the released structure is one simple span of the whole length and the redundants are reactions. The two released structures have nothing in common — different shapes, different deflections, different everything — and the diagrams they produce lie on top of each other to 9e-15 of the peak moment. Which restraints are released is a choice about the arithmetic and not about the structure, which is a fact worth trusting: it means a hand calculation can pick whichever release makes the sums easiest and be sure of the answer.

Choose what to take away

The other machine for a redundant structure works by removing restraints until what is left can be solved by statics, then putting back exactly enough force to close the gaps that opened. Which restraints are removed does not change the answer at all, and changes the arithmetic completely — one choice gives a tridiagonal matrix a person can solve on paper, and another gives a full one.

deflection · Force method One member's stiffness, scattered into the freedoms it touches. A member's own six-by-six stiffness matrix relates the forces at its two ends to the displacements there, and it is written in the member's own axes. Assembly is two operations and no physics: rotate it into the structure's axes, then add each of its thirty-six entries into the row and column of the global freedom that entry belongs to. Every member does the same, and the sum is the structure. The shaded rows and columns are the six freedoms this one member reaches; every other entry it contributes is exactly zero, and that is the whole reason a global stiffness matrix is sparse. Nothing here is an approximation — the result is the same equilibrium and the same compatibility a hand method writes, in an order a machine can follow.

The answer that depends on how it was divided

Every computed answer in this collection came out of a structure chopped into pieces — elements, strips, stations, trial positions. The chopping is invisible in the result and it is not neutral: some divisions give the exact answer, some give one that is always too stiff, and one of them changes nothing but the cost of getting there.

deflection · Discretisation A tenth of a per cent of the thrust is all of the moment. The thrust a two-hinged arch loses to its own axial shortening, against rise-to-span. The flexibility equation's denominator has two terms — ∫y²ds/EI for bending and ∫cos²θ ds/EA for shortening — and their ratio is about (15/8)(i/f)², the square of the radius of gyration over the RISE. The loss is a little more than that ratio, because the released rib also shortens under its own shear, and at the 10 per cent rise drawn it is 0.11 per cent of the thrust, which sounds like a rounding error and is not: a parabolic arch under a uniform load is funicular, so the rigid solution has NO crown moment at all, and the 0.11 per cent that the rib shortening removes from the thrust leaves 29 kNm behind. The correction that is a tenth of a per cent of the thrust is a hundred per cent of the bending. At a two per cent rise the loss is 2.6 per cent, because a shallow arch's thrust is enormous and its lever arm is not.

The arch that gets shorter

A parabolic arch under a uniform load is funicular, so the perfect solution gives it no bending at all. Then the rib shortens under its own thrust by a tenth of a per cent, and every kilonewton-metre of moment the arch will ever carry comes from that.

deflection · Rib shortening 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 The answer arrives in instalments. The hogging moment at support 1 of a three-span beam, cycle by cycle. It starts at the fixed-end moment of 98.0 kNm — the value with every joint clamped — and settles at 156.9 kNm against an exact 156.9. The error falls by about a factor of four per cycle: 21.03, 5.92, 1.54, 0.60 kNm after one, two, three and four. Two cycles is an engineering answer and nobody had to invert anything.

Why it converges, and how fast

Moment distribution is an iteration, and iterations do not always converge. This one always does, at a rate the beam's own proportions fix — about a factor of four per cycle on a regular beam and considerably worse on an irregular one, which is where the method's reputation for two cycles being enough comes from and where it stops being true.

deflection · Moment distribution Two quotients from one guessed shape. The critical load of a pin-ended column whose outer quarters keep 10% of the middle's flexural rigidity from four guessed shapes, each worked two ways: Rayleigh's quotient, strain energy of the guess's own curvature over the work of the load, and Timoshenko's, which uses the curvature the guess's bending moment would cause instead. The exact load is 3.225 EI₀/L². a half sine: 8.256 by Rayleigh and 3.745 by Timoshenko; its own sag shape: 8.041 by Rayleigh and 3.712 by Timoshenko; a mid-span sag: 8.875 by Rayleigh and 3.847 by Timoshenko; a parabola: 6.600 by Rayleigh and 3.492 by Timoshenko. Both are upper bounds, and from the same shape the second is never the worse of the two.

The bound from underneath

Rayleigh's quotient turns a guessed shape into a buckling load that is always too high. Divide the same guess differently — use the curvature its bending moment would cause instead of its own — and a parabola that was 21.6 per cent high is 1.3 per cent high. Add one more number that needs no guess at all and the load is caught from below as well, which is the only side of it an amplifier can safely use.

stability · Stability energy

Named alongside it

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

Virtual workCompatibilityForce methodMoment distributionReciprocitySuperpositionArchAxial shorteningBending momentConvergenceDeflectionEigenvalue

All concepts