[Notes on Diffy Qs home] [PDF version] [Buy cheap paperback version]

[next] [prev] [prev-tail] [tail] [up]

Note: 2 lectures, §3.4 in [EP], §3.7 in [BD]

Let us look at some applications of linear second order constant coefficient equations.

Our first example is a mass on a spring. Suppose we have a mass (in kilograms) connected by a spring with spring constant (in newtons per meter) to a fixed wall. There may be some external force (in newtons) acting on the mass. Finally, there is some friction measured by (in newton-seconds per meter) as the mass slides along the floor (or perhaps there is a damper connected).

Let be the displacement of the mass ( is the rest position), with growing to the right (away from the wall). The force exerted by the spring is proportional to the compression of the spring by Hooke’s law. Therefore, it is in the negative direction. Similarly the amount of force exerted by friction is proportional to the velocity of the mass. By Newton’s second law we know that force equals mass times acceleration and hence or

This is a linear second order constant coefficient ODE. We set up some terminology about this equation. We say the motion is

- (i)
- forced, if (if is not identically zero),
- (ii)
- unforced or free, if (if is identically zero),
- (iii)
- damped, if , and
- (iv)
- undamped, if .

This system appears in lots of applications even if it does not at first seem like it. Many real world scenarios can be simplified to a mass on a spring. For example, a bungee jump setup is essentially a mass and spring system (you are the mass). It would be good if someone did the math before you jump off the bridge, right? Let us give two other examples.

Here is an example for electrical engineers. Suppose that we have the pictured RLC circuit. There is a resistor with a resistance of ohms, an inductor with an inductance of henries, and a capacitor with a capacitance of farads. There is also an electric source (such as a battery) giving a voltage of volts at time (measured in seconds). Let be the charge in coulombs on the capacitor and be the current in the circuit. The relation between the two is . By elementary principles we have that . If we differentiate we get

This is an nonhomogeneous second order constant coefficient linear equation. Further, as , and are all positive, this system behaves just like the mass and spring system. The position of the mass is replaced by the current. Mass is replaced by the inductance, damping is replaced by resistance and the spring constant is replaced by one over the capacitance. The change in voltage becomes the forcing function. Hence for constant voltage this is an unforced motion.

Our next example behaves like a mass and spring system only approximately. Suppose we have a mass on a pendulum of length . We wish to find an equation for the angle . Let be the force of gravity. Elementary physics mandates that the equation is of the form

Let us derive this equation using Newton’s second law; force equals mass times acceleration. The acceleration is and mass is . So has to be equal to the tangential component of the force given by the gravity, that is in the opposite direction. So . The curiously cancels from the equation.

Now we make our approximation. For small we have that approximately . This can be seen by looking at the graph. In Figure 2.1 we can see that for approximately (in radians) the graphs of and are almost the same.

Therefore, when the swings are small, is always small and we can model the behavior by the simpler linear equation

Note that the errors that we get from the approximation build up. So after a very long time, the behavior of the real system might be substantially different from our solution. Also we will see that in a mass-spring system, the amplitude is independent of the period. This is not true for a pendulum. Nevertheless, for reasonably short periods of time and small swings (for example if the pendulum is very long), the approximation is reasonably good.

In real world problems it is often necessary to make these types of simplifications. We must understand both the mathematics and the physics of the situation to see if the simplification is valid in the context of the questions we are trying to answer.

In this section we will only consider free or unforced motion, as we cannot yet solve nonhomogeneous equations. Let us start with undamped motion where . We have the equation

If we divide by and let , then we can write the equation as

The general solution to this equation is

By a trigonometric identity, we have that for two different constants and , we have

It is not hard to compute that and . Therefore, we let and be our arbitrary constants and write .

Exercise 2.4.1: Justify the above identity and verify the equations for and . Hint: Start with and multiply by . Then think what should and be.

While it is generally easier to use the first form with and to solve for the initial conditions, the second form is much more natural. The constants and have very nice interpretation. We look at the form of the solution

We can see that the amplitude is , is the (angular) frequency, and is the so-called phase shift. The phase shift just shifts the graph left or right. We call the natural (angular) frequency. This entire setup is usually called simple harmonic motion.

Let us pause to explain the word angular before the word frequency. The units of are radians per unit time, not cycles per unit time as is the usual measure of frequency. Because one cycle is radians, the usual frequency is given by . It is simply a matter of where we put the constant , and that is a matter of taste.

The period of the motion is one over the frequency (in cycles per unit time) and hence . That is the amount of time it takes to complete one full cycle.

Example 2.4.1: Suppose that and . The whole mass and spring setup is sitting on a truck that was traveling at 1 ^{m}/_{s}. The truck crashes and hence stops. The mass was held in place 0.5 meters forward from the rest position. During the crash the mass gets loose. That is, the mass is now moving forward at 1 ^{m}/_{s}, while the other end of the spring is held in place. The mass therefore starts oscillating. What is the frequency of the resulting oscillation and what is the amplitude. The units are the mks units (meters-kilograms-seconds).

The setup means that the mass was at half a meter in the positive direction during the crash and relative to the wall the spring is mounted to, the mass was moving forward (in the positive direction) at 1 ^{m}/_{s}. This gives us the initial conditions.

So the equation with initial conditions is

We can directly compute . Hence the angular frequency is 2. The usual frequency in Hertz (cycles per second) is .

The general solution is

Letting means . Then . Letting we get . Therefore, the amplitude is . The solution is

A plot of is shown in Figure 2.2.

In general, for free undamped motion, a solution of the form

corresponds to the initial conditions and . Therefore, it is easy to figure out and from the initial conditions. The amplitude and the phase shift can then be computed from and . In the example, we have already found the amplitude . Let us compute the phase shift. We know that . We take the arctangent of 1 and get approximately 0.785. We still need to check if this is in the correct quadrant (and add to if it is not). Since both and are positive, then should be in the first quadrant, and 0.785 radians really is in the first quadrant.

Note: Many calculators and computer software do not only have the atan function for arctangent, but also what is sometimes called atan2. This function takes two arguments, and , and returns a in the correct quadrant for you.

Let us now focus on damped motion. Let us rewrite the equation

as

where

The characteristic equation is

Using the quadratic formula we get that the roots are

The form of the solution depends on whether we get complex or real roots. We get real roots if and only if the following number is nonnegative:

The sign of is the same as the sign of . Thus we get real roots if and only if is nonnegative, or in other words if .

When , we say the system is overdamped. In this case, there are two distinct real roots and . Both roots are negative: As is always less than , then is negative in either case.

The solution is

Since are negative, as . Thus the mass will tend towards the rest position as time goes to infinity. For a few sample plots for different initial conditions, see Figure 2.3.

Do note that no oscillation happens. In fact, the graph will cross the axis at most once. To see why, we try to solve . Therefore, and using laws of exponents we obtain

This equation has at most one solution . For some initial conditions the graph will never cross the axis, as is evident from the sample graphs.

Example 2.4.2: Suppose the mass is released from rest. That is and . Then

It is not hard to see that this satisfies the initial conditions.

When , we say the system is critically damped. In this case, there is one root of multiplicity 2 and this root is . Our solution is

The behavior of a critically damped system is very similar to an overdamped system. After all a critically damped system is in some sense a limit of overdamped systems. Since these equations are really only an approximation to the real world, in reality we are never critically damped, it is a place we can only reach in theory. We are always a little bit underdamped or a little bit overdamped. It is better not to dwell on critical damping.

When , we say the system is underdamped. In this case, the roots are complex.

where . Our solution is

or

An example plot is given in Figure 2.4. Note that we still have that as .

In the figure we also show the envelope curves and . The solution is the oscillating line between the two envelope curves. The envelope curves give the maximum amplitude of the oscillation at any given point in time. For example if you are bungee jumping, you are really interested in computing the envelope curve so that you do not hit the concrete with your head.

The phase shift just shifts the graph left or right but within the envelope curves (the envelope curves do not change if changes).

Finally note that the angular pseudo-frequency (we do not call it a frequency since the solution is not really a periodic function) becomes smaller when the damping (and hence ) becomes larger. This makes sense. When we change the damping just a little bit, we do not expect the behavior of the solution to change dramatically. If we keep making larger, then at some point the solution should start looking like the solution for critical damping or overdamping, where no oscillation happens. So if approaches , we want to approach 0.

On the other hand when becomes smaller, approaches ( is always smaller than ), and the solution looks more and more like the steady periodic motion of the undamped case. The envelope curves become flatter and flatter as (and hence ) goes to 0.

Exercise 2.4.2: Consider a mass and spring system with a mass , spring constant , and damping constant . a) Set up and find the general solution of the system. b) Is the system underdamped, overdamped or critically damped? c) If the system is not critically damped, find a that makes the system critically damped.

Exercise 2.4.3: Do Exercise 2.4.2 for , , and .

Exercise 2.4.4: Using the mks units (meters-kilograms-seconds), suppose you have a spring with spring constant 4 ^{N}/_{m}. You want to use it to weigh items. Assume no friction. You place the mass on the spring and put it in motion. a) You count and find that the frequency is 0.8 Hz (cycles per second). What is the mass? b) Find a formula for the mass given the frequency in Hz.

Exercise 2.4.5: Suppose we add possible friction to Exercise 2.4.4. Further, suppose you do not know the spring constant, but you have two reference weights 1 kg and 2 kg to calibrate your setup. You put each in motion on your spring and measure the frequency. For the 1 kg weight you measured 1.1 Hz, for the 2 kg weight you measured 0.8 Hz. a) Find (spring constant) and (damping constant). b) Find a formula for the mass in terms of the frequency in Hz. Note that there may be more than one possible mass for a given frequency. c) For an unknown object you measured 0.2 Hz, what is the mass of the object? Suppose that you know that the mass of the unknown object is more than a kilogram.

Exercise 2.4.6: Suppose you wish to measure the friction a mass of 0.1 kg experiences as it slides along a floor (you wish to find ). You have a spring with spring constant . You take the spring, you attach it to the mass and fix it to a wall. Then you pull on the spring and let the mass go. You find that the mass oscillates with frequency 1 Hz. What is the friction?

Exercise 2.4.101: A mass of kilograms is on a spring with spring constant newtons per meter with no damping. Suppose the system is at rest and at time the mass is kicked and starts traveling at 2 meters per second. How large does have to be to so that the mass does not go further than 3 meters from the rest position?

Exercise 2.4.102: Suppose we have an RLC circuit with a resistor of 100 miliohms (0.1 ohms), inductor of inductance of 50 millihenries (0.05 henries), and a capacitor of 5 farads, with constant voltage. a) Set up the ODE equation for the current . b) Find the general solution. c) Solve for and .

Exercise 2.4.103: A 5000 kg railcar hits a bumper (a spring) at 1 ^{m}/_{s}, and the spring compresses by 0.1 m. Assume no damping. a) Find . b) Find out how far does the spring compress when a 10000 kg railcar hits the spring at the same speed. c) If the spring would break if it compresses further than 0.3 m, what is the maximum mass of a railcar that can hit it at 1 ^{m}/_{s}? d) What is the maximum mass of a railcar that can hit the spring without breaking at 2 ^{m}/_{s}?