Probabilistic Analysis of Geometric Uncertainty Effects on Blade Modal Response

Author(s):  
Jeffrey M. Brown ◽  
Joeseph Slater ◽  
Ramana V. Grandhi

This paper investigates the effect of manufacturing variations on the modal response of a transonic low aspect ratio fan. A simulated set of coordinate measurement machine measurements from a single rotor, representative of actual manufacturing variations, are used to investigate geometric effects. Principal component analysis is used to statistically model spatial geometry variations and reduce variable space dimensionality. Statistics from this analysis are used with Monte Carlo sampling to generate random blades realizations that are used to predict response distributions for a simulated fleet of 1000 blades. An existing approach to approximate blade frequency response is extended to include modal displacement and stress. These approximations are based on eigensensitivity analysis and first order Taylor series approximations. An approximation error analysis is conducted to quantify accuracy. The effect of small geometry variations on blade natural frequency, mode shape, and modal stress is investigated with results showing that small variations on the order of mils can cause significant variations in both scale and location of free and forced response.

Author(s):  
Jeffrey M. Brown ◽  
Ramana V. Grandhi

This paper investigates the effect of manufacturing variations on the blade-alone forced response of a transonic low aspect ratio fan. A simulated set of coordinate measurement machine measurements from a single rotor, representative of actual manufacturing variations, are used to investigate geometric effects. A reduced order model is developed to rapidly solve for the forced response and is based on eigensensitivity analysis and dynamic response mode superposition. An approximation error analysis is conducted to quantify accuracy of the new tool and errors between approximate and full finite element analysis solutions are shown to be small for low order modes with some high order modes having moderate error. A study of the simulated measured blade results show a significant amount of forced response variation along the leading edge of the airfoil. Statistics from this simulated measured rotor are used with Monte Carlo sampling to generate random blades realizations that are solved with the reduced order model. This procedure allows the prediction of the variation across an entire fleet of blades from a small sample of blades. The large variations predicted, up to 40%, could have a significant impact of the blade design process including the procedures to account for foreign object damage damage tolerance, how non-intrusive stress measurement systems are used, and how mistuning prediction algorithms are validated.


Author(s):  
Marco Gambitta ◽  
Arnold Kühhorn ◽  
Sven Schrape

Abstract The present work focuses on the effect of the manufacturing geometrical variability on the high-pressure compressor of a turbofan engine for civil aviation. The deviations of the geometry over the axial compressor blades are studied and modeled for the representation in the computational models. Such variability is of particular interest for the forced response problem, where small deviations of the geometry from the ideal nominal model can cause significant differences in the vibrational responses. The information regarding the geometrical mistuning is extracted from a set of manufactured components surface scans of a blade integrated disk (blisk) rotor. The optically measured geometries are parameterized, defining a set of opportune variables to describe the deviations. The dimension of the variables domain is reduced using the principal component analysis approach and a reconstruction of the modeled geometries is performed for the implementation in CFD and FEM solvers. The generated model allows a stochastic representation of the variability, providing an optimal set of variables to represent it. The aeroelastic analyses considering geometry based mistuning is carried out on a test-rig case, focusing on how such variability can affect the modal forcing generated on the blades. The force generated by the unsteady pressure field over the selected vibrational mode shapes of the rotor blades is computed through a validated CFD model. The uncertainty quantification of the geometrical variability effect on the modal forcing is performed employing Monte Carlo methods on a reduced model for the CFD solution, based on a single passage multi-blade row setup. The amplitude shift of the unsteady modal forcing is studied for different engine orders. In particular the scatter of the main engine orders forcing amplitudes for the manufactured blades can be compared with the nominal responses to predict the possible amplification due to the geometrical variability. Finally the results are compared to a full assembly computational model to assess the influence of multiple variable blades.


Author(s):  
Nicolas Driot ◽  
Alain Berlioz ◽  
Claude-Henri Lamarque

The aim of this work is to apply stochastic methods to investigate uncertain parameters of rotating machines with constant speed of rotation subjected to a support motion. As the geometry of the skew disk is not well defined, randomness is introduced and affects the amplitude of the internal excitation in the time-variant equations of motion. This causes uncertainty in dynamical behavior, leading us to investigate its robustness. Stability under uncertainty is first studied by introducing a transformation of coordinates (feasible in this case) to make the problem simpler. Then, at a point far from the unstable area, the random forced steady state response is computed from the original equations of motion. An analytical method provides the probability of instability, whereas Taguchi’s method is used to provide statistical moments of the forced response.


2015 ◽  
Vol 137 (7) ◽  
Author(s):  
Joseph A. Beck ◽  
Jeffrey M. Brown ◽  
Alex A. Kaszynski ◽  
Charles J. Cross ◽  
Joseph C. Slater

New geometric mistuning modeling approaches for integrally bladed rotors (IBRs) are developed for incorporating geometric perturbations to a fundamental disk–blade sector, particularly the disk–blade boundary or connection. Reduced-order models (ROMs) are developed from a Craig–Bampton component mode synthesis (C–B CMS) framework that is further reduced by a truncated set of interface modes that are obtained from an Eigen-analysis of the C–B CMS constraint degrees of freedom (DOFs). An investigation into using a set of tuned interface modes and tuned constraint modes for model reduction is then performed, which offers significant computational savings for subsequent analyses. Two configurations of disk–blade connection mistuning are investigated: as-measured principal component (PC) deviations and random perturbations to the interblade spacing. Furthermore, the perturbation sizes are amplified to investigate the significance of incorporating mistuned disk–blade connections during solid model generation from optically scanned geometries. Free and forced response results are obtained for each ROM and each disk–blade connection type and compared to full finite element model (FEM) solutions. It is shown that the developed methods provide accurate results with a reduction in solution time compared to the full FEM. In addition, results indicate that the inclusion of a mistuned disk–blade connection deviations are small or conditions where large perturbations are localized to a small areas of the disk–blade connection.


Author(s):  
Vinod Vishwakarma ◽  
Alok Sinha ◽  
Yasharth Bhartiya ◽  
Jeffery M. Brown

Modified modal domain analysis (MMDA), a reduced order modeling technique, is applied to a geometrically mistuned integrally bladed rotor to obtain its natural frequencies, mode shapes, and forced response. The geometric mistuning of blades is described in terms of proper orthogonal decomposition (POD) of the coordinate measurement machine (CMM) data. Results from MMDA are compared to those from the full (360 deg) rotor Ansys model. It is found that the MMDA can accurately predict natural frequencies, mode shapes, and forced response. The effects of the number of POD features and the number of tuned modes used as bases for model reduction are examined. Results from frequency mistuning approaches, fundamental mistuning model (FMM) and subset of nominal modes (SNM), are also generated and compared to those from full (360 deg) rotor Ansys model. It is clearly seen that FMM and SNM are unable to yield accurate results whereas MMDA yields highly accurate results.


Author(s):  
Vsevolod Kharyton ◽  
Jean-Pierre Laine ◽  
Fabrice Thouverez ◽  
Olexiy Kucher

The study intends to simulate the process of the blade tip amplitude calculation by the tip-timing method. An attention is focused on tip-timing measurements for detection of a cracked blade from the bladed disk forced response. The cracked blade is considered within frameworks of the bladed disk dynamic model that takes into account mistuning presence. Nonlinear formulation of a crack behavior is done with the harmonic balance method in its combination with the contact analysis that allows simulation of crack breathing. In order to make the cracked blade detection process evident, the crack length and location are set in such a way as to produce the cracked blade frequency localization. Reconstruction of the blade tip amplitudes is attained with the arriving time of measured probes of the blade tips. The results are compared with the blade forced response obtained by the bladed disk dynamic model. A possibility is also considered how to reconstruct time-history of the bladed disk forced response with tip-timing data.


Author(s):  
Peter Ho¨nisch ◽  
Arnold Ku¨hhorn ◽  
Bernd Beirow

The effect of blade frequency mistuning on the forced response of integral radial turbines is studied by means of experimental and numerical analyses. Blade dominated frequencies representing the mistuning are identified based on blade by blade measurements using the example of a MTU ZR140 turbine blisk. Based on these results, numerical simulations of the blade by blade measurements are performed, aiming to update the originally ideal (tuned) finite element model. The damping information to be considered in the update process is taken from results of an experimental modal analysis. The quality of the model is proved by well correlated frequency response functions (FRF) of numerical and experimental analyses. Finally, the models are used to simulate the forced response due to travelling wave excitations. As a result, mode localization phenomena and response amplifications compared to tuned blisks are proved. In order to round off the contribution to a more enhanced understanding of the radial turbine blisk dynamics optically based geometry measurements are performed to assess the influence of geometrical deviations on frequency mistuning. It is shown that geometric imperfections can be the main driver causing a mistuned response characteristic.


Author(s):  
Bernd Beirow ◽  
Arnold Kühhorn ◽  
Robby Weber ◽  
Frederik Popig

Abstract The last stage bladed disk of a steam turbine is analyzed with respect to both flutter susceptibility and limitation of forced response. Due to the lack of variable stator vanes unfavorable flow conditions may occur which increases the risk of flutter at part load conditions. For this reason, intentional mistuning is employed with the objective to prevent any self-excited vibrations. A first step in this direction is done by choosing alternate mistuning, which keeps the manufactural efforts in limits. In this sense, two different series of blades have been made. However, small deviations from the design intention are unavoidable due to the manufacturing procedure, which could be proved by bonk tests carried out earlier. The influence of these additional deviations is considered in numerical simulations. Moreover, the strong dependence of blade frequencies on the speed is taken into account since centrifugal stiffening effects significantly attenuate the blade-to-blade frequency difference. Focusing on the first flap mode it could be shown that a mitigation of flutter susceptibility is achieved by prescribing alternate mistuning, which indeed evokes an increase of originally small aerodynamic damping ratios. Nevertheless, the occurrence of negative damping ratios could not be completely precluded at part load conditions. That is why optimization studies are conducted based on genetic algorithms with the objective function of maximizing the lowest aerodynamic damping ratios. Finally, mistuning patterns could be identified featuring a tremendous increase of aerodynamic damping ratios. The robustness of the solutions could be proved by superimposing additional random mistuning.


2009 ◽  
Vol 131 (4) ◽  
Author(s):  
Tomokazu Miyakozawa ◽  
Robert E. Kielb ◽  
Kenneth C. Hall

Most of the existing mistuning research assumes that the aerodynamic forces on each of the blades are identical except for an interblade phase angle shift. In reality, blades also undergo asymmetric steady and unsteady aerodynamic forces due to manufacturing variations, blending, mis-staggered, or in-service wear or damage, which cause aerodynamically asymmetric systems. This paper presents the results of sensitivity studies on forced response due to aerodynamic asymmetry perturbations. The focus is only on the asymmetries due to blade motions. Hence, no asymmetric forcing functions are considered. Aerodynamic coupling due to blade motions in the equation of motion is represented using the single family of modes approach. The unsteady aerodynamic forces are computed using computational fluid dynamics (CFD) methods assuming aerodynamic symmetry. Then, the aerodynamic asymmetry is applied by perturbing the influence coefficient matrix in the physical coordinates such that the matrix is no longer circulant. Therefore, the resulting aerodynamic modal forces in the traveling wave coordinates become a full matrix. These aerodynamic perturbations influence both stiffness and damping while traditional frequency mistuning analysis only perturbs the stiffness. It was found that maximum blade amplitudes are significantly influenced by the perturbation of the imaginary part (damping) of unsteady aerodynamic modal forces. This is contrary to blade frequency mistuning where the stiffness perturbation dominates.


Sign in / Sign up

Export Citation Format

Share Document