# On the Dynamic Analysis of a Beam Carrying Multiple Mass-Spring-Mass-Damper System

Authors: O. Barry; D. C. D. Oguamanam; J. W. Zu
Journal: Shock and Vibration
 (2014)

---

## Abstract

The exact natural frequencies, mode shapes, and the corresponding orthogonality relations are important in forced vibration analysis via modal expansion. In the present paper, a free vibration analysis is conducted to determine the exact natural frequencies and mode shapes of an axially loaded beam carrying several absorbers. An explicit expression is presented for the generalized orthogonality relations. These generalized orthogonality conditions are employed along with the assumed modes method to perform forced vibration analysis. The present approach is compared to other approximate methods in the literature with the classical orthogonality relations and different choice of mode shapes. The results indicate that the use of the generalized orthogonality relation with the exact mode shapes is required for a precise investigation of the dynamic response of a beam with mass-spring-mass-damper system.


---

## Body

## 1. Introduction

The study of beam vibrations is attracting continued interest all these years because of the wide range of engineering applications involving beams. Beam models are used to idealize buildings components, bridges, overhead transmission lines, microelectromechanical systems, and many more. In the case of overhead transmission lines a single conductor may be modeled as a beam that is subjected to a tensile load with the damper represented as a mass-spring-damper-mass system.Numerous authors have studied the free and/or forced response of beams carrying in-span mass or spring-mass systems. In most investigations on the free vibration, different authors have used analytical methods to determine the exact solutions for the natural frequencies and mode shapes (e.g., [1–13] and references mentioned therein).There are, however, few published exact solutions of the forced vibration problems. The most common analytical approach used in forced vibration analysis is based on the assumed modes method in conjunction with the classical orthogonality conditions which may be expressed as∫
0
L
Y
i
Y
j
d
x
=
δ
i
j, where L is the length of the beam, Y
i and Y
j are the ith and jth normal mode shapes, respectively, and δ
i
j is the Kronecker delta function. Hereinafter a beam without mass or spring-mass system is referred to as a bare beam, or else it is called a loaded beam.In [14–16], the authors employed the mode shapes of the bare beam and the classical orthogonality condition to study the forced response of loaded beams. Combinations of the mode shapes of the loaded beam with the classical orthogonality relation were employed in [17–23]. These two approaches are simple and easy to implement and, in some cases, give very good approximations. However, it was shown in [21, 24] that the mode shapes of the loaded beam and those of the bare beam are different. Further, Hassanpour et al. [25] showed that the mode shapes of the loaded beam were not orthonormal under the classical orthogonality conditions. They also presented the generalized orthogonality conditions for an axially loaded beam with in-span concentrated mass and springs. The study ignored the case where a suspended mass is connected to the beam via a spring-dashpot system.The majority of the studies reported in the literature focused on the study of beams with either an in-span mass or spring-mass system. A series combination of each type was found in some cases. However, works that simultaneously examined the free and forced vibration of a beam with multiple in-span mass-spring-damper-mass systems and subjected to an axial force were not found. The investigation of this problem is reported in the present paper; it is an extension of the work presented in [13]. The exact natural frequencies and mode shapes are determined and explicit expressions are presented for the orthogonality relations. Parametric studies are used to examine the effect of the magnitude and location of the in-span mass-spring-damper-mass system on the natural frequencies and vibrational response. The effect of the use of the classical orthogonality condition is also examined. The results are validated using the finite element method and results from the literature.


## 2. Equations of Motion

A schematic of a beam with several in-span mass-spring-damper-mass systems is depicted in Figure1. The system kinetic energy 𝒯 and potential energy 𝒱 may be expressed as
(1)
𝒯
=
∑
i
=
1
N
‍
(
1
2
m
∫
L
i
-
1
L
i
w
˙
i
2
(
x
,
t
)
d
x
+
1
2
M
d
i
z
˙
i
2
+
1
2
M
c
i
z
˙
0
i
2
)
,
𝒱
=
∑
i
=
1
N
‍
(
1
2
E
I
∫
L
i
-
1
L
i
w
i
′
′
2
(
x
,
t
)
d
x
+
1
2
c
d
i
(
z
˙
0
i
-
z
˙
i
)
2
+
1
2
k
i
(
z
0
i
-
z
i
)
2
+
1
2
T
∫
L
i
-
1
L
i
w
i
′
2
(
x
,
t
)
d
x
)
,
where z
0
i
(
t
)
=
w
i
(
L
i
,
t
), E
I is the beam flexural rigidity, m is the mass per unit length of the beam, and T is the beam pretension. The overdots and primes denote temporal and spatial derivatives, respectively. The subscript “i” denotes the position of the mass-spring-damper-mass system and “N” is the total number of attached mass-spring-damper-mass systems. The beam is assumed to be uniform and the pretension is constant.Figure 1
Schematic of an axial beam carrying several mass-spring-damper-mass systems.Introducing these energies equations into Hamilton’s principle and adding the forcing termF
(
x
,
t
) yield the equations of motion and continuity conditions:
(2)
E
I
w
i
′
′
′
′
+
m
w
¨
i
-
T
w
i
′
′
=
F
(
x
,
t
)
,
(3)
M
d
i
z
¨
+
k
i
(
z
i
-
z
0
i
)
+
c
d
i
(
z
˙
i
-
z
˙
0
i
)
=
0
,
(4)
w
i
(
L
i
,
t
)
=
w
i
+
1
(
L
i
,
t
)
,
(5)
w
i
′
(
L
i
,
t
)
=
w
i
+
1
′
(
L
i
,
t
)
,
(6)
w
i
′
′
(
L
i
,
t
)
=
w
i
+
1
′
′
(
L
i
,
t
)
-
M
c
i
z
¨
0
i
+
E
I
w
i
′
′
′
(
L
i
,
t
)
-
T
w
i
′
(
L
i
,
t
)
-
k
i
(
z
0
i
-
z
i
)
-
c
d
i
(
z
˙
0
i
-
z
˙
i
)
,
(7)
-
E
I
w
i
+
1
′
′
′
(
L
i
,
t
)
+
T
w
i
+
1
′
(
L
i
,
t
)
=
0
.


## 3. Frequency Equation and Mode Shapes

The deformationsw
i
(
x
i
,
t
) and displacements z
(
t
) are expressed as
(8)
w
i
(
x
i
,
t
)
=
L
W
i
(
ζ
)
e
ω
t
,
z
i
(
t
)
=
L
A
i
e
ω
t
,
where W
i
(
ζ
i
) and A
i are the respective nondimensional amplitudes of w
i
(
x
i
,
t
) and z
i
(
t
) and ω is the circular complex natural frequency of the system. Substituting (8) into (2)–(7) yields the following nondimensional system equations:
(9)
W
i
′
′
′
′
(
ζ
i
)
-
s
2
W
i
′
′
(
ζ
i
)
+
Ω
4
W
i
ζ
i
=
0
,
(10)
A
i
-
K
i
W
i
(
ξ
i
)
=
0
,
where
(11)
W
i
(
ζ
i
)
=
c
1
i
sin
⁡
α
ζ
i
+
c
2
i
cos
⁡
⁡
α
ζ
i
+
c
3
i
sinh
⁡
⁡
β
ζ
i
+
c
4
i
cosh
⁡
⁡
β
ζ
i
.
Substituting (8) into (3)–(7) yields
(12)
W
i
(
ξ
i
)
=
W
i
+
1
(
ξ
i
)
,
(13)
W
i
′
(
ξ
i
)
=
W
i
+
1
′
(
ξ
i
)
,
(14)
W
i
′
′
(
ξ
i
)
=
W
i
+
1
′
′
(
ξ
i
)
,
(15)
W
i
′
′
′
(
ξ
i
)
-
W
i
+
1
′
′
′
(
ξ
i
)
-
η
i
W
i
(
ξ
i
)
+
γ
i
A
i
=
0
,
where the following nondimensional variables are used:
(16)
ξ
i
=
L
i
L
,
A
i
=
K
W
i
(
ξ
i
)
,
K
i
=
k
i
+
ω
c
d
i
k
i
+
M
d
i
ω
2
+
ω
c
d
i
,
s
2
=
T
L
2
E
I
,
Ω
4
=
m
ω
2
E
I
L
4
,
η
i
=
k
i
+
ω
2
M
c
i
+
ω
c
d
i
E
I
L
3
,
γ
i
=
k
i
+
ω
c
d
i
E
I
L
3
,
ζ
=
x
L
,
α
=
-
s
2
2
+
s
4
4
+
Ω
4
,
β
=
s
2
2
+
s
4
4
+
Ω
4
.
Substituting (10) into (15) yields
(17)
W
i
′
′
′
(
ξ
i
)
-
W
i
′
′
′
(
ξ
i
)
+
W
i
(
ξ
i
)
(
K
i
γ
i
-
η
i
)
=
0
.
The use of any classical boundary conditions at each end of the beam along with (12)–(14) and (17) yields a set of 4
+
4
N algebraic homogeneous equations (4 equations from the boundary condition at the ends and 4
N equations from the continuity relations). These algebraic equations are linear in the unknown coefficients (C’s) and they can be presented in matrix format as
(18)
[
ℱ
]
(
4
+
4
N
)
X
(
4
+
4
N
)
{
C
}
(
4
+
4
N
)
X
(
1
)
=
{
0
}
(
4
+
4
N
)
X
(
1
)
,
where the elements of the matrix ℱ are listed in the Appendix. A nontrivial solution is obtained when matrix ℱ is singular. Hence, the characteristic or frequency equation is obtained as
(19)
det
⁡
⁡
(
[
ℱ
]
(
4
+
4
N
)
X
(
4
+
4
N
)
)
=
0
.
The mode shapes associated with each beam segment are obtained by substituting the integration constants from (18) into (11). An example of a matrix that yields the frequency equation is provided in the Appendix for the case of a cantilever beam with one in-span mass-spring-damper-mass support.


## 4. Generalized Orthogonality Relations

The solution of the equations of motion, (2) and (3), can be expressed as
(20)
w
i
(
x
,
t
)
=
Y
i
(
x
)
(
r
)
e
ω
t
,
z
i
(
t
)
=
Z
i
(
r
)
e
ω
t
,
where the superscript “r” denotes the mode number and
(21)
Y
i
(
x
)
=
L
W
i
(
ζ
)
,
Z
i
=
L
A
i
.
Substituting (20) into (3) and multiplying the resulting equation by Z
(
s
) yield
(22)
k
i
Z
i
(
r
)
Z
i
(
s
)
+
M
d
i
ω
r
2
Z
i
(
r
)
Z
i
(
s
)
+
ω
r
c
d
i
Z
i
(
r
)
Z
i
(
s
)
=
Y
i
*
(
r
)
Z
i
(
s
)
(
k
i
+
ω
r
c
d
i
)
,
where Y
i
*
(
r
)
=
Y
i
(
r
)
(
ξ
1
), Interchanging “r” and “s” in (22) and subtracting the resulting equation from (22) yield
(23)
(
ω
r
2
-
ω
s
2
)
M
d
i
Z
i
(
r
)
Z
i
(
s
)
+
c
d
i
(
ω
r
-
ω
s
)
Z
i
(
r
)
Z
i
(
s
)
=
k
i
(
Y
i
*
(
r
)
Z
i
(
s
)
-
Y
i
*
(
s
)
Z
i
(
r
)
)
+
c
d
i
(
ω
r
Y
i
*
(
r
)
Z
i
(
s
)
-
ω
s
Y
i
*
(
s
)
Z
i
(
r
)
)
.
With reference to the equations of motion of the beam, substituting (20) into (2), then multiplying the resulting equation by Y
(
s
), and integrating over the entire length of the beam, as well as applying the continuity conditions, (12)–(15), with any classical boundary conditions except those for free ends yield
(24)
ω
r
2
∑
i
=
1
N
‍
(
m
∫
0
L
i
Y
i
(
r
)
Y
i
(
s
)
d
x
+
M
c
i
Y
i
*
(
r
)
Y
i
*
(
s
)
)
+
ω
r
∑
i
=
1
N
‍
c
d
i
Y
i
*
(
r
)
Y
i
*
(
s
)
=
-
∑
i
=
1
N
‍
(
E
I
∫
0
L
i
Y
i
′
′
(
r
)
Y
i
′
′
(
s
)
d
x
k
k
k
k
k
k
k
k
+
T
∫
0
L
i
Y
i
′
(
r
)
Y
i
′
(
s
)
d
x
+
k
i
Y
i
*
(
r
)
Y
i
*
(
s
)
k
k
k
k
k
k
k
k
k
k
-
k
i
Z
i
(
r
)
Y
i
*
(
s
)
-
c
d
i
ω
r
Z
i
(
r
)
Y
i
*
(
s
)
∫
0
L
i
‍
)
.
Rewriting (24) and interchanging “r” and “s” yield
(25)
ω
s
2
∑
i
=
1
N
(
m
∫
0
L
i
Y
i
(
s
)
Y
i
(
r
)
d
x
+
M
c
i
Y
i
*
(
s
)
Y
i
*
(
r
)
)
+
ω
s
∑
i
=
1
N
c
d
i
Y
i
*
(
s
)
Y
i
*
(
r
)
=
-
∑
i
=
1
N
(
E
I
∫
0
L
i
Y
i
′
′
(
s
)
Y
i
′
′
(
r
)
d
x
k
k
k
k
k
k
k
k
+
T
∫
0
L
i
Y
i
′
(
s
)
Y
i
′
(
r
)
d
x
+
k
i
Y
i
*
(
r
)
Y
i
*
(
s
)
k
k
k
k
k
k
k
k
k
i
-
k
i
Z
i
(
s
)
Y
i
*
(
r
)
-
c
d
i
ω
s
Z
i
(
s
)
Y
i
*
(
r
)
∫
0
L
i
‍
)
.
Subtracting (25) from (24) and substituting (23) into the resulting equation yield
(26)
(
ω
r
2
-
ω
s
2
)
∑
i
=
1
N
‍
(
m
∫
0
L
i
Y
i
(
r
)
Y
i
(
s
)
d
x
+
M
c
i
Y
i
*
(
r
)
Y
i
*
(
s
)
k
k
k
k
k
k
k
k
k
k
i
+
M
d
i
Z
i
(
r
)
Z
i
(
s
)
∫
0
L
i
‍
)
+
(
ω
r
-
ω
s
)
∑
i
=
1
N
‍
c
d
i
(
Y
i
*
(
r
)
Y
i
*
(
s
)
+
Z
i
(
r
)
Z
i
(
s
)
)
=
0
.
From (26), the first set of orthogonality relation is obtained as
(27)
∑
i
=
1
N
‍
(
m
∫
0
L
i
Y
i
(
r
)
Y
i
(
s
)
d
x
+
M
c
i
Y
i
*
(
r
)
Y
i
*
(
s
)
+
M
d
i
Z
i
(
r
)
Z
i
(
s
)
)
=
δ
r
s
,
where δ
r
s is the Kronecker delta. The second set of orthogonality relation is expressed as
(28)
∑
i
=
1
N
‍
(
Y
i
*
(
r
)
Y
i
*
(
s
)
+
Z
i
(
r
)
Z
i
(
s
)
)
=
δ
r
s
.
The use of (22) and (25) with the aid of some algebraic manipulation yields the third and fourth set of orthogonality relation. This may be written as
(29)
∑
i
=
1
N
‍
(
E
I
∫
0
L
i
Y
i
′
′
(
r
)
Y
i
′
′
(
s
)
d
x
+
T
∫
0
L
i
Y
i
′
(
r
)
Y
i
′
(
s
)
d
x
)
+
∑
i
=
1
N
‍
k
i
(
Y
i
*
(
r
)
Y
i
*
(
s
)
-
Z
i
(
s
)
Y
i
*
(
r
)
-
Z
i
(
r
)
Y
i
*
(
s
)
+
Z
i
(
r
)
Z
i
(
s
)
)
=
δ
r
s
,
(30)
∑
i
=
1
N
‍
(
Y
i
*
(
r
)
Z
i
*
(
s
)
+
Z
i
(
r
)
Y
i
(
s
)
)
=
δ
r
s
.


## 5. Forced Vibration

Assume a harmonic forceF
(
x
,
t
) is arbitrarily applied along the span of the beam as depicted in Figure 1. The response of the loaded beam is now derived using the generalized orthogonality relations (27)–(30). Let the excitation force be applied at a location a
f from the left-hand end on the beam; the governing equations of motion are now
(31)
m
w
¨
i
+
E
I
w
i
′
′
′
′
-
T
w
i
′
′
=
F
(
x
,
t
)
δ
(
x
-
a
f
)
,
(32)
M
d
i
z
¨
i
+
k
i
(
z
i
-
w
i
(
L
i
)
)
+
c
d
i
(
z
˙
i
-
w
˙
i
(
L
i
)
)
=
0
.
Using the assumed mode method, the transverse displacement of the beam and the displacement of vibration absorber may be expressed as
(33)
w
i
=
∑
r
=
1
N
r
‍
q
r
(
t
)
Y
i
(
r
)
(
x
)
,
(34)
z
i
=
∑
r
=
1
N
r
‍
q
r
(
t
)
Z
i
(
r
)
,
where N
r is the number of retained modes, Y
i
(
r
)
(
x
) is the mode shape corresponding to the rth mode, Z
i
(
r
) is the displacement amplitude of the absorber, and q
r
(
t
) is the rth generalized coordinate. Substituting (33) and (34) into (31) and (32), respectively, yields
(35)
m
c
∑
r
=
1
N
r
‍
q
¨
r
Y
i
(
r
)
+
E
I
∑
r
=
1
N
r
‍
q
r
Y
i
′
′
′
′
(
r
)
-
T
∑
r
=
1
N
r
‍
q
r
Y
i
′
′
(
r
)
,
=
F
(
t
)
δ
(
x
-
a
f
)
M
d
i
∑
r
=
1
N
r
‍
q
¨
r
Y
m
i
(
r
)
+
k
i
∑
r
=
1
N
r
‍
q
r
(
Z
i
(
r
)
-
Y
c
i
(
r
)
(
L
i
)
)
+
c
d
i
∑
r
=
1
N
r
‍
q
˙
r
(
Z
i
(
r
)
-
Y
c
i
(
r
)
(
L
i
)
)
=
0
.
Multiplying (35) by Y
i
(
s
) and Z
i
(
s
), respectively, adding the two resulting equations, integrating over the entire length of the beam, and applying the continuity conditions, (12)–(15), with any classical boundary conditions (except for free ends) yield
(36)
∑
r
=
1
N
r
‍
q
¨
r
{
m
∫
0
L
i
Y
i
(
r
)
Y
i
(
s
)
d
x
+
M
c
i
Y
i
*
(
r
)
Y
i
*
(
s
)
+
M
d
i
Z
i
(
r
)
Z
i
(
s
)
}
+
c
d
i
∑
r
=
1
N
r
‍
q
˙
r
{
Y
i
*
(
r
)
Y
i
*
(
s
)
-
Z
i
(
s
)
Y
i
*
(
r
)
-
Z
i
(
r
)
Y
i
*
(
s
)
+
Z
i
(
r
)
Z
i
(
s
)
}
+
∑
r
=
1
N
r
‍
q
r
{
∫
0
L
i
(
E
I
Y
i
′
′
(
r
)
Y
i
′
′
(
s
)
+
T
Y
i
′
(
r
)
Y
i
′
(
s
)
)
d
x
+
k
i
(
Y
i
*
(
r
)
Y
i
*
(
s
)
-
Z
i
(
s
)
Y
i
*
(
r
)
-
Z
i
(
r
)
Y
i
*
(
s
)
+
Z
i
(
r
)
Z
i
(
s
)
)
∫
0
L
i
}
=
F
(
x
,
t
)
Y
i
(
s
)
(
a
f
)
.
Use of the orthogonality relations, (27)–(30), yields the following uncoupled differential equation:
(37)
[
M
r
r
]
{
q
¨
r
}
+
[
C
r
r
]
{
q
˙
r
}
+
[
K
r
r
]
{
q
r
}
=
{
F
r
}
,
where the matrices M
r
r, C
r, and K
r
r are expressed as
(38)
M
r
r
=
∑
i
=
1
N
‍
(
m
∫
0
L
i
Y
i
(
r
)
2
d
x
+
M
c
i
Y
i
*
(
r
)
2
+
M
d
i
Z
i
(
r
)
2
)
,
C
r
r
=
∑
i
=
1
N
‍
c
d
i
(
Y
i
*
(
r
)
-
Z
i
(
r
)
)
2
,
K
r
=
∑
i
=
1
N
‍
{
∫
0
L
i
(
E
I
Y
i
′
′
(
r
)
2
+
T
Y
i
′
(
r
)
2
)
d
x
+
k
i
(
Y
i
*
(
r
)
-
Z
i
(
r
)
)
2
}
,
F
r
=
F
(
x
,
t
)
Y
i
(
r
)
(
a
f
)
.
The amplitude of the vibration absorber can be readily expressed as
(39)
Z
i
(
r
)
=
κ
i
Y
i
*
(
r
)
,
κ
i
=
(
1
+
(
2
ς
i
r
i
)
2
(
1
-
r
i
2
)
2
+
(
2
ς
i
r
i
)
2
)
1
/
2
,
ς
i
=
c
d
i
2
M
d
i
ω
s
i
,
r
i
=
Ω
ω
s
i
,
ω
s
i
=
k
i
M
d
i
.


## 6. Numerical Simulation

The sets of parameters employed in the numerical examples are taken from [1, 16]; they are tabulated in Table 1. The validity of the free vibration of the present model is inferred from the results tabulated in Table 2. In the case of a cantilevered beam with three spring-mass-damper systems, the first five natural frequencies of the present model are compared to those of [16]. The results show very good agreement with a maximum error of 0.28%.Table 1
Material properties and parameters.
Parameter
Reference [1]
Reference [16]
Units
EI
63476.1
36.4583
Nm2
L
1.0
1.0
m
m
15.385
0.675
kg/m
m
1
3.0775
0.1
kg
m
2
4.614
0.1
kg
m
3
7.69
0.1
kg
m
4
9.997
—
kg
m
5
15.38
—
kg
k
1
190428
0.1
N/m
k
2
222166
0.1
N/m
k
3
285642.45
0.1
N/m
k
4
317380
—
N/m
k
5
380856
—
N/m
c
d
1, c
d
2, c
d
3
—
0.1
Ns/mTable 2
Natural frequency validation (rad/s).
Mode
Cantilevered
Pinned-pinned
3 spring-mass-damper
3 spring-mass
5 spring-mass
Reference [16]
Present
Reference [1]
Present
Reference [1]
Present
1
-
0.367
±
25.8769
i
-
0.367
±
25.8544
i
152.7339
152.7339
150.9571
150.9571
2
-
0.274
±
161.7655
i
-
0.275
±
161.9399
i
185.0949
185.0949
169.4728
169.4728
3
-
0.031
±
452.1367
i
-
0.031
±
452.4320
i
247.8313
247.8313
187.9146
187.9146
4
-
0.171
±
887.2157
i
-
0.174
±
887.3469
i
677.5959
677.5959
217.1278
217.1278
5
-
0.081
±
1467.5170
i
-
0.082
±
1468.8356
i
2548.6572
2548.6572
247.9867
247.9867With respect to a pinned-pinned beam with attached three and five spring-mass systems, the validation is done via [1] and shows excellent agreement with four decimal places. The corresponding mode shapes for a beam carrying the three in-span spring-mass systems are depicted in Figure 2. This figure is identical to Figure 3 from [1].Figure 2
Mode shapes of a simply supported beam carrying three in-span spring-mass systems using the same parameters from [1].For the forced vibration simulations, an excitation force ofF
(
t
)
=
10
sin
⁡
(
Ω
t
) is used with zero initial conditions. The force is first applied at the free end of a cantilevered beam with three identical spring-mass-damper systems located at 0.1, 0.5, and 0.9 m from the clamp end. The time history of the vertical displacement of the free end is presented to serve for comparison between the proposed approach (i.e., the combination of the generalized orthogonality conditions with the exact mode shapes of the loaded beam), the finite element method (FEM), and the approach used in [16]. The results are depicted in Figures 3 and 4 for an excitation frequency of 5 and 10 rad/s, respectively. All three methods are in agreement, but the FEM results agree better with the present approach than that of [16].Figure 3
Time history of the free-end vertical displacement for a cantilevered beam carrying three mass-spring-damper systems for a forcing frequencyΩ
=
5 rad/s, k
=
0.1 N/m, T
=
0, M
d
=
0.1 kg, M
c
=
0 kg, c
d
=
0.1 Ns/m, and L
1
;
2
;
3
=
0.1, 0.5, 0.9 m.Figure 4
Time history of the free-end vertical displacement for a cantilevered beam carrying three mass-spring-damper systems for a forcing frequencyΩ
=
10 rad/s, k
=
0.1 N/m, T
=
0, M
d
=
0.1 kg, M
c
=
0 kg, c
d
=
0.1 Ns/m, and L
1
;
2
;
3
=
0.1, 0.5, 0.9 m.The frequency response curves for a simply supported beam with one in-span mass-spring-damper-mass system attached at 0.1 m from one end using three approaches are depicted in Figure5. The first approach, which is the most accurate [25], employs the generalized orthogonality conditions along with the mode shapes of the loaded beam. The classical orthogonality relations and the mode shapes of the loaded beam are employed in the second approach, such as in [17–20]. The third method also uses the classical orthogonality conditions, but the mode shapes are those of the bare beam (see [14–16]). All three approaches yield very different midspan vertical displacements. The first approach is the exact solution and hence the error associated with the difference is defined with respect to the first approach. The frequency response amplitude errors are depicted in Figure 6 where the error can be as high as 104.Figure 5
Frequency response of a simply supported beam carrying one in-span mass-spring-damper-mass system fork
=
10 N/m, M
d
=
0.1 kg, M
c
=
0.05 kg, T
=
360 N, E
I
=
36.458 Nm2, c
d
=
0.1 Ns/m, and L
d
=
0.1 m.Figure 6
Frequency response amplitude error for a simply supported beam carrying one in-span mass-spring-damper-mass system fork
=
10 N/m, M
d
=
0.1 kg, M
c
=
0.05 kg, T
=
360 N, E
I
=
36.458 Nm2, c
d
=
0.1 Ns/m, and L
d
=
0.1 m.Figure7 is an identical plot to Figure 5 with the pinned-pinned configuration replaced by a guided-guided one. The corresponding frequency response amplitude error is shown in Figure 8. The maximum amplitude error for the guided-guided beam is higher than that for the simply supported beam. A similar observation with regard to the discrepancy associated with the use of the classical orthogonality was reported in [25] for a simply supported beam carrying a heavy mass. An error as high as 105 was reported.Figure 7
Frequency response of a guided-guided beam carrying one in-span mass-spring-damper-mass system fork
=
10 N/m, M
d
=
0.1 kg, M
c
=
0.05 kg, T
=
360 N, E
I
=
36.458 Nm2, c
d
=
0.1 Ns/m, and L
d
=
0.1 m.Figure 8
Frequency response amplitude error of a guided-guided beam carrying one in-span mass-spring-damper-mass system fork
=
10 N/m, M
d
=
0.1 kg, M
c
=
0.05 kg, T
=
360 N, E
I
=
36.458 Nm2, c
d
=
0.1 Ns/m, and L
d
=
0.1 m.The three approaches are also examined using the parameters taken from [1] for a simply supported beam with one in-span mass-spring-damper system and no tension. Figures 9 and 10 show the frequency response of the midspan displacement for a damping coefficient c
d
i
=
0.1 Ns/m and c
d
i
=
100 Ns/m, respectively. It can be observed that all three methods yield very similar results and the plots are barely distinguishable. The amplitude errors associated with the results presented in Figures 9 and 10 are depicted in Figures 11 and 12, respectively. But for excitation frequencies closer to 600 rad/s the results indicate very minimal error. The discrepancy is more pronounced in the case of the classical orthogonality with the combination of the bare beam mode shapes.Figure 9
Frequency response of a simply supported beam carrying one in-span mass-spring-damper-mass system fork
=
190428 N/m, M
d
=
3.0775 kg, M
c
=
0 kg, T
=
0 N, E
I
=
63476.1 Nm2, c
d
=
0.1 Ns/m, and L
d
=
0.75 m.Figure 10
Frequency response of a simply supported beam carrying one in-span mass-spring-damper-mass system fork
=
190428 N/m, M
d
=
3.0775 kg, M
c
=
0 kg, T
=
0 N, E
I
=
63476.1 Nm2, c
d
=
100 Ns/m, and L
d
=
0.75 m.Figure 11
Frequency response amplitude error of a simply supported beam carrying one in-span mass-spring-damper-mass system fork
=
190428 N/m, M
d
=
3.0775 kg, M
c
=
0 kg, T
=
0 N, E
I
=
63476.1 Nm2, c
d
=
0.1 Ns/m, and L
d
=
0.75 m.Figure 12
Frequency response amplitude error of a simply supported beam carrying one in-span mass-spring-damper-mass system fork
=
190428 N/m, M
d
=
3.0775 kg, M
c
=
0 kg, T
=
0 N, E
I
=
63476.1 Nm2, c
d
=
100 Ns/m, and L
d
=
0.75 m.


## 7. Conclusions

A common approach to studying the vibrational response of an elastically loaded beam involves the use of the classical orthogonality relations along with the mode shapes of the bare beam. It has been shown, however, that the mode shapes of the bare beam can be quite different from those of the loaded beam. In the present paper, the exact natural frequencies and mode shapes were presented for an axially loaded beam carrying several vibration absorbers. Explicit expressions were presented for the generalized orthogonality condition. The obtained generalized orthogonality relation was employed along with the assumed modes method to study the forced vibrational response. The numerical simulations indicated that using the common approach could produce erroneous results. The combination of the exact mode shapes of the loaded beam and the corresponding generalized orthogonality relation is necessary for accurate dynamic modeling of a beam carrying elastically mounted masses with dampers.



---

*Source: 485630-2014-02-23.xml*