Effect of Axial Load Mode Shapes A N D Frequencies Beams
Effect of Axial Load Mode Shapes A N D Frequencies Beams
Effect of Axial Load Mode Shapes A N D Frequencies Beams
R=19760006440 2018-11-04T02:12:00+00:00Z
ON MODE SHAPES A N D
FREQUENCIES OF BEAMS
Fmncis J. Shaker
Lewis Research Center
Cleveland, Ohio 44135
.. ~ .
1. Report No.
.
NASA
.
TN D-8109
4. Title and Subtitle
I 2. Government Accession No. 3. Recipient's Catalog No.
5. Report Date
December 1975
EFFECT O F AXIAL LOAD ON MODE SHAPES
AND FREQUENCIES O F BEAMS 6. Performing Organization Code
_ _ - .
16. Abstract
An investigation of the effect of axial load on the natural frequencies and mode shapes of uniform
beams and of a cantilevered beam with a concentrated m a s s a t the tip is presented. Character
istic equations which yield the frequencies and mode shape functions for the various cases are
given. The solutions to these equations a r e presented by a series of graphs so that frequency a s
a function of axial load can readily be determined. The effect of axial load on the mode shapes
are also depicted by another series of graphs.
.
17. Key Words (Suggested by Author(~)) 18. Distribution Statement
Structural dynamics -
Unclassified unlimited
Beam vibration under preload STAR category 39
-~ .~
19. Security Classif. (of this report) 20. Security Classif. (of this page)
Unclassified Unclassified
EFFECT OF AXIAL LOAD ON MODE SHAPES AND FREQUENCIES OF BEAMS
by Francis J. Shaker
SUMMARY
An investigation o the effect of axial load on the natural ,,,equencies and mode
shapes of uniform beams with various types of boundary conditions and of a cantilevered
beam with a concentrated mass at the tip is presented. This investigation yielded ex
pressions for the mode shapes and characteristic equations for the various cases con
sidered. For the uniform beams the characteristic equations were solved either numer
ically or in closed form, and the results a r e presented by a series of graphs showing the
effect of preload on the different types of beams. The effect of axial load on the mode
shapes is also shown in graphical form for several different loading conditions. For the
cantilevered beam with a tip mass, two types of axial loads w e r e considered: In the
first case the axial load vector remained constant, and in the second case the load was
directed through the root of the beam at all times. The results of this portion of the in
vestigation a r e presented in graphs that show the effects of both tip-mass variation and
, axial load on the fundamental frequency of the system.
INTRODUCTION
SYMBOLS
-
k nondimensional axial load parameter, )/Pl/EI
-
kcr nondimensional axial load parameter corresponding to critical load,
1 length of beam
axial load
2
shearing force distribution
time
beam displacement
nondimensional beam displacement, v/Z
mode shape of nth mode
X lengthwise coordinate
-
X nondimensional lengthwise coordinate, x/Z
characteristic values, eq. (9)
characteristic values for nth mode
frequency parameter, ?a
nondimensional frequency parameter, P 1
nondimensional frequency parameter corresponding to nth mode of
vibration, pwn
Superscript:
?
differentiation with respect to x
THEORETICAL ANALYSIS
(la)
From elementary beam theory the moment-curvature relationship for the sign convention
depicted in figures 1is
M(x,t) = E1 a2v(x, t)
n
ax
I
av(x,t)
(2 )
Finally, from equations (la) and (2) the equation of motion in terms of the displacement
variable can be written as
Equation (3) represents the governing partial differential equation describing the motion
of the beam, and equations (IC)and (2) represent the moment and shear distribution
along the beam.
Solution of equation and boundary conditions. - The time variable in equations (1) to
~
I
M(x,t) = m(x) sin Wt
4 2
a 4
)
+ k2 9- p4v(x) =0
dx dx2
m(x) = E1
d2v(x)
dx2
q(x) = -E1
d3v(x)
-1
2 dv(x)
+ k dx
where k2 = P/EI and p4 = pd2/EI. It will be convenient at this point to put equa
tion (5) into nondimensional form by letting x=
x/l, v=
v/l, k = k l , and = Pi?. Usp
ing these expressions allows equations (5) and (6) to become
-
T""(x) + k-2-ttv (x) - -4-
p V(X) =0 (7)
- EI-tt
m(x) = -V
1
The solution to equation (7), which can be readily verified by direct substitution, is
given by
I 5
(7
1
al= -- + / p- +) pl ' 2
CY2= -+ {-y-+ P
(10)
The frequency equation can now be determined from equation (9) once the boundary con
ditions are stipulated. The boundary conditions being considered at either end of the
Pinned
v =0 and m = 0 o r v f f = O
Clamped
-
v = O and v ' = O
Free
m = 0 or 5''= 0 and q = 0 or 7'" + -k 2-1v = 0
Guided -
The various combinations of these boundary conditions at both ends of the beam a r e
listed in table I. If the boundary conditions given by equations (11) a r e used at X = 0,
two of the unknown constants in equation (9) can be evaluated, and the resulting solutions
are for the beam free at E = 0
a,? +
@2 sin a2%)
( 12a)
"1
Note that each of equations (12) contains two unknowns. When the boundary conditions at
-
x = 1 are specified, these equations will yield two homogeneous, algebraic equations in
these two unknowns. The determinant of the coefficients of the unknowns must be equal
to zero for a nontrivial solution to exist. Equating this determinant to zero then yields
the characteristic equation from which the frequency parameter p can be determined a s
a function of axial load P or axial load parameter E. To illustrate, consider the case
where the beam is free at x
= 1. From equation ( l l c ) the boundary conditions are
From equations (l2c) and (13) the following two homogeneous, algebraic equations in A
and B are obtained after some manipulation:
cos a l ) + B b ;
CY; sinh cyl - (144
Equating the determinant of the coefficients of A and B in equations (14) to zero, ex
panding, and simplifying will then yield the following characteristic equation for the free-
f r e e beam:
-
2p6(1 cosh CY^ -2-4(k
cos a2)+ k + 3z4) sinh cy1 sin a2 = 0
Similarly, for the other boundary conditions at x = 1 the following characteristic equa
tions are obtained: for the free-guided beam
3 3
2 sinh a1 cos
CY CY
2 + a1 cosh cy1 sin CY^ = 0 (1%)
3 3
-a2 cosh al sin a2 + a1 sinh al cos a2 = 0
-4
2P + (2p4 + E4)cosh a1 cos a2 - -2-2
p k sinh a1 sin ai = 0
-2
2p (1 - cosh a1 cos a2) - -2k sinh al sin a2 = 0
sin a2 = 0
Equations (15g) and (15i) yield the following expressions for the frequency parameter:
for the pinned-pinned or guided-guided beam
(16)
Equations (16) and (17) exhibit the main characteristics of all the frequency equations
given by equations (15). Namely, as P (or E) Increase, the frequencies decrease until
P equals the Euler buckling load Per of the beam. At this point the lowest elastic fre
quency becomes zero. For P greater than Per, the fundamental frequency becomes
complex, and the corresponding mode shape is unstable (ref. 6, p. 4 5 1 and ref. 8, p.
302). Thus, the maximum value for P of practical interest is Per. The values of
Per and Ecr for beams with the different boundary conditions are tabulated in table E.
To determine the frequency parameter for the beams with other boundary conditions,
equations (15) for any value of less than the cor
- can be solved numerically for
responding kcr .
For the case of a cantilevered beam with a concentrated mass at the free end, two
types of axial loads will be considered. In the first case (fig. 2(a)) the axial load will be
taken as a constant in both magnitude and direction s o that this case corresponds to the
uniform clamped-free beam case except that a concentrated mass is added to the free
end. In the second case (fig. 3(a)) the axial load is assumed to be directed through the
root or fixed end of the beam at all times. For this case the axial load is not constant
since it changes direction during the motion. This type of loading will be referred to as
a directed axial load. For both cases the deflection is given by (eq. (12c) since the
boundary condition at the fixed end is the same as for a uniform cantilevered beam.
Also, the moment boundary condition at the free end is the same for both cases and is
given by
?'(I) =0
The difference between the two cases will be in the shear boundary condition at the free
end. For the case of constant axiai load, the shear boundary conditions can be found by
satisfging the equilibrium condition a t the free end which is depicted in figure 2@).
Thus
F"'(I) -2-1
-i- k
-4-
v (1)t qTP ~ ( 1=) 0 (20)
where qT is the ratio of the tip m a s s AT to total beam mass, AB = pQ. From equa
tions (12c), (18), and (20) the following two homogeneous equations in A and B are
obtained :
4cos a 2 ) +B ( 4
A{g sinh a1
A($ cosh al +
-zp4 + -2-2
p k sinh a1 sin a2 - (2B4 + E4)cosh a1 cos a2
For the case of a cantilevered beam with a tip mass and a directed axial load, the con
dition a t the free end is shown in figure 3@). From this figure the equilibrium condition
can be written as
(23)
Comparing this equation with equation (20) shows that the term 77TF4 - ii2 replaces
7Tp4. Thus multiplying equation ( 2 2 ) by p2 and replacing qTp4 with (qTp4 E 2 ) yields -
the characteristic equation for the beam with a directed axial load.
+ (qT? -x2)(4
+ ai)(al cosh al sin a2 - a2 sinh a1 cos a2)= 0
The mode shapes for a beam under preload can be determined from equations (12)
and one other boundary condition. For example, assume that the nth root of the fre
quency equation is known and is given by Fn.
Corresponding to there will be an &
10
h
cyln and a!
2n
which a r e obtained by substituting on
into equations (10). For the case
of a free-free beam, the displacement in the nth mode is given by equation (12a) as
a!
2n
2n
For the free-free beam case either of equations (14) can be used to determine a relation
between An and Bn. Using equation (14a) yields
Bn = -An F X )
(cash oln - COS CY )
(26)
@ln
Substituting equations (26) into (25) and arbitrarily setting An = 1 will then yield the
desired mode shape or eigenfunction for the free-free beam. In a similar manner the
mode shapes for the other boundary conditions can be determined. These functions are
summarized in the following equations: for free-free and free-pinned beams
2
C'*X
-+
@ln
n 2
CY
2n
11
- - -
vn(x) = (cash -
C Y ~ XCOS C Y ~ ~ X )
- - -
v,(x) = cosh -
C Y ~ Xcos C Y - @ln
~ 2n ~
-sin
CY C Y (27d)
~ ~
12
Vn@ =
It can readily be shown that the mode shapes for the cantilevered beams with a tip mass
a r e the same as those given by equation ( 2 7 ~ ) .
To determine the effect of axial load on the natural frequencies of a uniform beam,
the characteristic equations (given by eqs. (15a) to (15g) were solved numerically by the
method of bisection (ref. 9). In this investigation the axial load ratio P/Pcr was varied
between -1 and 1, and the first three lowest frequencies were determined. The results
for the fundamental elastic body frequencies a r e presented in figure 4. The pinned-
pinned, guided-guided, and pinned-guided curves in this figure were developed from
equations (16) and (17). In this figure the curves for the free-free and clamped-clamped
beams appear a s a single curve. Although the frequencies for these two cases were not
identically equal over the entire range of P/Pcr, the differences could not be shown for
the scale used in the figure. The same is true for the free-guided and clamped-guided
cases. The variation in frequency a s a function of axial load for the first few lower f r e
quencies for each uniform beam case is given in figure 5. It can be noted that the var
iation in the second and third frequency for all cases shown is almost linear for the
range of P/Pcr considered. These higher mode frequencies will go to zero when the
second and third buckling loads a r e reached. However, these portions of the curves a r e
of no practical interest unless the beam is prevented from buckling in its fundamental
mode.
To show the effect of axial loading on the mode shapes of the beam, equations (27)
were used to determine the modes for P/Pcr = -1.0, 0, and 0. 8. These calculations
were normalized, and the results a r e shown in figures 6 to 12 for the first three modes.
These figures indicate that the effect of axial load is greatest on the fundamental mode
and that the effect decreases rapidly as the mode number increases. The most
13
pronounced effect occurs for the free-pinned beam i n its fundamental mode (fig. 8(a)).
Note, however, that all the cases exhibit this effect as P/PCr approaches 1.
For a cantilevered beam with a tip m a s s equation (22) was solved for the case of a
constant axial load and equation (24) for an axial load directed through the root of the
beam. The method of solution of these equations was the same as that used for the uni
form beams. The effects of both tip mass ratio, qT = JT/.AB, and axial load ratio
P/PCr on the fundamental frequencies were determined, and the results a r e shown in
figures 13 to 16. Figure 13 shows the variation in the frequency parameters g; iK
as a function of tip mass ratio, qT, for values of P/PCr between -1 and 0. 8. Note that
8 4 s equals w1 d W E I , so that the frequency parameter for this case
is related to the total mass of the system and not just the beam mass. Figure 14 shows
the same data plotted as a function of axial load ratio with the tip mass ratio as a param
eter. Figures 15 and 16 are a similar set of graphs for the cantilever with a directed
axial load.
REFERENCES
14
6. Meirovitch, Leonard: Analytical Methods in Vibrations. The MacMillan Co., 1967.
7. Bisplinghoff, Raymond L. ; Ashley, Holt; and Hoffman, Robert L. : Aeroelasticity.
Addison-Wesley Publishing Co. , 1955.
8. Tong, Kin N. : Theory of Mechanical Vibrations. John Wiley & Sons., Inc., 1960.
9. Ralston, Anthony: A First Course in Numerical Analysis. McGraw-Hill Book Co.,
?
Inc., 1965.
15
TABLE I. - BOUNDARY CONDITIONS FOR UNIFORM BEAM UNDER CONSTANT AXIAL LOAD
- -~ .. ~ .
Free-free
-r Z
p
0 1
P ."(l) = 0
T'l'(1) + - E 2 7 (1) = I
F r e e - guided p PE 7'(0) = 0
-
?(l) = 0
v"'(0) +k%(O)= 0 v"'(1) + k % ( l ) = c
Guided- guided
-
k i d e d - pinned P P v'(0) = 0 G(1) = 0
;y'yo) +i27(0) = 0 v"(1) = 0
-
Xamped-guided J(0) = 0 v'(1) = 0
I
v'(0) = 0 v"'(1) + k 2 7 ( 1 ) = 0
~~ ~ .
-
Xamped- clamped J(0)= 0 v(1) = 0
v'(0) = 0 v'(1) = 0
16
-
Free-free sin kn = 0 nn n 2 E1
l2
-
-
Free guided COS kn = 0 (2n - 1)
-7r --
n 2 EI
2
l2
-
Free-pinned sin k = 0 nn n 2 E1 (Same as free-free)
l2
-
k i d e d - guided sin k =0 nn n 2 E1
l2 -
v1 = cos n Tz
-
- -
h i d e d - pinned COS k
n
=0 -n(2n
2
1) _
r2 _
E1
L2
>+I- v1 = cos n(x/2)
.
- (2n - 1) n 2 E1
Xamped-free COS k =0 -n --
2
l2
---
-
'inned- pinned sink = 0 nn n 2 E1 (Same as free-free)
..
l2
-
Ilamped-guided sink = 0
n nn n 2 E1
l2
-
k
lamped- clamped sin 2 = o 2nn 4n 2 -
E1
2
l2
Per
v1 = cos 27rz - 1
17
I
Ill
tv
...
^ l -
., Element 7
(a) Coordinate system f o r u n i f o r m vibrating beam (b) Sign convention f o r forces and moment
u n d e r constant axial load. acting o n a r b i t r a r y element.
Figure 1. - V i b r a t i n g beam w i t h arbitrary boundary conditions.
rv fli
(a) Coordinate system f o r cantilevered beam w i t h (b) Forces acting at free e n d of cantilevered beam
t i p mass u n d e r constant axial load. u n d e r constant axial load.
Figure 2. - Cantilevered beam w i t h t i p mass u n d e r constant axial load.
-
AV
- -*
X
(a) Coordinate system f o r cantilevered beam w i t h (b) Forces acting at free e n d of cantilevered beam
t i p mass u n d e r directed axial load. u n d e r directed axial load.
Figure 3. - Cantilevered beam w i t h t i p mass u n d e r directed axial load.
32 - Free-free o r clamped-clamped
--Free-pinned
- Clamped-pinned
---e
N
1 -
c1
i
aJ
-
EE 16
P
Y
L
-
-
----- --- ---_ ---- --_
--------
,-
4-
-------
- ----- ----
-1.0 -.8 -.6 -.4 -.2 0 .2 .4 .6 .8 1.0
Axial load ratio. PIP,
Figure 4 - Fundamental frequency as f u n c t i o n of axial load for u n i f o r m beam w i t h various boundary conditions.
19
-
- 02
6: 22.5
90
"i
15.0
12.5
40
P;
Lo.ol
7.5
20 h
2.5 -j
10
1 1 I I I I I
U
(b) For guided-pinned beam.
80 P3
140
6;
120
50
100
80
6;
60 -
10 - B:
20
0
-1.0 - , 8 -,6 -.4 -.2 0 .2 .4 .6 .8 1.0
0
-1.0 -.8 -.6 -.4 -.2 Axial load ratio, PIP,,
(di For free-guided beam.
(c)For free-free beam.
Figure 5. - Frequency as function of axial load
20
12or
I.
30
2 60, I$
v
3 20
* 40
10
20 B;
21
llllll I
-- -.a I
m
I - L I _ L L - I _ L _ l
E (a1 First mode.
-1. 0 - . 8 k - ‘ ” 1 1 1 I 1 I I I
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.0 0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.0
Nondimensional beam coordinate, x l l
22
it
,
I
'::h A xial
ratio,
load
PI Pc r
-_8
I l l l l l l l l u
(at First mode.
r
-.4
-. 6
-x
0 .1
I
.2
I
.3
I
.4
I
.5
I
.6
I
.7 .8
u .9 1.0 0
I I
.1
1.8'1
.2 .3
I
.4
I
,5 .6
I I
,7
Nondimensional beam coordinate, x l l
(bl Second mode. ( c l Third mode
Figure 7. - Effect of axial load on mode shape of free-guided beam.
23
-_2
-.4
-.6
- . a k l I: I ..I LI I 1-1
0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.0
Nondimensional beam coordinate, XU
(b) Second mode. (c)Third mode
Figure 8. - Effect of axial load on mode shape of free-pinned beam.
24
p
i,
1.0
.9
.6
.5
.4
.3
.2-
-
L
I I I
c
c
m
E
. 0' t l I I I I I \
la) First mode.
-.8
-1.0
0
coordinate, XU
(bl Second mode. (c)Third mode.
Figure 9. - Effect of axial load on mode shape of clamped-pinned beam.
25
Axial load
L
L k
26
beam
(bl Second mode. (c) Third mode.
Figure 11. - Effect of axial load on mode shape of clamped-clamped beam.
27
28
5. c
4.5
I\ to asb,;ea; of tip mass
Ratio 4'8L
4.4
Axial load
ratio,
P I Pc r
-1
-.6
-.4
-.2
0
.2
.4
.6
~ . .
Axial load ratio, PIPcr
L \
Figure 14. - Frequency as function of cantilever beam with tip mass.
I
0 1 2 3 4 5 6 7 8
Axial load
ratio,
L
a,
I
a,
1 -Q
N E
E
m
n
LL
L
0 I
lt I .
.8 ,
M
I
03
cp
w
w
- - . t-,
I . 1
N A T I O N A L AERONAUTICS A N D SPACE A D M I N I S T R A T I O N
WASHINGTON. D.C. 2 0 5 4 6 POSTAGE A N D FEES P A I D
N A T I O N A L AERONAUTICS A N D
O F F I C I A L BUSINESS SPACE ADMINISTRATION
451
PENALTY FOR PRIVATE USE $300 S P E C I A L FOURTH-CLASS R A T E USMAIL
BOOK
0 9 5 C U I Cl IJ D 7 5 1 1 2 8 S009g3DS
DEPT OF THE AIR FORCE
B F WEAPONS L A 3 O R A T O R Y
RTTN: T E C H N I C A L L I B R A R Y (SUL)
KIRTLAWD A P E N M 8 7 1 1 7