Structural Engineering & Statics
where the load goes
Concept

Stiffness matrix — where it appears

The array relating the forces at a structure's freedoms to the displacements there. It is assembled by adding each member's own matrix into the rows and columns of the freedoms that member touches. Nothing in the assembly is an approximation; it is the same equilibrium and compatibility a hand method writes, in an order a machine can follow.

Named by 7 essays across 3 fields — each of them below, with the objects they name alongside it.

Three modes of a five-storey frame. The 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.

A structure has more than one period

One mass on one spring has one period. A building with eight floors has eight, each with a shape of its own, and the second one bends the building into a curve nobody drew. They are not harmonics, they do not interfere, and each behaves as though the others were not there.

dynamics · Mode shapes The check that everything adds up, and the error it cannot see. Four versions of the same 2-bay, 3-storey frame, with the global equilibrium residual each one produces — the sum of the reactions against the sum of the applied loads, as a fraction of the applied total. It is the first thing every analysis prints and it is worth having: a lost restraint and a load entered in the wrong unit both show up immediately, at 6% and 24%, because both change what the structure is carrying. The fourth bar is the point. A member whose stiffness is wrong by a factor of ten redistributes the internal forces completely — the second bar shows the change in the member forces, 18% — and the global residual is exactly zero, because the wrong answer is still in equilibrium with the same loads. Equilibrium is one equation per degree of freedom of the whole body, and a stiffness error lives entirely in the many equations underneath it. A model can satisfy every equilibrium check ever devised and be a model of a different structure.

The stiffness the load takes away

Buckling is usually taught as an event — a critical load, a bifurcation, a mode. Written as a matrix it stops being an event at all. A compressive load subtracts a stiffness from the structure, the subtraction grows with the load, and the critical load is simply where what is left reaches zero.

stability · Geometric stiffness One member's stiffness, scattered into the freedoms it touches. A member's own six-by-six stiffness matrix relates the forces at its two ends to the displacements there, and it is written in the member's own axes. Assembly is two operations and no physics: rotate it into the structure's axes, then add each of its thirty-six entries into the row and column of the global freedom that entry belongs to. Every member does the same, and the sum is the structure. The shaded rows and columns are the six freedoms this one member reaches; every other entry it contributes is exactly zero, and that is the whole reason a global stiffness matrix is sparse. Nothing here is an approximation — the result is the same equilibrium and the same compatibility a hand method writes, in an order a machine can follow.

The answer that depends on how it was divided

Every computed answer in this collection came out of a structure chopped into pieces — elements, strips, stations, trial positions. The chopping is invisible in the result and it is not neutral: some divisions give the exact answer, some give one that is always too stiff, and one of them changes nothing but the cost of getting there.

deflection · Discretisation 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 (7 − 2√2)/2, and become a complex pair, which is oscillation that grows. With no damping that merge is also where the column goes unstable. 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.

The load it cannot buckle under

Every stability calculation on this site rests on an assumption nobody states: that the load has a potential, so a critical load is where a total potential energy stops being a minimum. A load that turns with the structure it is pushing has no potential, and the static analysis of such a column returns no critical load at all — a determinant that never vanishes, for a column that fails at a perfectly finite one.

stability · Follower force Three quantities converging at three different rates. How much of the exact answer a truncated modal analysis of a ten-storey building reaches, against how many modes it keeps. The mass count is the rule — ninety per cent, reached at two modes. The base shear is ahead of it: one mode carries 85 per cent of the mass and 98 per cent of the base shear. The force in the top storey is behind it, at 80 per cent on one mode and 94 on two. The rule is written in the quantity that converges fastest, and it is checked against a quantity nobody plots.

The modes that were left out

Nobody runs every mode a model has, and the rule for how many is a mass count — enough to account for ninety per cent of the structure. The rule is written in the one quantity that converges fastest. On a twenty-storey frame two modes give the base shear to within one per cent and the force in the top storey to within twenty.

dynamics · Mode shapes Two capacity curves that agree, from two answers that do not. Base shear against roof drift for the same eight-storey frame with a soft ground storey, pushed with a fixed triangular pattern and with one recomputed from the tangent stiffness at every step. At a roof displacement of 426 mm they differ by 8.7 per cent — which is the number a capacity-spectrum procedure reads, and it is the number that is nearly the same. The storey drifts underneath these two curves differ by a factor of 3.8 at the fourth floor, and the drift is what the analysis was run for.

The pattern that stopped describing the building

A pushover analysis pushes with a load pattern chosen to resemble the first mode, and by the time the structure has done anything worth analysing it no longer has that mode. Recomputing the pattern as the frame softens changes the base shear by nine per cent and the drift at the fourth floor by a factor of four.

dynamics · Pushover The modes of five masses joined by springs and held by nothing. Five equal masses joined in a line by equal springs, with nothing holding them to the ground. The first four modes, at zero, 0.62, 1.18, 1.62 times the frequency of one mass on one spring. The first is every mass moving together with no spring stretched at all: a mode at exactly zero frequency, a real solution of the eigenvalue problem, which carries all of the mass and none of the strain. Every other mode has the ends moving against each other, and carries none of the mass under a uniform acceleration.

The modes at zero frequency

A structure held by nothing — a span being launched, a segment on a crane, a pontoon — has a mode in which it moves as one body and stretches nothing, at a frequency of exactly zero. That is a real mode, not a glitch in the stiffness matrix. It carries every kilogram of the structure under a uniform acceleration and leaves the flexible modes none at all. A load that is not uniform is a different matter: pushed suddenly from one end, the structure has no static answer to give, only an acceleration with a vibration riding on it.

dynamics · Mode shapes

Named alongside it

The objects these essays reach for when they reach for this one.

EigenvalueDegrees of freedomMode shapeConvergenceModal analysisSoft storeyBucklingCompatibilityCritical loadModal massNatural periodOrthogonality

All concepts