Simply supported beam—moment-area method

Continuing our odyssey through various ways of calculating the deflection at the center of a simply supported beam with a uniform load, today we use the moment-area method.

There are two moment-area theorems used to calculate the slopes and deflections of a beam. Here’s how they’re given in the textbook I used as an undergrad, the 3rd edition of Elementary Structural Analysis by Norris, Wilbur, and Utku.1

The first theorem is

The change in the slope of the tangents of the elastic curve between two points A and B is equal to the area under the M/EI diagram between these two points.

The second theorem is

The deflection of point B on the elastic curve from the tangent to this curve at point A is equal to the static moment about an axis through B of the area under the M/EI diagram between points A and B.

The “elastic curve” is the curve made by the beam in its deflected position. The “static moment” of an area is the product of that area and the distance from the area’s centroid to the given point. (It can also be expressed as an integral, as we’ll see in a minute). This sort of thing is easier to describe with examples than with words.

The first theorem should be pretty obvious, given that we already know that

y=MEI

so

y(xA)y(xB)=xAxBMEIdx

The second theorem isn’t so obvious, but if we write it in terms of an integral,

xAxBMEI(xxb)dx=xBxAy(xxb)dx

you can get to the second moment-area theorem through integration by parts.

But we’re not here to derive the moment-area method—I did that last year—we’re here to use it. For our problem, let’s sketch the deflected shape of the beam and note that the slope at its center must be zero by symmetry.

Beam deflection with tangent line

We’ll call the center A and the support at the left end B. The deflection at A is equal to the deflection of B relative to the horizontal tangent at A, which we’ll call δ.

By the second moment-area theorem, we take the moment diagram, divide it by EI, and multiply the area under the left half of the curve by the distance from B to the centroid of that area.

Moment diagram with shaded half

The area under the shaded part of the parabola is ⅔ of the area of the enclosing rectangle,

23(wL28EI)(L2)=wL324EI

The distance from B to the centroid of the shaded area is ⅝ of the distance from B to A, so

δ=(wL324EI)(58L2)=5wL4384EI

which is the answer we were expecting.

How did I know the area under a parabola and the location of its centroid without looking it up? Lots and lots of practice doing problems like this 45 years ago. There’s a diagram with all the values in a post I wrote a few years ago.

One last thing: Because the slope at the center of the beam is zero, the slope at the left end should be equal to the shaded area calculated above. Click back to the previous article in this series, and you’ll see that it does match the value of θ0 calculated there.


  1. Commonly known as “the MIT book” because that’s where the two senior authors taught. 


Simply supported beam—fourth-order ODE

In the introductory post to this series, I mentioned that the shear is the derivative of the moment. It’s easy to see that for our specific problem, the simply supported beam with a uniform distributed load,

Shear and moment diagrams

but it’s also true in general. A further relationship is that the distributed load function—let’s call it q—is the derivative of the shear (with a minus sign according to the usual sign convention). Therefore

V=dMdxandq=dVdx

which can be combined to give

q=d2Mdx2

Again, for our specific problem with a uniform load, q is just the constant value w, and it’s easy to see that the slope of the shear diagram is w.

In yesterday’s post, we used the differential relationship between the moment and displacement,

M=EId2ydx2

and the fact that we already knew that the moment was a parabola to get a solution quickly. In general, though, people usually put the last two equations together to get

EId4ydx4=q

This is a non-homogeneous fourth-order linear differential equation with constant coefficients. The solution will have four constants of integration, which can be written in terms of the initial conditions, i.e., the displacement, slope, moment, and shear at x=0:

y=y0+θ0xM02EIx2V06EIx3+yp

(For small displacements, the slope and angle are the same, which is why the initial slope is called θ0.)

The final term, yp, is the particular solution to the original equation, i.e., four integrations of q. Since q in our problem is just the constant w,

yp=w24EIx4

We know three of the initial conditions right off the bat:

y0=0,M0=0,V0=wL2

Therefore,

y=θ0xwL12EIx3+w24EIx4

We solve for the fourth initial condition by noting that the displacement at x=L is zero:

y(L)=θ0LwL412EI+wL424EI=0

Solving for θ0 gives us

θ0=wL324EI

Using this, the displacement at the center of the beam is

y(L2)=wL4EI(148196+1384)=5wL4384EI

which is the formula we were looking for.


Simply supported beam—second-order ODE

Here’s the first of the derivations for the center deflection of a simply supported beam with a uniform load.

Beam and moment diagram

We start with the differential relationship between the bending moment, M, and the deflection, y:

M=EId2ydx2

The second derivative of y is the curvature of the beam (for small deflections, which is one of the fundamental assumptions of beam theory), and the negative sign is there to account for the usual sign conventions for moment and displacement.

M is a parabola that passes through 0 at each end of the beam and peaks at wL2/8 at the center. Its formula is

M=wL2xw2x2

Therefore,

y=wEI(12x2L2x)

where I’ve started using primes for differentiation.

Integrating once gives

y=wEI(16x314Lx2)+C1

Symmetry tells us the slope at the center of the beam is zero, so

y(L2)=wEI(148L3116L3)+C1=0

Which means

C1=wL324EI

Plugging this result in and integrating again gives us

y=wEI(124x4112Lx3+124L3x)+C2

Because the deflection is zero at x=0,

C2=0

and the deflection at the center of the beam is

y(L2)=wL4EI(1384196+148)=5wL4384EI

which is the answer we were expecting.


Areas of my expertise

A few years ago, I wrote a post describing how I asked ChatGPT to solve a couple of elementary beam bending problems and how its answers were persistently wrong, even after I told it the mistakes it had made. For the first problem, determining the deflection at the center of a simply supported beam with a uniform load, ChatGPT gave the correct formula—presumably because the correct formula was part of its training corpus—but couldn’t come up with the correct numerical solution. As I said in the post:

Strictly speaking, this wasn’t a good example of a structural analysis homework problem. Students don’t get asked to just look up formulas and plug in numbers. More likely, they’d be asked to derive the equation that ChatGPT started with by either solving the differential equation for beam deflection or using some simplified technique like the moment-area or conjugate beam method. I didn’t think asking ChatGPT to do something like that would be fair.

This got me wondering how many ways I could derive the formula. A handful of ways came to me immediately, and I kept thinking of other methods over the course of the next several weeks.

Here’s a sketch of the problem:

Simply supported beam with uniform load

where w is the intensity of the load, in units of force per length, L is the length of the beam, E is the modulus of elasticity of the beam’s material, and I is the moment of inertia of the beam’s cross-section. I’m not going to get into the details of these terms or the assumptions implicit in my derivations. Suffice it to say that I’m using the typical definitions and assumptions described in strength of materials and structural analysis textbooks.

I gave myself some rules for the derivations:

I scratched out the derivations in my notebook, eventually coming up with twelve ways. They were:

  1. Second-order differential equation
  2. Fourth-order differential equation
  3. The moment-area method
  4. The conjugate beam method
  5. The slope-deflection method
  6. The “myosotis” method
  7. Energy minimization with polynomials
  8. Energy minimization with Fourier series
  9. Castigliano’s second method
  10. Finite element analysis
  11. The dummy unit load method
  12. Newmark’s method

I thought about presenting the derivations here, but I dithered over the best way to organize them. Eventually, other parts of my life intruded, and I gave up on the idea. It wasn’t until I wrote about the definition of “kip” a couple of weeks ago that I decided to just do a brain dump of all the derivations, one post for each. That’s what you’ll see here for the next couple of weeks. I know most of you don’t care about this sort of stuff, but I don’t care that you don’t care. Forewarned is forearmed—you’ll know what each post is about from their titles and can skip as you see fit.

For those few who are interested, this post will serve as a table of contents. The items in the list above will be turned into links as the posts are written.

Let me put a couple of things here that will be common. First, the deflection at the center is

5wL4384EI

This is the formula each post is aiming towards.

Second, the upward reaction forces at each end of the beam are

wL2

which is, as you might expect, half the total applied load.

Third, the shear and moment diagrams for the beam are

Shear and moment diagrams

The moment diagram is very important to many of the derivations. It’s a parabola with a peak value of

wL28

We won’t be using the shear diagram1 directly in any of the derivations, but I tend to draw it whenever I draw a moment diagram. The mathematically inclined might notice that the shear is the derivative of the moment. It passes through zero when the moment is at its peak.

OK, that’s the setup. We’ll start zipping through the derivations next time.


  1. Shear is usually denoted V because it’s a vertical force in most beams.