Effect of Manufacturing Variability on Turbine Engine Performance: A Probabilistic Study

While turbine engine Original Equipment Manufacturers (OEMs) accumulated significant experience in the application of probabilistic methods (PM) and uncertainty quantification (UQ) methods to specific technical disciplines and engine components, experience with system-level PM applications has been limited. To demonstrate the feasibility and benefits of an integrated PM-based system, a numerical case study has been developed around the Honeywell turbine engine application. The case study uses experimental observations of engine performance such as horsepower and fuel flow from a population of engines. Due to manufacturing variability, there are unit-to-unit and supplier-to-supplier variations in compressor blade geometry. Blade inspection data are available for the characterization of these geometric variations, and CFD analysis can be linked to the engine performance model, so that the effect of blade geometry variation on system-level performance characteristics can be quantified. Other elements of the case study included the use of engine performance and blade geometry data to perform Bayesian updating of the model inputs, such as efficiency adders and turbine tip clearances. A probabilistic engine performance model was developed, system-level sensitivity analysis performed, and the predicted distribution of engine performance metrics was calibrated against the observed distributions. This paper describes the model development approach and key simulation results. The benefits of using PM and UQ methods in the system-level framework are discussed. This case study was developed under Defense Advanced Research Projects Agency (DARPA) funding which is gratefully acknowledged.

Rotordynamic Crack Diagnosis: Distinguishing Crack Depth and Location

The goal of this work is to establish a condition monitoring regimen capable of diagnosing the depth and location of a transverse fatigue crack in a rotordynamic system. The success of an on-line crack diagnosis regimen hinges on the accuracy of the crack model used. The model should account for the depth of the crack and the localization of the crack along the shaft. Negating the influence of crack location on system response ignores a crucial component of real cracks. Two gaping crack models are presented; the first simulates a finite-width manufactured notch, while the second models an open fatigue crack. An overhung rotordynamic system is modeled, imitating an available rotordynamic test rig. Four degree-of-freedom equations of motion for both crack models are presented and discussed, along with corresponding transfer matrix techniques. Free and forced response analyses are performed, with emphasis placed on results applicable to condition monitoring. It is demonstrated that two identifiers are necessary to diagnose the crack parameters: the 2X resonance frequency and the magnitude of the 2X component of the rotor angular response at resonance. First, a contour plot of the 2X resonant shaft speed versus crack depth and location is generated. The magnitude of the 2X component of the rotor’s angular response along the desired contour is obtained, narrowing the possible pairs of crack location/depth to either one or two possibilities. Practical aspects of the diagnosis procedure are then discussed.

Hybrid Analysis of Gas Annular Seals With Energy Equation

Annular seals are used in turbomachinery to reduce secondary flow between regions of high and low pressure. In a vibrating rotor system, the non-axisymmetric pressure field developed in the small clearance between the rotor and the seal generate reactionary forces that can affect the stability of the entire rotor system. Traditionally, two analyses have been used to study the fluid flow in seals, bulk-flow analysis and computational fluid dynamics (CFD). Bulk-flow methods are computational inexpensive, but solve simplified equations that rely on empirically derived coefficients and are moderately accurate. CFD analyses generally provide more accurate results than bulk-flow codes, but solution time can vary between days and weeks. For gas damper seals, these analyses have been developed with the assumption that the flow can be treated as isothermal. Some experimental studies show that the difference between the inlet and outlet temperature temperatures is less than 5% but initial CFD studies show that there can be a significant temperature change which can have an effect on the density field. Thus, a comprehensive analysis requires the solution of an energy equation. Recently, a new hybrid method that employs a CFD analysis for the base state, unperturbed flow and a bulk-flow analysis for the first order, perturbed flow has been developed. This method has shown to compare well with full CFD analysis and experimental data while being computationally efficient. In this study, the previously developed hybrid method is extended to include the effects of non-isothermal flow. The hybrid method with energy equation is then compared with the isothermal hybrid method and experimental data for several test cases of hole-pattern seals and the importance of the use of energy equation is studied.

Experimental Study on Lateral and Torsional Vibration of Cracked Rotor With Torsional Excitation

Crack failures in rotating machinery can result in catastrophic accidents, and they are are difficult to detect online. Condition monitoring is widely applied in field to detect changes of vibration, and form diagnostic features. However, effective features in vibration of the cracked rotor need more tests, especially validating the features with experiments. This work carried out an experimental study on cracked rotors in laboratory. The experiments are as following: (I) vibration of the rotor in normal condition is firstly tested, where lateral vibration and torsional vibration are measured; (II) torsional excitation is exerted on driven end of rotor system, and vibration characteristics of the rotor are tested; (III) cracked rotors are tested with transverse and slant cracks, respectively. With the measured signals, comparisons of vibrations in normal rotor and cracked rotors are carried out. The results show that, the transverse crack introduces more significant changes in 1X frequency and coupled frequency, while the slant crack employs larger changes in 2X frequency. And variation of phases of 1X frequency is presented. Also, the crack plays an impact on the torsional responses.

An Optimal Automated Zone Allocation Approach for Risk Assessment of Gas Turbine Engine Components

High-energy rotating components of gas turbine engines may contain rare material anomalies that can lead to uncontained engine failures. The Federal Aviation Administration and the aircraft engine industry have been developing enhanced life management methods to address the rare but significant threats posed by these anomalies. One of the outcomes of this effort has been a zone-based risk assessment methodology in which component fracture risk is estimated using groupings of elements called zones that are associated with 2D finite element (FE) stress and temperature models. Previous papers have presented processes for creation of zones either manually or via an automatic algorithm in which zones are assigned to each finite element in a component model. These processes may require significant human time and computer time. The focus of this paper is on the optimal allocation of multiple finite elements to zones that minimizes the total number of zones required to compute the fracture risk of a component. An algorithm is described that uses a relatively coarse response surface method to estimate the conditional risk value at each node in a finite element model. Zones are initially defined for each finite element in the model, and the algorithm identifies and merges zones based on minimizing the influence on component risk. The process continues until all of the zones have been merged into a single zone. The zone sequence is applied in reverse order to identify the minimum number of zones that satisfies component target risk or convergence threshold constraints. This solution provides the optimal allocation of finite elements to zones. The algorithm is demonstrated for a representative gas turbine engine component. The approach significantly improves the computational efficiency of the zone-based risk analysis process.

On Conservative CTOD Estimation Using Far Crack Deformation Field

In general engineering practice, crack tip opening displacement (CTOD) is very convenient approach for prediction of the components fracture mechanics (FM) lifetime. FM lifetime calculations are defined very well in industry and the lifetime prediction methods based on the CTOD resolve linear and nonlinear material behavior for monotonic and cyclic responses. The experiments confirm that under plasticity conditions the crack tip blunts for small scale or large scale yielding while, crack flanks open against each other only under elastic conditions. However, the CTOD application requires a very fine mesh in order to predict a crack tip deformation in reliable manner. Therefore, much more engineering work have to be involved in fine FE modeling. The crack tip flank deformation is crucial parameter responsible for reliable prediction of the nonlinear energy release rate, which is obtained from Hutchinson-Rice-Rosengren solution and the Shih rule. In accordance with design guidelines, the nonlinear energy release rate obtained from the CTOD must be evaluated conservatively to meet demands of RAM (Reliability, Availability and Maintainability). By using far crack deformation field, the paper proposes an engineering approach, which predicts the CTOD in a conservative manner under elastic-plastic conditions. This novel method is validated numerically by applying the well-known J-integral approach.

Structural Change Quantification in Rotor Systems Based on Measured Resonance and Antiresonance Frequencies

A structural change quantification methodology is proposed in which the magnitude and location of a structural alteration is identified experimentally in a rotor system. The resonance and antiresonance frequencies are captured from multiple frequency response functions and are compared with baseline data to extract frequency shifts due to these features. The resulting expression contains sufficient information to identify the dynamic characteristics of the rotor in both the frequency and spatial domains. A finite element model with carefully selected tunable parameters is iteratively adjusted using a numerical optimization algorithm to determine the source of the structural change. The methodology is experimentally demonstrated on a test rig with a laterally damaged rotor and the frequency response functions are acquired through utilization of magnetic actuators positioned near the ball bearings.

Damping Identification and its Comparison for Various Types of Blade Couplings

Regarding steam turbine blade vibrations, damping of blade as well as bladed disc mode shapes is one of the most important parameters in terms of steam turbine operation. A value of the parameter depends on properties of material used for manufacturing and construction elements of the blades and the discs such as blade roots, shrouds, tiebosses (snubbers) and dampers. This article deals with a comparison of damping of mode shapes for particular blade couplings and shows which methods are suitable for determination of the damping in individual cases. The whole identification procedure of the damping together with its specifics is also presented. At first, an identification technique of material damping ratio is introduced and its results are given for different materials. The material damping ratio is assessed as material strain dependent. Subsequently, damping ratio of bladed disc mode shapes under bladed disc rotation is identified taking into account two alternatives. The alternatives differ in such a way that blades have been free for the first time and then coupled with friction dampers. Outcomes presented in the article illustrate good agreement between damping ratio of bladed disc mode shapes with free blades and material used for manufacturing of the blades. On the other hand, damping ratio of bladed disc mode shapes with friction dampers is significantly different and strongly dependent on blade vibration amplitudes as well as nodal diameters of bladed disc mode shapes. Finally, nonlinear behavior of the bladed disc has been revealed along large blade vibration amplitudes and higher nodal diameters of the disc. The non-linear behavior manifests itself in such a way that values of natural frequencies of the disc have become dependent on blade vibration amplitudes.

Improvement of Compressor Blade Vibrations Spectral Analysis From Tip Timing Data: Aliasing Reduction

Tip-timing is a technique for measuring rotating blades vibrations in operation. Its concept exists since the early 70s but it has been more experimented in the last decade through improvements in hardware and software capabilities. It consists of a set of sensors mounted on stator casings that record blade passing times. Then, from this measurement, blade vibrations can be estimated. Tip-timing technique presents several advantages compared to usual mean of measurement: strain gages. Indeed, installation is easier, it is non-intrusive and all the blades can be monitored. However, resulting sampling depends on physical configuration i.e. number of sensors. In practice, as the number of sensors is limited, sampling rate is low in relation to the physical observed frequencies and do not respect the Shannon criterion. Thus, it generates important aliasing effects in spectrum, which makes the analysis difficult. In fact, such measurements aim to lead to pseudo-blind analysis, especially for asynchronous vibrations, when there is no hypothesis of underneath structural model. This main problem of aliasing is partially softened by using a minimum variance spectral estimator that iteratively reduces non-physical content in spectrum, but pseudo-blind analysis remains complicated. This paper presents a comparison of several methods based on different multisampling averaging concepts for reducing aliasing. Multisampling averaging consists in averaging spectrums from different sampling patterns of the same original signal, so that only stable physical content remains in the final spectrum. This study is presented on different industrial test cases of blade vibrations.