Analytical Analysis on Nonlinear Parametric Vibration of an Axially Moving String with Fractional Viscoelastic Damping

The nonlinear parametric vibration of an axially moving string made by rubber-like materials is studied in the paper.The fractional viscoelastic model is used to describe the damping of the string. Then, 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. Taking into consideration the fractional calculus law of Riemann-Liouville form, the principal parametric resonance is analytically investigated via applying the direct multiscale method. Numerical results are presented to show the influences of the fractional order, the stiffness constant, the viscosity coefficient, and the axial-speed fluctuation amplitude on steady-state responses. It is noticeable that the amplitudes and existing intervals of steady-state responses predicted by Kirchhoff ’s fractional material model are much larger than those predicted by Mote’s fractional material model.


Introduction
Axially moving structures are one of common elements in many mechanical systems, which is extensively used in engineering fields, such as paper sheets, magnetic tapes, power transmission belts, chains, fiber textiles, aerial tramways, pipes conveying fluids, and thread lines.In many cases, due to initial, parametric, and external excitations, the generation of unwanted transverse vibrations may limit their applications.Therefore, the dynamic behaviors of such devices have been widely investigated by numerous scholars for the past decades and are still of interesting today [1][2][3][4][5][6][7][8][9][10].Moreover, according to the demand of the actual engineering problems, some researchers have applied their theoretical results to design and optimize the axially moving structures [11,12].
In general, small imperfections induced by either geometrical or dynamic sources may bring about the occurrence of an unsteady axial speed or tension.For example, when an axially moving belt is installed on rotating pulleys, the torsional vibration of the pulleys would lead to a small fluctuation in the axial moving velocity, and consequently this system may exhibit more complicated dynamics behaviors like the parametric resonances.Thus, the parametric resonances of the axially moving structures caused by the pulsatile transport speed or tension have received extensive attention.Fung et al. [13] applied the Galerkin method and the numerical integration technique to investigate the parametric vibration of a viscoelastic string with the nonuniform transport speed.Pellicano et al. [14] used the approximate analytical and the experimental methods to study the primary and the parametric resonances of a power transmission belt with the fluctuation of the tension.Chen et al. [15] employed the averaging method to investigate the stability problems of an axially accelerating tensioned beam under the condition of the subharmonic and combination resonances.Pakdemirli and Öz [16] adopted a perturbation technique to discuss stable regions of a simply supported axially moving beams subjected to sum-and difference-type combination resonances.Ghayesh [17] utilized the method of multiple scales to study the stability characteristics for principal and combination parametric resonances of an axially moving string with the partial elastic support.
With the development of engineering technique, the various complicated materials are utilized to fabricate the axially moving structures.In order to better understand energy dissipation mechanism of such materials, some viscoelastic constitutive relations like the Kelvin-Voigt form have been used to describe mechanical behaviors of the materials, which is widely applied in axially moving continuums such as strings [18,19], beams [20,21], belts [22,23], and plates [24,25].In addition, some structures made of more complex viscoelastic materials are explored.Marynowski and Kapitaniak [26] presented a mathematical model of an axially moving viscoelastic beam with the three-parameter Zener element and investigated both regular and chaos motions using the Galerkin method.Chen et al. [27] analytically studied nonlinear parametric responses of an axially moving string composed of the complicated viscoelastic material based on the Boltzmann superposition principle.Wang and Chen [28] developed the differential quadrature scheme to determine the stable boundary of an axially moving viscoelastic beam with the standard linear solid in the case of the principal resonance.Ding and Chen [29] employed the material time derivative to characterize the viscoelastic property of an axially moving viscoelastic beam.They investigated the stability in principle resonance of the system by using analytical and numerical methods.
On the other hand, the axially moving structures made by rubber-like materials are widely applied in some engineering fields.To exactly describe the viscoelastic features, some scholars have adopted the fractional derivative theory to model the structures.For example, Chen et al. [30] established the fractional dynamic model of an axially moving string and analyzed the transient responses using the Galerkin and numerical methods.Yang et al. [31,32] employed the multiscale method to investigate the nonlinear free and parametric vibrations of an axially moving viscoelastic string with a fractional order damping.
The above-mentioned studies are concentrated on axially moving strings based on Mote's model.Besides, the nonlinear integro-partial-differential equation called Kirchhoff 's model is also used for describing the transverse motion of axially moving strings [33][34][35].Nevertheless, the applications of Kirchhoff 's model in axially moving structures obeying the fractional differentiation law are rather limited, and the difference of fractional parametric resonances between the Mote's model and the Kirchhoff 's model remains unclear.Therefore, the present paper further explores the nonlinear parametric vibration of an axially moving string with the fractional viscoelastic damping based on the literature [32].

Equations of Motion
Consider a uniform, axially moving viscoelastic string made by rubber-like materials, with linear density , the length , and the length cross-sectional area , traveling at timedependent axially speed ( t) between two fixed supports at both ends, shown in Figure 1.The symbol û( x, t) represents the transverse displacement, x is the coordinate along the axial direction, and t is the time.On the basis of the Euler beam theory and the Newton second law [34], the Kirchhoff 's mathematical model governing transverse motion can be obtained as where  0 represents the initial tension and σ( x, t) denotes the axial disturbed stress.According to the Lagrangian strain [35], the disturbed strain accounted for the geometric nonlinearity can be expressed as In the work, the viscoelastic damping of the string is supposed as obeying the fractional derivative Kelvin-Voigt model [32].Thus, the stress-strain relation can be given as where  0 and  1 represent the stiffness constant and viscosity coefficient of the string, respectively, and  is a constant which is used to describe the viscosity characteristic.Considering the fractional derivative operator defined by the Riemann-Liouville form [36,37], we write the expression   / t + ( ) in (3) as follows: where Γ is the Gamma function.The fractional derivative operator is applied to investigate an intermediate viscoelastic characteristic between the elastic string ( = 0) and Kelvin-Voigt viscoelastic string ( = 1).For calculating the fractional differentiation operator, we introduce the following property as where  denotes complex number √ −1.Substituting ( 2), (3), and ( 4) into ( 1), the mathematical model of transverse motion for the axially moving viscoelastic string can be obtained as For the axially moving string with two fixed supports, the boundary conditions are satisfied by Introduce the dimensionless variables where the book-keeping device  is a small parameter, indicating that both the stiffness constant and the viscosity coefficient are very small.Substituting ( 8) into ( 6) and ( 7), one can obtain the following dimensionless dynamic equation, and boundary conditions, In order to study parametric vibrations, we suppose that the axial moving speed is a small harmonically varying about the constant mean speed  0 : where  1 and Ω are the fluctuation amplitude and frequency, respectively.Substituting ( 11) into (9) yields

Analytical Treatments
Equation ( 12) is a fractional, nonlinear integro-partialdifferential equation; the direct multiscale method is applied to solve this equation for finding analytical solutions.Considering a first-order uniform approximation, one gives where  0 and  1 are the displacement functions at orders 1 and ,  0 =  and  1 =  are the fast and slow time scales, respectively.Substituting (13) into (12) and then equalizing the coefficients  0 and  lead to At  0 , the solution can be written as follows: where   is a complex function which is only related to the slow time scale  1 , the overbar represents the complex conjugate, and   and   () are the th natural frequency and mode function of the corresponding linear system, which is given by the following [1]: It is widely known that the principal parametric resonance may occur if Ω approaches 2  .Thus, the nearness relation can be supposed as follows: in which  is a detuning parameter.Substituting ( 16), (17), and ( 18) into ( 15), one yields the following solvability condition [38,39]: where where  1 and  2 are the constant coefficients depending on the fraction order, viscosity, stiffness constant, axial moving speed, and the linear frequency.Separating these coefficients into real and image parts, one yields Express the solution to (19) in the polar form: in which real functions   and   are the amplitude and phase of the th resonance, respectively.Substituting of ( 21) and ( 22) in (19) and then separating real and imaginary parts, we have where For steady-state responses, the amplitude   and the new phase   should be constants; thus,  Eliminating phase   from (25), we can yield the amplitude of steady-state responses: where the positive and the minus signs represent the bigger and the smaller amplitudes, respectively.From [27,28], it is proved that the bigger (smaller) response amplitude is stable (unstable).In addition, (26) is the analytical expression for steady-state responses of the string in principle parametric resonances, where the existing condition of the steady-state responses should be satisfied by

Numerical Results and Discussion
In the section, based on ( 26) and ( 27), the numerical results are presented to discuss the influences of physical parameters like the fractional order on the amplitudes and the existing conditions of steady-state responses in the principle parametric resonances of the first two modes.
Here we choose the parameters as  0 = 0.2,  1 = 0.5,  0 = 0.5,  = 1, and  = 0.5; the differences of steady-state responses for the string based on Kirchhoff 's and Mote's fractional model are shown in Figures 2(a) and 2(b).It can be seen that the response amplitudes and the upper boundaries of exiting conditions predicted by Kirchhoff 's model are much larger than those predicted by Mote's model.This implies that the application of the Kirchhoff 's fractional model in axially moving strings could yield greater security prediction in the engineering.Therefore, we study the parametric resonances of the fractional damping string based on Kirchhoff's model in the next section.
The influence of the fractional order on steady-state responses for Kirchhoff 's model is illustrated in Figures 3(a) and 3(b).It is revealed that the response amplitudes decrease with the fractional order increases.Moreover, decreasing the fractional order leads to a larger upper boundary of the existing conditions, but the lower boundary is independent of the fractional order.It is indicated that the larger fractional order would lead to the more evident influence of the damping and the decreasing in existing intervals of steady-state responses.Besides, the response amplitude in the first parametric resonance is much larger than that in the second one.
The changes of steady-state responses with the stiffness constant and the viscosity coefficient of the string for the Kirchhoff 's fractional model are plotted in Figures 4 and  5, respectively.It can be demonstrated that the response amplitudes decrease dramatically with the increasing of the stiffness constant and the viscosity coefficient.Moreover, the larger the stiffness constant, the larger the existing intervals of steady-state responses.However, a smaller viscosity coefficient can increase the existing intervals, which is in accordance with in the previous viscoelastic string.
Figure 6 depicts the effect of the axial-speed fluctuation amplitude on steady-state responses for Kirchhoff 's fractional string.It is observed that an increase in the axial-speed fluctuation amplitude causes the higher in both the response amplitudes and the existing intervals of steady-state responses.Therefore, we can conclude that the fractional damping string with a smaller axial-speed fluctuation amplitude has a smaller steady-state response.

Conclusions
In this study, on the basis of the fractional definition and calculus law of the Riemann-Liouville form, the nonlinear parametric vibration of an axially moving string with the fractional damping is investigated based on Kirchhoff 's   (2) In Kirchhoff 's model, the fractional derivative has a significant effect on steady-state responses; the smaller fractional order can cause more rapid increasing in the response amplitudes and existing intervals.(3) In Kirchhoff 's fractional model, the response amplitudes decrease with the increasing of the stiffness constant and viscosity coefficient, while they increase as the axial-speed fluctuation amplitude increases.In addition, the larger the stiffness constant, the higher the axial-speed fluctuation amplitude; and the smaller viscosity coefficient tends to enlarge existing intervals of steady-state responses.

Figure 1 :
Figure 1: The physical model of the axially moving strings made by rubber-like materials.
The second-order mode

Figure 2 :
Figure 2: Comparison the steady-state responses between Kirchhoff 's and Mote's fractional model.
The second-order mode

Figure 3 :
Figure 3: Influence of the fractional order on the steady-state responses based on Kirchhoff 's model.

Figure 4 :( 1 )
Figure 4: Influence of the stiffness constant on the steady-state responses based on Kirchhoff 's fractional model.
The second-order mode

Figure 5 :
Figure 5: Effect of the viscosity coefficient on the steady-state responses based on Kirchhoff 's fractional model.

Figure 6 :
Figure 6: Effect of the axial-speed fluctuation amplitude on the steady-state responses based on Kirchhoff 's fractional model.