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 itself. A 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 stiffness and geometric-stiffness matrices that give a tip-loaded column its Euler load. 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.
Fig. 1 A 30 m column carrying nothing but itself, with its buckled shape beside 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: it comes out at qcrL3/EI=7.835q_{cr}L^3/EI = 7.835, and it is equivalent to a tip load of 3.18 times as much total weight. 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 — and the height limit that follows puts a section of 80 mm radius of gyration on the ground at 52 m.

Which free body produced the number

Cut the column at height xx. The piece above it weighs q(L−x)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 curve. Failure 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.
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.467 EI/L22.467\,EI/L^2. Written the same way as the self-weight answer:

qcrL⋅L2EI=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 matters. Two 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.
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 utilisation matters beyond the eigenvalue itself, because it is the denominator of the amplification — a load that makes itself worse is what a mast at 60% of its critical load has, and the wind deflection it was designed for is multiplied by 2.5 along with the moment at its base. So the interaction being linear does not mean the consequences are: a mast that has spent 60% of its capacity on self-weight has spent rather more than 60% of its comfort.

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.837 E i2ρ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 Lmax∝i2/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

That is the same steel in a different shape deciding a height rather than a strength. Every section in that table can be rolled at the same weight per metre and therefore at the same qq; what separates them is the one length a section takes into a column, which is the only property of the cross-section that survives into the answer at all. Area cancels, depth does not enter on its own, and I/A\sqrt{I/A} is the whole of it.

The same column redrawn at the tube’s radius of gyration rather than the 80 mm of the first figure:

A column that has nothing on it but itself. A 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 stiffness and geometric-stiffness matrices that give a tip-loaded column its Euler load. 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 175 mm falls over on its own at 87 m and one of twice that reaches only 138.
Fig. 4 The identical calculation with one number changed. The eigenvalue is still 7.835 and the mode is the same shape, because neither depends on the section at all — but the height at which the column falls over on its own has gone from 52 m to 87 m, and one of twice that radius reaches 138. Nothing about the steel changed; only the distance of its material from its own centroid did.

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.

The claim is easier to believe when the same column is drawn twice with the material swapped out from under it. Aluminium first, at a third of steel’s modulus and a third of its density:

A column that has nothing on it but itself. A 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 stiffness and geometric-stiffness matrices that give a tip-loaded column its Euler load. 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 51 m and one of twice that reaches only 81.
Fig. 5 The same 30 m column in aluminium. The eigenvalue is unchanged at 7.835 — it is a property of the geometry and the load distribution, and no material enters it — and the height limit for an 80 mm radius of gyration comes out at 51 m against steel’s 52. A material a third as stiff and a third as heavy builds the same tower.

Then the harder case, because timber is not merely lighter and softer in the same proportion but nineteen times softer and seventeen times lighter, and there is no reason in advance that those two factors should nearly cancel.

A column that has nothing on it but itself. A 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 stiffness and geometric-stiffness matrices that give a tip-loaded column its Euler load. 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 50 m and one of twice that reaches only 79.
Fig. 6 And in structural timber, at E=11,000E = 11{,}000 N/mm² and a density of 450 kg/m³. The limit is 50 m: three per cent below aluminium’s and four per cent below steel’s, from a material with nothing else in common with either. Three drawings, three materials, one tower — because E/ρE/\rho is 27, 26 and 24 MN·m/kg and nothing else in the expression cares.

The ordinary column problem asks the same question from the other end, and it is worth naming which end. What slenderness costs is a curve of capacity against slenderness; the self-weight limit is that curve read at its far right, 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.

The same exponent argument in its other famous form is span to the fourth: deflection under self-weight goes as L4/EIL^4/EI, which under the same scaling goes 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 is really the question of whether stiffness or strength arrives first, asked about size instead of about a section: strength scales one way, stiffness another and stability a third, so a structure large enough fails by whichever of them scales worst, whatever its 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.

Drawn, it is the same column with the taper switched on:

A column that has nothing on it but itself. A 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 = 6.181, a number with no π in it, computed here as the smallest eigenvalue of the same stiffness and geometric-stiffness matrices that give a tip-loaded column its Euler load. It is equivalent to a tip load of 3.66 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 48 m and one of twice that reaches only 76.
Fig. 7 The bottom row of that table as a picture: 85% of the stiffness removed at the top, linearly. The eigenvalue falls from 7.835 to 6.181 — a fifth — and the tip-load equivalent rises from 3.18 to 3.66, because the taper has moved the weight further toward the base where it does least harm. The mode is visibly more bunched at the base than the prismatic one, which is the mechanism rather than a consequence of it.

This site has met the section that changes along the span in bending, where the question is which section along the taper is the worst one. Here the question 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. What a real column does instead of buckling is bend from the first increment of load, amplified by 1/(1−P/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 (2j−1/3,1/3)2\left(2j_{-1/3,1}/3\right)^2 where j−1/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.

There is a third route to 7.837 that needs neither Bessel functions nor sixty elements, and it is worth knowing because it is the one that can be done on paper. Guess the shape, equate the strain energy stored in that shape to the work the weight does in lowering itself, and the ratio is an upper bound on the eigenvalue; with a reasonable guess it lands within a per cent or two of the exact value, and the guess does not have to be good to be useful.

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.

A scale model cannot show it

There is a consequence for testing that follows from the exponents and that is worth having, because it removes a check somebody might otherwise think they have.

Take a structure and build a geometrically similar model at scale ss. Its second moment goes as s4s^4 and its area as s2s^2, so the group that sets the limit,

Hcr∝(EIρAg)1/3∝(s2)1/3=s2/3H_{cr} \propto \left(\frac{EI}{\rho A g}\right)^{1/3} \propto \left(s^2\right)^{1/3} = s^{2/3}

while the model’s actual height goes as ss. So the ratio the structure cares about scales as

HHcr∝ss2/3=s1/3\frac{H}{H_{cr}} \propto \frac{s}{s^{2/3}} = s^{1/3}

At a scale of 1:100 that is 0.2150.215 — the model sits at about a fifth of the prototype’s fraction of its own limit, which is to say nearly five times further from failing.

So a model that stands up proves nothing whatever about the full-size structure’s self-weight stability. It is the same conclusion the size effect reaches about strength, arriving from a completely different mechanism: a quantity with a length in it does not scale, and here the length is the height at which the structure’s own weight becomes its critical load.

Nor can the discrepancy be corrected by loading the model, because the effect being tested is distributed along the height rather than applied at the top. Adding weight at the tip reproduces the tip-load eigenvalue, not this one, and the two have different mode shapes and different coefficients.

Holding the top changes everything

Every number above is for a free-standing member — Greenhill’s cantilever, fixed at the base and free at the top, with the coefficient 7.837.

Restrain the top and the coefficient rises by a large factor, in exactly the way a column’s effective length does and for the same reason: the mode has to acquire a node, and the shape it is forced into stores far more strain energy for the same tip movement.

That matters most for the family of members the previous section identified as the ones this calculation is actually made about. A falsework tower leg or a temporary prop is nearly always in contact with the structure it is propping, so its top is restrained laterally by the thing it is holding up — and it is not free-standing at all. Treating it as a Greenhill cantilever is conservative by a wide margin.

But the restraint has to be real, and this is where the argument turns. A prop whose head bears on a soffit is restrained against lateral movement only by friction at the bearing, and friction against a horizontal force that is itself proportional to the lean is a restraint of unknown stiffness. A prop with a positive tie at its head is restrained; a prop resting under a slab is not, and the two look identical on site.

Drawing the restrained case beside the free-standing one is the whole of the difference:

A column that has nothing on it but itself. A 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 = 52.488, a number with no π in it, computed here as the smallest eigenvalue of the same stiffness and geometric-stiffness matrices that give a tip-loaded column its Euler load. It is equivalent to a tip load of 2.60 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 97 m and one of twice that reaches only 154.
Fig. 8 The same 30 m column, fixed at the base and now held laterally at the top. The eigenvalue rises from 7.835 to 52.488, a factor of 6.7, and the mode has acquired a node it did not have — which is where the extra strain energy comes from. The height limit for an 80 mm radius of gyration goes from 52 m to 97 m, and the tip-load equivalent falls from 3.18 to 2.60, because the buckle is no longer concentrated where the weight is least.

So the honest position for a temporary member is that it sits somewhere between two coefficients differing by a factor of nearly seven, and which one applies is decided by a detail at the head rather than by anything about the member. The stability of a prop is a property of how it is fixed at the top, and that is the one thing the calculation does not contain.

The same ambiguity runs up a falsework tower. Each lift is braced to the one below, so the legs are restrained at every node — but only against the tower’s own frame, which is itself free to sway unless something anchors it. A tower whose legs are perfectly braced to each other and whose whole frame is free to lean has legs with a short effective length and a frame with a long one, and the eigenvalue that governs belongs to the second. It is the same trap as a row of beams braced only to each other: a bracing system stiff internally and anchored to nothing restrains a mode the structure was not going to use.

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.

What links here

Every essay whose body links to this one.

The objects this essay names

Each one links to every other essay that touches it.

BucklingColumnCritical loadEigenvalueGeometric stiffnessInteractionMastRadius of gyrationScaleSelf-weightSlendernessSpecific stiffnessTaper