#### Abstract

In this paper, dynamic response analysis of a forced fractional viscoelastic beam under moving external load is studied. The beauty of this study is that the effect of values of fractional order, the effect of internal damping, and the effect of intensity value of the moving force load on the dynamic response of the beam are analyzed. Constitutive equations for fractional order viscoelastic beam are constructed in the manner of Euler–Bernoulli beam theory. Solution of the fractional beam system is obtained by using Bernoulli collocation method. Obtained results are presented in the tables and graphical forms for two different beam systems, which are polybutadiene beam and butyl B252 beam.

#### 1. Introduction

Theory and applications of beams are very important research area due to its wide usage areas in applied sciences. Especially after starting the space adventure of the mankind, the demand to more resistant structures has great importance. Beams are generally modeled based on Euler–Bernoulli beam theory, which is called classical beam theory. The background of beam theory goes on Newton’s second law and some different aspects of beams, such as modeling, analysis of bending-buckling, and reinforcement and control, are hot topics of research papers since the beginning of the nineteenth century. The books can provide a general overview about the Euler–Bernoulli beam theory, please see [1–3]. Some important studies related to beams modeled in the sense of classical beam theory are also summarized as follows, but not limited to [4–15]. The beam systems in [1–15] have the integer order derivatives of the state function. In the beginning of 1930s, fractional derivative was introduced for describing the constitutive relation of some beam materials [16], and after 1980s, since fractional order equations have good memory and can be used to describe material properties more accurately with fewer parameters, they are considered to be good mathematical models for describing the dynamic mechanical behavior of materials [17]. In [18], the dynamic behavior of the thin plates resting on a fractionally damped viscoelastic foundation subjected to a moving point load is investigated and results show that the damping of the foundation system increases with increasing the order of the fractional derivative, which leads to a decrease in the dynamic response. In [19], the dynamic response spectra of fractionally damped viscoelastic beams subjected to concentrated moving load are presented and results reveal that with an increase in the order of the fractional derivative, the system damping of the system increases and the dynamic amplification factor (DAF) decreases, especially in the dynamic zone of the sweep parameter. In [20], the precise integration method (PIM) is extended to numerically integrate the equation of motion with fractional terms, which offers high accuracy and obtained numerical results indicate the viscoelastic dampers can enhance the seismic performance of structures significantly. In [21], the nonstationary free vibration and nonlinear dynamic behavior of the viscoelastic nanoplates are analyzed. Obtained results show that the viscoelastic model-based vibration is nonstationary unlike the elastic model. Moreover, the damping mechanism of the viscoelasticity is amplitude dependent and the contribution of the viscoelastic damping terms at higher forcing conditions becomes noticeable. On the other hand, several numerical methods are developed and employed for better analyzing the fractional mechanical systems. Widely used methods for fractional systems are finite element method [22], Galerkin method [23], variational iteration method [24], and multiscale method [25, 26]. Especially, papers existing in the literature, which include a solution method for analyzing the dynamic response of a fractional order beam system, can be shortly listed as [19, 27–30]. In [19], the authors combined Galerkin method and Newton–Raphson method for analyzing the vibration of a fractional beam equation and they compared the results for only seeing the effects of fractional or integer derivatives. In [27], the author considered the dynamic response analyzing of a fractional order viscoelastic beam by means of green function method. In [27], the author only compared the results based on changes on the fractional derivative between . In [28], the authors employed the Adomian decomposition method for solving a fractional beam equation and they only observed the effect of the order of fractional derivative. In [29], the authors used the dynamic green function method for analyzing the dynamic response in a fractional beam equation and the beam equation does not include the damping term. Results are simulated for only indicating the effects of order of fractional derivative. In [30], the author employed the green function method for a fractional viscoelastic beam system subjected to a base excitation. After obtaining the solution, the author compared the results corresponding to different fractional order derivative. By comparing the present study with the studies existing in the literature, objectives of the present study are expressed as follows:(i)In this paper, Bernoulli collocation method is firstly employed for analyzing the fractional viscoelastic beam equation. In the literature, especially for the fractional beam systems, green function method, Galerkin method, Newton–Raphson method, Adomian decomposition method, and Bernoulli collocation method in this paper were used, but by comparing these five methods, it is clear that Bernoulli collocation method is new and has less computational process and less work.(ii)In the literature, the authors only considered and discussed the effects of order of fractional derivative on the dynamic response. But, we discussed both the effects of the order of fractional order derivative and the effects of damping coefficient term and the effect of density of moving force load. So, it is said that the present study has wider perspective than other studies.(iii)Also, in the literature, results are obtained for one beam system. In this paper, effects of order of fractional derivative, effects of damping coefficient term, and the effect of density of moving force load are observed and compared for two different beam systems which are polybutadiene beam and butyl B252 beam.

For theoretical and experimental review about the fractional Euler–Bernoulli beams, please see [31]. Specifically, in the present paper, displacement analysis of a forced fractional viscoelastic beam is studied. External moving force load perfectly moves on the beam with the velocity from the left edge to the right edge of the beam. The solution of the fractional beam system is obtained by means of Bernoulli collocation method. The main advantage of the Bernoulli collocation method is that employing the Bernoulli polynomials is easier than Chebyshev, Bessel polynomials, and Haar wavelets [32–34]. These advantages of Bernoulli polynomials provide us for obtaining the solution by making less computational process in shorter time. In the step of employing the Bernoulli collocation method, some external moving force loads having different load intensities are considered and also the effects of internal damping and fractional order of the derivative are searched for a fractional beam system. In the simulations, two different beam systems, which are polybutadiene beam and butyl B252 beam, are taken into account for being compared each other in the aspects of internal damping effects and resistance to effect of external moving force. Comparison results of the beam systems are presented in tables and graphics. The rest of the paper is organized as follows: in the next section, definition of the displacement analysis problem for a fractional viscoelastic beam is presented and scheme of the beam is overviewed. In the third section, short definition of the fractional derivative in the Caputo sense is introduced. In the fourth section, Bernoulli collocation method is explained and adopted to the present problem. In the fifth section, obtained results are given and discussions are made in the light of employing the Bernoulli collocation method to fractional viscoelastic beam system.

#### 2. Definition of the Problem

The motion equation of the fractional viscoelastic homogeneous beam is obtained by considering the Euler–Bernoulli beam theory by ignoring shear deformation factor and rotary inertia of the beam. The beam is considered as a uniform viscoelastic beam and mechanical energy dissipation inside the beam is modeled by fractional order differential equations. By taking into account the [35], stress-strain constitutive relation of a fractional viscoelastic beam is given as follows:in which is the Young’s modulus of the viscoelastic beam, is the damping coefficient, and is the fractional derivative operator with the order with respect to . The simply supported viscoelastic beam initially is at rest and nondeformed. The beam is subjected to a horizontally moving constant force load with the velocity from the left edge to right edge of the beam, respect to axis. In the light of [27], let us introduce the formulation of a fractional viscoelastic beam structure illustrated in Figure 1.in which is the deflection of the viscoelastic beam in , is the time variable, is the final time observed duration, is the space variable, is the length of the viscoelastic beam, is the cross-section area of the structure, is the material mass density of the viscoelastic beam, is the axial moment of inertia of the beam, is a constant showing intensity of the external moving force load, is the Dirac-delta function, and is the velocity of the moving force load with the condition . Equation (2) is subjected to the following boundary conditions:and the following initial conditions:in which . means to square-integrable functions space in the manner of Hilbert in the domain in the Lebesgue sense with the following norm and inner product:

Let us assume that

After substituting the equations (6) into (2) and multiplying both sides of equation (2) with , integrating on , we obtain the following ordinary differential equation as follows:

Equation (7) is subjected to the following initial conditions:

#### 3. The Fractional Derivative in the Caputo Sense

*Definition.* The Caputo definition of the fractional-order derivative iswhere is the order of the derivative and is the smallest integer greater than . For the Caputo derivative, we have

#### 4. Bernoulli Collocation Method

The recurrence relation of the Bernoulli polynomials is defined by

For , . The first few Bernoulli polynomials are

Our goal is to get the approximate solution as the truncated Bernoulli series defined bywhere denotes the Bernoulli polynomials; are the unknown coefficients for Bernoulli polynomial, and is any positive integer which possess . Let us assume that linear combination of Bernoulli polynomials equation (14) is an approximate solution of equation (7). Our purpose is to determine the matrix forms of equation (7) by using (14). Firstly, we can write Bernoulli polynomials (12) in the matrix formwhere , , , and

The matrix form of equation (14) by a truncated Bernoulli series is given by

By using equations (15) and (17), the matrix relation is expressed aswhere

By using equation (18), we obtain the following relation:

By substituting the Bernoulli collocation points given byinto equation (21), we obtainand the compact form of the relation (23) becomes

In this way, the unknown Bernoulli coefficients are obtained by solving the system. Then, these coefficients are substituted into (14), and the approximate solution is obtained. For more details, see [36].

#### 5. Simulation Results and Discussion

Bernoulli collocation method for obtaining the solution of fractional viscoelastic beam equation is employed. Hence, displacement analysis of a forced fractional viscoelastic beam is investigated by taking into account the different moving force loads, different values of internal damping coefficient, and different values of fractional order of derivative. Obtained results are simulated and presented in the tables and graphical forms. The velocity, from left to right, of the external moving force on beam is considered as . In order to observe the dynamic response of the viscoelastic beams under the different intensity of external moving force, the intensity constant of the external moving force load on the beam is involved to computation as . Also, the values in Tables 1–6 are computed on , which is the middle point of the fractional viscoelastic beams. Observed duration of time is . In the first case, forced displacement analysis of a polybutadiene beam is observed for different values of moving force load and results are presented in Table 1. The length and material density of the fractional viscoelastic beam are taken into account as and , respectively. The cross-sectional area is , moment of inertia is , and Young’s modulus is for a fractional viscoelastic polybutadiene beam. Also, the order of fractional derivative is evaluated as 0.528 for the results in Figure 2 and Tables 1 and 3. By observing Figure 2, it is concluded that while the intensity of the external moving load force increases, namely, is 1 to 25 and 50, the displacement of the fractional viscoelastic polybutadiene beam also increases. Also, parallel observation results to Figure 2 are obtained by taking into account Table 1. For example, on the moment , the amount of the displacement of the polybutadiene beam is measured as 0.001 1 for , 0.028 for , and 0.056 for . This observation is valid the entire time interval for polybutadiene beam. Also, the effect of internal damping on the displacement is presented in Table 3 for polybutadiene beam. The internal damping coefficient is evaluated from 0.2 to 1 and by examining Table 3; it reveals that when internal damping coefficient increases, the displacement of the polybutadiene beam decreases under the same conditions. The effect of the fractional order to system is observed from Table 5 and it can be concluded that while increasing the values of the fractional derivative, the value of the displacement is decreasing. In the second case, a butyl B252 beam is taken into account by the coefficients; the cross-sectional area is , moment of inertia is , and Young’s modulus is . The order of fractional derivative is considered as 0.519 for Figure 3 and Tables 2 and 4. By checking Figure 3, it is easy to see that displacements corresponding to much bigger intensity of moving force load are much bigger. For example, on the moment , while to , corresponding displacements are calculated as 0.00089, 0.022, and 0.044, respectively. This is effective along the observation duration. In Table 4, some results related to the effect of internal damping are presented and internal damping coefficient is included in the computation as 0.2 to 1. After looking at Table 4, it is concluded that while internal damping coefficient decreases, the displacement of the butyl B252 beam increases and relation between the effects of internal damping and displacements is inversely proportional. The relation between the displacement and fractional order in the system is vice versa. As understood from Table 6, while decreasing the values of the fractional derivative, the value of the displacement is increasing. These observation results of the present study are also compatible with the results existing in the literature. By taking into account Tables 1–6 and Figures 2 and 3 and comparing these two kinds of fractional viscoelastic beams, it is seen that the polybutadiene beam has more greater displacements than butyl B252 beam under same conditions. Also, the effect of internal damping coefficient is more visible on the butyl B252 beam according to polybutadiene beam. These observations make clear that butyl B252 beam is stronger and preferable than the polybutadiene beam.

**(a)**

**(b)**

**(c)**

**(a)**

**(b)**

**(c)**

#### 6. Conclusion

In this study, the Bernoulli collocation method as a new solution method for obtaining the approximate solution of a fractional viscoelastic beam model subjected to moving force load is employed. Dynamic response analysis of the fractional viscoelastic beam model is investigated for two different specific beams: polybutadiene beam and butyl B252 beam. Displacement analysis of a point on the fractional viscoelastic beams is studied for different moving force loads and also effect of the internal damping to displacement is observed for different internal damping coefficients. Moreover, dynamic response of the fractional viscoelastic beam is examined for different values of the fractional order. Obtained results are presented in tables and graphics and results reveal that Bernoulli collocation method is very effective and powerful solution method for obtaining the solution of fractional order viscoelastic beam models. After observing Figures 2 and 3, it is easy to conclude that as the moving force load increases, the displacement of a point on the beams also increases. Also, numerical results, presented in Tables 1–4, show that under the same moving force load with the same internal damping effect, the displacement of a point on the polybutadiene beam is greater than that corresponding to butyl B252 beam. Moreover, under the same moving force load, changes in the displacements of a point on the beams are examined in the aspect of different internal damping effects and observations made clear that butyl B252 beam better reflects the effect of internal damping to displacement of a point on the beam. By comparing polybutadiene beam and butyl B252 beam, it is concluded that polybutadiene beam is more open to destructive effects of vibrations under the same conditions with the butyl B252 beam.

#### Data Availability

The data used to support the findings of this study are included within the article.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Authors’ Contributions

The authors completed this study and wrote and approved the final version of the manuscript.