Updated 23 July 2022
Reading time: 19 mins

How to Analyse Indeterminate Beams using the Moment Distribution Method

Learn how to draw shear force & bending moment diagrams for indeterminate beams using the moment distribution method
[object Object]
by Dr Seán Carroll

In this tutorial we’ll explore the moment distribution method. This is an excellent technique for quickly determining the shear force and bending moment diagrams for indeterminate beam and frame structures. In this tutorial, we’ll focus on applying the moment distribution method to beams.

We’ll start by getting a clear understanding of the steps in the procedure before applying what we’ve learned to a more challenging worked example at the end. If you’d prefer to watch me explain the solution, you can watch video below.

This tutorial is based on my course, Indeterminate Structures and the Moment Distribution Method. After you finish this tutorial, you can go there to keep learning.

Indeterminate Structures & The Moment Distribution Method

Unlock indeterminate structures using the moment distribution structural analysis method.

After completing this course...

  • You will have a deep understanding of how structures respond to loading.
  • You will clearly understand the methods used to analyse statically indeterminate beams and frames.
  • You will be able to analyse both non-sway and sway frames – essential when designing real world structures.

If you’re reading this, I’m assuming you’re already comfortable drawing shear force and bending moment diagrams for statically determinate beams. If you’re not, work your way through this tutorial first.

1.0 Introduction to the Moment Distribution Method

Let’s start by summarising the key features of the moment distribution method; the technique seeks to identify the bending moments at internal joints through an iterative process of applying balancing and redistribution moments.

Iterations continue, successively reducing and moment imbalance, until moment equilibrium is achieved at all joints in the structure. With these internal moments established, span moments, shear forces and support reactions are determined using free body diagrams and simple statics.

We start by fixing all internal joints against rotation. This is often referred to as locking the structure. In a multi-span beam, this results in a series of beam segments, isolated from each other by locked joints. Next, we determine the bending moments that develop at each locked joint as a result of the loading on each beam segment. Typically there will be a moment imbalance at each joint. At this point we enter the iterative moment balancing process; for each joint in turn, we:

  • apply a balancing moment to eliminate the imbalance
  • distribute the balancing moment between the members meeting at the joint, in proportion to their flexural stiffnesses
  • carry over 50% of the distributed moment to the other end of each of the members meeting at the joint  -assuming the adjacent joint is capable of resisting moments – we’ll clarify this below)

This carry-over moment will now unbalance the joints in the structure again. This means we need to repeat the balance and distribution process. This process repeats over and over again, however, with each iteration the moment imbalances in the structure become smaller and smaller.

As usual, the only real way to make sense of this is to watch it in action, so let’s work our way through a simple example.

2.0 A Simple Example

Consider the simple two-span continuous indeterminate beam in Fig 1. below. This beam is statically indeterminate because there are more than three unknown reactions. So, we can’t determine their values with only 3 equations of statics. This makes the structure a prime candidate for a moment distribution analysis.

Moment distribution method 1 | EngineeringSkills.com

Fig 1. Two-span continuous indeterminate beam.

2.1 Locking the joints against rotation

The first step is to lock any joint not already fixed against rotation, so in this case, that’s joint B. The structure is now essentially two isolated spans, AB and BC with the loading on one span having no impact on the other.

Moment distribution method 2 | EngineeringSkills.com

Fig 2. Beam with all joints fixed against rotation.

Next, we can determine the moments that develop at each locked joint due to the span loading in each beam segment. The quickest way to do this is by using tables of fixed-end moments. A table of common fixed-end moments is provided, Fig. 3.

Moment distribution method 3 | EngineeringSkills.com

Fig 3. Common fixed-end moments.

We can see that in this case, the joint moments that develop due to the loading between AB are 100kNm100\:kNm. The joint moments induced by the loading on BC are 208.3kNm208.3\:kNm. Another way to think about the fixed-end moments is that they’re the reaction moments that would develop at the fixed ends, in response to the applied loading.

Moment distribution method 4 | EngineeringSkills.com

Fig 4. Fixed-end moments that develop at the locked joints.

2.2 Balancing the joint

If we now focus on joint B, we can see a moment imbalance due to the clockwise 100kNm100\:kNm from AB and the counter-clockwise 208kNm208\:kNm from BC. If we now imagine releasing this joint, it would, in theory, experience a net counter-clockwise moment of 108kNm108\:kNm. In response to this, we will add a clockwise balancing moment of 108.3kNm108.3\:kNm to eliminate the moment imbalance, Fig. 5.

Moment distribution method 5 | EngineeringSkills.com

Fig 5. Moment imbalance at B and the applied balancing moment.

2.3 Distributing the balancing moment

In adding the balancing moment to the joint, it must be distributed between all of the members that meet at the joint. This distribution is done in proportion to the flexural stiffnesses of the members meeting at the joint. We talk more about how these stiffnesses are determined below but for now, let’s just assume that beam segment AB is twice as stiff as segment BC.

This would mean that the balancing moment is distributed between AB and BC in a 2:1 ratio. So, 72.2kNm72.2\:kNm into member AB and 36.1kNm36.1\:kNm into BC. As a result, the joint is now balanced with moments of 172.2kNm172.2\:kNm balancing each other on either side of the joint.

Moment distribution method 6 | EngineeringSkills.com

Fig 6. Distributing the balancing moment.

2.4 The carry-over moment

We’re not finished just yet, the final step is to propagate the carry-over moment to the opposite ends of each beam segment. Since the structure has fixed supports at A and C, we can’t apply balancing moments at B without a corresponding moment developing at the outer supports.

As the supports are fully fixed, 50%50\:\% of the balancing moments will be carried over, Fig 7. Had A or B been a pin or roller support, which offers no resistance to rotation, the carry-over moments would be zero.

Moment distribution method 7 | EngineeringSkills.com

Fig 7. Carry-over moments.

Now we can add up the final moments at each joint, Fig. 8. This was a relatively simple distribution that only required a single balancing iteration. For structures with more than one internal locked joint, multiple balancing iterations are required – we’ll see an example of this below.

Moment distribution method 8 | EngineeringSkills.com

Fig 8. Adding up the final joint moments for the structure.

For more complex, multi-iteration structures, we can use a table to help keep track of the analysis. A simple table tracking this analysis is shown below, Table 1.

JointABBC
MemberABBABCCB
Distribution Factor0.660.33
Fixed-end moments-100100-208.3208.3
Distribution72.236.1
Carry-over36.118.1
Final-63.9 kNm172.2 kNm-172.2 kNm226.4 kNm

3.0 Relative Stiffnesses and Moment Distribution Factors

In the last example we assumed a ratio of stiffnesses between the beam segments. This was just so that we could work through an analysis from start to finish, without getting distracted.

We actually need to more carefully evaluate the relative stiffness of each member. We’ll discuss how to do this next. The flexural stiffness of a beam element depends on the following factors:

  • The beam material, captured using Young’s modulus, EE
  • The cross-section shape, captured with the second moment of area, II
  • The length of the beam between supports, LL
  • The degree of rotational restraint at supports.

In terms of rotational fixity, the beam segment will have one of two possible stiffnesses. If there is complete rotational fixity or resistance to rotation provided at both ends of the beam, the stiffness will be,

K=4EIL(1)\boxed{K=\frac{4EI}{L}} \tag{1}

However, if one end of the beam is pinned and therefore has no resistance to rotation at that end, the beam stiffness will be,

K=3EIL(2)\boxed{K=\frac{3EI}{L}} \tag{2}

We can more readily summarise the possible fixities and their associated stiffnesses graphically, Fig 9.

Moment distribution method 9 | EngineeringSkills.com

Fig 9. Beam element stiffnesses based on rotational fixity at supports.

With this additional information, let’s return to our example question and update it with the correct element stiffnesses. Considering AB first; it’s fixed against rotation at A and since it’s continuous over support B, it will have a stiffness of,

KAB=4EILKAB=4EI8KAB=0.5EI(3)\begin{align*} K_{AB} &= \frac{4EI}{L}\\\\ K_{AB} &= \frac{4EI}{8}\\\\ K_{AB} &= 0.5EI \end{align*} \tag{3}

The same reasoning applies to the segment BC and so its stiffness will be,

KBC=4EILKBC=4EI10KBC=0.4EI(4)\begin{align*} K_{BC} &= \frac{4EI}{L}\\\\ K_{BC} &= \frac{4EI}{10}\\\\ K_{BC} &= 0.4EI \end{align*} \tag{4}

Now we can calculate more accurate distribution factors governing the degree to which the balancing moment is distributed into each member. The distribution factor for member AB is given by,

DFAB=KABKAB+KBC=0.50.5+0.4=0.56(5)DF_{AB} = \frac{K_{AB}}{K_{AB}+K{BC}} = \frac{0.5}{0.5+0.4} = 0.56 \tag{5}

Since there is only two members meeting at joint B, we know that the distribution factor for BC is DFBC=1DFABDF_{BC} = 1-DF_{AB}. For completeness however, we’ll also state the complete expression,

DFBC=KBCKAB+KBC=0.40.5+0.4=0.44(6)DF_{BC} = \frac{K_{BC}}{K_{AB}+K{BC}} = \frac{0.4}{0.5+0.4} = 0.44 \tag{6}

Repeating the distribution analysis with the correct distribution factors reveals the actual final set of support moments for the structure, Table 2.

JointABBC
MemberABBABCCB
Distribution Factor0.560.44
Fixed-end moments-100100-208.3208.3
Distribution(0.56×108.3) =60.65(0.44×108.3)=47.65
Carry-over30.3323.83
Final-69.67 kNm160.65 kNm-160.65 kNm232.13 kNm

4.0 Reactions, Shear Forces and Bending Moments

Our analysis so far has revealed the bending moments at each joint. But to complete the analysis and produce shear force and bending moment diagrams, we have more work to do.

Fortunately, the missing information can be easily obtained using simple statics and free body diagrams. Consider a free body diagram of beam segment AB, Fig 10.

Moment distribution method 10 | EngineeringSkills.com

Fig 10. Free body diagram of beam segment AB.

We have two unknowns here, the shear forces at A and B. We can determine the shear force immediately to the left of B, VBLV_{B_L}, by evaluating the sum of the moments about A. We’ll assume clockwise moments are positive.

MA=0(+)(7)\sum M_A=0\:\:(\curvearrowright +) \tag{7}
69.7+160.7+(4×100)8VBL=0(8)-69.7+160.7+(4\times 100)-8V_{B_L} = 0 \tag{8}
VBL=61.38kN(9)\boxed{V_{B_L}=61.38\:kN} \tag{9}

Next we can evaluate the sum of the vertical forces to determine VAV_A.

Fy=0(+)(10)\sum F_y=0\:\:(\uparrow +) \tag{10}
VA100+61.38=0(11)V_A-100+61.38=0 \tag{11}
VA=38.62kN(12)\boxed{V_A=38.62\:kN} \tag{12}

Because the beam segment is subject to a single point load, we know there will be a peak moment under the point load and that the moment will vary linearly between this peak and the two support moments.

Similarly, we know that the shear force diagram will consist of two horizontal straight lines. At this stage in the process, you’ll need to have a good understanding of basic shear force and bending moment diagram construction – review this tutorial now, if needed.

We can determine the internal bending moment under the point load by making a cut in the structure at this location to reveal the internal moment, MDM_D, Fig 11.

Moment distribution method 11 | EngineeringSkills.com

Fig 11. Free body diagram of beam segment AD

Evaluating the sum of the moments about the cut yields,

MD=0(+)(13)\sum M_D=0\:\:(\curvearrowright +) \tag{13}
MD69.7+(4×38.62)=0(14)-M_D-69.7+(4\times 38.62)=0 \tag{14}
MD=84.78kNm(15)\boxed{M_D=84.78\:kNm} \tag{15}

We can repeat this process now for beam segment BC, Fig. 12.

Moment distribution method 12 | EngineeringSkills.com

Fig 12. Free body diagram of beam segment BC.

Evaluating the sum of the moments about B first.

MB=0(+)(16)\sum M_B=0\:\:(\curvearrowright +) \tag{16}
160.7+232.1(25×10×5)10VC=0(17)-160.7+232.1-(25 \times 10\times 5)-10V_C=0 \tag{17}
VC=132.14kNm(18)\boxed{V_C=132.14\:kNm} \tag{18}

Now evaluating the sum of the forces in the vertical direction.

Fy=0(+)(19)\sum F_y = 0 \:\:(\uparrow +) \tag{19}
VBR+132.14(25×10)=0(20)V_{B_R}+132.14-(25\times 10)=0 \tag{20}
VBR=117.86kN(21)\boxed{V_{B_R}=117.86\:kN} \tag{21}

This beam segment is subject to a uniformly distributed load (UDL) which means the shear force will vary linearly between VBL=117.86kNV_{B_L}=117.86\:kN on the left and VC=132.14kNV_C=132.14\:kN on the right.

The location of maximum moment will be where the shear force is zero. So, to determine this distance, xx to the right of B, we can divide the shear on the left by the rate at which it’s reducing, 25kN/m25\:kN/m, due to the UDL,

x=117.8625=4.71m(22)x=\frac{117.86}{25} = 4.71\:m \tag{22}

This is where we can make our cut to reveal the maximum span moment, Fig. 13

Moment distribution method 13 | EngineeringSkills.com

Fig 13. Free body diagram of beam segment identifying the location of the maximum moment.

Evaluating the sum of the moments about the cut location,

Mx=0(+)(23)\sum M_x=0\:\:(\curvearrowright +) \tag{23}
Mx160.7+(117.86×4.71)(25×4.7122)=0(24)-M_x-160.7+(117.86\times 4.71)-\left(25\times \frac{4.71^2}{2}\right)=0 \tag{24}
Mx=117.1kNm(25)\boxed{M_x = 117.1\:kNm} \tag{25}

Now we have all the information we need to sketch out the complete shear force and bending moment diagrams for this beam, Fig 14.

Moment distribution method 14 | EngineeringSkills.com

Fig 14. Sketch of the shear force and bending moment diagram.

This analysis walkthrough has demonstrated the complete moment distribution method. You should now have a good understanding of the various steps involved and recognise the potential this technique offers for manual analysis of both beam and frame structures. In the next section, we’ll work through a slightly more complex beam example that will require us to implement a multi-iteration solution.

5.0 Moment Distribution Method: Multi-iteration Worked Example

For this example, consider the multi-span continuous beam shown below, Fig 15. We will again assume EIEI is constant for this beam. With this example, we’ll see that a relatively complex span and loading arrangement can be handled quite quickly using the moment distribution process.

With some practice, the processing can become very streamlined and somewhat formulaic to implement. You can continue reading through the solution or watch the full solution video below.

Moment distribution method 15 | EngineeringSkills.com

Fig 15. Multi-span continuous beam.

5.1 Member Stiffnesses

Based on our earlier discussion of element stiffnesses, we can state the following:

KAB=4EIL=4EI8=0.5EI(26)K_{AB} = \frac{4EI}{L} = \frac{4EI}{8} = 0.5EI \tag{26}
KBC=4EIL=4EI6=0.67EI(27)K_{BC} = \frac{4EI}{L} = \frac{4EI}{6} = 0.67EI \tag{27}
KCD=3EIL=3EI4=0.75EI(28)K_{CD} = \frac{3EI}{L} = \frac{3EI}{4} = 0.75EI \tag{28}

We don’t calculate a stiffness for segment DE because no balancing moment will be transmitted into this segment. This is because we already know what the final moment will be at this location due to the fact that it’s a cantilever.

5.2 Distribution Factors

We calculate the distribution factors for each internal joint in the usual way, again noting specifically that no balancing moment will be distributed into DE and therefore the distribution factor for DC will be 1.

DFBA=0.50.5+0.67=0.43(29)DF_{BA} = \frac{0.5}{0.5+0.67} = 0.43 \tag{29}
DFBC=10.43=0.57(30)DF_{BC} = 1-0.43 = 0.57 \tag{30}
DFCB=0.670.67+0.75=0.47(31)DF_{CB} = \frac{0.67}{0.67+0.75} = 0.47 \tag{31}
DFCD=10.47=0.53(32)DF_{CD} = 1-0.47=0.53 \tag{32}
DFDC=1(due to cantilever)(33)DF_{DC} = 1\:\:(\text{due to cantilever}) \tag{33}
DFDE=0(due to cantilever)(34)DF_{DE} = 0\:\:(\text{due to cantilever}) \tag{34}

5.3 Fixed-end Moments

Considering segment AB first, Fig 16, the fixed-end moments are obtained as,

FEMAB=wL212=10×8212=53.33kNm(35)FEM_{AB} = \frac{-wL^2}{12} = \frac{-10\times 8^2}{12}=-53.33\:kNm \tag{35}
FEMBA=53.33kNm(36)FEM_{BA} =53.33\:kNm \tag{36}
Moment distribution method 16 | EngineeringSkills.com

Fig 16. Fixed-end moments, segment AB.

Next, we can evaluate beam segment BC, Fig 17.

FEMBC=Pab2L2=75×2×4262=66.7kNm(37)FEM_{BC} = \frac{Pab^2}{L^2} = \frac{75\times 2\times 4^2}{6^2}=66.7\:kNm \tag{37}
FEMCB=Pba2L2=75×4×2262=33.3kNm(38)FEM_{CB} = \frac{-Pba^2}{L^2} = \frac{-75\times 4\times 2^2}{6^2}=-33.3\:kNm \tag{38}
Moment distribution method 17 | EngineeringSkills.com

Fig 17. Fixed-end moments, segment BC.

When evaluating the fixed-end moments for segment CD we fix joint C as usual. However, since the beam is free to rotate to the right of D, i.e. it is not restrained due to the fact it’s a cantilever to the right of D, we leave this joint pinned. So, the fixed-end moment, FEMCDFEM_{CD} is calculated assuming a propped cantilever model, Fig 18.

FEMCD=wL2b=5×428=10kNm(39)FEM_{CD} = \frac{-wL^2}{b} = \frac{-5\times 4^2}{8} =-10\:kNm \tag{39}
FEMDC=0(40)FEM_{DC} = 0 \tag{40}
Moment distribution method 18 | EngineeringSkills.com

Fig 18. Fixed-end moments, segment CD.

Finally we can determine the fixed-end moments at DE and ED which are really just the cantilever moments evaluating from simple statics.

FEMDE=50kNm(41)FEM_{DE} = -50\:kNm \tag{41}
FEMED=50kNm(42)FEM_{ED} = 50\:kNm \tag{42}

5.4 Moment Distribution

Now that we’ve established all of the input information, we can construct the moment distribution table and process the distribution, Table 3. Unlike our previous example, because this structure has multiple internal joints, we’ll need to work through multiple balancing and distribution iterations to progressively reduce the moment imbalance at each support.

The first four rows of the table contain all of the structure and loading information we’ve just determined. Thereafter we are:

  1. Identifying the out of balance moment at each joint and calculating the balancing moment to distribute into each member based on the distribution factors.
  2. Once all joints are balanced, we carry over 50%50\% of the balancing moment to adjacent locked joints. Note again, that there is no carry-over moment applied to DC since the cantilever moment at this location is already determined.
  3. The carry-over moments have unbalanced the joints again, so we need to determine a new balancing moment for each joint, so we return to step one and repeat the process.

This process continues until the balancing moments being applied are sufficiently small. Practically, we can stop balancing once the moments have reduced to about 1 or 2 percent of the initial fixed-end moments. This leaves us with balanced internal joints.

JointABBCCDDE
MemberABBABCCBCDDCDEED
Distribution Factors0.430.570.470.5310
Fixed-end Moments-53.353.366.7-33.3-100-5050
Distribution(-120)-51.6(-120)-68.4(43.3)20.35(43.3)22.9550
Carry-over-25.810.18-34.225
Distribution-4.38-5.8(9.2)4.32(9.2)4.88
Carry-over-2.192.16-2.9
Distribution-0.93-1.231.361.54
Final-81.29 kNm-3.61 kNm3.61 kNm-44.37 kNm44.37 kNm50 kNm-50 kNm50 kNm

At this point, we can use free body diagrams and the equations of statics to evaluate the remaining unknown shear forces and bending moments. I won’t go through that step-by-step here because the process is pretty much the same as that demonstrated in the previous example. If you want to see every step, you can watch the solution video where I go through the complete process. The final shear force and bending moments diagrams are sketched below, Fig. 19.

Moment distribution method 19 | EngineeringSkills.com

Fig 19. Sketch of the final shear force and bending moment diagrams.

As a means of confirming this result, we can input this structure into our 2D beam and frame Python solver and obtain the same result, admittedly a lot quicker, Fig 20.

Moment distribution method 20 | EngineeringSkills.com

Fig 20. Graphical output from the 2D beam and frame solver we built in this course.

I hope this tutorial has given you a sense of how useful the moment distribution method can be. It’s a particularly handy tool for sub-frame analysis, allowing shear forces and bending moments to be quickly established.

In practice, you’re probably more likely to use a software programme to perform the analyse – not least because it makes analysis iterations faster, when for example, you need to alter the loading on the structure. Having said that, it’s still a good idea to have a manual analysis technique up your sleeve!

Although we’ve demonstrated the key features of the moment distribution method here, there’s more to cover in order to analyse frame structures, particularly sway frames. If you want to continue studying the moment distribution method with me, take a look at the following courses:

That’s all for now, see you in the next one.


If you want full access to the complete EngineeringSkills course (and code!) library, for one low annual fee, consider subscribing to the All Access Membership or Lifetime Access Membership if you never want to pay for another course again :)

All Access Annual Membership

Learn, revise or refresh your knowledge and master engineering analysis and design

Access Every Course and Tool

  • Over 833 lectures & over 170 hours of HD video content
  • Access all downloads, pdf guides & Python codes
  • Packed development roadmap of courses & tools 🏗️
  • Price Guarantee – avoid future price rises as we grow
  • Priority Q&A support
  • Course completion certificates
  • Early access to new courses
getting-started
Dr Seán Carroll
BEng (Hons), MSc, PhD, CEng MIEI, FHEA
Hi, I’m Seán, the founder of EngineeringSkills.com (formerly DegreeTutors.com). I hope you found this tutorial helpful. After spending 10 years as a university lecturer in structural engineering, I started this site to help more people understand engineering and get as much enjoyment from studying it as I do. Feel free to get in touch or follow me on any of the social accounts.

Dr Seán Carroll's latest courses.

Analytical Modelling of Plate and Shell Structures: Part 1 - Plates

Analytical Modelling of Plate and Shell Structures: Part 1 - Plates

A practical guide to the analysis of circular and rectangular plates under load, from first principles.

Fundamentals of Reinforced Concrete Design to Eurocode 2

Fundamentals of Reinforced Concrete Design to Eurocode 2

An introduction to ultimate limit state design for bending and shear with optional calculation automation using Python.

Modelling and Analysis of Non-linear Cablenet Structures using Python and Blender

Modelling and Analysis of Non-linear Cablenet Structures using Python and Blender

Learn how to combine parametric modelling, exploratory form-finding and iterative analysis techniques to simulate 3D tensile structures.

Non-linear Finite Element Analysis of 2D Catenary & Cable Structures using Python

Non-linear Finite Element Analysis of 2D Catenary & Cable Structures using Python

Build an iterative solution toolbox to analyse structures that exhibit geometric non-linearity due to large deflections.


Do you have some knowledge or expertise you'd like to share with the EngineeringSkills community?
Check out our guest writer programme - we pay for every article we publish.