ON THE ELLIPTIC HARMONIC BALANCE METHOD FOR MIXED PARITY NON-LINEAR OSCILLATORS

2000 ◽  
Vol 233 (5) ◽  
pp. 935-937 ◽  
Author(s):  
M. BELHAQ ◽  
F. LAKRAD
Author(s):  
Christian Frey ◽  
Graham Ashcroft ◽  
Hans-Peter Kersken

This paper compares various approaches to simulate unsteady blade row interactions in turbomachinery. Unsteady simulations of turbomachinery flows have gained importance over the last years since increasing computing power allows the user to consider 3D unsteady flows for industrially relevant configurations. Furthermore, for turbomachinery flows, the last two decades have seen considerable efforts in developing adequate CFD methods which exploit the rotational symmetries of blade rows and are therefore up to several orders of magnitude more efficient than the standard unsteady approach for full wheel configurations. This paper focusses on the harmonic balance method which has been developed recently by the authors. The system of equations as well as the iterative solver are formulated in the frequency domain. The aim of this paper is to compare the harmonic balance method with the time-linearized as well as the non-linear unsteady approach. For the latter the unsteady flow fields in a fan stage are compared to reference results obtained with a highly resolved unsteady simulation. Moreover the amplitudes of the acoustic modes which are due to the rotor stator interaction are compared to measurement data available for this fan stage. The harmonic balance results for different sets of harmonics in the blade rows are used to explain the minor discrepancies between the time-linearized and unsteady results published by the authors in previous publications. The results show that the differences are primarily due to the neglection of the two-way coupling in the time-linearized simulations.


Author(s):  
Vincent Jaumouille´ ◽  
Jean-Jacques Sinou

Aeronautical structures are commonly assembled with bolted joints in which friction phenomena provide damping on the dynamic behaviour. Some models, mostly non linear, have consequently been developed and the harmonic balance method (HBM) is adapted to compute non linear response functions in the frequency domain. The basic idea is to develop the response as a Fourier series and to solve equations linking Fourier coefficients. One specific HBM feature is that response accuracy improves as the number of harmonics increases, at the expense of larger computational time. Thus the aim of this study is to develop an adaptive HBM which appreciates numerically the contribution of each harmonic on the dynamic response. For a given precision, the number of retained harmonics is adapted by an algorithm which integrates a numerical criterion based on an approximate strain energy. The application case is an asymmetrical two cantilever beam system linked by a bolted joint represented by a nonlinear LuGre model. Condensation and continuation methods are used to accelerate calculation. Adaptive HBM shows that, for a given value of the criterion, the number of harmonics may increase on resonances indicating that non linear effects are predominant.


2012 ◽  
Vol 19 (6) ◽  
pp. 1281-1295 ◽  
Author(s):  
V. Jaumouillé ◽  
J.-J. Sinou ◽  
B. Petitjean

In the presented work, a non linear effect of rubber referred as Fletcher-Gent effect or Payne effect is investigated. It leads to a change in the rubber dynamic modulus with vibration amplitudes and, consequently, modifies resonance frequencies of mechanical systems including non linear elastomers. In this study a new methodology is developed to take into account Payne effect in a linear viscoelastic rubber material. Small vibration amplitudes around a no-preloaded state are predicted by considering frequency and amplitude dependencies of the material. This methodology has the advantage of using tabular experimental data from characterization tests which avoids the development of a complex model. In order to compute frequency responses, the non linear harmonic balance method is used and, for each iteration, new rubber properties are affected at each element according to its strain state. An equivalent strain measure is evaluated from the element strain energy density. This equivalent strain allows to associate dynamic properties of a material element subjected to multiaxial strain state with experimental dynamic properties of a material sample subjected to an uniaxial strain state. Practically, DMAP alter procedures are developed in order to evaluate energies in models defined with MSC.Nastran and the non linear solver is developed with Matlab. The method is applied on a satellite instrument isolator including four non linear rubber mounts. A non homogeneous spatial distribution of element equivalent strains is observed. Moreover, the maximum equivalent strain varies with frequency. These two observations validate the use of a specific methodology to deal with amplitude dependency of rubber.


Sign in / Sign up

Export Citation Format

Share Document