Abstract

In this paper, we study the stochastic P-bifurcation problem for axially moving of a bistable viscoelastic beam with fractional derivatives of high order nonlinear terms under Gaussian white noise excitation. First, using the principle for minimum mean square error, we show that the fractional derivative term is equivalent to a linear combination of the damping force and restoring force, so that the original system can be simplified to an equivalent system. Second, we obtain the stationary Probability Density Function (PDF) of the system’s amplitude by stochastic averaging, and using singularity theory, we find the critical parametric condition for stochastic P-bifurcation of amplitude of the system. Finally, we analyze the types of the stationary PDF curves of the system qualitatively by choosing parameters corresponding to each region within the transition set curve. We verify the theoretical analysis and calculation of the transition set by showing the consistency of the numerical results obtained by Monte Carlo simulation with the analytical results. The method used in this paper directly guides the design of the fractional order viscoelastic material model to adjust the response of the system.

1. Introduction

Fractional calculus is a generalization of integer-order calculus, it extends the order of calculus operation from the traditional integer order to the case of noninteger order, and it has a history of more than 300 years as so far. Due to the limitation of the definition of integer-order derivative, it cannot express the memory property of viscoelastic substances. The definition of fractional derivative contains convolution, which can well express the memory effect and show the cumulative effect over time. Compared with the traditional integer-order calculus, fractional calculus has more advantages and is a suitable mathematical tool for describing the memory characteristics [111] and in recent years, it has become the powerful mathematical tool in many disciplines, especially in the study of viscoelastic materials.

The fractional derivative can accurately describe the constitutive relation of viscoelastic materials with fewer parameters, so the studies of fractional differential equations on the typical mechanical properties and the influences of fractional order parameters on the system are very necessary and have important significance. In recent years, many scholars have done a lot of work and achieved fruitful results in this field: Li and Tang studied the nonlinear parametric vibration of an axially moving string made by rubber-like materials, a new nonlinear fractional mathematical model governing transverse motion of the string is derived based on Newton’s second law, the Euler beam theory, and the Lagrangian strain, and the principal parametric resonance is analytically investigated via applying the direct multiscale method [12]. Liu et al. introduced a transfer entropy and surrogate data algorithm to identify the nonlinearity level of the system by using a numerical solution of nonlinear response of beams, the Galerkin method was applied to discretize the dimensionless differential governing equation of the forced vibration, and then the fourth-order Runge-Kutta method was used to obtain the time history response of the lateral displacement [13]. Liu et al. investigated the stochastic stability of coupled viscoelastic system with nonviscously damping driven by white noise through moment Lyapunov exponents and Lyapunov exponents, obtained the coupled Itô stochastic differential equations of the norm of the response and angles process by using the coordinate transformation, and discussed the effects of various physical quantities of stochastic coupled system on the stochastic stability [14]. Nutting, Gemant and Scott-Blair et al. [1517] first proposed the fractional derivative models to study the constitutive relation of viscoelastic materials and the research on the viscoelastic materials with fractional derivative is also increasing, and so far, it is still a research hotspot [1825]. Rodr Guez et al. calculated the correlation function of transverse wave in linear and homogeneous viscoelastic liquid by the Generalized Langevin Equation (GLE) method and the influence of fractional correlation function on the dynamic behavior of the system is analyzed [26]. Bagley and Torvik used fractional calculus to study the dynamic behavior of viscoelastic damping structure and the responses of the system under general load as well as step load are analyzed respectively [27, 28]. Pakdemirli and Ulsoy studied the primary parametric resonance and combined resonance of the axial acceleration rope based on the discrete perturbation method and the multiscale method [29]. Zhang and Zhu analyzed the stability and dynamic response of viscoelastic belt under parametric excitation by the multiscale method [30, 31]. Chen et al. studied the dynamic behavior and steady-state response of axially accelerating viscoelastic beam by the Galerkin method [3235], derived the differential equation of nonlinear vibration for axially moving viscoelastic rope, and then pointed out that the damping of viscoelastic rope only exists in the nonlinear term [36, 37]. Leung et al. studied the steady-state response of a simply supported viscoelastic column under the axial harmonic excitation based on the fractional derivative constitutive model of cubic nonlinear and derived the generalized Mathieu-Duffing equation with time delay by the Galerkin discrete method, then the bifurcation behavior of the system caused by the order of the fractional derivative is analyzed [38]. Ghayesh and Moradian developed the Kelvin-Voigt viscoelastic model of the axially moving and the tensile belt, and then found the existence of nontrivial limit cycle in this system [39]. Liu et al. studied the dynamic response of an axially moving viscoelastic beam under random disorder periodic excitation, the first order expression of the solution is obtained by the multiscale method, and the stochastic jump phenomenon between the steady-state solutions is carried out [40]. Yang and Fang derived the system equation based on Newton’s second law and the fractional Kelvin constitutive relation and then studied the stability of the axially moving beam under the parametric resonance condition [41]. Leung et al. studied the single mode dynamic characteristics of the nonlinear arch with the fractional derivative, the steady-state solution of the system is obtained based on the residual harmonic homotopy method, and the influence of the parametric variation on the dynamic behaviors of the viscoelastic damping material is analyzed [42]. Galucio et al. obtained the fractional derivative model to describe the viscoelasticity of the system based on the Timoshenko theory and Euler-Bernoulli hypothesis and proposed a finite element formula for analyzing the sandwich beam of viscoelastic material with fractional derivative and the results were verified numerically [43].

Due to the complexity of fractional derivative, the analysis method of it becomes more difficult, the study on the vibration characteristics of the parameters can only be qualitatively analyzed, and the critical conditions of the parametric influences cannot be found, which affect the analysis and design of such systems, as well as the stochastic P-bifurcation of bistability for the viscoelastic beam with fractional derivatives of high order nonlinear terms under random noise excitation has not been reported. In view of the above situation, the nonlinear vibration of viscoelastic beam with fractional constitutive relation under Gaussian white noise excitation is taken as an example, the transition set curve of the fractional order system as well as the critical parametric condition for stochastic P-bifurcation of the system is obtained by the singularity theory, and then the types of stationary PDF curves of the system in each region in the parametric plane divided by the transition set are analyzed. By the method of Monte Carlo simulation, the numerical results are compared with the analytical results obtained in this paper, it can be seen that the numerical solutions are in good agreements with the analytical solutions, and thus the correctness of the theoretical analysis in this paper is verified.

2. Equation of Axially Moving Viscoelastic Beam

There are many definitions of fractional derivatives, and the Riemann-Liouville derivative and Caputo derivative are commonly used. The initial conditions corresponding to the Riemann-Liouville derivative have no physical meanings, however, the initial conditions of the systems described by the Caputo derivative have clear physical meanings and their forms are the same as the initial conditions for the differential equations of integer order. So in this paper, the Caputo-type fractional derivative is adopted as follows:where , , is the Euler Gamma function, and is the m order derivative of .

For a given physical system, due to the fact that the initial moment of the oscillator is , the following form of the Caputo derivative is often used:

In this paper, the transverse vibration of a viscoelastic simply supported beam under lateral excitation as shown in Figure 1 is considered; applying the d’Alambert principle, the governing equation can be written as [42]where is the mass density of the beam, is the area of cross-section, is the bending moment, is the lateral force, and is the horizontal force. From (3), we haveAssuming the material of the beam obeys a fractional derivative viscoelastic constitutive relation:where is the order of fractional derivative as is defined in (2), is the material modulus ratio, and is axial strain component.

When the deformation of the beam is small, the axial strain and lateral displacement satisfy the relationship as follows:Substituting (6) into (5) yieldsThe relationship between bending moment and axial stress of the beam can be expressed as follows:where is the thickness of the beam.

From (7) and (8), the expression of the bending moment can be obtained as follows:where .

The expression of the horizontal tension isSubstituting (9) and (10) into system (4), system (4) can be rewritten asThe boundary conditions areAccording to the boundary conditions (12), the solution of system (11) can be expressed as the Fourier series:Assume that the initial transverse vibration of the system is and the transverse load of the system satisfies the following form:where is the Gaussian white noise, which satisfies , , represents the noise intensity, and is the Dirac function.

By the discrete format based on Galerkin method, system (11) can be reduced to the fractional differential equation as follows:where is the   () order Caputo derivative of as is defined in (2), andFor convenience, system (15) can be represented as follows:where .

The fractional derivative term has contributions to both damping and restoring forces [4447], hence, introducing the following equivalent system:where , are the coefficients of equivalent damping and restoring forces of fractional derivative , respectively.

The error between system (17) and (18) isThe necessary conditions for minimum mean square error are [48]Substituting (19) into (20) yieldsAssume that the solution of system (18) has the following form:where is the natural frequency of system (17).

Based on (22), we can obtainSubstituting (22) and (23) into (21) yieldsTo simplify (24) and (25) further, asymptotic integrals are introduced as follows:Substituting (26) into (24) and (25) and averaging them across produce the ultimate forms of and as follows:Therefore, the equivalent system associated with system (18) can be expressed as follows:where and .

Next, we consider the stochastic P-bifurcation of system (28) which comprises the fractional derivatives of high order nonlinear terms and analyze the influence of parametric variation on the system response.

3. The Stationary PDF of Amplitude

For the system (28), the material modulus ratio is given as , coefficients of nonlinear terms are given as , , respectively, and nature frequency is given as . For the convenience to discuss the parametric influence, the bifurcation diagram of amplitude of the limit cycle along with variation of the fractional order is shown in Figure 2 when .

As can be seen from Figure 2, the solution corresponding to the solid line is almost completely coincided with the numerical solution, which proves the correctness and accuracy of the approximate analytical result of the deterministic system; at the same time, it shows that the solution corresponding to the solid line is stable and the solution corresponding to the dotted line is unstable. And it also can be seen that there is 1 attractor in the system when changes in the interval : equilibrium, as shown in Figure 3(a); there are 2 attractors when changes in the interval [0.6818, 1]: equilibrium and limit cycle, as shown in Figure 3(b).

In order to obtain the stationary PDF of the amplitude of system (28), the following transformation is introduced:where is the natural frequency of the equivalent system (28), and represent the amplitude process and the phase process of the system response, respectively, and they are all random processes.

Substituting (29) into (28), we can obtainin whichandEquation (30) can be regarded as a Stratonovich stochastic differential equation, by adding the corresponding Wong-

Zakai correction term: it can be transformed into the following Itô stochastic differential equation:where is a standard Wiener process andBy the stochastic averaging method, averaging (33) regarding , the following averaged Itô equation can be obtained as follows:whereEquations (35) and (36) show that the averaged Itô equation for is independent of and the random process is a one-dimensional diffusion process. The corresponding FPK equation associated with can be written asThe boundary conditions areThus, based on these boundary conditions (38), we can obtain the stationary PDF of amplitude as follows:where C is a normalized constant, which satisfiesSubstituting (36) into (39), the explicit expression for the stationary PDF of amplitude can be obtained as follows:where .

4. Stochastic P-Bifurcation

Stochastic P-bifurcation refers to the changes of the number of peaks in the PDF curve; in order to obtain the critical parametric condition for stochastic P-bifurcation, the influences of the parameters for stochastic P-bifurcation of the system are analyzed by using the singularity theory below.

For the sake of convenience, we can write as follows:in whichAccording to the singularity theory, the stationary PDF needs to satisfy the following two conditions:Substituting (42) into (44), we can obtain the following condition [49]:where represents the critical condition for the changes of the number of peaks in the PDF curve.

Substituting (43) into (45), we can get the critical parametric condition for stochastic P-bifurcation of the system as follows:where .

And amplitude satisfies

Taking the parameters as , , , and , according to (46) and (47), the transition set for stochastic P-bifurcation of the system with the unfolding parameters and can be obtained, as shown in Figure 4.

According to the singularity theory, the topological structures of the stationary PDF curves of different points () in the same region are qualitatively the same. Taking a point () in each region, all varieties of the stationary PDF curves which are qualitatively different could be obtained. It can be seen that the unfolding parameter plane is divided into two subregions by the transition set curve; for convenience, each region in Figure 4 is marked with a number.

Taking a given point () in each of the two subregions of Figure 4, the characteristics of stationary PDF curves are analyzed, and the corresponding results are shown in Figure 5.

As can be seen from Figure 4, the parametric region where the PDF curve appears bimodal is surrounded by an approximately triangular region. When the parameter () is taken in the region 1, the PDF curve has only a distinct peak, as shown in Figure 5(a); in region 2, the PDF curve has a distinct peak near the origin, but the probability is obviously not zero far away from the origin; there are both the equilibrium and limit cycle in the system simultaneously, as shown in Figure 5(b).

The analysis results above show that the stationary PDF curves of the system amplitude in any two adjacent regions in Figure 4 are qualitatively different. No matter the values of the unfolding parameters cross any line in the figure, the system will occur stochastic P-bifurcation behaviors, so the transition set curve is just the critical parametric condition for the stochastic P-bifurcation of the system, and the analytic results in Figure 5 are in good agreement with the numerical results by Monte Carlo simulation, which further verify the correctness of the theoretical analysis.

5. Conclusion

In this paper, we studied the stochastic P-bifurcation for axially moving of a bistable viscoelastic beam model with fractional derivatives of high order nonlinear terms under Gaussian white noise excitation. According to the minimum mean square error principle, we transformed the original system into an equivalent and simplified system and obtained the stationary PDF of the system amplitude using stochastic averaging. In addition, we obtained the critical parametric condition for stochastic P-bifurcation of the system using singularity theory; based on this, the system response can be maintained at the monostability or small amplitude near the equilibrium by selecting the appropriate unfolding parameters, providing theoretical guidance for system design in practical engineering, and avoiding the instability and damage caused by the large amplitude vibration or nonlinear jump phenomenon of the system. Finally, the numerical results by Monte Carlo simulation of the original system also verify the theoretical results obtained in this paper. We conclude that the order of fractional derivative and the noise intensity can both cause stochastic P-bifurcation of the system and the number of peaks in the stationary PDF curve of system amplitude can change from 2 to 1 by selecting the appropriate unfolding parameters (); it also shows that the method used in this paper is feasible to analyze the stochastic P-bifurcation behaviors of viscoelastic material system with fractional constitutive relation.

Data Availability

Data that support the findings of this study are included within the article (and its additional files). All other data are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the National 973 Project of China (Grant no. 2014CB046805) and the Natural Science Foundation of China (Grant no. 11372211).