Stability

Too tall for nothing but itself

Every critical load on this site so far has been applied at the top of a column. A mast carries a load that is zero at the top and largest at the base, the governing equation stops being harmonic, and the answer comes out as a number with no π in it — along with a maximum height that is almost the same for steel, aluminium and wood.

Assumes Strong enough and still falls over, The ends decide the length that matters and Span to the fourth, which is why spans are short.

Euler’s column has a force on top of it, constant all the way down. That single feature is what makes the governing equation EIv+Pv=0EIv'''' + Pv'' = 0 with constant coefficients, what makes the buckled shape a sine, and what puts π\pi into every answer in the subject.

A mast, a chimney, a lighting column and a tree do not have that. Their load is their own weight: zero at the tip and largest at the base, varying continuously in between. The coefficient in the differential equation varies with height, the solution stops being harmonic, and the critical quantity is no longer a force at all — it is a load per unit length.

qcrL3EI=7.837\frac{q_{cr}L^3}{EI} = 7.837

with no π\pi anywhere in it.

A column that has nothing on it but itselfA 30 m column carrying no load except its own weight, with its buckled shape and the axial force that produced it. The force is zero at the top and largest at the base, which is why the answer is not Euler's: the eigenvalue is a load intensity and comes out as q_cr L³/EI = 7.835, a number with no π in it, computed here as the smallest eigenvalue of the same matrices every other column on this site uses. It is equivalent to a tip load of 3.18 times as much total weight — a column carries its own weight better than it carries somebody else's, because most of the weight is near the base where the buckle is not. The height limit that follows is a cube root, so a section of radius of gyration 80 mm falls over on its own at 52 m and one of twice that reaches only 82.the buckleits own weightlargest at the basei = 40 mm32 mi = 60 mm43 mi = 80 mm52 mi = 120 mm68 mi = 200 mm95 mthe height it falls over atq L³/EI = 7.835a cube root — twice the shape buys 59%
Fig. 1 The column, its buckled shape, and the axial force that produced it. The force is largest at the base and nothing at the top, which is why the mode is bunched toward the base and why the eigenvalue is not Euler’s. The number is measured here rather than quoted: it is the smallest eigenvalue of the same matrices every other column on this site is built from.

Which free body produced the number

Cut the column at height xx. The piece above it weighs q(Lx)q(L-x), so that is the axial force at the cut — and it is a function of xx rather than a constant.

Feed that into the geometric stiffness. Each element of the column contributes a KgK_g proportional to the axial force it carries, so the destabilising matrix is built from a varying force rather than a uniform one, and the eigenvalue problem

(KλKg)v=0\left(K - \lambda K_g\right)\mathbf{v} = 0

is the same problem with a different second matrix. Solved on sixty elements it gives λL3/EI=7.8365\lambda L^3/EI = 7.8365, against the value 7.837 that Greenhill obtained in 1881 by recognising the differential equation as Bessel’s of order one-third.

The same machinery, handed a constant axial force instead, reproduces π2EI/L2\pi^2EI/L^2, π2EI/4L2\pi^2EI/4L^2 and 20.19EI/L220.19EI/L^2 to seven figures. That is what makes 7.8365 worth quoting rather than merely computed.

The column curveFailure load against slenderness, as a fraction of the squash load. A stocky column crushes; a slender one buckles at the Euler load; the crossover is where the two curves meet, and real columns fall below both near it.5010015020000.20.40.60.811.2slenderness (effective length ÷ radius of gyration)they cross at λ = 75squashingEuler bucklingreal columns, which are neither
Fig. 2 The constant-force case, for comparison. Everything about the mechanism is identical and everything about the answer differs, because a coefficient that was constant has become a linear function of height.

A column carries its own weight better than somebody else’s

Put the total weight qLqL at the tip instead and the column buckles at π2EI/4L2\pi^2EI/4L^2, which for the same column is 2.467EI/L22.467\,EI/L^2. Written the same way as the self-weight answer:

qcrLL2EI=7.837againstPcrL2EI=2.467q_{cr}L \cdot \frac{L^2}{EI} = 7.837 \qquad\text{against}\qquad P_{cr}\frac{L^2}{EI} = 2.467

so the column will carry 3.18 times as much total weight if the weight is spread along it as it will if the same weight is put on top.

The reason is where the buckle is. A cantilever’s first mode has its largest curvature at the base and its largest displacement at the tip; the destabilising work a load does is its force times the shortening of the column beneath it, which is largest for a load at the top. Self-weight puts most of its force near the base, where the shortening beneath is small, so most of it is doing very little harm.

The load’s distribution matters as much as its size, which is a statement no column curve or effective length can accommodate.

The ends decide the length that mattersTwo columns of identical height and section, buckling under two sets of end conditions. The effective length factor is the fraction of the column that behaves like a pin-ended one, and the buckling load goes as its inverse square.K = 0.5both ends fixedK = 0.7one fixed, one pinnedsame column, same section, two ways of holding the endsthe load at which each buckles goes as 1 ÷ K²
Fig. 3 The effective length device exists to stretch Euler’s answer to more end conditions. It cannot stretch to this: an effective length is a length at which a constant force gives the same critical load, and there is no constant force here to have one.

The two together, and the interaction is a straight line

Real masts have both — their own weight and something on top, an aerial, a tank, a lamp. Sweeping the reference load from all-tip to all-self-weight and reading the eigenvalue at each:

PPcr+qL(qL)cr=1.00  to  1.01\frac{P}{P_{cr}} + \frac{qL}{(qL)_{cr}} = 1.00 \;\text{to}\; 1.01

across the whole range, with a maximum departure of 1.4%.

That is a linear interaction, and linear interactions are rare enough in stability to be worth noticing. Two loads that buckle a member in the same mode and differ only in where they are applied add their utilisations, and nothing more complicated is needed.

The practical form: a lighting column that is at 60% of its self-weight buckling load has 40% of its tip capacity left, and no correction factor is required to say so.

The load that makes itself worseThe amplification of a deflection against the ratio of applied load to buckling load. A structure at half its buckling load deflects twice as far as first-order analysis predicts, and the curve runs away well before the load is reached.00.20.40.60.80246810applied load ÷ buckling load1.3×1.7×2.5×first-order analysis says the answer is always 1×one over one minus the ratio
Fig. 4 The amplification that follows from any of it. The critical load is the denominator; a mast at 60% of it has its wind deflection multiplied by 2.5, and the moment at its base multiplied with it.

The height limit, and the cube root that flattens everything

Write q=ρgAq = \rho g A and I=Ai2I = A i^2, where ii is the radius of gyration. The area cancels:

Lmax=(7.837Ei2ρg)1/3L_{max} = \left(\frac{7.837\, E\, i^2}{\rho g}\right)^{1/3}

Three things are worth reading off that expression.

Strength does not appear. Not the yield stress, not the ultimate stress, not any material property except the modulus and the density. A column that falls over under its own weight has not been overstressed anywhere; it has run out of stiffness.

It is a cube root, so nothing moves it much. Quadrupling EIEI buys 59%. Halving the density buys 26%. A material twice as stiff and half as dense — which does not exist — buys 59% again.

And the shape matters more than the material. ii enters squared inside a cube root, so Lmaxi2/3L_{max} \propto i^{2/3}: doubling the radius of gyration buys 22/3=1.5872^{2/3} = 1.587, exactly. That is more than any material substitution available, and it is free — a hollow tube of the same area has several times the radius of gyration of a solid bar.

radius of gyration height limit, steel
40 mm 32.5 m
80 mm 51.5 m
120 mm 67.5 m
175 mm (a 508 × 12.5 tube) 86.8 m
400 mm 150.7 m
The same material, three waysThree cross-sections of identical area, so identical weight and cost, with the second moment of area computed from each profile's own geometry. Only the arrangement differs, and the stiffest is many times the flattest.tall rectangleI = 64.86 × 10⁶1.0× the firstsquare hollowI = 98.17 × 10⁶1.5× the firstI-sectionI = 66.77 × 10⁶1.0× the firstevery section here has an area of 19458 — only the shape differsthe bar is the second moment of area, to scale
Fig. 5 The same steel in a different shape, which here decides a height rather than a strength. Every one of these sections has the same weight per metre and the same qq; their radii of gyration differ by a factor of three, and their height limits by a factor of two.
The one length a section carries into a columnFour profiles of equal area, with the radius of gyration r = √(I/A) drawn as the distance it is — a pair of lines either side of the centroid, at the depth the whole area would have to sit at to give the section the second moment it has. As a 4 m pin-ended column the same 19458 mm² of material carries between 1988 and 8649 kN, in the ratio of the squares of those radii and of nothing else.the same, laid flatr = 28.1 mmλ = 1421988 kNsquarer = 40.3 mmλ = 994087 kNtall rectangler = 57.7 mmλ = 698402 kNI-sectionr = 58.6 mmλ = 688649 kNthe dashed pair is ±r about the centroidthe bar is the Euler load at 4 m, to scale
Fig. 6 And the one length a section takes into a column, which is the only property of the cross-section that appears in the answer at all. Area cancels, depth does not enter, and I/A\sqrt{I/A} is the whole of it.

Three materials, one tower

Now the result that makes the cube root worth having. Put a radius of gyration of 80 mm into the expression for four materials:

EE (N/mm²) ρ\rho (kg/m³) E/ρE/\rho (MN·m/kg) LmaxL_{max}
steel 210,000 7,850 27 51.5 m
aluminium 70,000 2,700 26 51.0 m
structural timber 11,000 450 24 50.0 m
concrete 30,000 2,400 13 40.0 m

The first three agree to within three per cent, across a factor of nineteen in modulus and a factor of seventeen in density.

They agree because the answer depends on E/ρE/\rho alone — the specific stiffness — and the specific stiffness of nearly every structural material in common use is within a small factor of 25 MN·m/kg. Steel is not stiffer than wood in the sense that matters here; it is stiffer and heavier in almost exactly the same proportion.

That is not a coincidence about materials; it is a fact about atoms, and it belongs to physics rather than to this site. What belongs here is the consequence: a self-weight height limit is not a material choice. It is a shape choice, and the only material that changes it much is concrete, whose specific stiffness is half everybody else’s.

Length costs more than it looksThe same column section at four lengths, with the buckling capacity of each drawn as a bar. Capacity falls as the inverse square of the length, so a column three times as long carries a ninth as much.1× the length100% of the capacity1.5× the length44% of the capacity2× the length25% of the capacity3× the length11% of the capacityidentical section, identical material, identical end conditions
Fig. 7 What slenderness costs in the ordinary column problem. The self-weight limit is the same graph read at its far right-hand end, where the load has fallen to what the member weighs and the question is whether there is anything left over.

The scaling argument underneath it

There is a shorter way to reach the same conclusion, and it explains why the cube root had to be there.

Scale a column by a factor α\alpha in every dimension. Its weight per unit length goes as α2\alpha^2; its bending stiffness EIEI goes as α4\alpha^4; its length goes as α\alpha. So the group qL3/EIqL^3/EI goes as α2α3/α4=α\alpha^2 \cdot \alpha^3/\alpha^4 = \alpha.

The dimensionless measure of how close a column is to falling over under its own weight grows linearly with the scale. Double every dimension of a tower and it is twice as close to its own limit; there is no arrangement of material that avoids it, because the argument used no material property at all.

That is the same square-cube reasoning that runs through everything on this site about size, and it arrives here in its purest form: no stress, no strength, no factor of safety, just an exponent.

Deflection goes as the fourth power of the spanDeflection against span for a constant load intensity and section, with two slower relationships drawn faintly behind it for comparison: the load itself, which grows in proportion to the span, and the bending moment, which grows as its square. Doubling the span multiplies the deflection by sixteen, while the moment only quadruples.11.522.533.54050100150200250300span, relative to the first16×81×moment: the squareload: the first powerdeflection: the fourth
Fig. 8 The span dependence in its other famous form. Deflection under self-weight goes as L4/EIL^4/EI over the same scaling and therefore as α3\alpha^3 — so a scaled-up structure gets closer to buckling linearly and floppier cubically, and the second usually stops the design first.
Which limit arrives firstUtilisation of the strength limit and of the deflection limit, against span. Strength grows as the square of the span and deflection as the fourth power, so the two cross — and past the crossing a beam is sized by how far it moves rather than by what it can carry.0.60.811.21.41.61.8200.511.5span, relative to the firstthey cross heredeflection runs out at 1.40strength runs out at 1.41the limitstrengthdeflection
Fig. 9 Which failure arrives first, which is the question the scaling argument really answers. Strength scales one way with size, stiffness another, and stability a third; a structure large enough will fail by the one that scales worst, whatever the designer intended.

Taper, and why every mast has one

A real mast is not prismatic. It tapers, because the base carries the whole weight and the tip carries nothing, and a section sized for the base is wasted everywhere else.

Taper cuts both ways for stability. It removes stiffness from the top of the column, where the buckling mode has its largest displacement; but it also removes weight from the top, where the weight is doing most of the destabilising. Computed on the same column with the stiffness falling linearly to a fraction of its base value:

stiffness lost at the top qcrL3/EIbaseq_{cr}L^3/EI_{base}
none 7.836
30% 7.319
60% 6.741
85% 6.182

A column that has given up 85% of its tip stiffness has given up 21% of its buckling load. That is a very good trade, and it is why masts, chimneys and trees are all tapered — the second-moment penalty is small because the top of the column contributes little bending resistance to a mode whose curvature is at the base.

Two curves climbing together, and the one that catches up firstA 6 m member tapering from 200 to 600 mm, with the moment it carries and the moment it can carry drawn on the same scale below it. The demand rises linearly and the capacity as the square of the depth, so the gap between them closes and then opens again. It is narrowest at 3.00 m from the free end, where the member is 400 mm deep and 79% used, against 70% at the root where the moment is largest.30 kNwhat it can carrywhat it carriesworst at 3.00 mfree endroot79% used at the governing station · 70% at the root
Fig. 10 The tapered member, which this site has already met in bending. There the question is where the worst section is; here it is what the taper does to a mode shape, and the two answers are unrelated except in sharing a geometry.

Where the model stops

The column is elastic and perfectly straight. At the slendernesses this matters at — a 50 m steel tube of 175 mm radius of gyration has a slenderness of nearly 300 — it genuinely is elastic, and it is emphatically not straight. A real mast has an out-of-straightness of the order of L/500L/500, so it never bifurcates: it bends from the first increment of self-weight, and the eigenvalue is an asymptote rather than an event.

There is no wind. Every real mast is designed for a lateral load rather than for this, and the self-weight buckling load enters the calculation as the denominator in the amplification factor instead of as a capacity in its own right.

And the material is linear. For a tall timber pole it is not, and the modulus that matters for a load sustained for decades is the long-term one, which is perhaps 60% of the short-term value — worth 15% of the height limit.

What the pictures cannot show

The mode is drawn at an amplitude that means nothing. An eigenvector has no scale, so what the figure gives is the shape the column would take and no information whatever about how far it goes.

Nor can it show the growth. A real column approaches the critical load by bending more and more, and the relationship between load and deflection is a hyperbola that never reaches the asymptote — so a photograph of a mast just below its critical load and one at half of it look very much the same, and the difference is a factor of two in the moment at the base.

A column that was never straightLoad against lateral deflection at mid-height, for a column starting with an initial bow of 0.002. There is no critical value to reach: the deflection grows from the first increment, slowly at first and then without bound as the ratio approaches 1.00 — which is the Euler load, and which the column therefore never attains. The perfect column, drawn for comparison, sits on the vertical axis until it arrives there and then has no answer at all.00.0050.010.0150.020.02500.20.40.60.81lateral deflection at mid-heightload ÷ P꜀ᵣP ÷ P꜀ᵣ = 1.00, approached and never reachedinitial bow: δ₀ = 0.002
Fig. 11 What a real column does instead of buckling: a deflection that grows from nothing, amplified by 1/(1P/Pcr)1/(1-P/P_{cr}), with the eigenvalue visible only as the place the curve is heading.

The history, and the equation that was already solved

Greenhill gave the answer in 1881, in a paper about how tall a tree or a mast could be. He recognised that the governing equation, with an axial force linear in the coordinate, becomes Bessel’s equation of order one-third under a change of variable — and Bessel’s equation had been solved decades earlier for an entirely different reason.

The number 7.837 is (2j1/3,1/3)2\left(2j_{-1/3,1}/3\right)^2 where j1/3,1j_{-1/3,1} is the first zero of a Bessel function of order minus one-third. It is not a rounded value or a fitted one; it is exact, and the seven-thousandths agreement with the finite-element answer above is the finite element being checked rather than the constant.

What is worth carrying from that is the habit rather than the result. A varying coefficient turns a harmonic equation into a special-function one, and the special functions are almost always already tabulated — so a problem that looks intractable is often a problem in a different notation.

Four guesses at one buckling modeA 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.a half sine9.870 EI/L²exactits own sag shape9.882 EI/L²0.13% higha mid-span sag10.000 EI/L²1.32% higha parabola12.000 EI/L²21.59% highreference9.8696 EI/L²ten Ritz terms,as an eigenvalue problemevery guess is anupper boundP
Fig. 12 The energy route to the same number, which needs no Bessel functions at all. Assume a shape, equate the strain energy to the work the weight does in lowering itself, and get an upper bound — and with a reasonable guess it lands within a per cent or two of 7.837.

The assumption the figure rests on

The column is fixed at its base and free at the top, which is the mast’s condition and gives 7.837. Change it and the number changes a great deal: pinned at both ends gives 18.57, and fixed at the base with the top held laterally gives 52.5. Those are the same three-to-twenty spread that the end conditions produce for a tip-loaded column, and a self-weight problem is no less sensitive to them.

Where the limit actually gets close

It is fair to ask whether any of this ever governs, and the answer is that it governs in exactly one family of structures: the very slender ones that carry nothing.

A building is nowhere near it. A forty-storey core has a radius of gyration of several metres, so its self-weight height limit is measured in kilometres; what stops a building getting taller is drift, and after that lifts.

A chimney is closer. A 60 m reinforced concrete stack of 3 m diameter has a radius of gyration around 1.1 m and a limit near 250 m, so it is at a quarter of it — enough that the second-order amplification of its wind moment is a real term but not enough to threaten it.

A lighting column is closest of anything built. An 12 m aluminium column of 60 mm radius of gyration has a limit around 43 m: a factor of 3.6, which after imperfections and the amplification of the wind case is not the comfortable margin it sounds.

And a scaffold standard, a temporary prop, a falsework tower leg — those are the members where a designer meets this calculation for real, and they meet it because temporary works are the one place a structure is deliberately built as slender as it can possibly be.

The ladder from here

Later rungs on this anchor: the Bessel solution itself, and why an axial force linear in the coordinate produces an equation of order one-third. The tapered column solved as a variational problem rather than by finite elements. Trees, whose height limit is set by water transport rather than by buckling, and which sit at a comfortable factor below this one anyway. The self-weight buckling of a shell — a chimney or a silo — where the mode is a local ripple rather than a global lean. Guyed masts, where the guys are springs along the height and the eigenvalue becomes a function of their stiffness. And the same question for a hanging column, where the weight puts the member in tension and there is no eigenvalue at all.

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.

BucklingColumnCritical loadEigenvalueGeometric stiffnessInteractionMastRadius of gyrationScaleSelf weightSlendernessSpecific stiffnessTaper