Pendulum, let it go

This is a requested post. Clearly, I do requests. The idea here is that I am going to give all the details needed to determine the equation of motion (and then model it) for a basic pendulum. Warning: this post is a little more advanced than my normal posts. There are some prerequisites. You need to understand derivatives. I will assume that you do. Here is a pendulum. (and this time I will stick to my variables)

i-bececf3a25249f4632112bba54ff76c2-2009-12-08_untitled_94.jpg

Like I said before, this is a tricky problem unless I use some tricks. The problem is that the tension the string exerts on the mass changes. Here is my trick: think about a coordinate system that moves with the mass.

i-15cbd32b79cae39b565b8834c0011c70-2009-12-08_untitled.jpg

The mass only moves in the direction of the theta axis. So, I only care about forces in this direction. The tension from the string is always perpendicular to the direction of motion. There is a theta-direction-component of the gravitational force. It is:

i-7d951687043e1a58ec0f2bd9d57f3f84-2009-12-08_la_te_xi_t_1_7.jpg

Now I need the acceleration in the theta-direction. This would be related to the second derivative with respect to time of the angle:

i-4f3b156d0f1ec8eb6fbec6ecdcd69fe2-2009-12-08_la_te_xi_t_1_8.jpg

This is using the common relationship between angular and linear quantities for something moving in a circle:

i-32d584f5410aa7f0a804c939b095f893-2009-12-08_la_te_xi_t_1_9.jpg

So, now I can put this together in Newton's second law:

i-38be87fca56405078b9f924c28518e20-2009-12-08_la_te_xi_t_1_10.jpg

And the masses cancel (motion of this type of pendulum does not depend on mass). This leaves the following:

i-cb902e365a58cf1a289f7675b331dc5d-2009-12-08_la_te_xi_t_1_11.jpg

Looks good. I have a differential equation relating theta and time. I should be all set. However, this really isn't a very easy equation to solve. So, the trick is to only look for the cases when theta is small. Here is a plot of the sine of theta as a function of theta.

i-bfe7432382f1037c6719f1571f7cccad-2009-12-08_sintheta.jpg

Actually, the blue line is sine of theta and the red line is theta = theta. For theta less than 0.4 radians (22 degrees), these two functions are very similar. So, for that case, I can write the equation as:

i-6bdc0fec60d3fe276c443592eff6b99c-2009-12-08_la_te_xi_t_1_12.jpg

Now, this is a differential equation that I can solve. If you like, you can make this a homework problem for a diff-eq class. What method should I use to solve this differential equation? I will use the one I always use - guessing. Really, this is legit. If I can guess a solution and that solution works, I am finished. What function when I take the derivative with respect to time twice, do I get the same function back (with a negative constant)? There are two that easily work (there are actually more than two). Take a look at these two functions:

i-4dddf73353a79aa402a27459b6d01c6d-2009-12-08_la_te_xi_t_1_13.jpg

(I know I could add a phase - but I am not going to do that) Actually, if each of these are solutions, then the sum of these two are a solution.

i-4bbfb291c97288dba8fe5c10d53608f3-2009-12-08_la_te_xi_t_1_14.jpg

Let me show you that this is indeed a solution by taking the derivative (with respect to time) twice.

i-81ef23bef86b60e0d57c8a48481e6c1b-2009-12-08_la_te_xi_t_1_18.jpg

The only way this can be a solution is if:

i-277e784983cee1789be9060ac26cecfc-2009-12-08_la_te_xi_t_1_19.jpg

So, it is finished. If the angle is small, then the motion is sinusoidal with an angular frequency the depends on g and the length (which is your traditional textbook answer - except maybe they use L instead of R).

Oh wait. I just realized that I never solved for A and B. These depend on the initial conditions. I can uniquely define the initial conditions if I know the initial angle and the initial angular velocity. So, at t = 0 seconds:

i-69f58a805c7bc5b97a4c6a841948b95b-2009-12-09_la_te_xi_t_1_1.jpg

I know what you are thinking - but what if the angle is not small? Then I can go back to the original equation for which there is no simple solution. I can easily create a numerical solution for this (in a spreadsheet or python or something). For this case, I will use the Euler Method to solve for it. The basic idea is to break the problem into tiny time steps. During each step, I can calculate the angular acceleration (second derivative with respect to time of the angle) using the solution above (for the very first calculation, I can use the initial conditions)

i-5ba1bf57296ab1b8901db8261307c573-2009-12-08_la_te_xi_t_1_20.jpg

Now, during this time interval, the following is true regarding the rate of change of the angle and the second derivative of the rate of change. (I am using dot-notation where 1 dot means derivative with respect to time and two dots means a second time derivative).

i-7c0965cc23e7ee6660a7829165c5214c-2009-12-08_la_te_xi_t_1_21.jpg

So, if my time interval is small, I can pretend like the theta-double dot does not change during this interval (basically true). Then if I know one theta-dot, I can find the next one.

i-28df328045c11335c425812ce2e8c629-2009-12-08_la_te_xi_t_1_22.jpg

I can use this same trick to find theta.

i-6df901b48010ffbd5aaf06dd8b623c78-2009-12-08_la_te_xi_t_1_23.jpg

Yes, I know there are more elegant ways to do this, but my computer is fast enough to do it the rough way. If I just keep doing this in small steps, I can find the answer. Normally, I would do this in python (because it is awesome), but in this case I will do it in a spreadsheet. Here it is (feel free to play with it).

Now I am ready to put all of this in a spread sheet.

A couple of notes:

  • I also plotted the solution from the small angle approximation - so you could get a comparison
  • Apparently, google docs does not like plotting data in non-adjacent columns, so I put the small angle calculation right next to theta calculation
  • I also calculated x and y for the mass, but I didn't use that
  • I put dt as a small number so that the data would look ok, it should probably be a little smaller.
  • My angles are in radians

In case you don't want to play around with the spreadsheet, here is a plot of the two solutions for an initial angle of pi/4.

More like this

neat... if you do requests, I have a practical one I need answered.

can you show me a way to calculate forces on tightropes of varying degrees of tautness? trying to set up a tightrope/slackrope frame and would like to be able to calculate the forces of the rope on the frame for a given weight of rope walker.

Thanks for the refresher. The important thing that you derived is that the period is independent of the mass and proportional to the square root of the length. Practical application: I bought a really neat wind chime (gong actually). In moderate winds it didn't chime. But I noticed that it was a coupled pendulum (a more involved analysis case) with roughly equal lengths. In moderate winds the large gong and the knocker swung in unison. So I changed the lengths of one of the coupled pendulums by a few percent and now it chimes like a bandit in relatively mild wind.

@Rhett:
Nice article.

The small-angle version is a friendlier version of the usual textbook treatment, and I thought the numerical integration was really nice. It may be all the new modern technology (free online spreadsheets, fast computers) but that's the least intimidating example of numerical integration I've ever seen.

By Anonymous Coward (not verified) on 09 Dec 2009 #permalink

Thanks for taking up my request! I thought the whole large angle situation was impossible. But numerical solutions, while I like it less, works, so I am happy.