Analysis of the thermo-elastic vibration for axially moving Euler beam

The thermo-elastic vibration response of simple supported axially moving Euler beam is investigated. The differential equation of moving beam is established by recourse to Hamilton principle and the thermal effects is considered by introducing the equivalent thermal bending moment. A 2-D transient temperature field is calculated by the alternating-directional implicit (ADI) method and the equivalent thermal moment is calculated numerically. The dimensionless equation is discretized by Galerkin method and the modal analysis of gyroscopic system is used to calculate the forced vibration response. The time-history curve of the beam’s upper middle point is obtained for thermal or non-thermal situations.


Introduction
In the field of engineering applications, many structures are typical axially moving systems, such as aerial cable tramways, band saws and power transmission chains.Mote [1] studied the band saw vibrations by regarding it as an axially moving beam and the frequency curves were obtained.Simpson [2] analyzed the natural frequencies and mode curves of the axially moving beam with clamped boundary, his calculated results show that modes distorted violently as the moving speed increases.Moreover, Wickert [3] brought forth that axially moving beam belongs to gyroscopic system and researched the forced vibration of axially moving beam by modal analysis and Green's function method.The above articles conducted their investigations on the vibration performance and response of the axially moving beam under a certain uniform temperature field.However, in the actual situations, many structures are operated in transient thermal environment, and thus the thermal effect must be considered.Guo [4] considered the coupled thermo-elastic vibration characteristics of the axially moving beam, and acquired the curves of complex frequencies versus the axially moving speed.Wang and Li [5] computed the thermo-elastic responses of axially moving viscoelastic sandwich beams using pseudo-normal modes method, and found that the thermal impact has large influence on the beam's response.However, in the above investigations the vibration responses were mostly calculated under assumed non-constant temperature fields directly, while the others did not mention the specified algorithms they used.
As the problem involves two physical fields, it is important to find an appropriate method to ensure the problem being solved in a reasonable computational cost.In this paper, the temperature fields are calculated firstly by the ADI method, then the superposition method is used to calculate the vibration response.

The basic equations
An axially moving beam with uniform cross-section is considered as shown in Fig. 1.The thermal boundaries are: the upper one suddenly getting a certain temperature rise and remains unchanged, while the others remain their original temperatures unchanged.Consider the beam of density , length , width , height ℎ, cross-section area , moment of inertial , tension , Young's modulus , and coefficient of thermal expansion .The beam travels at a uniform constant velocity between two boundaries.The out-plane motion of the beam is specified by the transverse displacement , .and are the amplitude and frequency of the external excitation respectively.By inserting the thermal moment into the equation derived by Mote [1], we attain the differential equation: for the simple supports boundary conditions: where and is the initial temperature.
where = ⁄ is the velocity of the longitudinal wave.Omitting the bars in the dimensionless variables, Eq. ( 1) can be transformed into dimensionless equation as: When using the Galerkin method to discretize the equation, the displacement function can be expanded as: Introducing Eq. ( 5) into Eq.( 4) and moving the thermal terms into right part, the equation becomes: where: Multiplying the two sides of Eq. ( 6) by and integrating the equation in the dimensionless range [0, 1], let = 4, then we have: Introducing the state variables and general load vectors = , = 2 + , Eq. ( 7) can be arranged as: where = , = − .Eq. ( 9) is a standard form of the gyroscopic system, so it can be solved by the modal superposition method.
According to the research of Meirovitch [6,7], the solution to the eigenvalue problems: consists of pairs of pure imaginary complex conjugates = ± and pairs of complex conjugate eigenvectors and = 1,2. . . .Let and be the real part and imaginary part of and any arbitrary 2 -dimensional vector can be expressed as a linear combination form: For the load vector , the corresponding general coordinates are: = cos − − sin − + 0 cos − 0 sin .

Heat transfer model
For calculating the thermal load in right side, the 2D temperature field of the structure must be solved.The differential equation can be written as the following dimensionless form: , represents the thermal conductivity and the volume heat capacity.The thermal boundaries are written as: 0, , = 0, 1, , = 0, , 1 2 , = 0, , − 1 2 , = 99, , , 0 = 0. (14) The alternating-direction implicit method [8,9] is an effective approach in solving the heat-transfer problems which divides each iteration step into two steps for applying implicit format in and directions respectively, as shown in Fig. 2. Thus, the equation in each step can be transformed into a tri-diagonal form matrix which can be solved by Thomas algorithm efficiently.
At the first step, applying implicit format in direction and explicit format in direction, we have: then, applying formulation in the contrast directions, we have: where = ⁄ , = ⁄ .According to the equations above, the temperature field in any time can be calculated.After numerical integrating and differentiating, and ∂ ∂ ⁄ can be attained.Therefore, the thermal moment will be evaluated.

Numerical examples
Consider a beam with modulus of elasticity = 7.0×10 11 Pa, density = 2778 kg/m 3 .Let the length = 0.5 m, the width = 0.05 m, the height ℎ = 0.05 m, the axial speed = /100, the initial tension = 6000 N, the force amplitude = 500 N, the heat conductivity = 207 W/(m 2 •K), the specific heat = 2484000 J/(m 3 •K) and the coefficient of thermal expansion = 2.39×10 -5 m/K.The whole real-time response will be solved by dividing the heating process into short timesteps and the thermal load can be calculated numerically.At each step, assuming ∂ ∂ ⁄ is constant, the load can be written as: where = 2 2 ⁄ , 0, 2 3 ⁄ , 0 and is constant vector.In each time step period , , the general coordinates , can be written as: where = 2 0 2 3 0 0 0 0 0 .First, the eigenvalue problem can be solved and the frequencies of the system are presented as in Table 1.
Then, the temperature field can be drawn in Fig. 3.For the thermo-elastic issue, the time-region response under mechanical and thermal loads can be solved as follow.The beating phenomenon appears because the frequency of load is near the fundamental frequency of the system.Fig. 4(a) shows that the deflection at middle point presents the beating phenomenon indeed when the temperature is a constant, which shows the validity of the present numerical method.Fig. 4(b) is the result, when a sudden rise of temperature is applied on the upper boundary, which shows a large additional deflection.When the temperature condition disappears, the equilibrium position returns to the original straight line.

Conclusions
The modal superposition method of gyroscopic system has been used to compute the thermo-elastic response of axially moving beam subjected to mechanical and thermal loads.The 2D temperature field is calculated numerically by the ADI method and the differential equation is discretized by the Galerkin method.When the temperature field is a constant, the displacement appears a beating phenomenon under the harmonic load whose frequency is closed to fundamental frequency.The results show that thermal load has important influence on the beam's deflection which deviates from straight equilibrium line during a rise of temperature condition is applied on boundary.The method used herein allows one to obtain a fast solution of thermo-elastic vibration problem.

Fig. 3 .
The transient temperature field in a) 200 and b) 400 dimensionless time

Fig. 4 .
Fig. 4. Displacement in middle point for non-thermal and thermal situations

Table 1 .
The first four order frequencies of axially moving beam