Multiple Degree of Freedom Systems: Free Vibrations of Two Degree of Freedom Systems
Figure 8.1: Simple two degree of freedom system
Consider a system comprised of two masses and several springs as shown in Figure 8.1. This is a two degree of freedom system since two coordinates are required to completely specify the configuration of the system. One set of coordinates which can be used is and , measured from the equilibrium position, as shown. However, this is not the only set of coordinates which could be used. In fact, while and seem like a natural and logical choice, in some sense there is a better choice, for reasons to be discussed shortly.
Preliminary Physical Reasoning
Before finding the actual solution to this problem, consider two specific initial displacements from the equilibrium configuration as shown in Figure 8.2 and see what we can deduce about the resulting motion. In Figure 8.2, and are the initial displacements of the two masses.
Figure 8.2: Two specific initial displacements of a simple two degree of freedom system
In case (a), both masses are displaced an equal amount in the same direction. In such a displacement, the central spring is neither stretched nor compressed (relative to its equilibrium length). Due to the symmetry in the system, both masses will undergo the same motion so that the distance between them will remain constant. The central spring remains completely undeformed in the resulting motion and applies no force on the masses. It therefore has no effect on the resulting motion (i.e. for this specific initial displacement, the motion of the system would be the same if the central spring were not present as shown in Figure 8.3.
Figure 8.3: Simple two degree of freedom system with central spring removed
As a result, we expect that if the masses are initially displaced with and released from rest, the resulting motion would have and the oscillations would occur at a frequency of
In the second case (b), the masses are again displaced by equal amounts but in opposite directions. The central spring will extend and compress in this case. However, due to the symmetry in the system, the midpoint of the central spring remains stationary and does not move to the left or the right. As a result, the system behaves as if this were a fixed point, so that in this case we can think of the system as equivalent to that shown in Figure 8.4.
Figure 8.4: Simple two degree of freedom system with central point fixed
The fixed boundary in Figure 8.4 has the same effect on the system as the stationary central point. Since one half of the middle spring appears in each system, the effective spring constant in each system is (remember that, other factors being equal, shorter springs are stiffer). Each mass in Figure 8.4 therefore is supported by two springs in parallel so the effective stiffness of each system is . As a result, we expect that if the masses are initially displaced with and released from rest, the resulting motion would be such that and the oscillations would occur at a frequency of
Equations of Motion
Having selected the coordinates ( and ) to describe the configuration of the system, the next step is to find the equations of motion in terms of the chosen coordinates (and their derivatives). Since this is a two degree of freedom system, we will get two equations of motion. There is usually one equation of motion per degree of freedom in the system. There are many ways to find these equations, and most of these methods are beyond the scope of this course. Here we will focus on the use of Newton’s Laws. As with single degree of freedom systems, when drawing FBD/MAD’s we generally assume that all of the coordinates and their derivatives are positive to get the correct signs in the equations of motion. To make things easier, you may also assume if you wish that one of the coordinates is larger than the other, say assume , although it is not necessary to do so.
Figure 8.5 show the FBD/MAD for each of the masses in the system.
Figure 8.5: FBD/MAD’s for two degree of freedom system
Applying Newtons law’s for mass 1 gives
while for mass 2 we get
so that the two equations of motion for the system are
(8.1a)
(8.1b)
It will be convenient to write these equations in , introducing the concept of and matrices, as
(8.2)
An important feature of these equations is that they are , meaning that both and appear in each of these equations. As a result it is not possible to solve either (8.1a) or (8.1b) independently. They must be solved as a pair.
Solution to Equations of Motion
To solve ve (8.2) we look for solutions in which both masses are undergoing harmonic motion at the same frequency (known as s simple simultaneous harmonic motion). We assume the motion of the masses is of the form
(8.3)
where and are the amplitudes of motion of masses 1 and 2 respectively. With this assumption we see that,
so the equations of motion become
or
which must be satisfied at all times . The only way this can be true is if
(8.4)
There are two possibilities:

This implies that there is no motion (see equation (8.3)). This is called the solution and is not of interest here. The determinant is equal to zero.
This is a result from linear algebra and is required for equation (8.4) to have a nontrivial solution.
For a matrix , the determinant is so
or
(8.5)
This is known as the for this system and can be seen to be a quadratic equation for . %We can solve this equation to get two values for .
The two roots are
or
(8.6)
Equation {8.6) gives us the two natural frequencies for our two degree of freedom system. While the numbering of these frequencies is essentially arbitrary, it is customary to number them from lowest to highest. Therefore is understood to be lowest natural frequency of the system. This is also referred to as the natural frequency. These natural frequencies agree with those expected based on our earlier preliminary physical reasoning. When the system is vibrating at one of these frequencies, we say that it is vibrating in one of its . A two degree of freedom system will therefore have two modes of vibration.
Note that for our two degree of freedom system we obtained a second order polynomial (quadratic) characteristic equation from which we have found two natural frequencies. In general, if we have an degree of freedom system, we will find equations of motion. The solution of these leads to an order polynomial characteristic equation from which we can find natural frequencies.
Having found the natural frequencies, we would like to determine the amplitudes and that occur in each mode. To do so we again consider equation (8.4)
Since we have now found the two natural frequencies, we can use these in turn to find and . However, we cannot simply substitute the values of from (8.6) into ) into equation (8.4) and solve for and due to the fact that the matrix is singular for these two values of (that was the condition we used to find these natural frequencies originally).
Since the square matrix equation (8.4) will be singular, we cannot solve for and uniquely. However, we can find a ratio between and for each . To see this, consider the first natural frequency . Substituting this value into equation (8.4) gives
The first of these equations gives
or
(8.7)
Similarly, the second gives
or, once again,
which shows that it does not matter which of the equations in (8.4) is used.
Now consider the second natural frequency . Substituting this result into equation (8.4) gives
The first of these equations gives
or
(8.8)
The second equation gives
or
as before.
Therefore as we can see, there is a of the amplitudes of motion of the two masses when they are vibrating at one of the system’s natural frequencies. We have also shown that either of the equations in equation (8.4) can be used to determine this amplitude ratio.
In summary, we have shown that for this system
(8.9)
(8.10)
The terms
(8.11)
represent the ratio of the amplitudes of the masses when they execute simple simultaneous harmonic motion at the first and second natural frequencies respectively. These do not determine the values of and uniquely, only their ratio. These are known as the .
To use this information effectively, typically one of the amplitudes in the system is arbitrarily assigned a value of 1. The magnitude of the other amplitude can then be uniquely determined using the appropriate ratio. For example, if we arbitrarily assign , then we see from equation (8.9) that for ,
For the second natural frequency (i.e. in the second mode), equation (8.10) gives
We can represent these results as
(8.12)
With these results, the response of the system in each of the first and second modes as given in \eref{eq:AssumedSolution} can be written as
and
The terms
(8.13)
are known as the first and second respectively. These are also referred to as the modes of the system. Often, these mode shapes are normalized so that they represent a unit vector. We could then equivalently write them as
As we have seen there are several ways that the mode shapes can be expressed. Each of these is correct, although in this course we will most often use the form in equation (8.12). What is important is that that the ratio of the amplitudes in any given mode is as given in equations (8.7) and (8.8).
Response in Terms of Normal Modes
Normal modes can be used to conveniently express the response of the system. We have found two natural frequencies which each satisfy the equations of motion. For each of these frequencies there is an associated mode shape. Since the governing differential equations (8.2) are linear, the total solution can be expressed as a of these two solutions
(8.14)
where and represent the actual amplitude of the motion of mass 1 in the first and second modes respectively, which will depend on the initial conditions as we will discuss shortly. Note that equation (8.14) simply represents two equations
where
It can be verified by direct substitution that these equations are indeed the solutions to the equations of motion.
There are four in the solution: and . The (scalar) quantities and represent of each mode is present in the response while and represent the relative phase angle of each of the modes. Note that the equations in (8.14) can be rewritten as
(8.15)
where and are now the four arbitrary constants. These quantities are determined by the . With four arbitrary constants, we will need four initial conditions to determine them uniquely. The initial conditions used are typically the initial position and velocity of each of the masses. If for example
equation (8.15) at becomes
so that
(8.16)
Differentiating {8.15}
and substituting the above initial conditions leads to
so that
which can again be solved to find
(8.17)
As a result, the total response then becomes
(8.18)
If equation (8.14} was to be used then note that
To see how we can use these results, consider the two specific sets of initial conditions considered previously:
 Both masses displaced the same amount in the same direction and released from rest. The associated initial conditions are
The response from (8.18) is then
which shows that the response of the system is entirely in the first mode. The masses move together at the same frequency () with the same amplitude in the same direction () at all times. In this motion the middle spring remains undeformed.
 Both masses displaced the same amount in opposite directions and released from rest. The associated initial conditions are
The response from (8.18) is then
which shows that now the response of the system is entirely in the second mode. The masses again move with the same amplitude and the same frequency () but in opposite directions () at all times. In this motion the midpoint of the central spring does not move.
This highlights a general result, but a very important one. The response of a multiple degree of freedom linear system will be a linear combination of the normal mode responses of the system, with each mode responding at its associated natural frequency. The amount (amplitude) of each mode and the corresponding phase angle are determined solely by the initial conditions. This is the fundamental idea behind the approach known as . As we have seen above, if a system is started in one of its modes (both position and velocity must satisfy these conditions) then it will vibrate purely in that mode. If the initial conditions do not start the system in one of its modes, then the response will in general be a combination of the modes in the system.
These results completely agree with (and generalize) the previous results we predicted for certain special situations based on purely physical reasoning. Of course, the mathematical approach outlined here will usually be required when the systems become more complex.
Beatin
Beating is a phenomenon that occurs in MDOF systems when natural frequencies are “close” to each other. To see this, consider the simple \tdof system we have been investigating with initial conditions such that only one mass is displaced from the equilibrium position and released from rest
With these initial conditions we find from (8.16) and (8.17) that
so the response will be
or
Using appropriate trigonometric identities, these can equivalently be rewritten
(8.19)
If we define
(8.19) becomes
Now, if the two natural frequencies are close to each other, we can interpret these motions as vibrations occurring at a frequency of with an amplitude that varies slowly at a frequency . We can define the as
and the associated as
Figure 8.6 illustrates this response of the two masses with
so that in this case
Figure 8.6: Beating response of two degree of freedom system
Note that in the beating response, there is a transfer of energy between the two masses. As can be seen in Figure (8.6), when mass 1 is vibrating with a maximum amplitude, the amplitude of mass 2 is small so that mass 1 has a large kinetic energy while mass 2 has almost none. One half of a beating period later the situation is reversed and mass 2 will have a large kinetic energy while mass 1 will have very little. Note further that since the system under consideration is , the total mechanical energy remains constant at all times .