Fabrice Thouverez

A dynamic Lagrangian frequency-time method for the vibration of dry-friction-damped systems

Abstract A new frequency–time domain procedure, the dynamic Lagrangian mixed frequency–time method (DLFT), is proposed to calculate the non-linear steady state response to periodic excitation of structural systems subject to dry friction damping. In this formulation, the dynamic Lagrangians are defined as the non-linear contact forces obtained from the equations of motion in the frequency domain, with the adjunction of a penalization on the difference between the interface displacements calculate by the non-linear solver in the frequency domain and those calculated in the time domain from the non-linear contact forces, thus accounting for Coulomb friction and non-penetration conditions. The dynamic Lagrangians allow one to solve for the non-linear forces between two points in contact without using artifacts such as springs. The new DLFT method is thus particularly well suited to handling finite element models of structures in frictional contact, as it does not require a special model for the contact interface. Dynamic Lagrangians are also better suited to frequency-domain friction problems than the traditional time-domain method of augmented Lagrangians. Furthermore, a reduction of the non-linear system to relative interface displacements is introduced to decrease the computation time. The DLFT method is validated for a beam in contact with a flexible dry friction element connected to ground, for frictional constraints that feature two-dimensional relative motion. Results are also obtained for a large-scale structural system with a large number of one-dimensional dry-friction dampers. The DLFT method is shown to be accurate and fast, and it does not suffer from convergence problems, at least in the examples studied.

Qualitative analysis of forced response of blisks with friction ring dampers

A damping strategy for blisks (integrally bladed disks) of turbomachinery involving a friction ring is investigated. These rings, located in grooves underside the wheel of the blisks, are held in contact by centrifugal loads and the energy is dissipated when relative motions between the ring and the disk occur. A representative lumped parameter model of the system is introduced and the steady-state nonlinear response is derived using a multi-harmonic balance method combined with an AFT procedure where the friction force is calculated in the time domain. Numerical simulations are presented for several damper characteristics and several excitation configurations. From these results, the performance of this damping strategy is discussed and some design guidelines are given.

Numerical and Experimental Study of Friction Damping in Blade Attachments of Rotating Bladed Disks

In order to mitigate high cycle fatigue risks in bladed disks, the prediction of the vibration levels early in the design process is important. Therefore, the different sources of damping need to be modeled accurately. In this paper the impact of friction in blade attachments on forced response is investigated both numerically and experimentally. An efficient multiharmonic balance method is proposed in order to compute the forced response of bladed disks with contact and friction nonlinearities in blade roots. For experimental validation purposes, a rotating bladed disk was tested in a vacuum chamber, with excitation being provided by piezoelectric actuators. A model of the rig was built and numerical results were obtained with a normal load dependent coefficient of friction and a constant material damping ratio. Nonlinear behavior observed experimentally at resonances was well reproduced and an acceptable correlation was found with experimental resonant frequencies, amplitudes, and amount of damping throughout the spinning speed and excitation level range. The proposed numerical method can therefore serve to enhance the prediction of the alternating stresses in bladed disk assemblies.

Assessment of a semi-Hertzian method for determination of wheel–rail contact patch

Wheel–rail contact calculations are essential for simulating railway vehicle dynamic behavior. Currently, these simulations usually use the Hertz contact theory to calculate normal forces and Kalkers ‘FASTSIM’ program to evaluate tangential stresses. Since 1996, new methods called semi-Hertzian have appeared: 5 7 (STRIPES). These methods attempt to estimate the non-elliptical contact patches with a discrete extension of the Hertz theory. As a continuation of 2, a validation of the STRIPES method for normal problem computing on three test cases is proposed in this article. The test cases do not fulfill the hypothesis required for the Hertz theory. Then, the Kalkers FASTSIM algorithm is adapted to STRIPES patch calculus to perform tangential forces computation. This adaptation is assessed using Kalkers CONTACT algorithm.

Analysis of friction and instability by the centre manifold theory for a non-linear sprag-slip model

This paper presents the research devoted to the study of instability phenomena in non-linear model with a constant brake friction coefficient. Indeed, the impact of unstable oscillations can be catastrophic. It can cause vehicle control problems and component degradation. Accordingly, complex stability analysis is required. This paper outlines stability analysis and centre manifold approach for studying instability problems. To put it more precisely, one considers brake vibrations and more specifically heavy trucks judder where the dynamic characteristics of the whole front axle assembly is concerned, even if the source of judder is located in the brake system. The modelling introduces the sprag-slip mechanism based on dynamic coupling due to buttressing. The non-linearity is expressed as a polynomial with quadratic and cubic terms. This model does not require the use of brake negative coefficient, in order to predict the instability phenomena. Finally, the centre manifold approach is used to obtain equations for the limit cycle amplitudes. The centre manifold theory allows the reduction of the number of equations of the original system in order to obtain a simplified system, without loosing the dynamics of the original system as well as the contributions of non-linear terms. The goal is the study of the stability analysis and the validation of the centre manifold approach for a complex non-linear model by comparing results obtained by solving the full system and by using the centre manifold approach. The brake friction coefficient is used as an unfolding parameter of the fundamental Hopf bifurcation point.

The influence of crack-imbalance orientation and orbital evolution for an extended cracked Jeffcott rotor

Vibration peaks occurring at rational fractions of the fundamental rotating critical speed, here named Local Resonances, facilitate cracked shaft detection during machine shut-down. A modified Jeffcott-rotor on journal bearings accounting for gravity effects and oscillating around nontrivial equilibrium points is employed. Modal parameter selection allows this linear model to represent first mode characteristics of real machines. Orbit evolution and vibration patterns are analyzed, yielding useful results. Crack detection results indicate that, instead of 1x and 2x components, analysis of the remaining local resonances should have priority; this is due to crack-residual imbalance interaction and to 2x multiple induced origins. Therefore, local resonances and orbital evolution around 1/2, 1/3 and 1/4 of the critical speed are emphasized for various crack-imbalance orientations. To cite this article: J. Gomez-Mancilla et al., C. R. Mecanique 332 (2004).

Methods to reduce non-linear mechanical systems for instability computation

SummaryNon-linear dynamical structures depending on control parameters are encountered in many areas of science and engineering. In the study of non-linear dynamical systems depending on a given control parameter, the stability analysis and the associated non-linear behaviour in a near-critical steady-state equilibrium point are two of the most important points; they make it possible to validate and characterize the non-linear structures. Stability is investigated by determining eigenvalues of the linearized perturbation equations about each steady-state operating point, or by calculating the Jacobian of the system at the equilibrium points. While the conditions and the values of the parameters which cause instability can be investigated by using linearized equations of motion studies of the non-linear behaviour of vibration problems, on the other hand, require the complete non-linear expressions of systems. Due to the complexity of non-linear systems and to save time, simplifications and reductions in the mathematical complexity of the non-linear equations are usually required. The principal idea for these non-linear methods is to reduce the order of the system and eliminate as many non-linearities as possible in the system of equations.In this paper, a study devoted to evaluating the instability phenomena in non-linear models is presented. It outlines stability analysis and gives a non-linear strategy by constructing a reduced order model and simplifying the non-linearities, based on three non-linear methods: the centre manifold concept, the rational approximants and the Alternating Frequency/Time domain method. The computational procedures to determine the reduced and simplified system via the centre manifold approach and the fractional approximants, as well as the approximation of the responses as a Fourier series via the harmonic balance method, are presented and discussed. These non-linear methods for calculating the dynamical behaviour of non-linear systems with several degrees-of-freedom and non-linearities are tested in the case of mechanical systems with many degrees-of-freedom possessing polynomial non-linearities. Results obtained are compared with those estimated by a classical Runge-Kutta integration procedure.Moreover, an extension of the centre manifold approach using rational approximants is proposed and used to explore the dynamics of non-linear systems, by extending the domain of convergence of the non-linear reduced system and evaluating its performance and suitability.

Stability analysis of rotating beams rubbing on an elastic circular structure

This paper presents the stability analysis of a system composed of rotating beams on a flexible, circular fixed ring, using the Routh-Hurwitz criterion. The model displayed has been fully developed within the rotating frame by use of an energy approach. The beams considered possess two degrees of freedom (dofs), a flexural motion as well as a traction/compression motion. In-plane deformations of the ring will be considered. Divergences and mode couplings have thus been underscored within the rotating frame and in order to simplify understanding of all these phenomena, the dofs of the beams will first be treated separately and then together. The dynamics of radial rotating loads on an elastic ring can create divergence instabilities as well as post-critical mode couplings. Moreover, the flexural motion of beam rubbing on the ring can also lead to mode couplings and to the locus-veering phenomenon. The presence of rubbing seems to make the system unstable as soon as the rotational speed of the beams is greater than zero. Lastly, the influence of an angle between the beams and the normal to the rings inner surface will be studied with respect to system stability, thus highlighting a shift frequency phenomenon.

Forced Response Analysis of Integrally Bladed Disks With Friction Ring Dampers

This paper investigates a damping strategy for integrally bladed disks (blisks) based on the use of friction rings. The steady-state forced response of the blisk with friction rings is derived using the so-called dynamic Lagrangian frequency-time method adapted to cyclic structures with rotating excitations. In addition, an original approach for optimal determination of the number of Fourier harmonics is proposed. In numerical applications, a representative compressor blisk featuring several rings is considered. Each substructure is modeled using finite-elements and a reduced-order modeling technique is used for the blisk. The efficiency of this damping technology is investigated, and friction dissipation phenomena are interpreted with respect to frequency responses. It is shown that the friction damping effectiveness depends mainly on the level of dynamic coupling between blades and disk, and on whether the dynamics features significant alternating stick/slip phases. Through parameter studies, design guidelines are also proposed.

Dynamics of Multistage Bladed Disks Systems

This paper presents a new and original method for dynamical analysis of multistage cyclic structures such as turbomachinery compressors or turbines. Each stage is modeled cyclically by its elementary sector and the interstage coupling is achieved through a cyclic recombination of the interface degrees of freedom. This method is quite simple to set up; it allows us to handle the finite element models of each stages sector directly and, as in classical cyclic symmetry analysis, to study the nodal diameter problems separately. The method is first validated on a simple case study which shows good agreements with a complete 360 deg reference calculation. An industrial example involving two HP compressor stages is then presented. Then the forced response application is presented in which synchronous engine order type excitations are considered.