International Journal of Rotating Machinery

VolumeÂ 2008, Article IDÂ 752062, 10 pages

http://dx.doi.org/10.1155/2008/752062

## Free Vibration Analysis of a Rotating Composite Shaft Using the -Version of the Finite Element Method

Department of Mechanical Engineering, Faculty of Sciences Engineering, Abou Bakr Belkaid University, Tlemcen 13000, Algeria

Received 19 November 2007; Accepted 24 May 2008

Academic Editor: AgnesÂ Muszynska

Copyright Â© 2008 A. Boukhalfa et al. 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

This paper is concerned with the dynamic behavior of the rotating composite shaft on rigid bearings. A -version, hierarchical finite element is employed to define the model. A theoretical study allows the establishment of the kinetic energy and the strain energy of the shaft, necessary to the result of the equations of motion. In this model the transverse shear deformation, rotary inertia and gyroscopic effects, as well as the coupling effect due to the lamination of composite layers have been incorporated. A hierarchical beam finite element with six degrees of freedom per node is developed and used to find the natural frequencies of a rotating composite shaft. A program is elaborate for the calculation of the eigenfrequencies and critical speeds of a rotating composite shaft. To verify the present model, the critical speeds of composite shaft systems are compared with those available in the literature. The efficiency and accuracy of the methods employed are discussed.

#### 1. Introduction

The application of composite shafts has come a long way from early low-speed automotive driveshafts to helicopter tail rotors operating above the second critical speed. With operation at supercritical speeds, a substantial amount of payoffs and net system weight reductions are possible. At the same time, the rotordynamic aspects assume more importance, and detailed analysis is required. There are some technological problems associated with implementation, such as joints with bearings, affixing of lumped masses, couplings, provision of external damping, and so forth. The solutions proposed are just adequate, but require substantial refinements, which might explain some of the differing experiences of various authors.

Zinberg and Symonds [1] described a boron/epoxy composite tail rotor driveshaft for a helicopter. The critical speeds were determined using equivalent modulus beam theory (EMBT), assuming the shaft to be a thin-walled circular tube simply supported at the ends. Shear deformation was not taken into account. The shaft critical speed was determined by extrapolation of the unbalance response curve which was obtained in the subcritical region.

Dos Reis et al. [2] published analytical investigations on thin-walled layered composite cylindrical tubes. In part III of the series of publications, the beam element was extended to formulate the problem of a rotor supported on general eight coefficient bearings. Results were obtained for shaft configuration of Zinberg and Symmonds. The authors have shown that bending-stretching coupling and shear-normal coupling effects change with stacking sequence, and alter the frequency values. Gupta and Singh [3] studied the effect of shear-normal coupling on rotor natural frequencies and modal damping. Kim and Bert [4] have formulated the problem of determination of critical speeds of a composite shaft including the effects of bending-twisting coupling. The shaft was modelled as a Bresse-Timoshenko beam. The shaft gyroscopics have also been included. The results compare well with Zinberg's rotor [1]. In another study, Bert and Kim [5] have analysed the dynamic instability of a composite drive shaft subjected to fluctuating torque and/or rotational speed by using various thin shell theories. The rotational effects include centrifugal and Coriolis forces. Dynamic instability regions for a long span simply supported shaft are presented.

M.-Y. Chang et al. [6] published the vibration behaviours of the rotating composite shafts. In the model, the transverse shear deformation, rotary inertia, and gyroscopic effects, as well as the coupling effect due to the lamination of composite layers have been incorporated. The model based on a first-order shear deformable beam theory (continuum-based Timoshenko beam theory).

C.-Y. Chang et al. [7] published the vibration analysis of rotating composite shafts containing randomly oriented reinforcements. The Mori-Tanaka mean-field theory is adopted here to account for the interaction at the finite concentrations of reinforcements in the composite material.

Additional recent work on composite shafts dealing with both the theoretical and experimental aspects was reported by Singh [8], Gupta and Singh [3], and Singh and Gupta [9, 10]. Rotordynamic formulation based on equivalent modulus beam theory was developed for a composite rotor with a number of lumped masses and supported on general eight coefficient bearings. A layerwise beam theory was derived by Gupta and Singh [3] from an available shell theory, with a layerwise displacement field, and was then extended to solve a general composite rotordynamic problem. The conventional rotor dynamic parameters as well as critical speeds, natural frequencies, damping factors, unbalance response, and threshold of stability were analysed in detail and results from the formulations based on the two theories, namely, the equivalent modulus beam theory (EMBT) and layerwise beam theory (LBT) were compared [9, 10]. The experimental rotordynamic studies carried by Singh and Gupta [11, 12] were conducted on two filament wound carbon/epoxy shafts with constant winding angles ( and ). Progressive balancing had to be carried out to enable the shaft to traverse through the first critical speed. Inspite of the very different shaft configurations used, the authors' have shown that bending-stretching coupling and shear-normal coupling effects change with stacking sequence, and alter the frequency values.

Some practical aspects such as effect of shaft disc angular misalignment, interaction between shaft bow, which is common in composite shafts and rotor unbalance, and an unsuccessful operation of a composite rotor with an external damper were discussed and reported by Singh and Gupta [11]. The Bode and cascade plots were generated and orbital analysis at various operating speeds was performed. The experimental critical speeds showed good correlation with the theoretical prediction.

This paper deals with the -version, hierarchical finite element method applied to free vibration analysis of rotating composite shafts. The hierarchical concept for finite element shape functions has been investigated during the past 25 years. Babuska et al. [13] established a theoretical basis for -elements, where the mesh keeps unchanged and the polynomial degree of the shape functions is increased; however, in the standard -version of the finite element method the mesh is refined to achieve convergence and the polynomial degree of the shape functions remains unchanged. Since then, standard forms of the hierarchical shape functions have been represented in the literature elsewhere; see for instance [14, 15].

Meirovitch and Baruh [16] and Zhu [17] have shown that the hierarchical finite element method yields a better accuracy than the -version for eigenvalues problems. The hierarchical shape functions used by Bardell [18] are based on integrated Legendre orthogonal polynomials; the symbolic computing is used to calculate the mass and stiffness matrices of beams and plates. CĂ´tĂ© and Charron [19] give the selection of -version shape functions for plate vibration analysis.

In the presented composite shaft model, the Timoshenko theory will be adopted. It is the purpose of the present work to study dynamic characteristics such as natural frequencies, whirling frequencies, and the critical speeds of the rotating composite shaft. In the model, the transverse shear deformation, rotary inertia, and gyroscopic effects, as well as the coupling effect due to the lamination of composite layers have been incorporated. To determine the rotating shaft system's responses, the hierarchical finite element method with trigonometric shape functions [20, 21] is used here to approximate the governing equations by a system of ordinary differential equations.

#### 2. Equations of Motion

#### 2.1. Kinetic and Strain Energy Expressions

The shaft is modelled as a Timoshenko beam, that is, first-order shear deformation theory with rotary inertia and gyroscopic effect is used. The shaft rotates at constant speed about its longitudinal axis. Due to the presence of fibers oriented than axially or circumferentially, coupling is made between bending and twisting. The shaft has a uniform, circular cross-section.

The following displacement field of a rotating shaft is assumed by choosing the coordinate axis to coincide with the shaft axis: where , , and are the flexural displacements of any point on the cross-section of the shaft in the , , and directions (see Figure 1), the variables , , and are the flexural displacements of the shaft's axis, while and are the rotation angles of the cross-section, about the and axis, respectively. The is the angular displacement of the cross-section due to the torsion deformation of the shaft.

The strain components in the cylindrical coordinate system (as shown in Figures 1-2) can be written in terms of the displacement variables defined earlier as

Let us consider a composite shaft consists of layered (Figure 2) of fiber inclusion reinforced laminate. The constitutive relations for each layer are described by where are the effective elastic constants, they are related to lamination angle (as shown in Figures 3-4) and the elastic constants of principal axes.

The stress-strain relations of the th layer expressed in the cylindrical coordinate system (Figure 5) can be expressed as where is the transverse shear correction factor.

The formula of the strain energy is

The various components of strain energy come from the shaft: where and are given in the appendix.

The kinetic energy of the rotating composite shaft, including the effects of translatory and rotary inertia, can be written as where is the rotating speed of the shaft which is assumed constant, is the length of the shaft, the term accounts for the gyroscopic effect, and represent the rotary inertia effect. The mass moments of inertia , the diametrical mass moments of inertia , and polar mass moment of inertia of rotating shaft per unit length are defined in the appendix. As the term is far smaller than , it will neglected in further analysis.

#### 2.2. Hierarchical Beam Element Formulation

The spinning flexible beam is descretised into one hierarchical finite element is shown in Figure 6. The element's nodal d.o.f. at each node are , , , , , and . The local and nondimensional coordinates are related by The vector displacement formed by the variables , , , , , and can be written as And it can be expressed as where [] is the matrix of the shape functions, given by where , , , , , and are the numbers of hierarchical terms of displacements (are the numbers of shape functions of displacements).

In this work,

The vector of generalised coordinates given by where

The group of the shape functions used in this study is given by

The functions (, ) are those of the finite element method necessary to describe the nodal displacements of the element, whereas the trigonometric functions contribute only to the internal field of displacement and do not affect nodal displacements. The most attractive particularity of the trigonometric functions is that they offer great numerical stability. The shaft is modeled by only one element called hierarchical finite element.

By applying the Euler-Lagrange equations, the equations of motion of free vibration of spinning flexible shaft can be obtained:

The system obtained is linear of which equations are coupled once by the gyroscopic effect, represented by the matrix [].

[] and [] are the conventional hierarchical finite element mass and stiffness matrix, [] is the gyroscopic matrix (the different matrices of the system of equation are given in the appendix).

#### 3. Results

A program based on the formulation proposed was developed for the resolution of (16).

#### 3.1. Convergence

First, the material properties for boron-epoxy are listed in Table 1, and the geometric parameters are , , , 10 layers of equal thickness (, , , , ). The shear correction factor and the rotating speed .

The results of the five bending modes for various boundary conditions of the composite shaft as a function of the number of hierarchical terms are shown in Figures 7 and 8. Figures clearly show that rapid convergence from above to the exact values occur as the number of hierarchical terms is increased.

#### 3.2. Numerical Examples and Discussions

In all the following examples, .

In the first example, the properties of the boron-epoxy composite shaft are given by Table 1. The results obtained using the present model are shown in Table 2 together with those of referenced papers. As can be seen from the table, our results are close to those predicted by other beam theories. Since in the studied example the wall of the shaft is relatively thin, models based on shell theories [4] are expected to yield more accurate results. In the present example, the critical speed measured from the experiment however is still underestimated by using Sander shell theory, while overestimated by Donnell shallow shell theory. When the material of the shaft is changed to the graphite-epoxy given in Table 1 with other conditions left unchanged, the critical speed obtained from the present model is shown in Table 3. In this case, the result from the present model is compatible to that of the Bresse-Timoshenko beam theory of M.-Y. Chang et al. [6].

Next, comparisons are made with those of [6, 22] for different length to mean diameter ratios . The shafts being analysed are made of the graphite-epoxy material given in Table 1 and all have the same lamination [////]. The mean diameter and the wall thickness of the shaft remain the same as the previous examples. The shear correction factor being used is again 0.503. The results are listed in Table 4. Further results being compared are for generalized orthotropic composite tube of different lamination angles . The results are shown in Table 5.

In our work, the shaft is modelled by only one element with two nodes, but in the model of [6] the shaft is modelled by 20 finite elements of equal length (-version). The rapid convergence while taking only one element and a reduced number of shape functions shows the advantage of the method used.

From the thin-walled shaft systems studied above, the present shaft model yields result in all cases close to those of the model of Bert and Kim [22] based on the Bresse-Timoshenko theory. One should stress here that the present model is not only applicable to the thin-walled composite shafts as studied above, but also to the thick-walled shafts as well as to the solid ones.

In
the following example, the frequencies responses of a graphite-epoxy composite
shaft system are analysed. The material properties are those listed in Table 1.
The lamination scheme of the shaft remains the same as previous examples, while
its geometric properties, the Campbell
diagram containing the frequencies of the first five pairs of bending whirling
modes of the above composite system is shown in Figure 9. The intersection
point of the line (Î© = *Ď‰*) with the whirling frequency curves indicate the speed
at which the shaft will vibrate violently (i.e., the critical speed). The first
10 eigenvalues correspond to 5 forward (F) and 5 backward (B) whirling bending
modes of the shaft.

The gyroscopic effect causes a coupling of orthogonal displacements to the axis of rotation, and by consequence separates the frequencies in two branches: backward precession mode and forward precession mode. In all cases, the backward modes increase with increasing rotating speed however the forward modes decrease.

Figure 10 shows the variation of the bending fundamental frequency *Ď‰* as a function of the rotating speed Î© (Campbell diagram) for different boundary conditions and different rotating speeds.

The gyroscopic effect inherent to rotating structures induces a precession motion. The backward precession modes (1B) increase with increasing the rotating speed, however the forward precession (1F) modes decrease.

#### 4. Conclusion

This paper has presented the free vibration analysis of spinning composite shaft using the -version, hierarchical finite element method with trigonometric shape functions. Results obtained using the method has been evaluated against those available in the literature and the agreement has been found to be good. The main conclusions have emerged from this work, these are itemised below.(1) Monotonic and uniform convergence is found to occur as the number of hierarchical functions is increased.(2)The dynamic characteristics of rotating composite shaft are influenced significantly by varying the fiber orientation, the rotating speed, and boundary conditions.(3)The gyroscopic effect causes a coupling of orthogonal displacements to the axis of rotation, and by consequence separates the frequencies in two branches: backward and forward precession modes. In all cases, the backward modes increase with increasing rotating speed however the forward modes decrease.

#### Appendix

The terms , of (6) and , , of (7) are given as follows: where is the number of the layer, is the th layer inner radius of the composite shaft, and is the th layer outer of the composite shaft. is the length of the composite shaft and is the density of the th layer of the composite shaft.

Whereas the various matrices of (16) are expressed as follows: where , with ().

The terms of the matrices are a function of the integrals, () indicate the number of the shape functions used, where is the order of derivation.

*Nomenclature*

: | Displacement in direction |

: | Displacement in direction |

: | Displacement in direction |

: | Rotation angles of the cross-section on axis |

: | Rotation angles of the cross-section on axis |

: | Angular displacement of the cross-section due to the torsion deformation of the shaft |

: | Cartesian coordinates |

Axes of the Cartesian coordinates | |

: | Cylindrical coordinates |

Axes of the cylindrical coordinates | |

(1, 2, 3): | Principal axes of a layer of laminate |

: | Elastic constants |

: | Young modulus |

: | Shear modulus |

: | Poisson coefficient |

: | Masse density |

: | Shear correction factor |

: | Length of the shaft |

: | Mean radius of the shaft |

: | Wall thickness of the shaft |

: | Number of the layer of the composite shaft |

: | The th layer inner radius of the composite shaft |

: | The th layer outer radius of the composite shaft |

: | Inner radius of the composite shaft |

: | Outer radius of the composite shaft |

: | Lamination angle |

: | Circumferential coordinate |

: | Local and nondimensional coordinates |

: | Strain tensor |

: | Stress tensor |

: | Frequency, eigenvalue |

: | Rotating speed |

: | Shape functions |

: | Matrix of the shape functions |

: | Number of the shape functions or number of hierarchical terms |

: | Time |

: | Kinetic energy |

: | Strain energy |

: | Generalised coordinates, with () |

: | Masse matrix |

: | Stiffness matrix |

: | Gyroscopic matrix. |

#### References

- H. Zinberg and M. F. Symonds, â€śThe development of an advanced composite tail rotor driveshaft,â€ť in
*Proceedings of the 26th Annual Forum of the American Helicopter Society*, Washington, DC, USA, June 1970. - H. L. M. dos Reis, R. B. Goldman, and P. H. Verstrate, â€śThin-walled laminated composite cylindrical tubes—part III: critical speed analysis,â€ť
*Journal of Composites Technology and Research*, vol. 9, no. 2, pp. 58â€“62, 1987. View at Google Scholar - K. Gupta and S. E. Singh, â€śDynamics of composite rotors,â€ť in
*Proceedings of the Indo-US Symposium on Emerging Trends in Vibration and Noise Engineering*, pp. 59â€“70, New Delhi, India, March 1996. - C. D. Kim and C. W. Bert, â€śCritical speed analysis of laminated composite, hollow drive shafts,â€ť
*Composites Engineering*, vol. 3, pp. 633â€“643, 1993. View at Google Scholar - C. W. Bert and C. D. Kim, â€śDynamic instability of composite-material drive shaft subjected to
fluctuating torque and/or rotational speed,â€ť
*Dynamics and Stability of Systems*, vol. 10, no. 2, pp. 125â€“147, 1995. View at Google Scholar - M.-Y. Chang, J.-K. Chen, and C.-Y. Chang, â€śA simple spinning laminated composite shaft model,â€ť
*International Journal of Solids and Structures*, vol. 41, no. 3-4, pp. 637â€“662, 2004. View at Publisher Â· View at Google Scholar - C.-Y. Chang, M.-Y. Chang, and J. H. Huang, â€śVibration analysis of rotating composite shafts containing randomly oriented reinforcements,â€ť
*Composite Structures*, vol. 63, no. 1, pp. 21â€“32, 2004. View at Publisher Â· View at Google Scholar - S. P. Singh,
*Some studies on dynamics of composite shafts*, Ph.D. thesis, Mechanical Engineering Department, IIT, Delhi, India, 1992. - S. P. Singh and K. Gupta, â€śDynamic analysis of composite rotors,â€ť in
*Proceedings of the 5th International Symposium on Rotating Machinery (ISROMAC '94)*, pp. 204â€“219, Kaanapali, Hawaii, USA, May 1994. - S. P. Singh and K. Gupta, â€śDynamic analysis of composite rotors,â€ť
*International Journal of Rotating Machinery*, vol. 2, no. 3, pp. 179â€“186, 1996. View at Publisher Â· View at Google Scholar - S. E. Singh and K. Gupta, â€śExperimental studies on composite shafts,â€ť in
*Proceedings of the International Conference on Advances in Mechanical Engineering*, pp. 1205â€“1221, Bangalore, India, December 1995. - S. P. Singh and K. Gupta, â€śComposite shaft rotordynamic analysis using a layerwise theory,â€ť
*Journal of Sound and Vibration*, vol. 191, no. 5, pp. 739â€“756, 1996. View at Publisher Â· View at Google Scholar - I. Babuska, B. A. Szabó, and I. N. Katz, â€śThe $p$-version of the finite element method,â€ť
*SIAM Journal on Numerical Analysis*, vol. 18, no. 3, pp. 515â€“545, 1981. View at Publisher Â· View at Google Scholar - B. A. Szabó and G. J. Sahrmann, â€śHierarchic plate and shell models based on $p$-extension,â€ť
*International Journal for Numerical Methods in Engineering*, vol. 26, no. 8, pp. 1855â€“1881, 1988. View at Publisher Â· View at Google Scholar - B. A. Szabó and I. Babuska,
*Finite Element Analysis*, John Wiley & Sons, New York, NY, USA, 1991. - L. Meirovitch and H. Baruh, â€śOn the inclusion principle for the hierarchical finite element method,â€ť
*International Journal for Numerical Methods in Engineering*, vol. 19, no. 2, pp. 281â€“291, 1983. View at Publisher Â· View at Google Scholar - D. C. Zhu, â€śDevelopment of hierarchical finite element method at BIAA,â€ť in
*Proceedings of the International Conference on Computational Mechanics*, pp. 123â€“128, Tokyo, Japan, May 1986. - N. S. Bardell, â€śThe application of symbolic computing to the hierarchical finite element method,â€ť
*International Journal for Numerical Methods in Engineering*, vol. 28, no. 5, pp. 1181â€“1204, 1989. View at Publisher Â· View at Google Scholar Â· View at MathSciNet - A. Côté and F. Charron, â€śOn the selection of $p$-version shape functions for plate vibration problems,â€ť
*Computers & Structures*, vol. 79, no. 1, pp. 119â€“130, 2001. View at Publisher Â· View at Google Scholar - A. Houmat, â€śA sector fourier $p$-element applied to free vibration analysis of sectorial plates,â€ť
*Journal of Sound and Vibration*, vol. 243, no. 2, pp. 269â€“282, 2001. View at Publisher Â· View at Google Scholar - S. M. Hamza Cherif, â€śFree vibration analysis of rotating flexible beams by using the fourier $p$- version of the finite element,â€ť
*International Journal of Computational Methods*, vol. 2, no. 2, pp. 255â€“269, 2005. View at Publisher Â· View at Google Scholar - C. W. Bert and C.-D. Kim, â€śWhirling of composite-material driveshafts including bending-twisting
coupling and transverse shear deformation,â€ť
*Journal of Vibration and Acoustics*, vol. 117, no. 1, pp. 17â€“21, 1995. View at Publisher Â· View at Google Scholar - C. W. Bert, â€śThe effect of bending-twisting coupling on the critical speed of a driveshafts,â€ť in
*Proceedings of the 6th Japan-US Conference on Composites Materials*, pp. 29â€“36, Technomic, Orlando, Fla, USA, June 1992.