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. 


Unsound

I mentioned last week that some wording in a news article struck me as odd. A similar thing happened when I read this Scientific American piece (that’s an Apple News link—I don’t have a SciAm subscription). In this case, I didn’t even have to read the article; the oddity is right there in the headline: “A SpaceX rocket booster is on track to hit the Moon at several times the speed of sound.”

I know what the writer means, of course. He’s telling us that the SpaceX debris will hit the Moon at several thousand feet per second, 1,000 fps being about the speed of sound we’re used to here on Earth. But the Moon isn’t the Earth, so it’s kind of a weird comparison, don’t you think? It’s not as if the booster is going to create a sonic boom. I’d find this less odd if the story were in a normal newspaper or magazine instead of Scientific American.

Again, I realize that the writer is just giving his readers a point of comparison, but is the speed of sound (on Earth) really something most people have an intuitive feel for? I remember learning as a kid to count off the seconds between a lightning flash and the resulting thunder to get the distance to the lightning strike in thousands of feet, but that doesn’t mean I have a strong sense of the speed of sound. Other than “really fast,” which isn’t all that helpful.

It’s mentioned in the body of the article that the speed of the booster at impact, which will be on August 5, is estimated to be about 5,400 mph. That’s about 100 times the posted speed limit on a two-lane highway, which is something most Americans do have a feel for.

By the way, if you have any interest in the problem of space debris coming back and striking the Earth (where the speed of sound might be relevant), you should follow Sam Lawler on Mastodon. She’s a professor of astronomy at the University of Regina and also has a keen interest in the clogging of near-Earth orbital space by the huge number of satellites launched in recent years. Her Mastodon feed is also a great source of farm animal photos, most recently a crop of baby goats.


A Key followup

Earlier this week, I saw this article in Apple News. It’s from the San Francisco Chronicle and discusses a recent report on the safety of the Golden Gate Bridge. The report is one of several reports spurred by the Francis Scott Key Bridge collapse a couple of years ago.

Golden Gate Bridge from Wikipedia

Image from Wikipedia.

Spoiler: the report finds the Golden Gate Bridge safe—quite unlikely to suffer damage from the impact of a ship. This has to do with the Golden Gate’s two towers:

A striking omission from the Chronicle story is a link to the report itself. But that appears to be the fault of Apple News. I can’t read the story on the Chronicle’s website because I’m not a subscriber, but this reprint on Yahoo! News includes the report embedded as a PDF and available to download. So it looks like Apple removed the report itself, which is pretty poor form.

The nice thing about finding the actual report, written by HDR, Inc., is that it confirmed some suspicions I had regarding this paragraph in the Chronicle story:

The resilience of the Golden Gate Bridge partly comes from sheer strength. The south tower, on the San Francisco side, is described in the report as a “robust structural feature like no other” and is surrounded by a reinforced concrete protective shell up to 28 feet thick. It can withstand about 50,000 kips of force, or roughly 25,000 tons. In many cases, engineers found, a ship would crumple and absorb its own impact energy before it could seriously damage the structure.

First, it’s not “sheer strength,” it’s “shear strength.” HDR calculated the strength of the protective structure around the south tower to be at least 50,000 kips when that force is attempting to shear through the reinforced concrete wall. They specifically mention shearing capacity, shearing interfaces, and shearing area when discussing this calculation on pp. 59–60 of the report. While it’s true that engineers tend to be crummy writers, we definitely know the difference between “shear” and “sheer.”

Second, I suspect the writer, Brooke Park, doesn’t know what a kip is. When writing for a general audience1 most people wouldn’t use the word “kip” without saying it’s short for “kilopound.” Which is to say,

1kip=1000lbs=12ton

There’s no “roughly” about it.2

Still, this oddball paragraph doesn’t affect the overall story and provides the engineers who read it some amusement. Too bad about Apple’s redaction of a link to the report itself, though. That’s sheer incompetence.


  1. Which HDR isn’t, so I don’t blame them for using structural engineering terms without further explanation. 

  2. You could argue that Park’s “roughly” was meant to parallel the “about” in the previous clause. I think that’s an overly generous interpretation. Also, don’t write to me about long tons or metric tonnes—there’s no way Park was talking about those.