Deflection

Solved by passing it around

An indeterminate structure needs simultaneous equations, and for thirty years engineers solved them without writing any down. Clamp every joint, release one, share out what is left over, pass half of it along, and repeat — and the answer walks in, three figures correct after four cycles.

Assumes One support too many, and what it costs to know, The moment over the support, and what it buys and One deflection, without solving everything.

A continuous beam over four supports cannot be solved by statics. There are more reactions than equations, the answer depends on the beam’s stiffness, and the honest route is to write down a set of simultaneous equations and solve them. For a three-span beam that is three equations; for a twelve-storey frame it is several hundred, and in 1930 there was no way to solve several hundred simultaneous equations at all.

Hardy Cross’s answer was to not solve them.

The answer arrives in instalmentsThe hogging moment at support 1 of a three-span beam, cycle by cycle. It starts at the fixed-end moment of 53.3 kNm — the value with every joint clamped — and settles at 64.0 kNm against an exact 64.0. The error falls by about a factor of four per cycle: 14.00, 3.50, 1.95, 0.59 kNm after one, two, three and four. Two cycles is an engineering answer and nobody had to invert anything.012345678020406080cycles of distributionmoment at the support (kNm)exact: 64.0all joints clamped
Fig. 1 The hogging moment at the first interior support of a three-span beam under 10 kN/m, cycle by cycle. It starts at the fixed-end moment of 53.33 kNm — the value with every joint clamped — and settles at 64.00 against an exact 64.00 obtained from the stiffness matrix. The error falls by about a factor of four per cycle, and after four cycles it is 0.014 kNm.

Every number in that sequence was obtained by arithmetic a person can do: multiply by a fraction, halve, add. No equation was ever assembled, nothing was inverted, and the answer is exact in the limit.

The three rules

Clamp every joint. With every joint held against rotation, each span is a fixed-ended beam and its end moments are known in closed form — wL2/12wL^2/12 for a uniform load. The structure is now in equilibrium with the load and is not compatible: the joints are being held by moments nobody applied.

Release a joint and share out the imbalance. At a joint where two members meet, the two fixed-end moments do not cancel, and the difference is what the clamp is holding. Take the clamp off and the joint rotates until the members’ resisting moments balance. Since each member’s resisting moment is its rotational stiffness times the same rotation, the imbalance divides between them in proportion to stiffness:

DFi=kikwithk=4EIL\text{DF}_i = \frac{k_i}{\sum k} \quad\text{with}\quad k = \frac{4EI}{L}

Carry over half. A rotation applied at one end of a member with its far end fixed produces a moment at that far end of half the near-end value, with the same sign. So each share sent into a member deposits half of itself at the other end, which unbalances the joint there, which is why the process repeats.

What each member takes is decided before anything is distributedDistribution factors at every joint of a three-span beam of 8, 8, 8 m. At each joint the out-of-balance moment is shared between the members meeting there in proportion to 4EI/L, so a short span takes more of it than a long one — the factors are properties of the geometry and are written down once, before any arithmetic happens. The fixed-end moments the method starts from are wL²/12: 53.3, 53.3, 53.3 kNm.joint 01.000joint 10.5000.500joint 20.5000.500joint 31.0008 mFEM 53.3 kNm8 mFEM 53.3 kNm8 mFEM 53.3 kNmthe factors at a joint sum to one — nothing is created or lost in a distribution
Fig. 2 The distribution factors of the three-span beam, written down before any arithmetic happens. At each interior joint two equal members meet, so each takes half. At the pinned ends a single member takes all of the imbalance, which is why the end moment goes to zero on the first cycle and stays there. The fixed-end moments the method starts from are 53.33 kNm on every span.

That is the entire method. It is short enough to memorise and it is what an engineer of the 1930s to the 1960s used for every frame they ever designed.

The two numbers that make it work

The factor of a half deserves a second look, because it is the whole reason the process converges rather than wandering.

For a prismatic member with its far end fixed, the near-end stiffness is 4EI/L4EI/L and applying a unit rotation there produces 2EI/L2EI/L at the far end — a ratio of exactly 0.5. Every cycle therefore multiplies the residual imbalance by (a distribution factor below one) times (a half), so the residue is reduced by a factor of at least two per member per cycle, and in practice by three or four.

That is geometric convergence with a comfortable ratio. The table above shows it: 55.0, 28.8, 7.8, 2.4, 0.62, 0.16, 0.04 kNm of imbalance on successive cycles. Anyone can see when to stop, which is a property no matrix solution has — a matrix answer is either obtained or not, and a distribution answer arrives with its own error estimate visible in the last row.

The second number is the 4EI/L4EI/L itself, and it is worth stating what it hides. It is the moment needed at one end of a member to produce a unit rotation there while the far end is held, and it is the same quantity a stiffness matrix carries in its diagonal. Moment distribution is not an alternative to the stiffness method; it is Gauss–Seidel iteration on the stiffness equations, performed by somebody who never wrote the matrix down. Cross did not describe it that way, and the identification came later.

It agrees with a method that shares nothing with it

The claim that the method converges to the exact answer is checkable, and on this site it is checked by computing the same beam a second way.

The same diagram, from two entirely different arithmeticsBending moments on a three-span beam under 10 kN/m. The curve is the stiffness solution; the marked points are what moment distribution reached after eight cycles of hand arithmetic — -0.0, 64.0, 64.0, 0.0 kNm at the supports against an exact 0.0, 64.0, 64.0, 0.0. The two share no code and no equations, which is the only arrangement in which agreement is evidence of anything.-0.064.064.0the largest disagreement is 2.4e-3 kNm
Fig. 3 Bending moments on the three-span beam. The shaded curve is the stiffness solution — a matrix assembled, restrained freedoms struck out, and the system solved by elimination. The marked points are what moment distribution reached in eight cycles of the arithmetic above. The two share no code and no equations, and they agree to a hundredth of a kilonewton-metre.

The comparison is the point rather than the agreement. A method that is checked against its own last iteration is checking convergence and not correctness, and an iterative scheme can converge beautifully to the wrong answer if its distribution factors are wrong. Agreement between two routes that share no arithmetic is the only evidence that either is right, which is why this site’s gate computes both.

3 continuous spans against 3 simple onesThe bending moment in a continuous beam, solved by the stiffness method, drawn over the moment in the same spans made simply supported. The peak sagging moment falls from 30.6 to 19.6, and a hogging moment of 24.5 appears over the supports where there was none.moment19.6 sagging24.5 hogging30.6 if the spans were simplereactions 14.0 38.5 38.5 14.0 — the inner supports carry far more than a sharethe continuous case needed stiffness; the comparison did not
Fig. 4 The same kind of beam with its simply supported comparison drawn behind: continuity converts part of the midspan sagging moment into hogging over the supports, and the peak moment falls by a fifth. This is the result moment distribution exists to produce, and the reason engineers wanted it badly enough to iterate by hand for it.

The distribution factors are the design

The arithmetic gets more interesting when the spans are not equal, because the distribution factors stop being halves and start expressing a design decision.

The answer arrives in instalmentsThe hogging moment at support 1 of a three-span beam, cycle by cycle. It starts at the fixed-end moment of 30.0 kNm — the value with every joint clamped — and settles at 72.4 kNm against an exact 72.4. The error falls by about a factor of four per cycle: 11.98, 5.42, 0.71, 0.12 kNm after one, two, three and four. Two cycles is an engineering answer and nobody had to invert anything.012345678020406080cycles of distributionmoment at the support (kNm)exact: 72.4all joints clamped
Fig. 5 The same beam with spans of 6, 10 and 6 m. The distribution factors at the first interior joint are 0.625 to the short span and 0.375 to the long one — stiffness goes as 1/L, so the short member takes the larger share of every imbalance. The support moment converges to 72.38 kNm, against an exact 72.38.

That the short span attracts more moment is the practical content of the whole method, and it is a fact about stiffness rather than about load. A designer who deepens one member has attracted moment into it; one who lengthens a span has pushed moment away from it. In a hand method those consequences are visible in the second line of the table, before the calculation is finished — which is what a designer wants, because the question being asked is usually what should this be rather than what is this.

There is a lesson here that survived the method. Moment follows stiffness, and stiffness is chosen, so an indeterminate structure’s internal forces are partly a matter of design intent rather than of external circumstance. That statement is easy to say and hard to feel from a computer output; it is unavoidable when the sharing is done by hand.

What it is not good at

Sway. A frame free to move sideways has an extra unknown per storey that no amount of joint rotation addresses, and the treatment is a separate outer iteration: guess the sway, distribute, find the out-of-balance shear, correct the sway, repeat. It works and it is laborious, and it is where the method’s economy ends. A leaning portal frame is exactly this case.

Imposed displacements. A settlement or a temperature effect enters as a set of fixed-end moments computed from the movement, which is fine, but the physical story the method tells is much less natural — the engineer is now distributing moments that arrived from a support going down rather than from a load going on.

3 continuous spans against 3 simple onesThe bending moment in a continuous beam whose support 1 has settled by 0.012. Three curves: the moment the load causes, the moment the settlement causes on its own — dashed, peaking at 10.6, and in equilibrium with no applied load at all — and their sum, which is what the beam carries, peaking at 31.6 against 24.5 without the settlement. The settlement field is proportional to EI: a stiffer beam is punished harder for the same movement, which is the opposite of every intuition load-carrying gives.this support 0.012 lowmoment19.6 sagging24.5 hogging30.6 if the spans were simplereactions 14.0 38.5 38.5 14.0 — the inner supports carry far more than a sharethe continuous case needed stiffness; the comparison did not
Fig. 6 A support that has moved 12 mm, and the moments that follow with no load at all. Moment distribution handles this by converting the movement into fixed-end moments and proceeding as usual; what it cannot easily convey is that the answer is proportional to EI, so a stiffer beam is punished for its stiffness — the finding that makes settlement dangerous.

Anything with axial deformation in it. The method assumes members do not change length, which is excellent for beams and frames and useless for a truss. The competing hand method — virtual work with a unit load — has the opposite profile: it handles trusses beautifully and continuous frames laboriously.

The deflection at x = 3, by virtual workThree diagrams: the moment from the real load, the moment from a unit load placed where the answer is wanted, and their product. The area under the third, divided by EI, is the deflection — 197.50 here. No standard case was consulted, so the method works for any load pattern at all.real Mpeak 32.0a unit load, here and nowhere elseunit mM × marea ÷ EI = 197.50the unit load is the only place the question 'deflection where?' is asked
Fig. 7 The other hand method: a unit load placed where the answer is wanted, and an integral of the product of two moment diagrams. It gives one number for one question and needs a fresh calculation for the next one, where moment distribution gives every moment in the frame at once. The two are the flexibility and stiffness approaches, and which is less work depends entirely on which the structure has more of — redundancies or degrees of freedom.

The tables in every handbook are this method, frozen

Open any structural handbook at the page of continuous-beam coefficients — 0.080, 0.100, 0.025, the little grid of numbers for two, three, four and five equal spans — and what is printed there is the output of this method, run once by somebody in the 1930s and copied ever since.

The three-span case above is the entry 0.100wL20.100\,wL^2 at the interior supports, which is 64.0 kNm for w=10w = 10 and L=8L = 8, exactly as computed. The two-span case is 0.125wL20.125\,wL^2, which is wL2/8wL^2/8 — the same number a simply supported beam carries at midspan, arriving at a support instead. Those coefficients are exact for the case they describe and they are the reason a generation could design continuous beams without doing any of the arithmetic above.

One support too manyThe same uniformly loaded beam with three sets of restraints, and the bending moment in each. Adding restraint moves moment from mid-span to the supports and lowers the peak — but only the first case can be solved by statics.simply supportedstatics alonesag 32.0propped at one endneeds stiffnesssag 18.0hog 32.0built in at both endsneeds stiffnesssag 10.7hog 21.3the load never changes; only what is holding the endsthe built-in case peaks at two-thirds of the simple span's moment
Fig. 8 Where those coefficients come from: the same beam with restraints added one at a time. Statics gives the first case and nothing beyond it; the other two need the stiffness, and the peak moment falls each time a restraint is added because part of it has been moved to the support. Moment distribution is how those numbers were obtained before a matrix could be inverted.

The tables have two limitations that the method itself does not, and both are worth naming because the tables outlived the understanding. They assume equal spans, and the section on distribution factors above shows how quickly unequal ones change the answer. And they assume a pattern of loading — usually all spans loaded — where the governing case for a support moment is often alternate spans loaded, which is the same question an influence line answers and which the coefficients handle by having a second row.

Continuity is bought and paid for

It is worth being explicit about what all this arithmetic is for, because the answer is not obvious from the numbers.

A simply supported 8 m beam under 10 kN/m carries 80 kNm at midspan. Make three of them continuous and the interior spans carry 40 kNm at midspan and 64 at the supports — so the largest moment anywhere has fallen from 80 to 64, a reduction of a fifth, for no material and no change of span. That is the return on continuity, and it is why the method existed.

What was paid for it is threefold, and each cost is somewhere else in this collection. The structure is now indeterminate, so a support that settles produces moments with no load at all. The peak moment is at a support, where the beam is usually shallowest in a haunched design and where the connection has to carry it across a gap. And the beam is now sensitive to the pattern of loading rather than only to its magnitude, which doubles the number of cases to check.

Continuity is a trade of one large moment for two smaller ones plus a set of new sensitivities, and moment distribution is the instrument that made the trade computable by hand. That it also gave engineers an intuition for where moment goes was a side effect, and it was the part that was lost when the arithmetic moved to a machine.

What it was really for

Cross published the method in 1930 in a paper of ten pages, and it was in universal use within five years. The usual account is that it saved arithmetic, and it did — but the deeper reason it took over is the one visible in the tables above.

The method produces intermediate results that mean something. The first line is the structure with every joint clamped; the second is what one joint does when released; the third is the effect of that release travelling to its neighbours. An engineer following it is watching a moment move through a structure, and the number they end with is one they have watched arrive. It is a calculation with a narrative, and the narrative is the mechanism.

That is worth stating precisely because it is the property a matrix solution lacks, and the loss was real. A generation of engineers who had distributed moments by hand had an intuition for which member attracts what, and it came from the arithmetic rather than from a course. The replacement of the method by matrix analysis in the 1970s was correct on every technical ground — sway, three dimensions, thousands of members — and it removed the only part of the calculation a person could see.

There is one more thing the method quietly taught, and it is visible in the very first line of every table. The fixed-end moments are the answer for a structure with every joint clamped, which is the stiffest possible arrangement, and the distribution process only ever takes moment away from the supports and gives it to the spans. So the first line is a bound: the support moment of a real continuous beam is never larger than its fixed-end moment, and a designer in a hurry could stop at line one and be conservative at the supports. Half the rules of thumb in continuous-beam design are that observation, and the count of unknowns is what says why the bound exists — clamping every joint makes the structure maximally redundant, and every release moves it back toward statics.

What the picture cannot show

Everything above is elastic and linear. The carry-over factor of a half assumes a prismatic member behaving elastically; a member that has cracked, yielded or has a variable section has different stiffness and carry-over factors, which is arithmetic rather than a difficulty and which every textbook tabulates.

The convergence rate depends on the structure. The three-span beam above converges in four cycles because its distribution factors are around a half. A frame with many members at a joint, or with very unequal stiffnesses, converges more slowly, and a structure with a near-mechanism in it converges very slowly indeed — the same ill-conditioning a matrix solver would meet, showing up as patience rather than as a warning.

Nothing here is a check. The method produces a set of moments in equilibrium and compatible, and it has no opinion about whether they are acceptable. Everything about capacity, deflection and stability is a separate question asked afterwards.

Where the ladder goes

The first rung is sway, which is where the method’s outer iteration lives and which is the case every real frame has.

The second is the identification with Gauss–Seidel, and what it says about the structure: iterative methods converge fast when a system is diagonally dominant, and a structure’s stiffness matrix is diagonally dominant exactly when each joint is better connected to itself than to its neighbours. The method converges quickly on ordinary structures because ordinary structures are locally stiff, which is a statement about buildings rather than about arithmetic.

The third is the one this collection keeps returning to: what happens when the moments obtained above exceed what a section can carry. The elastic distribution is one equilibrium solution among many, and the plastic one redistributes it entirely — which means the careful sharing above describes the structure only until first yield, and the collapse load does not depend on 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.

Carry overContinuityDistribution factorFixed end momentIndeterminacyIterationMoment distributionStiffness