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

Authors: Ying Li; Ye Tang
Journal: Mathematical Problems in Engineering
 (2017)

---

## Abstract

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.


---

## Body

## 1. 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–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–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].


## 2. Equations of Motion

Consider a uniform, axially moving viscoelastic string made by rubber-like materials, with linear densityρ, the length L, and the length cross-sectional area A, traveling at time-dependent axially speed Vt^ between two fixed supports at both ends, shown in Figure 1. The symbol u^(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 (1)ρA∂2u^∂t^2+2V∂2u^∂x^∂t^+dVdt^∂u^∂x^+V2∂2u^∂x^2-P0∂2u^∂x^2-∂2u^∂x^21L∫0LAσ^x^,t^dx^=0,where P0 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 (2)εLx^,t^=12∂u^∂x^2.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 (3)σ^x^,t^=E0εLx^,t^+E1∂α∂t^αRL+εLx^,t^,where E0 and E1 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^αRL+ in (3) as follows:(4)∂α∂t^αRL+εL=1Γ1-αddt^∫0t^t^-τ-αεLτdτ,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(5)dαdt^αRL+eiωt^=iωαeiωt^+sin⁡παπ∫0∞uαu+iωe-ut^du,where i denotes complex number -1.Figure 1
The physical model of the axially moving strings made by rubber-like materials.Substituting (2), (3), and (4) into (1), the mathematical model of transverse motion for the axially moving viscoelastic string can be obtained as (6)ρA∂2u^∂t^2+2V∂2u^∂x^∂t^+dVdT∂u^∂x^+V2∂2u^∂x^2-P0∂2u^∂x^2=E0A2L∂2u^∂x^2∫0L∂u^∂x^2dx^+E1A2L∂2u^∂x^2∫0L∂α∂tα∂u^∂x^2dx^.For the axially moving string with two fixed supports, the boundary conditions are satisfied by(7)u0,t^=uL,t^=0.Introduce the dimensionless variables(8)u=u^L,x=x^L,t=t^LP0ρA,γ=VρAP0,εe0=E0AP0,εη=E1APα/2-1ρAL2α,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,(9)∂2u∂t2+2γ∂2u∂x∂t+dγdt∂u∂x+γ2-1∂2u∂x2=12εe0∂2u∂x2∫01∂u∂x2dx+12εη∂2u∂x2∫01∂α∂tα∂u∂x2dx,and boundary conditions,(10)u0,t=u1,t=0.In order to study parametric vibrations, we suppose that the axial moving speed is a small harmonically varying about the constant mean speed γ0:(11)γ=γ0+εγ1sin⁡Ωt,where γ1 and Ω are the fluctuation amplitude and frequency, respectively. Substituting (11) into (9) yields(12)∂2u∂t2+2γ0∂2u∂x∂t+dγdt∂u∂x+γ02-1∂2u∂x2=ε12e0∂2u∂x2∫01∂u∂x2dx+12η∂2u∂x2∫01∂α∂tα∂u∂x2dx-γ1Ωcos⁡Ωt∂u∂x-2γ1sin⁡Ωt∂2u∂x∂t-2γ0γ1sin⁡Ωt∂2u∂x2.


## 3. Analytical Treatments

Equation (12) is a fractional, nonlinear integro-partial-differential equation; the direct multiscale method is applied to solve this equation for finding analytical solutions. Considering a first-order uniform approximation, one gives (13)ux,t,ε=u0x,T0,T1+εu1x,T0,T1+⋯,where u0 and u1 are the displacement functions at orders 1 and ε, T0=t and T1=εt are the fast and slow time scales, respectively. Substituting (13) into (12) and then equalizing the coefficients ε0 and ε lead to (14)ε0:∂2u0∂T02+2γ0∂2u0∂T0∂x+γ02-1∂2u0∂x2=0,(15)ε1:∂2u1∂T02+2γ0∂2u1∂T0∂x+γ02-1∂2u1∂x2=-2∂2u0∂T0∂T1-2γ0∂2u0∂x∂T1+12e0∂2u0∂x2∫01∂u∂x2dx+12η∂2u0∂x2∫01∂∂T0RL+α∂u∂x2dx-γ1Ωcos⁡ΩT0∂u0∂x-2γ1sin⁡ΩT0∂2u0∂x∂T0-2γ0γ1sin⁡ΩT0∂2u0∂x2.At ε0, the solution can be written as follows:(16)u0x,T0,T1=AnT1ϕnxeiωnT0+A-nT1ϕ-nxe-iωnT0,where An is a complex function which is only related to the slow time scale T1, the overbar represents the complex conjugate, and ωn and ϕn(x) are the nth natural frequency and mode function of the corresponding linear system, which is given by the following [1]: (17)ωn=nπ1-γ02,ϕnx=2sin⁡nπxeinπγ0x.It is widely known that the principal parametric resonance may occur if Ω approaches 2ωn. Thus, the nearness relation can be supposed as follows:(18)Ω=2ωn+εσ,in which σ is a detuning parameter.Substituting (16), (17), and (18) into (15), one yields the following solvability condition [38, 39]:(19)dAndT1+χ1An2A-n+χ2A-neiσT1=0,where(20)χ1=-2e0∫01ϕ-n′′ϕ-ndx∫01ϕn′ϕ-n′dx+e0+η2iωnα∫01ϕ-n′′ϕ-ndx∫01ϕn′2dx4iωn∫01ϕnϕ-ndx+γ0∫01ϕn′ϕ-ndx,χ2=-ωn-Ω/2∫01ϕ-n′ϕ-ndx+2γ0∫01ϕ-n′′ϕ-ndx2iωn∫01ϕnϕ-ndx+γ0∫01ϕn′ϕ-ndxγ1,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(21)χ1=Reχ1+iImχ1,χ2=Reχ2+iImχ2.Express the solution to (19) in the polar form:(22)An=12anT1eiϕnT1,in which real functions an and ϕn are the amplitude and phase of the nth resonance, respectively. Substituting of (21) and (22) in (19) and then separating real and imaginary parts, we have(23)dandT1=anImχ2sin⁡θn-Reχ2cos⁡θn-14Reχ1an3,dθndT1=σ+2Reχ2sin⁡θn-Imχ2cos⁡θn+12Imχ1an2,where(24)θn=σT1-2ϕn.For steady-state responses, the amplitude an and the new phase θn should be constants; thus, (25)anImχ2sin⁡θn-Reχ2cos⁡θn-14Reχ1an3=0,σ+2Reχ2sin⁡θn-Imχ2cos⁡θn+12Imχ1an2=0.Eliminating phase θn from (25), we can yield the amplitude of steady-state responses: (26)an=-2σImχ1±24χ12χ22-σ2Reχ12χ1,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(27)∓2χ1≤σ≤2χ1χ2Reχ2.


## 4. 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, e0=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.Figure 2
Comparison the steady-state responses between Kirchhoff’s and Mote’s fractional model.
(a)
The first-order mode
(b)
The second-order modeThe influence of the fractional order on steady-state responses for Kirchhoff’s model is illustrated in Figures3(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.Figure 3
Influence of the fractional order on the steady-state responses based on Kirchhoff’s model.
(a)
The first-order mode
(b)
The second-order modeThe 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 Figures4 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 4
Influence of the stiffness constant on the steady-state responses based on Kirchhoff’s fractional model.
(a)
The first-order mode
(b)
The second-order modeFigure 5
Effect of the viscosity coefficient on the steady-state responses based on Kirchhoff’s fractional model.
(a)
The first-order mode
(b)
The second-order modeFigure6 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.Figure 6
Effect of the axial-speed fluctuation amplitude on the steady-state responses based on Kirchhoff’s fractional model.
(a)
The first-order mode
(b)
The second-order mode


## 5. 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 model. Consider that the viscoelastic behavior of the string is described by the fractional Kelvin-Voigt constitutive relation; a new nonlinear fractional mathematical model governing transverse vibration of the string is derived from the Euler beam theory and Newton’s second law. The direct multiscale method is applied to investigate the steady-state responses of the fractional damping string undergoing the principle parametric resonances. The investigation yields the following conclusions:(1)
The amplitudes and existing intervals of steady-state responses predicted by Kirchhoff’s fractional model are much larger than those predicted by Mote’s fractional model.(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.



---

*Source: 1393954-2017-11-16.xml*