Stability

The buckling load with no compression in it

Twist a straight bar hard enough and it snaps into a helix, with no axial load on it anywhere. The load at which that happens is 2πEI/L — first power of the length, and no shear modulus in it at all, so how stiff the bar is in torsion has nothing to do with the torque that buckles it in torsion.

Assumes Strong enough and still falls over, The internal force with no diagram and The stiffness the load takes away.

Every buckling load in this collection so far has had a compression in it somewhere. Euler’s column is a compression along the member. A beam that fails sideways has a compression flange. A plate that ripples is in compression across its width, and a ring that goes out of round is under a hoop compression. The mechanism is always the same: a force that was doing something useful, once the member has moved, starts doing something else.

There is a case where nothing is in compression at all and the same thing happens.

A buckling load with no compression anywhere in the member. A tube of radius 81.6 mm and wall 5 mm under pure torque — no axial load at all — buckles into a helix at T = 2πEI/L, which for the section drawn is 1877 kNm. Two things are absent from that expression and both are surprising: the length appears to the first power rather than the second, and the shear modulus does not appear at all, so the torsional stiffness of the member has nothing to do with the torque at which it buckles in torsion. The comparison is the torque that yields the same tube, 42.9 kNm — a factor of 43.8 below it. The mode arrives first only past a length of 263 m, which is πE/τ_y = 3219 radii and contains no thickness whatever.
Fig. 1 A tube under pure torque — no axial load anywhere on it — buckling into a helix at T = 2πEI/L, which for a 168 × 5 section 6 m long is 1,832 kNm. Below it, the same expression against length, with the torque that yields the same tube for comparison. Two things are absent from 2πEI/L and both are surprising: the length appears to the first power, and the shear modulus does not appear at all.

Why a torque can buckle something

Take the bar straight, and apply a torque TT about its own axis at each end. The torque vector points along the bar. Nothing about that configuration is unstable, and nothing in the bar is in compression.

Now let the bar bow, by a very small amount, into a shape y(z)y(z). The ends have rotated, because the ends of a bowed bar are no longer perpendicular to the original axis, and the torque vector — which is applied about a fixed direction, or about the bar’s own end tangent, depending on the mechanism applying it — is no longer parallel to the deflected bar. It now has a component transverse to the bar.

A transverse component of a torque is a bending moment about an axis at right angles to the plane of the bow, and it therefore bends the bar in the perpendicular plane. That bow acquires its own transverse torque component, which bends the bar back in the first plane. The two bending planes drive each other, and the deflected shape they sustain between them is a helix.

Writing the two coupled equations,

EIy=Tz,EIz=Ty,EI\,y'''' = -T\,z''', \qquad EI\,z'''' = T\,y''',

and combining them in the complex variable u=y+izu = y + iz gives EIu=iTuEI\,u'''' = iT\,u''', whose non-trivial solution for a bar pinned at both ends against bending but free to rotate about its axis is

Tcr=2πEIL.T_{cr} = \frac{2\pi EI}{L}.

Greenhill wrote it in 1883, the year after he did the same for a column under its own weight — the problem this site has already drawn — and the second result is much less known than the first.

What is not in the expression

There is no GG and no JJ. The critical torque of a member in torsion contains only its bending stiffness. That is worth stopping on, because it is entirely counter-intuitive and it is exactly right: what buckles is a bending mode, and the torque is merely the agent driving it. The bar’s resistance to being twisted is irrelevant, because it is not being asked to twist any more than it already is.

One slit, and the torsional stiffness falls by a factor of hundreds. A 200 by 200 box of 8 mm wall, drawn closed and then slit along its length. Closed, the torque runs round the wall as a shear flow and the torsion constant is 5.66×10⁷ mm⁴; slit, the loop is broken and only each wall's own thickness resists, giving 1.31×10⁵ mm⁴. The ratio is 432 to one, so the same torque twists the slit section 432 times as far and raises a peak shear stress 36 times as high. Nothing about the material changed.
Fig. 2 The consequence, put beside the factor of several hundred that closing a section is worth in torsional stiffness. A tube and the same tube slit along its length differ by two orders of magnitude in J, and by nothing whatever in the torque at which they buckle. Every intuition about open and closed sections in torsion is about a different question from this one.

The length appears to the first power. Euler’s load goes as 1/L21/L^2 and this goes as 1/L1/L, which means the two do not scale together: the ratio Tcr/(PEL)T_{cr}/(P_E L) is a constant of the section, and a longer member loses axial capacity faster than it loses torsional capacity. That is the opposite of what the shared word “buckling” suggests.

And the constant is 2π2\pi. No square roots, no numerical fitting, no Poisson’s ratio. The end conditions change it in a very simple way — 2 for a bar pinned at both ends, 4 for one clamped at both, and 1.4303 for a cantilever with a torque at its free end.

The ends decide the length that matters. Four columns of identical height and section, buckling under four 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 familiar version of the same statement for a column, where the four end conditions multiply the critical load by 4, 1, 0.5 and 0.25. Under torque the same four conditions multiply it by 2, 1, and 0.72 relative to the pinned case — a much narrower spread, because a helix is a mode that cares less about its ends than a half-sine does.

Which free body produced the number

The free body is a length of the bar cut at one end, with the applied torque on the cut face and the bar deflected. Two things are on it: the torque, drawn as a vector along the original axis, and the internal bending moments on the cut. The whole instability is contained in the observation that those two are not parallel once the bar has moved, and the component of the first about the second’s axis is what the equations balance.

That is the same structure of argument as second-order analysis in general: a force whose direction relative to the member changes as the member moves. What is different here is that the force is a couple rather than a push, and that the equilibrium being disturbed is one in which no member was under any threat.

Two frequencies that meet, and a determinant that never moves. The two natural frequencies of Ziegler's two-bar column against the follower load Pℓ/k, with the determinant of its stiffness matrix drawn along the top. The determinant is k² at every load — it varies over this whole axis by 1.1e-16 of itself, which is round-off — so a static buckling analysis of this structure finds no critical load whatever and reports it as stable everywhere. The frequencies say otherwise: they approach, meet at Pℓ/k = 2.0858 — the closed form is (7 − 2√2)/2 = 2.0858 — and become a complex pair, which is oscillation that grows. Adding any internal damping at all drops the load at which that happens to 1.4643, which is 41/28 and 30% below the undamped value; the limit of the damped system is not the undamped system, which is the paradox Ziegler found in 1952 and which was taken for an arithmetic error for a decade.
Fig. 4 The other member of the same family — a load that turns to follow the member it is destabilising, which cannot be treated by an eigenvalue at all. Torsional buckling sits between the two: the answer depends on how the torque is applied, and the three end conditions above are three different answers to what is the torque attached to.

The mode it is not

There is a second thing called torsional buckling in this subject and it is a different phenomenon with an unfortunately similar name.

A channel has three critical loads, not one. The three critical loads of a channel in compression, against its length, with the load it actually buckles at drawn over them. At 3000 mm the flexural loads are 19014 kN about the major axis and 3139 kN about the minor, while twisting about the shear centre takes 1962 kN. The lowest root is 1879 kN, and the column twists. The shear centre sits 109.7 mm from the centroid, so the modes cannot happen separately: the lowest root of the coupled problem is 4.2 per cent below the lowest of the three, and the Wagner coefficient β is 0.60. The governing mode changes at 5813 mm: below that length the column twists, above it, it bends — because the torsional resistance keeps a term that does not grow when the member is shortened, and the flexural loads have none.
Fig. 5 A column that twists instead of bending, which is the other object the words describe. There the load is an axial compression, the mode is a rotation of the cross-section about the shear centre, and the critical load contains GJ and the warping constant — everything the result on this page does not have. The two are unrelated except in vocabulary.

Telling them apart is easy once the question is what is applied. In flexural-torsional buckling an axial force finds a twisting mode. Here a torque finds a bending mode. The first is a property of the section’s shape — it is why channels and angles are awkward and hollow sections are not — and the second is a property of nothing but EIEI and LL.

The number that ends it

A critical load is only a design case if it arrives before something else. Compare the critical torque with the torque that simply yields the same section, Ty=τyJ/rT_y = \tau_y J/r:

TcrTy=2πEπr3t/Lτy2πr3t/r=πErLτy.\frac{T_{cr}}{T_y} = \frac{2\pi E \pi r^3 t / L}{\tau_y \cdot 2\pi r^3 t / r} = \frac{\pi E r}{L\,\tau_y}.

Every dimension has cancelled except the radius, and the thickness has gone entirely. Setting the ratio to one gives the length at which the mode becomes reachable:

L=πEτyr3,142rL^* = \frac{\pi E}{\tau_y}\,r \approx 3{,}142\,r

for any structural steel, since E/τyE/\tau_y is about a thousand for all of them. The 168 mm tube drawn would have to be 256 m long before it buckled under torque rather than yielding — and at 6 m the factor is 42.7.

Length costs more than it looks. The 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.
Fig. 6 For scale: an axially loaded member of the same section reaches its Euler load at a length of about 25 radii. The two failure modes that share the word buckling are separated by a factor of more than a hundred in the slenderness required to reach them, and the reason is the power of the length.

So this is a real instability, exactly computed, that governs nothing anybody would build. It matters in exactly one place, and there the length is not a problem: a drill string is a tube several kilometres long, driven by torque from the surface, and helical buckling is its central design difficulty.

The half that does matter

The result that a structural engineer can use is not the critical torque but its effect on something else.

The load that makes itself worse. The 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.
Fig. 7 Because a torque and an axial load destabilise the same bending mode, they compete for the same capacity. A member carrying half its critical torque has lost part of its axial capacity even though nothing about its compression has changed, and the interaction is very nearly a circle in the two ratios.

(TTcr)2+PPE=1\left(\frac{T}{T_{cr}}\right)^2 + \frac{P}{P_E} = 1

is the linear-theory result for pinned ends, and it is worth reading for what it says at small TT: because the torque term is squared, a torque of a third of the critical value costs about 11 per cent of the axial capacity and one of a tenth costs 1 per cent. The interaction is flat where it matters, which is the reason nobody checks it.

That flatness is the real conclusion of this essay. The mode exists, the expression is exact and elegant, the coupling is genuine — and the practical instruction it produces is ignore it, arrived at by computing it rather than by never having asked.

There is one place the flatness stops helping, and it is where a member is already close to its axial capacity. At P/PE=0.9P/P_E = 0.9 the remaining torque capacity is 0.1=0.32\sqrt{0.1} = 0.32 of TcrT_{cr} — so a nearly-buckled column has lost two thirds of a capacity it never knew it had. Whether that matters depends on whether anything is applying a torque, which in a building frame it very rarely is.

A column that was never straight. Load 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.
Fig. 8 And the reason the interaction is drawn as a bifurcation curve rather than as a design curve. Every real member has an imperfection, so nothing ever reaches a bifurcation load of any kind; what a real member does is amplify what it started with. A torque-and-compression interaction curve is a boundary in a plane no real member’s path reaches — it is a statement about where the amplification becomes unbounded, not about where anything fails.

What it shares with the rest of the subject

Strip away the torque and the helix, and this result belongs to the same family as everything else in this field.

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. 9 The column curve, which is the shape every one of these results takes once yield is admitted. Below a certain slenderness the material decides and the elastic critical load is irrelevant; above it the elastic load decides and the material is irrelevant; and the interesting members are in the transition. Torsional buckling has the same two branches — the yield torque and 2πEI/L2\pi EI/L — and the transition sits at 3,142 radii instead of about 90.

The general statement is that a load applied to a member does two things at once: it produces stresses, and it removes stiffness. The removal is the geometric stiffness, a term proportional to the force already present, and a buckling load is the load at which it has removed all of the structural stiffness there was. Which mode it removes the stiffness from depends on how the load moves when the member does, and a torque moves in a way that attacks bending.

That framing also explains the absent GG in one line. Geometric stiffness is subtracted from whatever stiffness resists the mode being driven; the mode here is bending; so the stiffness that appears is the bending one. The torsional stiffness resists a mode nobody is driving.

The drill string, which is where it is not academic

Every structure in this collection is short in the sense that matters here. A drill string is not: it is a continuous tube two to eight kilometres long, driven from the surface by a torque large enough to turn a bit through rock, and carrying almost no axial compression because most of its weight is held up in tension from above.

Put 3,000 m into L=3,142rL^* = 3{,}142\,r and the radius at which the mode arrives is 955 mm — vastly larger than a drill pipe, which is about 60 mm. So a drill string is not merely past the threshold; it is past it by more than an order of magnitude, and helical buckling is not a check but the normal operating state of the lower part of the string.

Two things change once the mode is not merely reachable but unavoidable. The helix has a borehole wall to lie against, so the deflection is bounded and the problem becomes one of contact rather than of stability — what is wanted is the pitch of the helix and the friction it generates, not a critical load. And the friction the helix generates against the wall absorbs torque that was meant for the bit, which is why the torque at the bit is a fraction of the torque at the surface and why the fraction is not known.

Neither of those questions is a structural-statics question, and both are downstream of an expression written in 1883 about a bar with nothing pushing on it.

Where the model stops

The bar is prismatic, straight and elastic. Every one of those is doing work. A bar with an initial bow has no bifurcation at all; the helix grows from the first newton-metre, and the amplification is one over one minus the load ratio, exactly as an imperfect column’s deflection does.

The torque is applied in a stated way. “Semi-tangential”, “quasi-tangential” and “axial” torques are three different loadings with three different critical values, differing by up to a factor of two — because a torque is a couple, and a couple applied through a universal joint is not the same load as one applied through a gear. That sensitivity is the same one a follower force has, and it means the answer depends on the machine at the end of the member.

The section is doubly symmetric. For anything else the bending and twisting are coupled before the torque arrives, and the problem becomes a general flexural-torsional one with the applied torque as an extra term.

The comparison with yield uses a thin tube. For an open section the elastic torsional stress at a given torque is far higher, so TyT_y falls and the reachable length falls with it. An open section still buckles at 2πEI/L2\pi EI/L and yields very much sooner, so the mode is further out of reach rather than nearer.

And the drawing is of the mode, not of the path. The helix has an undetermined amplitude, as every eigenvector does, and no stress can be computed from it. What happens past the critical torque — whether the helix grows stably or the bar wraps itself into a lock-up — is a large-displacement question this linear result cannot be asked.

The ladder from here

Later rungs on this anchor: the drill string proper, where the tube is confined inside a borehole so the helix has a wall to lie against, and the buckling is a contact problem whose answer is a pitch rather than a load. Coiled tubing, which is the same problem with residual curvature in the member before it starts. The combined torque, axial load and internal pressure case, which is the real drilling condition. Torsional buckling of a bar with an initial helical imperfection, and the amplification factor it obeys. Lateral-torsional buckling read as a member of this family, where the destabilising agent is a bending moment rather than a torque and the mode is the same coupled bending-and-twisting. And the historical thread: Greenhill’s two 1883 papers between them opened the whole subject of non-conservative and geometrically-coupled instability, and everything on this site about second-order effects is downstream of them.

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.

BucklingCritical loadEigenvalueGeometric stiffnessInteractionSlendernessTorsionWhich failure arrives first