2 Modelling free mechanical oscillations

Consider the following schematic diagram of a shock absorber:

Figure 6

No alt text was set. Please request alt text from the person who provided you with this resource.

The equation of motion can be described in terms of the vertical displacement x of the mass.

Let m be the mass, k d x d t the damping force resulting from the dashpot and n x the restoring force resulting from the spring. Here, k and n are constants.

Then the equation of motion is

m d 2 x d t 2 = k d x d t n x .

Suppose that the mass is displaced a distance x 0 initially and released from rest. Then at t = 0 , x = x 0 and d x d t = 0 . Writing the differential equation in standard form gives

m d 2 x d t 2 + k d x d t + n x = 0.

We shall see that the nature of the oscillations described by this differential equation depends crucially upon the relative values of the mechanical constants m , k and n . This will be explored in subsequent Tasks.

Task!

Find and solve the auxiliary equation of the differential equation

m d 2 x d t 2 + k d x d t + n x = 0.

Putting x = e λ t , the auxiliary equation is m λ 2 + k λ + n = 0.

Hence λ = k ± k 2 4 m n 2 m .

The value of k controls the amount of damping in the system. We explore the solution for various values of k .

2.1 Case 1: No damping

If k = 0 then there is no damping. We expect, in this case, that once motion has started it will continue for ever. The motion that ensues is called simple harmonic motion . In this case we have

λ = ± 4 m n 2 m , that is, λ = ± n m i where i 2 = 1.

and the solution for the displacement x is:

x = A cos n m t + B sin n m t where A , B are arbitrary constants.

Task!

Impose the initial conditions x = x 0 and d x d t = 0 at t = 0 to find the unique solution to the ODE:

d x d t = n m A sin n m t + n m B cos n m t

When t = 0 , d x d t = 0 so that n m B = 0 so that B = 0.

Therefore x = A cos n m t .

Imposing the remaining initial condition: when t = 0 , x = x 0 so that x 0 = A and finally:

x = x 0 cos n m t .

2.2 Case 2: Light damping

If k 2 4 m n < 0 , i.e. k 2 < 4 m n then the roots of the auxiliary equation are complex:

λ 1 = k + i 4 m n k 2 2 m λ 2 = k i 4 m n k 2 2 m

Then, after some rearrangement:

x = e k t 2 m A cos p t + B sin p t in which p = 4 m n k 2 2 m .

Task!

If m = 1 , n = 1 and k = 1 find λ 1 and λ 2 and then find the solution for the displacement x .

λ = 1 + i 4 1 2 = 1 2 ± i 3 2 . Hence x = e t 2 A cos 3 2 t + B sin 3 2 t .

Impose the initial conditions x = x 0 , d x d t = 0 at t = 0 to find the arbitrary constants and hence find the solution to the ODE:

Differentiating, we obtain

d x d t = 1 2 e t 2 A cos 3 2 t + B sin 3 2 t + e t 2 3 2 A sin 3 2 t + 3 2 B cos 3 2 t

At t = 0 ,

x = x 0 = A (i)

d x d t = 0 = 1 2 A + 3 2 B (ii)

Solving (i) and (ii) we obtain

A = x 0 B = 3 3 x 0 then x = x 0 e t 2 cos 3 2 t + 3 3 sin 3 2 t . The graph of x against t is shown in Figure 7. This is the case of light damping. As the damping in the system decreases (i.e. k 0 ) the number of oscillations (in a given time interval) will increase. In many mechanical systems these oscillations are usually unwanted and the designer would choose a value of k to either reduce them or to eliminate them altogether. For the choice k 2 = 4 m n , known as the critical damping case, all the oscillations are absent.

Figure 7

No alt text was set. Please request alt text from the person who provided you with this resource.

2.3 Case 3: Heavy damping

If k 2 4 m n > 0 , i.e. k 2 > 4 m n , then there are two real roots of the auxiliary equation, λ 1 and λ 2 :

λ 1 = k + k 2 4 m n 2 m λ 2 = k k 2 4 m n 2 m

Then

x = A e λ 1 t + B e λ 2 t .

Task!

If m = 1 , n = 1 and k = 2.5 find λ 1 and λ 2 and then find the solution for the displacement x .

λ = 2.5 ± 6.25 4 2 = 1.25 ± 0.75

Hence λ 1 , λ 2 = 0.5 , 2 and so x = A e 0.5 t + B e 2 t Impose the initial conditions x = x 0 , d x d t = 0 at t = 0 to find the arbitrary constants and hence find the solution to the ODE.

Differentiating, we obtain

d x d t = 0.5 A e 0.5 t 2 B e 2 t

At t = 0 ,

x = x 0 = A + B (i)

d x d t = 0 = 0.5 A 2 B (ii)

Solving (i) and (ii) we obtain A = 4 3 x 0 B = 1 3 x 0 then x = 1 3 x 0 ( 4 e 0.5 t e 2 t ) .

The graph of x against t is shown below. This is the case of heavy damping.

No alt text was set. Please request alt text from the person who provided you with this resource.

Other cases are dealt with in the Exercises at the end of the Section.