NASTRAN Analysis of a Turbine Blade and Comparison with Test and Field Data

Similar documents
EFFECT OF TAPER AND TWISTED BLADE IN STEAM TURBINES

FATIGUE LIFE PREDICTION OF TURBOMACHINE BLADING

Structural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian

Deflections and Strains in Cracked Shafts due to Rotating Loads: A Numerical and Experimental Analysis

Frequency Failure Investigation on Shrouded Steam Turbine Blade through Dynamic Analysis

Modal Analysis: What it is and is not Gerrit Visser

Static and Dynamic Analysis of mm Steel Last Stage Blade for Steam Turbine

Failure Analysis of Last Stage Low Pressure Steam Turbine Blade

Vibration Analysis Of Rotating Cantilever Beam By Simulating Foreign Object Damage Nandish R.V 1, Somashekar H 2, S. Paul Vizhian 3, S.

Blade Group Fatigue Life Calculation under Resonant Stresses

Deflections and Strains in Cracked Shafts Due to Rotating Loads: A Numerical and Experimental Analysis

Transient Vibration Prediction for Rotors on Ball Bearings Using Load-Dependent Nonlinear Bearing Stiffness

Modal and Harmonic analysis of L.P. Turbine of a small Turbo- Fan engine using Finite Element Method

Some effects of large blade deflections on aeroelastic stability

Nonlinear Rolling Element Bearings in MADYN 2000 Version 4.3

Numerical Study on Performance of Curved Wind Turbine Blade for Loads Reduction

Response Spectrum Analysis Shock and Seismic. FEMAP & NX Nastran

SAMCEF For ROTORS. Chapter 1 : Physical Aspects of rotor dynamics. This document is the property of SAMTECH S.A. MEF A, Page 1

AEROELASTICITY IN AXIAL FLOW TURBOMACHINES

Vibration Analysis. with SOLIDWORKS Simulation 2018 SDC. Paul M. Kurowski. Better Textbooks. Lower Prices.

Article. Kamlesh Purohit and Manish Bhandari*

Structural Analysis of Wind Turbine Blades

NX Nastran 10. Rotor Dynamics User s Guide

What is the crack propagation rate for 7075-T6 aluminium alloy.

Vibration characteristics analysis of a centrifugal impeller

Propagation of Uncertainty in Stress Estimation in Turbine Engine Blades

DYNAMIC STRESS MEASUREMENT OF CENTRIFUGAL COMPRESSOR IMPELLER AND STUDY FOR STRENGTH CRITERIA BASED ON CORRELATION BY UNSTEADY CFD

Analytical Strip Method for Thin Isotropic Cylindrical Shells

International Journal on Emerging Technologies (Special Issue on NCRIET-2015) 6(2): (2015)

Review of the Master SN Neuber Rule in the ASME Division 2 Rewrite Project

PLEASURE VESSEL VIBRATION AND NOISE FINITE ELEMENT ANALYSIS

STRAIN ASSESSMENT USFOS

Prediction of Propeller Blade Stress Distribution Through FEA

Centrifugal Compressor Root Cause Analysis. (Pulsation and Vibration Issues)

Fatigue Life Prediction of Turbomachine Blading

Workshop 8. Lateral Buckling

4.5 The framework element stiffness matrix

A STUDY ON THE FRACTURE A SIROCCO FAN IMPELLER

202 Index. failure, 26 field equation, 122 force, 1

Aeroelastic effects of large blade deflections for wind turbines

DETC98/PTG-5788 VIBRO-ACOUSTIC STUDIES OF TRANSMISSION CASING STRUCTURES

FE-Analysis of Stringer-to-floor-beam Connections in Riveted Railway Bridges

Improving the Accuracy of Dynamic Vibration Fatigue Simulation

Basics of Finite Element Analysis. Strength of Materials, Solid Mechanics

CHAPTER 4 DESIGN AND ANALYSIS OF CANTILEVER BEAM ELECTROSTATIC ACTUATORS

Hardened Concrete. Lecture No. 16

New Representation of Bearings in LS-DYNA

0000. Finite element modeling of a wind turbine blade

FREE VIBRATIONS ANALYSIS OF SHROUDED BLADED DISCS WITH ONE LOOSE BLADE

CHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES

COPYRIGHTED MATERIAL. Index

STRUCTURAL PROPERTIES OF HELICOPTER ROTOR BLADE SECTIONS

D && 9.0 DYNAMIC ANALYSIS

A beam reduction method for wing aeroelastic design optimisation with detailed stress constraints

Optimum Height of Plate Stiffener under Pressure Effect

ROLLER BEARING FAILURES IN REDUCTION GEAR CAUSED BY INADEQUATE DAMPING BY ELASTIC COUPLINGS FOR LOW ORDER EXCITATIONS

3D Finite Element Modeling and Vibration Analysis of Gas Turbine Structural Elements

Figure 1 Lifting Lug Geometry with Weld

Quintic beam closed form matrices (revised 2/21, 2/23/12) General elastic beam with an elastic foundation

CRITERIA FOR SELECTION OF FEM MODELS.

MITOCW MITRES2_002S10nonlinear_lec20_300k-mp4

8/1/2009. CAE 7962 Presentation

RESONANCE IDENTIFICATION FOR IMPELLERS

AERO 214. Lab II. Measurement of elastic moduli using bending of beams and torsion of bars

4. SHAFTS. A shaft is an element used to transmit power and torque, and it can support

Investigation on fluid added mass effect in the modal response of a pump-turbine runner

Part D: Frames and Plates

Finite Element Analysis of Compression of Thin, High Modulus, Cylindrical Shells with Low-Modulus Core

Effects of Structural Forces on the Dynamic Performance of High Speed Rotating Impellers.

ARTICLE A-8000 STRESSES IN PERFORATED FLAT PLATES

HISTORY OF EVOLUTION, PROGRESS AND APPLICATION OF SAFE DIAGRAM FOR TUNED AND MISTUNED SYSTEMS. Murari Singh, Ph.D

Methods of Analysis. Force or Flexibility Method

DYNAMIC CHARACTERSTIC ESTIMATION OF STRUCTURAL MATERIALS BY MODAL ANALYSIS USING ANSYS

2. Determine the deflection at C of the beam given in fig below. Use principal of virtual work. W L/2 B A L C

Presented By: EAS 6939 Aerospace Structural Composites

#SEU16. FEA in Solid Edge and FEMAP Mark Sherman

Applied Modal Analysis of Wind Turbine Blades

Numerical Study on Performance of Innovative Wind Turbine Blade for Load Reduction

AERSYS KNOWLEDGE UNIT

UNIT- I Thin plate theory, Structural Instability:

MITOCW MITRES2_002S10linear_lec07_300k-mp4

THERMOWELL VIBRATION INVESTIGATION AND ANALYSIS

Prediction of Axial Compressor Blade Vibration by Modelling Fluid-Structure Interaction

Design of the blisk of an aircraft turbojet engine and verification of its resonance free operation

Torsion of Solid Sections. Introduction

History of WRC 107. Using WRC 107 and NozzlePRO FEA. Presented by: Ray Delaforce

(Refer Slide Time: 1: 19)

INTERMEDIATE MECHANICS OF DEFORMABLE BODIES (58:150/51:151/53:140) Fall 2003

FORCED RESPONSE COMPUTATION FOR BLADED DISKS INDUSTRIAL PRACTICES AND ADVANCED METHODS

Chapter 2: Rigid Bar Supported by Two Buckled Struts under Axial, Harmonic, Displacement Excitation..14

Fatigue Crack Analysis on the Bracket of Sanding Nozzle of CRH5 EMU Bogie

Direct Strength Method for Steel Deck

EFFECT ON THE AERODYNAMIC FORCING IN AN AXIAL TURBINE STAGE ON VARYING STATOR BLADE SHAPE

CAEFEM v9.5 Information

Elastic Properties of Solids (One or two weights)

Procedure for Performing Stress Analysis by Means of Finite Element Method (FEM)

Sabah Shawkat Cabinet of Structural Engineering Walls carrying vertical loads should be designed as columns. Basically walls are designed in

Bird Impact Analysis of Pre-Stressed Fan Blades Using Explicit Finite Element Code

COUPLED USE OF FEA AND EMA FOR THE INVESTIGATION OF DYNAMIC BEHAVIOUR OF AN INJECTION PUMP

Section 6: PRISMATIC BEAMS. Beam Theory

Transcription:

75-GT-77 I Copyright 1975 by ASME $3. PER COPY author(s). $1. TO ASME MEMBERS The Society shall not be responsible for statements or opinions advanced in papers or in discussion at meetings of the Society or of its Divisions or Sections, or printed in its publications. Discussion is printeo only if the paper is published in an ASME journal or Proceedings. Released for general publication upon presentation. Full credit should be given to ASME, the Technical Division, and the NASTRAN Analysis of a Turbine Blade and Comparison with Test and Field Data J. M. ALLEN Fellow Engineer L. B. ERICKSON Engineer Westinghouse Electric Corp., Gas Turbine Systems Division, Engineering Department, Eddystone, Pa. A NASTRAN finite element analysis of a free standing gas turbine blade is presented. The analysis entails calculation of the first four natural frequencies, mode shapes, and relative vibratory stresses, as well as deflections and stresses due to centrifugal loading. The stiffening effect of the centrifugal force field was accounted for by using NASTRAN's differential stiffness option. Natural frequencies measured in a rotating test correlated well with computed results. Areas of maximum vibratory stress (fundamental mode) coincided with the three zones of crack initiation observed in a metallographic examination of a fatigue failure. Airfoil stress distributions were found to be significantly different from that predicted by generalized beam theory, especially near the airfoil-platform junction. Contributed by the Gas Turbine Division of The American Society of Mechanical Engineers for presentation at the Gas Turbine Conference & Products Show, Houston, Texas, March 2-6, 1975. Manuscript received at ASME Headquarters December 2, 197. Copies will be available until December 1, 1975. THE AMERICAN SOCIETY OF MECHANICAL ENGINEERS, UNITED ENGINEERING CENTER, 35 EAST 7th STREET, NEW YORK, N.Y. 117

NASTRAN Analysis of a Turbine Blade and Comparison with Test and Field Data J. M. ALLEN L. B. ERICKSON INTRODUCTION One of the most significant advances in solid mechanics in recent years is the development of finite element computer programs. One of the most powerful and best known codes is NASTRAN (1, 2), l It is especially well suited for turbine blade analysis because of its capacity to handle dynamic problems with many degrees of freedom, and its capability to include the stiffening effect of centrifugal force (). NASTRAN became a part of the Westinghouse computer system in 197. Since then it has been used extensively in the design and analysis of steam and gas turbine blades. In this paper, results of an analysis of a free standing gas turbine blade are reported, and correlations made using test data and field failure experience. A number of comparisons between NASTRAN and generalized beam theory are made to show some of the limitations in the traditional method. CF STRESSES AND DEFLECTIONS Stresses and deflections were calculated at the design speed of 36 rpm, To account for the stiffening effect of rotation, the differential stiffness feature in NASTRAN was em- FINITE ELEMENT MODEL The subject of this analysis is the last low blade of the Westinghouse W51 family of gas turbines. The airfoil is approximately 18 in. (5.7 cm) high, and has a nearly constant chord of 6.5 in. (16.5 cm), Fig. 1. A proprietary, pre-processer, mesh generator program was used to generate the finite element model. 6 elements, connected at 26 grid (nodal) points, were used to approximate the blade geometry, Fig. 2. The element employed is identified in NASTRAN as CTRIA2 (1) and is a triangular, flat shell element with both inplane (membrane) and bending stiffness. Fig. 2 is a tengential projection of the middle surface of the elements. Fig. 3 is a cross-sectional view of the model near the base of the airfoil. 1 Numbers in parentheses designate References at end of paper. Fig. 1 W51 last row turbine blade 2

X-COORDINATE, CM -8 6 - -2 2 6 8 E AIRFOIL X COORDINATE, INCHES Fig, 3 Element model of airfoil hub cross section } -PLATFORM SHANK Fig. 2 NASTRAN finite-element model mid surface projection ployed. Differential stiffness is a first order approximation to large deformation effects, such as those considered in beam-column action. Briefly, the procedure is to first solve the problem assuming no effect of deflection on loading. Then, using the resulting displacements normal to the plate elements and the membrane forces in each element, a linearized first approximation of the stiffening influence of the membrane forces is obtained in the form of a differential stiffness matrix. It is assumed that the applied loads from which the differential stiffness is derived remain fixed in magnitude and direction during motion of the structure, and that the points of application move with the structure. The elements of the differential stiffness matrix are added to the corresponding elements in the basic stiffness matrix and the problem solved again. This is the only iteration made, though some increase in accuracy could be obtained by additional repetitions of the process. However, it is doubtful that this increase in accuracy is enough to justify the increase in computer costs. All displacements, stresses, and natural frequencies referred to in this paper are based on the differential stiffness method. Boundary conditions are such that the six degrees of freedom are suppressed for the eleven nodal points along the bottom row of the model, with the exception of the translation displacement along this line. Here only the nodal point on the leading edge was fixed from moving in this direction. In other words, the blade is clamped where the blade root and disk steeple first come into contact, but the root is free to move in the direction of the groove from a fixed point on the leading edge, The remaining grid points in the model are each allowed six degrees of freedom. Airfoil displacements due to cf loading are shown in Fig.. The overall distortion may be characterized as a general bending and untwisting of the airfoil. While the untwisting is relatively small in this instance, 1. deg at the tip, this untwisting and/or local distortion of airfoil shape can be an important design consideration in long, thin, airfoils under high cf fields. Cf maximum principal stress patterns are shown in Fig. 5. These are Calcomp plots created by running the NASTRAN stress output through a proprietary, post-processor, plotting program. In these plots, the curved middle surface is plotted flat, which explains the difference in appearance between Figs. 2 and 5. It can be seen that the stress distribution at any given cross section is far from uniform, as would be indicated by elementary beam theory (load divided by crosssectional area). This is especially true at the base of the airfoil, Fig. 6. Because of the untwisting action described previously, the hub section cannot remain plane after deformation, and at the same time, stay compatible with the heavy platform. Thus, a high degree of warping occurs in sections near the airfoil-platform junction, as indicated by the stress pattern exhibited in Fig. 6, The blade was also analyzed using a generalized beam theory (3, 5). For comparison 3

1 9 8 7 6 5 3 MAX MIN 2 1 CONVEX SURFACE SURFACE 1 1 KSI 6.895 MPz Fig, 5 Maximum principal stress (ksi) due to centrifugal force Fig. Deflection due to centrifugal force purposes, the hub section stress results are plotted in Fig. 6. It can be seen that there is no similarity between the NASTRAN and generalized beam patterns. NASTRAN predicts that the maximum stress occurs on the suction side near the leading edge, and is nearly twice the nominal stress (68 versus 37). As we move away from the airfoilplatform junction, the NASTRAN calculated stress distribution changes dramatically and becomes somewhat more uniform. For example, at 2 percent airfoil height, the maximum stress occurs on the concave surface at mid-chord, and is now only 28 percent greater than the nominal, Fig. 7. Here the stress pattern is fairly similar to that predicted by generalized beam theory, Fig. 7. It should be noted that the beam analysis also included the stiffening effect of rotation, and assumed the blade to be clamped at the root. No measurement of strains due to cf loading of this blade are available to compare with these computations. However, it has been observed that a NASTRAN analysis of an axial flow compressor blade yielded cf stress patterns very similar in appearance to those measured by htsuka in his investigation of compressor blade untwist [Reference (5), Figs. 13(b) and 13(d)]. NATURAL FREQUENCIES AND MODE SHAPES The first four natural frequencies were calculated using both NASTRAN and generalized beam theory. The results are presented in Table 1. Also shown are natural frequencies obtained from strain gage output in a rotating test of this blade. Blade excitation was supplied by steam jets in an evacuated chamber. Evacuation was achieved by condensing steam. Correlation between test and NASTRAN results is good. Deviation between average test values and NASTRAN is less than six percent for all modes. Generalized beam theory is less consistent showing deviations of 17 and 27 percent for the second and fourth modes. The main reason NASTRAN predictions are high is because no root-groove flexibility was provided in the model. As pointed out in the discussion on boundary conditions, the blade shank was considered to be clamped along the bottom row of elements. In reality, this condition can never by fully attained. Any cantilever support will elastically distort and allow some rotation at the base. It is important to remember that the same clamping condition was assumed for both beam and NASTRAN models. Mode shapes are displayed in Fig. 8. These are plots of the deformed structure mid-

w z I X COORDINATE, CM 8 8 12 7 e X COORDINATE, INCHES 8 çdnve X NA ST A A N 5 ti y CO NOMINAL 3 N 2 N 2 1 GENERALIZED BEAM THEORY 6 3 NOMINAL - - --- - 1 w CON VE% 2 y 1 ~ C D KE Fig. 6 Stress due to cf hub airfoil section 6 3 X COORDINATE, CM -12-8 - 8 ti 2 3 I w t 3 8 >.1 g 2 3-2 -1 1 2 3 X COORDINATE, INCHES 9 X MIN. MAX. I so RAX, 8 NASTRAN 5 6 NOMINAL --- - - --- 2 y CONVEX D 2 B NOMINAL GENERALIZED BEAM THEORY surface superimposed on the original geometry midsurface. These plots were made using NASTRAN's plotting option (2). Since this is a free vibration solution, only relative displacements can be 3 1 CONVEX 1 Fig. 7 Stress due to of 2 percent height 3 D CONVEX SURFACE SURFACE 1 KXI E835.i1Pa Fig. 9 Relative vibratory stress (ksi) radial direction computed. Hence, one nodal displacement component must be arbitrarily pre-set and all others calculated relative to it. In this case, the largest deflection component in the global coordinate system (the cylindrical coordinate system was used) was set equal to 1 in. (2.5 cm), Fig. 8. The first mode shown in Fig. 8 is, of course, the fundamental mode which is predominently a beam type bending about the minimum inertia axis of the airfoil. The third mode exhibits a high degree of twisting, and has the appearance o.f a 5

Table 1 Calculated and Measured Natural Frequencies MODE AVERAGE TEST FREQUENCY NASTRAN GENERALIZED BEAM CALCULATED % ERROR CALCULATED % ERROR 1 179 Hz 189 Hz +5.6 18 Hz + 2.8 2 391 +3.3 59 +17. 3 535 559 +.5 531 -.75 677 713 +5.3 859 +26.9 T 12 T X COORDINATE, CM 8 8 12 8 X COORDINATE CM 12 8 8 12 2 Q' 2 ^ 1 1 _ O f -1 2-3 -2-1 1 2 3 X COORDINATE, INCHES - 3 2 1 1 2 3 X COORDINATE. INCHES 1-5 6 3 NA ST RAN 1 fi q 5 3 N -5 CONVEX 3 1-6 5 I CONVEX -1 ` 36 fi 15 - GENERALIZED BEAM THETAY 9 15 GENERALIZED BEAM THEORY 9 1 6 1 6 5 3 N x. 5 3 N N y ti 5 3 CONVEX Fig. 1 Mode 1 vibratory stress hub airfoil section 5 5B ^ 3 1 cdnvex Fig. 11 Mode 1 vibratory stress 2 percent height Boo torsional mode. This is probably why the beam program predictions are so good for the first and third modes. A considerable amount of twisting and plate-like bending is present in the second and fourth modes, particularly the fourth mode. This explains why the beam frequency predictions are so poor for these modes. VIBRATORY STRESSES Relative vibratory stresses were also calculated for each of the four modes. Since we 6

Fig. 1 Typical trailing edge crack (1.5 x). Arrow indicates origin ^g. Fig. 13 Typical leading edge crack (1.5 x). Arrow indicates origin only have data on the first mode to compare with the analytical results, only the first mode stress patterns are included. Fig. 9 shows the isostress contours for the stress in the radial direction (the maximum or minimum principal and radial stresses are very nearly the same for this mode). The relative stresses shown are quite large because of the arbitrarily chosen maximum deflection of 1 in. (2.5 cm). Fig. 1 is a plot of surface stress versus axial length of the base airfoil section. Surface stresses computed by the generalized beam program are also plotted on the same graph. Here again, we see that plane sections do not remain plane after deformation. This is especially noticeable on the leading edge where the maximum stress calculated by NASTRAN occurs on the convex side, which is opposite to that which is demanded by beam theory. The other peak stress points are on the convex surface at the thickest part of the airfoil and on the trailing edge con- Fig. 15 Typical convex side crack (leading edge crack on left). Arrow indicates origin cave side. This agrees with the beam prediction, but NASTRAN shows a larger gradient of stress across the trailing edge. This comes about because of local (plate-like) bending. Away from the airfoil-platform junction, the stress pattern changes and becomes more like the beam distribution. -Phis is illustrated in Fig. 11 for the 2 percent height section. Note, however, that some plate-like bending occurs here also. Direct measurements of stress patterns in this blade have never been made, but field failure data does exist to support the NASTRAN vibratory stress results. In 1972, a W51 unit experienced a severe fuel nozzle clogging problem, which resulted in a large third engine order stimulus at design speed (36 rpm). This set up a fundamental F

mode resonant vibration. Two blades fractured and five others developed cracks near the airfoil base. Metallographic examination of the cracked blades revealed that the failures were due to high cycle fatigue, and that there were three points of crack initiation, Figs. 12, 13, 1, and 15. The most interesting point to note is that the leading edge crack, Figs. 12 and 13, did indeed start on the convex side as predicted by NASTRAN; not on the concave side as required by beam theory, Fig. 1. The other cracks started on the concave side of the trailing edge, Fig. 1, and on the convex side at the thickest part of the airfoil, Fig. 15, as predicted by both NASTRAN and beam theory. COMPUTER COSTS The results reported herein were obtained using a CDC 66 version of NASTRAN, level 15. The cost of the analysis, including mesh generation, was about 13, based on a costing rate of.2 per cs second. Since these calculations were made, Westinghouse has installed CDC 76 computers for its engineering and scientific computations. Experience with the 76 version indicates that the cost of this analysis would be about one-fourth as much, or X36O. A factor of one-fourth cannot be used as a general rule to estimate 66 versus 76 costs because it depends on the program, the problem type, and the size of the problem. The important point is that the cost of using NASTRAN as a design tool is reasonable, especially if a CDC 76 version is available to the user. CONCLUSIONS It has been demonstrated that NASTRAN is a powerful tool for the calculation of elastic stresses, displacements, and natural frequencies in turbine blades. Although strains were not measured, field failure experience shows that peak vibratory stress locations predicted by NASTRAN agree with the three crack initiation zones observed in failed blades. Of special interest is the fact that NASTRAN correctly predicts that the fundamental mode peak stress occurs at the airfoil base section on the suction side of the leading edge, rather than on the pressure side, as predicted by beam theory. In view of this correlation, one would expect the of stress and deflection patterns to be equally good. Also, a close similarity was noted between calculated of stress patterns for a Westinghouse axial flow compressor blade, and the patterns measured by Ohtsuka in his investigation of compressor blade untwist. The cf stress distribution at any given cross section was found to be far from uniform, especially at the base of the airfoil. Here the peak stress occurred on the suction side of the leading edge and is nearly twice the nominal stress. The correlation between calculated and measured natural frequencies (first four modes) was good. Deviation between average test frequency and calculation was less than six percent for all modes, whereas generalized beam theory was less consistent, showing 17 and 27 percent deviations for the second and fourth modes. NASTRAN is uniquely suited for turbine blade analysis because of its capacity to handle large dynamic problems, and its differential stiffness feature, which takes into account the approximate effect of deflection onloading. REFERENCES 1 MacNeal, R. H., ed., "NASTRAN Theoretical Manual," NASA SP-221, Sept. 197, 2 McCormick, C. W., ed., "NASTRAN Users' Manual," NASA SP-222, Sept. 197. 3 Houbolt, J. C., and Brooks, G. W., "Differential Equations of Motion for Combined Flapwise Bending, Chordwise Bending, and Torsion of Twisted Non-Uniform Rotor Blades," NACA Report 136, 195 8. Filstrup, A. W., "Finite Element Analysis of a Gas Turbine Blade," ASME Paper 7-- WA/GT-11, Nov. 197. 5 Ohtsuka, M., "Untwist of Rotating Blades," ASME Paper 7-GT-2, Apr. 197. 8