0
Research Papers

Numerical Tracking of Limit Points for Direct Parametric Analysis in Nonlinear Rotordynamics

[+] Author and Article Information
Lihan Xie

Université de Lyon,
CNRS,
INSA-Lyon,
LaMCoS, UMR5259,
Villeurbanne F-69621, France;
CEA-Saclay, DEN, DANS, DM2S, SEMT, DYN
Gif sur Yvette F-91191, France
e-mail: lihan.xie@insa-lyon.fr

Sébastien Baguet

Université de Lyon,
CNRS,
INSA-Lyon,
LaMCoS, UMR5259,
Villeurbanne F-69621, France
e-mail: sebastien.baguet@insa-lyon.fr

Benoit Prabel

CEA-Saclay, DEN, DANS, DM2S, SEMT, DYN
Gif sur Yvette F-91191, France;
IMSIA,
UMR9219, CNRS-EDF-CEA-ENSTA,
Université Paris Saclay,
Palaiseau F-91762, France
e-mail: benoit.prabel@cea.fr

Régis Dufour

Université de Lyon,
CNRS,
INSA-Lyon,
LaMCoS, UMR5259
Villeurbanne F-69621, France
e-mail: regis.dufour@insa-lyon.fr

1Corresponding author.

Contributed by the Technical Committee on Vibration and Sound of ASME for publication in the JOURNAL OF VIBRATION AND ACOUSTICS. Manuscript received June 10, 2015; final manuscript received November 13, 2015; published online January 20, 2016. Assoc. Editor: Izhak Bucher.

J. Vib. Acoust 138(2), 021007 (Jan 20, 2016) (9 pages) Paper No: VIB-15-1210; doi: 10.1115/1.4032182 History: Received June 10, 2015; Revised November 13, 2015

A frequency-domain approach for direct parametric analysis of limit points (LPs) of nonlinear dynamical systems is presented in this paper. Instead of computing responses curves for several values of a given system parameter, the direct tracking of LPs is performed. The whole numerical procedure is based on the harmonic balance method (HBM) and can be decomposed in three distinct steps. First, a response curve is calculated by HBM combined with a continuation technique until an LP is detected. Then this starting LP is used to initialize the direct tracking of LPs which is based on the combination of a so-called extended system and a continuation technique. With only one computation, a complete branch of LPs is obtained, which provides the stability boundary with respect to system parameters such as nonlinearity or excitation level. Several numerical examples demonstrate the capabilities and the performance of the proposed method.

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

References

Jacquet-Richardet, G. , Torkhani, M. , Cartraud, P. , Thouverez, F. , Baranger, T. N. , Herran, M. , Gibert, C. , Baguet, S. , Almeida, P. , and Peletan, L. , 2013, “ Rotor to Stator Contacts in Turbomachines. Review and Application,” Mech. Syst. Signal Process., 40(2), pp. 401–420. [CrossRef]
Duran, C. , Manin, L. , Andrianoely, M.-A. , Bordegaray, C. , Battle, F. , and Dufour, R. , 2015, “ Effect of Rotor-Stator Contact on the Mass Unbalance Response,” 9th IFToMM International Conference on Rotor Dynamics (Mechanisms and Machine Science), P. Pennacchi, ed., Vol. 21, Springer International Publishing, pp. 1965–1975.
Chen, C. , Dai, L. , and Fu, Y. , 2007, “ Nonlinear Response and Dynamic Stability of a Cracked Rotor,” Commun. Nonlinear Sci. Numer. Simul., 12(6), pp. 1023–1037. [CrossRef]
Cheng, L. , Li, N. , Chen, X. , and He, Z. , 2011, “ The Influence of Crack Breathing and Imbalance Orientation Angle on the Characteristics of the Critical Speed of a Cracked Rotor,” J. Sound Vib., 330(9), pp. 2031–2048. [CrossRef]
Khonsari, M. M. , and Chang, Y. J. , 1993, “ Stability Boundary of Non-Linear Orbits Within Clearance Circle of Journal Bearings,” ASME J. Vib. Acoust., 115(3), pp. 303–307. [CrossRef]
Brown, R. D. , Addison, P. , and Chan, A. , 1994, “ Chaos in the Unbalance Response of Journal Bearings,” Nonlinear Dyn., 5(4), pp. 421–432. [CrossRef]
Dakel, M.-Z. , Baguet, S. , and Dufour, R. , 2014, “ Nonlinear Dynamics of a Support-Excited Flexible Rotor With Hydrodynamic Journal Bearings,” J. Sound Vib., 333(10), pp. 2774–2799. [CrossRef]
Sundararajan, P. , and Noah, S. , 1997, “ Dynamics of Forced Nonlinear Systems Using Shooting/Arc-Length Continuation Method—Application to Rotor Systems,” ASME J. Vib. Acoust., 119(1), pp. 9–19. [CrossRef]
Kim, Y.-B. , and Noah, S. , 1996, “ Quasi-Periodic Response and Stability Analysis for a Non-Linear Jeffcott Rotor,” J. Sound Vib., 190(2), pp. 239–253. [CrossRef]
Von Groll, G. , and Ewins, D. J. , 2001, “ The Harmonic Balance Method With Arc-Length Continuation in Rotor/Stator Contact Problems,” J. Sound Vib., 241(2), pp. 223–233. [CrossRef]
Nayfeh, A. H. , and Mook, D. T. , 1995, Nonlinear Oscillations, Wiley, New York.
Berlioz, A. , Dufour, R. , and Sinha, S. , 2000, “ Bifurcation in a Nonlinear Autoparametric System Using Experimental and Numerical Investigations,” Nonlinear Dyn., 23(2), pp. 175–187. [CrossRef]
Seydel, R. , 2009, Practical Bifurcation and Stability Analysis (Interdisciplinary Applied Mathematics), Springer, New York.
Moore, G. , and Spence, A. , 1980, “ The Calculation of Turning Points of Nonlinear Equations,” SIAM J. Numer. Anal., 17(4), pp. 567–576. [CrossRef]
Wriggers, P. , and Simo, J. C. , 1990, “ A General Procedure for the Direct Computation of Turning and Bifurcation Points,” Int. J. Numer. Methods Eng., 30(1), pp. 155–176.
Battini, J.-M. , Pacoste, C. , and Eriksson, A. , 2003, “ Improved Minimal Augmentation Procedure for the Direct Computation of Critical Points,” Comput. Methods Appl. Mech. Eng., 192(16–18), pp. 2169–2185. [CrossRef]
Petrov, E. P. , 2007, “ Direct Parametric Analysis of Resonance Regimes for Nonlinear Vibrations of Bladed Disks,” ASME J. Turbomach., 129(3), pp. 495–502. [CrossRef]
Petrov, E. P. , 2009, “ Method for Sensitivity Analysis of Resonance Forced Response of Bladed Disks With Nonlinear Contact Interfaces,” ASME J. Eng. Gas Turbines Power, 131(2), p. 022510.
Liao, H. , and Wang, J. , 2013, “ Maximization of the Vibration Amplitude and Bifurcation Analysis of Nonlinear Systems Using the Constrained Optimization Shooting Method,” J. Sound Vib., 332(16), pp. 3781–3793. [CrossRef]
Baguet, S. , and Cochelin, B. , 2002, “ Stability of Thin-Shell Structures and Imperfection Sensitivity Analysis With the Asymptotic Numerical Method,” Rev. Eur. Elém. Finis, 11(2–4), pp. 493–509.
Cameron, T. M. , and Griffin, J. H. , 1989, “ An Alternating Frequency/Time Domain Method for Calculating the Steady-State Response of Nonlinear Dynamic Systems,” ASME J. Appl. Mech., 56(1), pp. 149–154. [CrossRef]
Crisfield, M. , 1981, “ A Fast Incremental/Iterative Solution Procedure That Handles Snap-Through,” Comput. Struct., 13(1), pp. 55–62. [CrossRef]
Peletan, L. , Baguet, S. , Torkhani, M. , and Jacquet-Richardet, G. , 2013, “ A Comparison of Stability Computational Methods for Periodic Solution of Nonlinear Problems With Application to Rotordynamics,” Nonlinear Dyn., 72(3), pp. 671–682. [CrossRef]
Keller, H. B. , 1977, “ Numerical Solution of Bifurcation and Nonlinear Eigenvalue Problems,” Applications of Bifurcation Theory, P. H. Rabinowitz , ed., Academic Press, New York, pp. 359–384.
Felippa, C. A. , 1987, “ Traversing Critical Points by Penalty Springs,” NUMETA'87 Conference, Nijhoff Pubs, Dordrecht, Holland.
Jiang, J. , 2009, “ Determination of the Global Responses Characteristics of a Piecewise Smooth Dynamical System With Contact,” Nonlinear Dyn., 57(3), pp. 351–361. [CrossRef]
Peletan, L. , Baguet, S. , Torkhani, M. , and Jacquet-Richardet, G. , 2014, “ Quasi-Periodic Harmonic Balance Method for Rubbing Self-Induced Vibrations in Rotor-Stator Dynamics,” Nonlinear Dyn., 78(4), pp. 2501–2515. [CrossRef]
Lalanne, M. , and Ferraris, G. , 1998, Rotordynamics Prediction in Engineering, 2nd ed., Wiley, New York.
CAST3M, 2015, “ Finite Element Software,” CEA (French Atomic Energy Commission), http://www-cast3m.cea.fr/

Figures

Grahic Jump Location
Fig. 1

Frequency response of the Duffing system for α = 0.02

Grahic Jump Location
Fig. 2

Frequency response of the Duffing system for α = 10

Grahic Jump Location
Fig. 3

LPs tracking for α ∈ [0 10]

Grahic Jump Location
Fig. 4

Zoom of LPs tracking

Grahic Jump Location
Fig. 5

Projections of LPs tracking

Grahic Jump Location
Fig. 6

Forced response of the Jeffcott rotor for μ = 0

Grahic Jump Location
Fig. 7

Forced response of the Jeffcott rotor for μ = 0.11

Grahic Jump Location
Fig. 8

Forced response of the Jeffcott rotor for μ = 0.2

Grahic Jump Location
Fig. 9

Forced response of the Jeffcott rotor for μ = 0.11, comparison with time integration for increasing and decreasing frequency sweeps

Grahic Jump Location
Fig. 10

LPs tracking of the Jeffcott rotor with varying friction coefficient μ

Grahic Jump Location
Fig. 11

Multi-DOFs FE rotor [23]

Grahic Jump Location
Fig. 12

FE rotor: forced response at node 6 for μ = 0

Grahic Jump Location
Fig. 13

FE rotor: forced response at node 6 for μ = 0.03

Grahic Jump Location
Fig. 14

FE rotor: forced response at node 6 for μ = 0.03, comparison with time integration for increasing sweep

Grahic Jump Location
Fig. 15

FE rotor: LPs tracking with varying friction coefficientμ

Tables

Errata

Discussions

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