Dynamics

The period nobody chose

Every structure has a natural period, it decides the answer to every question in this field, and no drawing anywhere records it. It is a consequence of a mass picked for one reason and a stiffness picked for another — and it has already been computed, by the serviceability check.

Assumes The deflection that arrives three years late and Twice the deflection, for the same load.

The last essay’s answer depended on one number about the structure and nothing else: how long it takes to swing once. A quarter of a second was nearly instantaneous for that platform, and it would have been leisurely for a bracket and geological for a suspension bridge.

That number is not on any drawing. Nobody specified it, no client asked for it, and it does not appear in any calculation done for strength. It is a consequence — of a mass fixed by what the building is made of and what goes in it, and a stiffness fixed by a deflection limit that was checked for an entirely different reason.

Rayleigh's method: the frequency read off the deflection that was computed anywayA simply supported beam of 8 m, sagging 13.08 mm under its own weight, sampled at 21 stations. Rayleigh's quotient over those deflections gives 4.91 Hz against the exact 4.91 Hz — 0.119% high, and high rather than low because an assumed shape is a constraint and a constraint stiffens. The rule of thumb, 18 divided by the square root of the deflection in millimetres, gives 4.98 Hz.8 m simple span · sag 13.08 mmthe static shape, not an eigenvectorω² = g · Σ(w·δ) ⁄ Σ(w·δ²)Σ(w·δ) = 197.3 · Σ(w·δ²) = 2.030Rayleigh, from the static shape: 4.915 Hzexact, from the characteristic equation: 4.909 Hz18/√δ with δ in mm: 4.977 Hz
Fig. 1 An 8 m floor beam sagging 13.08 mm under its own weight — a calculation done for the span/360 check and then filed. Rayleigh’s quotient over the same deflections gives 4.915 Hz. Solving the beam’s characteristic equation properly gives 4.909 Hz. The difference is 0.12%, in the safe direction for a reason that is a theorem rather than luck.

The frequency was in the deflection all along. It takes one line of arithmetic to get it out, and the line is the same for a beam, a floor, a frame or a mast.

Why mass over stiffness, and nothing else

Draw the free body once more: a mass, a spring, and no applied load at all. The spring pulls back with kuku and the mass resists acceleration with mu¨m\ddot u, so

mu¨+ku=0m\ddot u + ku = 0

and the only motion that satisfies it is a sine wave at ω=k/m\omega = \sqrt{k/m}. Two properties of the structure go in and one number comes out; nothing else can enter, because nothing else is in the equation.

That is a stronger statement than it looks. The amplitude is not in it. A structure swings at the same rate whether it is nudged or shoved, which is why a period can be measured at all — hit the thing however hard is convenient and time the swings. It is also why the whole of this field can be organised around a single number per structure, in a way that no strength calculation ever can.

The square root is the part that misleads. Frequency responds to stiffness and mass weakly, and both directions have been measured on the same beam:

A beam's frequency against its span, which falls as one over the squareThe fundamental frequency of a simply supported beam of fixed section and fixed mass per metre, against its span, from 2.5 to 20 m. At 8 m it is 4.91 Hz; At 16 m it is 1.23 Hz. Doubling the span quarters the frequency, because the frequency goes as βL squared over the span squared and the flexibility it is competing with goes as the fourth power.46810121416182001020304050span (m)fundamental frequency (Hz)8 m — 4.91 Hz16 m — 1.23 Hz
Fig. 2 The same section and the same mass per metre, at spans from 2.5 to 20 m. The frequency falls as one over the square of the span: 19.6 Hz at 4 m, 4.91 Hz at 8 m, 1.23 Hz at 16 m. Doubling a span quarters its frequency — the strongest lever in this field, and it is geometry rather than material.

Doubling the stiffness of a beam raises its frequency by a factor of 1.4142, which is √2 and is the entire return on doubling the material working in bending. Adding half again to the mass drops it by 0.8165, which is 1/√1.5. Doubling the span, by contrast, quarters the frequency, because span enters both terms at once — the stiffness falls as the cube of it and the mass rises in proportion to it.

So the honest summary of how to change a structure’s period is: change the span. The rest is square roots.

Rayleigh’s method, which is free

The exact frequency of a beam comes from a characteristic equation with hyperbolic functions in it. The approximate one comes from a quantity every serviceability check already produces.

Rayleigh’s argument is an energy statement and it takes three lines. Assume the structure vibrates in some shape; at the extremes of the motion all the energy is strain energy, and at the middle all of it is kinetic. Equate the two, and the frequency drops out:

ω2=gwδwδ2\omega^2 = \frac{g\sum w\delta}{\sum w\delta^2}

with ww the weight at each point and δ\delta its deflection under those weights. Every quantity in it is in the static analysis.

Two things about this are worth knowing, and both are provable rather than empirical.

It is always an over-estimate. Assuming a shape is imposing a constraint, a constraint stiffens, and a stiffer structure has a higher frequency. So Rayleigh’s answer bounds the true one from above, always — 4.915 against 4.909 here, and never the other way round. A method that is guaranteed to err in one direction is worth far more than one that is more accurate on average.

It is astonishingly close. The static deflected shape under self-weight is not the first mode, but it is very nearly the first mode, because the loading that produces it is distributed the same way the inertia forces are. A tenth of a per cent is typical, and the error is second-order in the difference between the assumed shape and the real one — which is the general property of a Rayleigh quotient and the reason it was worth a name.

The rule of thumb every engineer carries falls straight out of it. Take a simply supported beam under self-weight, put the closed-form deflection into the quotient, and the result is

f18δf \approx \frac{18}{\sqrt{\delta}}

with δ\delta the mid-span deflection in millimetres. On the beam above that gives 4.977 Hz against 4.909 — 1.4% high, and high for the same reason. A floor that sags 13 mm is a 5 Hz floor, and that sentence is all most problems need, since what a floor does about footfall turns on which side of 5 Hz it lands.

The ends matter more than anything in the middle

The characteristic equation is where the end conditions enter, and they enter hard.

What the ends are worth: four support conditions at one spanThe fundamental frequency of the same beam — 8 m, the same section, the same mass per metre — under four support conditions. simply supported: 4.91 Hz; cantilever: 1.75 Hz; fixed at both ends: 11.13 Hz; fixed one end, pinned the other: 7.67 Hz. The range is a factor of 6.36, and nothing about the beam itself changed.8 m span, identical sectionthe ends alone4.91 Hz1.00× the simple beamsimply supported1.75 Hz0.36× the simple beamcantilever11.13 Hz2.27× the simple beamfixed at both ends7.67 Hz1.56× the simple beamfixed one end, pinned the other
Fig. 3 One beam, one span, one section, four sets of end conditions. The cantilever is at 1.75 Hz and the beam fixed at both ends at 11.13 Hz — a factor of 6.36, with nothing about the beam itself changed. Fixing both ends of a simply supported beam multiplies its frequency by 2.27.

That factor of 6.36 is the same phenomenon as the factor of sixteen between a column’s end conditions, and it arises the same way: the boundary conditions decide the shape, the shape decides the curvature, and the curvature is what the stiffness acts on. In both cases the number that changes is a root of a transcendental equation, and in both cases the practical consequence is that the least well-known part of a structure is the part that decides the answer.

Because the frequency goes as the square of that root, small differences in end fixity are magnified. A beam whose ends are neither pinned nor fixed — which is every beam, as the connections field spends its time establishing — has a frequency somewhere between 4.91 and 11.13 Hz, and a measurement on the finished floor is worth more than any calculation.

Higher modes exist, and are not evenly spaced

The characteristic equation has infinitely many roots, and each is a different way the beam can vibrate.

The first three modes of a simply supported beamThree modes of a simply supported beam of 8 m span, drawn from the general solution with the constants fixed by the support conditions rather than assumed to be sines. Mode 1 is at 4.91 Hz with βL = 3.1416; Mode 2 is at 19.63 Hz with βL = 6.2832; Mode 3 is at 44.18 Hz with βL = 9.4248. The frequencies go as the square of βL, so the threeth mode is 9.0 times the first. The marked points are the nodes.Simply supported, 8 m spaneach shape scaled to its own peakmode 1: 4.91 HzβL = 3.1416mode 2: 19.63 HzβL = 6.2832mode 3: 44.18 HzβL = 9.4248
Fig. 4 The first three modes of the simply supported beam, drawn from the general solution rather than assumed to be sines. The frequencies are 4.91, 19.6 and 44.2 Hz — in the ratio 1 : 4 : 9, because βL is nπ and the frequency goes as its square. The marked points are the nodes: the second mode has one, the third has two, and a beam struck at a node cannot excite that mode at all.
The first three modes of a cantilever beamThree modes of a cantilever beam of 8 m span, drawn from the general solution with the constants fixed by the support conditions rather than assumed to be sines. Mode 1 is at 1.75 Hz with βL = 1.8751; Mode 2 is at 10.96 Hz with βL = 4.6941; Mode 3 is at 30.69 Hz with βL = 7.8548. The frequencies go as the square of βL, so the threeth mode is 17.5 times the first. The marked points are the nodes.Cantilever, 8 m spaneach shape scaled to its own peakmode 1: 1.75 HzβL = 1.8751mode 2: 10.96 HzβL = 4.6941mode 3: 30.69 HzβL = 7.8548
Fig. 5 The same three modes for a cantilever, where the spacing is not 1 : 4 : 9. The roots are 1.875, 4.694 and 7.855, so the ratios are 1 : 6.27 : 17.55. The even spacing of the simple beam’s overtones is a property of that particular set of end conditions and of nothing else — which is the whole difference between a structure and a musical instrument.

The ratio 1 : 4 : 9 is worth pausing on because it is the reason a simply supported beam is a poor musical analogy. A string’s overtones are 1 : 2 : 3, which is why a string sounds like a note; a beam’s are 1 : 4 : 9, which is why a struck beam sounds like a clank. The difference is that a string’s restoring force comes from tension and a beam’s from bending, so one is governed by a second derivative and the other by a fourth.

A frame is not a beam, and the period is still a division

For a building the mass sits at the floors and the stiffness is in the columns between them, and the same square root applies at a different scale.

Three modes of a five-storey frameThe first three mode shapes of a five-storey shear frame, from the eigenvalue problem rather than sketched. Mode 1 has a period of 0.6 s, no node and carries 88.0% of the mass; Mode 2 has a period of 0.21 s, one node and carries 8.7% of the mass; Mode 3 has a period of 0.13 s, two nodes and carries 2.4% of the mass. The nth mode crosses the axis n−1 times, which is a theorem rather than a drawing convention.five floors · 300 t eachmode 10.6 s · 1.65 Hz88.0% of the massno nodemode 20.21 s · 4.83 Hz8.7% of the massone nodemode 30.13 s · 7.61 Hz2.4% of the masstwo nodes
Fig. 6 A five-storey frame with 300 tonnes at each floor and 400 MN/m of storey stiffness. Its fundamental period is 0.605 s. The second and third modes are at 0.207 s and 0.130 s, and each has a shape of its own — which is the subject of the next rung and not of this one.

Across a range of heights that model reproduces the rule of thumb every seismic engineer uses. Three storeys gives 0.387 s, five gives 0.605, eight gives 0.932, twelve gives 1.370, twenty gives 2.246 — which is 0.11 to 0.13 seconds per storey, against the standard estimate of about a tenth of a second per storey. The rule is not folklore; it is what m/k\sqrt{m/k} produces when mass and stiffness both scale with height in the ways that buildings make them.

It also shows where the rule stops.

Three modes of a five-storey frameThe first three mode shapes of a five-storey shear frame whose ground storey is 35% as stiff as the others, from the eigenvalue problem rather than sketched. Mode 1 has a period of 0.79 s, no node and carries 96.3% of the mass; Mode 2 has a period of 0.24 s, one node and carries 3.3% of the mass; Mode 3 has a period of 0.14 s, two nodes and carries 0.4% of the mass. The nth mode crosses the axis n−1 times, which is a theorem rather than a drawing convention.five floors · soft ground storeymode 10.79 s · 1.27 Hz96.3% of the massno nodemode 20.24 s · 4.15 Hz3.3% of the massone nodemode 30.14 s · 7.08 Hz0.4% of the masstwo nodes
Fig. 7 The same frame with its ground storey softened to 35% of the others — an open frontage, a car park, a hotel lobby. The fundamental period goes from 0.605 s to 0.801 s, a third longer for a change to one storey out of five, and the first mode’s shape changes with it: the ground-floor drift, 0.285 of the roof’s before, is 0.570 now. Half the building’s deformation has moved into one storey.

That is the soft-storey failure written as an eigenvalue rather than as a photograph. Nothing in the frame got weaker — every column is the one that was there before — and the deformation the earthquake will demand has been concentrated into the storey with the fewest columns to share it.

Which mass is moving

There is a quantity hidden in Rayleigh’s quotient that is worth taking out and looking at, because it answers a question that otherwise sounds unanswerable: when a beam vibrates, how much of it is moving?

Not all of it, and not at the same speed. The mid-span of a simply supported beam travels the full amplitude and the ends do not move at all, so the kinetic energy is less than it would be if the whole mass moved together. Do the integral for a half-sine shape and the answer comes out at half the total mass — that is the mass which, placed at mid-span on a spring of the beam’s stiffness, would have the same period as the beam.

For a cantilever in its first mode the number is about a quarter of the total, because most of the beam is near the root and hardly moves.

This is not bookkeeping. It is the quantity that decides whether a person walking on a floor can shift it at all, and it is why a long-span floor is worse than a short one for reasons beyond its lower frequency: the walker’s own weight is a larger fraction of the mass that is actually participating. A structure’s dynamic size is not its size, and the ratio between the two is a property of the mode rather than of the object.

Rayleigh, Dunkerley, and the shafts that came first

The method is Lord Rayleigh’s, from The Theory of Sound in 1877, and it belongs to a period when the subject was acoustics rather than structures — the questions being asked were about bells, plates and organ pipes, and the structural applications came afterwards.

The engineering pressure came from rotating machinery. A turbine shaft has a speed at which it whirls violently, and it is the shaft’s own bending frequency; running through it on the way to operating speed is a design problem with an expensive failure mode. Dunkerley’s rule, published in the 1890s, gave a way to estimate the critical speed of a shaft carrying several rotors by adding the reciprocals of the squares of the frequencies each rotor would give on its own — an approximation that errs low, where Rayleigh’s errs high, so the pair of them bracket the answer.

The two rules bracketing from opposite sides is the useful part, and it is a pattern this site meets in other fields: the plastic collapse load is bracketed by an upper-bound mechanism and a lower-bound stress field, and both are worth more together than either is alone.

What is striking, from here, is that this apparatus was fully developed before it was needed for buildings at all. Structures were not analysed dynamically until the twentieth century — not because the mathematics was missing but because the loads that need it were rare when buildings were heavy and stiff, and because nobody had a way to measure the one number the theory cannot supply.

What the picture cannot show, and what the number assumes

The figures above give a frequency to three decimal places for a beam whose real frequency nobody knows to two.

The mass is a guess about occupancy. A floor’s period depends on what is standing on it, and the load that will be there when the vibration matters is not the load the strength calculation used. Design codes ask for a fraction of the imposed load in the vibrating mass, and the fraction is a judgement.

The stiffness is a guess about cracking. A concrete floor’s second moment of area depends on how cracked it is, which depends on how it has been loaded since it was cast. A factor of two between the uncracked and fully cracked values is ordinary, and it moves the frequency by 1.4.

The end conditions are the largest uncertainty and are drawn as though they were exact. The figure above shows four discrete cases. The real beam is between two of them, and the answer moves by a factor of two across that gap.

So the honest reading is that a computed natural frequency is a number with a band around it of perhaps ±20%, and the design decisions that turn on it should survive that band. This is exactly why a floor’s response is presented as a curve against frequency rather than as a single number: what a designer needs to know is not the frequency but whether the neighbourhood of it is a bad place to be.

The decay measures what nothing else can

One property of a structure genuinely cannot be calculated from the drawing, and the same test that measures the period gives it away.

Pulled to 20 mm and let go, on a structure of 0.500 s periodDisplacement against time for a single-degree-of-freedom structure of natural period 0.500 s and 2.0% damping, under pulled to 20 mm and let go. The static deflection under the same peak force is 12.67 mm and the peak response is 20 mm — a factor of 1.58.024681012-20-101020time (s)displacement (mm)pulled to 20 mm and let goeach cycle is 0.882 of the one beforepeak 20 mm at 0.000 s
Fig. 8 The same structure, pulled 20 mm out of position and released. The period is read off the spacing of the peaks; the damping is read off how fast they shrink — each cycle here is 0.882 of the one before. Everything about the first quantity was predictable from a drawing. Nothing about the second was.

Which is the subject of the next essay, and the reason it is a separate one: the period comes from the structure’s geometry and its material, and the damping comes from neither.

Where the ladder goes

Three directions lead out of here, and each is an essay of its own.

The first is that a real structure has many periods and not one, and the question of which of them matters is the question of where the mass is — the next rung.

The second is what happens when something drives a structure at the period it has, which is where the number computed here turns into a number of millimetres.

The third is the practical one, and it is the reason engineers compute this at all: the period decides which of the loads a structure meets are static and which are not. A 5 Hz floor is dynamic under footfall and static under wind. A 0.2 Hz tower is dynamic under wind and static under footfall. Same two loads, same physics, opposite conclusions — and the only thing that decides which is which is a number nobody wrote down.

What this makes readable

Essays that name this one as a prerequisite.

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.

Characteristic equationDeflectionModal massMode shapeNatural frequencyNatural periodRayleigh methodStiffness