Dynamics

The motion that feeds itself

A steady wind contains no frequency at all, and it can destroy a bridge. The force that does it is manufactured by the structure's own movement, so there is no excitation to resonate with — there is a wind speed above which the equilibrium is unstable, and below which nothing happens.

Assumes The wind that brings its own frequency and The bridge that was pushed by its own sway.

An iced power line in a steady winter wind swings through metres, slowly, for hours. A rectangular signal gantry in a steady breeze oscillates until its welds crack. A suspension bridge deck in a steady 19 m/s wind twists itself apart.

None of those is resonance, and the reason is the word steady. A steady wind contains no frequency for a structure to resonate with. The oscillation has to be manufactured, and the structure manufactures it.

The damping a wind leaves behind, and the speed that uses it upTotal damping ratio against wind speed, for a structure with 0.60% of its own and a modal mass of 25 kg per metre at 1.2 Hz. The aerodynamic damping of a section whose lift falls with angle of attack is negative and grows with speed, so at 4.14 m/s what is left is nothing. Below the threshold a disturbance dies away. Above it, the structure feeds itself and the motion grows out of nothing at all — which is what makes this a different mechanism from resonance rather than a severe case of it.024681012-0.01-0.0050.005wind speed (m/s)total damping ratio0.60% damping · 1.2 Hz4.14 m/s — nothing leftbelow: a disturbance dies awayabove: the structure drives itself
Fig. 1 The total damping of a bluff section against wind speed. The structure’s own 0.6% is used up at 4.1 m/s, and above that the total is negative: any disturbance grows, with nothing oscillating at any frequency anywhere in the problem. The horizontal axis is a speed, and the curve crosses zero rather than peaking.

Den Hartog’s criterion

Consider a section moving across a steady wind at velocity y˙\dot y. Relative to the section, the wind is no longer horizontal: it arrives at an angle αy˙/U\alpha \approx \dot y / U, because the section’s own motion has added a component.

If the section’s lift changes with angle of attack, that induced angle produces a change in lift. And if the lift falls as the angle rises — which is what happens on a bluff section past a certain angle, and never on a well-behaved aerofoil — the change in force is in the same direction as the motion.

Writing the aerodynamic force as a damping term gives the equation of motion

my¨+(c+12ρUD[dCLdα+CD])y˙+ky=0m\ddot y + \left(c + \tfrac{1}{2}\rho U D \left[\frac{dC_L}{d\alpha} + C_D\right]\right)\dot y + ky = 0

and the bracket is Den Hartog’s criterion. When

dCLdα+CD<0\frac{dC_L}{d\alpha} + C_D < 0

the aerodynamic damping is negative, and since it grows with UU there is a speed at which it overcomes the structure’s own:

Ucrit=4mζωρDdCL/dα+CDU_{\text{crit}} = \frac{4 m \zeta \omega}{\rho D \left|dC_L/d\alpha + C_D\right|}

Three things about that expression deserve attention, and each contradicts an intuition.

The criterion is about the SHAPE. Whether a section can gallop at all is decided by two aerodynamic coefficients and nothing else. A circular cylinder cannot — it has no angle of attack, so dCL/dαdC_L/d\alpha is zero and the criterion is satisfied by the drag alone. A square section can, a D-section can, and an iced conductor can, because the ice makes an aerofoil of an object that was a circle.

The frequency does not appear except through the damping. There is no coincidence to avoid, no band to stay out of, and nothing to be tuned away from. This is what most distinguishes it from every other problem in this field.

Mass and damping appear as a product, exactly as in the Scruton number, and for the same reason: what matters is how much the structure can dissipate against how much the air can supply.

The damping a wind leaves behind, and the speed that uses it upTotal damping ratio against wind speed, for a structure with 0.60% of its own and a modal mass of 200 kg per metre at 1.2 Hz. The aerodynamic damping of a section whose lift falls with angle of attack is negative and grows with speed, so at 33.09 m/s what is left is nothing. Below the threshold a disturbance dies away. Above it, the structure feeds itself and the motion grows out of nothing at all — which is what makes this a different mechanism from resonance rather than a severe case of it.05101520253035400.0020.0040.006wind speed (m/s)total damping ratio0.60% damping · 1.2 Hz33.09 m/s — nothing leftbelow: a disturbance dies awayabove: the structure drives itself
Fig. 2 The same section eight times heavier per metre. The critical speed goes from 4.1 m/s to 33.1 — eight times, in exact proportion — and moves from a breeze that occurs weekly to a wind that occurs in a serious storm. Mass buys critical speed linearly, which is the strongest lever available and the most expensive.
The damping a wind leaves behind, and the speed that uses it upTotal damping ratio against wind speed, for a structure with 0.60% of its own and a modal mass of 25 kg per metre at 1.2 Hz. The aerodynamic damping of a section whose lift falls with angle of attack is negative and grows with speed, so at 0 m/s what is left is nothing. Below the threshold a disturbance dies away. Above it, the structure feeds itself and the motion grows out of nothing at all — which is what makes this a different mechanism from resonance rather than a severe case of it.05101520253035400.010.020.030.04wind speed (m/s)total damping ratio0.60% damping · 1.2 Hzstable at every value drawnbelow: a disturbance dies awayabove: the structure drives itself
Fig. 3 And the same structure with a section whose lift RISES with angle of attack. The criterion is positive, the aerodynamic damping is positive at every speed, and the curve never crosses zero: this section does not gallop at any wind speed whatever. Shape is the only variable in this problem that can remove the instability rather than postponing it.

What growth without a driver looks like

Either side of the threshold, from the same 2 mm nudgeThe same deck, displaced 2 mm and released, at 3.33 and 4.76 m/s — one below the threshold of 4.14 and one above. The total damping ratios are 0.117% and -0.090%. After a minute the first has fallen to 1.18 mm and the second has reached 3 mm and is still growing. Nothing is forcing either of them.0102030405060-3-2-1123time (s)displacement (mm)3.33 m/s — total damping 0.117%4.76 m/s — total damping -0.090%
Fig. 4 Either side of the threshold, from the same 2 mm disturbance, with a steady wind and no oscillating force anywhere in the integration. Below the critical speed the motion dies as any structure’s does. Above it the same disturbance grows and keeps growing. Nothing is being applied at any frequency; the structure is extracting energy from a wind that is not varying at all.

That figure is the whole distinction, and it is worth stating as a pair of tests a reader can apply to any reported failure.

A resonant response requires an excitation at a matching frequency; is bounded by damping at 1/2ζ1/2\zeta; grows to a steady amplitude and stays; and disappears when the excitation moves off the frequency.

A self-excited response requires only that a parameter exceed a threshold; is not bounded by anything linear; grows exponentially from any disturbance; and cannot be removed by detuning.

Pulled to 2 mm and let go, on a structure of 0.833 s periodDisplacement against time for a single-degree-of-freedom structure of natural period 0.833 s and -0.6% damping, under pulled to 2 mm and let go. The static deflection under the same peak force is 35.15 mm and the peak response is 9.03 mm — a factor of 0.26.051015202530-10-5510time (s)displacement (mm)pulled to 2 mm and let gonothing takes energy outpeak 9.03 mm at 33 s
Fig. 5 The growth on its own, as a free vibration with the damping ratio negative. Forty cycles from a 2 mm disturbance and the amplitude has multiplied several times over, with the equation containing no forcing term at all. Every self-excited problem in this field reduces to this picture; only the mechanism supplying the negative sign differs.
Pulled to 2 mm and let go, on a structure of 0.833 s periodDisplacement against time for a single-degree-of-freedom structure of natural period 0.833 s and 0.6% damping, under pulled to 2 mm and let go. The static deflection under the same peak force is 35.15 mm and the peak response is 2 mm — a factor of 0.06.051015202530-2-112time (s)displacement (mm)pulled to 2 mm and let goeach cycle is 0.963 of the one beforepeak 2 mm at 0.000 s
Fig. 6 And the same structure at the same damping magnitude with the sign the right way round — a wind just below the critical speed. The difference between the two figures is one sign in one coefficient, and it is the difference between a structure that is fine and a structure that is not.

Tacoma Narrows, and why the textbooks are wrong

The bridge that collapsed on 7 November 1940 is the most famous illustration of resonance in engineering, and it was not resonance.

The wind that day was around 19 m/s and steady. The motion that destroyed the deck was torsional — the deck twisting about its own axis, one edge up and the other down — and it developed after the bridge had been oscillating vertically for some hours, following a change in the deck’s behaviour when a mid-span tie slipped. The torsional motion grew over about an hour to angles beyond 30 degrees.

The mechanism was aeroelastic: the deck’s own rotation changed the flow around it in a way that added energy to the rotation. It is usually classified as torsional flutter, which is the same family as galloping with an angular coordinate rather than a linear one, and complicated by the deck’s H-shaped section shedding vortices at the same time. The deck was stiffened by shallow plate girders rather than by a truss, on a main span of 853 m, giving a depth-to-span ratio of about one in three hundred and fifty and a torsional stiffness lower than any comparable bridge then standing — which is the property every subsequent long-span deck has been designed to have more of.

What it certainly was not is a vortex-shedding resonance. The shedding frequency of the deck at that wind speed was around 1 Hz, and the torsional mode was at 0.2 Hz. The numbers do not match, and they were checked at the time.

The reason the wrong explanation survived is instructive, and it is not carelessness. The film shows a structure oscillating at a steady frequency with growing amplitude, and that is exactly what resonance looks like. The two mechanisms produce visually identical evidence, and the only way to distinguish them is by asking what the excitation was — which requires knowing that a steady wind has no frequency in it, and that a growing oscillation therefore has to be self-excited.

This site’s rule about naming the free body has a version for this field: name the excitation, and say what its frequency is. If the answer is “the structure’s own”, the problem is not resonance.

What it took to establish that

The correction to the Tacoma story was not a matter of anyone rereading the film. It came from the wind tunnel work that followed the collapse and from the aeroelastic theory that the aircraft industry had been building since the 1920s, when flutter began destroying aeroplanes for exactly the same reason.

That transfer is worth noticing. The mechanism that destroyed a bridge had been understood, analysed and designed against by another profession for fifteen years, in a field that shared no journals, no conferences and no vocabulary with structural engineering. The bridge engineers’ apparatus for wind was a static pressure; the aeronautical engineers’ was a set of motion-dependent forces measured on a model.

It is the strongest argument this collection can offer for the figure-first habit it is built on: the two professions were drawing different free bodies. One drew a deck with a pressure on it; the other drew a section with forces that depended on how the section was moving. Everything else followed from that choice, and the choice was invisible to both of them until a bridge fell down.

What flutter needs that galloping does not

Galloping is the one-degree-of-freedom member of the family and the only one with a criterion this simple. Classical flutter — of an aircraft wing, and of a bridge deck — needs two degrees of freedom that couple.

The mechanism is that vertical motion and rotation feed each other: rotating changes the lift, which drives the vertical motion, which changes the angle of attack, which drives the rotation. Below the critical speed the two motions are out of phase in a way that dissipates energy and above it they are in phase in a way that absorbs it, and the transition is a speed at which two of the system’s modes coalesce.

Two consequences worth carrying:

Flutter needs the two frequencies to be close, which makes it the one member of this family for which a mode list is the first thing to look at. A deck whose torsional frequency is well above its vertical one is much less susceptible, which is why modern long-span decks are torsionally stiff box girders rather than open trusses or plate girders, and why the ratio of the two frequencies is a headline number in a long-span design.

It cannot be found by any calculation that treats the wind as a load. The critical speed comes out of an eigenvalue problem in which the aerodynamic forces are functions of the motion — the flutter derivatives, measured in a wind tunnel on a section model. That measurement is a routine part of long-span bridge design and exists because of 1940.

The shedding frequency against wind speed, and where it stops following itVortices leave a 1.2 m cylinder at St·U/D, which is the faint straight line, so there is always a wind speed at which they match the structure's own 0.9 Hz — here 6 m/s, which is a breeze. Inside a band from 5.1 to 7.8 m/s the shedding abandons the line and locks on to the structure, driven by the structure's own motion.024681012141618200123wind speed (m/s)shedding frequency (Hz)1.2 m across · St 0.18 · 0.9 Hzlock-in6 m/s
Fig. 7 The third member of the family, for comparison: vortex lock-in, where the structure’s motion captures the shedding frequency over a band of wind speeds. It is the mild case — the amplitude is limited by the flow’s own nonlinearity and the band is finite — and it shares the essential feature, which is that the excitation is not independent of the response.
How much a harmonic force is magnified, at three damping ratiosDisplacement amplitude divided by the deflection the same force would produce if it were applied slowly, against the ratio of the forcing frequency to the structure's own, at 1%, 2%, 5% of critical damping. At the natural frequency the magnification is 83.33, 25, 10 respectively — one over twice the damping ratio, and nothing else in the problem enters it.00.511.522.50510152025forcing frequency ÷ natural frequencyamplitude ÷ static deflection1% damping — 83.33× at the peak2% damping — 25.01× at the peak5% damping — 10.01× at the peak
Fig. 8 For comparison, what a resonant response looks like when it is bounded: a peak, at a particular frequency ratio, of a height set by the damping. Every self-excited problem in this essay has no such curve at all. There is no frequency axis, because there is no forcing frequency — the picture is a threshold in a speed, and past it the amplitude is not a number the linear theory can supply.

The comparison is worth making explicit because it is diagnostic. If a structure’s problem can be drawn as a peak on a magnification curve, damping bounds it and a tuned mass can reduce it. If it can only be drawn as a crossing, then damping does not bound the response — it sets the speed at which the response begins, which is a different thing to buy and a different thing to check.

The family, and why it belongs with buckling

Three problems, in three industries, with three names: a crowd on a footbridge, a wind on a bluff section, and a flow over a coupled pair of modes. In each of them the force is a function of the motion; in each there is a threshold in a parameter; and in each, below the threshold nothing happens at all.

That is not the structure of a load problem. It is the structure of a stability problem, and the resemblance to buckling is exact rather than poetic:

buckling self-excitation
what is lost stiffness damping
driving parameter axial load wind speed, or crowd size
below the threshold any disturbance decays any disturbance decays
above it any disturbance grows any disturbance grows
size of the disturbance irrelevant irrelevant
what a linear analysis gives the critical value the critical value
what bounds the answer nonlinearity and imperfections nonlinearity

A column that is strong enough and falls over anyway and a deck that is strong enough and tears itself apart anyway are the same failure of the same kind of thinking: a check on the response to a given load says nothing whatever about whether the equilibrium the response is measured from is stable.

What the picture cannot show

The coefficients are measurements. dCL/dαdC_L/d\alpha and CDC_D come from wind tunnel tests on the section, they vary with Reynolds number and with turbulence, and for an iced conductor they vary with the shape of the ice. The criterion is exact and its inputs are not.

The growth stops somewhere. The linear model grows without bound. In practice galloping settles into a large-amplitude limit cycle, because the aerodynamic coefficients change once the amplitude is big enough that the induced angle is no longer small. The structure oscillates hugely and steadily rather than accelerating to destruction — which is why iced conductors gallop for hours and are usually still there afterwards.

One mode is assumed. The criterion above is written for a section translating in one direction. A real member gallops in whichever mode is least stable, and for a long cable or conductor that is a mode with several half-waves along the span — which changes the effective mass and the damping, and therefore the critical speed, without changing the criterion.

Turbulence matters and is not in the model. Atmospheric turbulence disrupts the coherent flow that galloping needs, so a section in rough flow is more stable than in smooth. The dangerous conditions are again steady ones — which is the third time in this field that the worst case has turned out to be a mild, steady, unremarkable day rather than a storm.

Where the ladder goes

This is the last rung of the field’s foundation and it points at the whole of it. The essays before it built up a way of asking what a structure does about a load; this one is where that question stops being the right one.

The immediate direction is the flutter analysis proper — the two-degree-of-freedom eigenvalue problem with measured aerodynamic derivatives, which is how every long-span bridge deck is now checked and how the aircraft industry has worked since the 1930s.

The wider direction is the classification. This field opened with a load that arrives suddenly and a factor of two, and it ends with a load that does not arrive at all until the structure invites it. Between the two lie every intermediate case: a load with a frequency, a load with a spectrum, a load that is a displacement of the ground, and a load supplied by the people using the structure. What organises them is not what causes the load but how the load is related to the response — independent of it, statistically related to it, or manufactured by it — and that is the ordering worth carrying out of this field into any other.

Read that way, the field has a shape. The first four essays are about a structure meeting a load it does not affect. The middle five are about loads described statistically, where the structure and the excitation are still independent but the answer is a distribution. And the last three are about loads that would not exist if the structure held still — which is the point at which the free body has to include the thing applying the force, and the point at which dynamics stops being statics with an extra term.

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.

AeroelasticityCritical velocityDampingFlutterGallopingNegative dampingSelf excitationStability threshold