Bending Vibrations
Bending Vibrations
Bending Vibrations
com/
Structures
Journal of Intelligent Material Systems and
http://jim.sagepub.com/content/22/4/307
The online version of this article can be found at:
DOI: 10.1177/1045389X10379661
2011 22: 307 originally published online 11 August 2010 Journal of Intelligent Material Systems and Structures
Norman M. Wereley, Gang Wang and Anirban Chaudhuri
Demonstration of Uniform Cantilevered Beam Bending Vibration Using a Pair of Piezoelectric Actuators
Published by:
http://www.sagepublications.com
can be found at: Journal of Intelligent Material Systems and Structures Additional services and information for
http://jim.sagepub.com/subscriptions Subscriptions:
http://www.sagepub.com/journalsReprints.nav Reprints:
http://www.sagepub.com/journalsPermissions.nav Permissions:
http://jim.sagepub.com/content/22/4/307.refs.html Citations:
What is This?
1
n1
X
n
x T
n
t :
Power amplifier Function generator
+
PZT
Beam Strobe
light
Figure 2. Schematic of experimental setup.
Table 3. Dimensions and properties of beam and
piezoactuator.
Quantity Symbol Value Value
(SI Units) (English units)
Length of beam L
b
61.2775 cm 24.125 in.
Width of beam b
b
2.54cm 1.0 in.
Thickness of beam h
b
1.0583 cm 1/24 in.
Youngs modulus of Al E 69 GPa
Density of aluminum 2700 kg/m
3
Length of piezo L
c
5.08 in. 2.0 in.
Width of piezo b
c
2.54cm 1.0 in
Thickness of piezo h
c
0.03048 cm 12 mils
Youngs modulus
of piezoelectric material
E
c
62 GPa
Density of piezoelectric
material
c
7500 Kg/m
3
x0
0,
d
2
w x, t
d
2
x
xL
b
0 and
d
3
w x, t
d
3
x
xL
b
0:
Using the boundary conditions, the shape function for
the n-th mode is given by:
X
n
x cosh
n
x cos
n
x
n
sinh
n
x sin
n
x ,
2
where
n
is the solution of characteristic equation:
cosh
n
l cos
n
l 1, 3
and
n
sinh
n
l sin
n
l
cosh
n
l cos
n
l
: 4
For a uniform beam with rectangular cross-section, the
natural frequencies (in radians) of vibration are given
by !
n
n
EI=m
p
, n 1, 2, 3, . . . . The key issue here
is to solve the characteristic equation and find the
solution for
n
from (3). In this article, the function
fzero in MATLAB software was used for solving the
characteristic equation; the n-th natural frequency and
mode shape were calculated by substituting the corre-
sponding value of
n
. The computational code is given
in Table 4.
Numerical Methods
The numerical approaches to solve the beam trans-
verse bending vibration problem will be outlined in
this section. Instead of directly solving the governing
equation for beam bending vibration as shown in the
exact method, the beam bending vibration energies are
utilized (Lagrangian approach). Based on the varia-
tional principle, solving the boundary value problem
for the governing equation is equivalent to finding the
solution of for minimum system energy during bending
(Meirovitch, 1986; Inman, 1994). Therefore, we can
develop very useful approximate methods to solve the
dynamic problem, such as lumped parameter method,
assumed modes method and finite element method.
Figure 3. Picture of actual experimental setup.
Figure 4. Instantaneous pictures of bending mode shapes of a uni-
form cantilever beam vibrating. Excitation was provided by a pair of
piezoelectric actuators. A strobe light was used to capture the mode
shapes at maximum positive and negative amplitudes.
Demonstration of Uniform Cantilevered Beam Bending Vibration 311
at NANYANG TECH UNIV LIBRARY on September 7, 2014 jim.sagepub.com Downloaded from
These methods discretize a system represented by a
second-order differential equation:
M x Kx f x, t , 5
where M is the total mass matrix, K is the stiffness
matrix, and f(x, t) is applied force. Vector x is the
unknown displacement. For the free vibration problem
( f x, t 0), we need to solve a general eigenvalue and
eigenvector problem. Assuming a harmonic response
of the form x Xe
j!t
and substituting into the discre-
tized equation (4), we obtain the general eigenvalue
problem:
KX !
2
MX, 6
where ! is the natural frequency (in radian), and X is the
corresponding eigenvector. The mode shape function(s)
can be constructed from the eigenvector(s). For the uni-
form beam bending vibration problem, the application
of the lumped parameter, assumed modes method, and
finite element method is straightforward and well pre-
sented in the textbooks by Inman (1994) and Meirovitch
(1986). We will briefly present results of elemental mass
and stiffness matrices for each method; MATLAB codes
were developed for each method and are tabulated in
Tables 57.
(1) Lumped Parameter
The influence coefficients, C(i, j), and the corresponding
entries of the stiffness matrix, K(i, j), are calculated
using the following expressions:
C i, j
1
EI
x
j
x
2
i
2
x
2
i
6
_ _
for x
i
x
j
1
EI
x
i
x
2
j
2
x
2
j
6
_ _
for x
i
4x
j
_
_
:
K i, j C i, j
1
, M i, i m x
i
x
i1
(2) Assumed Modes
The entries of the stiffness and mass matrices are calcu-
lated from:
K i, j
_
l
0
EI x
d
2
X
i
x
dx
2
d
2
X
j
x
dx
2
dx,
M i, j
_
l
0
m x X
i
x X
j
x dx
where X
i
x is the assumed i-th mode shape function.
(3) Finite Element Methods
The element stiffness matrix and mass matrix are
(Inman, 1994):
K
EI
l
3
12 6l 12 6l
6l 4l
2
6l 2l
2
12 6l 12 6l
6l 2l
2
6l 4l
2
_
_
_
_
_
_
and
M
ml
420
156 22l 54 13l
22l 4l
2
13l 3l
2
54 13l 156 22l
13l 3l
2
22l 4l
2
_
_
_
_
_
_
:
Table 4. MATLAB code for exact solution of a uniform beam bending vibration.
% Physical properties
Lb = 24.125*0.0254; bb = 1*0.0254; hb = (1/24) *0.0254;
Eb = 69e9; rhob = 2700;
Lc = 2*0.0254; bc = 1*0.0254; hc = 12e-3*0.0254;
Ec = 62e9; rhoc = 7500;
Nmodes = 10; % number of modes desired
betainit = 0.5*pi : pi : (Nmodes-0.5)*pi; % initial estimates
dx = Lb / 100;
x = 0:dx:Lb; % points along beam length
for i = 1 : Nmodes,
BetaL = fzero(cos(x)*cosh(x)+1,betainit(i));
Beta(i) = BetaL / Lb;
al(i) = (sinh(BetaL)-sin(BetaL)) / (cosh(BetaL)+cos(BetaL));
xx = Beta(i)*x;
phi(:,i) = (cosh(xx)-cos(xx)) - al(i)*(sinh(xx)-sin(xx));
end
% frequencies from exact solution
fexact = Beta.*Beta*sqrt(EI/m)/(2*pi);
312 N. M. WERELEY ET AL.
at NANYANG TECH UNIV LIBRARY on September 7, 2014 jim.sagepub.com Downloaded from
TYPICAL RESULTS ASSUMING UNIFORM BEAM
Modal Frequencies
In this analysis, the effects of the PZT wafers are
neglected and their mass and stiffness contributions
are not considered. The beam is assumed to be uniform
and PZT only works as an exciter. Table 8 shows the
results from all the methods; the measured natural fre-
quencies agree well with our prediction results from the
exact solution, lumped parameter, assumed modes, and
finite element methods. Note that the results of the
exact solution and the assumed modes method are
exactly the same, which is because the exact mode
shapes of a cantilever beam were used for the numer-
ical method.
Table 6. MATLAB code for solution of uniform beam bending vibration using assumed modes method.
Nmodes = 10; % number of modes desired
betainit = 0.5*pi : pi : (Nmodes-0.5)*pi; % initial estimates
dx = Lb / 100;
x = 0:dx:Lb; % points along beam length
for i = 1 : Nmodes,
BetaL = fzero(cos(x)*cosh(x)+1,betainit(i));
Beta(i) = BetaL / Lb;
al(i) = (sinh(BetaL)-sin(BetaL)) / (cosh(BetaL)+cos(BetaL));
xx = Beta(i) * x;
phi(:,i) = (cosh(xx)-cos(xx)) - al(i)*(sinh(xx)-sin(xx));
% numerical differentiation
dphi(:,i) = [0 ; diff(phi(:,i))] / dx;
ddphi(:,i) = [0 ; diff(dphi(:,i))] / dx;
end
Kam = zeros(Nmodes); Mam = zeros(Nmodes); % initialisation
for i = 1 : Nmodes,
for j = 1 : Nmodes,% numerical integration
Kam(i,j) = trapz(x,EI*ddphi(:,i).*ddphi(:,j));
Mam(i,j) = trapz(x,m*phi(:,i).*phi(:,j));
end
end
[evecam,evalam] = eig(Kam,Mam);
% frequencies from assumed modes method
fam = sqrt(diag(evallp))/(2*pi);
Table 5. MATLAB code for solution of uniform beam bending vibration using lumped parameter method.
Nlump = 10; % number of lumps
Le = Lb/Nlump; % length of each lump
x = [Le/2:Le:Lb-Le/2]; % mid-point of each lump
Klp = zeros(Nlump); Clp = zeros(Nlump);
Mlp = m*Le*eye(Nlump);
for i = 1 : Nlump, % influence coefficients
Clp(i,i) = x(i)^3/EI/3;
for j = i+1 : Nlump,
Clp(i,j) = (x(j)-x(i)/3)*x(i)^2/EI/2;
Clp(j,i) = Clp(i,j);
end
end
Klp = inv(Clp);
[eveclp,evallp] = eig(Klp,Mlp);
% frequencies from lumped parameter method
flp = sqrt(diag(evallp))/(2*pi);
Demonstration of Uniform Cantilevered Beam Bending Vibration 313
at NANYANG TECH UNIV LIBRARY on September 7, 2014 jim.sagepub.com Downloaded from
Mode Shapes
The mode shape functions of beam bending vibra-
tion were extracted from the experimental measure-
ments. The data was then tabulated for each mode
and the resulting measured mode shapes were com-
pared to the mode shapes predicted using the
EulerBernoulli uniform beam analysis. The measured
and predicted mode shapes agree quite well, as shown
in Figure 5.
ACCOUNTING FOR PIEZOACTUATOR
In the above analyses, the inertia and stiffness effects
of the PZT patches were neglected. If the mass and
stiffness of PZTs are incorporated into our analysis,
then the beam is not uniform anymore and it is diffi-
cult to obtain the exact solutions. However, the
problem can be numerically solved if the region on
the beam with PZTs is treated as an element with mod-
ified stiffness and additional mass in the lumped
parameter, assumed modes, and finite element meth-
ods; the results are tabulated in Table 9. When com-
pared with the results in Table 8, the predicted natural
frequencies increase for each mode compared to the
uniform beam case, which imply that the added stiff-
ness effects of the PZTs are more dominant than the
mass or inertia contributions.
CONCLUSIONS
The cantilever beam test outlined in this article is
simple to construct and provides physical insight into
the actual manifestation of mode shapes and their
associated modal frequencies. The analysis of the
cantilevered beam is straightforward, and the analysis
Table 7. MATLAB code for solution of uniform beam bending vibration using finite element method.
Ne = 10; % number of elements
Le = Lb/Ne;
Ke = [12 6*Le -12 6*Le; 6*Le 4*Le^2 -6*Le 2*Le^2;
-12 -6*Le 12 -6*Le; 6*Le 2*Le^2 -6*Le 4*Le^2] / Le^3;
Me = [156 22*Le 54 -13*Le; 22*Le 4*Le^2 13*Le -3*Le^2;
54 13*Le 156 -22*Le; -13*Le -3*Le^2 -22*Le 4*Le^2] * Le/420;
Mg = zeros(2*Ne+2); Kg = zeros(2*Ne+2);
for k = 1 : Ne,
kbgn = 2*k-1; kend = 2*k+2;
Mg(kbgn:kend,kbgn:kend) = Mg(kbgn:kend,kbgn:kend) + m*Me;
Kg(kbgn:kend,kbgn:kend) = Kg(kbgn:kend,kbgn:kend) + EI*Ke;
end
Ktot = Kg(3:2*Ne+2,3:2*Ne+2);
Mtot = Mg(3:2*Ne+2,3:2*Ne+2);
[evecfem,evalfem] = eig(Ktot,Mtot);
% frequencies from fem
ffem = sqrt(diag(evallp))/(2*pi);
Table 8. Measured and predicted modal frequencies of a uniform beam neglecting effects of stiffness
and mass of the piezoelectric actuator (uniform beam approach).
Modal frequencies (Hz)
Mode number Measured
EulerBernoulli
uniform beam
Lumped
parameter
analysis
(10 elements)
Assumed
modes
analysis
(10 modes)
Finite element
analysis
(10 elements)
1 2.22 2.30 2.31 2.30 2.30
2 13.64 14.42 14.54 14.42 14.42
3 36.46 40.39 40.92 40.39 40.40
4 72.42 79.14 80.54 79.14 79.22
5 117.40 130.83 133.52 130.83 131.16
314 N. M. WERELEY ET AL.
at NANYANG TECH UNIV LIBRARY on September 7, 2014 jim.sagepub.com Downloaded from
can be easily made more challenging by introducing
additional effects into the analysis. For the undergrad-
uate with a typical structures course at the senior
level, the beam can be analyzed, neglecting the piezo-
actuators for simplicity, using the exact solution of the
EulerBernoulli uniform beam equation, and a finite
element analysis can also be developed and validated.
For the first year graduate student in an introductory
structural dynamics course, the mass and stiffness
loading of the piezoelectric actuator can also be
added to the analysis, and additional analysis meth-
ods, such as lumped parameter or assumed modes
analyses, can be attempted. Computer codes that can
be implemented in MATLAB have also been provided
in this article.
SMART BEAM TEACHING AID
The analyses in this article have been coded in
MATLAB and included in a CD-ROM teaching aid,
which can be provided along with a beam specimen as
described in this article, for easy application in the
classroom.
ACKNOWLEDGMENT
Supported by National Science Foundation (NSF)
Career Award to NMW, CAREER: Passive and
Active Damping Control for Large Civil Structures,
Award Number 9734244.
REFERENCES
Chen, P.C. and Chopra, I. 1996. Induced Strain Actuation of
Composite Beams and Rotor Blades with Embedded
Piezoceramics Elements, Smart Materials and Structures,
5:3548.
Chopra, I. 2002. Review of State of Art of Smart Structures and
Integrated Systems, AIAA Journal, 40:21452187.
Crawley, E.F. and Anderson, E.H. 1990. Detailed Models of Piezo-
ceramic Actuation of Beams, Journal of Intelligent Material
Systems and Structures, 1:425.
Crawley, E.F. and de Luis, J. 1987. Use of Piezoelectric Actuators as
Elements of Intelligent Structures, AIAA Journal,
25:13731385.
Epps, J.J. and Chopra, I. 2001. In-flight Tracking of Helicopter
Rotor Blades using Shape Memory Alloy Actuators, Smart
Materials and Structures, 10:104111.
f
1
(x/L)
f
2
(x/L)
f
3
(x/L)
f
4
(x/L)
f
5
(x/L)
1
0.5
0
1
0
1
1
0
1
1
0
1
1
0
1
Mode 1: w
1
=2.22 Hz
Mode 2: w
2
=13.64 Hz
Mode 3: w
3
=37.54 Hz
Mode 4: w
4
=72.42 Hz
Mode 5: w
5
=117.40 Hz
x/L
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Figure 5. Comparison of experimentally measured mode shapes (circle) with the exact solutions of the EulerBernoulli beam equation (solid
lines) for a uniform beam.
Table 9. Estimates of modal frequencies for a non-
uniform beam including piezoelectric actuator.
Modal frequencies (Hz)
Mode
number Measured
Lumped
parameter
analysis
(10 elements)
Assumed
modes
analysis
(10 modes)
Finite
element
analysis
(10 elements)
1 2.22 2.32 2.64 2.58
2 13.64 14.58 16.25 15.88
3 36.46 40.94 44.62 43.74
4 72.42 80.36 85.59 84.16
5 117.40 132.73 138.43 136.43
Demonstration of Uniform Cantilevered Beam Bending Vibration 315
at NANYANG TECH UNIV LIBRARY on September 7, 2014 jim.sagepub.com Downloaded from
Giurgiutiu, V., Chaudhry, Z. and Rogers, C.A. 1995. Engineering
Feasibility of Induced Strain Actuators for Rotor Blade Active
Vibration Control, Journal of Intelligent Material Systems and
Structures, 6:583597.
Inman, D.J. 1994. Engineering Vibrations, Prentice Hall, Englewood
Cliffs, NJ.
Kamath, G.M., Wereley, N.M. and Jolly, M.R. 1999.
Characterization of Magnetorheological Helicopter Lag
Dampers, Journal of the American Helicopter Society,
44:234248.
Kudva, J.N. 2004. Overview of the DARPA Smart Wing Project,
Journal of Intelligent Material Systems and Structures,
15:261267.
Meirovitch, L. 1986. Elements of Vibration Analysis, McGraw-Hill,
New York.
Niezrecki, C., Brei, D., Balakrishnan, S. and Moskalik, A. 2001.
Piezoelectric Actuation: State of the Art, The Shock and
Vibration Digest, 33:269280.
Tummers, B. 2006. DataThief III. Available at: http://datathief.org/
(accessed date July 2010).
Vishay Measurement Group. Surface Preparation for Strain Gage
Bonding, Application Note B-129-8, Available at: http://
www.vishay.com (accessed date July 2010).
Wereley, N.M. 2000. The Smart Beam Teaching Aid CD-ROM,
Department of Aerospace Engineering, University of Maryland.
316 N. M. WERELEY ET AL.
at NANYANG TECH UNIV LIBRARY on September 7, 2014 jim.sagepub.com Downloaded from