FREE VIBRATION ANALYSIS OF EULER-BERNOULLI BEAMS UNDER

NON-CLASSICAL BOUNDARY CONDITIONS

João Fernandes da Silva, jfernandes.mecanica@gmail.com1
Lucas Allende Dias do Nascimento, allendezim@gmail.com1
Simone dos Santos Hoefel, simone.santos@ufpi.edu.br1

1
Universidade Federal do Piauí, Campus Ministro Petrônio Portella, s/n - Ininga, Teresina-PI, 64049-550

Abstract: The flexible beams carrying attachments and non-classical boundary conditions often appear in engineering
structures, modal analysis of those structures is important and necessary in structural design. The analysis of vibrating
beams with ends elastically restrained against rotation and translation or with ends carrying concentrated masses or
rotational inertias is of great interest in a variety of practical cases. In this analysis, the governing differential equations
of the beam, which is a partial differential equation with variable coefficients, and that of the mass-spring system, which
is an ordinary differential equation, are found. The exact solution of the problem is then obtained using classical and
non-classical boundary conditions and the eigenvalues and eigenfunctions are found. Finally, some cases with available
results in the literature are presented and analyzed.
Keywords: Euler-Bernoulli beam, Non-classical boundary conditions, Natural frequencies

1. INTRODUCTION

The interest in non-classical mechanical boundary conditions is originated in the study of stability of elastic structures,
and a convenient reference for the various mechanical loadings considered until 1963 is the book by Bolotin (1963).
The dynamic analysis of beams with ends elastically restrained against rotation and translation or with ends carrying
concentrated masses or rotational inertia is a classical structural problem, which nowadays is becoming more and more
important, even in mechanical engineering and in aeronautic engineering. Numerous authors have approached the analysis
assuming that the beam is sufficiently slender to be considered as an Euler-Bernoulli beam, and trying to analytically solve
the resulting fourth-order differential equation with variable coefficients. In exact method, difficulties arises in solving
roots of the characteristic equation, except for very simple boundary conditions, that one has to go for numerical solution
and in determining the normal modes of the system.
Grossi and Arenas (1996) employed both the classical Rayleigh-Ritz method and the optimized Rayleigh-Schmidt
method to find the frequencies with varying width and height. A lot of numerical results were given, for various non-
classical boundary conditions. Craver Jr. and Jampala (1993) examined the free vibration frequencies of a cantilever
beam with variable cross-section and constraining springs. De Rosa and Auciello (1996) gave the exact free frequencies
of a beam with linearly varying cross-section, in the presence of generic non-classical boundary conditions, so that all the
usual boundary conditions can be treated as particular cases. Chun (1972) studied the free vibration of a Bernoulli-Euler
beam hinged at one end by a rotational spring with constant spring stiffness and with the other end free. Wang and Lin
(1996) utilized the Fourier series to investigate the dynamic analysis of beams having arbitrary boundary conditions. The
Rayleigh-Ritz approach was used by Zhou and Cheung (2000) to find the first free vibration frequencies of three different
tapered beams with various boundary conditions and truncation factors. Maurizi et al. (1976) studied the problem of
free vibration of a uniform beam hinged at one end by a rotational spring and subjected to the restraining action of a
translational spring at the other end using exact expression of trigonometric and hyperbolic functions. Maurizi et al.
(1990) studied ends elastically restrained against rotation and translation in Timoshenko beams. Frequency equation of
the cantilever beams carrying springs and masses by use the normal summation technique, was presented by Gürgöze
(1984). Gürgöze (1998) also presented two formulations for the frequency equation of a clamped-free Euler-Bernoulli
beam to which several spring-mass systems are attached in span. Recently, Liu and Gurram (2009) used his variational
iteration method to calculate the natural frequencies and mode shapes of an Euler-Bernoulli beam under various supporting
conditions.
 I X C o n g r e s s o N a c i o n a l d e E n g e n h a r i a M e c â n i c a, 2 1 a 2 5 d e a g o s t o d e 2 0 1 6, F o r t a l e z a - C e a r á

This article focuses on the free vibration analysis of Euler-Bernoulli beams under non-classical boundary conditions.
The governing differential equations of the beam are presented and the exact solution of the problem is then obtained
using the pertinent boundary conditions. The eigenvalues and eigenfunctions of the problem which are frequencies and
mode shapes of the system are derived for various properties of the system such as stiffness of springs and attached mass.
Some numerical examples with available results in the literature are analyzed and its results discussed.

2. EULER-BERNOULLI BEAM THEORY

The partial differential equation of motion for free vibration of a Euler-Bernoulli beam is given by (Euler, 1773):

∂ 4 v(x, t) ∂ 2 v(x, t)
EI + ρA = 0, (1)
∂x4 ∂t2
where A is the cross-sectional area, I is the moment of inertia of the cross-sectional area, E is the modulus of elasticity,
ρ is the mass per unit of volume and v(x, t) is the transverse deflection at the axial location x and at the time t. Since the
equation of the motion involves a second order derivative with respect to time and a fourth order derivative with respect to
position, two initial equations and four boundary conditions are needed for finding a unique solution for v(x, t). Equation
of motion can be solved by using the methods of separation of variables and eigenfunctions expansion. The transverse
deflection v(x, t) are separated into two functions as follows:

v(x, t) = V (x) T (t), (2)

where V (x) is know as the normal mode or characteristic function of the beam. Therefore, the Eq. (1) is rewritten as (Han
et al., 1999):

1 ∂ 4 V (x) 1 ∂ 2 T (t)
c2 4
=− = ω2 , (3)
V (x) ∂x T (t) ∂t2
p
knowing that c = EI/ρA and ω is the natural frequency. Note that the left side of the Eq. (3) is relation to position and
the right side to time and both can be splitted into two equations (Avcar, 2014):
d4 W (x)
− β 4 W (x) = 0, (4)
dx4

d2 T (t)
+ ω 2 T (t) = 0, (5)
dt2
where β is the eigenvalue of the normal mode V (x) and
r
ω
β= . (6)
c

Solutions of T (t) and V (x) in Eq. (4) and Eq. (5), respectively, are given by:
T (t) = D cos (ωt − φ), (7)
V (x) = C1 sin(βx) + C2 cos(βx) + C3 sinh(βx) + C4 cosh(βx). (8)

C1 , C2 , C3 and C4 are constants which can be found by using the appropriate boundary conditions, D is a constant
which can be obtained by using initial conditions and φ is the phase angle.

3. BOUNDARY CONDITIONS AND FREQUENCY EQUATIONS

Each one of the boundary conditions provides two equations necessary to determine the unknown constants C1 to C4
and the values of β in Eq. (8). They are usually classified in two groups: the classical and the non-classical boundary
conditions. The classical boundary conditions are the most common in the vibration analysis in beams and are generally
classified in free, clamped and supported ends and their equations are showed in Tab. 1. Non-classical boundary conditions
are less common, but their importance are as the same as the classical conditions. These boundary conditions consists
in points connected to masses, springs, rotational inertia or, even, dampers. All those conditions are presented in Tab. 2,
where a = −1 and b = 1 for the left end and a = 1 and b = −1 for the right end. km is the stiffness of the linear spring,
kr is the stiffness of the torsional spring, cm is the stiffness of the damper, cr is the stiffness of the torsional damper, mc
is the value of the concentrated mass and the Ic is the value of the rotational inertia of the concentrated mass.
 I X C o n g r e s s o N a c i o n a l d e E n g e n h a r i a M e c â n i c a, 2 1 a 2 5 d e a g o s t o d e 2 0 1 6, F o r t a l e z a - C e a r á

Table 1: Classical Boundary Conditions.

Boundary Condition Equations

∂ 2 v(x, t) ∂ 3 v(x, t)
Free end EI =0 and EI =0
∂x2 ∂x3

∂ 2 v(x, t)
Supported end v(x, t) = 0 and EI =0
∂x2
∂v(x, t)
Clamped end v(x, t) = 0 and EI =0
∂x

Table 2: Non-classical Boundary Conditions.

Boundary Condition Equations

∂ 2 v(x, t) ∂ 3 v(x, t)
Linear Spring EI =0 and EI = akm v(x, t)
∂x2 ∂x3

∂ 3 v(x, t) ∂ 2 v(x, t) ∂v(x, t)

Torsional Spring EI =0 and EI = bkr
∂x3 ∂x2 ∂x

∂ 2 v(x, t) ∂ 3 v(x, t) ∂ 2 v(x, t)

Concentrated Mass EI =0 and EI = am c
∂x2 ∂x3 ∂t2

∂ 3 v(x, t) ∂ 2 v(x, t) ∂ 3 v(x, t)

Rotational Inertia EI =0 and EI = bI c
∂x3 ∂x2 ∂x ∂t2

∂ 2 v(x, t) ∂ 3 v(x, t) ∂v(x, t)

Damper EI =0 and EI = acm
∂x2 ∂x3 ∂t

∂ 3 v(x, t) ∂ 2 v(x, t) ∂ 2 v(x, t)

Torsional Damper EI =0 and EI = bcr
∂x3 ∂x2 ∂x ∂t

The application of appropriate boundary conditions presented in Tab. 2 into the Eq. (8) and its derivatives yields for
each type of beam a set of four homogeneous linear algebraic equations in four constants C1 to C4 . A non-trivial solution
exists only if the determinant of the coefficients vanishes. The solving of the determinant yields the frequency equation or
the characteristic equation for each type of beam. Table 3 show frequency equations for some beams under non-classical
boundary conditions presented in Tab. 2, where:

Q1 = [ cos( βn l ) cosh( βn l ) − 1 ],

Q2 = sin( βn l ) sinh( βn l ),

Q3 = sin( βn l ) cosh( βn l),

Q4 = cos( βn l ) sinh( βn l ) and

Q5 = cos( βn l ) cosh( βn l ).
 I X C o n g r e s s o N a c i o n a l d e E n g e n h a r i a M e c â n i c a, 2 1 a 2 5 d e a g o s t o d e 2 0 1 6, F o r t a l e z a - C e a r á

Table 3: Frequency equations for some beams under non-classical boundary conditions.

Boundary Condition Frequency Equations

Linear Spring - Linear Spring (EI)2 β 6 Q1 − 2km

2
Q2 + 2EIkm β 3 [Q3 − Q4 ] = 0.

Linear Spring - Concentrated Mass ρAEIβ 3 Q1 + EImc β 4 [Q4 − Q3 ] + km ρA[Q3 − Q4 ] + 2km βQ2 = 0.

Linear Spring - Torsional Spring EIβkm [Q3 − Q4 ] + EIβ 3 kr [Q4 + Q3 ] − 2km kr Q5 + (EI)2 β 4 Q1 = 0.

Torsional Spring - Torsional Spring 2kr2 Q2 + E 2 I 2 β 2 Q1 + 2βkr EI[Q4 + Q3 ] = 0.

These frequency equations are highly transcendental and need to be solved numerically. In this paper was used the
false position root finding method.

4. NUMERICAL EXAMPLES

This section presents numerical examples for Euler-Bernoulli beams subjected to three kind of boundary conditions.
The same parameters values are considered for all examples and are given by EI = 8000 N.m2 , ρA = 3000 kg/m,
L = 1 m, mc = 100 kg, km = 1000 N/m, kr = 1000 N.m/rad and the frequency equation are showed in Tab. 3 for
each case analyzed.

4.1 Beam connected to linear springs at both ends

Consider a beam connected to linear springs at both ends shown in Fig. 1, the natural frequencies values and the
eigenvalues obtained from the frequency equation are presented in Tab. 4.

km km
Figure 1: Beam connected to linear springs at both ends.

Table 4: Natural frequencies values ω and eingenvalues βn of a Euler-Bernoulli beam connected to linear springs
at both ends.

n=1 n=2 n=3 n=4 n=5 n=6 n=7

Eigenvalues βn 4, 7324 7, 8537 10, 9958 14, 1373 17, 2788 20, 4204 23, 5620
Natural frequencies ωn 36, 5719 100, 7245 197, 4412 326, 3732 487, 5419 680, 9451 906, 5825

The normal mode V (x) and the slope θ(x), which is the derivative of V (x), are as follow:
V (x) = σ1 [cos(βx) + cosh(βx)] + σ2 sin(βx) + sinh(βx), (9)
θ(x) = β{σ1 [−sin(βx) + sinh(βx)] + σ2 cos(βx) + cosh(βx)}, (10)
where
EIβ 3 [sinh(βn l) − sin(βn l)]
σ1 =
−EIβ 3 cosh(β 3
n l) + 2km sin(βn l) + EIβ cos(βn l)

and
2km sinh(βn l) − EIβ 3 cos(βn l) + EIβ 3 cosh(βn l)
σ2 = .
−EIβ 3 cosh(βn l) + 2km sin(βn l) + EIβ 3 cos(βn l)
Figure 2 showed mode shape of the first four eigenvalues and Fig. 3 showed the plots of slope. We can observed that
the major deflection always occurs at the ends of the beam. Note that the modal shapes and graphs of slope of this beam
are very similar to those of a free-free beam, as showed in Fig. 4 and Fig. 5. Also, if the stiffness km is equal to zero, the
resulting frequency equation is one for a free-free beam as follows
cos(βn l) cosh(βn l) − 1 = 0. (11)
 I X C o n g r e s s o N a c i o n a l d e E n g e n h a r i a M e c â n i c a, 2 1 a 2 5 d e a g o s t o d e 2 0 1 6, F o r t a l e z a - C e a r á

2
1st Modal Shape 2nd Modal Shape 3rd Modal Shape 4th Modal Shape
1.5

0.5

0
V(x)

−0.5

−1

−1.5

−2

−2.5
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
x
Figure 2: Modal Shapes for a Euler-Bernoulli beam connected to linear springs at both ends.

30
β1 β2 β3 β4

20

10
θ(x)

−10

−20

−30
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
x
Figure 3: Slope of Euler-Bernoulli beam connected to linear springs at both ends associate to the first four eigen-
values βn .

2
1st Modal Shape 2nd Modal Shape 3rd Modal Shape 4th Modal Shape
1.5

0.5

0
V(x)

−0.5

−1

−1.5

−2

−2.5
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
x
Figure 4: Modal Shapes for a free-free Euler-Bernoulli beam.

Furthermore, the natural frequency values found are very closer from the values for a free-free beam, which are
presented in Tab. 5. For a free-free beam, the normal mode V(x) and the slope θ(x) are:
V (x) = σ[cos(βx) + cosh(βx)] + sin(βx) + sinh(βx), (12)
θ(x) = β{σ[−sin(βx) + sinh(βx)] + cos(βx) + cosh(βx)}, (13)
where
sin(βn l) − sinh(βn l)
σ= .
cosh(βn l) − cos(βn l)

4.2 Beam connected to a linear spring at left end and to a concentrated mass at right end

Now, a beam connected to a linear spring at left end and to a concentrated mass at right end shown in Fig. 6 is analyzed.
The normal mode V (x) and the slope θ(x) are:
V (x) = σ1 [cos(βx) + cosh(βx)] + σ2 sin(βx) + sinh(βx), (14)
 I X C o n g r e s s o N a c i o n a l d e E n g e n h a r i a M e c â n i c a, 2 1 a 2 5 d e a g o s t o d e 2 0 1 6, F o r t a l e z a - C e a r á

30
β1 β2 β3 β4

20

10
θ(x)

−10

−20

−30
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
x
Figure 5: Slope of free-free Euler-Bernoulli beam associate to the first four eigenvalues βn .

Table 5: Natural frequencies values ω and eigenvalues βn of a free-free Euler-Bernoulli beam.

n=1 n=2 n=3 n=4 n=5 n=6 n=7

Eigenvalues βn 4, 7300 7, 8532 10, 9956 14, 1372 17, 2788 20, 4204 23, 5619
Natural frequencies ωn 36, 5354 100, 7113 197, 4344 326, 3691 487, 5391 680, 9431 906, 5811

θ(x) = β{σ1 [−sin(βx) + sinh(βx)] + σ2 cos(βx) + cosh(βx)}, (15)

where
EIβ 3 [sinh(βn l) − sin(βn l)]
σ1 =
−EIβ 3 cosh(βn l) + 2km sin(βn l) + EIβ 3 cos(βn l)
and
2km sinh(βn l) − EIβ 3 cos(βn l) + EIβ 3 cosh(βn l)
σ2 = .
−EIβ 3 cosh(βn l) + 2km sin(βn l) + EIβ 3 cos(βn l)

L
mc
km
Figure 6: Beam connected to a linear spring at left both end and to a concentrated mass at right end.

The natural frequencies values ω and the eigenvalues β are presented in Tab. 6. As the previous beam, the modal
shapes in Fig. 7 and the slopes in Fig. 8 the are similar to those of a free-free beam, except for the deflection at left end
be greater than at the right end. Besides, considering mc and km equals to zero, the frequency equation becomes equal to
frequency equation of free-free beam.

Table 6: Natural frequencies values ω and eigenvalues βn of a Euler-Bernoulli beam connected to spring at left
end and to a concentrated mass at right end.

n=1 n=2 n=3 n=4 n=5 n=6 n=7

Eigenvalues βn 4, 5988 7, 6530 10, 7380 13, 8315 16, 9325 20, 0395 23, 1512
Natural frequencies ωn 34, 5368 95, 6418 188, 2911 312, 4104 468, 1967 655, 7769 875, 2458

4.3 Beam connected to a linear spring at left end and to a torsional spring at right end

In this example, a beam connected to a linear spring at left end and to a torsional spring at right end is considered,
as showed in the Fig. 9. The natural frequencies values and the eigenvalues are presented in Tab. 7 and the normal mode
V (x) and the slope θ(x) are:

V (x) = σ1 [−cos(βx) − cosh(βx)] + σ2 sin(βx) + sinh(βx), (16)

θ(x) = β{σ1 [sin(βx) − sinh(βx)] + σ2 cos(βx) + cosh(βx)}, (17)

 I X C o n g r e s s o N a c i o n a l d e E n g e n h a r i a M e c â n i c a, 2 1 a 2 5 d e a g o s t o d e 2 0 1 6, F o r t a l e z a - C e a r á

2
1st Modal Shape 2nd Modal Shape 3rd Modal Shape 4th Modal Shape
1.5

0.5

0
V(x)

−0.5

−1

−1.5

−2

−2.5
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
x
Figure 7: Modal Shapes for a Euler-Bernoulli beam connected to a linear spring at left end and to a concentrated
mass at right end.

30
β1 β2 β3 β4

20

10
θ(x)

−10

−20

−30
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
x
Figure 8: Slope of Euler-Bernoulli beam connected to a linear spring at left end and to a concentrated mass at
right end associate to the first four eigenvalues βn .

where
EIβ 3 [cos(βn l) − cosh(βn l)]
σ1 =
2km cos(βn l) − EIβ 3 sinh(βn l) − EIβ 3 sin(βn l)
and
2km cosh(βn l) − EIβ 3 sinh(βn l) − EIβ 3 sin(βn l)
σ2 = .
2km cos(βn l) − EIβ 3 sinh(βn l) − EIβ 3 sin(βn l)

km kr

Figure 9: Beam connected to a linear spring at left both end and to a torsional spring at right end.

Modal shapes presented in Fig. 10 and the plots of slopes in Fig. 11 are still very similar to those of a free-free beam
and to those of the previous beams. Furthermore, if km and kr are zero, the result is similar to the frequency equation for
a free-free beam.

Table 7: Natural frequencies values ω and eingenvalues βn of a Euler-Bernoulli beam connected to a linear spring
at left end and to a torsional spring at right end.

n=1 n=2 n=3 n=4 n=5 n=6 n=7

Eigenvalues βn 4, 7559 7, 8691 11, 0069 14, 1460 17, 2860 20, 4264 23, 5672
Natural frequencies ωn 36, 9365 101, 1199 197, 8412 326, 7757 487, 9457 681, 3498 906, 9878
 I X C o n g r e s s o N a c i o n a l d e E n g e n h a r i a M e c â n i c a, 2 1 a 2 5 d e a g o s t o d e 2 0 1 6, F o r t a l e z a - C e a r á

2
1st Modal Shape 2nd Modal Shape 3rd Modal Shape 4th Modal Shape
1.5

0.5

0
V(x)

−0.5

−1

−1.5

−2

−2.5
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
x
Figure 10: Modal Shapes for a Euler-Bernoulli beam connected to a linear spring at left end and to a torsional
spring at right end.

30
β1 β2 β3 β4

20

10
θ(x)

−10

−20

−30
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
x
Figure 11: Slope of Euler-Bernoulli beam connected to a linear spring at left end and to a torsional spring at right
end associate to the first four eigenvalues βn .

4.4 Beam connected to torsional springs at both ends

Finally, consider a beam connected to a torsional spring at left end and to a torsional spring at right end as showed in
Fig. 12 and with normal mode V (x) and the slope θ(x) given by:

V (x) = σ1 [sin(βx) + sinh(βx)] + σ2 cos(βx) + cosh(βx), (18)

θ(x) = β{σ1 [cos(βx) + cosh(βx)] − σ2 sin(βx) + senh(βx)}, (19)

where
EIβ[sinh(βn l) + sin(βn l)]
σ1 =
EIβ[cos(βn l) − cosh(βn l)] + 2βkr sin(βn l)
and
EIβ[cos(βn l) − cosh(βn l)] − 2βkr sinh(βn l)
σ2 = .
EIβ[cos(βn l) − cosh(βn l)] + 2βkr sin(βn l)

kr kr

Figure 12: Beam connected to torsional springs at both ends.

According Tab. 8, the highest values of the natural frequencies values ω and the eigenvalues β were found for the beam
connected to torsional springs. For this beam, the modal shapes showed in the Fig. 13 and the graphs of slope in Fig. 14
are the inverse for those of each one of the preceding beams and for a free-free beam, but they still have a similarity.
Besides, the frequency equation for a free-free beam is obtained if the value of kr is zero.
 I X C o n g r e s s o N a c i o n a l d e E n g e n h a r i a M e c â n i c a, 2 1 a 2 5 d e a g o s t o d e 2 0 1 6, F o r t a l e z a - C e a r á

Table 8: Natural frequencies values ω and eingenvalues βn of a Euler-Bernoulli beam connected to torsional
springs at both ends.

n=1 n=2 n=3 n=4 n=5 n=6 n=7

Eigenvalues βn 4, 7793 7, 8845 11, 0180 14, 15470 17, 2931 20, 4325 23, 5725
Natural frequencies ωn 37, 3004 101, 5146 198, 2408 327, 1779 488, 3494 681, 7544 907, 3931

2.5
1st Modal Shape 2nd Modal Shape 3rd Modal Shape 4th Modal Shape
2

1.5

0.5
V(x)

−0.5

−1

−1.5

−2

−2.5
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
x
Figure 13: Modal Shapes for a Euler-Bernoulli beam connected to torsional springs at both ends.

30
β1 β2 β3 β4
20

10

0
θ(x)

−10

−20

−30

−40
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
x
Figure 14: Slope of Euler-Bernoulli beam connected to torsional springs at both ends associate to the first four
eigenvalues βn .

5. CONCLUSION

In this paper, free vibration of beams considering different boundary conditions are analyzed. The natural frequencies
and the eigenvalues for Euler-Bernoulli beams are obtained and their modal shapes are showed in graphs. The eigenvalues
and the eigenfunctions, for each case analyzed, are very similar to the available results in the literature, like the studies
conducted by Morelatto (2000).
It was noted that the major deflection and slope of the beams in study always happens at their ends. That occurs
because the spring-mass-beam system presents an initial deflection in order to present an static equilibrium. Since that
does not have any excitation after the static equilibrium, the springs or masses will not deflect anymore. Meanwhile, in
free vibration, the entire beam is vibrating trying to restore to the position of equilibrium, always finding a maximum
deflection and slopes at the ends of the beams. For this reason, the behavior of these beams are similar to a free-free beam,
even their frequency equation can become that of a free-free beam, if the values for the stiffness of the springs and for the
masses are equals to zero.

6. REFERENCES

Avcar, M., 2014. “Free vibration analysis of beams considering different geometric characteristics and boundary condi-
tions”. Journal of mechanics and applications, Vol. 4(3), pp. 94–100.
Bolotin, V.V., 1963. Nonconservative Problems of the Theory of Elastic Stability. Macmillan Co., New York.
Chun, K.R., 1972. “Free vibration of a beam with one end spring-hinged and the other free”. Journal of Applied
Mechanics, Vol. 39, pp. 1154–1155.
Craver Jr., L. and Jampala, P., 1993. “Transverse vibrations of a linearly tapered cantilever beam with constraining
 I X C o n g r e s s o N a c i o n a l d e E n g e n h a r i a M e c â n i c a, 2 1 a 2 5 d e a g o s t o d e 2 0 1 6, F o r t a l e z a - C e a r á

springs”. Journal of Sound and Vibration, Vol. 166, pp. 521–529.

De Rosa, M.A. and Auciello, N.M., 1996. “Free vibrations of tapered beams with exible ends”. Computers Structures,
Vol. 60, pp. 197–202.
Euler, L., 1773. “De motu vibratorio laminarum elasticarum, ubi plures novae vibrationum species hactenus non pertrac-
tatae evolvuntur.” Novi Commentarii academiae scientiarum Petropolitanae, Vol. 17, pp. 449–487.
Grossi, R.O. and Arenas, B.d.V., 1996. “A variational approach to the vibration of tapered beams with elastically restrained
ends”. Journal of Sound and Vibration, Vol. 195, pp. 507–511.
Gürgöze, M., 1984. “A note on the vibrations of restrained beams and rods with point masses”. Journal of Sound and
Vibration, Vol. 96(4), pp. 461–468.
Gürgöze, M., 1998. “On the alternative formulations of the frequency equation of a bernoulli-euler beam to which several
spring-mass systems are attached in-span”. Journal of Sound and Vibration, Vol. 217, pp. 585–595.
Han, S.M., Benaroya, H. and Wei, T., 1999. “Dynamics of transversely vibrating beams using four engineering theories”.
Journal of Sound and Vibration, Vol. 225(5), pp. 935–988.
Liu, Y. and Gurram, C.S., 2009. “The use of he’s variational iteration method for obtaining the free vibration of an
euler–bernoulli beam”. Mathematical and Computer Modelling, Vol. 50, pp. 1545–1552.
Maurizi, M.J., Rosi, R.E. and Reyes, J.A., 1976. “Vibration frequencies for a beam with one end spring-hinged and
subjected to a translational restraint at the other end”. Journal of Sound and Vibration, Vol. 48, pp. 565–568.
Maurizi, M., Rossi, R. and Belles, P., 1990. “Free vibrations of uniform timoshenko beams with ends elastically restrained
against rotation and translation”. Journal of Sound and Vibration, Vol. 141(2), pp. 359–362.
Morelatto, T., 2000. Cálculo de Modos Vibratórios no Modelo Estrutural de Euler-Bernoulli com Condições de Contorno
Não-Clássicas. Master’s thesis, Instituto de Matemática, Universidade Federal do Rio Grande do Sul, Porto Alegre.
2000.
Wang, J.T.S. and Lin, C.C., 1996. “Dynamic analysis of generally supported beams using fourier series”. Journal of
Sound and Vibration, Vol. 196(3), pp. 285–293.
Zhou, D. and Cheung, Y.K., 2000. “The free vibrations of a type of tapered beams”. Computer methods in applied
mechanics and engineering, Vol. 188, pp. 203–216.

7. RESPONSIBILITY NOTICE

The authors are the only responsible for the printed material included in this paper.