Refine
Year of publication
Document Type
Keywords
Institute
BTU
This paper describes preparation, execution and evaluation of a comprehensive bladed disk spin test series. At the example of an turbine impeller the effects of rotation and temperature are analyzed with special focus on mistuning and damping. The forced response is measured synchronously via 13 identical positioned strain gauges on each blade as well as via blade tip-timing. Subsequently it is possible to compare the results of both systems. During the test series rotational speed varies in the range from 10.000 up to 19.000 RPM. Simultaneously, the wheel is heated up to 820 K by an oven. A number of pre-selected natural frequencies, damping ratios and operating deflection shapes are evaluated and compared with respect to different rotational speeds and impeller temperatures.
Copyright © 2015 by ASME
To achieve reverse objectives in engine design, advanced modelling and analysis methods are among the key research technologies. In the presented work, a robust design optimization of a first stage high pressure turbine blade has been carried out. This blade derives from a current production of a Rolls-Royce aero engine.
The motivation of this work is to show that the methodology of robust design optimization can be applied to high pressure turbine blades. A fully automated workflow, which encapsulated the integral blade design and analysis process, has been used. The main workflow objective is a representative life value of the external surface of the blade. In addition, the workflow enables the engineering uses to consider sub objectives like mass, efficiency and life at critical locations of the blade. These can also be taken into account in the multi-objective robust design optimization. This research also focuses on the use of surrogate models, with attention to the delivery of a physically correct result. For this purpose, the validation of the applied methods has a huge significance and a toolbox was created to generate and evaluate the quality of the surrogate models.
In the present case sixteen geometry parameters were considered. In order to show that this methodology is not limited to geometry variation, parameters for material specification and for boundary conditions were varied in addition. The surrogate model was trained by the workflow generated DoE-data and could be used for different kinds of optimization.
As a conclusion, it has been demonstrated that the methodology can be used for the engineering design process of turbine blades, while delivering physically correct results. The different techniques for surrogate modelling were examined and compared. With the help of these surrogate models, an optimization of life, mass and efficiency with 22.5 million evaluations was possible. Finally, an overview of the methodology for the case of a real world turbine blade could be given, and an improved blade in the sense of multi-objective robust design was found.
Copyright © 2015 by Rolls-Royce Deutschland Ltd & Co KG
The forced response of an E3E-type high pressure compressor blisk front rotor is analyzed with regard to intentional mistuning and its robustness towards additional random mistuning. Both a chosen alternating mistuning pattern and artificial mistuning patterns optimized concerning the forced response are considered. Focusing on three different blade modes, subset of nominal system mode-based reduced order models are employed to compute the forced response. The disk remains unchanged while the Young’s modulus of each blade is used to define the particular mistuning pattern. The well established aerodynamic influence coefficient technique is employed to model aeroelastic coupling and hence to consider the strongly mode- and inter blade phase angle-dependent aerodynamic damping contribution.
It has been found that a reduction of the maximum forced response beyond that of the tuned reference can be achieved for particular mistuning patterns and all modes considered. This implies an exciting engine order which would cause a low nodal diameter mode in case of a tuned blisk. At best a nearly 50% reduction of maximum response magnitudes is computed for the fundamental bending mode and large mistuning. The solution proved to be robust towards additional random mistuning of reasonable magnitude, which is of particular interest with regard to a potential technical realization. In case of small mistuning as assumed for the first torsion and the longitudinal bending mode the advantage of achieving response magnitudes beyond the tuned reference gets lost indeed, if random mistuning is superimposed. However, mostly a lower response level is calculated compared to responses obtained from models adjusted to mistuning determined by experiment.
Copyright © 2015 by ASME
With increasing demands for reliability of modern turbomachinery blades the quantification of uncertainty and its impact on the designed product has become an important part of the development process. This paper aims to contribute to an improved approximation of expected vibration amplitudes of a mistuned rotor assembly under certain assumptions on the probability distribution of the blade’s natural frequencies. A previously widely used lumped mass model is employed to represent the vibrational behavior of a cyclic symmetric structure. Aerodynamic coupling of the blades is considered based on the concept of influence coefficients leading to individual damping of the traveling wave modes. The natural frequencies of individual rotor blades are assumed to be normal distributed and the required variance could be estimated due to experiences with the applied manufacturing process. Under these conditions it is possible to derive the probability distribution of the off-diagonal terms in the mistuned equations of motions, that are responsible for the coupling of different circumferential modes. Knowing these distributions recent limits on the maximum attainable mistuned vibration amplitude are improved. The improvement is achieved due to the fact, that the maximum amplification depends on the mistuning strength. This improved limit can be used in the development process, as it could partly replace probabilistic studies with surrogate models of reduced order. The obtained results are verified with numerical simulations of the underlying structural model with random mistuning patterns based on a normal distribution of individual blade frequencies.
Copyright © 2015 by ASME
In the drive for higher cycle efficiencies in gas turbine engines, turbine blades are seeing an increasingly high heat load. This in turn demands improvements in the internal cooling system and a better understanding of both the level and distribution of the internal heat-transfer. A typical approach to enhance the internal cooling of the turbine blade is by casting angled ‘low blockage’ ribs on the walls of the cooling channels. The objective of the present paper is to determine the detailed Nusselt number distribution in rectangular internal channels with ribs. This knowledge can be used to guide the overall design e.g. to achieve high levels of heat-transfer where required. The effects of rotation as well as the interaction effects of the position and direction of ribs on opposite walls of the cooling channel have been investigated.
Numerical calculations have been carried out using the commercial CFD code Fluent to investigate the local Nusselt number enhancement factor in rectangular ducts of different aspect ratios (0.5, 1 and 2) which have 45° or 90° angled ribs located on two opposite walls. This has been studied for different Rotation number Ro (0–0.45) and with a Reynolds number >30000.
The first series of studies has been carried out with the same experimental setup as by Han [1]. The geometry was slightly changed to avoid the effect of high heat transfer at the entry.
This study identifies important vortical structures, which are dependent on the direction and the position of the ribs. This has a profound effect on the distribution of heat-transfer within the passage. It is shown that the two smooth walls of the duct have different average Nusselt number ratio Nu/NuFD enhancement depending on the rib angle.
In addition, based on numerical investigations, simple correlations have been developed for the rotational influence of the internal Nusselt number distribution. A major finding is that the effect of rotation is dominant for low aspect ratio channels and the local enhancement due to the rib position and angle is more dominant for high aspect ratio channels.
The forced response of an E3E-type HPC-blisk front rotor is analyzed with regard to varying mistuning and the consideration of the fluid-structure interaction (FSI). For that purpose, a reduced order model is used in which the disk remains unchanged and mechanical properties of the blades namely stiffness and damping are adjusted to measured as well as intentional blade frequency mistuning distributions. The aerodynamic influence coefficient technique is employed to model the aeroelastics. Depending on the blade mode, the exciting engine order and aerodynamic influences it is sought for the worst mistuning distributions with respect to the maximum blade displacement based on optimization analyses. Genetic algorithms using blade alone frequencies as design variables are applied. The validity of the Whitehead-limit is assessed in this context. In particular, the question is addressed if and how far aeroelastic effects, mainly caused by aerodynamic damping, combined with mistuning can even cause a reduction of the forced response compared to the ideally tuned blisk. It is shown that the strong dependence of the aerodynamic damping on the inter-blade phase angle is the main driver for a possible response attenuation considering the fundamental as well as a higher blade mode. Furthermore, the differences to the blisk vibration response without a consideration of the flow and an increase of the disk’s stiffness are discussed. Closing, the influence of pure damping mistuning is analyzed again using optimization.
The application of high pressure compressor (HPC) rotors manufactured as blisk (Blade Integrated Disk) is ever-expanding in modern jet engine designs. Despite the major advantages of less mass and higher efficiency, the most challenging problem is lower mechanical damping due to the loss of damping between blades root’s and the disk. Mistuning is induced by material inhomogeneities, manufacturing tolerances or wear during use and leads to amplitude magnification and mode localization. From the experimental point of view mistuning can be evaluated via experimental vibration analysis in terms of frequency deviations. Furthermore optical measurements can be evaluated in terms of geometrical deviations between the real and designed geometry. From the structural point of view a mistuned blisk model can be obtained by morphing the nodes of the geometrical tuned FE model or by performing blade individual stiffness mistuning due to modification of Young’s modulus. The following work is focused on the numerical prediction of mistuned blisk vibrations. Therefore, the research blisk of the 4 stage research compressor, manufactured as job-production, is analyzed. For this research blisk optical measurement data as well as experimentally obtained frequency patterns are available. In a first part mistuning identification in terms of experimental vibration analysis and Proper Orthogonal Decomposition of the geometrical deviations is presented. In a second part mistuning modeling in terms of stiffness mistuning and geometrical mistuning is applied to the tuned FE-model and the numerical results are evaluated against experimental data regarding accuracy. Furthermore, the impact of geometrical deviations on mistuning is analyzed.
Copyright © 2015 by Rolls-Royce Deutschland Ltd & Co KG
Die Verbesserung von Kalibriermethoden für Hochdruckverdichterstufen in ausgeführter Integralbauweise verlangt von den Ingenieuren ein detaillierteres Verständnis der strukturdynamischen Zusammenhänge dieses komplexen Bauteils. Als Ausgangsbasis werden zunächst für das perfekte Normdesign bei Ausnutzung der vorliegenden zyklischen Rotationssymmetrie numerische Ergebnisse vorgestellt. Dabei stehen bezüglich des Schwingverhaltens die Kopplungseffekte zwischen Scheibe und Schaufeln im Mittelpunkt, wobei experimentelle Untersu-chungen am realen Bauteil ergänzend hinzugezogen werden. Zukünftig wird in einem weiteren Schritt im Rahmen erheblich aufwendigerer Berechnungen mit der Berücksichtigung von Imperfektionen hinsichtlich der filigranen Schaufeln die Rotationssymmetrie aufgegeben und somit die Basis für eine gezielte Beeinflussung des Schwingverhaltens infolge von Mistuning - Effekten geschaffen.
Considerations About the Necessary Mesh Density of Bearings in Detailed Finite Element Models
(2014)
In order to predict properly the behavior of very complex mechanical models, much focus needs to be done on the level of model details as well as the techniques used to predict component movements based on mechanical and thermal loads. In particular bearings and joints play a significant role in big model assemblies such as aero-engines where most of the structures are modeled by solid elements rather than using traditional shell or beam elements. This approach is very complex and expensive in terms of computational effort however can lead to increased simulation accuracy when done properly.
This publication will answer questions related to the modeling techniques required for detailed bearings (not only rotor bearings) and joints of an aero-engine application where these components are not idealized or simplified, i.e. the ball in a ball joint is meshed by solid elements. The investigation will discuss the question about the number of elements needed to mesh a ball joint in such a way that a correct rotation is possible. If the mesh is too coarse, a rotation would be impossible due to the faceting of the rotating parts. It will be shown that a proper rotation depends on a correct relation between the radius of the inner ball, the gap size of the bearing and the number of nodes in the circumferential direction. For this relation an analytical formula has been derived and successfully tested at some examples. Furthermore, some finite element codes offer the possibility of smoothed contact surfaces. These options improve the situation but they might cause issues too. For this reason these smoothed contact options are also considered under the aspects mentioned above.
Copyright © 2014 by Rolls-Royce Deutschland Ltd & Co KG
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.
Modal and Aeroelastic Analysis of a Mistuned Compressor Blisk Using an Equivalent Blisk Model
(2011)
The oil and gas, chemical, and process industries employ centrifugal compressors for a wide range of applications. Due to this, the conditions under which centrifugal compressors have to operate, vary significantly from case to case. Gas pipeline compressors, for example, may feature discharge pressures well over 100 bar. During the last decades, comprehensive research was conducted on the impact of high pressure operating conditions on the vibrational behavior of centrifugal compressors. Nowadays, it is well-known that an increase in gas pressure levels leads to a more pronounced interaction between the side cavities and the impeller, which results in a frequency shift of the acoustic and structural modes.
For the safe operation of compressors, it is necessary to predict these coupled natural frequencies accurately. The state-of-the-art approach to achieve this objective is the finite element method. While this technique provides high-quality results, it incurs high computational costs and is, therefore, time-consuming.
The authors of the current paper propose a generalized model to overcome this challenge. It uses the uncoupled modes of the impeller and side cavities in a modal superposition to approximate the coupled system's natural frequencies. In this way, the intended design geometries are considered while reducing the computational effort significantly. In a numerical study, the generalized model is applied to different systems of increasing complexity, and the results are compared to a finite element analysis. Finally, the paper concludes with a discussion of the limitations and benefits of all employed numerical methods.
The manufacturing geometrical variability is a source of uncertainty, which cannot be avoided in the realization of machinery components. Deviations of a part geometry from its nominal design are inevitably present due to the manufacturing process. In the case of the aeroelastic forced response problem within axial compressors, these uncertainties may affect the vibration characteristics. For this reason, the impact of geometrical uncertainties due to the manufacturing process onto the modal forcing of axial compressor blades is investigated in this study. The research focuses on the vibrational behavior of an axial compressor rotor blisk. In particular, the amplitude of the forces acting as a source of excitation on the vibrating blades is studied. The geometrical variability of the upstream stator is investigated as input uncertainty. The variability is modeled starting from a series of optical surface scans. A stochastic model is created to represent the measured manufacturing geometrical deviations from the nominal model. A data reduction methodology is proposed in order to represent the uncertainty with a minimal set of variables. The manufacturing geometrical variability model allows to represent the input uncertainty and probabilistically evaluate its impact on the aeroelastic problem. An uncertainty quantification is performed in order to evaluate the resulting variability on the modal forcing acting on the vibrating rotor blades. Of particular interest is the possible rise of low engine orders due to the mistuned flow field along the annulus. A reconstruction algorithm allows the representation of the variability during one rotor revolution. The uncertainty on low harmonics of the modal rotor forcing can be therefore identified and quantified.
Objective of this paper is to analyze the consequences of borescope blending repairs on the aeroelastic behavior of a modern high pressure compressor (HPC) blisk. To investigate the blending consequences in terms of aerodynamic damping and forcing changes, a generic blending of a rotor blade is modeled. Steady-state flow parameters like total pressure ratio, polytropic efficiency, and the loss coefficient are compared. Furthermore, aerodynamic damping is computed utilizing the aerodynamic influence coefficient (AIC) approach for both geometries. Results are confirmed by single passage flutter (SPF) simulations for specific interblade phase angles (IBPA) of interest. Finally, a unidirectional forced response analysis for the nominal and the blended rotor is conducted to determine the aerodynamic force exciting the blade motion. The frequency content as well as the forcing amplitudes is obtained from Fourier transformation of the forcing signal. As a result of the present analysis, the change of the blade vibration amplitude is computed.
The manufacturing geometrical variability is a source of uncertainty, which cannot be avoided in the realization of a machinery. Deviations of a component geometry from its nominal design are inevitably present due to the manufacturing process. In the case of the aeroelastic forced response problem within axial compressors, these uncertainties may affect the vibration characteristics. For this reason, the impact of geometrical uncertainties due to the manufacturing process onto the modal forcing of axial compressor blades is investigated in this study.
The research focuses on the vibrational behavior of an axial compressor rotor blisk (blade-integrated disk) and in particular the amplitude of the forces acting as source of excitation on the vibrating blades (modal forcing). Within this context, the geometry of the upstream stator plays an important role as in general the main harmonics of the rotor excitation forces are produced by its wake. Therefore, small variations of the upstream stators geometries, such as the ones caused by the manufacturing process, may affect the resulting forcing. In particular, the geometrical variability of the upstream stator implies that the hypothesis of a cyclic-symmetrical flow is no longer valid. This may cause the introduction of lower harmonic components in the modal forces, generally referred to as Low Engine Orders (LEO).
The geometrical variability is modelled starting from a series of optical surface scans. A set of optical measurements of manufactured stator blades originating from the same nominal design constitutes the baseline dataset on which the geometrical model is built. The measured blades as well as the relative nominal geometry are parametrized to describe the individual blades surfaces. The parameterization is accomplished by slicing the surfaces in radial sections and describing each of these with a set of NACA-like parameters [1]. The measured geometrical deviations from the nominal model can therefore be described as an offset of such parameters. A reduced representation of the variables representing the input uncertainty (noise variables) is obtained via Principal Components Analysis. Afterwards a sampling on the reduced noise variables domain can be done to represent the modelled uncertainty and perform an Uncertainty Quantification (UQ) on the relative quantities of interest, in this case the modal forcing.
The computation of the modal forcing is done through a CFD solver, computing the unsteady flow field around the rotor blades. The domain considered in this case is a 1.5 stage of the axial compressor, including the rotor and the up- and down-stream stators. The solutions are initialized from a validated steady state solution of the considered compressor rig. The time-dependent pressure field calculated on the rotor blades is projected onto the relative vibrational mode shapes of interests (from structural modal analyzes). The resulting forces are analyzed by means of their spectrum, evaluating the amplitudes for the present engine orders (higher harmonics of the shaft mechanical speed).
The UQ uses Monte Carlo methods to evaluate the impact of the geometrical variability onto the modal forcing. The modelled uncertainty on the geometries is introduced into the CFD solver to compute the deviations on the quantities of interest. A reconstruction of the forces acting on the rotor during one revolution is obtained. This allows to evaluate the uncertainty on the present engine orders as well as the possible rise of LEO for the rotor blades in presence of a mistuned upstream stator.
[1]: Lange A., Vogeler K., Gümmer V., Schrapp H. and Clemen C. (2009). “Introduction of a Parameter Based Compressor Blade Model for Considering Measured Geometry Uncertainties in Numerical Simulation.” Proceedings of ASME Turbo Expo. GT2009-59937
Simulation fluidgedämpfter Strukturschwingungen mittels partitioniertem Kopplungssatz via MpCCI
(2006)
Moderne Verbrennungskraftmaschinen müssen ein stetig wachsendes Anforderungsprofil in Bezug auf Wirtschaftlichkeit, Leistung und Umweltfreundlichkeit erfüllen. In diesem Zusammenhang hat die Turboaufladung von Verbrennungsmotoren an Bedeutung gewonnen. Bei Turboladern kleiner und mittlerer Baugröße, deren Turbinen einen Durchmesser zwischen 30 und 250 mm aufweisen, kommen vornehmlich gegossene Laufräder zum Einsatz. Am Institut für Verkehrstechnik der Brandenburgischen Technischen Universität Cottbus-Senftenberg wurde im Rahmen eines FVV-Forschungsvorhabens der Einfluss charakteristischer Betriebsgrößen eines Turboladers auf das strukturdynamische Verhalten solcher Radialturbinenräder untersucht.
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 can lead to flow separation in some circumstances. Consequently, there is the risk of flutter in principle, particularly at nominal speed under part load conditions. For this reason, intentional mistuning is employed by the manufacturer 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 since only two different blade designs are allowed. In this sense, two different series of blades have been made. However, it is well known that 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 but unwanted deviations is considered in numerical simulations. Moreover, the strong dependence of blade frequencies on the speed is taken into account since it significantly attenuates the blade to blade frequency difference in this particular case.
Within an academic study the turbine wheel is modelled as blade integrated disk in order to demonstrate fundamental effects of intentional mistuning on flutter susceptibility and forced response. For that purpose, reduced order models are built up by using the subset of nominal system mode approach introduced by Yang and Griffin [1], which conveniently allows for taking into account both differing mistuning patterns and the impact of aeroelastic interaction. Focusing on the first flap mode it could be shown that a mitigation of flutter susceptibility is achieved by prescribing alternate mistuning, which indeed affects 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. Again only two different blade designs are admitted. Finally, mistuning patterns could be identified causing a tremendous increase of aerodynamic damping ratios. The robustness of the solutions found could be proved by superimposing additional random mistuning. Another study is focused on the impact of mistuning strength.
Further analyses are addressing the forced response at part speed conditions, where different resonance crossings are becoming apparent in the Campbell plot. An increase of the forced response compared to the tuned counterpart is partly unpreventable because of unfavorable aerodynamic damping curves. Independently, the maximum forced response has to be limited also in case of applying large intentional mistuning.
[1] Yang, M. T., Griffin, J. H., „A Reduced-Order model of Mistuning Using a Subset of Nominal System Modes“. J Eng Gas Turb Power, 123, pp. 893-900 (2001).
This paper aims at contributing to a better understanding of the effect of Tyler–Sofrin Modes (TSMs) on forced vibration responses by analysing a 4.5-stage research axial compressor rig. The first part starts with a brief review of the involved physical mechanisms and necessary prerequisites for the generation of TSMs in multistage engines. This review is supported by unsteady CFD simulations of a quasi 2D section of the studied engine. It is shown that the amplitude increasing effect due to mistuning can be further amplified by the presence of TSMs. Furthermore, the sensitivity with respect to the structural coupling of the blades and the damping as well as the shape of the expected envelope is analysed.
The second part deals with the Rotor 2 blisk of the research compressor rig. The resonance of a higher blade mode with the engine order of the upstream stator is studied in two different flow conditions realised by different variable stator vane (VSV) schedules which allows to separate the influence of TSMs from the impact of mistuning. A subset of nominal system modes representation of the rotor is used to describe its mistuned vibration behaviour, and unsteady CFD simulations are used to characterise the present strength of the TSMs in the particular operating conditions. Measured maximum amplitude vs blade pattern and frequency response functions are compared against the predictions of the aeromechanical models in order to assess the strength of the TSMs as well as its influence on vibration levels.
In order to prepare an advanced 4-stage high-pressure compressor rig test campaign, details regarding both accomplishment and analysis of preliminary experiments are provided in this paper. The superior objective of the research project is to contribute to a reliable but simultaneously less conservative design of future high pressure blade integrated disks (blisk). It is planned to achieve trend-setting advances based on a close combination of both numerical and experimental analyses. The analyses are focused on the second rotor of this research compressor, which is the only one being manufactured as blisk. The comprehensive test program is addressing both surge and forced response analyses e.g. caused by low engine order excitation. Among others the interaction of aeroelastics and blade mistuning is demanding attention in this regard. That is why structural models are needed, allowing for an accurate forced response prediction close to reality. Furthermore, these models are required to support the assessment of blade tip timing (BTT) data gathered in the rig tests and strain gauge (s/g) data as well. To gain the maximum information regarding the correlation between BTT data, s/g-data and pressure gauge data, every blade of the second stage rotor (28 blades) is applied with s/g. However, it is well known that s/g on blades can contribute additional mistuning that had to be considered upon updating structural models.
Due to the relevance of mistuning, efforts are made for its accurate experimental determination. Blade-by-blade impact tests according to a patented approach are used for this purpose. From the research point of view, it is most interesting to determine both the effect s/g-instrumentation and assembling the compressor stages on blade frequency mistuning. That is why experimental mistuning tests carried out immediately after manufacturing the blisk are repeated twice, namely, after s/g instrumentation and after assembling. To complete the pre-test program, the pure mechanical damping and modal damping ratios dependent on the ambient pressure are experimentally determined inside a pressure vessel. Subsequently the mistuning data gained before is used for updating subset of nominal system mode (SNM) models. Aerodynamic influence coefficients (AICs) are implemented to take aeroelastic interaction into account for forced response analyses. Within a comparison of different models, it is shown for the fundamental flap mode (1F) that the s/g instrumentation significantly affects the forced response, whereas the impact of assembling the compressor plays a minor role.
Especially for rotating structures like bearings non-smooth contact surfaces, as they appear in classical finite element modeling, lead to various problems during the analysis involving mesh interlocking and spurious oscillations in contact forces. In order to eliminate these issues, very fine meshes and additional smoothing strategies are employed in case of the finite element method (FEM).
But also Non-Uniform Rational B-splines (NURBS) based isogeometric analysis (IGA) can be very useful for the contact analysis due to the inherent higher order continuity of NURBS basis functions. In this contribution, the contact behavior using classical FEA and IGA approaches is studied by means of an example of a pendulum under gravitational load. In addition, a more practical example of a coupled IGA–FEM problem with a cylindrical roller bearing is also reported in this paper. This research study of contact analysis has been carried out for the above mentioned examples using LS-DYNA and illustrates that contact surfaces of coarsely meshed geometry lock the rotation of the parts in case of classical FEM. On the contrary, when the contact surface is represented by NURBS elements it allows the rotation of
the parts and effectively alleviates the contact force oscillation.
Prediction of Geometrically Induced Localization Effects Using a Subset of Nominal System Modes
(2019)
It is widely known that the vibration characteristics of blade integrated discs can dramatically change in the presence of manufacturing tolerances and wear. In this context, an increasing number of publications discuss the influence of the geometrical variability of blades on phenomena like frequency splitting and mode localization. This contribution is investigating the validity of a stiffness modified reduced order model for predicting the modal parameters of a geometrically mistuned compressor stage. In detail, the natural frequencies and mode shapes, as well as the corresponding mistuning patterns, are experimentally determined for an exemplary rotor. Furthermore, a blue light fringe projector is used to identify the geometrical differences between the actual rotor and the nominal blisk design. With the help of these digitization results, a realistic finite element model of the whole compressor stage is generated. Beyond that, a reduced order model is implemented based on the nominal design intention. Finally, the numerical predictions of the geometrically updated finite element model and the stiffness modified reduced order model are compared to the vibration measurement results. The investigation is completed by pointing out the benefits and limitations of the SNM-approach in the context of geometrically induced mistuning effects.
A novel approach for the identification of tuned aerodynamic damping, system frequencies, forcing and mistuning has been introduced in the first part of this paper. It is based on the forced response equations of motion for a blade dominated mode family. A least squares formulation allows to identify the system’s parameters directly from measured frequency response functions (FRFs) of all blades recorded during a sweep through a resonance. The second part has dealt with its modification and application to experimental modal analyses of blisks at rest.
This 3rd part aims at presenting the application of the approach to blade tip timing (BTT) data acquired in rig tests. Therefore, blisk rotors of two different engines are studied: a single stage fan rig and a 4.5 stage high pressure compressor (HPC) rig. The rig test campaign of the fan blisk included also an intentional mistuning experiment that allows to study the performance of the identification approach for a similar rotor with two different mistuning levels. It is demonstrated that the approach can identify aerodynamic damping curves, system frequencies, mistuning pattern and forced travelling wave modes (TWMs) from state of the art BTT data monitored during rig or engine tests. All derived mistuning patterns could be verified with reference measurements at standstill. The derived aerodynamic damping curves and system frequencies show a reasonable agreement with simulations. For the HPC case a multitude of excited TWMs could be identified which also lines up with previous simulations.
An axial turbine blisk for turbocharger applications is analyzed with respect to the effect of intentional mistuning on the forced response. Originally, the intentional mistuning pattern has been designed by employing a genetic algorithm optimization in order to reduce the forced response caused by low engine order excitation (LEO) of the fundamental flap mode. The solution found has been implemented in a prototype of that blisk. For the purpose of comparison, a second reference blisk has been manufactured without intentional mistuning. The actual mistuning distributions of the blisks have been identified by employing blade-by-blade impact testing. Alternatively, a new inverse approach has been employed, which is based on a least squares formulation and benefits from less experimental effort. Based on the information gained by the aforementioned testing procedures, subset of nominal systems (SNM)-models have been updated, which allow for considering the aeroelastic coupling by means of aerodynamic influence coefficients (AIC). Despite of small but unavoidable deviations from the design intention it could be proved within numerical simulations that the intended 70 per cent reduction of the maximum forced response is nevertheless achieved.
In addition, the paper is addressing the effect of the aforementioned intentional mistuning pattern on a higher mode, which is relevant for the durability as well. Hence, new SNM-models have to be updated in order to calculate the forced response due to EO-excitation caused by the nozzle guide vane. Although the original mistuning pattern has been optimized solely for reducing the forced response of the fundamental flap mode, it hardly affects the higher mode forced response in a negative manner.
Comparative Evaluation of Isogeometric Analysis and Classical FEM with Regard to Contact Anaylsis
(2019)
Isogeometric analysis represents a newly developed technique that offers the application of Computer Aided Designs (CAD) concept of Non-uniform Rational B-Splines (NURBS) tool to describe the geometry of the computational domain. The simplified transition of CAD models into the computational domain eliminates the problems arising from the geometrical discontinuities induced by the faceted approximation of the mesh. Moreover, numerical analysis directly on NURBS objects significantly reduces the design-to-analysis time compared to traditional FEA approach. In the field of contact mechanics, when finite elements are applied to geometry with curved surfaces, the result is a non-smooth geometrical representation of interface surfaces which may lead to mesh interlocking, high jumps and spurious oscillations in contact forces. To eliminate these issues, various surface smoothening strategies are to be employed in case of FEM. Isogeometric based analysis alleviates these issues without employing any additional smoothening strategy due to inherent higher order continuity of NURBS basis functions and much more accurate results are obtained compared to conventional FE approach. In the current study, LS-DYNA is used to demonstrate the capabilities and advantage of an isogeometric analysis though an example of pendulum under gravitational load. The numerical simulation results are analytically validated and the comparison of NURBS surfaces with faceted surfaces is carried out to investigate the accuracy.
Turbocharging is known to be a well-established technology for an engine's efficiency and power output by forcing extra compressed air into the combustion chamber. The centrifugal loads, necessary flow deflections, unsteady pressure fluctuations, and structural temperature gradients put a high strain on rotating components. Additionally, those components are prone to high-cycle fatigue. The Chair of Structural Mechanics and Vehicle Vibrational Technology at the BTU Cottbus-Senftenberg investigated the impact of manufacturing tolerances on the vibrational behavior of several turbine and compressor impellers. Finally, it is shown that intentional mistuning can lead to significantly lower stresses.