There’s nothing you and I won’t do

I’ve talked about area moments of inertia in the past couple of posts. Today, thanks to Phil Plait, resident astronomer at Scientific American, I have an excuse to talk about mass moments of inertia, the kind you learned about in first semester physics class.

Plait’s article today takes on the Modern English song “I Melt with You” and looks into the energy it would take to stop the world.1 By “stopping,” he initially assumes that means stopping the rotation of the Earth about its axis. He gives the value, 2×1029joules, but doesn’t explain how he calculated that number. Kind of disappointing for SciAm and definitely not acceptable here at ANIAT. Especially since the calculation is easy.

The kinetic energy of a rotating body is

12Iω2

where I is the mass moment of inertia of the body about its axis of rotation and ω is its rotational speed.

The value of I is kind of tricky because, as Plait points out, the Earth’s density varies with depth. But we can make an initial estimate by assuming the Earth is a sphere of uniform density. The moment of inertia for that kind of body is

I=25mr2

We get the mass and (average) radius of the Earth from the sadly defunct NASA Earth Fact Sheet, currently available through the Internet Archive. The mass is 5.9722×1024kg and the mean radius is 6371km. That gives us a moment of inertia of

I=25(5.9722×1024kg)(6.371×106m)2=9.6964×1037kgm2

Relative to the stars, the Earth rotates once every 23.9345 hours (that’s the “sidereal rotation period” on the Fact Sheet), so the rotational speed is

ω=2π(23.9345×60×60)s=7.29211×105s1

Putting these together, we get a kinetic energy of

12(9.6964×1037kgm2)(7.29211×105s1)2=2.578×1029J

This is higher than Plait’s value, but he may have rounded down because he knew the density of the Earth is higher near the center. Also, I suspect he was mostly after an order-of-magnitude answer.

I found a better estimate of I on this page from Eric Weisstein’s World of Science site. It was determined by Kurt Lambeck to be

I=8.034×1037kgm2

(Lambeck’s book came out just a couple of years before “I Melt with You,” which I’m pretty sure is a coincidence.)

Using this value, we get a rotational kinetic energy of

12(8.034×1037kgm2)(7.29211×105s1)2=2.136×1029J

which is a better match for Plait’s energy value.

Plait also considers another interpretation of “I’ll stop the world”: stopping the Earth in its orbit. I won’t get into the details of this, first because it doesn’t involve the moment of inertia, and second because stopping the Earth this way would have it fall into the Sun. Not much doubt about melting in that case.


  1. He also gets into the question of melting the Earth’s crust, which I won’t be addressing. After all, the song says “I’ll stop the world and melt with you.” As far as I’m concerned, only the singer and his girl melt, not the whole goddamned world. 


Axial vs. bending stiffness

At the end of yesterday’s post, I mentioned that axial deformations in framework members are typically two orders of magnitude less than the lateral (bending) deformations. It’s easy to see why.

Consider a column with length L, cross-sectional area A, and moment of inertia I. We’ll calculate the tip deflection of this column under two loads of the same magnitude, P: an axial compression load and a horizontal load to the right.

Column with axial and lateral loads

Under the compression load, the tip of the column will move down by

PLEA

Under the horizontal load, the tip will move to the right by

PL33EI

The ratio of lateral deflection to axial deflection is therefore

13AL2I

When calculating buckling loads of columns, it’s often useful to consider the radius of gyration, r, which is defined this way:

I=Ar2r=IA

We’re not doing a buckling calculation here, but we can use the radius of gyration to help us get a sense of the deflection ratio, which we can rewrite as

13(Lr)2

The radius of gyration is some fraction of a characteristic dimension of the cross-section. For example, if the cross-section is a square with side length b, the area is b2 and the moment of inertia is b4/12, so

r=b4/12b2=b230.289b

In typical structures, the stick-like members—beams, columns, braces, etc.—tend to be long and thin, with aspect ratios of about 10. The deflection ratio is roughly the square of the aspect ratio, so lateral deflections are typically two orders of magnitude greater than axial deflections. That’s why it’s common to ignore the axial deflections in analyses like yesterday’s.

This doesn’t mean axial deflections can always be ignored. In tall buildings, the columns carry the weight of all the floors above them, and column shortening has to be accounted for, especially in the lower floors. Also, when using modern numerical techniques for structural analysis, there’s basically no extra work involved in calculating the axial deflections along with the bending deflections, so the axial deflections pop out of the computer as a matter of course.


Frame sidesway again

I’ve been meaning to follow up on the “sense of structure” post I wrote just after Thanksgiving. I don’t have anything profound to say, but I did want to give complete solutions for both the off-center vertical load problem and the horizontal load problem.

First, though, I want to mention a nice email conversation I had with Karl Hoitsma, proprietor of the Mathpax blog. He pointed me to this post of his from a few years ago in which he describes being asked essentially the same frame sidesway question I was. The only real differences were that his oral exam was about five years after mine and that he was at Texas A&M instead of U of I. Oh, and he hadn’t been prepped for the question.

In that post, Karl gives a few solutions to the problem, which I found fun to go through. I had already solved the horizontal load problem for a simpler frame, and it was interesting to compare our approaches to the problem. As I’ve been trying to get better at Mathematica, I decided to do the off-center load problem and redo the horizontal load problem using the Wolfram Language. I won’t go through the code here, just the results.

Let’s start with the simpler problem: a portal frame with fixed column bases in which the beam and columns are identical in both cross-section and length.

Side load

E is the modulus of elasticity of the material, and I is the moment of inertia of the cross-section.1

This is a statically indeterminate system. You cannot determine the reactions at the bases of the columns through statics alone; you have to account for the stiffnesses of the members. But once you’ve done that, the reactions are these:

Side load reactions

It’s a decent exercise to confirm that these reactions are in equilibrium. It should be immediately obvious that the horizontal and vertical forces are in balance. It takes only a little more effort to show that the clockwise and counter-clockwise moments about any point are also in balance. (As a practical matter, I’d choose the base of one of the columns to take moments about, as that eliminates three of the forces from the equation.)

What I really wanted to show was the deflected shape of the frame under this loading. If you recall from that earlier post, I drew a deflected shape in which the left half of the beam goes up and the right half goes down. Knowing this was key to using reciprocity to determine which way the frame sways when it’s subjected to an off-center vertical load. Here’s the deflected shape plotted from a formal solution instead of just drawn on the basis of my sense of the structure’s behavior:

Side load deflection

The leftward displacement at the top of each column is

Δ=5PL384EI

and the counter-clockwise rotation at the top of each column is

θ=PL228EI

It’s the CCW rotation that lifts the left half of the beam upward and pushes the right half downward. As you can see, though, the vertical deflection of the beam is pretty small compared to the horizontal deflection of the columns.

A note on units: In these equations, P has units of force, L has units of length, E has units of force per length squared, and I has units of length to the fourth power. So Δ has units of length, as expected, and θ is a pure number, which means we’re getting the rotation in radians.

OK, now let’s look at the frame with a downward off-center load:

Off-center load

I’ve put the load at the quarter point of the beam for convenience, but it doesn’t matter where it is as long as it’s to the left of the centerline. Here are the reactions:

Off-center load reactions

Again, you can confirm that these forces and moments are in equilibrium. Here’s the deflected shape:

Off-center load deflection

For this loading, the leftward displacement at the top of each column is

Δ=3PL3896EI

The minus sign means the displacement is actually to the right, which matches what we figured out a few weeks ago. Now we have the magnitude in addition to the direction.

As I said above, I don’t want to go through the details of the code that got me to these answers, but if you’re interested, you can see it using the following Wolfram Cloud links:

I should mention that in both of these analyses I’m ignoring the axial deformation of the members. This is a common practice in frame analysis, as axial deformations are typically about two orders of magnitude less than lateral deformations. Leaving out the axial terms cuts out a lot of complexity in symbolic solutions like this while retaining good accuracy.


  1. In physics class, you used a different moment of inertia when you were learning about rotating bodies. Strictly speaking, the I in this problem should be called the second moment of area, but structural engineers always call it the moment of inertia. It’s a reasonable reuse of the phrase, as the second moment of area is a measure of the resistance of a cross-section to rotation as a structural member is bent. 


Be careful where you click

I ran into an example of especially bad design on a web page yesterday. Not “how it looks” design; “how it works” design.

Here’s an excerpt from an online payment page:

Excerpt from payment web page

Looks straightforward, right? Just click in the field, type in the amount to pay, and off you go. But unlike almost every other example you’ll see nowadays, this field doesn’t automatically select its dummy contents when you click in it. So what you type does not replace the dummy contents. Here’s a video showing what can happen.

There are three interactions in the video.

  1. I clicked to the right of the “$0.00” default amount and typed “12345.” You can see the zeros disappear as I type, and the amount at the end is $123.45. Let’s assume this is what I wanted to enter.
  2. I clicked just to the right of the decimal point and typed “12345.” The cursor first jumps to the beginning of the amount. After I type the “1,” it jumps to the end, and the “2345” fills in from the right. But two of the zeros never disappear, and the amount at the end is $10,023.45.
  3. I clicked just to the right of the dollar sign and typed “12345.” After I type “123,” the cursor jumps to the end and the “45” fills in from the right. None of the zeros disappear, and the amount at the end is $123,000.45.

I want to emphasize that I never moved the cursor myself after the initial click in the field. All I did was type.

As you might guess, I discovered this when I nearly made an enormous overpayment.

It’s true that a standard HTML <input type="text" /> form field with a value attribute will put the cursor where you click, so a user could screw up the entry in that kind of field, too. But the value attribute is meant to be a default, a common legitimate entry. No one ever goes online to pay $0.00. A more appropriate pure HTML solution would be to use the placeholder attribute, which puts gray text in the field to suggest what the field is for. As soon as you start typing in such a field, the placeholder text disappears.

For the designers of this website, though, standard HTML just wasn’t good enough. They wanted to show off their elite JavaScript skills to keep the dollar sign and decimal point in place and add commas as the user types. In the process, they created a field in which the cursor jumps around haphazardly and the user can easily create a mistaken payment.

Am I the only user who’s noticed this? Impossible.1 But we’ve gotten so used to crap like this that we just sigh, figure out how to fix things on our end, and move on with our lives. I did point this problem out to the vendor, who was surprised to hear it, but payment processing is a service they buy from someone else. It’s unlikely that the people responsible for the poor user interface will ever be told to fix it.


  1. I probably am the only person who tried out all the different places to click. Not included in the video: clicking just to the left of the decimal point, which works the same as clicking just to right of it; and clicking between the zeros after the decimal point, which leads to an entry of $1,023.45