**4. Numerical analysis of the model**

All the parameters concerning the chosen models of beam, of foundation and of the aerodynamic force are presented in Ref. [21]. Theses parameters clearly help to calculated the dimensionless coefficients defined in Eq. (7). It is well know that the validation of the results obtained through analytical investigation is guarantee by the perfect match with the results obtained through numerical simulations. thus, The numerical scheme used in this chapter is based on the Grunwald-Letnikov definition of the fractional order derivative Eq. (39). [37–40] and the Newton-Leipnik algorithm [37, 38]

$$D\_{\pi}^{a}\left[\chi\left(\boldsymbol{\tau}\_{\pi\_{f}}\right)\right] \approx h^{-a} \sum\_{l=0}^{n\_{f}} \mathbf{C}\_{l}^{a} \chi\left(\boldsymbol{\tau}\_{\pi\_{f}-l}\right),\tag{39}$$

comparison between the results from the mathematical analysis (curve with dotted line) and the results obtained from numerical simulation of Eq. (9) using Eq. (39). The match between the results shows a good level of precision of the approximation made in obtaining Eq. (24). This figure also reveals that the vibration amplitude of the beam decreases and the resonance frequency of the system increases as the number of bearings increases. **Figure 2(b)**, shows the effect of loads number on the beam response. It is observed that as the value of *Nv* increases, the amplitude of

*Vibrations of an Elastic Beam Subjected by Two Kinds of Moving Loads and Positioned on…*

Looking at the effects of the order of the fractional derivative *α* on the amplitude of the beam, we obtain the graph of **Figure 3**. Small value of *α* leads to large value of the maximum vibration amplitude. It is also clearly shown that the system is more stable for the highest order of the derivative. The multivalued solution appears for small order and disappears progressively as the order increases. The resonance (a

The good match between the analytical and the numerical results gives a validation

*Effect of the fractional order on the amplitude of the beam A when the dimensionless stiffness k*<sup>0</sup>

*Analytical curves (-), numerical curve (o). All The parameters are given in [21].*

Also, the stochastic analysis has allowed to estimate the probabilistic distribution as a consequence of the wind random effects. The beam response, and more specifically the stationary probability density function *Ps*ð Þ *a* of its amplitude *a*, can also be retrieved (**Figures 4** and **5**). This type of analysis indicates that as the additive wind turbulence parameter increases, the peak value of the probability density function

<sup>0</sup> increases, see **Figure 3(a)-(d)**.

*<sup>o</sup> varying:*

vibration at the resonant state merely increases.

*DOI: http://dx.doi.org/10.5772/intechopen.96878*

peak of the amplitude) appears as the parameter *k*<sup>0</sup>

of the approximations made.

**Figure 3.**

**23**

where *h* is the integration step and the coefficients *C<sup>α</sup> <sup>l</sup>* satisfy the following recursive relations:

$$\mathbf{C}\_0^a = \mathbf{1}, \quad \mathbf{C}\_l^a = \left(\mathbf{1} - \frac{\mathbf{1} + a}{l}\right) \mathbf{C}\_{l-1}^a. \tag{40}$$

Now we display in some figures the effects of the main parameters of the proposed model. For example, **Figure 2(a)** shows the effect of the number of the bearings on the amplitude of vibration of the beam. This graph also shows a

#### **Figure 2.**

*Effects of the number of the bearings Np (a) and number of moving loads Nv (b) on the amplitude response of the beam when driving frequency Ω: Analytical curves (-), numerical curves (o). All The parameters are given in [21].*
