Under, Over and Critical Damping

Download as pdf or txt
Download as pdf or txt
You are on page 1of 6

Under, Over and Critical Damping 1.

Response to Damping
As we saw, the unforced damped harmonic oscillator has equation mx + bx + kx = 0, with m > 0, b 0 and k > 0. It has characteristic equation ms2 + bs + k = 0 with characteristic roots

..

(1)

b2 4mk 2m

(2)

There are three cases depending on the sign of the expression under the square root: i) b2 < 4mk (this will be underdamping, b is small relative to m and k). ii) b2 > 4mk (this will be overdamping, b is large relative to m and k). iii) b2 = 4mk (this will be critical damping, b is just between over and underdamping. Mathematically, the easiest case is overdamping because the roots are real. However, most people think of the oscillatory behavior of a damped oscillator. Since this is connected to underdamping we start with that case. Case (i) Underdamping (non-real complex roots) If b2 < 4mk then the term under the square root is negative and the characteristic roots are not real. In order for b2 < 4mk the damping constant b must be relatively small. First we use the roots (2) to solve equation (1). Let d = |b2 4mk|/2m. Then we have b characteristic roots: i d . leading to 2m complex exponential solutions: e(b/2m+id )t , e(b/2mid )t . The basic real solutions are ebt/2m cos(d t) and ebt/2m sin(d t). The general real solution is found by taking linear combinations of the two basic solutions, that is: x (t) = c1 ebt/2m cos(d t) + c2 ebt/2m sin(d t)

Under, Over and Critical Damping

OCW 18.03SC

or x (t) = ebt/2m (c1 cos(d t) + c2 sin(d t)) = Aebt/2m cos(d t ). (3)

Lets analyze this physically. When b = 0 the response is a sinusoid. Damping is a frictional force, so it generates heat and dissipates energy. When the damping constant b is small we would expect the system to still oscillate, but with decreasing amplitude as its energy is converted to heat. Over time it should come to rest at equilibrium. This is exactly what we see in (3). The factor cos(d t ) shows the oscillation. The exponential factor ebt/2m has a negative exponent and therefore gives the decaying amplitude. As t , the exponential goes asymptotically to 0, so x (t) also goes asympotically to its equilibrium position x = 0. We call d the damped angular (or circular) frequency of the system. This is sometimes called a pseudo-frequency of x (t). We need to be careful to call it a pseudo-frequency because x (t) is not periodic and only periodic functions have a frequency. Nonetheless, x (t) does oscillate, crossing x = 0 twice each pseudo-period. Example 1. Show that the system x + 1x + 3x = 0 is underdamped, nd its damped angular frequency and graph the solution with initial conditions . x (0) = 1, x (0) = 0. Solution. Characteristic equation: s2 + s + 3 = 0. Characteristic roots: 1/2 i 11/2. Basic real solutions: et/2 cos( 11 t/2), et/2 sin( 11 t/2). General solution: x (t) = et/2 (c1 cos( 11 t/2) + c2 sin( 11 t/2)) = Aet/2 cos( 11 t/2 ). Since the roots have nonzero imaginary part, the system is underdamped. The damped angular frequency is d = 11/2. The initial conditions are satised when c1 = 1 and c2 = 1/ 11. So, 1 t/2 x (t) = e cos( 11 t/2) + sin( 11 t/2) 11 12 = et/2 cos( 11 t/2 ), 11 where = tan1 (1/ 11).

..

Under, Over and Critical Damping

OCW 18.03SC

Figure 1: The damped oscillation for example 1. Case (ii) Overdamping (distinct real roots) If b2 > 4mk then the term under the square root is positive and the characteristic roots are real and distinct. In order for b2 > 4mk the damping constant b must be relatively large. One extremely important thing to notice is that in this case the roots are both negative. You can see this by looking at the formula (2). The term under the square root is positive by assumption, so the roots are real. Since b2 4mk < b2 the square root is less than b and therefore the root b + b2 4mk < 0. The other root is clearly negative. Now we use the roots to solve equation (1) in this case. b + b2 4mk b b2 4mk Characteristic roots: r1 = , r2 = . 2m 2m Exponential solutions: er1 t , er2 t . General solution: x ( t ) = c 1 e r1 t + c 2 e r2 t . Lets analyze this physically. When the damping is large the frictional force is so great that the system cant oscillate. It might sound odd, but an unforced overdamped harmonic oscillator does not oscillate. Since both exponents are negative every solution in this case goes asymptotically to the equilibrium x = 0. At the top of many doors is a spring to make them shut automatically. The spring is damped to control the rate at which the door closes. If the damper is strong enough, so that the spring is overdamped, then the door just settles back to the equilibrium position (i.e. the closed position) without oscillating which is usually what is wanted in this case. Example 2. Show that the system x + 4x + 3x = 0 is overdamped and . graph the solution with initial conditions x (0) = 1, x (0) = 0. Which root controls how fast the solution returns to equilibrium? 3

..

Under, Over and Critical Damping Solution. Characteristic equation: s2 + 4s + 3 = 0. Characteristic roots: (this factors) 1, 3. Exponential solutions: et , e3t . General solution: x ( t ) = c 1 e t + c 2 e 3t .

OCW 18.03SC

Because the roots are real and different, the system is overdamped. The intial conditions are satised when c1 = 3/2, c2 = 1/2. So, x (t) = 3et /2 e3t /2.

x 0.0 0 0.2 0.4

0.6

0.8

1.0

2 t

Figure 2: The overdamped graph for example 2. goes to 0 more slowly than e3t/2 it controls the rate at which x goes to 0. (Remember, it is the term that goes to zero slowest term that controls the rate.) Case (iii) Critical Damping (repeated real roots) If b2 = 4mk then the term under the square root is 0 and the characteristic polynomial has repeated roots, b/2m, b/2m. Now we use the roots to solve equation (1) in this case. We have only one exponential solution, so we need to multiply it by t to get the second solution. Basic solutions: ebt/2m , tebt/2m . General solution: x (t) = ebt/2m (c1 + c2 t). As in the overdamped case, this does not oscillate. It is worth noting that for a xed m and k, choosing b to be the critical damping value gives the fastest return of the system to its equilibrium position. In engineering design this is often a desirable property. This can be seen by considering the roots, but we will not go through the algebra that shows this. (See gure (4).) 4 Because et

Under, Over and Critical Damping

OCW 18.03SC

Example 3. Show that the system x + 4x + 4x = 0 is critically damped and . graph the solution with initial conditions x (0) = 1, x (0) = 0. Solution. Characteristic equation: s2 + 4s + 4 = 0. Characteristic roots: (this factors) 2, 2. Exponential solutions: (only one) e2t . General solution: x ( t ) = e 2t ( c 1 + c 2 t ). Because the roots are repeated, the system is critically damped. The intial conditions are satised when c1 = 1, c2 = 2. So, x ( t ) = e 2t ( 1 + 2 t ).

..

x 0.0 0 0.2 0.4

0.6

0.8

1.0

2 t

Figure 3: The critically damped graph for example 3. Notice that qualitatively the graphs for the overdamped and critically damped cases are similar. The following gure shows plots for solutions to x + bx + x = 0 with . initial conditions x (0) = 1, x (0) = 0. The three plots are b = 1 underdamped; b = 2 critically damped (dashed line); b = 3 overdamped. Notice that the critically damped curve has the fastest decay.

..

x1 -0.2 0 0.2

0.6

1.0

4 t

Figure 4: Plots of solutions to x + bx + x = 0.

..

MIT OpenCourseWare http://ocw.mit.edu

Fall 2011

18.03SC Differential Equations

For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.

You might also like