Response Of Bladed Disks Research Articles

Intentional Blade Mistuning to Emulate the Dynamic Influence of Split Ring Damper

Abstract A novel methodology to predict responses of bladed disks incorporated with a split ring damper is proposed. To account for the dynamic characteristics of a split ring damper, intentional blade stiffness mistuning is presented. It is assumed that the mistuning is proportional to the blade stiffness matrix and uniformly distributed to all blades. Therefore, responses of bladed disks integrated with a split ring damper are captured by the detuned bladed disk integrated with a continuous ring damper. The proposed method is demonstrated through a lumped parameter model and validated using a complete finite element model. A key advantage of the proposed approach is that the dynamic characteristics of a split ring damper are emulated by a continuous ring damper in the detuned bladed disk. Moreover, the approach holds the potential to facilitate reduced-order modeling (ROM) by employing sector-level calculations. This results in substantial reductions in computational efforts while maintaining accurate predictions of forced responses.

Numerical Investigation of Friction Induced Interaction of Flutter Modes in a Realistic Low Pressure Turbine Rotor

Abstract Flutter response of bladed disks is a quite an important problem in modern designs of low pressure turbines, but, even for tuned configurations, it is not fully understood yet. Some research work suggests that, in a tuned rotor, the flutter induced vibration state that sets in consists of just a single traveling wave mode, while others show flutter states composed of several traveling waves active at the same time. Many of these studies were performed using conceptual models, such as mass-spring systems or asymptotically reduced models. In this work, a high fidelity bladed disk model is integrated into the time domain for the case of a realistic aerodynamically unstable low pressure turbine. In this configuration, the only nonlinear effect is the friction present at the blade disk attachment, which is responsible for the saturation of the flutter growth, and it is also the only mechanism that allows the interaction among the different unstable traveling wave modes. A case with multiple unstable modes of the same modal family is studied first. After that, the study is extended to a case where modes of different modal families are unstable and have similar aerodynamic instability levels. The interaction of unstable modes from different modal families is studied, we believe, for the first time, as previous flutter analyses were restricted to the interaction of unstable modes from the same family.

Study on the vibration characteristics of bladed disks (Effect of damping variation on the resonant responses of bladed disks)

Study on the vibration characteristics of bladed disks (Effect of damping variation on the resonant responses of bladed disks)

Validation of a Modal Work Approach for Forced Response Analysis of Bladed Disks

The paper describes a numerical method based on a modal work approach to evaluate the forced response of bladed disks and its validation against numerical results obtained by a commercial FEM code. Forcing functions caused by rotor–stator interactions are extracted from CFD unsteady solutions properly decomposed in time and space to separate the spinning perturbation acting on the bladed disk in a cyclic environment. The method was firstly applied on a dummy test case with cyclic symmetry where the forcing function distributions were arbitrarily selected: comparisons for resonance and out of resonance conditions revealed an excellent agreement between the two numerical methods. Finally, the validation was extended to a more realistic test case representative of a low-pressure turbine bladed rotor subjected to the wakes of two upstream rows: an IGV with low blade count and a stator row. The results show a good agreement and suggest computing the forced response problem on the finer CFD blade surface grid to achieve a better accuracy. The successful validation of the method, closely linked to the CFD environment, creates the opportunity to include the tool in an integrated multi-objective procedure able to account for aeromechanical aspects.

Comparative Studies of Surrogate Models for Response Analysis of Mistuned Bladed Disks

The forced responses of bladed disks are highly sensitive to inevitable random mistuning. Considerable computational efforts are required for the sampling process to assess the statistical vibration properties of mistuned bladed disks. Therefore, efficient surrogate models are preferred to accelerate the process for probabilistic analysis. In this paper, four surrogate models are utilized to construct the relation between random mistuning and forced response amplitudes, which are polynomial chaos expansion (PCE), response surface method (RSM), artificial neural networks (ANN) and Kriging interpolation, respectively. A bladed disk with 2-degrees-of-freedom (2-DOF) each sector is used to validate the effectiveness of the surrogate models. The effects of number of training samples on the surrogate model accuracy are discussed. The responses results of one blade (single output) and maximum response of all blades (multi-output) indicate that PCE and Kriging interpolation could yield accurate and stable predictions of the statistical characteristics of the forced responses. PCE is recommended for the mistuned response predictions due to its accuracy and efficiency.

Sensitivity and Forced Response Analysis of Anisotropy-Mistuned Bladed Disks With Nonlinear Contact Interfaces

AbstractA new method has been developed for the analysis of nonlinear forced response of bladed disks mistuned by blade anisotropy scatter and for the forced response sensitivity to blade material anisotropy orientations. The approach allows for the calculation of bladed disks with nonlinear friction contact interfaces using the multiharmonic balance method. The method uses efficient high-accuracy model reduction method for the minimization of the computational effort while providing required accuracy. The capabilities of the developed methods are validated and demonstrated using a two-blade model. A thorough study of the influence of the material anisotropy mistuning and its sensitivity on the characteristics of the forced response is carried out using finite element (FE) modes of anisotropy mistuned realistic bladed disk with nonlinear friction joints of blade roots and shroud contacts. The dependency of the nonlinear forced response on excitation level and contact pressure values has been carried out for anisotropy mistuned bladed disks.

Transient amplitude amplification of mistuned structures: An experimental validation

Until now, only few works have dealt with the transient vibration response of bladed disks of turbomachinery. However, transient vibration of turbines has become more and more important due to the changes on the energy production sector. Recently, it has been theoretically demonstrated that the maximum amplitude of the resonance passage of mistuned blisks may even exceed the maximum amplitude of the steady-state response. This effect is called transient amplitude amplification of mistuned structures (TAMS). The phenomenon is investigated within this article. By using a lumped mass model, it is shown how the maximum amplitude of the steady-state solution may be exceeded by the maximum amplitude of the resonance passage. Further, the occurrence of transient amplitude amplification of mistuned structures is experimentally validated.A new experimental set-up is designed to perform a resonance passage. Hence, an accelerated traveling wave type excitation is created. Loudspeakers are used to perform a contactless excitation. With the designed test rig, the occurrence of TAMS can be demonstrated. The results are the first experimental proof that the maximum of a transient vibration response can exceed the maximum of the stationary vibration response of the same dynamic system.

Vibration localization for a rotating mistuning bladed disk with the Coriolis effect by a state-space decoupling method

Slight deviations of blades due to manufacturing tolerances can cause mistuning of bladed disk leading to localized vibration, which can accelerate fatigue. Moreover, the rotating blades are subjected to the Coriolis effect and the influence of the Coriolis force on the natural frequencies of high-speed rotational bladed disks such as those of an aero-engine become more apparent. In this paper, the effect of Coriolis force on the forced response localization of a mistuned bladed disk are investigated, for conditions where the natural frequency located in the first and second modal families of the bladed disk. Mistuning is introduced by varying the Young’s modulus of each blade. Due to the asymmetric Coriolis matrix, it is not possible to directly decouple the system. A state-space decoupling method is developed to decouple the system to effectively calculate the forced response of bladed disk with the consideration of the Coriolis effect. The results show that response localization factor is increased by 13.09% considering the Coriolis force compared to the system without considering the Coriolis effect, in the case of where the first modal family is considered. In addition, the response localization factor with the consideration of the Coriolis force is decreased by 30.85% compared to the system without considering the Coriolis force, when the second modal family is considered. It indicates that the forced response localization with the consideration of the Coriolis effect will be changed obviously with the rotational speed increasing, when the concerning natural frequency is located in the first and second modal families. Furthermore, the effect of Coriolis force causes changes in the resonant frequencies and resonant amplitude, but does not introduce additional resonant peaks for the case of the mistuned bladed disk.

Reduced-Order Modeling of Bladed Disks With Friction Ring Dampers

An efficient methodology to predict the nonlinear response of bladed disks with a dry friction ring damper is proposed. Designing frictional interfaces for bladed-disk systems is an important approach to dissipate vibration energy. One emerging technology uses ring dampers, which are ringlike substructures constrained to move inside a groove at the root of the blades. Such rings are in contact with the bladed disk due to centrifugal forces, and they create nonlinear dissipation by relative motion between the ring and the disk. The analysis of the dynamic response of nonlinear structures is commonly done by numerical integration of the equations of motion, which is computationally inefficient, especially for steady-state responses. To address this issue, reduced-order models (ROMs) are developed to capture the nonlinear behavior due to contact friction. The approach is based on expressing the nonlinear forces as equivalent nonlinear damping and stiffness parameters. The method requires only sector-level calculations and allows precalculation of the response-dependent equivalent terms. These factors contribute to the increase of the computational speed of the iterative solution methods. A model of a bladed disk and damper is used to demonstrate the method. Macro- and micro-slip are used in the friction model to account for realistic behavior of dry friction damping. For validation, responses due to steady-state traveling wave excitations are examined. Results computed by ROMs are compared with results from transient dynamic analysis (TDA) in ansys with the full-order model. It is found that the steady-state responses predicted from the ROMs and the results from ansys are in good agreement, and that the ROMs reduce computation time significantly.

On the choice of contact parameters for the forced response calculation of a bladed disk with underplatform dampers

AbstractUnderplatform dampers (UPDs) are still in use to reduce the vibration amplitude of turbine blades and to shift the position of resonant frequencies. The dynamics of blades with UPDs is nonlinear and the analysis is challenging from both the experimental and the numerical point of view. A key point in obtaining a predictive numerical tool is the choice of the correct contact parameters (contact stiffness and friction coefficient) that are required as input to the contact model. The paper presents different approaches to choose these parameters: the contact stiffness in normal and tangential direction are both calculated and measured. The calculation is based on the analytical models in literature, the measurements are carried out on a dedicated test rig. The friction coefficient is also measured. Test results of the forced response of the same bladed disk with UPDs are available for each blade, they come from an experimental campaign under controlled excitation and centrifugal force. The forced response of the bladed disk is not used as a mean to tune the contact parameters, but rather as a validation tool: the effect of the different choices of contact parameters in the code is highlighted by the comparison of the calculated and experimental forced response of the bladed disk.

Vibration Characteristics of a Mistuned Bladed Disk considering the Effect of Coriolis Forces

To investigate the influence of Coriolis force on vibration characteristics of mistuned bladed disk, a bladed disk with 22 blades is employed and the effects of different rotational speeds and excitation engine orders on the maximum forced response are discussed considering the effects of Coriolis forces. The results show that if there are frequency veering regions, the largest split of double natural frequencies of each modal family considering the effects of Coriolis forces appears at frequency veering region. In addition, the amplitude magnification factor considering the Coriolis effects is increased by 1.02% compared to the system without considering the Coriolis effects as the rotating speed is 3000 rpm, while the amplitude magnification factor is increased by 2.76% as the rotating speed is 10000 rpm. The results indicate that the amplitude magnification factor may be moderately enhanced with the increasing of rotating speed. Moreover, the position of the maximum forced response of bladed disk may shift from one blade to another with the increasing of the rotational speed, when the effects of Coriolis forces are considered.

Optimization of Intentional Mistuning Patterns for the Mitigation of the Effects of Random Mistuning

This paper focuses on the optimization of intentional mistuning patterns for the reduction of the sensitivity of the forced response of bladed disks to random mistuning. Intentional mistuning is achieved here by using two different blade types (denoted as A and B) around the disk. It is thus desired to find the arrangement of these A and B blades (A/B pattern) that leads to the smallest 99th percentile of the amplitude of blade response in the presence of random mistuning. It is first demonstrated that there usually is a large number of local minima and further that the cost of a function evaluation is high. Accordingly, two novel, dedicated optimization algorithms are formulated and validated to address this specific problem. Both algorithms proceed in a two-step fashion. The first step, which consists of an optimization in a reduced space, leads to a series of good initial guesses. A local search from these initial guesses forms the second step of the methods. A detailed validation effort of this new procedure was next achieved on a single-degree-of-freedom-per-blade model, a reduced order model of a blisk, and that of an impeller considered in an earlier study. In all validation cases, the two novel algorithms were found to converge to the global optimum or very close to it at a small computational cost. Finally, the results of these optimization efforts further demonstrate the value of intentional mistuning to increase the robustness of bladed disks to random mistuning.

A Reduced-Order Meshless Energy Model for the Vibrations of Mistuned Bladed Disks—Part II: Finite Element Benchmark Comparisons

This paper draws upon the theoretical basis and applicability of the three-dimensional (3-D) reduced-order spectral-based “meshless” energy technology presented in a companion paper (McGee et al., 2013, “A Reduced-Order Meshless Energy Model for the Vibrations of Mistuned Bladed Disks—Part I: Theoretical Basis,” ASME J. Turbomach., to be published) to predict free and forced responses of bladed disks comprised of randomly mistuned blades integrally attached to a flexible disk. The 3-D reduced-order spectral-based model employed is an alternative choice in the computational modeling landscape of bladed disks, such as conventionally-used finite element methods and component mode synthesis techniques, and even emerging element-free Hamiltonian–Galerkin, Petrov–Galerkin, boundary integral, and kernel-particle methods. This is because continuum-based modeling of a full disk annulus of mistuned blades is, at present, a steep task using these latter approaches for modal-type mistuning and/or rogue blade failure analysis. Hence, a considerably simplified and idealized bladed disk of 20 randomly mistuned blades mounted to a flexible disk was created and modeled not only to analyze its free and forced 3-D responses, but also to compare the predictive capability of the present reduced-order spectral-based “meshless” technology to general-purpose finite element procedures widely-used in industry practice. To benchmark future development of reduced-order technologies of turbomachinery mechanics analysts may use the present 3-D findings of the idealized 20-bladed disk as a new standard test model. Application of the 3-D reduced-order spectral-based “meshless” technology to an industry integrally-bladed rotor, having all of its blades modally mistuned, is also offered, where reasonably sufficient upper-bounds on the exact free and forced 3-D responses are predicted. These predictions expound new solutions of 3-D vibration effects of modal mistuning strength and pattern, interblade mechanical coupling, and localized modes on the free and forced response amplitudes.

A Reduced-Order Meshless Energy Model for the Vibrations of Mistuned Bladed Disks—Part I: Theoretical Basis

In this paper, a reduced order model for the vibrations of bladed disk assemblies was achieved. The system studied was a 3D annulus of shroudless, “custom-tailored,” mistuned blades attached to a flexible disk. Specifically, the annulus was modeled as a spectral-based “meshless” continuum structure utilizing only nodal data to describe the arbitrary volume in which the system's dynamical energy was minimized. An extended Ritz variational procedure was used to minimize this energy, subjected to constraints imposed by an assumed 3D displacement field of mathematically complete, orthonormal “blade-disk” polynomials multiplied by generalized coefficients. The coefficients were determined by constraining the polynomial series to satisfy the extended Ritz stationary equations and essential boundary conditions of the bladed disk. From this, the governing equations of motion were generated into their usual dynamical forms to calculate upper-bounds on the actual free and forced responses of bladed disks. No conventional finite elements and element connectivity or component substructuring data were needed. This paper, Part I, outlines the theoretical foundation of the present model, and through extensive Monte Carlo simulations, establishes the analytical basis, predictive accuracy, and re-analysis efficiency of the present technology in the prediction of 3D maximum response amplitude of mistuned bladed disks having increasing numbers of nodal diameter excitations. Further applications validating the 3D approach against conventional finite element procedures of free and forced response prediction of a mistuned Integrally-Bladed Rotor used in practice is presented in a companion paper, Part II (Fang, McGee, and El-Aini, 2013, “A Reduced-Order Meshless Energy Model for the Vibrations of Mistuned Bladed Disks—Part II: Finite Element Benchmark Comparisons, ASME J. Turbomach., to be published.

Mistuning Identification of Integrally Bladed Disks With Cascaded Optimization and Neural Networks

Mistuning affects forced response of bladed disks drastically; therefore, its identification plays an essential role in the forced response analysis of bladed disk assemblies. Forced response analysis of mistuned bladed disk assemblies has drawn wide attention of researchers but there are a very limited number of studies dealing with identification of mistuning, especially if the component under consideration is an integrally bladed disk (blisk). This paper presents two new methods to identify mistuning of a bladed disk from the assembly modes via utilizing cascaded optimization and neural networks. It is assumed that a tuned mathematical model of the blisk under consideration is readily available, which is always the case for today’s realistic bladed disk assemblies. In the first method, a data set of selected mode shapes and natural frequencies is created by a number of simulations performed by mistuning the tuned mathematical model randomly. A neural network created by considering the number of modes, is then trained with this data set. Upon training the network, it is used to identify mistuning of the rotor from measured data. The second method further improves the first one by using it as a starting point of an optimization routine and carries out an optimization to identify mistuning. To carry out identification analysis by means of the proposed methods, there are no limitations on the number of modes or natural frequencies to be used. Thus, unlike existing mistuning identification methods they are suitable for incomplete data as well. Moreover, since system modes are used rather than blade alone counterparts, the techniques are ready to be used for analysis of blisks. Case studies are performed to demonstrate the capabilities of the new methods by using two different mathematical models to create training data sets a lumped-parameter model and a relatively realistic reduced order model. Throughout the case studies, the effects of using incomplete mode families and random errors in assembly modes are investigated. The results show that, the proposed method utilizing cascaded optimization and neural networks can identify mistuning parameters of a realistic blisk system with an exceptional accuracy even in the presence of incomplete and noisy test data.

Improved Modal Localization and Excitation Factors for Understanding Mistuned Bladed Disk Response

The nodal diameter spectrum is defined based on the discrete Fourier transform of the mode shapes of bladed disks. The dynamic characteristics analysis and assessment for a simplified bladed disk is performed using the concept of the nodal diameter spectrum. For free vibration, themode localization factor is defined based on the nodal diameter spectrum, and the probabilistic modal properties are obtained for the bladed disk with randommistuning of blade stiffness. For the forced vibration, the modal excitation factor is defined based on the nodal diameter spectrum to evaluate the similarity of engine-order excitation and themistunedmodes.The threshold phenomenon of random mistuning strength and the existence of optimal intentional mistuning are well explained through the combined effects of the modal excitation factor and the corresponding mode localization level of the mistuned mode shape on the response amplitude magnification. The numerical results indicate that the nodal diameter spectrum is an important parameter describing the spatial distributions of the mistuned mode shapes and could be widely applied in the analysis of the dynamic characteristics of mistuned bladed disks.

Dissipated energy and boundary condition effects associated to dry friction on the dynamics of vibrating structures

In turbomachines, dry friction devices (under platform dampers, shrouds, and tie-wire) are usually introduced to reduce resonant responses of bladed disks. Dry friction between rubbing elements induces a highly nonlinear dynamic behaviour which flattens the frequency response functions. It is clear that such behaviour requires an optimisation process to find the optimum parameters that lead to the minimum forced response amplitudes. However, different interpretations still remain concerning the explanation of the physical origin of this type of flattening. The most common one is based on dissipated energy. In this case, heat resulting from the relative frictional motion between contacting surfaces is supposed to bring sufficient dissipation to flatten response functions. On the other hand, a different approach considers that a decrease in vibrational amplitudes is explained by changes in boundary conditions induced by a stick/slip behaviour. In this study, a single degree-of-freedom system is used and analysed both in time and in frequency domains (Harmonic Balance Method) in order to show the contribution of respectively energy dissipation and change of contact state on peak levels.

The effect of underplatform dampers on the forced response of bladed disks by a coupled static/dynamic harmonic balance method

Friction contacts are often used in turbomachinery design as passive damping systems. In particular, underplatform dampers are mechanical devices used to decrease the vibration amplitudes of bladed disks. Numerical codes are used to optimize during designing the underplatform damper effectiveness in order to limit the resonant stress level of the blades. In such codes, the contact model plays the most relevant role in calculation of the dissipated energy at friction interfaces. One of the most important contact parameters to consider in order to calculate the forced response of blades assembly is the static normal load acting at the contact, since its value strongly affects the area of the hysteresis loop of the tangential force, and therefore the amount of dissipation. A common procedure to estimate the static normal loads acting on underplatform dampers consists in decoupling the static and the dynamic balance of the damper. A preliminary static analysis of the contact is performed in order to get the static contact/gap status to use in the calculation, assuming that it does not change when vibration occurs. In this paper, a novel approach is proposed. The static and the dynamic displacements of the system (bladed disk+underplatform dampers) are coupled together during the forced response calculation. Static loads acting at the contacts follow from static displacements and no preliminary static analysis of the system is necessary. The proposed method is applied to a numerical test case representing a simplified bladed disk with underplatform dampers. Results are compared with those obtained with the classical approach.

A test rig for the investigation of the dynamic response of a bladed disk with underplatform dampers

The paper presents the planning, structure and preliminary results of a new test rig for the validation of numerical models that calculate the dynamic response of bladed disks with underplatform dampers (UPDs).

A sparse preconditioned iterative method for vibration analysis of geometrically mistuned bladed disks

In this paper, the Static Mode Compensation method to predict geometrical mistuning effects on the response of bladed disks is reviewed and its limitations are analyzed. This method proved to be effective only for narrow clusters of modes under localized low rank perturbation. A new algorithm is introduced to address its deficiencies that draws on optimal preconditioned methods for generalized eigenvalue problem featuring sparse matrix vector multiplications, being more attractive under limited memory constraints of multi-millon DOF FEM models. The central idea of the SMC is to correct nominal eigenspace using modal acceleration method. It has been extended here by replacing the quasi-static set of modes with inexact solutions of the linear Jacobi–Davidson correction equations. Some heuristic strategies are discussed to lower the computational effort given the block-circulant structure of the nominal system. Numerical experiments on an industrial scale bladed disk model show that this leads to a very competitive tool. Computational performance and accuracy of both methods is compared in two areas of spectrum. The study demonstrates low accuracy of SMC method in the modal interaction zone, while validating efficiency and accuracy of the new algorithm in both areas.