Abstract
This paper deals with a viscoelastic beam obeying a fractional differentiation constitutive law. The governing equation is derived from the viscoelastic material model. The equation of motion is solved by using the method of multiple scales. Additionally, principal parametric resonances are investigated in detail. The stability boundaries are also analytically determined from the solvability condition. It is concluded that the order and the coefficient of the fractional derivative have significant effect on the natural frequency and the amplitude of vibrations.
Keywords:
perturbation method; fractional derivative; method of multiple scales; linear vibrations1 Introduction
Many researchers have demonstrated the potential of viscoelastic materials to improve the dynamics of fractionally damped structures. Fractional derivatives are practically used in the field of engineering for describing viscoelastic features in structural dynamics [1]. Namely, linear or nonlinear vibrations of axially moving beams have been studied extensively by many researchers [2]. Fractional derivatives are used in the simplest viscoelastic models for some standard linear solid. It can be seen that the vibrations of the continuum are modeled in the form of a partial differential equation system [3]. These damping models involve ordinary integer differential operators that are relatively easy to manipulate [4]. On the other hand, fractional derivatives have more advantages in comparison with classical integerorder models [5].
The partial differential equations of fractional order are increasingly used to model problems in the continuum and other areas of application. The field of fractional calculus is of importance in various disciplines such as science, engineering, and pure and applied mathematics [6]. The numerical solution for the time fractional partial differential equations subject to the initialboundary value is introduced by Podlubny [5]. The finite difference method for a fractional partial differential equation is presented by Zhang [7]. Galucio et al. developed a finite element formulation of the fractional derivative viscoelastic model [4]. Chen et al. studied the transient responses of an axially accelerating viscoelastic string constituted by the fractional differentiation law [8]. Applications of the method of multiple scales to partial differential systems arising in nonlinear vibrations of continuous systems were considered by Boyacı and Pakdemirli [9]. The method of multiple scales is one of the most common perturbation methods used to investigate approximate analytical solutions of dynamical systems. The dynamic response of the continuum is analyzed by using this method.
In this paper, longitudinal vibrations of the beam with external harmonic force are studied. The model developed is used to show the applicability of the fractional damped model and to find an approximate solution of the problem. The RiemannLiouville fractional operator is emphasized among several definitions of a fractional operator [10,11]. On the other hand, the approximate solution of the beam modeled by a fractional derivative is obtained and an application of the fractional damped model is also given. Additionally, the effects of a fractional damping term on a dynamical system are investigated. Finally, it is seen that the fractional derivative also has an effect on damping as a result of the previous studies in the literature.
2 The equation of motion
The problem of giving the longitudinal vibration of a harmonic external forced beam is given by
where
where Ω is the frequency of a beam [12]. Introducing the dimensionless parameters as
we have the new dimensionless parameters
where ρ is density, A is the crosssectional area, and L is the length of the beam. Thus, the equation in the nondimensional form is presented as
where
3 The method of multiple scales
In this section, an approximate solution will be searched by using the method of multiple scales. This method is known as the directperturbation method which can be applied directly to the partial differential equation. In higherorder schemes and for finite mode truncations, the method yields better approximations to the real problem [13]. Let us consider the expansion
where
where
At order one, the solution is obtained as
where
and
Thus, the solutions of Eqs. (17) and (19) are
where
Thus, different cases arise depending on the numerical value of variation frequency. These cases will be treated in the following sections.
4 Case studies
In this section, we assume that one dominant mode of vibrations exists. As a result of the previous studies in the literature, it is seen that the results are the same in the finite mode analysis and in the infinite mode analysis [3,14]. Therefore, we consider one dominant mode of vibration in this study.
4.1
Ω
1
close to 0,
Ω
2
away from
ω
n
(
Ω
1
≅
0
,
Ω
2
≠
ω
n
)
For this case, we consider the case of the nearness of
where
where cc and NST denote complex conjugates and nonsecular terms, respectively. Thus, the solution of Eq. (25) is
where the first term is related to the secular terms and the second term is related to the nonsecular terms. Now, substituting Eq. (26) into Eq. (25), we obtain the equation
with the boundary conditions
Using the solvability condition [15], we then find
Thus, by the normalization given as
where
Then, the amplitude solution for the first order of the problem is as follows:
and the displacement is also obtained:
where
As seen in Figure 1, the fractional derivative αorder has an effect on the displacementtime curves. In Figure 2 and Figure 3, the effect of the variation of the coefficient λ is observed for the different functions on displacementtime curves.
Figure 1. Displacementtime curves for different values of the order of the fractional derivative
for
Figure 2. Displacementtime graph for different values ofλfor
Figure 3. Displacementtime graph for various values ofλfor
4.2
Ω
1
close to
2
ω
n
,
Ω
2
away from
ω
n
(
Ω
1
≅
2
ω
n
,
Ω
2
≠
ω
n
)
If we consider the parametric resonance, then
Hence, the solvability condition requires that
where
where
Substituting Eq. (38) into Eq. (37) and also obtaining the result placed into Eq. (36) (and separating into real and imaginary parts), we get
For a nontrivial solution (
Here, λ also must be zero for the steadystate condition. Thus, the stability boundaries are determined as follows:
Inserting
for the external excitation frequency. Thus, the two different values of
The variation of an unstable region for different values of λ is observed in Figure 4. Since the rigidity of the system is increased by decreasing the value of λ, the unstable region reduces expeditiously for smaller values of λ.
Figure 4. Stability boundaries for different values ofλfor
The variation of an unstable region with some different values of α for
Figure 5. Stability boundaries for different values ofαfor
Figure 6. Critical value ofaversus the value ofηfor various fractional orders (
Figure 7. Critical value ofaversus the value ofαfor various damping coefficients (
4.3
Ω
1
away from
2
ω
n
and 0,
Ω
2
away from
ω
n
(
Ω
1
≠
2
ω
n
, 0,
Ω
2
≠
ω
n
)
This case corresponds to the absence of any resonances. Then, Eq. (23) turns into
where cc is a complex conjugate and NST denotes nonsecular terms. Substituting Eq. (26) into Eq. (43), we obtain the equation
with the boundary conditions
Using the solvability condition [15], we find
By the normalization, then Eq. (46) becomes
Thus, the displacement is obtained as follows:
On the other hand, the amplitude is
The displacementtime variation for different values of α is seen in Figure 8. Also, it is shown that the damping increases while the value of coefficient λ diminishes in Figure 9.
Figure 8. Displacementtime graph for the different fractional order for
Figure 9. Displacementtime graph for different values ofλfor
4.4
Ω
1
away from
2
ω
n
and
Ω
2
close to
ω
n
(
Ω
1
≠
2
ω
n
≠
0
,
Ω
2
≅
ω
n
)
This case deals with the primary resonance
and substituting Eq. (50) into the equation below,
where
we then obtain
where
into Eq. (53), we find
By the same mathematical manipulation, the stability boundaries are calculated as follows:
4.5 Sum type of resonance (
Ω
1
+
Ω
2
≅
ω
n
+
ε
σ
n
)
In this case, we consider the sum or difference of internal and external forced frequency
since
where
Substituting Eq. (50) into Eq. (58) and also separating the equation into real and imaginary parts, we get
Inserting Eq. (54) into Eqs. (60) and (61), then we have
Therefore, the stability boundaries are obtained as follows:
5 Conclusion
In this study, the effects of the damping term modeled with a fractional derivative on the dynamic analysis of a beam having viscoelastic properties subject to the harmonic external force are investigated. The parametric or primary resonances in simple supported beams, the governing equation of which involves a fractional derivative, are also analyzed. It is concluded that the value of the natural frequency of the beam modeled with a fractional damper is greater than that of the beam modeled with a classical damper. The fractional derivative has no effect on the static behavior, but it has a significant impact on the dynamic behavior. Furthermore, it is seen that the unstable region in the resonance case diminishes when the order of the fractional derivative increases.
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
All authors read and approved the final manuscript.
References

Rossikhin, YA, Shitikova, MV: Application of fractional calculus for dynamic problems of solid mechanics: novel trends and recent results. Appl. Mech. Rev.. 63, (2010) Article ID 010801

Fung, RF, Huang, JS, Chen, YC, Yao, CM: Nonlinear dynamic analysis of the viscoelastic string with a harmonically varying transport speed. Comput. Struct.. 66(6), 777–784 (1998). Publisher Full Text

Pakdemirli, M, Boyacı, H: The directperturbation methods versus the discretizationperturbation method: linear systems. J. Sound Vib.. 199(5), 825–832 (1997). Publisher Full Text

Galucio, AC, Deu, JF, Ohayon, R: A fractional derivative viscoelastic model for hybrid activepassive damping treatments in time domain  application to sandwich beams. J. Intell. Mater. Syst. Struct.. 16, 33–45 (2005). Publisher Full Text

Podlubny, I: Fractional Differential Equations, Academic Press, San Diego (1999)

Agrawal, OP: Formulation of EulerLagrange equations for fractional variational problems. J. Math. Anal. Appl.. 272, 368–379 (2002). Publisher Full Text

Zhang, Y: A finite difference method for fractional partial differential equation. Appl. Math. Comput.. 215, 524–529 (2009). Publisher Full Text

Chen, LQ, Zhao, WJ, Zu, JW: Transient responses of an axially accelerating viscoelastic string constituted by a fractional differentiation law. J. Sound Vib.. 278, 861–871 (2004). Publisher Full Text

Boyacı, H, Pakdemirli, M: A comparison of different versions of the method of multiple scales for partial differential equations. J. Sound Vib.. 204(4), 595–607 (1997). Publisher Full Text

Shimizu, N, Zhang, W: Fractional calculus approach to dynamic problems of viscoelastic materials. JSME Int. J. Ser. C Mech. Syst. Mach. Elem. Manuf.. 42(4), 825–837 (1999)

Rossikhin, YA, Shitikova, MV: Application of fractional calculus for analysis of nonlinear damped vibrations of suspension bridges. J. Eng. Mech.. 124, 1029–1036 (1998). Publisher Full Text

Zhang, L, Zu, JW: Nonlinear vibration of parametrically excited moving belts, part I: dynamic response. J. Appl. Mech.. 66(2), 396–403 (1999). Publisher Full Text

Öz, HR, Pakdemirli, M, Boyacı, H: Nonlinear vibrations and stability of an axially moving beam with time dependent velocity. Int. J. NonLinear Mech.. 36, 107–115 (2001). Publisher Full Text

Pakdemirli, M, Boyacı, H: Comparison of directperturbation methods with discretizationperturbation methods for nonlinear vibrations. J. Sound Vib.. 186, 837–845 (1995). Publisher Full Text

Nayfeh, AH: Introduction to Perturbation Techniques, WileyInterscience, New York (1981)