Vibrations of Cantilever Beams

Vibrations of Cantilever Beams:
1 of 9
http://em-ntserver.unl.edu/Mechanics-Pages/Scott-Whitney/325hweb/Be...

Deflection, Frequency, and Research Uses
For: Dr. Negahban

EngrM 325H
Scott Whitney
April 23, 1999
Introduction
Measurements of thin film properties are difficult when compared to bulk materials. One method for finding
the modulus of elasticity of a thin film is from frequency analysis of a cantilever beam. A straight, horizontal
cantilever beam under a vertical load will deform into a curve. When this force is removed, the beam will
return to its original shape; however, its inertia will keep the beam in motion. Thus, the beam will vibrate at
its characteristic frequencies. If a thin film is sputtered onto the beam, the flexural rigidity will be altered.
11/12/2006 4:17 PM
2 of 9
This change causes the frequency of vibrations to shift. If the frequency shift is measured, the films elastic
modulus can be calculated.
Thin Cantilever Beam Setup

Beams studied in this paper are long, thin, cantilever beams. Figure 1, below, shows such a beam. One end of
the beam is fixed, while the other end is free. The origin of the coordinate axis is at the fixed end, point A. A
typical beam, used in this study, is L = 30 mm long, w = 5 mm wide, and t = 0.5 mm thick. The beam material
must be chosen so that its stiffness differs from the thin films stiffness so that the frequency shift is
significant. This study involved zinc oxide films (about 5 m thick), so a good substrate is aluminum.
Figure 1: Typical cantilever beam studied
Deflection of Cantilever Beam

If the free end of a cantilever beam is subjected to a point load, P, the beam will deflect into a curve. See
Figure 2 below. The larger the load, the greater the deflection, (x).
Figure 2: Cantilever beam deflection under load at fixed end

Assuming the beam undergoes small deflections, is in the linearly elastic region, and has a uniform
cross-section, the following equations can be used (Gere, p. 602).
The curvature of the beam, , is equal to the second derivative of the deflection
The curvature can also be related to the bending moment, M, and the flexural rigidity, EI,
11/12/2006 4:17 PM
3 of 9
where E is the elastic modulus of the beam and I is the moment of inertia. The bending moment in a beam
can be related to the shear force, V, and the lateral load, q, on the beam. Thus,
(1a,b,c)
For the load shown in Figure 2, the distributed load, shear force, and bending moment are:
Thus, the solution to Equation (1a) is
(2a)
At the free end of the beam, the displacement is:
(2b)
Vibrations of Beams
-Determination of Equations
When the force, P, is removed from a displaced beam, the beam will return to its original shape. However,
inertia of the beam will cause the beam to vibrate around that initial location. Assuming the elastic modulus,
inertia, and cross sectional area (A) are constant along the beam length, the equation for that vibration is
(Volterra, p. 310)
(3)
where
is the linear mass density of the beam. Equation (3) is best solved by separation of variables
(Atkins, p. A29). Assume that the displacement can be separated into two parts, one depends on position and
another on time.
(x,t) = X(x)f(t)
(4)
where X is independent of time, and f is independent of position. Then equation (3) becomes
Dividing equation (3) by
X(x)f(t):
11/12/2006 4:17 PM
4 of 9
(5)
Since the left side of equation (5) does not change as t varies, the right side must be a constant. Similarly
since the right half of equation (5) does not change as x varies, the left half must be a constant. Because each
side is constant, equation (4) is valid and the method of separation of variables may be used. Let this constant
be denoted
. It can be shown that
is a real quantity, and that
are natural frequencies of the beam.
Equation (5) can now be written as two differential equations (Volterra, p. 311),
(6a,b)
where
(6c)
In order to solve equation (6a), the following boundary conditions for a cantilever beam are needed
These boundary conditions come from the supports of a cantilever beam. The fixed end must have zero
displacement and zero slope due to the clamp. The free end cannot have a bending moment or a shearing
force.
The general solution to equation (6a) is a linear combination of trigonometric equations (Volterra, p.312)
(7)
Using the first boundary condition, C1 may be found
=>
Taking the first derivative of equation (7) and using the second boundary condition, C3 may be found
=>
The higher derivatives of equation (7), and remaining boundary conditions give
11/12/2006 4:17 PM
5 of 9
(8a)
(8b)
Equations (8a) and (8b) can be combined to give
(9)
Thus for a cantilever beam, equation (7) reduces to
(10)
According to Volterra, p. 312, the constants Cn are arbitrary. However, in order for the dynamic solution for
the displacement to be equal to the static solution (at time t=0), C2 must be equal to . With this value for
C2, Xn(0)=0 and Xn(L)=1.
Plugging equation (9) into either (8a) or (8b) will lead to the frequency equation for a cantilever beam,
(11)
The frequency equation can be solved for the constants, knL; the first six are shown below in Figure 3 (note,
kn=0 is ignored since it implies that the bar is at rest because =0). These constants along with equation (6c)
can be used to find the natural frequencies of a cantilever beam.
Figure 3: Constants for a cantilever beam vibrations. Note, since cosh(x) is large when x is large, knL
needs to be found with high precision.
For each frequency, there exists a characteristic vibration (Volterra, p. 319)
(12)
where An depends on the initial position at time t=0, and Bn depends on the initial velocity. In this study, the
beam starts its vibration when displaced and at rest. Thus Bn=0 and
(13)
The initial displacement, (x,0), was found above, equation (2a). Equation (13) can be solved analytically by
11/12/2006 4:17 PM
6 of 9
a computer math program,
-Examples of Vibration Modes

Typical values of An and frequency (for beams used in this study) are shown in Figure 3, above. Notice that
only the first few modes of vibration have significantly large values for the constant An. Thus, the higher
order vibration modes can be ignored. To get the total displacement of the cantilever beam, simply add all the
displacements found in equation (12) for each mode. Below, Figure 4 shows the displacement caused by each
mode at t=0; also included is the total initial displacement of the beam. From these displacements it is
apparent that even the third characteristic mode has little effect on the total displacement of the beam.
Figure 4: Initial displacement caused by each mode.

As time progresses, each mode will vibrate around the zero displacement line with the frequency listed in
Figure 3, above. Figure 5, below, shows this vibration for the first two modes, higher modes act similarly.
The total beam motion is complex; each characteristic mode vibrates with a different size, shape, and
frequency.
The following files can be downloaded to view an animation of the first few modes and total vibration of a
cantilever beam. To view these, the Lotus Screencam Player must be downloaded: Scplayer.exe
Total motion: 325Htot.scm
1st mode:
325h1st.scm
2nd mode:
325h2nd.scm
3rd mode:
325h3rd.scm
4th mode:
325h4th.scm
Figure 6, below, shows a reading from a piezoelectric sensor located near the fixed end of the beam; this
sensor shows the interaction of each vibration mode. Also shown is the vibration predicted from equation
(12) at that location. The actual and predicted signals are nearly identical; the major difference between them
is damping (both viscous and Kelvin-Voigt) of the higher order oscillations in the actual signal.
11/12/2006 4:17 PM
7 of 9
Figure 5: Vibrations of first two characteristic modes; other modes vibrate similarly. Time, t, is in
microseconds.
Figure 6: ZnO signal from a vibrating cantilever beam compared to theoretical vibrations.
Measurement of Modulus of Elasticity for Thin Films

-Reason for Using Frequency Analysis
If a cantilever beam is sputter coated with a thin film, then the flexural rigidity will change. A change in
stiffness will directly affect the frequency of the beams vibrations. Thus, the elastic modulus of the film can
be determined from this frequency shift. See Figure 7 for a schematic drawing of a sputtered beam; the
subscript s refers to the substrate and f refers to the film. Note, the height of the film is greatly
exaggerated, since
.
11/12/2006 4:17 PM
8 of 9
Figure 7: Cantilever Beam Sputter Coated with Thin Film

Theoretically, the stiffness could be determined from a stationary beam under a load, P, see equations (2a,b).
Suppose the flexural rigidity was increased by adding a thin film to the beam. Then, assuming the applied
load was constant, the deflection at the free end will decrease. For a typical beam described above, this
displacement change is extremely small. For example, an aluminum beam with initial displacement (at the
free end) of 0.55 mm will have an initial displacement of 0.49 mm with a relatively thick (8 m) ZnO film.
The 0.06 mm change is too small to be measured with any conventional tools, and a thinner ZnO film (1 m)
will give an immeasurable difference in displacement. Thus, stationary measurements of the films modulus
of elasticity are difficult, if not impossible.
However, the same beam coated with ZnO will have its fundamental frequency, , shift from 329.7 Hz to
342.0 Hz. Any oscilloscope can be used to measure the beams vibrations, see Figure 6, above. If enough data
points are taken, a Fourier transform of the signal can measure frequency shifts even as small as 0.2 Hz.
Therefore, the frequency shift can then be known with high accuracy, and the films elastic modulus can be
calculated.
-Determination of Equations
All of the equations derived in the previous sections (1-13) can still be applied to a two component beam if
several of the constants are changed: EI is replaced with EsIs+EfIf and is corrected as shown below. The
assumptions stated previously must still be fulfilled: small deflections, linearly elastic, and uniform
cross-section. The latter assumption implies that the film thickness cannot vary along the beam (difficult to do
when sputter coating long beams).
The linear mass density of the coated beam is now
The center of mass as well as the inertia of the coated beam will shift
Equation (6c) is used to measure the modulus of elasticity of the film from the frequency shift.
(14)
(15)
11/12/2006 4:17 PM
9 of 9
Since the amplitudes of higher vibration modes are so small, see Figure 4, equation (15) is solved using the
fundamental frequency, n=1.
Conclusion
Due to the small size of thin films, conventional methods of measuring their properties often do not work.
These thin film properties may differ from the bulk material properties. Therefore, alternative measurement
methods must be developed. Vibration of cantilever beams is one of these methods. It is not limited to just
determining the modulus of elasticity; other useful information, such as piezoelectric constants, can be
determined from cantilever beams.
Works Cited
Atkins, P. W. Physical Chemistry. 5th ed. New York: W. H. Freeman and Company, 1994.
Gere, J. M., Timoshenko, S. P. Mechanics of Materials. 4th ed. Boston: PWS Publishing Company, 1997.
Voltera, E., Zachmanoglou, E. C. Dynamics of Vibrations. Columbus, Charles E. Merrill Books, Inc., 1965.
11/12/2006 4:17 PM

Vibrations of Cantilever Beams

Uploaded by

Copyright:

Available Formats

Vibrations of Cantilever Beams

Uploaded by

Document Information

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Vibrations of Cantilever Beams

Uploaded by

Copyright:

Available Formats

Vibrations of Cantilever Beams: