You must log in to edit PetroWiki. Help with editing

Content of PetroWiki is intended for personal use only and to supplement, not replace, engineering judgment. SPE disclaims any and all liability for your use of such content. More information

Message: PetroWiki content is moving to OnePetro! Please note that all projects need to be complete by November 1, 2024, to ensure a smooth transition. Online editing will be turned off on this date.


Vibrating systems with multiple sources of motion

PetroWiki
Jump to navigation Jump to search

From an analytical perspective, systems with two or more sources of motion are modeled with multiple degrees of freedom (DOFs). At the basic level, systems with two and more DOFs are similar to single DOF systems. In a 2DOF system, two independent equations of motion are required to define a system (e.g., a double classical linear oscillator (CLO) or a double pendulum), but the sources need not be the same. For example, a system could have a translational and a rotational DOF. As long as the equations of motion are independent of one other, it is a 2DOF system. These equations must be solved simultaneously. (See a refresher on differential calculus if needed)

Free undamped 2DOF system

As in the SDOF system, an undamped (c1 and c2 = 0) system will be developed first. Fig. 1 shows such a system. As before, Newton’s second law can determine the equations of motion. In the SDOF system, a solution in the form of sine and cosine was used. For the first 2DOF system, another valid form of a solution, a sine with a phase angle, Φ, is used to show another solution form. The equation for mass 1 is

RTENOTITLE....................(1)

and for mass 2 is

RTENOTITLE....................(2)

  • Fig. 1 – 2DOF free and undamped CLO system.

    Fig. 1 – 2DOF free and undamped CLO system.

These two equations of motion must be solved simultaneously because they are coupled through the displacement terms, x1 and x2. As before, a solution is assumed and substituted back into the equations of motion. The assumed solution is

RTENOTITLE....................(3)

and

RTENOTITLE....................(4)

for which the second differentiation is

RTENOTITLE....................(5)

and

RTENOTITLE....................(6)

Substituting back into the equations of motion and collecting terms, the result is:

RTENOTITLE....................(7)

and

RTENOTITLE....................(8)

The only way not to have a trivial solution (C1 and C2 = 0) to these equations is to have the determinate of the coefficients be zero:

RTENOTITLE....................(9)

Using linear algebra, the determinate is:

RTENOTITLE....................(10)

which is a quadratic equation in terms of the square of the natural frequencies, ω2. This solves to:

RTENOTITLE....................(11)

There are two roots to this equation, which means that there are two natural frequencies. The solution therefore must be in terms of the two frequencies. Substituting the first natural frequency back into the equation of motion solution gives the result of the first natural frequency:

RTENOTITLE....................(12)

and

RTENOTITLE....................(13)

The ratio of C11 to C21, known as the mode shape, for the first natural frequency is:

RTENOTITLE....................(14)

therefore,

RTENOTITLE....................(15)

The first natural frequency motions then are:

RTENOTITLE....................(16)

and

RTENOTITLE....................(17)

The mode shape for the second natural frequency is:

RTENOTITLE....................(18)

therefore,

RTENOTITLE....................(19)

The second natural frequency motions then are:

RTENOTITLE....................(20)

and

RTENOTITLE....................(21)

The complete solution then is both displacements added into one equation:

RTENOTITLE....................(22)

and

RTENOTITLE....................(23)

The values of C11, C12, C21, C22, Φ1, and Φ2 depend on the initial conditions and the mode shapes, λ1 and λ2.

Ex. 1 is a 2DOF free and undamped CLO system.

  • Example 1 – 2DOF free and undamped CLO system.

    Example 1 – 2DOF free and undamped CLO system.

Free damped 2DOF system

Adding damping complicates the equations considerably, but the procedure remains the same. In this case, it is easier to use linear algebra. Using the same model as before (see Fig. 2), but adding viscous dampers (c1 and c2 ≠ 0), the equations of motion for the independent DOFs are determined for mass 1 as

RTENOTITLE....................(24)

and for mass 2 as

RTENOTITLE....................(25)

This can be written in matrix form as

RTENOTITLE....................(26)

or as

RTENOTITLE....................(27)

where M = the mass matrix, C = the damping matrix, K = the stiffness matrix, Ẍ = the acceleration vector, Ẋ = the velocity vector, and X = the displacement vector, which are given as:

RTENOTITLE....................(28)

RTENOTITLE....................(29)

RTENOTITLE....................(30)

RTENOTITLE....................(31)

RTENOTITLE....................(32)

and

RTENOTITLE....................(33)

If the solution is assumed to be of the form

RTENOTITLE....................(34)

then

RTENOTITLE....................(35)

and

RTENOTITLE....................(36)

Substituting back into the equation of motion, the result in matrix form is:

RTENOTITLE....................(37a)

which can be rewritten as

RTENOTITLE....................(37b)

  • Fig. 2 – 2DOF free and damped CLO system.

    Fig. 2 – 2DOF free and damped CLO system.

Because time is always positive and a nontrivial solution is desired, the only way this equation is true is if the determinate of the coefficient of Ceωt is zero; that is, if

RTENOTITLE....................(38)

The determinate of Eq. 38 is a fourth-order polynomial in terms of a, which means that there are four roots. These roots can be:

  1. Four real and negative roots.
  2. Two sets of complex conjugates with negative real parts.
  3. Two real and negative roots and one set of complex conjugates.

If number one is the case, then the result is an exponentially decaying motion without oscillation. It is similar to the overdamped case for an SDOF system. If number two is the case, then the motions will be exponentially decaying oscillations for both DOFs. This is similar to the underdamped case for an SDOF system. Finally, for case number three, either condition can occur.

Ex. 2 is a 2DOF free and underdamped CLO system.

  • Example 2 – 2DOF free and damped CLO system.

    Example 2 – 2DOF free and damped CLO system.


More information on linear algebra can be found on the mathematics of fluid flow page.

Forced damped 2DOF system

Adding forcing complicates the equations considerably yet again. The procedure is the same, however, and it is a matter of keeping the mathematics straight. Many texts are available to delve more deeply into this subject. For more information, please refer to the literature.[1][2][3][4][5][6][7][8][9]

Multiple-DOF systems

The previous discussion of 2DOF systems points out how to handle any DOF system. The last example used matrix notation to define the system for the solution process. Multiple-DOF systems are solved similarly with the primary difference being the degree of the defining matrices is greater, as is the degree of difficulty in solving the system. The matrix will have the same number of rows and columns as the degree of freedom. There are other methods (e.g., finite-element modeling) that can be used to tackle the complexity of multiple-DOF systems.

Nomenclature

c = axial damping coefficient, mL/t, lbf-ft/sec
C = constant of integration, various
C = damping matrix, mL/t, lbf-ft/sec
k = spring constant, m/t2, lbf/in.
K = stiffness matrix, m/t2, lbf/in.
m = mass, m, lbm
M = mass matrix, m, lbm
t = time, seconds
x = displacement, L, in.
ẋ = first derivative with respect to time of displacement (velocity), L/t, ft/sec
ẍ = second derivative with respect to time of displacement (acceleration), L/t 2 , ft/sec 2
X = displacement vector
Ẋ = velocity vector, L/t, ft/sec
Ẍ = acceleration vector, L/t2, ft/sec2
λ = mode shape, dimensionless
Φ = phase angle, rad
ω = frequency, 1/t, Hz

References

  1. ↑ Achenbach, J.D. 1973. Wave Propagation in Elastic Solids, seventh edition. Amsterdam: North-Holland Publishing Co.
  2. ↑ Clough, R.W. and Penzien, J. 1975. Dynamics of Structures. New York: McGraw-Hill Book Co.
  3. ↑ Den Hartog, J.P. 1934. Mechanical Vibrations. New York: Dover Publications.
  4. ↑ Meirovitch, L. 1986. Elements of Vibration Analysis, second edition. New York: McGraw-Hill Book Co.
  5. ↑ Nashif, A.D., Jones, D.I.G., and Henderson, J.P. 1985. Vibration Damping. New York: John Wiley & Sons.
  6. ↑ Seto, W. 1964. Mechanical Vibrations. New York: McGraw-Hill Book Co.
  7. ↑ Shabana, A.A. 1991. Theory of Vibration, Volume I: An Introduction. New York: Springer-Verlag.
  8. ↑ Shabana, A.A. 1991. Theory of Vibration, Volume II: Discrete and Continuous Systems. New York: Springer-Verlag.
  9. ↑ Elmore, W.C. and Heald, M.A. 1969. Physics of Waves. New York: Dover Publications.

Noteworthy papers in OnePetro

Use this section to list papers in OnePetro that a reader who wants to learn more should definitely read

External links

Use this section to provide links to relevant material on websites other than PetroWiki and OnePetro

See also

Vibration theory

Basic vibration analysis

Continuous vibration systems

Differential calculus refresher

PEH:Mathematics_of_Vibrating_Systems

Retrieved from "https://petrowiki.spe.org/w/index.php?title=Vibrating_systems_with_multiple_sources_of_motion&oldid=47957"