Vibration Control of Fractionally-Damped Beam Subjected to a Moving Vehicle and Attached to Fractionally-Damped Multiabsorbers

This paper presents the dynamic response of Bernoulli-Euler homogeneous isotropic fractionally-damped simply-supported beam. The beam is attached to multi single-degree-of-freedom (SDOF) fractionally-damped systems, and it is subjected to a vehicle moving with a constant velocity. The damping characteristics of the beam and SDOF systems are described in terms of fractional derivatives. Three coupled second-order fractional differential equations are produced and then they are solved by combining the Laplace transformwith the decompositionmethod.The obtained numerical results show that the dynamic response decreases as (a) the number of absorbers attached to the beam increases and (b) the damping-ratios of used absorbers and beam increase. However, there are some critical values of fractional derivatives which are different from unity at which the beam has less dynamic response than that obtained for the full-order derivatives model. Furthermore, the obtained results show very good agreements with special case studies that were published in the literature.


Introduction
Many physical problems were adequately described by differential equations of fractional order.The list of these problems is long and the areas of applications are broad.Polymeric damping, fluid mechanics, and the theory of viscoelasticity [1][2][3][4][5][6][7][8] are some of these applications.The solution in the closed form to a class of SDOF fractional oscillators subjected to an impulse response is investigated [9].Relationship between fractional oscillator processes and the corresponding fractional Brownian motion processes is well established in research.However, the increment process of fractional Brownian motion can be described by a random function with long-range dependence fractional Gaussian noise model [10].A fractal time series is taken as the solution of differential equation of a fractional order, or a response of a fractional system, or a fractional filter driven with a white noise in the domain of stochastic processes [11,12].A general approach for approximating ideal filters that are based on fractional calculus from the point of view of systems of fractional order was introduced [13,14].A new direct operational inversion method is introduced for solving coupled linear systems of ordinary fractional differential equations, where the obtained solutions are expressed explicitly in terms of multivariate Mittag-Leffler functions [9,15,16].The solutions of fractional Langevin equation of two different orders which are known as the fractional Ornstein-Uhlenbeck processes, based on Weyl and Riemann-Liouville fractional derivatives were obtained [9,17].Power-laws-type data may further, in earthquake and wind engineering applications, be governed by stochastically differential equations of fractional order [18,19].Analysis of signals generated by fractional-order process and applications of such model in a wide range of physical, mechanical, and biological systems are discussed in [20].However, as noted from the literature, a clear physical interpretation of the fractional derivative has been elusive.Good examples are (1) transverse vibration of homogeneous beams whose damping behavior is described by a fractional derivative of predefined Advances in Mathematical Physics constant orders was an important subject of matter [11,[21][22][23][24]; (2) diffusion process where fractional differential equations have been employed to describe an anomalous diffusion regime including both subdiffusion and superdiffusion [25]; and (3) many linear viscoelastic damping materials exhibit a macroscopic constitutive behavior which has been the subject of many investigations involving fractional order derivative [26].However, in the fractional-derivative model of viscoelastic materials over extended ranges of time and frequency, the deformation work corresponding to springs and losses corresponding to dashpots (described by fractional derivatives) have both energy types (stored and dissipated) at any point of the utilized materials [27,28].
The dynamic response characteristics of beams resulting from the passage of different moving load systems, and with the existence of dampers and oscillators, have been studied within many applications.Some of the applications consist of the vibrations existing in bridges, aircraft carriers, and railroad tracks due to moving vehicles, and vibrations existing in pipeline systems due to fluid flows [29][30][31].One of the main objectives of studying such problems is to minimize or to suppress the dynamic response characteristics or vibration levels of beams when subjected to different load systems.The simplest and most economical way to achieve this objective is by having different lumped appendages attached to the beam at different locations.In investigating the dynamic response of beams, if the travelling force is constant, then the time-dependent coefficients can be obtained in analytical forms [32] which make it easy to find a solution to any particular problem.However, in many cases, it is not just the response of the beam (which is a function of the spatial variable and time and depends on several parameters such as beam characteristics, the magnitude of the force, its velocity, number of appendages attached to the beam, and their locations) that the designer is interested in.However, it has long been observed that, as a structure is subjected to moving loads, the induced dynamic deflection and stress can be significantly higher than those observed in the static case [33].In this regard, the majority of the literature has been devoted to the study of the so-called moving force [34], moving mass [35][36][37][38], moving oscillator problems [39][40][41], and multiple moving oscillators [42].
The vibration of a beam excited by a moving unmaintained oscillator is discussed in [43].In this study, for certain values of system parameters and travelling speeds of the moving oscillator, separation of the oscillator from the beam can occur frequently and even more than once during the travel of the load [43].The problem of an oscillator traversing on an elastically-supported continuum is studied, whereas the flexibility in the boundaries of the continuum is modeled by linear-transverse springs [33].Green functions formulation to analyze the free vibration of a linear undamped beam subjected to grounded SDOF oscillators appendages at different discrete points is used [44], where researchers compared their results with the approximate solution obtained by Galarkin's method.
Searching the relevant literature shows only few papers dealing with the vibration of beams with constant values for the fractional damping behavior [23,24,45].However, in this study, the fractional damping behaviors of the beam and absorbers are taken with arbitrary orders.Therefore, this study is important for the following key reasons: the damping behavior of the dynamic response of the beam can be analyzed over a continuous wide range by means of fractional-derivative parameters.Moreover, the fractional derivative damping behavior may better converge to realistic experimental models than the well-known first order damping model.Hence, the fractional derivative model may precisely prescribe a nonlinear damping behavior than other existing damping models.
In this paper, a homogeneous isotropic Euler beam with different appendages attached at different locations subjected to a uniformly moving vehicle is investigated.The solution to the handled problem is based on decomposing the generalized transverse displacements of the beam, attached oscillators, and the moving vehicle into infinite series components.Then, the generalized displacements are recursively solved by in the Laplace transform domain.The solutions in this work mainly focus on the following characteristics: the maximum response of the beam for a given velocity of the vehicle, the maximum response over all possible velocities and the velocity at which it occurs, the optimal beam and absorbers fractional-damping characteristics, and the effect of the number of attached absorbers.

Formulation of the Problem
Transverse vibration of Bernoulli-Euler homogeneous isotropic fractionally-damped simply-supported beam is investigated.The beam is assumed to be subjected to a vehicle moving from left-to-right with a constant velocity (V).The governing equation of the beam under investigation, whose damping characteristics are described by a fractional derivative of order , can be written as BC's:  (0, ) = 0,  (, ) = 0,   (0, ) = 0, IC's:  (, 0) = 0,  (, 0)  = 0, where (⋅) refers to the unit step Heaviside function and (⋅) refers to the Kronecker delta function.The equations corresponding to the -SDOF absorbers that are appendaged to the beam are given as Advances in Mathematical Physics IC's: The equations corresponding to the moving vehicle can be expressed as IC's: where the beam is assumed to be initially at rest and {, , , , } are the length, modulus of elasticity, moment of inertia of cross-sectional area (), mass per unit length, and the coefficient of external damping of the beam, respectively.Furthermore, these values are assumed to be constants.The -rigid-body SDOF fractionally-damped absorbers of order , masses M , stiffnesses K , and damping constants Ĉ are attached to the beam, as shown in Figure 1 (for  = 1, 2, . . ., ).These elements are located at distances   measured from the left-end of the beam.However, the parameters { V ,  V ,  V } refer to mass, stiffness, and damping constants of the vehicle, respectively.Furthermore, the set {(, ), ẑ (), ()} represents the beam transverse deflection at location () and time (), the transverse dynamic displacements of the masses M , and vehicle's transverse displacement, respectively.Equations that are similar to (1a), (1b), (1c), (2a), (2b), (3a), and (3b) are presented in [42], but with the conventional first order derivative damping behaviors.The fractional derivatives of order  or , presented in (1a) and (2a), may be defined in many ways [47].Among these are the two most frequently encountered, the Riemann-Liouville and Caputo fractional derivatives of order  which are defined by the following convolution integrals, respectively: where  is a positive integer such that ( − 1 ≤  < ).
In modal form, the transverse deflection of the simplysupported beam can be written as [31,32]: where   () is the generalized displacement or the modal response of the beam and   () are the normal modes of the undamped free vibration of the simply-supported beam written as [46]   () = sin (  ) , where (  = /) is the frequency parameter associated with the simply-supported beam.Equation ( 6) should satisfy the following fourth-order partial differential equation: Substituting ( 5) into (1a), (2a), and (3a) leads to the following coupled-system of fractional differential equations: Multiplying (8a) by   () and integrating over the domain (0 The dynamic response of the beam represented by (1a) is analyzed after projecting the 4th order partial derivative into a complete orthogonal basis.For the present problem, the eigenfunctions of the linear operator representing the simply-supported beam without appendages given in (6) can be used.These eigenfunctions should satisfy the following orthogonality conditions [32,42]: By considering these orthogonality conditions, the differential equation of the th mode of the generalized beamdisplacement, (9), can be rewritten as Substituting ( 6) into (8a), (8b), (8c), and (11) and noting that ( = ) lead to the following coupled system of fractional differential equations: where are the undamped natural circular frequency (  ), beam damping ratio (  ), th oscillator damping ratio ( ζ ), vehicle damping ratio ( V ), the natural frequency of the vehicle ( V ), the natural frequency of the th oscillator ( ω ), and the load frequency (Ω  ), respectively, where ( 1 = 2/).

Results and Discussion
In this Section, two examples are given to verify the validity of the formal technique used in this paper.Moreover, three more important case studies are introduced to show the effect of fractional-damping behaviors for both of the beam and the SDOF absorbers and the effect of the number of SDOF absorbers attached to the beam ().In all case studies, the dynamic response of the beam is mainly governed by the fundamental mode ( = 1) of the beam because this mode has the significant effect on the beam dynamic response.

Verification Problem I.
This example verifies the validity of the technique used in this paper.In this example, an undamped beam with a moving oscillator is investigated [49].The parameters adopted for the beam under investigation are given in Table 1.The displacement response for the beam versus the normalized time is shown in Figure 2; that is, when ( = 0) the load enters the beam from left-hand side ( = 0) and when ( = /V) the moving load system leaves the righthand side of the beam ( = ).The curves of Figure 2 show the dynamic mid-span displacement response of the beam for the five cases of traveling velocities: V = 10 m/s, V = 27.778m/s, V = 40 m/s, V = 60 m/s, and V = 80 m/s, respectively.
In this example, it is shown that the dynamic displacement response increases as the traveling speed increase.This follows from the fact that the energy content of the elastic beam increases when the speed increases.Furthermore, Figure 2 shows that the amplitude of the displacement response for the force and mass moving systems increases as speed increases.The maximum displacement when V = 27.778m/s is ( max = 2.17453 m/s) which satisfies the results obtained by Duhamel integral [32] and is almost the same as that presented in [49].

Verification Problem II.
This verification problem presents a damped beam attached to a single absorber and subjected to a moving load of intensity 9.81 N; that is, the moving vehicle of unit mass approaches to the moving load ( V = 0 N/m,  V = 1 Ns/m,  V = 1 kg).The data presented in [50,51] are used for solving this problem, see Table 1.In this problem, the response of an initially stationary beam subjected to a moving load with a constant velocity is Problem in Section 3.5  considered.The beam is assumed to be free or appendaged to one, two, or three absorbers ( = 0, 1, 2, 3) and ( = 1, . . ., ) at different locations of the beam, see Figures 3, 4, and 5.
In this problem, ( = 0) refers to a beam without any appendages; that is, all parametric properties of the oscillators are null.The case for ( = 1) was studied in [51] as  a two-SDOF system, while in [50], it was analyzed using the Euler-Bernoulli theory.In [50], the Galerkin-Bubnov variation method is used to solve the beam problem with only one absorber attached at the mid-span of the beam.The order of the fractional derivative for both of the beam and the oscillator is assumed to be unity; that is,  =  = 1.Figures 3-5 show the maximum response of the beam at the mid-span versus the velocity of the moving load for three different values of the beam damping ratio; that is,  = 0,  = 0.075, and  = 0.15, respectively.The upper two curves in Figure 3 fit precisely the published results in [50,51].The curve for ( = 0) in Figure 3 fits exactly the well-known exact solution which is obtained by the Duhamel integral [32].Furthermore, we note from Figures 3-5 that, as both the number of absorbers which are attached to the beam and as the damping beam-ratio increases, the dynamic response decreases.This is clearly seen from the lowest curve in Figure 5 where  = 3 and  = 0.15.

Optimization of the Vehicle Dynamic Response.
The moving load problem that was solved in Section 3.2 is solved here for the more general case of a moving vehicle ( V = 700 N/m,  V = 0.1,  V = 1 kg).Figures 6, 7, and 8 illustrate the maximum response of the beam at the mid-span versus the velocity of the moving vehicle for three different values of the damping ratio of the beam; that is,  = 0,  = 0.075, and  = 0.15, respectively.Moreover, Figure 9 shows the maximum central beam displacement versus the velocity of the moving vehicle for different numbers of absorbers ( = 1,2).The curves in these figures show that the maximum central displacement of the beam decreases as the number of absorbers, beam damping ratio, and vehicle velocity increase.By comparing Figures 3-5 obtained for a moving load with Figures 6-8 obtained for a moving vehicle, it is clearly seen that the maximum displacement of the beam decreases for the moving vehicle case and this is due to the effect of the damping ratio of the vehicle.Moreover, the       of fractional derivatives for both of the beam and absorbers (0 <  < 2, 0 <  < 2).The 3D surface shown in Figure 10 (for the data given in Table 1) represents the variation of the maximum deflection of the beam at the mid-span (where the absorber is attached) with the beam fractional derivative () and the fractional derivative of the absorber ().The obtained surface in Figure 10 leads to very interesting critical values of  and .From this surface, the optimal values obtained are  = 1.245 and  = 1.190, at which the minimum beam deflection occurs.This reduction is due to the effect of the damper added to the vehicle.The optimal values obtained from the surface shown in Figure 10 further show the physical meaning of fractional derivative, which is equivalent to increasing the damping characteristics of the beam.

Optimization of the Number of Fractionally-Damped
Oscillators.In this section, the optimal values that were obtained in Sections 3.  The effect of the previously obtained optimal values is observed clearly when both Figures 5 and 8 are compared with Figure 13.The percent relative reduction in the maximum central deflection (when the beam subjected to a moving load is compared to that subjected to a moving    vehicle) is computed for ( = 0,  = 1,  = 2, and  = 3) as (14.5%, 13.85%, 12.74%, and 11.63%), respectively, see Figures 5 and 8 Figure 11: The dynamic response ((/2, )) for the beam with  = 0, subjected to a moving vehicle and attached to many absorbers. = 1,  = 2, and  = 3) as (27.58%, 41.27%, 48.99%, and 54.54%), respectively, see Figures 8 and 13.This reduction in the dynamic response of the beam is due to the effect of the damping of the beam, fractional damping, number of absorbers, and the vehicle-damping.However, the influence of the order of fractional derivatives, beam damping ratio,  oscillators damping ratios, and the number of the attached oscillators on the beam dynamic response is clearly observed.Figure 14 shows the dynamic response of the beam versus time for both of the moving load and vehicle for the beam damping ratio ( = 0.15) for which ( = 0,  = 1,  = 2,  = 3).The fractional derivatives for both the beam and absorbers are taken as  = 1.245 and  = 1.19, respectively.The effect of these optimal values on the beam response is clearly observed.The reduction in the dynamic response of the beam is due to the effect of the damping of the beam, fractional damping, number of attached absorbers, and the vehicle-damping.

Conclusions and Future Work
In this paper, an analytical model is introduced to solve the transverse vibration of Bernoulli-Euler beam attached to -SDOF absorbers with damping characteristics that are described in terms of fractional derivatives.The damping characteristics of the beam and absorbers are assumed to be of fractional-order derivatives.The method employed uses the Laplace transform with decomposing the displacement into infinite series to find the solution.However, the solution may be written in a closed form for some special cases; otherwise, it may be truncated and produced by using appropriate mathematical software such as Maple.
The numerical results obtained in this paper, illustrate that (1) the dynamic response decreases as the number of absorbers attached to the beam increases, (2) as both the damping-ratios of absorbers and beam increase, the dynamic response decreases, (3) the type of load further affects the dynamic response of the beam; hence, a moving vehicle reduces the response due to its damping property, and (4) there are some critical values of  and  at which the beam has less deflection than that which is obtained from the fullorder derivative models (i.e.,  = 1 and  = 1).Finally, even though the fractional derivative model provides better damping models of fractionally-damped structures and materials, a limited number of papers have been presented in this field.This is probably due to the fact that the underlying mathematics of the fractional derivative models are sophisticated and not well-developed.Hence, developments in this field will lead to the development of the applications of fractional derivative models in several other engineering problems.
The approach of the fractional derivative damping, used in this paper, will allow researchers to choose suitable mathematical models that may better converge to realistic experimental models than the well-known conventional firstorder damping models.Hence, the fractional derivative models may precisely represent a nonlinear damping behavior better than other existing damping models.
Future work will involve the following items for the application of the fractional calculus-based approach that is used in this paper within: (1) quantum and nano computing systems [52], (2) nanobeams, materials and medical applications [52][53][54], (3) fractal-based formulations to the microquasistatic thermoviscoelastic creep for rough surfaces in contact [55], and (4) other beam configurations for various types of nonlinearities, damping characteristics, and moving vehicles.

Figure 1 :
Figure 1: A Simply-supported beam attached to -absorbers and subjected to a moving vehicle with a constant speed.

Figure 2 :
Figure 2: The dynamic response of an elastic beam at mid-span for many load speeds.

Figure 3 :
Figure 3: The maximum displacement of the beam at mid-span versus the velocity of the moving load for  = 0.

Figure 4 :
Figure 4: The maximum displacement of the beam at mid-span versus the velocity of the moving load for  = 0.075.

Figure 5 :
Figure 5: The maximum displacement of the beam at mid-span versus the velocity of the moving load for  = 0.15.
vehicle velocity (m/s)

Figure 6 :
Figure 6: The maximum displacement of the beam at mid-span versus the velocity of the moving vehicle for  = 0.

Figure 7 :
Figure 7: The maximum displacement of the beam at mid-span versus the velocity of the moving vehicle for  = 0.075. .

3. 4 . 3
Optimization of the Fractional Derivatives.In this Section, one example is analyzed to show the effect of the order Advances in Mathematical Physics 9 Absolute maximum central beam displacement (mm) Moving vehicle velocity (m/s)

Figure 8 :
Figure 8: The maximum displacement of the beam at mid-span versus the velocity of the moving vehicle for  = 0.15.

3 and 3 . 4 (
= 1.245,  = 1.19,V = 25 m/s) are used to determine the dynamic response of the beam attached to many SDOF systems.The other parameters considered in this example are shown in Table 1.Figures 11, 12, 13, and 14 show the dynamic response of the beam for various values of beam damping ratio ( = 0,  = 0.075,  = 0.15), respectively.The effects of the beam damping ratio and the number of oscillators are clearly observed from the curves in Figures 11-14 .

Figure 9 :
Figure 9: The maximum displacement of the beam at mid-span versus the velocity of the moving vehicle with one and two ( = 1, 2) absorbers attached.

Figure 10 :
Figure 10: The maximum central deflection of the beam against the order of the beam's fractional derivative and the fractional derivative of the absorber (, ).
, while the percent relative reduction in the maximum central deflection when Figure8( = 1 and  = 1) is compared with Figure13( = 1.245 and  = 1.19) of the beams subjected to a moving vehicle is computed for ( = 0,

Figure 12 :
Figure 12: The dynamic response ((/2, )) for the beam with  = 0.075, subjected to a moving vehicle and attached to many absorbers.

Figure 13 :
Figure 13: The dynamic response ((/2, )) for the beam with  = 0.15, subjected to a moving vehicle attached to many absorbers.