Dynamics

The liquid has a period of its own

Shake a tank and its contents do not all go with it. Part of the liquid moves as though it were bolted to the wall and part sloshes at a period fixed by gravity and the radius, which the tank's stiffness has no influence over whatever. The split is decided by one proportion, and the two parts then take entirely different amounts of the earthquake.

Assumes Most of the mass moves together, The spectrum is not a load and The force that is only a radius.

Every dynamics figure on this site so far has put the mass inside the structure. A column has mass, a floor plate has mass, and when the ground moves they move with it because they are attached.

A liquid is not attached to anything. It is contained, which is a different relationship, and the difference is the whole of this essay. Shake a tank and the liquid near the bottom goes along as though it were solid, while the liquid near the free surface does something else entirely — it sloshes, at a period that has nothing to do with the tank’s stiffness and everything to do with gravity and the shape of the surface.

A broad tank sloshes and a tall one does notThe liquid's division into the part that moves with the wall and the part that sloshes, against the tank's proportion. The convective masses come from the potential-flow solution and the impulsive mass is whatever is left, so the two sum to the liquid's mass exactly at every proportion rather than approximately over part of the range. A tall tank at H/R = 3 is 84% impulsive and behaves almost like a solid; a shallow one at H/R = 0.5 is 72% convective and most of its contents never notice the earthquake. This tank sits at H/R = 1.33, which is 65% impulsive.0123400.20.40.60.81liquid depth ÷ radiusfraction of the liquidwith the wallslosheslever ÷ Hthis tank: H/R = 1.33, 65% impulsive
Fig. 1 The liquid’s division into the part that moves with the wall and the part that sloshes, against the tank’s proportion. The two sum to the liquid’s mass exactly at every proportion.

Which free body produced the number

The split comes from the potential-flow solution for a rigid circular tank, and it is worth saying which half of it is computed and which half is a remainder.

The convective modes are the sloshing modes of the free surface, and each has a mass, a height and a frequency that fall out of the flow solution:

mc,nm=2tanh(λnH/R)(λnH/R)(λn21),ωn2=λngRtanh ⁣(λnHR)\frac{m_{c,n}}{m} = \frac{2\tanh(\lambda_n H/R)}{(\lambda_n H/R)(\lambda_n^2 - 1)}, \qquad \omega_n^2 = \frac{\lambda_n g}{R}\tanh\!\left(\frac{\lambda_n H}{R}\right)

with λn\lambda_n the roots of a Bessel derivative. The impulsive mass is then whatever is left, which is what makes the accounting exact: the parts sum to the whole by construction rather than by an approximation that stops being true for a tall tank.

The height the impulsive mass acts at is not tabulated either; it comes from a statics argument. Accelerate a rigid tank of rigid liquid sideways at aa and the pressure on the wall is ρax\rho a x — independent of depth — so the resultant is mama acting at H/2H/2. The modal masses and heights therefore have to reproduce that:

mihi+mc,nhc,n=mH2m_i h_i + \sum m_{c,n} h_{c,n} = m\,\frac{H}{2}

and hih_i falls out. It is the one place in this field where a free body of the whole body settles a modal quantity, and it is worth having because it makes the decomposition checkable: a set of masses and heights that does not reproduce a resultant at mid-height is wrong, whatever it was derived from.

Tall tanks and broad ones are different objects

The split depends on one number, H/RH/R, and it depends on it strongly.

A tall tank at H/R=3H/R = 3 is 84% impulsive: most of its contents are too deep to know there is a surface, and the tank behaves almost as though it were carrying a solid. A broad shallow one at H/R=0.5H/R = 0.5 is 72% convective: most of its contents slosh, and most of them therefore never notice the earthquake at all.

That is a genuinely useful design fact and it runs the opposite way to intuition. A broad reservoir attracts far less seismic force than its weight suggests, and a tall slender tank attracts nearly all of it — the same total mass, twice the base shear, decided by a proportion.

The proportion enters through a hyperbolic tangent, which is the signature of a decay with depth. Below about one radius from the surface the sloshing motion has effectively died away, so the convective mass is roughly whatever liquid lies within a radius of the top and the impulsive mass is the rest. That reading makes the whole figure predictable without any of the algebra: a tank deeper than about two radii has a fixed convective mass and an impulsive one that grows with every further metre of depth, which is exactly the shape the curves take.

Two periods, orders apart

The reason the split matters is what the two masses meet when they arrive at a spectrum.

Two periods on one spectrum, and one of them is off the endThe design spectrum with the tank's two periods marked. The impulsive mass — 65% of the liquid, moving with the wall — has a period of 0.250 s and sits on the plateau at 7.36 m/s². The first convective mode has a period of 4.47 s, fixed by gravity and the radius and by nothing structural, and sits far out on the falling branch at 0.37 — a factor of 20 lower. So 35% of the liquid supplies only 3% of the base shear. The sloshing mass is nearly weightless as a force and decides the freeboard, which is a displacement.01234567802468period (s)spectral acceleration (m/s²)impulsive 65%7.36 m/s²convective0.37 m/s²35% of the liquid supplies 3% of the shear
Fig. 2 The two periods on one design spectrum. The impulsive mass sits on the plateau at 7.36 m/s² and the convective one far out on the falling branch at 0.37.

The impulsive period is a structural quantity — the tank wall’s flexibility with the impulsive mass on it — and it is a quarter of a second, which is on the plateau of any ordinary design spectrum. The convective period is 2π/ω12\pi/\omega_1, gravity and radius only, and it is 4.47 seconds, which is well past the corner where a spectrum has begun to fall as 1/T21/T^2.

The two accelerations therefore differ by a factor of twenty, and 35% of the liquid supplies 3% of the base shear.

That is the spectrum not being a load doing real work: the same ground motion is enormous for one part of the same tank and nearly absent for another, because the two parts ask it different questions.

The response spectrum of that record, at one damping ratioThe peak pseudo-acceleration of a single-degree-of-freedom structure against its natural period, for the record above, at 5% damping. Each curve is 44 complete time integrations — one structure per point. At zero period the structure is rigid and rides the ground, so the curve starts at the peak ground acceleration of 3.5 m/s²; it peaks at 6.06 m/s² at a period of 0.39 s, an amplification of 1.73.,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,0123456780123456natural period (s)spectral acceleration (m/s²)the ground's own peak, 3.5 m/s²5% damping — peak 6.06 m/s² at 0.39 s
Fig. 3 A spectrum computed rather than quoted — each point one complete time integration of one oscillator. At zero period the structure rides the ground; the peak is at 0.39 s and the long-period end falls away.

The overturning, which needs two heights

The base shear is one of two answers a tank needs, and the other one is where it acts.

Each mass has its own height, and the two are very different. The impulsive mass sits at about 41% of the liquid depth; the first convective mass sits at 66%. So the overturning moment is not the base shear times any single height, and the shares that make it up are different from the shares that made the shear.

The convective mass supplies 3% of the shear and rather more of the moment, because it is acting nearly two thirds of the way up. That is a small effect here and it is the general shape of a result worth carrying: a mass that is unimportant for force can be less unimportant for moment, because a moment weights everything by its height.

There is a further subtlety in what the moment is being computed for. The moment just above the base plate decides the wall’s own vertical compression and its anchorage; the moment below the base, including the pressures on the tank’s floor, decides the foundation. The two differ, and the modal heights for them are different numbers derived from the same flow field — which is why published tables give two sets and why using the wrong one is a common and entirely invisible error.

None of that is a refinement. On a slender tank the overturning moment is what governs, because it decides whether the shell lifts off its foundation on the windward side — and a shell that lifts and sets back down is a shell doing something no elastic analysis of it contains.

The wave, which is a displacement

The sloshing mass takes almost no force and it decides something else.

The wave stops growing with the tankSloshing wave height against the tank's radius, at one ground motion and one liquid depth. It rises as the square root of the radius while the convective period is short enough to sit on the spectrum's constant-velocity branch, and then flattens: past the corner the spectral acceleration falls as 1/T² while the period rises as the square root of the radius, so the product stops depending on the radius at all. A 12 m tank asks for 0.27 m of freeboard and a tank twice the size asks for 0.21. Nothing on this axis overtops the 0.50 m provided.0510152000.10.20.30.40.50.6tank radius (m)wave height (m)freeboard 0.50 m0.27 m herea bigger tank does not get a bigger wave
Fig. 4 Sloshing wave height against the tank’s radius. It rises as the square root of the radius and then flattens, because past the spectrum’s corner the acceleration falls as fast as the period rises.

A wave height is proportional to RSa(Tc)/gR \cdot S_a(T_c)/g, and in the falling branch SaS_a goes as 1/Tc21/T_c^2 while Tc2T_c^2 goes as RR — so the two dependencies cancel exactly and the wave stops growing. A tank twice the size gets the same wave.

The consequence for design is that freeboard is not scaled with the tank. It is a roughly constant height, set by the ground motion’s long-period content rather than by the vessel, and a large tank therefore needs proportionally less of it.

Where the freeboard is insufficient, the wave hits the roof — and then the convective mass, which had been contributing almost nothing, is stopped by a structure. The force it delivers over that short contact is impulsive by every argument about a load that is over before anything has moved, and a roof designed for a snow load meets a mass of liquid arriving with a velocity.

The pressure on the wall

The two masses put pressures on the wall with opposite shapes, which is the clearest picture of what the split means physically.

Two pressures with opposite shapes, on the same wallThe pressure on the windward meridian of a 18 m tank holding liquid 12 m deep, with the hydrostatic pressure drawn for scale. The impulsive part is uniform with depth — a rigid body accelerated sideways produces a pressure that does not know how deep it is — and reaches 43.1 kN/m², which is 37% of the hydrostatic pressure at the base. The convective part is the opposite shape: nearly nothing at the bottom and 2.78 kN/m² at the surface, because a wave is a disturbance of the free surface and has to die away with depth. Each is drawn at the amplitude that returns its own mass's base shear rather than to a shape.020406080100120024681012pressure (kN/m²)height above the base (m)impulsive acts at 41%convective at 66%hydrostaticuniform with depth, and 37% of the hydrostatic peak
Fig. 5 The pressure on the windward meridian, with the hydrostatic pressure drawn for scale. The impulsive part is uniform with depth and the convective part is nearly nothing at the bottom.

The impulsive part is uniform with depth, because a rigid body accelerated sideways produces a pressure that does not know how deep it is, and it reaches 43.1 kN/m² — 37% of the hydrostatic pressure at the base. The convective part is the opposite: a disturbance of the free surface, dying away downward, and largest at the top.

So the two masses load different parts of the wall, and a wall designed for the total force in the wrong distribution is wrong in both halves. It is also why the two cannot simply be added: their peaks occur at different times, being responses at periods twenty times apart, and the accepted treatment is to combine them as the square root of the sum of their squares rather than arithmetically — which for this tank makes almost no difference to the shear and rather more to the wave.

What the wall does with it

A tank wall carries a horizontal pressure as a hoop force, and that is where this essay meets the static half of the collection.

The free body that makes a hoop force a pressure times a radiusHalf a ring cut along a diameter, with the pressure drawn normal to the wall wherever the wall is. Vertical equilibrium of the half ring is the whole derivation: the pressure acts over the projected width 2R whatever the shape of the arc, the two cut faces carry N each, so N = pR — 1200 kN per metre here at 0.5 MPa on a 2.4 m radius. The result contains no wall thickness, no second moment, and no length along the pipe, which is why a hoop force is the one internal force in this collection that arrives with no lever arm attached to it. The stress does contain the thickness — 150 MPa at 8 mm — but the force does not, and a thicker wall carries exactly the same force at a lower stress.NN2R = 4.80 m — the projected widthp = 0.5 MPap · 2R = 2N ⇒ N = pR = 1200 kN/mno thickness in it, no second moment, no length
Fig. 6 The free body that makes a hoop force a pressure times a radius. Vertical equilibrium of a half ring is the whole derivation, and the result contains no wall thickness at all.

The force that is only a radius is N=pRN = pR, and it arrives with no lever arm — the one internal force in this collection with none. The seismic pressure adds to the hydrostatic one and the hoop force follows directly.

The reinforcement a tank wants most is not at the bottomHoop force up the wall of a 18 m tank holding 12 m of liquid, with the wall cast into its base slab. The dashed line is the membrane answer — γ(H − x)R, a triangle with its peak at the base — and it is what every hand calculation starts from. The solid line is what the wall actually carries. At the base the hoop force is **zero**, because a wall held there cannot move outwards and there is no hoop strain to go with a hoop force; the pressure is carried in vertical bending instead, at a fixing moment of 83.8 kNm/m. The membrane solution is recovered about 3.96 m up, and in between the two exchange the load. The peak is 805 kN/m at 25% of the height — 76% of the triangle's peak, and a third of the way up the wall.02004006008001000024681012hoop force (kN/m)height up the wall (m)peak 805 at 25% of the heightzero, at the basemembraneγ(H − x)Rheld atthe base
Fig. 7 Hoop force up the wall of a tank cast into its base. The membrane answer is a triangle peaking at the base; the real wall carries zero hoop force there and a peak a quarter of the way up.

The complication is at the base. A wall cast into its slab cannot move outward there, so it has no hoop strain and therefore no hoop force, and the pressure is carried in vertical bending instead at a fixing moment of 83.8 kNm/m. The peak hoop force is a quarter of the way up. Reinforcement placed by the membrane triangle is in the wrong place and misses the moment.

What the anchorage is being asked

The last part of a tank’s seismic design is at its perimeter, and it is a connection problem of exactly the kind this collection has already met.

A tank shell resists overturning by a couple: compression on the leeward side, tension on the windward. The tension has to be delivered into the foundation by holding-down bolts or by the weight of the shell and its contents over the uplifted zone — and where the bolts do the work, their capacity is a cone of concrete rather than anything about the steel.

Many tanks are not anchored at all. An unanchored tank resists overturning by lifting one side of its base plate, at which point the base plate acts as a membrane, the shell’s compression on the far side concentrates over a short arc, and the whole thing rocks. That is a perfectly acceptable design and it is not one any linear analysis describes: the mechanism is a large displacement, the restoring force comes from the liquid’s own weight over the uplifted area, and the period lengthens as the rocking grows.

The failure it protects against is worth naming. A rigidly anchored tank that cannot rock puts its whole overturning moment into the shell as vertical compression, and a thin cylinder under vertical compression buckles into the diamond pattern known as elephant’s foot at a stress far below anything the material would suggest. Letting the tank lift is a way of limiting the force it can attract — made weaker on purpose, applied to a vessel.

Where the model stops

The tank was assumed rigid. The convective solution is exact for a rigid wall, and a steel tank’s wall is not — its flexibility raises the impulsive period, and for a thin shell it can also couple with the shell’s own breathing modes, which the two-mass model has no representation of at all.

Only the first convective mode matters, and only mostly. The higher sloshing modes carry a few per cent of the mass between them, and they matter for the same reason higher modes matter anywhere.

How much of the mass each mode carries, over five modesThe effective modal mass of each of the five modes of a frame of five storeys, as a percentage of the total. Mode 1 carries 88.0% and mode 2 8.7%; two modes are needed to reach 90% of the mass. The masses sum to exactly the total, which is a property of the eigenvectors rather than a normalisation applied afterwards.five modes · 1500 ttwo modes reach 90% of itmode 1 · 0.6 s88.0%88.0% cumulativemode 2 · 0.21 s8.7%96.7% cumulativemode 3 · 0.13 s2.4%99.1% cumulativemode 4 · 0.1 s0.8%99.8% cumulativemode 5 · 0.09 s0.2%100.0% cumulative
Fig. 8 How much of the mass each mode carries in a different system, for comparison. The first mode carries most of it and the masses sum to exactly the total, which is a property of the eigenvectors.

Nothing here is about the tank’s own modes. A thin shell has breathing modes of its own, and a structure has more than one period applies to the container as much as to the contents.

A tank on the ground is not a tank on legs. Everything above assumes the tank’s base is on the ground and the impulsive period is short. An elevated tank has a long structural period, the impulsive mass sits somewhere quite different on the spectrum, and the two periods can approach each other — at which point the two masses interact rather than being independent.

Three modes of a seven-storey frameThe first three mode shapes of a seven-storey shear frame, from the eigenvalue problem rather than sketched. Mode 1 has a period of 0.72 s, no node and carries 86.2% of the mass; Mode 2 has a period of 0.24 s, one node and carries 9.0% of the mass; Mode 3 has a period of 0.15 s, two nodes and carries 2.9% of the mass. The nth mode crosses the axis n−1 times, which is a theorem rather than a drawing convention.seven floors · 300 t eachmode 10.72 s · 1.39 Hz86.2% of the massno nodemode 20.24 s · 4.1 Hz9.0% of the massone nodemode 30.15 s · 6.63 Hz2.9% of the masstwo nodes
Fig. 9 The general case of that difficulty. Modes exist in a sequence with the mass distributed among them, and treating a structure as having one is a decision about how much of the mass the first one carries.

What the picture cannot show

The two masses are a fiction. There is no part of the liquid that is impulsive and no part that is convective; there is one flow field, decomposed into modes for convenience, and the decomposition is exact in what it predicts and imaginary in what it depicts.

Nor does any figure show what a partially full tank does. Everything here is a function of H/RH/R, so a tank draws a path across the first figure as it empties, and its convective period lengthens while its impulsive mass falls — so the worst case is at some fill level between empty and full, and which one it is depends on the spectrum rather than on the tank.

The first three modes of a simply supported beamThree modes of a simply supported beam of 12 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 2.18 Hz with βL = 3.1416; Mode 2 is at 8.73 Hz with βL = 6.2832; Mode 3 is at 19.63 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, 12 m spaneach shape scaled to its own peakmode 1: 2.18 HzβL = 3.1416mode 2: 8.73 HzβL = 6.2832mode 3: 19.63 HzβL = 9.4248
Fig. 10 The familiar version of a modal decomposition. A beam’s modes are drawn from a general solution with the constants fixed by the supports, and their frequencies go as the square of the mode number.

The generalisation

The habit worth carrying is to ask, of any mass in a structure, whether it is attached.

A liquid is the extreme case and it is not the only one. Grain in a silo moves partly with the walls and partly not, and it also has friction, so its answer is different again. A suspended ceiling, a raised floor, a stack of stored goods, a crane’s load hanging on a rope: each is a mass that goes with the structure at some frequencies and not at others, and each is included at its full weight in almost every model ever built.

The general principle is that a mass contributes to a seismic force only in proportion to how much of it moves with the structure at the structure’s own period. Most of the mass moves together is the essay about how much of a building’s mass participates in its first mode; this is the same question asked of the contents rather than of the frame, and the answer is that a great deal of what a structure weighs is not, dynamically, part of it 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.

Base shearConvective massEigenvalueFree bodyFreeboardHoop tensionImpulsive massLoad pathMembrane actionModal massMode shapeNatural periodOverturningResponse spectrumSloshing