Abstract
The free vibration characteristics of the axially moving Timoshenko beam under compressive load are investigated. The differential governing equation of transverse vibration of the axially moving beam under axial load is established based on the Timoshenko beam theory and Hamilton's principle. The beam system is sorted by the positive definiteness of boundary conditions and the characteristic equations are deduced. The dimensionless complex frequencies of the axially moving Timoshenko beam with different boundary conditions under compressive load are calculated by the proposed modified Galerkin method and the differential quadrature method. The characteristics of the lower order real parts and imaginary parts of the complex frequencies are observed. The effect of compressive load factor, the moving speed and acceleration on the stability of system are discussed.
1. Introduction
Axially moving continuum is involved in industry of engineering and aerospace, such as bandsaws in engineering devices, magnetic bands and moving rockets with high speed, missiles and rigid pipeline in aircrafts and some other components, etc. Most of those devices can be simplified as axially moving beams. Dues to effects of axially moving and excitations of axial load, longitudinal and torsional vibrations can be ignored in contrast with transverse vibration of those systems. Therefore, Study of transverse vibration and dynamic behavior of axially moving beam is of great significance for performance improving and optimization design for these systems.
Many researches have been carried out on axially moving beam. Swope and Ames [1] established a mathematical model for the oscillation of a string as it is traversed and wound on a bobbin, and noted the importance of the ratio of the winding speed to wave velocity. Mote [2] found that the flexural natural frequencies of bandsaw decrease continuously with increasing band velocity and large amplitude equilibrium may be induced by small periodic band tension variations as in the case of the classical string. Lengoc [3] studied the free vibrations of axially moving bandsaws by eigenvalue method and its stability under parametric excitations. Wickert [4] analyzed the free nonlinear vibration of axially moving tensioned beam in sub and supercritical transport speed ranges for the pinnedpinned boundary condition. Öz [5, 6] calculated the natural frequencies for different flexural stiffness with time varying velocity using multiscale approach. Kong and Parker [7] used perturbation method to obtain algebraic equations of natural frequencies and studied the free vibration of the axially moving beam with small flexural stiffness. Yang and Chen [8] also used multiscale approach to determin the natural frequencies of axially moving beam. Lee and Jang [9] studied the dynamic characteristics and stability of an axially moving beam with infinite parts under different boundary conditions by using the spectral element method.
Most of these beam models are based on the EulerBernoulli beam theory, and it appears not so accurate for beams with low lengththickness ratio as a result of without considering shear deformation and rotary inertia of cross section. Simpson [10] firstly investigated the shear effect and analyzed the natural frequencies of the axially moving Timoshenko beam with clampedclamped boundary condition, in which the eigenvalue method was used and the divergence of frequency was shown. Chonan [11] used Laplace transform method to analyze the steadystate response of axially moving Timoshenko beam.
However, no axial load was considered in above two beam models. Lee et al. [12] formulated the spectral element model for the axially moving Timoshenko beam under a uniform axial tension and used exact dynamic stiffness matrix to provide high accuracy of solutions. By applying complex mode approach, Tang et al. [13] analyzed natural frequencies, modes and critical speeds of axially moving Timoshenko beams on different supports under uniform tension. On the other hand, there are some works for static beam with compressive axial load. Bokaina [14, 15] had investigated axial tension and compressive load effect for EulerBernoulli beam; Zhang et al. [16] studied the effect of compressive axial load on transverse vibrations of a doublebeam system and found that the properties of the forced transverse vibrations of the system are significantly dependent on the compressive axial load. Christian and Shu [17] analyzed the influence of axial compressive load on the vibration characteristics of a multipledelaminated beam. These are all for static beam without considering axially moving effect, and the boundary conditions of support are investigated for clamped and pinned situations, which is belong to the positive definite system (PDS), but for the positive semidefinite system (PSDS), such as freefree boundary conditions, there were few researches. Pourtakdoust and Assadian [18] investigated the thrust effect on the vibrational characteristics of flexible guided missiles, but axially motion was not involved. Therefore, to date no publication can be found for the study of moving beam with compressive load, which is still of key engineering significance.
In this paper, the governing equation of transverse vibration of the axially moving beam under compressive axial load is established based on the Timoshenko beam theory and Hamilton's principle at first, then the dimensionless complex frequencies of the axial moving Timoshenko beam with different boundary conditions of PDS and PSDS are calculated numerically by the proposed modified Galerkin method and the differentical quadrature method (DQM). The effect of compressive load factor, the moving speed and acceleration on the stability of system are analyzed, and the critical speeds and loads on different boundary conditions are predicted.
2. Governing equations
For an axially moving Timoshenko beam travelling with speed $V\left(T\right)$ under axial compressive constant load $N$, the stiffness of its cross section is set as $EI\left(X\right)$ and the shear stiffness is $\kappa GA$. $I\left(X\right)$ is the moment of inertia of cross section along length and $E$ is the modulus of elasticity. $\kappa $, $G$, $A$, $L$ and $\rho $ represent the shape factor, shearing modulus, area of cross section, length of beam and density of beam material, respectively. The transverse vibration can be described by the axial coordinate $X$ and time $T$, and the displacements in the axial and transverse coordinates of axial and transverse are $U(X,T)$ and $W(X,T)$, repectively. For uniform beams, rotatory inertia of cross section can be represented by $\rho I$. The model of the current moving Timoshenko beam under compressive load is shown in Fig. 1.
Fig. 1Model of axially moving Timoshenko beam under compressive load
Where $\mathrm{\Theta}$ is the angle of rotation of cross section and can be expressed as $\mathrm{\Theta}=\partial W/\partial X+\mathrm{\Psi}$ for Timoshenko beam, $\mathrm{\Psi}$ is the shear angle. For EulerBernoulli beam, $\mathrm{\Theta}=\partial W/\partial X$.
Without considering the coupling of longitudinal and transverse vibrations, the kinetic energy of the beam involves transverse vibration energy, rotational kinetic energy of cross section and the kinetic energy of rigid motion:
The components of potential energy involves potential energy due to compressive load, and strain energy of bending and shear:
According to Hamilton’s principle, the difference between the kinetic and potential energy is zero:
Substituting Eqs. (1) and (2) into Eq. (3), the expression of governing equations for axially moving Timoshenko beam under compressive load can be established based on the variational analysis in terms of axially moving speed and acceleration, but without derivative to load with respect to the axial coordinate and time:
where:
Substituting Eq. (4) into Eq. (5) yields the partial differential equations which contains transverse vibration only:
By introducing the dimensionless parameters:
Equation (7) can be expressed as dimensionless form:
where the dimensionless parameters ${k}_{1}$, ${k}_{2}$ and ${k}_{5}$associated with $G$ account for the effect of shear distortion, ${k}_{3}$ and ${k}_{4}$ represent the effect of rotary inertia and stiffness of the beam, respectively. Bending moment $M\left(X\right)$ and shear force $Q\left(X\right)$ can be defined as $M\left(X\right)=EI\partial \Theta /\partial X$, ${\mathrm{}Q\left(X\right)=M}_{X}=\kappa GA(\mathrm{\Theta}\partial W/\partial X)$.
Equation (8) contains 1st to 4th order partial differential components of dimensionless axial coordinate and time, and the general solution in the form of separation of variables can be expressed as:
where $\varphi =w\left(x\right)$ is the mode shape function, ${i}^{2}=\text{\u20131}$ and $\lambda $ denotes the natural frequency of beam.
The dimensionless form of different boundary conditions are shown in Table A1 in Appendix, and the static critical loads of each kind of support under compressive load are also given. It should be noted that the loads are true value for positive definite system (pinnedpinned, clampedclamped, clampedfree) and for positive semidefinite system (freefree and pinnedfree), critical load is assumed to be equal to that of pinnedpinned boundary condition for reference. Moreover, the positive semidefinite systems are assumed to be kept at steady state like flight vehicles, rigid motion of instability is restricted.
Boundary conditions used in reference [13] is a simplified form without taking into account of the impact of timevarying characteristics, which means the coefficient of $\lambda $ is zero. Though it makes the equations easy to solve, the result is not so accurate. In the present work, the items of timevarying are included in the boundary conditions as shown in Table A1.
3. Solving methods
3.1. Analytical solution of axial compressive load effect
The governing equation of transverse vibration of EulerBernoulli beam under axial compressive load can be translated to a dimensionless form for static beam as [14]:
where $\lambda $ is the eigenvalue which has a series of solution with the form ${\lambda}_{i}^{2}={\stackrel{}{\omega}}_{i}{L}^{2}\sqrt{\rho A/EI}$, $i=\mathrm{1,2},\dots ,$ and ${\stackrel{}{\omega}}_{i}$ is the $i\text{th}$ natural frequency of beam. ${w}_{xx}/{k}_{4}$ represents the effect of compressive load and the general solution of Eq. (10) is:
where ${C}_{1}\text{,}{C}_{2}\text{,}{C}_{3}$ and ${C}_{4}$ are constant coefficients. The characteristic equations of each boundary conditions are obtained by substituting equations in Table A1 into Eq. (11), as shown in Table A2 in Appendix.
Defining the variable of load ${k}_{s}={k}_{4}^{0.5}=\sqrt{N{L}^{2}/EI}\text{,}$ we have ${{\mathrm{\alpha}}_{1}=\left(0.5{k}_{s}^{2}+\sqrt{0.5{k}_{s}^{4}+{\lambda}^{4}}\right)}^{1/2}$ and ${{\mathrm{\alpha}}_{2}=\left(0.5{k}_{s}^{2}+\sqrt{0.5{k}_{s}^{4}+{\lambda}^{4}}\right)}^{1/2}\text{.}$ Moreover, the dimensionless load factor is defined as the ratio of true load to critical load:
where $\mu $ is a coefficient depending on the boundary condition, it takes ${\pi}^{2}$ for pinnedpinned, pinnedfree and freefree, ${\pi}^{2}/4$ for clampedfree and $4{\pi}^{2}$for clampedclamped, as shown in Table A1. There are no analytical solutions for equations in Table A2 except pinnedpinned and pinnedfree boundary conditions, for those cases the first order eigenvalue ${\lambda}_{1}$ can be obtained with the variation of axial compressive load by numerical solution, as shown in Fig. 2.
It is shown in Fig. 2 that ${\lambda}_{1}$ for clampedclamped and freefree beams are the same. For all boundary conditions ${\lambda}_{1}$ decreases slowly with the increasing axial load in the range of ${k}_{n}=\text{0.01}\text{0.8}$; when the load is near to critical load as ${k}_{n}=\text{0.8}\text{1}$, ${\lambda}_{1}$ decreases rapidly for all boundary conditions, when the load reaches to the critical load, ${\lambda}_{1}$ decreases to 0. Thus, the axial load causes the beam stiffness decreasing and the first order frequency of elastic mode reduces to 0 when it reach to the critical load while the beam undergoes instability state. Translating the natural frequencies to dimensionless forms and the relationship between dimensionless frequencies and axial compressive load can be obtained:
It is shown in Eq. (13) that ${\omega}_{i}$ is proportional to ${\lambda}_{i}^{2}$, but is inversely proportional to load $\sqrt{{k}_{n}}$, and is also dependent on factor$\mu $ of boundary condition. Fig. 3 presents the relationship between the first order dimensionless frequencies ${\omega}_{1}$ and the load factor ${k}_{n}$. Since ${k}_{n}\to 0$, ${\omega}_{i}\to +\infty $, load factor is defined in the range of ${k}_{n}=\text{0.01}\text{1}$ in this paper. ${\omega}_{1}$ is different from ${\lambda}_{1}$ and it decreases rapidly when ${k}_{n}=\text{0.01}\text{0.1}$ and more slowly when ${k}_{n}=\text{0.1}\text{1}$. ${\omega}_{1}$ decreases to 0 when ${k}_{n}=\text{1}$. The results of analytical solution can be used to verify with the numerical methods.
Fig. 2λ1 vs. kn of different boundary conditions
Fig. 3ω1 vs. kn of different boundary conditions
3.2. Modified Galerkin method (MGM)
The Galerkin method can be used to solve Eq. (8), when its trial functions satisfy both force and displacement boundary conditions. The traditional Galerkin method takes mode functions of a static beam without axial load as its trail functions and thus leads to some inaccuracy for solving such problems. In this work, the Glaerkin method is modified by taking the mode functions of beam with load as its trail functions, and thus the present algorithm is called the Modified Galerkin method (MGM). Variation of the first four order roots with load factor ${k}_{n}$ can be obtained by solving equations in Table A2 numerically, as shown in Fig. 4.
Fig. 4First four order eigenvalues vs. load of different boundary conditions
a) Pinnedpinned
b) Freefree
The 2nd, 3rd and 4th order eigenvalues decrease slowly with increasing compressive load and without significant drop like the 1st order eigenvalue when the load is close to the critical load, which indicates that the axial load has little effect on the higherorder eigenvalues other than the first order. Stiffness coefficient ${k}_{4}$ reaches its minimum value when the axial load closes to the critical one, and the contribution of the 1^{st} order eigenvalue caused by compressive load reaches to a maximum value. Obviously, eigenvalues of nonaxial load beam are just a special state of preloaded beam when ${k}_{n}=0$.
Substituting equations of boundary condition in Table A1 into Eq. (11) yields the mode functions of beam under compressive load, as shown in Table A3 in Appendix. Equation (9) is rewritten as $w\left(x,t\right)={\mathbf{\Phi}}^{\mathrm{T}}\mathbf{q}\text{,}\mathbf{}\mathbf{\Phi}=\left[{\varphi}_{1}\left(x\right)\mathbf{}\mathbf{}\mathbf{}{\varphi}_{2}\left(x\right)\dots \mathbf{}{\varphi}_{n}\left(x\right)\mathbf{}\right]\text{,}$$\mathbf{q}=\left[{q}_{1}\left(t\right)\mathbf{}\mathbf{}\mathbf{}{q}_{2}\left(t\right)\dots \mathbf{}{q}_{n}\left(t\right)\mathbf{}\right]\text{,}$${q}_{j}\left(t\right)={e}^{i{\lambda}_{j}t}$. ${\varphi}_{i}\left(x\right)$ is the Mode shape functions $\varphi \left(x\right)$ in Table A3 which is used as the trail functions. From Eq. (8) we have:
where:
${\mathbf{B}}_{3}={\int}_{0}^{1}\mathbf{\Phi}{\mathbf{\Phi}}_{xxx}^{T}dx/{\int}_{0}^{1}\mathbf{\Phi}{\mathbf{\Phi}}^{\mathrm{T}}dx,\mathit{}\mathit{}\mathit{}\mathit{}{\mathbf{B}}_{4}={\int}_{0}^{1}\mathbf{\Phi}{\mathbf{\Phi}}_{xxxx}^{T}dx/{\int}_{0}^{1}\mathbf{\Phi}{\mathbf{\Phi}}^{\mathrm{T}}dx,$
are all square matrix of ${n}_{G}$ th order, ${n}_{G}$ is the number of items of $\varphi \left(x\right)$.
Equation (14) can be expressed in matrix form as:
where $\mathbf{u}={\left[{\mathbf{q}}^{\mathbf{T}}\mathbf{}{\mathbf{q}}_{t}^{\mathbf{T}}\mathbf{}{\mathbf{q}}_{tt}^{\mathbf{T}}\mathbf{}{\mathbf{q}}_{ttt}^{\mathbf{T}}\right]}^{\mathbf{T}}$, $\mathbf{A}$ is the coefficient matrix:
where $0$ and $\mathbf{I}$ are the $n$ order zero matrix and unit matrix, respectively, and as a typical gyroscopic matrix, the translation exists for $\mathbf{A}$ as follow:
The eigenvalue ${\omega}_{i}$ of matrix $\mathbf{A}$ are pairs of complex numbers. According to Eqs. (14) and (15), the number of items for MGM truncation must be larger than the number of natural frequency order to be calculated, therefore, ${n}_{G}$ is determined to be 4 for the accurate calculation for the 1st and 2nd natural frequencies.
3.3. Differential quadrature method
Substituting Eq. (9) into Eq. (8) yields:
According to the differential quadrature method, $r$th order derivative value of function $f\left(x\right)$ at point ${x}_{i}$ can be expressed as a weighted sum of the value of all the nodes:
where ${A}_{ij}^{\left(r\right)}$presents the $r$th order weight coefficients of derivative. ${f}_{j}$ is the function value at ${x}_{j}$, $n$ is the number of nodes. Therefore, the first four order derivatives of mode functions of the nodes can be obtained by interpolation of the function values of the nodes as:
${\varphi}_{xxx}\left({x}_{i}\right)=\sum _{j=1}^{n}{A}_{ij}^{\left(3\right)}{\varphi}_{j},{\varphi}_{xxxx}\left({x}_{i}\right)=\sum _{j=1}^{n}{A}_{ij}^{\left(4\right)}{\varphi}_{j}.$
Nonuniform node distribution is used to discretize the beam and the $\delta $ method is introduced to deal with the boundary conditions, which is to add two nodes of $\delta $ distance from the two ends and renumber the nodes:
$i=\mathrm{3,4}\cdots ,n2,{x}_{n1}=1\delta ,{x}_{n}=1,$
where $\delta ={\text{10}}^{\text{4}}{\text{10}}^{\text{6}}$ [19]. The weight coefficients of DQM can be obtained through the Lagrange interpolation polynomial according to the interpolation principle. Expressions of the coefficients are:
Substituting Eqs. (21) and (22) into Eq. (17) yields:
Also, Eq. (23) can be translated to the matrix form as:
where, ${\mathbf{}\mathbf{B}}^{\left(\mathit{i}\right)}$, $i=\text{0}\text{4}$ are all $n$th order square matrix, $\mathbf{\Phi}={\left[{\varphi}_{1}\mathbf{}\mathbf{}{\varphi}_{2}\mathbf{}\cdots \mathbf{}{\varphi}_{n}\right]}^{T}$ is the model function vector. Expression of elements in ${\mathbf{}\mathbf{B}}^{\left(\mathit{i}\right)}$ are:
${B}_{ij}^{\left(1\right)}=\left(2v{A}_{ij}^{\left(1\right)}+4{k}_{2}v{v}_{t}{A}_{ij}^{\left(2\right)}2{k}_{1}v{A}_{ij}^{\left(3\right)}\right)i,$
${B}_{ij}^{\left(2\right)}={\delta}_{ij}5{k}_{2}{v}_{t}{A}_{ij}^{\left(1\right)}+\left({k}_{1}+{k}_{3}{k}_{2}{k}_{2}{v}^{2}\right){A}_{ij}^{\left(2\right)},$
${B}_{ij}^{\left(3\right)}={k}_{2}vi{A}_{ij}^{\left(1\right)},$
${B}_{ij}^{\left(4\right)}={k}_{2}{\delta}_{ij},$
where ${\delta}_{ij}=\left\{\begin{array}{c}0,i\ne j,\\ 1,i=j.\end{array}\right.$
For nodes at the two ends, i.e. $i,j=\mathrm{1,2},n1,n\text{,}$ since that the introduction of the timevarying items, not only elements in ${\mathbf{}\mathbf{B}}^{\left(0\right)}$ are affected, but also in ${\mathbf{}\mathbf{B}}^{\left(1\right)}$ and ${\mathbf{}\mathbf{B}}^{\left(2\right)}$. According to Table A1, elements in ${\mathbf{}\mathbf{B}}^{\left(0\right)}\text{,}$${\mathbf{B}}^{\left(1\right)}$ and ${\mathbf{}\mathbf{B}}^{\left(2\right)}$ are improved after introducing the boundary conditions with the axial compressive load. The DQM form of boundary conditions are shown in Table A4 in Appendix.
Then solving of Eq. (24) is transformed to the generalized eigenvalue problem. The eigenvalues of Eq. (24) represents the dimensionless natural frequency of beam system.
4. Verification and discussion
An uniform thinwalled beam with rectangular crosssection is used for solutions and discussions. The beam’s configuration and material parameters are: $L=\text{1}\text{}\text{m}$, $\mathrm{A}=\text{0.02}\text{}{\text{m}}^{\text{2}}\text{,}$$I=\text{1.67\xd7}{\text{10}}^{\text{5}}{\text{m}}^{\text{4}}\text{,}$$A=\text{0.02}\text{}{\text{m}}^{\text{2}}\text{,}$$E=\text{100}\text{}\text{GPa}\text{,}$$\mu =\text{0.3,}$$\rho =\text{1.89\xd7}{\text{10}}^{\text{3}}\text{kg}/{\text{m}}^{\text{3}}$.
4.1. Verification of numerical methods
The beam in engineering can be classified into two types according to the positive definiteness property of boundary condition as mentioned before. Take one from each type of boundary conditions for example, i.e. the pinnedpinned beam for PDS, and the freefree beam for PSDS. The numerical results of MGM and DQM are compared with the analytical ones under $v=\mathrm{}{v}_{t}=\text{0}$, as shown in Fig. 5.
Fig. 5 shows that numerical results of MGM and DQM both agree with analytical results in the range of ${k}_{n}=\text{0.01}\text{}1$, thus these two numerical methods are both applicable for different boundary conditions. Fig. 6 shows the results of DQM and MGM for the pinnedpinned beam in the range of $v=\text{0}\text{4}$, the results for 1st and 2nd natural frequencies of these two methods are shown to agree with each other well. The DQM is used in the following numerical analysis while the number of nodes $n$ is set to 13 for convergence.
Fig. 5Comparison of DQM, MGM and analytical solutions ω1 vs. kn
Fig. 6Comparison of DQM, MGM ω1,2 vs. v
4.2. Effect of axial compressive load
Fig. 7 indicates the variation of the 1st dimensionless complex frequency of PDS with axial compressive load. The complex frequencies of PDS are all real numbers with imaginary parts equal to 0 when axially moving speed is 0. The real part decreases rapidly in the range of ${k}_{n}=\text{0.01}\text{}\text{0.2}$ and slowly when ${k}_{n}=\text{0.2}\text{}1\text{.}$ When the load reaches to the critical load as ${k}_{n}=\text{1}$, $Re\left({\omega}_{1}\right)$ reduce to 0, which indicates the critical load factor ${k}_{ncr}=\text{1}$ when$v=\text{0}$. But for $v=\text{1}$, the imaginary part of clampedfree beam is negative. With increasing moving speed for all boundary conditions the real part of the 1st frequency reduces to zero and it becomes unstable by the divergence instability of their imaginary parts when the load reaches the critical loads, which decreases to around half of the value of static beam.
Fig. 7The 1st dimensionless complex frequency of PDS vs. load vt=0, v=0, 1. v=0 ; v=1∙
a)
b)
Fig. 8 indicates the variation of the 1st dimensionless complex frequency of PSDS with axial compressive load. As the cases in PDS, the frequency reduces with the increasing of moving speed. When $v=\text{1}$, the critical loads of freefree and pinnedfree beams reduce to $\text{0.8}{N}_{cr}$ and $\text{0.6}{N}_{cr}$, repectively. The imaginary parts $Im\left({\omega}_{1}\right)$ of those two conditions keep 0 before critical load when $v=\text{0}$. When $v=\text{1}$, it branches for freefree beam and the bifurcation point is the critical load point, but $Im\left({\omega}_{1}\right)$ of pinnedfree beam keeps rising with increasing load and it intersects the horizontal axis at the critical load point.
Fig. 8The 1st dimensionless complex frequency of PSDB vs. load vt=0, v=0, 1. v=0 ; v=1∙
a)
b)
Fig. 9 indicates the variation of the imaginary parts of the fundamental modes of PSDS with axial compressive load. There exists two fundamental modes for freefree beam, one for in plane translational motion and the other for in plane rotation, while the pinnedfree beam only has one for rotation. The real part of fundamental mode keep 0 for the two boundary conditions, and only the imaginary parts need to be discussed. For freefree beam, the imaginary frequency of the first mode of translational motion keeps 0, while that of the second mode of in plane rotation has slight increasing with ${k}_{n}$ for $v=\text{0}$ and has an obvious increasing when ${k}_{n}>\text{0.5}$ for $v=\text{1}\text{.}$ For pinnedfree beam, the imaginary frequency of in plane rotation keeps increasing in a small range with respect to ${k}_{n}$ for $v=\text{0}\text{,}$ but it has a complex variation of bifurcationcombinationrebifuration phenomenon with the increasing of ${k}_{n}$ for $v=\text{1}$.
Fig. 9Imaginary part of fundamental frequencies of PSDB vs. load vt=0, v=0, 1. v=0 ; v=1∙
4.3. Effect of axially moving
Figs. 1012 show the variations of the first three dimensionless complex frequencies of PDSs boundary cases with axially moving speed when ${v}_{t}=\text{0}$, ${k}_{n}=\text{0.1,}\text{}\text{0.5}$. For all the three boundary cases, the real parts of $\omega $, i.e. $Re\left({\omega}_{1}\right)$ decrease with the increasing of speed. The $Re\left({\omega}_{1}\right)$ reduces to 0 when the moving speed rises to the critical value ${v}_{cr}$ and the supercritical value of speed ${v}_{scr}$ The ${v}_{cr}$ of three cases are 6.5, 2.8 and 3.9 when ${k}_{n}=\text{0.1}$, the values decrease to 1.0, 1.1 and 1.3 when ${k}_{n}=\text{0.5}$, which indicates that axial the compressive axial load affects the critical moving speed a lot. For ${v}_{cr}<v<{v}_{scr}$, the 2nd and 3rd mode remain stable, but 1st mode is unstable with divergence instability. Special attention should be paid to the left parts of Fig. 10 for ${k}_{n}=\text{0.5}$ and Fig. 11 for ${k}_{n}=\text{0.1,}\text{}\text{0.5}$. The 1st order regains stability when $v$ reaches ${v}_{scr}$, after that, the real parts of 1st and 2nd complex frequencies merge to each other while the imaginary parts bifurcate with positive and negative values, which indicates that the two order modes couple with each other and the beam undergoes coupledmode flutter. For the clampedfree beam shown in Fig. 12, the real parts of 1st and 2nd complex frequencies get close but not merge together when $v$ reaches ${v}_{scr}$, and The $Im\left(\omega \right)$ of clampedfree beam keep negative as moving speed increasing and is different from the other two PDSs as a result of its asymmetric boundary condition.
Fig. 10The 1st, 2ndand 3rd dimensionless complex frequencies of pinnedpinned vs. speed vt=0, kn=0.1, 0.5
a)
b)
Fig. 11The 1st, 2ndand 3rd dimensionless complex frequencies of clampedclamped vs. speed vt=0, kn=0.1, 0.5
a)
b)
Fig. 12The 1st, 2ndand 3rd dimensionless complex frequencies of clampedfree vs. speed vt=0, kn=0.1, 0.5
a)
b)
Fig. 13The 1st, 2ndand 3rd dimensionless complex frequencies of freefree vs. speed vt=0, kn=0.1, 0.5
a)
b)
Figs. 13 and 14 show the variations of the first three order dimensionless complex frequencies of PSDSs boundary cases with axially moving speed when ${v}_{t}=\text{0}$, ${k}_{n}=\text{0.1,}\text{}\text{0.5}$. For freefree beam shown in Fig. 13, ${v}_{cr}$ reduces from 5.9 to 2.2 when the compressive axial load rises from $\text{0.1}{N}_{cr}$ to $\text{0.5}{N}_{cr}\text{.}$ For ${k}_{n}=\text{0.5}$, when $\text{2.2}<v<\text{2.6}$, the positive parts of 1st complex frequency keeps stable with value 0 while imaginary parts branches. When $v>\text{2.6}$, $Re\left({\omega}_{1}\right)$ reappears and rises up to couple with $Re\left({\omega}_{2}\right)$ at $v=\text{3.9}$. $Re\left({\omega}_{3}\right)$ decreases with load increasing while $Im\left({\omega}_{2}\right)$ and $Im\left({\omega}_{3}\right)$ both remains 0.
Fig. 14The 1st, 2ndand 3rd dimensionless complex frequencies of pinnedfree vs. speed vt=0, kn=0.1, 0.5
a)
b)
For the pinnedfree beam shown in Fig. 14, ${v}_{cr}$ is 6.6 under the load of $\text{0.1}{N}_{cr}$, and $Im\left({\omega}_{1}\right)$ branches and is asymmetric. The beam undergoes supercritical state when$v$ reaches to 6.9, 1st and 2nd frequencies are close to each other but do not couple. ${v}_{cr}$ reduces to 2.6 when the compressive load increases to $\text{0.5}{N}_{cr}\text{.}$ For all cases, $Re\left({\omega}_{2}\right)$ and $Re\left({\omega}_{3}\right)$ are positive while $Im\left({\omega}_{2}\right)$ and $Im\left({\omega}_{3}\right)$ are both negative and decrease with increasing speed.
Fig. 15 show the variations of the fundamental modes of PSDSs boundary cases with axially moving speed when ${v}_{t}=\text{0}$, ${k}_{n}=\text{0.1,}\text{}\text{0.5}$. The 1st fundamental frequency of freefree beam remains 0, while the 2nd order frequency decreases slowly with increasing moving speed but rises sharply when $v$ near to ${v}_{cr}$ for ${k}_{n}=\text{0.1.}$ The $Im\left({\omega}_{0}\right)$ branches when $v>\text{0}$ and merges at $v=\text{4.6}$. The law of $Im\left({\omega}_{0}\right)$ for those two boundary cases are similar, when ${k}_{n}=\text{0.5}$ the turning point moves forward.
Fig. 15Fundamental frequencies of PSDB vs. speed vt=0, kn=0.1, 0.5
a) Freefree
b) Pinnedfree
Figs. 16 and 17 show the variations of the first three dimensionless frequencies of pinnedpinned and pinnedfree beams with axially moving acceleration when ${k}_{n}=\text{0.1, 0.5}$, $v=\text{0,}\text{}\text{3}$. Fig. 16 indicates that the real parts of each order frequencies rise with the increasing moving acceleration when axial load is $\text{0.1}{N}_{cr}$. When the load increase to $\text{0.5}{N}_{cr}$, $Re\left({\omega}_{1}\right)$ turns to decrease for ${v}_{t}>\text{400}$ and reaches 0 when ${v}_{t}=\text{67}0$, while the 2nd and 3rd modes couple each other which indicates that the beam undergoes coupledmode flutter. For pinnedfree beam, fundamental mode merges with the 1st elastic mode when ${v}_{t}=\text{530}$, and the 2nd and 3rd modes couple at ${v}_{t}=\text{670}$. Thus, the relationship between elastic modes and axially moving acceleration is complex under a nonlinear coupling with the axial compressive load.
From Fig. 17, one can see that the real parts of each order rise with increasing moving acceleration. For ${v}_{t}>\text{500}$$Re\left({\omega}_{1}\right)$ of pinnedfree beam is larger at $v=\text{3}$ than $v=\text{1}$, which indicates that the effect of coupling speed and acceleration is not always linear.
Fig. 16Real part of the first three frequencies vs. acceleration kn=0.1, 0.5, v= 0
a) Pinnedpinned
b) Pinnedfree
Fig. 17Real part of the first three frequencies vs. acceleration kn=0.1, v=0, 1
a) Pinnedpinned
b) Pinnedfree
Fig. 18Fundamental frequencies of pinnedfree beam vs. acceleration kn=0.1, 0.5, v=0, 3
Fig. 18 shows the variations of the fundamental frequencies of pinnedfree beam with axially moving acceleration when ${k}_{n}=\text{0.1, 0.5}$, $v=\text{0}$ and $v=\text{0,}\text{}\text{3}$, ${k}_{n}=\text{0.1}$. The $Im\left({\omega}_{0}\right)$ branches when ${v}_{t}=\text{240}$ and this point agrees with $Re\left({\omega}_{1}\right)$. The $Im\left({\omega}_{0}\right)$ remains 0 when $v=\text{0}$ and decreases with the increasing of acceleration when $v>\text{0}$.
5. The critical loads and speeds
When the first elastic frequency ${w}_{1}$ reduces to zero, the beam reaches a critical state, and Eq. (8) becomes a timeindependent equilibrium linear equation:
It can be simplified without taking into account of axially acceleration i.e. ${v}_{t}=0$:
The roots of the characteristic equation of Eq. (27) is:
And the general solution of Eq. (27) is:
Equation (29) contains two parts, the first part ${C}_{1}+{C}_{2}x$ denotes two types of foundational motion, i.e. lateral rigid motion, and the inplane rotation. For PDS, the two foundational motions are limited. For PSDS, freefree beam have both foundational motions and pinnedfree beam only has inplane rotation. The second part denotes elastic modes.
Equation (27) can be expressed to a similar form of Eq. (10), where ${\alpha}_{1}^{2}=\frac{{v}^{2}+1}{{k}_{4}{k}_{1}{k}_{1}{v}^{2}}$, ${\alpha}_{2}=0$. Then $n$th eigenvalue can be obtained ${\alpha}_{n}$ are obtained: ${\alpha}_{n}=n\pi $, for pinnedpinned, pinnedfree and freefree beams, ${\alpha}_{n}=2n\pi $ for clampedclamped beam, ${\alpha}_{n}=(2n1)\pi /2$ for clampedfree beams. By defining a new parameters ${k}_{6}=\frac{E{I}^{2}}{\kappa G{A}^{2}{L}^{4}}$, then ${k}_{1}={k}_{6}\mu {k}_{n}$, ${k}_{4}=1/\mu {k}_{n}$. Taking these items into the expression of ${\alpha}_{1}\text{,}$ the relationship between critical moving speed and axial compressive load can be obtained as:
And the variation of critical load to moving speed can also be obtained at the same time.
Fig. 19 shows the relationship between critical load and speed of clampedfree beam. It indicates that with the effect of axially moving and axial compressive load, the beam reaches critical state more easily. The critical value of speed or load can be approximately predicted from these results.
Fig. 19Critical load and speed of clampedfree beam
a) Critical speed vs. load
b) Critical load vs. speed
6. Conclusions
In this paper, the effect of compressive load on complex frequencies and stabilities of axially moving Timoshenko beam with different supports are studied，expressions of the critical load and speed for various boundary conditions are determined. Results of the study can be summarized as follows:
The real part of first three complex frequencies of five different boundary conditions decrease with the increasing compressive load and axially moving speed. The first two order modes of PDS conditions undergo coupledmode flutter with increasing critical moving speed by divergent instability, while those of PSDS get closer but not merge together also. Rigid mode of PSDS undergoes instability with the increasing of compressive load and speed. The moving acceleration causes increasing of elastic and foundational frequencies.
Critical moving speed decreases with the increasing of the axial load, the critical load decreases with the increasing of axially moving speed, the coupling of these two effects causes the beam to reach instability state more easily.
References

Swope R. D., Ames W. F. Vibrations of a moving threadline. Journal of the Franklin Institute, Vol. 275, Issue 1, 1963, p. 3655.

Mote Jr. C. D. A study of band saw vibrations. Journal of the Franklin Institute, Vol. 279, Issue 6, 1965, p. 430444.

Lengoc L. Vibration of bandsaws. University of Canterbury, 1990.

Wickert J. A. Nonlinear vibration of a traveling tensioned beam. International Journal of NonLinear Mechanics, Vol. 27, Issue 3, 1992, p. 503517.

Öz H. R., Pakdemirli M. Vibrations of an axially accelerating beam with small flexural stiffness. Journal of Sound and Vibration, Vol. 234, 2000, p. 521535.

Öz H. R., Pakdemirli M., Boyacm H. NonLinear vibrations and stability of an axially moving beam with timedependent velocity. International Journal of NonLinear Mechanics, Vol. 36, Issue 3, 2001, p. 107115.

Kong L., Parker G. Approximate eigensolutions of axially moving beams with small flexural stiffness. Journal of Sound and Vibration, Vol. 276, Issue 2, 2004, p. 459469.

Yang X. D., Chen L. Q. Determination of the natural frequencies of axially moving beams by the method of multiple scales. Journal of Shanghai University, Vol. 11, Issue 3, 2007, p. 251254.

Lee U., Jang I. On the boundary conditions for axially moving beams. Journal of Sound and Vibration, Vol. 306, 2007, p. 675690.

Simpson A. Transverse modes and frequencies of beams translating between fixed end supports. Journal of Mechanical Engineering Science, Vol. 15, 1973, p. 159164.

Chonan S. Steady state response of an axially moving strip subjected to a stationary lateral load. Journal of Sound and Vibration, Vol. 107, Issue 1, 1986, p. 155165.

Lee U., Kim J. H., Hyungmi O. Spectral analysis for the transverse vibration of an axially moving Timoshenko beam. Journal of Sound and Vibration, Vol. 271, Issue 3, 2004, p. 685703.

Tang Y. Q., Chen L. Q., Yang X. D. Natural frequencies modes and critical speeds of axially moving Timoshenko beams with different boundary conditions. International Journal of Mechanical Sciences, Vol. 50, Issue 3, 2008, p. 14481458.

Bokaian A. Natural frequencies of beams under compressive axial loads. Journal of Sound and Vibration, Vol. 126, Issue 1, 1988, p. 4965.

Bokaian A. Natural frequencies of beams under tensile axial loads. Journal of Sound and Vibration, Vol. 142, Issue 3, 1990, p. 481498.

Zhang Y. Q., Lu Y., Ma G. W. Effect of compressive axial load on forced transverse vibrations of a doublebeam system. International Journal of Mechanical Sciences, Vol. 50, 2008, p. 299305.

Christian N. D., Shu D. W. Free vibration analysis of multiple delaminated beams under axial compressive load. Journal of Reinforced Plastics and Composites, Vol. 28, Issue 11, 2009, p. 13651381.

Pourtakdoust S. H., Assadian N. Investigation of thrust effect on the vibrational characteristics of flexible guided missiles. Journal of Sound and Vibration, Vol. 272, 2004, p. 287299.

Sung K. J., Charles W. B., Alfred G. S. Application of differential quadrature to static analysis of structural components. International Journal for Numerical Methods in Engineering, Vol. 28, Issue 3, 1989, p. 561577.