Structural Engineering & Statics
where the load goes
Generator

Four guesses at one buckling mode

Rendered here at the parameters it defaults to, with every essay that calls it — which is the same list as the blast radius of changing it.
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.

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.

5 essays call rayleigh-strut. The drawing above is what it returns with no arguments at all; every call below passes it something, because a placement that passes nothing draws whichever member of the family the generator happens to default to rather than the one its essay argues about.

Where it is called

Changing this generator changes every one of these figures.

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. Stability

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.

 The check that everything adds up, and the error it cannot see. Four versions of the same 2-bay, 3-storey frame, with the global equilibrium residual each one produces — the sum of the reactions against the sum of the applied loads, as a fraction of the applied total. It is the first thing every analysis prints and it is worth having: a lost restraint and a load entered in the wrong unit both show up immediately, at 6% and 24%, because both change what the structure is carrying. The fourth bar is the point. A member whose stiffness is wrong by a factor of ten redistributes the internal forces completely — the second bar shows the change in the member forces, 18% — and the global residual is exactly zero, because the wrong answer is still in equilibrium with the same loads. Equilibrium is one equation per degree of freedom of the whole body, and a stiffness error lives entirely in the many equations underneath it. A model can satisfy every equilibrium check ever devised and be a model of a different structure. Stability

The stiffness the load takes away

Buckling is usually taught as an event — a critical load, a bifurcation, a mode. Written as a matrix it stops being an event at all. A compressive load subtracts a stiffness from the structure, the subtraction grows with the load, and the critical load is simply where what is left reaches zero.

 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. Deflection

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.

 The average is not the answer, and it is unsafe. Critical load of a pinned column whose middle third has been given a different stiffness, against the whole-column Euler load, with the two numbers a hand check reaches for beside it. The eigenvalue is taken from K − P·Kg over 24 elements, so nothing here is a formula for a stepped column — it is the same computation the uniform case gets. At a middle third of 0.50 times the rest the true load is 0.612 of Euler's, the arithmetic average says 0.832 and the weakest segment says 0.496. The average is high by 36% and it is high on the unsafe side, because the third of the column it is averaging over is the third where the mode has all its curvature. The weakest-segment answer is safe everywhere and wasteful by about as much. Stability

An average stiffness is not a safe stiffness

Euler's load belongs to a column of one EI. Give the same column two, and the temptation is to average them — which is wrong, and wrong in the unsafe direction by a quarter. Buckling weights stiffness by the square of the curvature of the mode, so the middle of a pinned column decides everything and the ends decide almost nothing.

 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. Stability

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.

The library, page 4 of 7 — where rayleigh-strut sits