Modeling of a Mass-Spring-Damper System by Fractional Derivatives with and without a Singular Kernel

In this paper, the fractional equations of the mass-spring-damper system with Caputo and Caputo–Fabrizio derivatives are presented. The physical units of the system are preserved by introducing an auxiliary parameter σ. The input of the resulting equations is a constant and periodic source; for the Caputo case, we obtain the analytical solution, and the resulting equations are given in terms of the Mittag–Leffler function; for the Caputo–Fabrizio approach, the numerical solutions are obtained by the numerical Laplace transform algorithm. Our results show that the mechanical components exhibit viscoelastic behaviors producing temporal fractality at different scales and demonstrate the existence of Entropy 2015, 17 6290 material heterogeneities in the mechanical components. The Markovian nature of the model is recovered when the order of the fractional derivatives is equal to one.


Introduction
In recent years, fractional calculus (FC) has been increasingly applied in different fields of science [1][2][3][4][5][6][7].Physical phenomena related to electromagnetism, propagation of energy in dissipative systems, thermal stresses, models of porous electrodes, relaxation vibrations, viscoelasticity and thermoelasticity are successfully described by fractional differential equations [8][9][10][11][12][13][14][15][16][17][18][19].FC allows for the investigation of the nonlocal response of mechanical systems, this is the main advantage when compared to the classical calculus.Some research concerning classical mechanics introduces FC, for example in [20], the authors investigated the relationship of the Fourier transform of fractional order to harmonic oscillation.Ryabov in [21] studied the fractional oscillator engaging the Riemann-Liouville fractional derivative.In [22], the authors found the analytical solution of the fractional damped oscillator equation using the Caputo derivative.In [23], Stanislavsky dealt with generalization of the mechanical oscillator using fractional derivatives.Tarasov in [24] considered the oscillator of fractional order as a system with memory; the author suggested independent harmonic oscillators to define the environment.Caputo in [25] studied the classic second order differential equation of the damped oscillator; using this model, the author quantifies the response and the decaying oscillations of a seismograph and applied the FC approach; this model introduces the mathematical memory operator represented by the fractional order derivative in order to model the response curves of more complex instruments realistically.Recently, in [26], the mass-spring and the spring-damper systems without the source term were analyzed; the authors considered fractional time derivatives of the Caputo type.Other applications of FC to mechanical oscillators are given in [27][28][29].
In the literature, a number of definitions of the fractional derivatives have been introduced, namely the Hadamard, Erdelyi-Kober, Riemann-Liouville, Riesz, Weyl, Grünwald-Letnikov, Jumarie and the Caputo representation [30][31][32][33].A thorough analysis of fractional dynamical systems is necessary to achieve an appropriate definition of the fractional derivative.For example, the Riemann-Liouville definition entails physically unacceptable initial conditions (fractional order initial conditions) [34]; conversely, for the Caputo representation, the initial conditions are expressed in terms of integer-order derivatives having direct physical significance [35]; this definition is mainly used to include memory effects.Recently, Michele Caputo and Mauro Fabrizio in [36] presented a new definition of the fractional derivative without a singular kernel; this derivative possesses very interesting properties, for instance the possibility to describe fluctuations and structures with different scales.Furthermore, this definition allows for the description of mechanical properties related to damage, fatigue and material heterogeneities.The properties of this new fractional derivative are reviewed in detail in the paper [37].Recently, Abdon Atangana and Badr Alkahtani [38,39] studied the resistor-inductor-capacitor electrical circuit and the Keller-Segel model based on the Caputo-Fabrizio definition; in both works, the authors showed the application of the new proposed fractional derivative without a singular kernel.
In the present work, we use both the Caputo and the Caputo-Fabrizio fractional derivatives to present the analytical and numerical solutions of the mass-spring-damper systems for different source terms; the idea proposed in [26] is applied to construct the fractional differential equations, and our representation preserves the physical units of the physical system for any value by the fractional derivative exponent.The main reason for this work is to investigate and model the displacement of the oscillator in fractal geometries using the fractional derivative with a singular kernel (Caputo approach) and the fractional derivative without a singular kernel (Caputo-Fabrizio approach).
The manuscript is organized as follows: Section 2 explains the basic concepts of the FC; Section 3 presents the analytic solution of the systems and the simulations with different source terms; and the conclusions are given in Section 4.

Basic Concepts
The Caputo derivative (CD) is defined as follows [32]: where t is a CD with respect to t, ϕ R is the order of the fractional derivative and Γ(•) represents the gamma function.
The Laplace transform of the CD has the form [32]: Some common Laplace transforms are: The Mittag-Leffler function [40] is defined by a power as: The Caputo-Fabrizio derivative (CF) is defined as follows [36,37]: where in this definition, the derivative of a constant is equal to zero, but unlike the usual Caputo definition (1), the kernel does not have a singularity at t = α.
The Laplace transform of ( 7) is defined as follows [36,37]: for this representation in the time domain, it is suitable to use the Laplace transform [36,37].
From this expression, we have:

Mass-Spring-Damper System
According to [26], to be consistent with the dimensionality of the physical equation, an auxiliary parameter σ is introduced into the fractional temporal operator: and: where γ represents the order of the fractional temporal operator and σ has the dimension of seconds.
The auxiliary parameter σ is associated with the temporal components in the system (these components change the time constant of the system) [41]; for the case γ = 1, Expressions ( 13) and ( 14) become ordinary temporal operators.Following this idea, the equation of the mass-spring-damper system represented in Figure 1 is given by: m and for the CF, we have: where the mass is m, the damping coefficient is β, the spring constant is k and F (t) represents the forcing function.

Mass-Spring System
Case 1: Consider a constant source, F (t) = f 0 , x(0) = x 0 , (x 0 > 0), ẋ(0) = 0; Equation ( 17) may be written as follows: where: is the fractional angular frequency and η 2 0 = k m is the frequency in the classical case.The solution of ( 21) is: where E 2γ is given by (6).Concerning the classical case, from (22), we have η 2 = η 2 0 = k m ; then, the solution of Equation ( 21) is: In this case, a physical relationship exists between Equation ( 17), the order γ and σ: For Expression (25), the solution ( 23) is given by: The plots for different values of the fractional order γ are shown in Figure 2.For Equation ( 18), we have: the Laplace transform ( 10) of ( 27) yields: applying the numerical inverse Laplace transform algorithm [42] to (28), we obtain the time response.The plots for different values of the fractional order γ are shown in Figure 3.  Case 2: Consider a periodic source, F (t) = f 0 cos(ωt), x(0) = x 0 , (x 0 > 0), ẋ(0) = 0; Equation (17) may be written as follows: where η 2 is given by (22).The solution of Equation ( 29) is: where E 2γ is given by (6).For the classical case, from (22), we have η 2 = η 2 0 = k m ; the solution of Equation ( 29) is: Taking into account Expression (25), the solution (30) is: The plots for different values of the fractional order γ are shown in Figure 4.For Equation (18), we have: the Laplace transform (10) of (33) yields: applying the numerical inverse Laplace transform algorithm [42] to (34), we obtain the time response.The plots for different values of the fractional order γ are shown in Figure 5.
For the classical case, Expression (37) becomes: In this case, there is a physical relation between Equation ( 19), the order γ and σ: considering Expression (39); the solution ( 37) is: The plots for different values of the fractional order γ are shown in Figure 6.For Equation ( 20), we have: the Laplace transform ( 10) of ( 41) yields: applying the numerical inverse Laplace transform algorithm [42] to (42), we obtain the time response.The plots for different values of the fractional order γ are shown in Figure 7. Case 2: Consider a periodic source, F (t) = f 0 cos(ωt), x(0) = x 0 , (x 0 > 0); Equation ( 19) may be written as follows: the solution for Equation ( 43) is: where E γ is given by (6).
For the classical case, Expression (44) becomes: considering Expression (39), the solution (44) is: the plots for different values of the fractional order γ are shown in Figure 8.For Equation ( 20), we have: the Laplace transform (10) of ( 47) is: applying the numerical inverse Laplace transform algorithm [42] to (48), we obtain the time response.
The plots for different values of the fractional order γ are shown in Figure 9. .Damper-spring system with periodic source, Caputo-Fabrizio derivative approach.

Mass-Spring-Damper System
Case 3: Consider F (t) = 0 and initial conditions equal to zero; Equation (15) may be written as follows: where η 2 and τ are given by ( 22) and (36), respectively.The solution of Equation ( 49) is: where E γ is given by (6).
For the classical case, Expression (50) becomes: In this case, there is a physical relation between Equation (15), the order γ and σ: considering Expression (52), the solution (50) is: The plots for different values of the fractional order γ are shown in Figure 10.For Equation ( 16), we have: and the Laplace transform ( 10) of (54) yields: where η 2 and τ are given by ( 22) and (36), respectively.Applying the numerical inverse Laplace transform algorithm [42] to (55), we obtain the time response.The plots for different values of the fractional order γ are shown in Figure 11.

Conclusion
In this work, we present alternative fractional differential equations using Caputo and Caputo-Fabrizio fractional derivatives for modeling oscillating systems.The proposed equations provide a universal behavior (see Equations ( 25), ( 39) and ( 52)); the solutions presented preserve the dimensionality of the studied system for any value of the exponent of the fractional derivative.The advantage of this alternative representation when compared to the models presented in the literature is the physical compatibility of the solutions.
For the classical Caputo approach, the solutions incorporate and describe long-term memory effects (attenuation or dissipation); these effects are related to an algebraic decay related to the Mittag-Leffler function; when γ is less than one, the fractional differentiation with respect to time represents the non-local displacement effects of dissipation of energy (internal friction) represented by the fractional order γ, which is related to the displacement of the oscillator in fractal geometries.For the Caputo-Fabrizio fractional derivative, the numerical solutions show a change of the amplitude variation and the phase; these behaviors depend on the fractional derivative order and modify the constant time; particularly, when the fractional operator is less than one, the systems exhibit a fast stabilization.This new fractional definition allows describing the relaxation phenomena characteristic of viscoelastic materials, and the numerical solution exhibits temporal fractality at different scales, as well as the existence of material heterogeneities in the mechanical components.The two definitions of fractional derivatives must apply conveniently depending on the materials treated; the choice of the fractional derivative depends upon the problem studied and on the phenomenological behavior of the system.For both definitions, when γ = 1, the displacement has the ordinary behavior.
The results gathered in Figures 2-11 briefly show that when γ = 1, the system displays the Markovian nature.However, for values of γ < 1, the equations describe systems that are non-conservative (non-local in time); in this context, the different γ values are modified the time constant of the systems and exhibit fractional structures (components that show an intermediate behavior between a system that is conservative and dissipative).In the range γ ∈ (0.85, 1), the figures show that the system presents dissipative effects that correspond to the nonlinear situation of the physical process (realistic behavior that is non-local in time).In these cases, the systems modified the damping capacity, for example when γ = 0.85, the damping capacity is bigger than when γ = 0.95.Furthermore, the figures show that the Caputo fractional derivative is more affected by the past compared to the Caputo-Fabrizio fractional derivative, which shows a rapid stabilization.
The resulting equations represent a generalization of the classical mass-spring-damper model; the proposed representation can be used to describe a wide variety of systems, which had not been addressed due to the limitations of the classical calculus.In this context, we consider that these results are useful to understand the behavior of dynamical complex systems, mechanical vibrations, control theory, relaxation phenomena, viscoelasticity, viscoelastic damping and oscillatory processes.

Figure 2 .
Figure2.Mass-spring system with a constant source, Caputo derivative approach.

Figure 3 .
Figure 3. Mass-spring system with a constant source, Caputo-Fabrizio derivative approach.

Figure 4 .
Figure 4. Mass-spring system with a periodic source, Caputo derivative approach.

Figure 6 .
Figure 6.Damper-spring system with a constant source, Caputo derivative approach.

Figure 7 .
Figure 7. Damper-spring system with a constant source, Caputo-Fabrizio derivative approach.

Figure 8 .
Figure 8. Damper-spring system with a periodic source, Caputo derivative approach.