Next Article in Journal
Identification Approach for Nonlinear MIMO Dynamics of Closed-Loop Active Magnetic Bearing System
Previous Article in Journal
An Overview of Non-Isolated Hybrid Switched-Capacitor Step-Up DC–DC Converters
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Buckling Induced Strongly Nonlinear Vibration of Supercritical Axially Moving Beam

1
Tianjin Key Laboratory for Civil Aircraft Airworthiness and Maintenance, Civil Aviation University of China, Tianjin 300300, China
2
School of Mechanical Engineering and Automation, Northeastern University, Shenyang 110819, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(17), 8555; https://doi.org/10.3390/app12178555
Submission received: 23 July 2022 / Revised: 21 August 2022 / Accepted: 24 August 2022 / Published: 26 August 2022
(This article belongs to the Section Acoustics and Vibrations)

Abstract

:
In this paper, the postbuckling problem and the free linear and nonlinear vibration surrounding postbuckling configuration are accurately analyzed. The placement of translating beams that undergo significant overall buckling deformations in this high-speed region, leading to strongly nonlinear vibration, is studied. Using the Galerkin method, we found that the partial differential equation is reduced to an ordinary differential equation with cubic and quadratic terms, which differs from the previous literature. The supercritical natural frequency is analytically obtained, the strongly nonlinear vibration problem is solved using the homotopy analysis method (HAM), and the solution is examined using numerical results. Results show that the obtained analytical solution is valid for both weakly and strongly nonlinear cases.

1. Introduction

Axially moving systems, such as threads, textiles, wires, magnetic tapes, belts and cables, are present in a wide class of engineering structures. The vibration of such devices has been studied extensively [1,2,3,4,5,6]. In general, one-dimensional string or beam theory can be used to model axially moving elements, depending on flexural stiffness [7,8,9,10,11,12,13,14].
It is known that the linear frequency of a moving beam decreases as its axial speed increases. At certain critical speeds, the linear frequency vanishes and buckling occurs [3,7]. In previous studies, the linear and nonlinear vibration of axially moving beams were studied in the subcritical regime, and referred to as subcritical axially moving beams. Herein, we refer to beams traveling in the supercritical regime as supercritical axially moving beams. Öz [3] studied the linear vibration of an axially moving beam with time-dependent velocity below critical value. The multi-scales method was used to investigate the resulting parametric resonance. Li and Lu [15,16] described three-dimensional parametric vibrations of a simply-supported pipe conveying pulsating fluid, and found that the jump phenomenon for the parametric resonance in the post-buckling state could be clearly observed. Chen and Yang [17,18] derived a new governing equation for the traveling beam using geometric nonlinearity, which was used to account for finite stretching of the beam; the nonlinear vibration problem was studied using the perturbation method. Ghayesh et al. [19,20,21] contributed to the field of forced dynamics of a subcritical axially moving viscoelastic beam and investigated the beam’s bifurcation and chaotic motions.
The literature regarding supercritical axially moving beams is rather limited. Hwang and Perkins [22,23] studied the influence of the initial curvature caused by the support wheel and pulley on the bifurcation and stability in the supercritical regime. Ghayesh [24,25] included the nonlinear stretching effects for a moving beam and analyzed both the subcritical and supercritical speed range. Ding et al. [26,27,28,29] utilized the Galerkin method and other numerical methods to calculate the supercritical natural frequency and equilibrium configurations for different boundary conditions. Zhang et al. [30] studied the stability of an axially moving beam subjected to harmonic excitation. They employed the multi-scales method to calculate the steady-state response amplitude for the first single mode; the strongly nonlinear state was shown in some situations.
In this paper, we focus on the dynamics of supercritical axially moving beams. In this high-speed region, the translating beam undergoes large overall buckling deformations leading to strongly nonlinear vibration. This situation is common; however, the solution is not available in the literature. Using Galerkin truncation, we found that the partial differential equation can be simplified and that the dynamic process can be described by an ordinary differential equation, which differs from the previous literature. Without using transformation, we directly used the homotopy analysis method to address the strongly nonlinear problem.

2. Governing Equations and Solutions

The equation of motion of an axially moving beam takes the form [3,18]
ρ A ( 2 U T 2 + 2 V 2 U X T + V 2 2 U X 2 ) P 0 2 U X 2 + E I 4 U X 4 = E A L 0 L ( 1 2 ( U X ) 2 ) d X 2 U X 2
where ρ is linear density, U is transverse displacement, A is cross-sectional area, P0 is initial tension, V is transport speed, and EI is flexural stiffness. The beam travels between two simple supports, which are separated by L.
The case V = 0 corresponds to a static beam with deformation, and the governing equation is reduced to
ρ A 2 U T 2 P 0 2 U X 2 + E I 4 U X 4 = E A L 0 L ( 1 2 ( U X ) 2 ) d X 2 U X 2
Introduce the following non-dimensional quantities
x = X L , u = U L , t = T L P 0 ρ A , γ = V ρ A P 0 , v f 2 = E I P 0 L 2 , k f 2 = E A P 0
and Equation (4) is obtained by substituting Equation (3) into Equation (1).
u ¨ + 2 γ u ˙ + ( γ 2 1 ) u + v f 2 u ( 4 ) = 1 2 k f 2 u 0 1 u 2 d x
The pinned–pinned supported and the fixed–fixed boundary conditions are
u = 0   and   u = 0   at   x = 0 ,   1  
u = 0   and   u = 0   at   x = 0 ,   1  
Equation (7) is the supercritical equilibrium equation induced from Equation (4) and ψ is the supercritical configuration.
( γ 2 1 ) ψ + v f 2 ψ ( 4 ) = 1 2 k f 2 ψ 0 1 ψ 2 d x
The integral of Equation (7) is a constant for a given ψ. Let
Π = 1 2 k f 2 0 1 ψ 2 d x
Next, Equation (7) can be rewritten as
ψ ( 4 ) + λ 2 ψ = 0
where
λ 2 = ( γ 2 1 ) Π v f 2
The characteristic equation is shown in Equation (11).
sin λ = 0   or   λ = m π ,   m = 1 ,   2 ,   3    
The equilibrium configurations for the pinned–pinned supported case is
ψ ( x ) = c sin ( λ x )
where
c = ± 2 λ k f ( γ 2 1 ) v f 2 λ 2
Note that the equilibrium configurations (12) were also obtained by Wickert [7] using a different method.
Figure 1 illustrates subcritical and supercritical equilibrium configurations. The blue dashed line indicates the equilibrium position of a beam in the subcritical regime. If the beam speeds up and exceeds critical speed, it will stand vertically as an arch, and be symmetrical about the x axis, as shown by the solid black curve.
It is noted that the configuration given by (12) is symmetric for the pinned–pinned case. Yang et al. [31] found two special configurations for the fixed–fixed boundary condition, which are shown below.
(1)
Symmetric configuration:
ψ ( x ) = c [ 1 cos ( 2 k π x ) ] ,   λ = 2 π ,   4 π  
(2)
Anti-symmetric configuration:
ψ ( x ) = c [ 1 2 x cos ( λ x ) + 2 λ sin ( λ x ) ] ,   λ = 8.9868 ,   15.4505
In Equations (14) and (15), the coefficient c is defined by Equation (13).
To obtain the natural frequency of a supercritical axially moving beam subjected to buckling, we introduce a small disturbance as follows
u = ψ ( x ) + v ( x , t )
where v(x,t) is a small dynamic disturbance around the buckling configuration ψ ( x ) . Substituting Equation (16) into Equation (4) leads to
v ¨ + 2 γ v ˙ + ( γ 2 1 ) v + v f 2 v ( 4 ) = k f 2 ( ψ 0 1 ψ v d x + 1 2 ψ 0 1 v 2 d x + v 0 1 ψ v d x + 1 2 v 0 1 v 2 d x )
Note that the first term, k f 2 ψ 0 1 ψ v d x , on the right side of Equation (17) is linear. Using the configuration relation (7) and dropping the nonlinear terms of Equation (17), we obtain the linear vibration equation
v ¨ + 2 γ v ˙ + v f 2 λ 2 v + v f 2 v ( 4 ) 1 2 k f 2 ψ 0 1 ψ v d x = 0
Equation (19) corresponds to the case V = 0.
v ¨ + v f 2 λ 2 v + v f 2 v ( 4 ) 1 2 k f 2 ψ 0 1 ψ v d x = 0
in which ψ ( x ) is the linear equilibrium configuration defined by Equations (12), (14) and (15), depending on different boundary conditions. Supercritical linear vibration properties can be studied using Equation (18), which significantly differs from previous linear subcritical models [3,5].
We let
v ( x , t ) = j = 1 N ϕ j ( x ) e i ω t
where ϕ j ( x ) is the j-th supercritical linear vibration mode around the buckled shape and ω is the corresponding natural frequency. Substituting Equation (20) into Equation (18) leads to
v f 2 ϕ j ( 4 ) + 2 γ i ω ϕ j + v f 2 λ 2 ϕ j 1 2 k f 2 ψ 0 1 ψ ϕ j d x = ω 2 ϕ j
Multiplying Equation (21) by ϕ i ( x ) and integrating the result over [0, 1] yields
0 1 L ( ϕ j ) ϕ i d x = ω 2
where
L ( ϕ j ) = v f 2 ϕ j ( 4 ) + v f 2 λ 2 ϕ j + 2 γ i ω ϕ j 1 2 k f 2 ψ 0 1 ψ ϕ j d x
which denotes the linear operator. Evaluating Equation (23) yields
v f 2 λ i 4 δ i j v f 2 λ i 4 C i j + i ω i 2 B i j 1 2 k f 2 λ i 4 c 2 = ω i 2
where
B i j = 0 1 ϕ i ϕ j d x ,   C i j = 0 1 ϕ i ϕ j d x
Evaluations of B and C are shown in Appendix A.
If m = 1, the first supercritical natural frequency for a pinned–pinned case is
ω 1 2 1 2 k f 2 c 2 λ 4 = 1 2 k f 2 4 π 4 k f 2 ( ( γ 2 1 ) π 2 v f 2 ) 2 π 4 ( ( γ 2 1 ) π 2 v f 2 )
Equation (26) shows that if V = 0, ω will become a constant.
ω 1 2 2 π 2 v f 2
The first supercritical natural frequency for a fixed–fixed case is
ω 1 2 = v f 2 λ 1 4 v f 2 λ 1 4 C 11 1 2 k f 2 λ 1 4 c 2
Ding and Chen [32] adopted the finite difference method to calculate the natural frequency of an axially moving beam in the supercritical regime. They numerically showed that supercritical natural frequency for a pinned–pinned case is independent of the nonlinear coefficient, kf. The supercritical natural frequency increases regardless of the magnitude of the nonlinear coefficient, kf. In this paper, Equations (25) and (27) are the analytical approximation of a supercritical axially moving beam, in which kf disappears. Moreover, we analytically obtained the supercritical natural frequency for a fixed–fixed case and found that, in contrast, it is highly dependent on kf.
First, to facilitate comparison with classical results, we plot the natural frequencies in the subcritical region in Figure 2 and Figure 3, where ω is the axial speed and γ is the natural frequency. Note that the subcritical natural frequency decreases as the axial speed increases for both the pinned–pinned and fixed–fixed cases. If V = 0, these series of curves in Figure 2 and Figure 3 will reduce to a series of points on the y-axis, which means that the subcritical natural frequency becomes a fixed value.
In Figure 4 and Figure 5, variations in supercritical natural frequency with axial speed increases are depicted by Equations (25) and (27). Note that the supercritical natural frequencies increase as axial speeds increase, which is quite different from a subcritical case. Moreover, a beam with greater flexural stiffness has a larger vibration natural frequency and starts to vibrate at a higher axial speed. This effect of flexural stiffness on the supercritical natural frequency is the same as that of a subcritical case. The basis for this finding is quite reasonable; similar to the one-degree-of-freedom system, stiffness is proportional to natural frequency.
After analyzing nature frequencies, the nonlinear vibration Equation (17) is studied. We used standard Galerkin-type projections to approximate Equation (17). Assume the displacement expansion is
v = n = 1 N ϕ n ( x ) η n ( q )
where ϕ n ( x ) is the eigenfunction of the buckled beam. For a pinned–pinned case, ϕ n ( x ) = sin ( n π x ) . For a fixed–fixed case,
ϕ n ( x ) = cos ( β n x ) cosh ( β n x ) cos ( β n x ) cosh ( β n x ) sin ( β n x ) sinh ( β n x ) [ sin ( β n x ) sinh ( β n x ) ]
where βn is the n-th root of the equation
cos ( β n x ) cosh ( β n x ) = 1
In this section, we only take the pinned–pinned boundary condition as an example. ηn(t) is the generalized coordinates of the discretized system, and N is the number of retained modes. Substituting Equation (28) into Equation (17) yields
η ¨ m + ω m 2 η = k f 2 [ c i , j = 1 N A m i j + i , j = 1 N B m i j ]
where
A m i j = 1 2 n 2 π 2 ( 0 1 ϕ i ϕ m d x ) ( 0 1 ϕ i ϕ m d x ) + n π ( 0 1 ϕ i ϕ m d x ) ( 0 1 ϕ i ϕ m d x )
B m i j = 1 2 ( 0 1 ϕ i ϕ m d x ) ( 0 1 ϕ i ϕ m d x )
The single-mode approximation can be obtained by letting N = 1
η ¨ + ω 1 2 η + k f 2 ( c a 2 η 2 + a 3 η 3 ) = 0
Equation (32) is the key equation in this paper, which contains both a quadratic term and a cubic term at the same time. Evaluating a2 and a3 for N = 1, the results are
a 2 = 1 2 π 2 ( 0 1 sin π x ϕ 1 d x ) ( 0 1 ( ϕ 1 ) 2 ) π ( 0 1 cos π x ϕ 1 d x ) ( 0 1 ϕ 1 ϕ 1 d x ) = 1 8 π 4
a 3 = 1 2 ( 0 1 ϕ 1 ϕ 1 d x ) ( 0 1 ( ϕ 1 ) 2 d x ) = 1 2 ( 0 1 ( ϕ 1 ) 2 d x ) 2 = 1 8 π 4
The nonlinear free vibration problem of Equation (32) will be investigated in the next section. First, the parameters are discussed in detail. Recalling the non-dimensional flexural stiffness and nonlinear coefficient in Equation (3)
v f 2 = E I P 0 L 2 , k f 2 = E A P 0
it follows from Equation (35) that
0 k f 2 = 12 v f 2 ( L h ) 2
In general, note that the non-dimensional flexural stiffness is in the range [0, 1], and for a slender beam is in the range L / h 1 . The nonlinear coefficient, k f 2 , is not a small parameter; therefore, Equation (32) exhibits as strongly nonlinear, and the traditional perturbation method is not valid. In the next section, we employ the homotopy analysis method to solve Equation (32).

3. Principle of Homotopy Analysis

Liao [33] introduced the basic principle of homotopy analysis. Assume that Equation (34) is a nonlinear differential equation
F u ( x , t ) = g ( x , t )
where F is a nonlinear ordinary or partial differential operator, u(x,t) is an unknown function, and g is the source term.
F = L + N
where L and N indicate the linear and nonlinear terms, respectively. Equation (39) is the zero-order deformation equation:
( 1 q ) L [ Φ ( t ; q ) u 0 ] = q [ F Φ ( t ; q ) g ]
The convergence range can be well controlled by , and the constraint conditions for constructing auxiliary linear operator, L, are relatively loose.
When q = 0 and q = 1, it holds that
Φ ( t ; 0 ) = u 0 , Φ ( t ; 1 ) = u ( q , t )
As q changes from 0 to 1, the solution of Equation (39) also changes from u0 to u(x,t). Equations (41) and (42) are the Taylor series expansion of Φ(t;q), the exact solution of Equation (37).
Φ ( t ; q ) = u 0 ( t ; q ) + m = 1 + u m ( t ; q ) q m
u m ( t ; q ) = 1 m ! m Φ ( t ; q ) q m
If we perform differentiation on Equation (41) m times, the corresponding deformation equation can be obtained as Equation (43) and represented as a set of linear ordinary equations.
L [ u m ( t ; q ) χ m u m 1 ( t ; q ) ] = R m ( u m 1 )
where u m 1 is a vector defined as
u m = { u 0 ( x , t ) , u 1 ( x , t ) , u 2 ( x , t ) , u m ( x , t ) }
and
R m ( u m 1 ) = 1 ( m 1 ) ! m 1 F [ Φ ( t ; q ) ] q m 1 | q = 0
χ m = { 0 m = 1 1 m > 1

4. Application of the Homotopy Analysis Method

Under the transformation τ = Ω t , Equation (32) can be written as
Ω 2 η + ω 1 2 η + k f 2 ( c a 2 η 2 + a 3 η 3 ) = 0
with the following initial conditions:
η ( 0 ) = A ,   η ( 0 ) = 0
Thus, with these initial conditions, the initial guess of η(τ) for a zero-order deformation equation is chosen as follows.
η 0 = A cos τ
According to (46), we should choose the auxiliary linear operator
L [ η ( τ ; q ) ] = Ω 2 ( 2 η ( τ ; q ) τ 2 + η ( τ ; q ) )
which has the property
L ( C 1 sin τ + C 2 τ ) = 0
for any integration constants C1 and C2. Define the following nonlinear operator
N [ η ( τ ; q ) , Ω ] = Ω 2 2 η ( τ ; q ) τ 2 + ω 1 2 η ( τ ; q ) + k f 2 [ c a 2 η 2 ( τ ; q ) + a 3 η 3 ( τ ; q ) ]
The zero-order deformation is constructed as
( 1 q ) L [ Φ ( τ ; q ) η 0 ( τ ; q ) ] = q N [ Φ ( τ ; q ) ]
with the following initial conditions:
Φ ( τ ; 0 ) = η 0 ( t ) ,   Φ ( τ ; 1 ) = η ( τ ; q )
The functions Φ ( τ ; q ) and Ω ( q ) can be expanded as power series of the q using Taylor’s theorem as
Φ ( τ ; q ) = η 0 ( τ , q ) + m = 1 + η m ( τ , q ) q m
Ω ( q ) = Ω 0 + m = 1 + Ω m ( τ , q ) q m
where ηm and Ωm are the m-order deformation derivatives.
N m ( η m 1 , ω m 1 ) = 1 ( m 1 ) ! d m 1 N [ Φ , ω ] d q m 1 | q = 0
For free vibration without damping, periodic motion can be expressed by the following base functions
{ cos ( m τ ) , m = 1 , 2 , 3 , ... }
For the first-order approximation of HAM, N1 is as follows:
N 1 = ( ω 1 2 A Ω 2 A + 3 4 k f 2 a 3 A 3 ) cos τ + 1 2 k f 2 c a 2 A 2 cos 2 τ + 1 4 k f 2 a 3 A 3 cos 3 τ + 1 2 k f 2 c a 2 A 2
Thus, Ω0 can be written as
Ω 0 = ω 1 2 + 3 4 k f 2 a 3 A 2
Equation (59), does not contain the coefficient of quadratic term k f 2 c a 2 , which means that the nonlinear frequency is independent of a quadratic term. The cubic term of Equation (32) plays the dominant role in the nonlinear vibration process.
Solving Equation (52) for the first order solution, η1, is obtained as
η 1 = 1 96 Ω 0 2 [ ( 48 Ω 0 2 A 48 ω 1 2 A 32 k f 2 c a 2 A 2 33 k f 2 a 3 A 3 ) cos τ 16 k f 2 c a 2 A 2 cos 2 τ 3 k f 2 a 3 A 3 cos 3 τ + 48 k f 2 c a 2 A 2 ]
Finally, the approximate solution of Equation (32) is
η = η 0 + η 1 = A cos Ω t + 96 Ω 0 2 [ ( 48 Ω 0 2 A 48 ω 1 2 A 32 k f 2 c a 2 A 2 33 k f 2 a 3 A 3 ) cos τ 16 k f 2 c a 2 A 2 cos 2 τ 3 k f 2 a 3 A 3 cos 3 τ + 48 k f 2 c a 2 A 2 ]

5. Numerical Comparison

In order to verify the effectiveness of the proposed method, numerical comparison and discussion are presented in this section. In Figure 6, the curves at different axial speeds are plotted for vf = 0.8 kf = 1, where Ω is nonlinear frequency and A is amplitude. It is shown that the increasing trend of frequency is consistent with the axial speed. A lower axial speed yields a smaller nonlinear frequency.
Figure 7 shows the dependence of the nonlinear frequency on the amplitude at a fixed axial speed, γ = 3. The dimensionless stiffness is chosen as vf = 0.8. Moreover, we plot the nonlinear frequency for kf =1 and kf = 100, respectively. It is seen that the nonlinear coefficient, kf, is independent of nonlinear frequency. This finding is the same as that for linear frequency, which is independent of nonlinear coefficient kf.
Figure 8 depicts the relationship between nonlinear frequency and flexural stiffness, vf, at a fixed axial speed for γ = 3 and kf = 1. Increasing flexural stiffness causes an increase in nonlinear frequency.
When we compare the analytical solution (61) in Figure 9, let ħ = −1, k f 2 = 1 , γ = 3, A = 0.1, and vf = 0.8, which indicates a weakly nonlinear case. In this figure, the exact and approximate solutions are represented by a solid line and a dashed line, respectively, and they are calculated using the fourth-order Runge Kutta method and the average method. They are in good agreement, indicating that solution (61) is accurate enough for the weakly nonlinear case.
Figure 10 can prove the accuracy of the solutions (61) and the effectiveness of the homotopy analysis method. The parameters chosen in this section are γ = 3, A = 0.1, vf = 0.8, and the nonlinear coefficient is chosen as a larger value, k f 2 = 10 . The approximation, (61), is consistent with the numerical results. Note that the method in this paper is valid for weak and relatively large cases; it does not work for the extreme case where vf approaches infinity.

6. Conclusions

In this paper, we present analytical solutions of natural frequency and response for beams traveling in the supercritical regime. Postbuckling shapes were derived for pinned–pinned and fixed–fixed boundary conditions. The linear vibration around the equilibrium configuration (buckling shape) was analyzed and the corresponding natural frequency and modes were analytically derived. We found that supercritical natural frequency increases as axial speed increases. The equation of motion for a buckling axially moving beam processes quadratic and cubic nonlinearity. Homotopy analysis was employed to solve the nonlinear governing equation of a buckled axially moving beam. Numerical results show that the obtained solution is valid for both weakly and strongly nonlinear vibration cases.

Author Contributions

Conceptualization, Y.J. and T.Y.; methodology, Y.J. and T.Y.; validation, T.Y.; writing—original draft preparation, Y.J.; All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Scientific Research Project of the Tianjin Education Commission grant number 2019KJ121.

Acknowledgments

This study was supported by the Scientific Research Project of the Tianjin Education Commission (No. 2019KJ121).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. The constants Bnk and Cnk.
Table A1. The constants Bnk and Cnk.
Simple SupportsFixed Supports
Bnk λ n λ k λ n 2 λ k 2 [ ( 1 ) n + k 1 ] 2 λ n 2 λ k 2 λ n 4 λ k 4 [ ( 1 ) n + k 1 ]
Bnn00
Cnk0 2 λ n 2 λ k 2 λ n 4 λ k 4 ( λ n σ n λ k σ k ) [ ( 1 ) n + k 1 ]
Cnn λ n 2 λ n σ n ( 2 λ n σ n )
where σ n = sinh λ n sin λ n cosh λ n + cos λ n .

References

  1. Wickert, J. Non-linear vibration of a traveling tensioned beam. Int. J. Non-Linear Mech. 1992, 27, 503–517. [Google Scholar] [CrossRef]
  2. Parker, R.G. Supercritical speed stability of the trivial equilibrium of an axially-moving string on an elastic foundation. J. Sound Vib. 1999, 221, 205–219. [Google Scholar] [CrossRef]
  3. Öz, H.; Pakdemirli, M. Vibrations of an axially moving beam with time-dependent velocity. J. Sound Vib. 1999, 227, 239–257. [Google Scholar] [CrossRef]
  4. Zhu, W.; Mote, C. Free And Forced Response Of An Axially Moving String Transporting A Damped Linear Oscillator. J. Sound Vib. 1994, 177, 591–610. [Google Scholar] [CrossRef]
  5. Chen, L.-Q.; Tang, Y.-Q. Combination and principal parametric resonances of axially accelerating viscoelastic beams: Recognition of longitudinally varying tensions. J. Sound Vib. 2011, 330, 5598–5614. [Google Scholar] [CrossRef]
  6. Mote, C. A study of band saw vibrations. J. Frankl. Inst. 1965, 279, 430–444. [Google Scholar] [CrossRef]
  7. Wickert, J.A.; Mote, C.D. Classical Vibration Analysis of Axially Moving Continua. J. Appl. Mech. 1990, 57, 738–744. [Google Scholar] [CrossRef]
  8. Lin, W.; Qiao, N. Vibration and stability of an axially moving beam immersed in fluid. Int. J. Solids Struct. 2008, 45, 1445–1457. [Google Scholar] [CrossRef]
  9. Chen, L.Q.; Yang, X.D. Vibration and stability of an axially moving viscoelastic beam with hybrid supports. Eur. J. Mech.-A/Solids 2006, 25, 996–1008. [Google Scholar] [CrossRef]
  10. Chen, L.-Q. Analysis and Control of Transverse Vibrations of Axially Moving Strings. Appl. Mech. Rev. 2005, 58, 91–116. [Google Scholar] [CrossRef]
  11. Chen, L.-Q.; Zhao, W.-J.; Zu, J.W. Transient responses of an axially accelerating viscoelastic string constituted by a fractional differentiation law. J. Sound Vib. 2004, 278, 861–871. [Google Scholar] [CrossRef]
  12. Fung, R.-F.; Huang, J.-S.; Chen, Y.-C. The transient amplitude of the viscoelastic travelling string: An integral constitutive law. J. Sound Vib. 1997, 201, 153–167. [Google Scholar] [CrossRef]
  13. Öz, H.; Pakdemirli, M.; Özkaya, E. Transition behaviour from string to beam for an axially accelerating material. J. Sound Vib. 1998, 215, 571–576. [Google Scholar] [CrossRef]
  14. Pakdemirli, M.; Ulsoy, A.G. Stability analysis of an axially accelerating string. J. Sound Vib. 1997, 203, 815–832. [Google Scholar] [CrossRef]
  15. Li, Q.; Liu, W.; Lu, K.; Yue, Z. Flow-induced buckling statics and dynamics of imperfect pipes. Arch. Appl. Mech. 2021, 91, 4553–4569. [Google Scholar] [CrossRef]
  16. Li, Q.; Liu, W.; Lu, K.; Yue, Z. Three-dimensional parametric resonance of fluid-conveying pipes in the pre-buckling and post-buckling states. Int. J. Press. Vessel. Pip. 2021, 189, 104287. [Google Scholar] [CrossRef]
  17. Chen, L.-Q.; Yang, X.-D. Steady-state response of axially moving viscoelastic beams with pulsating speed: Comparison of two nonlinear models. Int. J. Solids Struct. 2005, 42, 37–50. [Google Scholar] [CrossRef]
  18. Chen, L.-Q.; Yang, X.-D. Nonlinear free transverse vibration of an axially moving beam: Comparison of two models. J. Sound Vib. 2007, 299, 348–354. [Google Scholar] [CrossRef]
  19. Ghayesh, M.H. Nonlinear forced dynamics of an axially moving viscoelastic beam with an internal resonance. Int. J. Mech. Sci. 2011, 53, 1022–1037. [Google Scholar] [CrossRef]
  20. Ghayesh, M.H.; Amabili, M. Three-dimensional nonlinear planar dynamics of an axially moving Timoshenko beam. Arch. Appl. Mech. 2013, 83, 591–604. [Google Scholar] [CrossRef]
  21. Ghayesh, M.H. Subharmonic dynamics of an axially accelerating beam. Arch. Appl. Mech. 2012, 82, 1169–1181. [Google Scholar] [CrossRef]
  22. Hwang, S.-J.; Perkins, N. Supercritical stability of an axially moving beam part I: Model and equilibrium analysis. J. Sound Vib. 1992, 154, 381–396. [Google Scholar] [CrossRef]
  23. Hwang, S.-J.; Perkins, N. Supercritical stability of an axially moving beam part II: Vibration and stability analyses. J. Sound Vib. 1992, 154, 397–409. [Google Scholar] [CrossRef]
  24. Ghayesh, M.H.; Amabili, M. Nonlinear dynamics of axially moving viscoelastic beams over the buckled state. Comput. Struct. 2012, 112–113, 406–421. [Google Scholar] [CrossRef]
  25. Ghayesh, M.H.; Kafiabad, H.A.; Reid, T. Sub- and super-critical nonlinear dynamics of a harmonically excited axially moving beam. Int. J. Solids Struct. 2012, 49, 227–243. [Google Scholar] [CrossRef]
  26. Ding, H.; Chen, L.-Q. Galerkin methods for natural frequencies of high-speed axially moving beams. J. Sound Vib. 2010, 329, 3484–3494. [Google Scholar] [CrossRef]
  27. Ding, H.; Chen, L.-Q. Equilibria of axially moving beams in the supercritical regime. Arch. Appl. Mech. 2011, 81, 51–64. [Google Scholar] [CrossRef]
  28. Ding, H.; Zhang, G.-C.; Chen, L.-Q. Supercritical equilibrium solutions of axially moving beams with hybrid boundary conditions. Mech. Res. Commun. 2011, 38, 52–56. [Google Scholar] [CrossRef]
  29. Ding, H.; Zhang, G.; Chen, L. Supercritical vibration of nonlinear coupled moving beams based on discrete Fourier transform. Int. J. Non-Linear Mech. 2012, 47, 1095–1104. [Google Scholar] [CrossRef]
  30. Zhang, G.-C.; Ding, H.; Chen, L.-Q.; Yang, S.-P. Galerkin method for steady-state response of nonlinear forced vibration of axially moving beams at supercritical speeds. J. Sound Vib. 2012, 331, 1612–1623. [Google Scholar] [CrossRef]
  31. Yang, T.Z.; Yang, X.D. Exact solution of supercritical axially moving beams: Symmetric and anti-symmetric configurations. Arch. Appl. Mech. 2013, in press. [Google Scholar] [CrossRef]
  32. Ding, H.; Chen, L.-Q. Natural frequencies of nonlinear vibration of axially moving beams. Nonlinear Dyn. 2010, 63, 125–134. [Google Scholar] [CrossRef]
  33. Liao, S.; Sherif, S. Beyond Perturbation: Introduction to the Homotopy Analysis Method; Chapman and Hall/CRC: New York, NY, USA, 2003; p. 336. [Google Scholar]
Figure 1. The supercritical equilibrium configuration of an axially moving beam.
Figure 1. The supercritical equilibrium configuration of an axially moving beam.
Applsci 12 08555 g001
Figure 2. The subcritical natural frequencies for simple supports.
Figure 2. The subcritical natural frequencies for simple supports.
Applsci 12 08555 g002
Figure 3. The subcritical natural frequencies for fixed supports.
Figure 3. The subcritical natural frequencies for fixed supports.
Applsci 12 08555 g003
Figure 4. The first supercritical natural frequencies for a pinned–pinned case.
Figure 4. The first supercritical natural frequencies for a pinned–pinned case.
Applsci 12 08555 g004
Figure 5. The first supercritical natural frequencies for a fixed–fixed case.
Figure 5. The first supercritical natural frequencies for a fixed–fixed case.
Applsci 12 08555 g005
Figure 6. The variation in nonlinear frequencies with amplitudes above different supercritical axial speeds.
Figure 6. The variation in nonlinear frequencies with amplitudes above different supercritical axial speeds.
Applsci 12 08555 g006
Figure 7. The variation in nonlinear frequencies with amplitudes for different nonlinear coefficient values.
Figure 7. The variation in nonlinear frequencies with amplitudes for different nonlinear coefficient values.
Applsci 12 08555 g007
Figure 8. The variation in nonlinear frequencies with amplitudes for different flexural stiffness values.
Figure 8. The variation in nonlinear frequencies with amplitudes for different flexural stiffness values.
Applsci 12 08555 g008
Figure 9. Comparison between the supercritical analytical solution (57) and the numerical solution for weakly nonlinear cases.
Figure 9. Comparison between the supercritical analytical solution (57) and the numerical solution for weakly nonlinear cases.
Applsci 12 08555 g009
Figure 10. Comparison between the supercritical analytical solution (57) and the numerical solution for strongly nonlinear cases.
Figure 10. Comparison between the supercritical analytical solution (57) and the numerical solution for strongly nonlinear cases.
Applsci 12 08555 g010
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jin, Y.; Yang, T. Buckling Induced Strongly Nonlinear Vibration of Supercritical Axially Moving Beam. Appl. Sci. 2022, 12, 8555. https://doi.org/10.3390/app12178555

AMA Style

Jin Y, Yang T. Buckling Induced Strongly Nonlinear Vibration of Supercritical Axially Moving Beam. Applied Sciences. 2022; 12(17):8555. https://doi.org/10.3390/app12178555

Chicago/Turabian Style

Jin, Yang, and Tianzhi Yang. 2022. "Buckling Induced Strongly Nonlinear Vibration of Supercritical Axially Moving Beam" Applied Sciences 12, no. 17: 8555. https://doi.org/10.3390/app12178555

APA Style

Jin, Y., & Yang, T. (2022). Buckling Induced Strongly Nonlinear Vibration of Supercritical Axially Moving Beam. Applied Sciences, 12(17), 8555. https://doi.org/10.3390/app12178555

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop