PEH:Mathematics of Vibrating Systems

Revision as of 16:32, 26 April 2017 by Denise Watts (Denisewatts) (talk | contribs)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)

Publication Information

Vol1GECover.png

Petroleum Engineering Handbook

Larry W. Lake, Editor-in-Chief

Volume I – General Engineering

John R. Fanchi, Editor

Chapter 1 – Mathematics of Vibrating Systems

Alfred W. Eustes, III, PhD., PE, SPE, Colorado School of Mines

Pgs. 1-43

ISBN 978-1-55563-108-6
Get permission for reuse



Introduction


Many mathematical tools can be used to analyze vibrational systems. One of the first mathematical tools a neophyte engineer learns is calculus. The basics of limits, differentiation, and integration permeate all of engineering mathematics. This chapter offers a cursory review of these topics and uses the mathematics of vibrations to demonstrate how the concepts operate. For more specific information on all these topics, consult relevant sections of this Handbook.

Many of the mathematical tools engineers use to evaluate and predict behavior, such as vibrations, require equations that have continuously varying terms. Often, there are many terms regarding the rate of change, or the rate of change of the rate of change, and so forth, with respect to some basis. For example, a velocity is the rate of change of distance with respect to time. Acceleration is the rate of change of the velocity, which makes it the rate of change of the rate of change of distance with respect to time. Determining the solutions to these types of equations is the basis of differential calculus.

An equation with continuously varying terms is a differential equation. If only one basis is changing, then it is an ordinary differential equation (ODE); however, if two or more bases are changing, then it is a partial differential equation (PDE). An ODE uses the notation "d" and a PDE uses ∂ to refer to change.

Understanding differentiation starts with an understanding of limits.

Limits

A graph is a useful method for determining how an equation behaves. The independent variable t in Eq. 1.1 determines how the dependent variable y behaves. The operators and constants in an equation specify this behavior. Fig. 1.1 shows the graph of Eq. 1.1, the distance of freefall over time with an initial velocity of zero. Down is considered negative in this equation:

RTENOTITLE....................(1.1)

The x-axis (abscissa) usually is the independent variable, and the y-axis (ordinate) usually is the dependent variable; however, many drilling charts hold an exception to this generality, in that their ordinate often is the independent variable, and their abscissa is the dependent one. An example of such a drilling chart is the depth vs. time graph.


In Fig. 1.1, at the time of 3 seconds, the distance is –96.522 ft. A tangent line to the graph at 3 seconds is known as the slope (A) of the graph at that point. To quickly estimate the slope of the tangent, divide the rise (Δy) by the run (Δt), as shown in Eq. 1.2:

RTENOTITLE....................(1.2)

In this case, the tangent y value at 2 seconds is –48.261 ft and at 4 seconds is –241.305 ft. The slope then is:

RTENOTITLE....................(1.3)

Because the units in this case are ft/sec, this slope gives the velocity at that point. It is the rate of change of the distance with respect to time.

A limit is defined as the value of a function at a given point as that point is approached from either higher or lower values (often referred to as approaching from the left or right, respectively). The limit (Y) of Eq. 1.1 at 3 seconds is:

RTENOTITLE....................(1.4)

Y is known as the limit of the function. In this simple case, Y is the same regardless of whether t approaches 3 from the left or the right. This is not true in all cases, however (e.g., with a discontinuous function). In these cases, the limit can be determined analytically. One can also determine the limit using a graph such as in Fig. 1.1.

Limits have the following properties:

RTENOTITLE....................(1.5)

RTENOTITLE....................(1.6)

RTENOTITLE....................(1.7)

and

RTENOTITLE....................(1.8)

Derivatives

As noted earlier, the slope of graph of Eq. 1.1 at 3 seconds = –96.522 ft/sec and is the velocity (v) of free-fall at 3 seconds from release. This value is known as the first derivative of Eq. 1.1 at the value of 3. It is written as:

RTENOTITLE....................(1.9)

and is defined as:

RTENOTITLE....................(1.10)

As the limit of the value of Δt approaches zero, the solution converges to the first derivative.

Derivatives have the following properties (r = constant).

RTENOTITLE....................(1.11)

RTENOTITLE....................(1.12)

RTENOTITLE....................(1.13)

RTENOTITLE....................(1.14)

RTENOTITLE....................(1.15)

RTENOTITLE....................(1.16)

In the case of Eq. 1.7, where Q = 0, L’Hopital’s rule can help find the limit. This is shown in Eq. 1.17:

RTENOTITLE....................(1.17)

Other rules regarding differentials are the following.

The linear superposition rule:

RTENOTITLE....................(1.18)

The product rule:

RTENOTITLE....................(1.19)

The quotient rule:

RTENOTITLE....................(1.20)

The chain rule (or function of a function):

RTENOTITLE....................(1.21)

Multiple differentiations can be shown by

RTENOTITLE....................(1.22)

and continued differentiations can be shown by

RTENOTITLE....................(1.23)

A useful point to recognize is where a slope equals zero, which can correspond to a maximum, a minimum, or an inflection. To determine these points, determine a first derivative of an equation. Then, set this first-derivative equation to equal zero and solve for the basis (the unknown). To determine whether this point is a maximum, a minimum, or an inflection, determine the second derivative of that equation. If that value is negative, the point is a maximum; if it is positive, the point is a minimum; and if it is zero, the point is an inflection.

The graph of Eq. 1.24 ( Fig. 1.2 ) is an example of this process:

RTENOTITLE....................(1.24)

The first derivative of Eq. 1.24 is:

RTENOTITLE....................(1.25)

which, when set equal to zero, is a quadratic equation with two roots, t = 3 and 1/3 . These two points correspond to the maximum and minimum points on the graph. To prove which is which, a second derivative is taken:

RTENOTITLE....................(1.26)

which at t = 3 and 1/3 is equal to 8 and –8, respectively. This means that at t = 3, the function is at a minimum and at t = 1/3, the function is at a maximum.


The first differentiation of the equation of the position of a free-falling object starting at rest (Eq. 1.1) gives the slope of the graph, which, as noted, is the velocity:

RTENOTITLE....................(1.27)

A second differentiation gives the change of the slope with respect to time (acceleration), and is:

RTENOTITLE....................(1.28)

which is the acceleration caused by Earth’s gravity.

Differential-Equation Solutions

Solutions to differential equations solved in closed form can range from trivial to impossible. Numerical methods often are required. Nevertheless, some general strategies have been developed to solve differential equations.[1][2][3]

An ODE with only first derivatives is known as a first-order ODE. A second-order ODE has second and possibly first derivatives. The same reasoning applies to third order and beyond. Likewise, when a PDE has only first derivatives, it is a first-order PDE. The second and third orders and beyond are defined on the basis of their highest-order derivative.

This section has covered some of the basics of ODE and PDE mathematics. The reader is urged to review mathematical texts and handbooks for more details on this subject.

To illustrate some of the aforementioned principles of ODE and PDE mathematics, the rest of this chapter discusses examples within the subject of vibrations. Vibration mathematics has been chosen because vibrations are notorious for instigating drilling and production problems. The next section covers some basic information on the nature of vibrations, with subsequent sections devoted to some aspect of the mathematical models of vibrations.

Introductory Thoughts About Vibrations


The fundamental theories of vibration are not new. Indeed, Saint-Venant[4] published his theory on the vibrations of rods in 1867, and Love[5] published an entire treatise on vibration theory in 1926. The mathematics of vibration theory involves infinite series, complex functions, and Fourier integral transforms, and its physics involves Newtonian mechanics and stress analyses. Until recently, except under relatively simple conditions, the complexity of such mathematics had restrained the application of vibration theory to solving simple common problems. Now, however, state-of-the-art computers can perform these complex calculations in a reasonable time frame, making possible a wave of new studies.

Vibration Theory


A vibration is a fluctuating motion about an equilibrium state. There are two types of vibration: deterministic and random. A deterministic vibration is one that can be characterized precisely, whereas a random vibration only can be analyzed statistically. The vibration generated by a pumping unit is an example of a deterministic vibration, and an intermittent sticking problem within the same system is a random vibration.

In mechanical systems, deterministic vibrations are excitations that elicit a response from a system, as shown schematically in Fig. 1.3. In theory, as long as two of the three variables (excitation, system, and response) are known, the third one can be determined; however, the mathematics might be challenging. Most often, the response function is sought, so that the excitation function and the system must be known.


Vibration systems can be linear or nonlinear, and discrete or continuous (Fig. 1.4). In all cases, a vibration system can be in one, two, or three mutually orthogonal dimensions. A linear system is a system in which proportionality (Eq. 1.29) and superposition (Eq. 1.30) are true, that is, in which:

RTENOTITLE....................(1.29)

and

RTENOTITLE....................(1.30)

When proportionality and superposition are not true, then the system is nonlinear.


A discrete system is one having a finite number of independent coordinates that can describe a system response. These independent coordinates are known as degrees of freedom (DOFs). If the motion of mass, either translational or rotational, of a vibrating system is a function of only one independent coordinate, then the system has one DOF. If two or more independent coordinates are required to describe one or both types of motion, then the system has two or more DOFs. If a system is continuous (an infinite set of independent coordinates is needed to describe the system response), it has an infinite number of DOFs. Because material structures all have a continuous nature, all systems have an infinite number of DOFs. Most systems have dominant DOFs; some even have a single dominant DOF. Such systems therefore can be characterized as discrete systems, which makes the mathematics more tractable.

If a system has a single DOF or set of DOFs in only one direction, it is a 1D system. If there are two mutually orthogonal directions for the DOF, it is a 2D system; and if there are three mutually orthogonal directions for the DOF, it is a 3D system.

As Fig. 1.5 shows, the excitation function can be periodic or transient, and absent or present. A periodic vibration is one that can be characterized mathematically as an indefinite repetition. A transient vibration is of finite length and is composed of waves that have a definite beginning and that eventually die out. These waves can be of extremely short duration or last for some time.


A standing wave is a vibration whose wave profile appears to be standing still, though actually the particles that make up the material are oscillating about an equilibrium position. Because of the geometry and boundary conditions of the material through which they are traveling, the waves and the reflected waves cancel and reinforce themselves over the same location in the material, which makes the wave profile appear not to be moving. The point at which no motion is occurring is a nodal point, or node. The point of maximum amplitude is the antinode.

In reality, all waves are transient in some way. If a wave is repeated over a longer time than it takes for a single wave to propagate through a material, then this series of waves can be called a vibration. All vibrations are transient, as well. If the vibration lasts longer than the time under analysis, then it can be characterized as infinite in length.

When the excitation is present and is actively affecting the system within the analysis time frame, the response is called a forced vibration. The response of a system with an absent excitation function—one that is not present within the analysis time frame—is called a free vibration. As such, the system can be responding to the removal of an excitation function. For example, if the response of a mass and spring system is sought after the system has been pulled down and released, the original excitation function (the pulling force) is considered absent because the analysis is being performed after the release.

Wave Propagation


The method by which a vibration travels through a system is known as wave propagation. When an external force is impressed on a real-world elastic body, the body does not react instantly over its entire length. The point immediately under the external force reacts first, and then the section just under that point reacts to the previous section’s reaction, and so on. This series of reactions is called wave propagation because the reactions propagate through the body over a period of time at a specific velocity. If the rate of change of the external force is slow enough, static equilibrium analysis can model the reactions adequately for most engineering applications. This is called rigid-body analysis. If the external force changes rapidly, however, wave-propagation analysis is necessary to model the reactions effectively.

There are many types of elastic waves. Some listed in this section are longitudinal, lateral, and bending waves.[6] Some of these are shown in Fig. 1.6.


In longitudinal waves (also variously called compression/tension, axial, dilatational, and irrotational waves), the particles that make up the elastic medium are forced directly toward and away from each other, and the direction of the particles’ motion is parallel to that of the wave motion. In most steels, longitudinal waves travel at ≈16,800 ft/sec. Longitudinal waves are not dispersive. This means that all the wave components that make up a longitudinal wave travel at the same velocity and, hence, do not separate (disperse).

In lateral waves (also known variously as shear, torsional, transverse, equivoluminal, and distortional waves), the particles slip beside each other, and move perpendicular to the direction of the wave motion. Because slipping uses more energy, lateral waves are slower. In steel, for example, they travel at 10,400 ft/sec. A rapidly changing torsional force on a section of pipe will cause a lateral wave to propagate from the point of application to all other parts of the pipe. It propagates as an angular twist. Lateral waves are nondispersive and have a similar solution method as the longitudinal waves; however, shear or transverse waves are dispersive (i.e., the wave components that make up a shear wave travel at different velocities). Their wave components will disperse and "smear" the initial wave profile. This complicates the analysis significantly.

Bending waves (flexural waves) travel as a bend in a bar or plate and have longitudinal and lateral components. Rotary shears and moments of inertia complicate bending-wave analysis.

Wave-propagation studies in petroleum engineering areas generally have been confined to longitudinal, torsional, and lateral waves. Of these, longitudinal waves generally are easiest to model and are considered in this section. A compression wave is a stress wave in which the propagated stress is in compression. Likewise, a tension wave is a stress wave in which the propagated stress is in tension.

Wave Behavior


Wave velocity depends primarily on density and modulus of elasticity but also is affected by damping and frequency. For example, hitting one end of a long steel rod with a hammer will generate a longitudinal wave that compresses the particles of the steel. The wave’s length is set by the length of time that the hammer is in contact with the end of the rod, whereas its magnitude is set by the force of the hammer blow. As the wave moves along the rod, the steel within the length of the wave is compressed. After the wave passes, the steel returns to its unstressed state, though not necessarily in the same location as before the wave passage.

As another example, twisting (shearing) a steel rod will generate a shear wave. A shear wave moves along the rod more slowly than the longitudinal wave does. Similarly to the longitudinal wave discussed above, its length is set by the duration of the twisting action, whereas its magnitude is set by the torque from the twisting action.

Waves act independently, but the stresses they create can be additive. For example, two equal compression waves that are generated simultaneously by hammer blows at each end of a long steel rod will meet in the center of the rod, pass through one another, and then each continue along the rod as if the other never existed (independence). While the waves are passing each other, however, the compression in the steel will be twice (additive) that of either wave.

Natural Frequencies and Resonance

Everything has a natural frequency, a frequency at which it would vibrate were it given the energy to vibrate and left alone. For instance, the human body has a natural frequency of ≈5 cycles/sec. All drill and rod strings have a natural frequency that depends on the material properties and geometry. The material properties determine the wave velocity, and the geometry determines how waves are reflected and refracted.

During wave propagation, the wave eventually reaches an end of the material. Some of the wave will reflect back to its source. If the reflection reaches the source at the same time a new wave is generated, the two waves will combine and be synchronized in phase. Later, if those two waves’ reflections return to the source at the same time the next new wave is generated, all three waves will combine. This will continue for as long as waves are generated under these conditions, and the resultant wave will increase in amplitude, theoretically to infinity. This is called resonance. The frequency at which resonance occurs is the natural frequency or an integer multiple of that frequency (called a harmonic). If this wave reinforcement is allowed to continue, the system eventually will either self-destruct or fatigue to failure.

A continuous system contains an infinite number of natural frequencies, whereas a discrete single-degree-of-freedom (SDOF) system (e.g., a point mass on a massless spring) has only one natural frequency. If two point masses are connected using two springs, then there are two natural frequencies in this 2DOF system. In general, the number of DOFs in a system determines the number of natural frequencies it has, which means that any discrete system will have a finite number of natural frequencies; however, in reality, there is an infinite number of natural frequencies because all systems are continuous. Some frequencies will have higher amplitudes than others. Such continuous systems with discrete higher-amplitude responses can be modeled with a discrete methodology.

Damping

Resonance energy does not reach an infinite value because of damping, the dissipation of energy over time or distance.[7] Without damping, or friction, the energy from vibrations would build until there is more energy than the structure can sustain, which can cause structural failure.

A wave propagating into a system adds energy to a system, whereas damping removes it. Generally, the dissipated energy from the vibration is converted to heat, and if damping does not take enough energy out of a system, the system can self-destruct from energy overload. The amount of energy in a system at a given time is reflected in the system’s stress/strain level. The more stresses/strains in the system, the higher the energy level. Once the stresses reach a value greater than the yield strength of the system, yield failure is imminent. If the stresses are greater than the ultimate strength of the material, failure is immediate.

In the borehole, three distinctive types of damping occur: viscous, Coulomb, and hysteretic. Viscous damping occurs when the damping force generated is proportional to the velocity of the particles. Coulomb damping (also called dry friction) is the force generated by the movement of materials past one another, and it usually is proportional to the force normal to the materials’ surfaces. The dynamic and static coefficients of friction are the proportionality constants. Hysteretic damping is the friction force generated by the relative motion of the internal planes of a material as a wave causes particle motion. Although this is true of all materials, some materials are viscoelastic (i.e., they show a much larger hysteretic effect than do others).

Viscous Damping. As noted above, viscous damping occurs when the damping force is proportional to the velocity of the particles. Viscous damping is shown by:

RTENOTITLE....................(1.31)

One way that viscous damping arises in jarring analysis is from the interaction of a solid and liquid at their interface, such as where the steel contacts the liquid mud along the sides of a drillstring.

One method for determining the damping involves noting the decrement of acceleration over one vibration cycle.[8] An impulse is impressed on the drillstring to produce a wave. While the wave is decaying, the acceleration is measured and recorded multiple times at one location on the string and at the same phase (i.e., crest to crest). The time between recordings also is noted. These values are used in Eq. 1.32 to compute the damping coefficient (c). Unfortunately, though, this method gives the total damping and does not distinguish between viscous and Coulomb damping.

RTENOTITLE....................(1.32)

Coulomb damping is the friction that occurs when two dry surfaces slide over each other, and its force is a constant value that is independent of particle velocity and displacement, but dependent on the friction factor (μ) and the force normal to the friction surface. This value is:

RTENOTITLE....................(1.33)

The Coulomb damping force always is of the opposite sign from that of the particle velocity, so that the damping force reverses when the particle velocity changes signs. This discontinuity makes it a nonlinear damping force, shown as:

RTENOTITLE....................(1.34)

Nonlinearity makes a closed-form solution to an equation of motion difficult.

Hysteretic Damping. Hysteretic damping also is called structural damping because it arises from internal friction within a structure. A wave moves through a material because the atomic structure is reacting to an applied force. As the atoms of the structure move, energy is lost through the interaction of these atoms with their neighboring atoms. Hysteretic damping is the energy lost when atoms move relative to each other.

If a material had a perfectly linear stress/strain relationship, hysteretic damping would not occur. In reality, though, there is no such thing as a perfectly linear stress/strain curve. Two curves develop on the stress/strain diagram while a material is stressed and relieved. The center area between these two curves represents the energy lost to internal friction. (This hysteresis loop is the reason for the name of this damping type.) This variation can be small, but the amount of energy dissipated can be large because high-frequency vibrations can cause this loop to be repeated many, many times over a given time period.[9]

The hysteretic-damping value is highly dependent on a number of factors. One factor is the condition of the material (i.e., chemical composition, inhomogeneities, and property changes caused by thermal and stress histories). Another is the state of internal stress from initial and subsequent thermal and stress histories. Also, the type and variation of stress—axial, torsional, shear, and/or bending—affect the hysteretic-damping value.

A way of looking at hysteretic-damping force is to set it proportional to the particle velocity divided by the wave frequency. This is shown in Eq. 1.35.[10]

RTENOTITLE....................(1.35)

Equivalent Springs


Many systems can be modeled as multiple springs. Such springs can be combined into a single, equivalent spring (Fig. 1.7. For parallel springs, the sum of the spring constants is equal to the equivalent spring constant (Eq. 1.36). For series springs, the reciprocal of the sum of the reciprocals of the spring constants is equal to the equivalent spring constant (Eq. 1.37). A linear spring oscillates in a single translational direction. A torsional spring oscillates with an angular twist (Eq. 1.38).

Boundary and Initial Conditions


The boundary conditions (how the ends of a system are attached) and initial condition (condition of the system at the start in time) are extremely important in vibration and wave propagation analysis. The specific solution of any ODE or PDE requires a set of boundary and/or initial conditions. Usually, a displacement (boundary condition) and an initial velocity (initial condition) are specified.

In wave propagation, the boundary conditions also dictate wave behavior. For example, a compression wave is reflected from a free end as a tension wave and from a fixed end as a compression wave. If two rods are connected at their ends and are of different geometry or material, then a fraction of the energy of the wave is reflected and the remaining portion of the energy is refracted at their connection. Other types of boundaries direct the system response by limiting the DOF. This includes boundary conditions of pinned, revolute, translational, translational and rotational; forcing function; mass spring and/or damper, and a semi-infinite connection. In addition, changes in material properties will affect the various constants and will cause wave-propagation reflections and refractions at the boundary between the properties. Fig. 1.8 shows some typical boundary conditions.

Mechanical Vibration Analysis


There are three components to mechanical vibration analysis: to determine the geometric compatibilities; to determine the constitutive (material properties) equations; and to determine the equilibrium condition.

The geometric compatibilities are the displacement constraints and connections. They also include the continuous properties, which state that the system does not separate into individual pieces. (If it does, that is another problem altogether.) The constitutive equations represent the material properties, which include mass, damping, and spring coefficients. These constitutive equations include stress/strain relationships and Hooke’s law (Eqs. 1.39a and 1.39b):

RTENOTITLE....................(1.39a)

or, in another form,

RTENOTITLE....................(1.39b)

The coefficient of Δl in Eq. 1.39b often is called the spring constant or stiffness constant.

The equilibrium condition is based on both static and dynamic conditions. A static equilibrium states that the sum of the forces acting on an object is equal to zero:

RTENOTITLE....................(1.40)

A dynamic equilibrium is based on Newton’s second law and is the basis of many vibration analysis methods. The sum of the forces acting on an object is equal to its mass times the acceleration of the object. Other dynamic-equilibrium analysis includes virtual work methods and energy-balance methods (Hamilton’s principle).

Newton’s second law for a translational system is

RTENOTITLE....................(1.41)

and for torsional systems is

RTENOTITLE....................(1.42)

Newton’s second law can be rewritten in a form known as D’Alembert’s principle:

RTENOTITLE....................(1.43)

in which mẌ is treated as a force and is called an inertial force.

Some basic equations of vibration analysis are shown in Table 1.1.

SDOF System


An SDOF system is the most basic vibration analysis. The typical SDOF system is the classical linear oscillator (CLO), as shown in Fig. 1.9. It consists of a point mass, spring, and damper.

Free Vibration Without Damping

The first analysis is free vibration without damping. Using Newton’s second law and D’Alembert’s principle, the equation that describes free vibration without damping (c = 0) is:

RTENOTITLE....................(1.52)

which gives RTENOTITLE, which when rearranged becomes:

RTENOTITLE....................(1.53)

The solution to this differential equation is:

RTENOTITLE....................(1.54)

The constant multiplying the t is the natural frequency of the system and is:

RTENOTITLE....................(1.55)

in radians/unit time. Multiplying by 1/2π gives the natural frequency in cycles/unit time. When Eq. 1.55 is substituted into Eq. 1.54, the result is:

RTENOTITLE....................(1.56)

The constants C1 and C2 are based on the initial and boundary conditions. If at time 0, x = x0 and =v0, the initial location and velocity, respectively, the first coefficient is:

RTENOTITLE....................(1.57)

which is the initial location. Differentiating once gives:

RTENOTITLE....................(1.58)

which, when t = 0, gives the other coefficient:

RTENOTITLE....................(1.59)

which is based on the initial velocity. The entire equation then is:

RTENOTITLE....................(1.60)

Ex. 1.1 is an SDOF free and undamped CLO system.

Free Vibration With Damping

The second analysis of free vibration is with damping (Fig. 1.10). Using Newton’s second law, the equation that describes free vibration with damping (c ≠ 0) is:

RTENOTITLE....................(1.61)

which is rearranged as before to get:

RTENOTITLE....................(1.62)

The general solution to this differential equation is:

RTENOTITLE....................(1.63)

although the specific solution depends on the value under the square root. When c2 - 4mk = 0, the system is critically damped. Another way to look at this critical damping point is:

RTENOTITLE....................(1.64)

Often, the damping coefficient is divided by the critical damping coefficient to get the critical damping ratio:

RTENOTITLE....................(1.65)

If ξ > 1, the system is underdamped. When disturbed, the system will experience an oscillating decay. If ξ < 1, the system is overdamped and, when disturbed, will die out without oscillating. If ξ = 1, the system is critically damped and also will not oscillate.


A useful simplifying equation is the "damped" natural frequency, ωd. It is:

RTENOTITLE....................(1.66)

If the system is underdamped, that is, if 0 ≤ ξ < 1 , the solution is:

RTENOTITLE....................(1.67)

or

RTENOTITLE....................(1.68)

where RTENOTITLE. If the initial and boundary coefficients are the same as before, then the solution is:

RTENOTITLE....................(1.69)

or

RTENOTITLE....................(1.70)

Ex. 1.2 is an SDOF free and underdamped CLO system.


If the system is overdamped, that is, if ξ>1 , the solution is:

RTENOTITLE....................(1.71)

but in this case, RTENOTITLE. The order changed because it was an imaginary number. With the same initial and boundary conditions as before, the solution is as before:

RTENOTITLE....................(1.72)

Ex. 1.3 is an SDOF free and overdamped CLO system.


If the system is critically damped, that is, if ξ=1 , the solution with the initial and boundary conditions is:

RTENOTITLE....................(1.73)

Ex. 1.4 is an SDOF free and critically damped CLO system.

Forced Vibration Without Damping

The next sets of systems have a forcing function driving the vibration. The first of these is a CLO without damping (c = 0), as shown in Fig. 1.11. The equation of motion for this system with F = F0 sin ωft, a sinusoidally varying force, using Newton’s second law, is:

RTENOTITLE....................(1.74)

In this case, there are two terms in the solution, the homogenous or transient term, and the particular or steady-state term. The homogenous term is the same as in a free-vibration case and is solved by setting the forcing function to zero (that is, the free-vibration case, Eq. 1.60). If the same initial and boundary conditions are applied as before, the solution for the homogenous case is the same as before:

RTENOTITLE....................(1.75)

The second term is the effect of the forcing function on the system. This is solved by assuming a particular solution and deriving it back:

RTENOTITLE....................(1.76)

RTENOTITLE....................(1.77)

RTENOTITLE....................(1.78)

Substituting the above equations into Eq. 1.74 gives:

RTENOTITLE RTENOTITLE....................(1.79)

Collecting the terms gives:

RTENOTITLE....................(1.80)

Equating coefficients shows that:

RTENOTITLE....................(1.81)

and

RTENOTITLE....................(1.82)

Therefore, the particular solution is:

RTENOTITLE....................(1.83)

which can be rewritten as:

RTENOTITLE....................(1.84)

where the reciprocal term in the parentheses sometimes is called the magnification factor. The total solution is:

RTENOTITLE
RTENOTITLE....................(1.85)

Ex. 1.5 is an SDOF forced and undamped CLO system.


Note that when the forcing frequency, ωf, matches the natural frequency, ωn, the value of the coefficient is infinity. This is the resonance condition, and it can lead to excessively large displacements (see Fig. 1.12).


When the forcing frequency is close to but not at the natural frequency, a beating phenomenon occurs. This appears as a low frequency impressed over the frequency of the system. When the engines of a twin-engine aircraft are not quite synchronized, for example, one can hear a beating sound as a low-frequency pulse (the "wow-wow" throb). Fig. 1.13 illustrates this beating phenomenon.

Forced Vibration With Damping

The second system with a forcing function driving the vibration is a CLO with damping (c ≠ 0), as shown in Fig. 1.14. The equation of motion for this system with the same force as before, F = F0 sin ωt, is:

RTENOTITLE....................(1.86)

The solution has two parts, as before. It is similar to the last example, except for an additional damping term. The particular solution is solved similarly to the last example in Eqs. 1.76,1.77, and 1.78. Differentiating and substituting into Eq. 1.86 gives:

RTENOTITLE
RTENOTITLE....................(1.87)

Rearranging gives:

RTENOTITLE
RTENOTITLE....................(1.88)

Equating coefficients as before yields:

RTENOTITLE....................(1.89)

and

RTENOTITLE....................(1.90)

Solving for the constants gives:

RTENOTITLE....................(1.91)

and

RTENOTITLE....................(1.92)

which gives the particular solution:

RTENOTITLE
RTENOTITLE....................(1.93)

The total solution is homogenous and the particular solutions added together. In this case, if the same initial and boundary conditions are applied as before, the homogenous solution is the same as in the free-vibration case and is Eq. 1.69. The particular solution is Eq. 1.93 and is the effect of the forcing function on the system. As noted before, the critical damping coefficient dictates the behavior of the homogenous part of the solution.


Ex. 1.6 is an SDOF forced and underdamped CLO system.

Two-Degree-of-Freedom (2DOF) Systems


At the basic level, systems with two and more DOFs are similar to SDOF systems. In a 2DOF system, two independent equations of motion are required to define a system (e.g., a double CLO or a double pendulum), but the DOFs 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.

Free Undamped 2DOF System

As in the SDOF system, an undamped (c1 and c2 = 0) system will be developed first. Fig. 1.15 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.94)

and for mass 2 is

RTENOTITLE....................(1.95)


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....................(1.96)

and

RTENOTITLE....................(1.97)

for which the second differentiation is

RTENOTITLE....................(1.98)

and

RTENOTITLE....................(1.99)

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

RTENOTITLE....................(1.100)

and

RTENOTITLE....................(1.101)

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....................(1.102)

Using linear algebra, the determinate is:

RTENOTITLE....................(1.103)

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

RTENOTITLE....................(1.104)

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....................(1.105)

and

RTENOTITLE....................(1.106)

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

RTENOTITLE....................(1.107)

therefore,

RTENOTITLE....................(1.108)

The first natural frequency motions then are:

RTENOTITLE....................(1.109)

and

RTENOTITLE....................(1.110)

The mode shape for the second natural frequency is:

RTENOTITLE....................(1.111)

therefore,

RTENOTITLE....................(1.112)

The second natural frequency motions then are:

RTENOTITLE....................(1.113)

and

RTENOTITLE....................(1.114)

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

RTENOTITLE
RTENOTITLE....................(1.115)

and

RTENOTITLE....................(1.116)

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

Ex. 1.7 is a 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. 1.16), but adding viscous dampers (c1 and c2 ≠ 0), the equations of motion for the independent DOFs are determined for mass 1 as

RTENOTITLE....................(1.127)

and for mass 2 as

RTENOTITLE....................(1.128)

This can be written in matrix form as

RTENOTITLE....................(1.129)

or as

RTENOTITLE....................(1.130)

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....................(1.131)

RTENOTITLE....................(1.132)

RTENOTITLE....................(1.133)

RTENOTITLE....................(1.134)

RTENOTITLE....................(1.135)

and

RTENOTITLE....................(1.136)

If the solution is assumed to be of the form

RTENOTITLE....................(1.137)

then

RTENOTITLE....................(1.138)

and

RTENOTITLE....................(1.139)

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

RTENOTITLE....................(1.140a)

which can be rewritten as

RTENOTITLE....................(1.140b)


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....................(1.141)

The determinate of Eq. 1.141 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. 1.8 is a 2DOF free and underdamped CLO system.


More information on linear algebra can be found in Chap. 2 of this section of the Handbook.

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.[11][12][13][14][15][16][17][18][19]

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.

Continuous Systems


If one continues to add DOFs, the limit at an infinite DOF defines a continuous system. The result becomes a PDE. The following is a brief description of the separation of variables method for solving a PDE.

Fig. 1.17 shows a freebody diagram for axial and torsional systems. The axial system equations will be used to determine the solution of the equations of motion. Eq. 1.162 is the axial equation of motion:

RTENOTITLE....................(1.162)

where RTENOTITLE = the inertial force, RTENOTITLE = the rate of strain change, mgc = the static weight of the element, and RTENOTITLE = the force from viscous damping. This PDE, Eq. 1.162, can be solved using the separation of variables method. This is shown as:

RTENOTITLE....................(1.164)


The following solution assumption is made concerning the time function:

RTENOTITLE....................(1.165)

This equation is substituted back into the assumed solution, which then is appropriately differentiated and substituted back into the equation of motion. The equation becomes

RTENOTITLE....................(1.166)

which is of the form

RTENOTITLE....................(1.167)

The standard solution of this equation is:

RTENOTITLE....................(1.168)

The constants of integration, C1 and C2, are determined by the initial and boundary conditions, and φ is a collection of the constants and is given by:

RTENOTITLE....................(1.169)

Therefore, the total solution is:

RTENOTITLE....................(1.170)

The solution to the torsional equation of motion is derived similarly to the axial equation, with the substitution of the appropriate variables and noting that there is no initial strain from gravity. The variables u, m, Ac, E, c, ω, vs, and φ are replaced by θ, I, J, G, cθ, ωθ, vθ, and η, respectively. The torsional equation of motion is:

RTENOTITLE....................(1.163)

This gives the solution as:

RTENOTITLE....................(1.171)

Constants C1 and C2 are based on the initial and boundary conditions, and η is a collection of the constants and is given by:

RTENOTITLE....................(1.172)

where

RTENOTITLE....................(1.173)

and

RTENOTITLE....................(1.174)

Wave Reflection From Various Geometric Boundaries


As discussed early in the chapter, wave propagation is the movement of a distinct group of waves through some material in response to an external force.

A key point in wave-propagation studies is how waves interact with geometric discontinuities. What happens as a wave meets a fixed or free boundary condition? Also, what happens to a wave as it encounters a geometrical area change or a change in material properties?

There are two limiting boundary conditions for wave propagation: a fixed (pinned) end (zero displacement) and a free end (zero stress). A fixed end is a boundary condition in which there is zero displacement. According to wave theory, during a wave encounter with a fixed end, the stress at the fixed end doubles during the passage of the wave. A reflection of a stress wave will simply bounce back with the same sign. A compression wave will reflect as a compression wave and a tension wave will reflect as a tension wave. At a fixed end, because the displacement is zero, the particle velocity will be zero. The wave particle velocity amplitude is inverted during a reflection from a fixed end.

A free end is defined as a boundary condition free to move. The stress at the free end is always zero. The effects on stress and particle velocity caused by a free end are opposite of the effects on stress on a fixed end. A compression wave encountering a free end reflects as a tension wave, and a tension wave reflects as a compression wave. The wave particle velocity values double during an encounter with a free end and reflect with the same sign.

As a wave encounters a change in cross-sectional area, some of the wave is reflected and some is transmitted (refracted). The amplitudes and sign of the waves depend on the relative change in cross-sectional area. The equation that describes the effect on the incident force, Fi, of a cross-sectional-area, density, or modulus-of-elasticity change for the transmitted wave is

RTENOTITLE....................(1.175)

and for the reflected wave is

RTENOTITLE....................(1.176)

If an incident wave encounters a junction where the relative change in cross-sectional area is greater than 1 (a smaller area to a larger area), most of the wave will transmit through the junction. Some of the wave will reflect from the junction and will keep the same sign. For example, a compression wave will transmit through the junction and keep going as a somewhat-diminished compression wave. The part of the wave that is reflected is still a compression wave, but its amplitude is less than that of the wave that transmitted though the junction.

On the other hand, if an incident wave encounters a junction where the relative change in cross-sectional area is less than 1 (a larger area to a smaller area), most of the wave will reflect off the junction, but some of it will transmit through the junction and will keep the same sign. For example, a compression wave will transmit through the junction and keep going as a diminished compression wave. The reflected part of the wave is a tension wave whose absolute amplitude is greater than that of the compression wave that is transmitted through the junction.

As with most drillstrings, there are many geometric discontinuities (changes in cross-sectional area) that will cause part of the wave to refract and part to reflect. For example, drill collars to heavyweight drillpipe to drillpipe all are geometric discontinuities. Sometimes, too, there are material discontinuities—changes in material density or modulus of elasticity—that cause refractions and reflections. A third possible type of discontinuity is when there are different endpoints. For example, if the pipe is stuck, one end can be modeled as stuck. If the pipe is hanging freely, such as with casing running, then the end is free.

More-detailed information on wave propagation can be found in several sources.[20][21][22][23][24][25][26][27]

Nomenclature


a = acceleration, L/t2, ft/sec2
A = slope, dimensionless
Ac = cross-sectional area, L2, in.2
c = axial damping coefficient, mL/t, lbf-ft/sec
ccrit = critical damping coefficient, dimensionless
cθ = torsional damping coefficient, mL/t, lbf-sec/rad
C = constant of integration, various
C = damping matrix, mL/t, lbf-ft/sec
di = inner diameter, L, in.
do = outer diameter, L, in.
E = modulus of elasticity, m/Lt2, psia
E(t) = excitation function
f = cyclic frequency, 1/t, cycle/sec
f(t) = function of t
f'(t) = first differential function
F = axial force, mL/t2, lbf
Fd = damping force, mL/t2, lbf
Ff = friction force, mL/t2, lbf
Fh = hysteretic force, mL/t2, lbf
Fi = incident force, mL/t2, lbf
Fn = normal force, mL/t2, lbf
Fr = reflected force, mL/t2, lbf
Ft = transmitted force, mL/t2, lbf
F0 = initial force, mL/t2, lbf
gc = gravitational constant, L/t2, 32.174 ft/sec2
g(t) = function of t
g'(t) = first derivative of function of g(t)
G = shear modulus, m/Lt2, psia
h = hysteretic factor, dimensionless
i = imaginary operator
I = second moment of inertia, L4, in.4
j = iteration index
J = polar moment, L3, in.3
k = spring constant, m/t2, lbf/in.
keq = equivalent spring constant, m/t2, lbf/in.
kw = wave number, 1/L, 1/ft
K = stiffness matrix, m/t2, lbf/in.
L = total length, L, ft
m = mass, m, lbm
mθ = mass polar moment of inertia, mL, lbf-sec2
M = mass matrix, m, lbm
n = exponent
P = generic value
Q = generic value
r = generic constant
rj = iterated generic constant
R(t) = response function
t = time, seconds
T = torque, mL2/t2, lbf-ft
T(t) = displacement function in terms of time, t
u = displacement, L, in.
U(x,t) = continuous displacement function, L, in.
v = velocity, L/t, ft/sec
vg = group speed, L/t, ft/sec
vs = sonic velocity, L/t, ft/sec
vw = phase velocity, L/t, ft/sec
v0 = initial velocity, L/t, ft/sec
vθ = torsional sonic velocity, L/t, ft/sec
x = displacement, L, in.
xh = homogeneous displacement, L, in.
xp = particular displacement, L, in.
x0 = initial displacement, L, in.
= first derivative with respect to time of displacement (velocity), L/t, ft/sec
p = particular velocity, L/t, ft/sec
= second derivative with respect to time of displacement (acceleration), L/t 2 , ft/sec 2
p = particular acceleration, L/t 2 , ft/sec 2
X(x) = displacement function in terms of location x
X = displacement vector
= velocity vector, L/t, ft/sec
= acceleration vector, L/t2, ft/sec2
y = dependent variable, various
yi = distance at point i, L, in.
Y = limit value, various
z = arbitrary constant, various
Δl = change in length, L, in.
Δt = change in time, t, seconds
Δy = change in dependent variable, various
α = amplitude, various
ε = strain, L/L, in./in.
η = convenient coefficient, 1/L, 1/ft
θ = twist, rad
Ӫ = second derivative with respect to time of twist (acceleration) rad/sec2
λ = mode shape, dimensionless
λω = wavelength, L, in.
μ = friction factor, dimensionless
ξ = critical damping ratio, dimensionless
ρ = density, m/L3, lbm/in.3
σ = stress, m/Lt2, psia
Φ = phase angle, rad
ω = frequency, 1/t, Hz
ωd = damped natural frequency, 1/t, Hz
ωf = forcing frequency, 1/t, Hz
ωn = natural frequency, 1/t, Hz
ωθ = twist natural frequency, 1/t, Hz
φ = convenient coefficient, 1/L, 1/ft
τ = period, t, seconds


References


  1. Fanchi, J.R. 1997. Math Refresher for Scientists and Engineers. New York: John Wiley & Sons.
  2. Leithold, L. 1972. The Calculus with Analytic Geometry. New York: Harper and Row.
  3. Bird, J.O. 2001. Newnes Engineering Mathematics Pocket Book, third edition. Oxford, UK: Newnes.
  4. Barré de Saint-Venant, A.J.C. 1867. Sur le choc longitudinal de deux barres élastiques. Journal de Mathématiques Pures et Appliquées 2 (12): 237.
  5. Love, A.E.H. 1926. A Treatise on the Mathematical Theory of Elasticity, fourth edition. New York: Dover Publications.
  6. Meyers, M.A. 1994. Dynamic Behavior of Materials. New York: John Wiley & Sons.
  7. Harris, C. and Crede, C. ed. 1988. Shock and Vibration Handbook, third edition. New York: McGraw-Hill Book Co.
  8. Dareing, D.W. and Livesay, B.J. 1968. Longitudinal and Angular Drill-String Vibrations With Damping. Journal of Engineering for Industry
  9. Kolski, H. 1963. Stress Waves in Solids. New York: Dover Publications.
  10. Doyle, J.F. 1991. Static And Dynamic Analysis Of Structures With an Emphasis on Mechanics and Computer Matrix Methods. Dordrecht, The Netherlands: Kluwer Academic Publishers.
  11. Achenbach, J.D. 1973. Wave Propagation in Elastic Solids, seventh edition. Amsterdam: North-Holland Publishing Co.
  12. Clough, R.W. and Penzien, J. 1975. Dynamics of Structures. New York: McGraw-Hill Book Co.
  13. Den Hartog, J.P. 1934. Mechanical Vibrations. New York: Dover Publications.
  14. Meirovitch, L. 1986. Elements of Vibration Analysis, second edition. New York: McGraw-Hill Book Co.
  15. Nashif, A.D., Jones, D.I.G., and Henderson, J.P. 1985. Vibration Damping. New York: John Wiley & Sons.
  16. Seto, W. 1964. Mechanical Vibrations. New York: McGraw-Hill Book Co.
  17. Shabana, A.A. 1991. Theory of Vibration, Volume I: An Introduction. New York: Springer-Verlag.
  18. Shabana, A.A. 1991. Theory of Vibration, Volume II: Discrete and Continuous Systems. New York: Springer-Verlag.
  19. Elmore, W.C. and Heald, M.A. 1969. Physics of Waves. New York: Dover Publications.
  20. Hudson, J.A. 1980. The Excitation and Propagation of Elastic Waves. Cambridge, UK: Cambridge University Press.
  21. Mal, A. and Singh, S.J. 1991. Deformation of Elastic Solids. Englewood Cliffs, New Jersey: Prentice-Hall.
  22. Sharman, R.V. 1963. Vibrations and Waves. London: Butterworth.
  23. Chin, W.C. 1994. Wave Propagation in Petroleum Engineering: Modern Applications to Drillstring Vibrations, Measurement-While-Drilling, Swab-Surge, and Geophysics. Houston, Texas: Gulf Publishing.
  24. Tolstoy, I. 1973. Wave Propagation. New York: McGraw-Hill Book Co.
  25. Doyle, J.F. 1989. Wave Propagation in Structures: An FFT-Based Spectral Analysis Methodology. New York: Springer-Verlag.
  26. Karl, J.H.: An Introduction To Digital Signal Processing, Academic Press Inc., San Diego, California (1989).
  27. Ramirez, R.W.: The FFT: Fundamentals and Concepts, Prentice-Hall, Englewood Cliffs, New Jersey (1985).

Suggested Reading


Aarrestad, T.V. and Kyllingstad, A. 1993. Measurements and Theoretical Models on Rig Suspension and the Effect on Drillstring Vibrations. SPE Drill & Compl 8 (3): 201-206. SPE-19553-PA. http://dx.doi.org/10.2118/19553-PA.

Aarrestad, T.V., Tonnesen, H.A., and Kyllingstad, A. 1986. Drillstring Vibrations: Comparison Between Theory and Experiments on a Full-Scale Research Drilling Rig. Presented at the SPE/IADC Drilling Conference, Dallas, Texas, 9-12 February 1986. SPE-14760-MS. http://dx.doi.org/10.2118/14760-MS.

Besaisow, A.A. and Payne, M.L. 1988. A Study of Excitation Mechanisms and Resonances Inducing Bottomhole-Assembly Vibrations. SPE Drill Eng 3 (1): 93-101. SPE-15560-PA. http://dx.doi.org/10.2118/15560-PA.

Booer, A.K. and Meehan, R.J. 1993. Drillstring Imaging - An Interpretation of Surface Drilling Vibrations. SPE Drill & Compl 8 (2): 93-98. SPE-23889-PA. http://dx.doi.org/10.2118/23889-PA.

Brett, J.F. 1992. The Genesis of Torsional Drillstring Vibrations. SPE Drill Eng 7 (3): 168-174. SPE-21943-PA. http://dx.doi.org/10.2118/21943-PA.

Broch, J.T. 1984. Mechanical Vibration and Shock Measurement. Denmark: Brüel & Kjær.

Dareing, D.W. 1984. Drill Collar Length is a Major Factor in Vibration Control. J Pet Technol 36 (4): 637–644. SPE-11228-PA. http://dx.doi.org/10.2118/11228-PA.

Doyle, J.F. 1988. A Spectrally Formulated Finite Element For Longitudinal Wave Propagation. The Intl. J. of Analytical and Experimental Modal Analysis 1 (January).

Doyle, J.F. 1986. Application Of The Fast-Fourier Transform (FFT) To Wave-Propagation Problems. The Intl. J. of Analytical and Experimental Modal Analysis 18 (October).

Dubinsky, V.S.H., Henneuse, H.P., and Kirkman, M.A. 1992. Surface Monitoring of Downhole Vibrations: Russian, European, and American Approaches. Presented at the European Petroleum Conference, Cannes, France, 16-18 November 1992. SPE-24969-MS. http://dx.doi.org/10.2118/24969-MS.

Dufeyte, M.-P. and Henneuse, H. 1991. Detection and Monitoring of the Slip-Stick Motion: Field Experiments. Presented at the SPE/IADC Drilling Conference, Amsterdam, 11–14 March. SPE-21945-MS. http://dx.doi.org/10.2118/21945-MS.

Dunayevsky, V.A., Abbassian, F., and Judzis, A. 1993. Dynamic Stability of Drillstrings Under Fluctuating Weight on Bit. SPE Drill & Compl 8 (2): 84–92. SPE-14329-PA. http://dx.doi.org/10.2118/14329-PA.

Eustes, A.W. 1996. A Frequency Domain Approach to Drillstring Jarring Analysis. PhD dissertation, Colorado School of Mines, Golden, Colorado.

Eustes, A.W., Long R.C., and Mitchell, B.J. 1995. Core Bit Frequency Signatures. Proc., ASME Energy-Sources Technology Conference & Exhibition (January) 17.

Eustes, I.A.W., Mitchell, B.J., and Stoner, M.S. 1994. Selection of Slim Hole Core Rods by Vibratory Analysis. J. Energy Resour. Technol. 116 (4): 251-257. http://dx.doi.org/10.1115/1.2906450.10.1115/1.2906450

Fertis, D.G. 1993. Nonlinear Mechanics. Boca Raton, Florida: CRC Press.

Fung, Y.C. 1965. Foundations of Solid Mechanics. Englewood Cliffs, New Jersey: Prentice-Hall.

Gieck, K. and Gieck, R. 1990. Engineering Formulas, sixth edition. New York City: McGraw-Hill Book Co. Inc.

Henneuse, H. 1992. Surface Detection of Vibrations and Drilling Optimization : Field Experience. Presented at the SPE/IADC Drilling Conference, New Orleans, Louisiana, 18-21 February 1992. SPE-23888-MS. http://dx.doi.org/10.2118/23888-MS.

Karasudhi, P. 1991. Foundations of Solid Mechanics. Dordrecht, The Netherlands: Kluwer Academic Publishers.

Kyllingstad, A. and Halsey, G.W. 1988. A Study of Slip/Stick Motion of the Bit. SPE Drill & Compl 3 (4): 369–373. SPE-16659-PA. http://dx.doi.org/10.2118/16659-PA.

Lai, W.M., Rubin, D., and Krempl, E. 1978. Introduction to Continuum Mechanics. Oxford, U.K.: Pergamon Press.

Macpherson, J. 1993. Vibration Signature of Roller Cone and PDC Bits. Baker Hughes INTEQ, Drilling Dynamics Technical Note (May).

Mitchell, B. 1995. Advanced Oilwell Drilling Engineering Handbook, tenth edition. Golden, Colorado: Mitchell Engineering.

Moon, F.C. 1987. Chaotic Vibrations: An Introduction for Applied Scientists and Engineers. New York City: John Wiley & Sons.

Paslay, P.R., Jan, Y.-M., Kingman, J.E.E. et al. 1992. Detection of BHA Lateral Resonances While Drilling With Surface Longitudinal and Torsional Sensors. Presented at the SPE Annual Technical Conference and Exhibition, Washington, D.C., 4-7 October 1992. SPE-24583-MS. http://dx.doi.org/10.2118/24583-MS.

Pestel, E. and Leckie, F.A. 1963. Matrix Methods in Elasto-Mechanics. New York City: McGraw-Hill Book Co. Inc.

Randall, R.B. 1987. Frequency Analysis. Denmark: Brüel & Kjær.

Rappold, K. 1993. Drilling Optimized with Surface Measurement of Downhole Vibrations. Oil & Gas J. 91 (7): 58.

Rappold, K. 1993. Drilling Dynamics—Conclusion Drillstring Vibration Measurements Detect Bit Stick-Slip. Oil & Gas J. (1 March): 66.

Reddy, J.N. 1984. An Introduction To The Finite Element Method. New York City: McGraw-Hill Book Co. Inc.

Skudrzyk, E. 1968. Simple and Complex Vibratory Systems. University Park, Pennsylvania: Pennsylvania State U. Press.

Timoshenko, S. and Gere, J. 1990. Mechanics of Materials, fourth edition, 1-912. Boston: PWS Publishing Co.

Young, W.C. 1989. Roark’s Formulas for Stress and Strain, sixth edition. New York City: McGraw-Hill Book Co. Inc.

Zienkiewicz, O.C. and Taylor, R.L. 1989. The Finite Element Method, Vol 1: Basic Formulation and Linear Problems, fourth edition. London, UK: McGraw-Hill Book Co. Inc.

SI Metric Conversion Factor


cycle/sec × 1.0* E + 00 = Hz
ft × 3.048* E – 01 = m
ft2 × 9.290 304* E – 02 = m2
in. × 2.54* E + 00 = cm
in.2 × 6.451 6* E + 00 = cm2
in.3 × 1.638 706 E + 01 = cm3
in.4 × 4.162 314 E – 07 = m4
lbf × 4.448 222 E + 00 = N
lbf-ft × 1.355 818 E + 00 = N•m
lbf-sec2 4.448 E + 00 = N•s2
lbm × 4.535 924 E – 01 = kg
psia × 6.894 757 E + 00 = kPa
rad × 1.00* E + 00 = rad


*

Conversion factor is exact.

Appendix – Glossary of Vibration Theory Terms


The period of a vibration is the time taken for a motion to repeat.

Frequency is the number of repeats per unit of time. This also is called a cyclic frequency. An angular frequency (sometimes called circular frequency) (ω) is measured in radians per unit time. Angular frequency is the product of 2π and the cyclic frequency.7

The amplitude of a frequency is the difference between the maximum or minimum values of a sinusoidally varying quantity and its mean. A wave’s crest is its maximum amplitude, and the trough is its minimum amplitude.

The wave period is the length of time is takes a wave to travel from crest to crest. The wave frequency is the inverse of the wave period.

The wavelength is the distance covered by one wave period.

The phase angle (abbreviated to "phase") of a wave is the fraction of the distance the wave has traveled relative to an arbitrary reference. The points at which the various waves cross the average amplitude line do not coincide unless the phase is zero. The phase is negative if the crossing occurs before the reference wave crossing, and positive if it occurs after the reference wave crossing.

Phase velocity is the velocity of a given point, as defined by the phase angle on a sinusoidal wave. When the phase velocity and the frequency remain constant, it is a nondispersing wave. If the phase velocity and the frequency are not constant, the wave disperses.

The wave number is the ratio of the wavelength to 2π.10 If this ratio is a real number, the wave will be a nondispersing wave (i.e., it will have a constant phase velocity for all its components and will retain its shape). If the wave number is complex, then the wave will be a dispersing wave (i.e., its components will have different phase velocities and it will not retain its shape and will disperse over time). All waves encountered in macroscopic reality are dispersing to some extent.

The group speed is the response of a wave’s group of sinusoidal components, and is different from phase velocity. Group speed is important; the interaction of all the wave components in terms of wave number, amplitude, and frequency causes the appearance of a carrier wave of some wave number, amplitude, and frequency. This carrier wave is modulated by group waves that propagate at group speeds. In a nondispersing wave, the group speeds are the same as the phase velocity of the wave components. In a dispersing wave, however, these group speeds can be different. If the group speeds are greater than the phase velocities, it will appear that a carrier wave originates at the back of the wave, propagates forward until it reaches the front of the wave, vanishes, and reappears at the back of the wave. If the group speeds are less than the phase velocities, it will appear that a carrier wave originates at the front of the wave, propagates backward until it reaches the back of the wave, vanishes, and reappears at the front of the wave. This can be visualized by thinking of a caterpillar’s motion. The caterpillar’s ripples appear to move backward, yet the caterpillar moves forward.24 This is one reason the analysis of dispersing waves can be very complicated.

The natural frequency of a system is the frequency at which the system vibrates when free of any friction or forcing functions.

A damped natural frequency is a natural frequency with friction.

Damping is the dissipation of energy with time or distance. Several forms of damping are discussed in this chapter.

A harmonic is a frequency that is an integer multiple of a given frequency.

A mode is a particular harmonic frequency.

The mode shape is a pattern of the system assumed during a harmonic.

A frequency spectrum is a band of frequencies.

A forcing function is an external force that is acting on a vibrating system. This has the effect of modifying the vibrations that are experienced in the system, sometimes significantly. For example, in a rotating-drillstring system, the drill bit will impart a forcing vibration to the lower end of the string.

A critical frequency is a frequency at which amplitudes become unbounded. This occurs when the frequency of the system matches one of the natural frequencies.

The critical rotary speed is a rotary speed that coincides with one of the critical frequencies of the system.

Resonance occurs when the frequency of the forcing function is equal to one of the natural frequencies of the system.

A stable system returns to its former position after the removal of the force that disturbed the system. In pipe work, stability is involved with buckling, bending, and straightness of the pipe.

An incident wave is a wave going into a geometric or material discontinuity.

A transmitted wave is the portion of a wave that passes through a geometric or material discontinuity.

A reflected wave is the portion of a wave that does not pass through a geometric or material discontinuity.

Impedance is the ratio of a force-like quantity to a velocity-like quantity when the arguments of the real or imaginary parts of the quantities linearly increase with time. The reciprocal of impedance is mobility. If the mechanical impedance of two materials and their geometries are equal, then the maximum amount of energy will be transmitted at their intersection.

The beating phenomenon occurs when two harmonic waves of slightly different frequencies are impressed on a body. They are a periodic variation in vibration at a frequency that is the difference between to two frequencies.

Self-excited vibrations occur when the forcing function is a function of the displacement, velocity, or acceleration of the system mass. If the energy gained from the self-excited vibrations is more than the system can absorb, the system will self-destruct.

A standing wave is a wave that appears to be standing still. There is no discernible beginning or end in time to the wave.

A transient wave is a wave that has a definite beginning and eventually dies out.

Strain energy is the energy gained by a mass as it is deformed.