Abstract
This study presented the transverse vibration of an axially moving beam with an intermediate nonlinear viscoelastic foundation. Hamilton’s principle was used to derive the nonlinear equations of motion. The finite difference and state-space methods transform the partial differential equations into a system of coupled first-order regular differential equations. The numerical modeling procedures are utilized for evaluating the effects of parameters, such as axial translation velocity, flexure rigidities of the beam, damping, and stiffness of the support on the transverse response amplitude and frequencies. It is observed that the dimensionless fundamental frequency and magnitude of axial speed had an inverse correlation. Furthermore, increasing the flexure rigidity of the beam reduced the transverse displacement, but at the same instant, fundamental frequency rises. Vibration amplitude is found to be significantly reduced with higher damping of support. It is also observed that an increase in the foundation damping leads to lower fundamental frequencies, whereas increasing the foundation stiffness results in higher frequencies.
1. Introduction
Axially moving beams have numerous applications in various branches of engineering, such as manufacturing, civil, mechanical, chemical, and aerospace engineering [1–4]. Typical applications explored in previous studies include conveyor belts, medical nanorobots, power saw, band saw, and magnetic tape [5–9]. At the nano/micron level, the motion of the substrate is significant, and small fluctuations can result in a low-quality product [10–13]. Transverse vibrations triggered by external disturbances during a particular process may be lowered by providing end supports beneath the moving beam in the form of rotating rolls [14]. Though end supports are beneficial, they are not sufficient to stabilize the motion of moving substrate. End support controller was used to stabilize the transverse vibration of an axially moving Kirchhoff beam [15]. Additional intermediate support in the moving span is often required [16]. Intermediate support between the rolls offers superior web processing quality in comparison to the traditional end support. Vibration of an axially moving beam with an intermediate foundation has become a burning research topic since the last decade. Several studies have attempted to analyze axially moving beams with elastic foundations. Dynamic analysis of coupled moving beams through a Winkler elastic support was investigated by Gaith and Müftü [17]. Yang et al. [18] studied the free vibration of axially moving elastic beams resting on an elastic foundation. In their study, they utilized the Galerkin method to discretize the equation of motion. Fundamental frequencies of transverse fluctuations in a moving beam were studied by utilizing different numerical methods [19]. Liu et al. [20] explored the dynamics and instability of the Euler–Bernoulli beam as a function of several variables such as axially moving speed, position, and weight of the lumped mass. The study results revealed that any increase in lumped mass is inversely associated with resulting systems’ fundamental frequencies. Dynamics of axially moving beam attached with energy sink and a piezoelectric device was studied in [21, 22]. A significant effect of the attachment on the response was observed. Norouzi and Younesian [23] investigated the chaotic behavior of an axially moving beam supported by a stiff foundation. Cosine-cosine function was used to obtain the approximate solution of the equation of motion. Chaos indication was obtained by setting Lyapunov exponent as a criterion. Mohammadzadeh and Mosayebi [24] investigated the vibration in rail supports and joints by considering it as an axially moving beam resting on an elastic foundation. The effect of rail support and joints in the locality of a bridge in railway track was investigated.
Kural and Ozkaya [25] investigated the performance of a microbeam fluid-carrying media by utilizing the modified couple stress theory. In another study, Yan et al. [26] employed the theoretical model to examine the stability and dynamics of a thin cantilever beam. The beam was attached with a moving station submerged in a fluid. Calim [27] explored the patterns of free and forced fluctuation in a functionally graded beam assembly supported with an intermediate elastic foundation. Tang et al. [28] explored the transverse vibration of an axially accelerating beam having time-dependent axial tension. Analysis revealed that variation in axial speed and axial tension has a direct influence on the resonance. Mirzabeigy and Madoliat [29] studied the large amplitude-free vibration of beams resting on the supporting foundation. The Winkler model and Euler–Bernoulli method were used to study the induced vibrations in elastic foundation and beam, respectively. Ghannadiasl [30] proposed the application of dynamic Green function (DGF) to investigate the dynamic response of railways under varying loading conditions. The load was varying in nature with different speeds and accelerations. Variation in moving load, as well as the elastic parameters of the foundation, was investigated. It was shown by the modeling results that the maximum deflection depends on the increasing or decreasing acceleration of the moving load.
Similarly, Zhang et al. [31] utilized the complex modal analysis to investigate the dynamic response of a moving beam resting on intermediate support. Ding et al. [32] investigated the axially moving beam traveling with super-critical speed. The impact of rotary inertia and shear deformation on the transverse fluctuations was investigated. Zhao et al. [33] analytically investigated an axially moving microbeam. Coupled thermoelastic vibration and heat transfer process was presented. Mohamed et al. [34] presented a novel numerical procedure to forecast the forced steady-state and nonlinear free response of curved beam under clamped-clamped boundary conditions. An and Su [35] presented the generalized integral transform technique as a numerical approach for the dynamic analysis of Timoshenko beam. Esen [36] considered functionally graded Timoshenko beam to investigate the dynamic response using the modified finite element method (FEM). The beam under consideration was resting on elastic support and was exposed to a mass moving with a variable speed. For the first time, Ding et al. [37] proposed a generalized boundary conditions approach for axially moving beams supported by vertical and torsional springs at both ends. Shao et al. [38] investigated the nonlinear dynamic behavior of a moving membrane with fluctuating speed. It was concluded that the nonlinear vibration characteristics of a membrane are sensitive to the initial motion conditions.
A thorough review of the literature reveals a significant lack of comprehensive investigations on the problem of an axially moving beam with nonlinear viscoelastic intermediate support. To the best of author’s knowledge, problem of transverse vibration in an axially moving beam is tackled by providing supports at the ends of the moving beam. In authors’ opinion, end supports are likely to result in lower quality process. Therefore, intermediate viscoelastic supports are presented here as a possible solution to mitigate the vibration in moving beams. Also, in view of the nature of involved equations of motion, authors also feel the need to investigate the application of the finite difference method (FDM) on the nonlinear dynamics of beams under consideration. In this study, transverse vibrations in axially traveling beams with a roll-to-roll configuration supported by the viscoelastic intermediate foundation were studied. Nonlinear equations of motion for transverse vibrations were derived using Hamilton’s principle. A model based on the finite difference method coupled with the state-space approach was presented. The effect of parameters, such as axial translation velocity, flexure rigidities of the beam, damping, and stiffness of the support on the transverse response amplitude and frequencies of the system in the subcritical region, were investigated.
2. Problem Formulation
An axially moving beam on a viscoelastic foundation for analyzing the vibrations characteristics is presented in Figure 1. This model mainly includes stiffness, damping, length, and end supports. The beam has flexural rigidity EI and axial tension T and moves with an axial velocity (). Transverse displacement of the beam is (x, t).

Free body diagram of an infinitesimal element of the considered beam is shown in Figure 2. Here, S is the shear force and Mb is the bending moment. Viscous damping force, spring force, and inertial force are represented by Fd, Fs, and FI, respectively. First, the total kinetic energy and potential energy of the moving beam were calculated by considering as the density and A as the cross-sectional area of the beam. Hamilton’s principle was then applied to derive the equation of motion. Transverse vibration of the beam is described bywhere E is the modulus off elasticity, is the density of the beam material, A is the cross-sectional area, and I is the moment inertia of the beam. Axial speed and axial tension in the beam are represented by and T, respectively. The linear stiffness, nonlinear stiffness, and damping of the foundation are represented by , , and , respectively. Introducing the following dimensionless variables transforms the equation of motion into a normalized equation:where represents the dimensionless transport speed (a ratio between physical and wave velocity in the beam), while parameter is the dimensionless flexural rigidity.

2.1. Finite Difference Formulation
Transverse vibrations of the moving beam on a nonlinear viscoelastic foundation are numerically solved, and essential mathematical formulation is presented in equation (2). An analysis approach based on the finite difference method is used to obtain second-order ordinary differential equations (ODEs) by transforming the fourth-order HPDE. Discretization of the equation is then achieved by the central difference scheme given by the following equations:
The following second-order ordinary differential equations are produced:where the total number of spatial points is represented by n and step size is represented by . The subsequent matrix form of equation (6) is given bywhere M, G, , and are matrices, which have the following structure:
System of second-order ordinary differential equations was then transformed into a system of first-order ordinary differential equations by applying state-space approach:where
The initial conditions used in the system of ODE’s (8) arewhere is the initial transverse deflection and represents the initial transverse speed.
3. Numerical Illustrations and Discussion
Numerical simulations were performed to investigate the effect of different factors such as flexural rigidity of the beam, axial speed of the beam, damping of the foundation, and stiffness of the foundation on the vibration characteristics of the axially traveling beam. Values of the parameters used in the simulations are given in Table 1. Simulations were performed for the initial transverse displacement of “” and initial transverse speed condition of “”. An in-house matlab code was used for the simulations.
First of all, the solution method used in the current analysis (Method 1) was verified by comparing it with the method (Method 2) presented by An and Su [35]. Comparison of the midspan transverse displacement for the case of an axially moving beam having nondimensional speed of “” and nondimensional flexure rigidity of “” is shown in Figure 3.

Response obtained from both solution methods was found to be in good agreement. It was also observed that the solution method used in this study is an efficient approach having long-term stability and excellent convergence.
Figures 4 and 5 show the comparison of transverse vibration as a function of axial speed in an axially moving beam. Transverse displacement at different locations () of the beam at different axial speeds, , 0.3, and 0.5, (equals to 3 m/s, 9 m/s, and 15 m/s, respectively) is presented in Figure 4. Figure 5 plots the trend of the fundamental frequency at different axial speeds. It may be noted from Figure 5 that an increased value of axial speed is accompanied by an exponential decrease in dimensionless fundamental frequency. For the stationary beam, it has a positive value of 7.3, which ultimately reaches zero for an axial speed approaching a value of 1.18 (equals to 35 m/s). The critical speed of the beam can also be interpreted as the speed at which the frequency of the vibration becomes zero; therefore, the critical speed of an axially moving beam under the abovementioned conditions is 35 m/s. Although the fundamental frequency of the vibration is experiencing a downward trend, the amplitude of vibration is unaffected, as evident from Figure 4. Therefore, it is concluded that the critical speed at flexure rigidity of “0.2” is 1.18. Similarly, Figures 6 and 7 show the impact of flexure rigidity on the transverse fluctuations of moving beams resting on the viscoelastic foundation.




Figure 6 shows the effect of flexural rigidity, , on the midspan transverse displacement of axially moving beam. Based on Figure 6, increasing the flexure rigidity of the beam reduces the transverse displacement (24% decrease in displacement magnitude by increasing from 0.2 to 0.3, and 47% decrease in magnitude by increasing from 0.3 to 0.5) but at the same time will increase the nondimensional fundamental frequency from 7.3 to 9.8, as shown in Figure 7. Flexure rigidity also increases the critical speed of the axially moving beam (from 1.18 to 1.65). The higher value of flexure rigidity will somehow increase the internal damping (in terms of lower vibration amplitude) but at the cost of higher fundamental frequency.
Figures 8 and 9 represent the effect of damping of viscoelastic support on the vibration of the beam. Transverse displacement of beam at different locations () for different values of the damping factor is given in Figure 9. A sharp decline in the response is observed throughout the beam length with an increase in the damping of the viscoelastic support.


Transverse displacement of midpoint () of the beam is plotted for different axial speeds () at increasing values of damping factor () (Figure 9). As shown in Figure 9, with an increase in damping of the support, there is a quite visible reduction in the vibration amplitude. Interestingly, the vibration amplitude reduction is more visible when beam is moving with higher axial speed. A rapid reduction in the vibration magnitude for the case of is observed when the axial speed is
The effect of the support damping on the nondimensional fundamental frequency of the beam is presented in Figure 10. Nondimensional frequency and the critical speed of the beam decrease with an increase in the damping. Nondimensional frequency decreases from 7.3 to 6.4, and critical speed decreases from 1.18 to 1.15 when the damping is increased from 0.05 to 3.5.

Figure 11 presents the variation of nondimensional fundamental frequency as a function of nondimensional axial speed at various magnitudes of support stiffness. An opposite effect to damping on the nondimensional frequency and critical speed is observed in the case of stiffness of the foundation. The frequency of the beam increases from 7.3 to 8.8, and the critical speed increases from 1.18 to 1.43 when the linear stiffness of the beam increases from 0.6 to 1.5 (corresponding to 800 N/m and 6000 N/m, respectively).

A case of a forced vibration was also considered with external excitation force . The magnitude of the excitation force is “F,” and the frequency of the excitation is Ω. Frequency and displacement of the axially moving beam were measured under various excitation frequencies and are presented in Figure 12.

(a)

(b)

(c)

(d)
Figures 12(a)–12(c) show the frequency domain response of the axially moving beam. Starting with lower excitation frequency (Ω = 1.0), frequency domain response, Figure 12(a) indicates that the fundamental frequency of the beam is higher than the excitation frequency (). As the excitation frequency is increased, fundamental frequency of the beam and the external excitation frequency match at “Ω = 7.4,” and as a result, the response increases exponentially (Figure 12(d)).
4. Conclusions
In this study, transverse vibrations in axially traveling beam, supported by the viscoelastic intermediate foundation, was studied. Nonlinear equations of motion were derived using Hamilton’s principle. The finite difference method coupled with the state-space approach was presented as an efficient solution technique. The effect of parameters, such as axial translation velocity, flexure rigidities of the beam, damping, and stiffness of the support on the transverse response amplitude and frequencies of the system in subcritical region, was investigated. The major findings of this study may be summarized as follows:(i)In the absence of effective damping, the dimensionless fundamental frequency was observed to decrease with the increasing value of axial speed. Having the positive value of 7.3 for the case of stationary beam ultimately reaches to zero for an axial speed approaching a value of 1.18. Although, the fundamental frequency of the vibration is decreasing, the amplitude of vibration is unaffected.(ii)Higher value of flexure rigidity increases the internal damping (in terms of lower vibration amplitude) but at the cost of higher fundamental frequency. A decrease of 24% in the peak magnitude of transverse displacement is observed with increasing from 0.2 to 0.3 and 47% decrease in the peak magnitude by increasing from 0.3 to 0.5. In contrast, increasing the flexural rigidity, nondimensional frequency, and critical speed of the beam increases from 7.3 to 9.8 and 1.18 to 1.65, respectively.(iii)With an increase in the damping of support, there is a visible reduction in the vibration amplitude. Interestingly, the vibration amplitude reduction is found to be more significant when the beam is moving with a higher axial speed. A significant reduction in the peak vibration magnitude is observed for the case of beam having and axial speed of Nondimensional frequency decreases from 7.3 to 6.4, and critical speed decreases from 1.18 to 1.15 when the damping is increased from 0.05 to 3.5.(iv)Higher value of the foundation stiffness results in the increase of both nondimensional frequencies as well as the critical speed of the beam. Critical speed increases from 1.18 to 1.43, whereas the nondimensional frequency of the beam increases from 7.3 to 8.8 when the linear stiffness of the beam is increased from 0.6 to 1.5.
Abbreviations
A: | Cross-sectional area of the beam |
: | Dimensionless axial speed |
E: | Modulus of elasticity of beam |
: | Initial transverse speed |
T: | Tension in beam |
: | Density of the beam |
: | Axial speed of beam |
L: | Beam length |
: | Foundation damping |
: | Transverse displacement |
: | Linear stiffness of foundation |
: | Dimensionless flexural rigidity |
: | Nonlinear stiffness of foundation |
: | Initial transverse deflection. |
Data Availability
The underlying data supporting the results of this work are found within the article.
Conflicts of Interest
The authors have no conflicts of interest to disclose.
Acknowledgments
The authors appreciate and acknowledge the research support provided by the Deanship of Scientific Research, Imam Abdulrahman Bin Faisal University (IAU), for conducting this study.