scholarly journals Dual Time Stepping Algorithms With the High Order Harmonic Balance Method for Contact Interfaces With Fretting-Wear

Author(s):  
Loi¨c Salles ◽  
Laurent Blanc ◽  
Fabrice Thouverez ◽  
Alexander M. Gouskov ◽  
Pierrick Jean

Contact interfaces with dry friction are frequently used in turbomachinery. Dry friction damping produced by the sliding surfaces of these interfaces reduces the amplitude of bladed-disk vibration. The relative displacements at these interfaces lead to fretting-wear which reduces the average life expectancy of the structure. Frequency response functions are calculated numerically by using the multi-Harmonic Balance Method (mHBM). The Dynamic Lagrangian Frequency-Time method is used to calculate contact forces in the frequency domain. A new strategy for solving non-linear systems based on dual time stepping is applied. This method is faster than using Newton solvers. It was used successfully for solving Nonlinear CFD equations in the frequency domain. This new approach allows identifying the steady state of worn systems by integrating wear rate equations a on dual time scale. The dual time equations are integrated by an implicit scheme. Of the different orders tested, the first order scheme provided the best results.

Author(s):  
Loïc Salles ◽  
Laurent Blanc ◽  
Fabrice Thouverez ◽  
Alexander M. Gouskov ◽  
Pierrick Jean

Contact interfaces with dry friction are frequently used in turbomachinery. Dry friction damping produced by the sliding surfaces of these interfaces reduces the amplitude of bladed-disk vibration. The relative displacements at these interfaces lead to fretting-wear which reduces the average life expectancy of the structure. Frequency response functions are calculated numerically by using the multi-harmonic balance method (mHBM). The dynamic Lagrangian frequency-time method is used to calculate contact forces in the frequency domain. A new strategy for solving nonlinear systems based on dual time stepping is applied. This method is faster than using Newton solvers. It was used successfully for solving Nonlinear CFD equations in the frequency domain. This new approach allows identifying the steady state of worn systems by integrating wear rate equations a on dual time scale. The dual time equations are integrated by an implicit scheme. Of the different orders tested, the first order scheme provided the best results.


2011 ◽  
Vol 50-51 ◽  
pp. 583-588
Author(s):  
Yong Qiang Shi ◽  
Qing Zhen Yang ◽  
Xin Hai Zhou

A harmonic balance approach has been developed to compute nonlinear viscous unsteady flows around oscillating blades. The computed results using two orders harmonic balance method are compared with those by conventional dual-time stepping method. Results obtained with the present method agree well with those from dual-time stepping method, which demonstrate the ability of the present analysis method to model accurately the unsteady flow. Furthermore,the present method is highly efficient. It is about 36 times fast than conventional dual-time stepping method in the present computation. Then the effects of oscillation amplitude and reduced frequency on unsteadiness of flows are studied. The analysis exploits the fact that, (1) the hysteresis effect of unsteady flow is hardly affected by oscillation amplitude, but the first harmonic unsteady pressure across the blade is proportional to oscillation amplitude; (2) the higher the reduced frequency, the wider the range of unsteady aerodynamic forces, the more intense the hysteresis effect.


2014 ◽  
Vol 2014 ◽  
pp. 1-8
Author(s):  
Hongliang Yao ◽  
Qian Zhao ◽  
Qi Xu ◽  
Bangchun Wen

The efficiency and accuracy of common time and frequency domain methods that are used to simulate the response of a rotor system with malfunctions are compared and analyzed. The Newmark method and the incremental harmonic balance method are selected as typical representatives of time and frequency domain methods, respectively. To improve the simulation efficiency, the fixed interface component mode synthesis approach is combined with the Newmark method and the receptance approach is combined with the incremental harmonic balance method. Numerical simulations are performed for rotor systems with single and double frequency excitations. The inherent characteristic that determines the efficiency of the two methods is analyzed. The results of the analysis indicated that frequency domain methods are suitable single and double frequency excitation rotor systems, whereas time domain methods are more suitable for multifrequency excitation rotor systems.


Author(s):  
Shangguan Bo ◽  
Zili Xu ◽  
Qilin Wu ◽  
XianDing Zhou ◽  
ShouHong Cao

To understand the mechanism of interfacial damping of axial loosely assembled dovetail to suppress blade vibration, a dry friction force model is presented by the Coulomb friction law and the macroslip model, and the mathematical expression of the friction force is derived. The nonlinear friction force is linearized as an equivalent stiffness and an equivalent damping through the one-term harmonic balance method. The effect of centrifugal force on the equivalent stiffness and the equivalent damping is studied. The forced response of one simplified blade with loosely assembled dovetail attachment is predicted by the harmonic balance method, in which the blade is described by the lumped mass and spring model, and the friction contact joints is simplified as a ideal friction damper. The results show that the equivalent stiffness of loosely assembled dovetail attachment increases with blade centrifugal force, gradually reaches a certain value, and there exists the maximum value for the equivalent stiffness. The equivalent damping increases at the beginning and then decreases with blade centrifugal force increasing, there exists a maximum too. The resonant frequency of blade rises with blade centrifugal force, but it no longer increases when the centrifugal force exceed a certain value. There exists a special centrifugal force on which the effect of dry friction damping is the best.


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.


1985 ◽  
Vol 52 (4) ◽  
pp. 958-964 ◽  
Author(s):  
C. Pierre ◽  
A. A. Ferri ◽  
E. H. Dowell

A multi-harmonic, frequency domain analysis of dry friction damped systems is presented which uses an incremental harmonic balance method. When compared with time domain solution methods, it is found that the incremental harmonic balance method can yield very accurate results with some advantages over the time domain methods. Both one and two degree-of-freedom systems are studied.


Author(s):  
Yann Colaïtis ◽  
Alain Batailly

Abstract In this study, a frequency-domain approach based on the harmonic balance method coupled to a predictor-corrector continuation algorithm is implemented for the qualitative analysis of blade-tip/casing contacts in aircraft engines. Unilateral contact and dry friction are taken into account through a regularized penalty law. To enhance the robustness of the methodology, particular attention is paid to the mitigation of the Gibbs phenomenon. To this end, the employed Alternating Frequency/Time scheme features a Lanczos σ-approximation so that spurious oscillations of the computed nonlinear contact forces become negligible. This approach is applied in combination with a model reduction technique on an industrial compressor blade: NASA rotor 37. In order to assess the influence of both the contact law regularization and the Lanczos σ-approximation, obtained results are thoroughly compared to an existing time integration-based numerical strategy relying on a Lagrange multiplier-based approach for contact treatment and that was previously confronted to experimental results. Presented results underline the very good agreement between the proposed methodology and the reference time integration numerical strategy. The proposed developments thus complement existing results on blade-tip/casing contact adding a much needed qualitative understanding of the interaction and an accurate assessment of the contact stiffening phenomenon.


2021 ◽  
Author(s):  
Yann Colaïtis ◽  
Alain Batailly

Abstract In this study, a frequency-domain approach based on the harmonic balance method coupled to a predictor-corrector continuation algorithm is implemented for the qualitative analysis of blade-tip/casing contacts in aircraft engines. Unilateral contact and dry friction are taken into account through a regularized penalty law. To enhance the robustness of the methodology, particular attention is paid to the mitigation of the Gibbs phenomenon. To this end, the employed Alternating Frequency/Time scheme features a Lanczos σ-approximation so that spurious oscillations of the computed nonlinear contact forces become negligible. This approach is applied in combination with a model reduction technique on an industrial compressor blade: NASA rotor 37. In order to assess the influence of both the contact law regularization and the Lanczos σ-approximation, obtained results are thoroughly compared to an existing time integration-based numerical strategy relying on a Lagrange multiplier-based approach for contact treatment and that was previously confronted to experimental results. Presented results underline the very good agreement between the proposed methodology and the reference time integration numerical strategy. The proposed developments thus complement existing results on blade-tip/casing contact adding a much needed qualitative understanding of the interaction and an accurate assessment of the contact stiffening phenomenon.


Sign in / Sign up

Export Citation Format

Share Document