Advances in Acoustics and Vibration

Volume 2014 (2014), Article ID 410851, 14 pages

http://dx.doi.org/10.1155/2014/410851

Research Article

Vibration Analysis of Hollow Tapered Shaft Rotor

P. M. G. Bashir Asdaque and R. K. Behera

Department of Mechanical Engineering, National Institute of Technology Rourkela, Odisha 769008, India

Received 27 December 2013; Accepted 6 March 2014; Published 28 April 2014

Academic Editor: Abdelkrim Khelif

Copyright © 2014 P. M. G. Bashir Asdaque and R. K. Behera. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Shafts or circular cross-section beams are important parts of rotating systems and their geometries play important role in rotor dynamics. Hollow tapered shaft rotors with uniform thickness and uniform bore are considered. Critical speeds or whirling frequency conditions are computed using transfer matrix method and then the results were compared using finite element method. For particular shaft lengths and rotating speeds, response of the hollow tapered shaft-rotor system is determined for the establishment of dynamic characteristics. Nonrotating conditions are also considered and results obtained are plotted.

1. Introduction

Shaft is a major component of any rotating system, used to transmit torque and rotation. Hence the study shaft-rotor systems has been the concern of researchers for more than a century and will continue to persist as an active area of research and analysis in near future. Geometry of shaft is of the main concern during the study of any rotating system. Most papers related to shaft-rotor systems consider cylindrical shaft elements for study and analysis of rotating systems. The first idea of transfer matrix method (TMM) was compiled by Holzer for finding natural frequencies of torsional systems and later adapted by Myklestad [1, 2] for computing natural frequencies of airplane wings, coupled in bending and torsion. Gyroscopic moments were first introduced by Prohl [3] for rotor-bearing system analysis. Lund [4] used complex variables as the next significant advancement in the method. An improved method for calculating critical speeds and rotor stability of turbo machinery was investigated by Murphy and Vance [5]. Whalley and Abdul-Ameer [6] used frequency response analysis for particular profiled shafts to study dynamic response of distributed-lumped shaft rotor system. They studied the system behavior in terms of frequency response for the shafts with diameters which are functions of their lengths. They derived an analytical method which uses Euler-Bernoulli beam theory in combination with TMM. On the other hand, there are large numbers of numerical applications of finite element techniques for the calculation of whirling and the computation of maximum dynamic magnitude. In this regard, Ruhl and Booker [7] modeled the distributed parameter turbo rotor systems using finite element method (FEM). Nelson and McVaugh [8] reduced large number of eigenvalues and eigenvectors identified, following finite element analysis, and the erroneous modes of vibration predicted were eliminated. Nelson [9] again formulated the equations of motion for a uniform rotating shaft element using deformation shape functions developed from Timoshenko beam theory including the effects of translational and rotational inertia, gyroscopic moments, bending and shear deformation, and axial load. Greenhill et al. [10] derived equation of motion for a conical beam finite element form Timoshenko beam theory and include effects of translational and rotational inertia, gyroscopic moments, bending and shear deformation, axial load, and internal damping. Genta and Gugliotta [11] also analyzed element with annular cross-section based on Timoshenko beam theory having two degrees of freedom at each node. Mohiuddin and Khulief [12] derived a finite element model of a tapered rotating cracked shaft for modal analysis and dynamic modeling of a rotor-bearing system, based on Timoshenko beam theory, that is, included shear deformation and rotary inertia. Rouch and Kao [13] presented numerically integrated formulation of a tapered beam element for rotor dynamics.

In this era of machines, tapered shafts are widely used for rotating systems. Using the approach of Whalley and Abdul-Ameer [6], the dynamic analysis of hollow tapered shaft-rotor has been done. Later the results obtained were compared with that obtained from finite element method. The effect of length and speed on the dynamic analysis of the hollow tapered shaft-rotor system is also clearly shown. Frequency response of the rotor system for an impulse of unit force at the free end is determined in terms of critical speeds for various rotor speeds and shaft lengths.

2. Transfer Matrix Method

2.1. Shaft Model

The shaft model is derived in matrix form Whalley and Abdul-Ameer [6] as

where

The complete derivation is present in [6].

2.2. Rigid Disk

The output vector from the shaft will become the input for the rigid disk model, as shown in Figure 1; that is, for disk model, we have

where , , , and are the deflections, slopes, bending moments, and shear forces at the free and fixed end, respectively.

Figure 1: Vibrating shaft element.

Hence, writing in matrix form, we have

where

Form transfer matrix method [6]

where

and input-output vectors relationship is given by

After applying the boundary conditions for cantilever beam, deflection at the free end is obtained and hence leads to transfer function.

3. Finite Element Method

The vector of nodal displacements is given by

So, each element is having eight degrees of freedom.

3.1. Rigid Disc

Rigid disk is having two translations and two rotations in and direction, respectively (considering coordinate in axial direction). For constant spin condition, the Lagrangian equation of motion is given by

where

The forcing term may include mass unbalance and other external forces.

3.2. Finite Shaft-Rotor Element

The rotor-shaft element considered here has eight degrees of freedom, that is, four degrees of freedom per node as in Nelson and McVaugh [8]. For constant spin condition, the Lagrangian equation of motion is given by

where

Except skew-symmetric gyroscopic matrix , others are symmetric matrices. Since the element is linearly tapered, area and inertias are the function of the shaft-rotor length. The translational shape function is given by

where

The rotational shape function is given by

The element matrices are assembled together to get the equation of motion for the complete system.

4. Numerical Results

4.1. Rotating Condition

By using the transfer matrix approach as in the paper of Whalley and Abdul-Ameer, we will ultimately get the transfer function which will be plotted. The FEM will be applied and then is compared with the TMM approach of Whalley and Abdul-Ameer [6].

Example 1. Let us consider a cantilever tubular shaft with uniform thickness and a disc at the free end with downward unit force, , on the disc as shown in Figure 2. The default values of various parameters are tabulated in Table 1.

Table 1: Various parameters of the system shown in Figure 2 and their default values.

Parameters
Values

Length of the shaft rotor, (m)
0.1

Mass of the disk, (Kg)
0.7443

Diameter of the disk, (m)
0.09

Young’s modulus of elasticity, (GPa)
209

Density of the material, (Kg/m³)
7800

Rotational speed, (rpm)
10000

Thickness of the hollow shaft, (m)
0.002

Beginning radius, (m)
0.0050

End radius, (m)
0.0037

Figure 2: Hollow tapered shaft disc with uniform thickness.

Transfer function for tubular shaft with constant thickness as shown in Figure 2 for default values is given by

Bode plots for different lengths and rotating speed have been plotted using MATLAB software, as shown in Figures 3 and 4.

Figure 3: Bode plot for different lengths with TMM.

Figure 4: Bode plot for various speeds with TMM.

Applying FEM on the same system, we get mass, gyroscopic, and stiffness matrices. A finite hollow tapered shaft element is shown in Figure 5.

Figure 5: Hollow tapered shaft finite element with uniform thickness.

The stiffness matrix for hollow tapered shaft element with uniform thickness is given by is given by

where elements of the stiffness matrices are

Translational mass matrix is given by

where elements of translational mass matrix are given by

Rotational mass matrix is given by

where elements of rotational mass matrix are given by

Gyroscopic matrix is given by

where the elements of the gyroscopic matrix are given by

Discretizing the tapered shaft into six elements as shown in Figure 6 and then assembling, we get the assembled equation of motion

where is the assembled mass matrix containing both the translational and rotational mass matrices.

Figure 6: Discretized shaft element.

The assembled equation of motion is arranged in the first order state vector form

where

The shaft rotor has been discretized into six elements of equal length. Hence the order of assembled matrices, after applying the fixed-free boundary condition, is .

MATLAB program is used to find the bode plot for different values of shaft length and rotor speed as shown in Figures 7 and 8 and are found to be in good agreement with bode plots found using TMM as shown in Figures 3 and 4.

Figure 7: Bode plot for different shaft lengths (FEM).

Figure 8: Bode plot for various rotor speeds (FEM).

Example 2. Let us consider a cantilever hollow shaft with uniform bore and a disc at the free end, as shown in Figure 9. The values of various parameters are tabulated in Table 2.

Table 2: Various parameters of the system shown in Figure 9 and their default values.

Parameters
Values

Length of the shaft rotor, (m)
0.1

Mass of the disk, (Kg)
0.7443

Diameter of the disk, (m)
0.09

Young’s modulus of elasticity, (GPa)
209

Density of the material, (Kg/m³)
7800

Rotational speed, (rpm)
10000

Inner radius of hollow shaft, (m)
0.001

Beginning radius, (m)
0.0050

End radius, (m)
0.0037

Figure 9: Hollow tapered shaft disc with uniform bore and vertically downward force on the disc.

The transfer function for hollow shaft with constant thickness for default values is given by

The bode plots for different lengths and rotating speeds have been plotted using MATLAB software, as shown in Figures 10 and 11.

Figure 10: Bode plot for changing length for hollow tapered shaft rotor with uniform bore.

Figure 11: Bode plot for different speeds for hollow tapered shaft rotor with uniform bore.

Applying FEM in the same system, we get mass, gyroscopic, and stiffness matrices. A finite hollow tapered shaft element with uniform bore is shown in Figure 12.

Figure 12: Hollow tapered shaft finite element with uniform bore.

The stiffness matrix for hollow tapered shaft with uniform bore is given by

where elements of the stiffness matrices are

Translational mass matrix is given by

where elements of translational mass matrix are given by

Rotational mass matrix is given by

where the elements of the rotational mass matrix are given by

The gyroscopic matrix is given by

where the elements of the gyroscopic matrix are

As proceeded in Example 1, bode plots are obtained for various shaft lengths and rotor speeds as shown in Figures 13 and 14 and are found to be in good agreement with bode plots found using TMM as shown in Figures 10 and 11.

Figure 13: Bode plot for different shaft lengths (FEM).

Figure 14: Bode plot for different rotating speeds (FEM).

4.2. Nonrotating Conditions

Bode plot for nonrotating (1 or 2 rpm) tapered shaft-rotor system is slightly different from rotating conditions in terms of amplitude. Hollow shaft with uniform thickness is considered. Bode plots are obtained for zero rpm as shown in Figure 15 with TMM.

Figure 15: Bode plot at zero rpm with TMM for uniform thickness hollow tapered shaft-rotor system.

Applying FEM, then for zero rpm we get the bode plot as shown in Figure 16.

Figure 16: Bode plot at zero rpm with FEM for uniform thickness hollow tapered shaft-rotor system.

5. Conclusions

Shaft geometry plays one of the important roles in dynamic characteristics of rotating systems. Vibration analysis with the help of bode plots has been done for hollow tapered shaft-rotor system. Both TMM and FEM have been used for the purpose. The equation of motion for a tapered beam finite element has been developed using Euler-Bernoulli beam theory. Mass, stiffness, and gyroscopic matrices are found and values of all these elements are stated in a systematic manner for ease of understanding. The results obtained from both methods are compared and are found to be in good agreement. However, the above procedures show that the method of TMM is simpler in calculations. Two types of hollow tapered shafts have been analyzed, that is, one with uniform thickness and another with uniform bore.

The length of the shaft is a vital parameter that affects the frequency response of the shaft-rotor system. As shown in Figures 3, 7, 10, and 13 there is increase in amplitude of vibration for increased value of shaft length, while reducing the whirling speed of shaft. The system exhibits this behaviour due to the bending effect of the shaft and stiffness change. For changing lengths, there are large differences in frequencies even for small increase in lengths. Bode plots obtained in Figures 4, 8, 11, and 14 show that rotating speeds have very little effect on the critical frequencies; however with increasing speed, the amplitude is lowered due to gyroscopic couple. The phase angle changes abruptly for lower value of shaft speeds than higher speeds. This is due to the fact that there is reduction in gyroscopic couple as the rotational speed of the shaft decreases thereby giving larger amplitude of vibration. Nonrotating conditions are also shown in Figures 15 and 16 as bode plots which show that any rotating system at a very low speed vibrates with high amplitudes due to lack of gyroscopic couple.

Effects of bearing may be included in the problem. Geared systems and other rotary elements can be mounted instead of discs and further calculations can be made. Multidiscs and other complex problems can be solved using these methods.

Notations

| Compliance per unit of length (function)
| Inertia per unit of length (function)
| Length of shaft
| Modulus of elasticity
| System model matrix
| Bending moment in plane (function)
| Shear force (function)
| Vertical deflection of shaft (function)
| Slope of the shaft (function)
| Mass moment of inertia (function)
| Material density
| Shaft polar moment of inertia
| Whirling frequency
| Wave propagation factor (function)
| Shaft-rotor rotational speed
| Rigid rotor model matrix
| Polar moment of inertia of disc
| Beginning radius of the shaft element
| Inner radius of the shaft element
| End radius of the shaft element
| Mass of the disc attached at free end
| Translational mass matrix for disc
| Rotational mass matrix for disc
| Gyroscopic matrix for disc
| Translational mass matrix
| Rotational mass matrix
| Gyroscopic matrix for element
| Stiffness matrix for shaft element
| External force matrix
| Cross-sectional area (function)
| Diametral inertia (function)
| Polar inertia (function)
| Rotational speed in rad/s
| Translational shape function matrix
| Rotational shape function matrix.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

N. O. Myklestad, “A new method for calculating natural modes of uncoupled bending vibration of airplane wings,” Journal of the Aeronautical Sciences, vol. 11, no. 2, pp. 153–162, 1944.

N. O. Myklestad, “New method of calculating natural modes of coupled bending-torsion vibration of beams,” Transactions of the ASME, vol. 67, no. 1, pp. 61–67, 1945.

M. A. Prohl, “A general method for calculating critical speeds of flexible rotors,” Transactions of the ASME, Journal of Applied Mechanics, vol. 66, pp. A-142–A-148, 1945.

J. W. Lund, “Stability and damped critical speeds of a flexible rotor in fluid-film bearings,” Transactions of the ASME, Journal of Engineering for Industry Series B, vol. 96, no. 2, pp. 509–517, 1974.

B. T. Murphy and J. M. Vance, “An improved method for calculating critical speeds and rotor dynamic stability of turbo machinery,” Transactions of the ASME, Journal of Engineering for Power, vol. 105, no. 3, pp. 591–595, 1983.

R. Whalley and A. Abdul-Ameer, “Whirling prediction with geometrical shaft profiling,” Applied Mathematical Modelling, vol. 33, no. 7, pp. 3166–3177, 2009.

R. L. Ruhl and J. F. Booker, “A finite element model for distributed parameter turborotor systems,” Journal of Engineering for Industry, vol. 94, no. 1, pp. 126–132, 1972.

H. D. Nelson and J. M. McVaugh, “The dynamics of rotor-bearing systems using finite elements,” Transactions of the ASME, Journal of Engineering, vol. 98, no. 2, pp. 593–600, 1976.

H. D. Nelson, “A finite rotating shaft element using timoshenko beam theory,” Journal of Mechanical Design, vol. 102, no. 4, pp. 793–803, 1980.

L. M. Greenhill, W. B. Bickford, and H. D. Nelson, “A conical beam finite element for rotor dynamics analysis,” Journal of Vibration, Acoustics, Stress, and Reliability in Design, vol. 107, no. 4, pp. 421–430, 1985.

G. Genta and A. Gugliotta, “A conical element for finite element rotor dynamics,” Journal of Sound and Vibration, vol. 120, no. 1, pp. 175–182, 1988.

M. A. Mohiuddin and Y. A. Khulief, “Modal characteristics of cracked rotors using a conical shaft finite element,” Computer Methods in Applied Mechanics and Engineering, vol. 162, no. 1–4, pp. 223–247, 1998.

K. E. Rouch and J.-S. Kao, “A tapered beam finite element for rotor dynamics analysis,” Journal of Sound and Vibration, vol. 66, no. 1, pp. 119–140, 1979.