Three-Dimensional Tailored Vibration Response and Flutter Control of High-Bypass Shroudless Aeroengine Fans

2013 ◽  
Vol 135 (2) ◽  
Author(s):  
O. G. McGee III ◽  
C. Fang

A new reduced-order design synthesis technology has been developed for vibration response and flutter control of cold-stream, high-bypass ratio, shroudless, aeroengine fans. To simplify the design synthesis (optimization) of the fan, a significant order reduction of the mechanical response and stiffness-shape design synthesis has been achieved. The assumed cyclic symmetric baseline fan is modeled as a cascade of tuned, shroudless, arbitrarily shaped, wide-chord laminated composite blades, each with a reduced order of degrees of freedom using a three-dimensional (3D) elasticity spectral-based energy model (McGee et al., 2013, “A Reduced-Order Meshless Energy Model for the Vibrations of Mistuned Bladed Disks—Part I: Theoretical Basis, ASME J. Turbomach., in press; Fang et al., 2013, “A Reduced-Order Meshless Energy Model for the Vibrations of Mistuned Bladed Disks—Part II: Finite Element Benchmark Comparisons, ASME J. Turbomach., in press). The uniqueness of the mechanical analysis is that the composite fan was modeled as a “meshless” continuum, consisting of nodal point data to describe the arbitrary volume. A stationary value of energy within the arbitrarily shaped composite fan annulus was achieved using an extended spectral-based Ritz procedure to obtain the dynamical equations of motion for 3D free vibration response of a rotating composite high-bypass fan. No additional kinematical constraints (as in beam, plate, or shell theories) were utilized in the 3D elasticity-based energy formulation. The convergence accuracy of the spectral-based 3D free vibration response predictions was nearly one percent upper-bounds on the exact mechanical response of the baseline composite fan, particularly in the lowest five modes studied closely in this work, as typically seen with spectral-based Ritz procedures employed in the analysis. The spectral-based 3D predictions was validated against those predicted using a general purpose finite element technology widely used by industry. In off-design operation, the frequency margins of the lower flex-torsion modes of a fan may be dangerously close to integral-order resonant and empirical stall flutter boundaries. For a given baseline composite fan, it is proposed that to reduce the likelihood of resonant response and flutter on a Campbell diagram, design analysts can efficiently unite the newly developed reduced-order 3D spectral-based energy reanalysis within a novel reduced-order spectral-based Kuhn–Tucker optimality design synthesis procedure to fairly accurately restructure the Campbell diagram of a composite high-bypass ratio fan using stiffness optimization (by means of proper choices of angle-ply orientations of the blade laminates) and mass-balancing (shape) optimization (by way of blade thickness variation tuning of the lower aerodynamic loading portion of the blades between the dovetail root section and the midradial height section of the composite fan annulus). Fan design optima is summarized that (1) achieves multiple frequency margins and satisfies multiple empirical stall flutter constraints, (2) controls the twist-flex vibratory response in the lowest (fundamental) mode, and (3) ensures the mechanical strength integrity of the optimized angle-ply lay-up under steady centrifugal tension and gas bending stresses. Baseline and optimally restructured Campbell diagrams and design sensitivity calculations are presented, comparing optimum solution accuracy and validity of the proposed reduced-order spectral-based design synthesis technology against optimum solutions generated from open-source nonlinear mathematical programming software (i.e., NASA’s general-purpose sequential unconstrained minimization technique, Newsumt-A) (Miura and Schmit, Jr., 1979, ”NEWSUMT–A, Fortran Program for Inequality Constrained Function Minimization—Users Guide,“ NASA CR-159070). Design histories of fan stiffness and mass balancing (or shape) along with nondimensional constraints (i.e., frequency margins, reduced frequencies, twist-flex vibratory response, first-ply failure principal stress limits, and dovetail-to-midblade height thickness distribution) show that a proper implementation of fan stiffness tailoring (via symmetric angle-ply orientations) and mass-balancing (thickness) optimization of the fan assembly produces a feasible Campbell diagram that satisfies all design goals. An off-design analysis of the optimized fan shows little sensitivity to twist-flex coupling response and flutter with respect to small variability or errors in optimum design construction. Industry manufacturing processes may introduce these small errors known as angle-ply laminate construction misalignments (Graham and Guentert, 1965, “Compressor Stall and Blade Vibration,” Aerodynamic Design of Axial-Flow Compressors, Chap. XI, NASA SP-36; Meher-Hornji, 1995, “Blading Vibration and Failures in Gas Turbines, Part A: Blading Dynamics and the Operating Environment,” ASME Paper 95-GT-418; Petrov et al., 2002, “A New Method for Dynamic Analysis of Mistuned Bladed Disks Based on the Exact Relationship Between Tuned and Mistuned Systems,” ASME J. Eng. Gas Turbines Power, 124(3), pp. 586–597; Wei and Pierre, 1990, “Statistical Analysis of the Forced Response of Mistuned Cyclic Assemblies,” ASME J. Eng. Gas Turbines Power, 28(5), pp. 861–868; Wisler, 1988, “Advanced Compressor and Fan Systems,” GE Aircraft Engines, Cincinnati, Ohio (also 1986 Lecture to ASME Turbomachinery Institute, Ames Iowa)).

Author(s):  
O. G. McGee ◽  
C. Fang

A new reduced-order design synthesis technology has been developed for vibration response and flutter control of cold-stream, high-bypass ratio, shroudless fans. To simplify the design synthesis (optimization) of the fan, a significant order reduction of the mechanical response and stiffness-shape design synthesis has been achieved. The assumed cyclic symmetric baseline fan is modeled as a cascade of tuned, shroudless, arbitrarily-shaped, wide-chord blades, each with a reduced-order of degrees of freedom using a three-dimensional (3-D) elasticity-based, meshless energy model [Fang, (2002); McGee et al (2008); Fang et al (2008)]. The convergence accuracy and mechanical response error of the present 3-D predictions were estimated nearly one percent above the exact mechanical response of the baseline fan. In off-design operation, the frequency margins of the lower flex-torsion modes of a fan may be dangerously close to integral-order resonant and empirical stall flutter boundaries. It is shown that an optimized mechanical stiffness through material properties (via. symmetric angle-ply orientations) and an optimized fan shape (via. variation of blade thickness from hub-to-mid-radial height) can be found to reduce the likelihood of resonant response and flutter on a Campbell diagram. A baseline fan is numerically optimized using a first-of-its-kind reduced-order design synthesis technology involving a solution of simultaneous nonlinear partial differential equations determining the necessary and sufficient Kuhn-Tucker conditions of optimality of constrained minimization. Solution accuracy and validity of the reduced-order design synthesis technology is benchmarked against a widely-used conventional method of nonlinear programming (via. sequential unconstrained minimization technique). Fan design optima is obtained that (1) achieves multiple frequency margins and satisfies multiple empirical stall flutter constraints, (2) controls the twist-flex vibratory response in the lowest (fundamental) mode, and (3) ensures the mechanical strength integrity of the optimized angle-ply lay-up under steady centrifugal tension and gas bending stresses. Baseline and optimally-restructured Campbell diagrams and design sensitivity calculations are presented. Design histories of fan stiffness and shape and nondimensional constraints (i.e., frequency margins, reduced frequencies, twist-flex vibratory response, first-ply failure principal stress limits, and hub-to-mid-blade height thickness distribution) show that a proper implementation of fan stiffness tailoring and shape (thickness) optimization of the fan assembly produces a feasible Campbell diagram that satisfies all design goals. An off-design analysis of the optimized fan shows little sensitivity to twist-flex coupling response and flutter with respect to small variability or errors in optimum design construction. Industry manufacturing processes may introduce these small errors. Finally, the developed reduced-ordered technology of fan design is shown to be highly cost-effective and accurate, when its predictive mechanical response capability is compared to general-purpose finite element technology widely-used by industry.


2018 ◽  
Vol 25 (2) ◽  
pp. 298-309
Author(s):  
Hao Zhou ◽  
Sheng Meng ◽  
Chunhong Mo ◽  
Lujun Wang ◽  
Xiukui Hu ◽  
...  

Thermoacoustic oscillation occurs in modern industrial furnaces, gas turbines, and liquid rockets. However, the thermoacoustic prediction tools for furnaces vibration are less developed. This paper presents a one-dimensional (1D) linear acoustic approach to analyze the three-dimensional acoustic modes of a 660 MWe oil-fuel furnace. The interaction between the flame and acoustic field is described with the flame transfer function. The global time delay is evaluated through a Reynolds averaged simulation. The results of the 1D acoustic approach are compared with real furnace test data. The unstable modes are close to the natural modes of the furnace, and the 30 Hz in the longitudinal mode is the strongest vibration frequency. The effects of inlet length reduction and separation plate removal are also examined. When the separation plates are removed, the time lag of flame in response to inlet flow decreases from 52.5 milliseconds (ms) to 43.8 ms. The results of the 1D approach and finite element method (FEM) show a same safe operation window. The reduced-order procedure and FEM adopted in this study give us a solution to mitigate the vibration in the furnace.


2013 ◽  
Vol 135 (6) ◽  
Author(s):  
C. Fang ◽  
O. G. McGee ◽  
Y. El Aini

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.


Author(s):  
Jean de Cazenove ◽  
Moustapha Mbaye ◽  
Jean-Philippe Ousty ◽  
Scott Cogan

This paper proposes an original approach to the reduced-order modelling of integrally bladed disks. It is proposed to build a reduction basis which is independent from the rotational speed, from only one modal, cyclic-symmetry calculation performed at rest, and a few static computations. Based on previous works, a polynomial expansion is used which leads to a parametric approximation of the stiffness matrix for the entire operating range. Furthermore, the Kirsch Combined Approximation (KCA) method is used for building the final reduction basis. This method is based on successive approximations of the negative binomial expansion applied to the reanalysis eigenproblem. After giving a general overview of the main theoretical aspects, the paper focuses on the reanalysis problem based on the combined approximations method. Finally, the application of the extended reduction method to the case of a real bladed disk is presented. It is shown that the use of combined approximations provides a very accurate estimation of a Campbell diagram, and allows substantial computational time savings.


Author(s):  
C C Siew ◽  
M Hill ◽  
R Holmes ◽  
M J Brennan

This paper uses two reduced-order methods to calculate the overall non-linear vibration response of a multi-mode rotor-squeeze film damper (SFD) bearing assembly. Good agreement has been found between their computed results and those calculated by the conventional Runge-Kutta-Merson method (RKMM), yet they require less than 5 per cent of the computing time consumed by the RKMM. They are applied to compute the vibration response, three-dimensional deflection shape and the overall vibrational kinetic energy of the rotor-bearing assembly. The assembly is simulated under various operating conditions including different unbalances, oil viscosities, static misalignments and damper factors. It is concluded that effective vibration control of a flexible rotor-SFD bearing assembly can be achieved when the vertical static eccentricity ratio of the SFD is set to a certain limit.


2013 ◽  
Vol 135 (6) ◽  
Author(s):  
O. G. McGee ◽  
C. Fang ◽  
Y. El-Aini

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.


Author(s):  
D. L. Callahan

Modern polishing, precision machining and microindentation techniques allow the processing and mechanical characterization of ceramics at nanometric scales and within entirely plastic deformation regimes. The mechanical response of most ceramics to such highly constrained contact is not predictable from macroscopic properties and the microstructural deformation patterns have proven difficult to characterize by the application of any individual technique. In this study, TEM techniques of contrast analysis and CBED are combined with stereographic analysis to construct a three-dimensional microstructure deformation map of the surface of a perfectly plastic microindentation on macroscopically brittle aluminum nitride.The bright field image in Figure 1 shows a lg Vickers microindentation contained within a single AlN grain far from any boundaries. High densities of dislocations are evident, particularly near facet edges but are not individually resolvable. The prominent bend contours also indicate the severity of plastic deformation. Figure 2 is a selected area diffraction pattern covering the entire indentation area.


2002 ◽  
Vol 124 (4) ◽  
pp. 988-993 ◽  
Author(s):  
V. Esfahanian ◽  
M. Behbahani-nejad

An approach to developing a general technique for constructing reduced-order models of unsteady flows about three-dimensional complex geometries is presented. The boundary element method along with the potential flow is used to analyze unsteady flows over two-dimensional airfoils, three-dimensional wings, and wing-body configurations. Eigenanalysis of unsteady flows over a NACA 0012 airfoil, a three-dimensional wing with the NACA 0012 section and a wing-body configuration is performed in time domain based on the unsteady boundary element formulation. Reduced-order models are constructed with and without the static correction. The numerical results demonstrate the accuracy and efficiency of the present method in reduced-order modeling of unsteady flows over complex configurations.


Author(s):  
Stephan Uhkoetter ◽  
Stefan aus der Wiesche ◽  
Michael Kursch ◽  
Christian Beck

The traditional method for hydrodynamic journal bearing analysis usually applies the lubrication theory based on the Reynolds equation and suitable empirical modifications to cover turbulence, heat transfer, and cavitation. In cases of complex bearing geometries for steam and heavy-duty gas turbines this approach has its obvious restrictions in regard to detail flow recirculation, mixing, mass balance, and filling level phenomena. These limitations could be circumvented by applying a computational fluid dynamics (CFD) approach resting closer to the fundamental physical laws. The present contribution reports about the state of the art of such a fully three-dimensional multiphase-flow CFD approach including cavitation and air entrainment for high-speed turbo-machinery journal bearings. It has been developed and validated using experimental data. Due to the high ambient shear rates in bearings, the multiphase-flow model for journal bearings requires substantial modifications in comparison to common two-phase flow simulations. Based on experimental data, it is found, that particular cavitation phenomena are essential for the understanding of steam and heavy-duty type gas turbine journal bearings.


2003 ◽  
Vol 2 (1) ◽  
pp. 59-73 ◽  
Author(s):  
Yaxin Song ◽  
C. J. Hartwigsen ◽  
Lawrence A. Bergman ◽  
Alexander F. Vakakis

Sign in / Sign up

Export Citation Format

Share Document