Research Papers

An Efficient Approach for the Frequency Analysis of Nonaxisymmetric Rotating Structures: Application to a Coupled Bladed Birotor System

[+] Author and Article Information
Cécile Dumartineix

École Centrale de Lyon,
Laboratoire de Tribologie et Dynamique des
UMR-CNRS 5513,
Écully cedex 69134, France;
Safran Aircraft Engines,
Rond Point René Ravaud—Réau,
Moissy-Cramayel 77550, France
e-mail: cecile.dumartineix@doctorant.ec-lyon.fr

Benjamin Chouvion, Fabrice Thouverez

École Centrale de Lyon,
Laboratoire de Tribologie et Dynamique des
UMR-CNRS 5513,
Écully cedex 69134, France

Marie-Océane Parent

Safran Aircraft Engines,
Rond Point René Ravaud—Réau,
Moissy-Cramayel 77550, France

1Corresponding author.

Manuscript received June 26, 2018; final manuscript received July 4, 2018; published online December 10, 2018. Editor: Jerzy T. Sawicki.

J. Eng. Gas Turbines Power 141(4), 041033 (Dec 10, 2018) (8 pages) Paper No: GTP-18-1371; doi: 10.1115/1.4040865 History: Received June 26, 2018; Revised July 04, 2018

The improvement of efficiency in the design of turbomachines requires a reliable prediction of the vibrating behavior of the whole structure. The simulation of blades vibrations is decisive and this is usually based on elaborated finite element model restricted to the bladed-disk. However, the blades dynamic behavior can be strongly affected by interactions with other parts of the engine. Global dynamic studies that consider these other parts are required but usually come with a high numerical cost. In the case of a birotor architecture, two coaxial rotors with different rotating speeds can be coupled with a bearing system. The mechanical coupling between the shafts generates energy exchange that alters the dynamic behavior of the blades. The equations of motion of the whole structure that take into account the coupling contain periodic time-dependent coefficients due to the difference of rotational speed between both rotors. Equations of this kind, with variable coefficients, are typically difficult to solve. This study presents a preprocessing method to guarantee the elimination of time-dependent coefficients in the birotor equations of motion. This method is tested with a simplified finite element model of two bladed-disks coupled with linear stiffnesses. We obtain accurate results when comparing frequency analysis of preprocessed equations with time-integration resolution of the initial set of equations. The developed methodology also offers a substantial time saving.

Copyright © 2019 by ASME
Your Session has timed out. Please sign back in to continue.


Legrand, M. , Pierre, C. , Cartraud, P. , and Lombard, J.-P. , 2009, “ Two-Dimensional Modeling of an Aircraft Engine Structural Bladed Disk-Casing Modal Interaction,” J. Sound Vib., 319(1)pp. 366–391. [CrossRef]
Thomas, D. L. , 1979, “ Dynamics of Rotationally Periodic Structures,” Int. J. Numer. Methods Eng., 14(1), pp. 81–102.
Laxalde, D. , Thouverez, F. , and Lombard, J. P. , 2007, “ Dynamical Analysis of Multi-Stage Cyclic Structures,” Mech. Res. Commun., 34(4), pp. 379–384. [CrossRef]
Battiato, G. , Firrone, C. M. , Berruti, T. M. , and Epureanu, B. I. , 2016, “ Reduced Order Modeling for Multi-Stage Coupling of Cyclic Symmetric Structures,” Proc. ISMA, pp. 3447–3462.
Ewins, D. J. , 2000, Modal Testing: Theory, Practice and Application, 2nd ed., Research Studies Press Wiley, New York.
Shapiro, V. E. , 2001, “ Rotating Class of Parametric Resonance Processes in Coupled Oscillators,” Phys. Lett. A, 290(5–6), pp. 288–296. [CrossRef]
Pradeep, S. , and Shrivastava, S. K. , 1988, “ On the Stability of the Damped Mathieu Equation,” Mech. Res. Commun., 15(6), pp. 353–359. [CrossRef]
Mond, M. , Cederbaum, G. , Khan, P. B. , and Zarmi, Y. , 1993, “ Stability Analysis of the Non-Linear Mathieu Equation,” J. Sound Vib., 167(1), pp. 77–89. [CrossRef]
Lalanne, M. , and Ferraris, G. , 2000, Rotordynamics Prediction in Engineering, 2nd ed., Wiley, Chichester, UK.
Gruin, M. , Thouverez, F. , Blanc, L. , and Jean, P. , 2011, “ Nonlinear Dynamics of a Bladed Dual-Shaft,” Eur. J. Comput. Mech., 20(1–4), pp. 207–225. [CrossRef]
Salvat, N. , Batailly, A. , and Legrand, M. , 2015, “ Two-Dimensional Modeling of Unilateral Contact-Induced Shaft Precessional Motions in Bladed-Disk/Casing Systems,” Int. J. Non-Linear Mech., 78, pp. 90–104. [CrossRef]
Gmür, T. , 1997, Dynamique Des Structures: Analyse Modale Numérique, Presses polytechniques et universitaires romandes, Lausanne, Switzerland.
Timoshenko, S. , 1937, Vibration Problems in Engineering, 2nd ed., D.Van Nostrand Company INC, New York.
Rao, S. S. , 2007, Vibration of Continuous Systems, Wiley, New York.
Lazarus, A. , Prabel, B. , and Combescure, D. , 2010, “ A 3D Finite Element Model for the Vibration Analysis of Asymmetric Rotating Machines,” J. Sound Vib., 329(18), pp. 3780–3797. [CrossRef]
Guskov, M. , Sinou, J. J. , and Thouverez, F. , 2008, “ Multi-Dimensional Harmonic Balance Applied to Rotor Dynamics,” Mech. Res. Commun., 35(8), pp. 537–545. [CrossRef]
Zhou, B. , Thouverez, F. , and Lenoir, D. , 2015, “ A Variable-Coefficient Harmonic Balance Method for the Prediction of Quasi-Periodic Response in Nonlinear Systems,” Mech. Syst. Signal Process., 64–65 pp. 233–244. [CrossRef]
Feldmann, P. , and Roychowdhury, J. , 1996, “ Computation of Circuit Waveform Envelopes Using an Efficient, Matrix-Decomposed Harmonic Balance Algorithm,” International Conference on Computer Aided Design, San Jose, CA, Nov. 10–14, pp. 295–300.
Laxalde, D. , 2007, “ Étude D'amortisseurs Non-Linéaires Appliqués Aux Roues Aubagées Et Aux Systèmes Multi-Étages,” Ph. D. Thèse de doctorat, Ecole Centrale de Lyon, Écully, France.
Kundert, K. S. , Sorkin, G. B. , and Sangiovanni-Vincentelli, A. , 1988, “ Applying Harmonic Balance to Almost-Periodic Circuits,” IEEE Trans. Microwave Theory Tech., (2), pp. 366–378.
Guédeney, T. , Gomar, A. , Gallard, F. , Sicot, F. , Dufour, G. , and Puigt, G. , 2013, “ Non-Uniform Time Sampling for Multiple-Frequency Harmonic Balance Computations,” J. Comput. Phys., 236, pp. 317–345. [CrossRef]


Grahic Jump Location
Fig. 1

Sectional view of a two-spool turbofan. Low pressure rotor () and high pressure rotor (), stator ().

Grahic Jump Location
Fig. 2

Studied birotor system: (a) simplified model of the bladed birotor system and (b) bearing model

Grahic Jump Location
Fig. 3

Fixing area of the bearing (circular line) on one bladed-disk for our studied case

Grahic Jump Location
Fig. 4

Finite element model of one sector of each bladed-disk of the studied birotor system

Grahic Jump Location
Fig. 5

Modal deformation of the first 2-diameter mode of the rotor R1

Grahic Jump Location
Fig. 6

Time radial displacements of a node of each rotor in the bearing area. Time integration of the initial equation of motion (), frequency resolution of the preprocessed equation (): (a) radial displacements of rotor R1 over time. Zoom on the stationary state () and (b) radial displacements of rotor R2 over time.

Grahic Jump Location
Fig. 7

Maximum radial displacements of a node of each rotor in the bearing area computed for a range of rotation speeds with a corotative configuration: (a) radial displacements of rotor R1 (left) and rotor R2 (right) at the bearing and (b) zoom on the crossing of resonance frequencies of both rotors



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In