Acta Polytechnica CTU Proceedings doi:10.14311/APP.2018.20.0056 Acta Polytechnica CTU Proceedings 20:56–64, 2018 © Czech Technical University in Prague, 2018 available online at http://ojs.cvut.cz/ojs/index.php/app THERMO-MECHANICAL FATIGUE ANALYSIS OF A STEAM TURBINE SHAFT Martin Nesládeka, ∗, Josef Jurenkaa, Michal Bartošáka, Milan Růžičkaa, Maxim Lutovinova, Jan Papugaa, Radek Procházkab, Jan Džuganb, Petr Měšťánekc a Czech Technical University in Prague, Faculty of Mechanical Engineering, Technická 4, 166 07 Prague 6, Czech Republic b COMTES FHT a.s., Průmyslová 995, 334 41 Dobřany, Czech Republic c Doosan Škoda Power s.r.o., Tylova 1/57, 301 28 Plzeň, Czech Republic ∗ corresponding author: martin.nesladek@fs.cvut.cz Abstract. Increasing demands on the flexibility of steam turbines due to the use of renewable energy sources substantially alters the fatigue strength requirements of components of these devices. Rapid start-ups as well as the increased number of the load cycles applied to the turbines must be handled by design methodologies. The goal of the work presented in this paper was to provide a computational framework applicable to the thermo-mechanical fatigue (TMF) prediction of steam turbine shafts. The so-called Damage Operator Approach by Nagode et al. has been implemented to the software codes and applied to fatigue analysis of the thermo-mechanical material response computed numerically by the finite element analysis. Experimental program conducted in order to identify the material thermo-mechanical behavior and to verify numerical simulations is introduced in the paper. Some results of TMF prediction of a sample steam turbine shaft are shown. Keywords: Fatigue prediction, thermo-mechanical fatigue, finite element analysis, Chaboche elastic- plastic model. 1. Introduction The energy industry has been experiencing signifi- cant changes recently. The overall share of renewable sources in the European Union countries has reached 17 percent [1]. The leading country in this statistic is Sweden having more than 50 percent of energy from renewable sources. In Czech Republic it is approxi- mately 15 percent and the overall trend is that this share will continue to grow significantly in the future - the EU strategy is to reach at least 27% by 2030. However, instability of the power output due to the nature of the renewable sources may cause problems in energy supplies if it is not compensated by some back-up sources. Currently, there are no reliable ways to save energy for power shortages in the network in a sufficiently large scale. This role must be taken over by traditional thermal power sources such as gas-fired power plants providing "just-in-time" energy supplies. The problem is that these plants already have to or will have to cover shortages in the daily regime. Increased number of start-ups and shut-downs as long as increasing demands on speed-up of the start-ups set new challenges in the gas and steam turbine design. Components of steam turbines are often exposed to thermo-mechanical loading due to variable thermal and mechanical loading conditions induced by the turbine start-ups and shut-downs. The main damage mechanisms, which if they occur together constitute the therm "thermo-mechanical faigue" (TMF), are often referred to as oxidation, mechanical fatigue and creep [2, 3]. Constitutive material model derived from the material testing, finite element analysis and suitable fatigue criterion is required for the reliable damage assessment of components subjected to cyclic thermo-mechanical loading conditions. This work is focused on the TMF assessment of a sample case of a steam turbine shaft. Motivation for this work originates from the need for revising the current steam turbine design strategies due to the above-mentioned reasons. Classical low-cycle fatigue (LCF) criteria as Smith- Watson-Topper (SWT) [4] and the Nihei [5] methods in their original form can only cope with isothermal loading. For the case of thermo-mechanical loading, cycle reference temperature, which could be used to select fatigue curves for the isothermal models, is somewhat arbitrary and error-prone. The so-called Damage Operator Approach (DOA) by Nagode et al. [6, 7] overcomes this issue by intro- ducing a procedure for continuous damage parameter calculation. In the original form it is based on the SWT method, however the damage induced by non- stationary thermal and mechanical loading is evalu- ated in a chronological manner, i.e. damage is com- puted and accumulated in any time point no matter if it is load reversal or not. Such approach allows for handling arbitrary thermal and mechanical load 56 http://dx.doi.org/10.14311/APP.2018.20.0056 http://ojs.cvut.cz/ojs/index.php/app vol. 20/2018 TMF analysis of a steam turbine shaft sequences without the need for evaluating the above- mentioned cycle reference temperature. The work presented here was done within the Flex- turbine project, which aims to revise the design of key turbine components and propose new design methodologies so that the turbines can cope with the changing situation in the transmission grids. The paper describes the approach to prediction of thermo- mechanical fatigue of steam turbine shafts. Especially the high-pressure sections of turbine shafts are ex- posed to complex thermo-mechanical loading due to variable service conditions. Within the experimental part of the work, experi- ments for identifying parameters of the selected consti- tutive model as well as parameters of damage model were carried out. The assumed constitutive model of the material mechanical response is temperature- dependent elastic-plastic and for the purpose of the structural finite element analysis (FEA) parameters of the Chaboche non-linear kinematic hardening were identified. The DOA allows for processing stress histo- ries computed by FEA in order to perform correction on the material viscous behavior, if needed. Strain controlled fatigue tests with hold period may be used for finding the parameters of the Maxwell model of vis- coplasticity, which may then be applied for evaluating the viscous strain and subsequent stress relaxation. A set of TMF tests for verifying the prediction methodology have been scheduled. Various phasing of temperature and mechanical load is tested. The experiments are conceived as strain-controlled with two channels of strain applied to tubular specimens - axial and torsional. TMF test setup is briefly intro- duced in the paper and some preliminary results are presented. Thermo-mechanical stress-strain response from the tests and numerical simulations is compared. The elastic-plastic material model is used in the structural FEA of the steam turbine shaft, which is performed on the basis of the simulated tempera- ture fields that are obtained from the previous tran- sient thermal FEA. Finite element calculations for the steam turbine shaft are performed in the com- mercial finite element software. Subsequently, TMF is predicted for the turbine shaft on the basis of the calculated stress, strain and temperature histories. 2. Experimental work As mentioned within the introduction section, the ex- perimental work may be divided into two main groups depending on the general purpose of the individual test - on the one hand basic material tests to identify parameters of constitutive models and on the other hand tests intended for verifying the proposed model- ing methodology. The material subjected to testing was a chrome-molybdenum creep resistant steel appli- cable to steam turbine shafts. Figure 1. The temperature dependent Young’s mod- ulus. Figure 2. The temperature dependent cyclically stable cyclic stress-strain curves. Figure 3. The temperature dependent Manson-Coffin and Basquin curves. 57 M. Nesládek, J. Jurenka, M. Bartošák et al. Acta Polytechnica CTU Proceedings (a) . Detail of grips with the test specimen and in- strumentation (b) . Overall view of the TMF test setup (c) . Tubular test specimen with surface speckle pattern and welded thermocouple Figure 4. TMF test setup and specimen. 1 - servohydraulic axial-torsional test machine, 2 - resistance heating, 3 - camera of the DIC system, 4 - infrared camera, 5 - test specimen. 2.1. Basic mechanical material tests Static tensile tests as well as LCF tests under the isothermal conditions have been carried out. The as- sumed temperature range of the tests was 20 – 600 ◦C, which correspond to operating conditions of the sam- ple steam turbine. Circular specimens with 10 mm diameter in the gauge section have been tested. Ra- diant heating of specimens was controlled by a split furnace. As a minimum, ten samples have been tested on each curve. Fig. 1, 2 and 3 show the main results of the isother- mal material tests, i.e. temperature dependent nor- malized Young’s modulus, cyclic stress-strain curves (CSSCs) and the Manson-Coffin and Basquin curves, respectively. The Young’s modulus and the CSSCs have been used for identifying the elastic-plastic ma- terial model based on the Chaboche’s procedure - section 3. The temperature dependent Manson-Coffin and Basquin curves are used as an input to the DOA. Numerical procedure based on this approach is de- scribed in section 4. LCF strain-controlled tests with hold period have been performed in order to identify temperature de- 0 500 1000 1500 0.2 0.4 0.6 0.8 1 1.2 t [s] N o rm a li z e d s tr e s s [ − ] 450 ° C 550 ° C 600 ° C Figure 5. The temperature dependent stress relax- ation curves. pendent stress relaxation curves - Fig. 5. This is an important input to the DOA allowing for separating elastic-plastic and viscous strain and modeling stress relaxation - section 4. Note that the curves in Fig. 5 are shown for illustrative purposes and will be fur- ther refined with continuation of the experimental program. 58 vol. 20/2018 TMF analysis of a steam turbine shaft Figure 6. Time dependence of axial strain and tem- perature during the in-phase TMF test. Furthermore, creep tests should be conducted in order to identify constants of the Larson-Miller pa- rameter that is applied to evaluating creep damage. 2.2. Verification tests Experimental setup for TMF tests is shown in Fig. 4a and 4b. The test machine is equipped with a servohydraulic actuator allowing for combined axial and torsional loading. Resistance heating is applied to specimens and the induced temperature is con- trolled by a thermocouple welded in the center of the specimen gauge section. Online axial and torsional deformation control is done by an optical measure- ment system based on the digital image correlation (DIC). Tubular specimens with 8.5 mm in diameter and 0.75 mm wall thickness of the gauge section are sprayed with temperature-stable paint creating a ran- dom speckle pattern - Fig. 4c. Design of the test specimens conforms to ASTM E2207 standard. Triangular strain and temperature in-phase and out-of-phase loading is assumed within the tests - Fig. 6 and 7. Until now, tests with axial strain loading at εa = 0.007 and R = −1 have been conducted. The maximum test temperature is 600 ◦C. More tests combining axial and torsional loading with various phases are scheduled in order to check the prediction quality of the selected model. 3. Numerical simulations 3.1. Constitutive material model A constitutive material model must be able to de- scribe the cyclic material behavior for the service temperature range of the steam turbine shaft. The temperature-dependent Chaboche non-linear kine- matic hardening model [8, 9] is selected for the elastic- plastic FEA. The yield criterion may be expressed as f(σ−α) = k, (1) where k is the material temperature dependent yield stress, σ is the stress tensor and α denotes the back- Figure 7. Time dependence of axial strain and tem- perature during the out-of-phase TMF test. stress tensor. The function f() denotes the von Mises stress norm, which may be in this case evaluated as f(σ−α) = √ 3 2 (S −α′) : (S −α′), (2) where S and α′ are the deviatoric parts of the stress and the backstress tensors, respectively. The : oper- ator denotes the tensor contraction. Assuming the non-linear kinematic hardening model, the total back- stress is composed of multiple backstress components as follows: α = ∑ i αi (3) The non-isothermal evolution of the backstress com- ponents is defined as: α̇i = 2 3 Ciε̇p −γiαiλ + 1 Ci ∂Ci ∂T αiṪ. (4) Ci and γi are temperature dependent material param- eters, which may be derived from the material me- chanical cyclic stress-strain response. εp is the plastic strain tensor and λ is accumulated plastic strain. In our case, the model with three backstress terms was selected, where γ3 is assumed constant and equal to zero. Plots in Fig. 8 and 9 show comparison of the mea- sured temperature-dependent material mechanical re- sponse with the one simulated by finite elements using the above-mentioned constitutive model. The plots show good agreement of simulated and experimentally determined material response. 3.2. FEA of a steam turbine shaft The non-linear transient elastic-plastic FEA of a sam- ple steam turbine shaft has been carried out in order to investigate local thermo-mechanical material response due to the turbine load regimes. In this paper, results of the cold-start load scenario are presented. At first, transient thermal analysis has to be conducted in or- der to analyze the time-dependent temperature spatial 59 M. Nesládek, J. Jurenka, M. Bartošák et al. Acta Polytechnica CTU Proceedings Figure 8. Comparison of measured and simulated me- chanical response due to in-phase uniaxial mechanical and thermal loading. Figure 9. Comparison of measured and simulated mechanical response due to out-of-phase uniaxial me- chanical and thermal loading. distribution. The cold-start regime is characterized by the turbine start-up from ambient temperature and after shut-down it cools down back to almost 20 ◦C. So the service temperature range corresponds to ∼ 20 − 600 ◦C and especially during the start-up process the material exhibits severe temperature gra- dients in the subsurface layers since it is exposed to superheated steam. The heating rates are the matter of the individual turbine design that has to ensure that the material fatigue strength corresponding to the expected turbine life is not exceeded. Three basic turbine load regimes are assumed in the design - cold, warm and hot starts differing, from a material point of view, in the initial temperature and allowable heating rates. However, regarding the necessary computation strategy, these regimes are similar. Temperature distribution in the turbine is a func- tion of complex time evolution of heat flux governed by the material surface and steam flow difference and heat transfer coefficients. The induced thermal gradi- ents result in the constrained material expansion with pronounced thermal stresses. Furthermore, mechani- cal loads due to shaft rotation and steam flow pressure on blade airfoils have to be taken into account. Fig. 10 shows time and subsurface distribution of thermo-mechanical material response within a single cold-start cycle in the domain marked in Fig. 10a. The observable discontinuities in temperature are present due to the step changes in the boundary conditions, which partly reflect the nature of the service condi- tions. On the other hand, this also partly results from modeling imperfections caused by higher anal- ysis incrementation than would be suitable to cover some steep transients. However, this is a reasonable compromise between accuracy and solution cost. The plots in Fig. 10 illustrate how complex the material thermo-mechanical response could be in such structure. For instance, in the R3 notch root the cyclic plastic stress-strain state under variable temperature is present. This justifies selection of the DOA as a method to assessing the TMF damage. 4. Thermo-mechanical fatigue prediction Fig. 11 graphically presents the main steps, inputs and outputs of the DOA, which was adopted and implemented to software codes for the purpose of this work. As mentioned in the introduction section, the main advantage of this approach is that it treats thermo-mechanical material response in step-by-step manner (i.e. in increment-by-increment manner if translated into the finite element terminology). This allows for handling arbitrary stress-strain and temper- ature histories in TMF prediction, while assigning the proper parameters of the material damage model to the individual load time points. In general, it requires local elastic-plastic stress and temperature from FEA as an input. The material inputs are temperature-dependent CSSCs, Manson- Coffin and Basquin curves, Larson-Miller parameter and stress relaxation curves for cases, where the ma- terial viscous effect is substantial. Note that d = 1/N denotes mechanical fatigue damage per single cycle. Prandtl hysteresis operator is employed in order to model the continuous fatigue mechanical damage in- fluenced by temperature changes. The total damage is evaluated as summation of the mechanical part and the creep part computed separately. The oxidation effect is treated indirectly since the material fatigue and creep tests are performed under environmental conditions. Palmgren-Miner’s rule [10] and Robnson’s life-fraction rule [11] are used for accumulating the mechanical fatigue and creep damage, respectively. This procedure has been coded to an in-house C++ procedure allowing for processing thermo-mechanical response from FEA. Refer to [12] for more details. 5. Results and discussion TMF prediction results obtained by the DOA are shown in Fig. 12 and Tab. 1. The contour plot in Fig. 12b shows distribution of 6 √ D quantity, where D denotes the TMF damage due to the cold-start turbine load regime. It may be observed that the critical place is the R3 domain. Note that these re- sults were obtained by the DOA without taking the 60 vol. 20/2018 TMF analysis of a steam turbine shaft (a) . The R3 domain and direction of the nodal path from which the results below have been extracted. (b) . Time and spatial distribution of thermo-mechanical material response in the R3 domain. Figure 10. FEA results in the R3 domain of the sample steam turbine shaft. R2 R3 R4 R5 HP1 upper fillet HP1 lower fillet TMF damage 1.31 · 10−4 2.66 · 10−4 1.69 · 10−4 2.11 · 10−4 2.25 · 10−5 1.81 · 10−4 TMF life 7633 3759 5917 4739 44 444 5525 Table 1. TMF fatigue prediction results by the DOA in the selected domains of the steam turbine. 61 M. Nesládek, J. Jurenka, M. Bartošák et al. Acta Polytechnica CTU Proceedings Strain controlled nonlinear Maxwell model to compute the viscoplastic strain εvp Thermo-mechanical elastic-plastic FEA by ANSYS Mechanical Strain response modelling by using the Prandtl hysteresis operators. Assumes kinematic strain hardening. T1 T2 T3 Cyclic stress-strain curves σa εa Tinterp ˙εvp = ( σ − k(T) K(T) )N(T ) T1 T2 T3 Stress relaxation curves σ t Tinterp Elastic-plastic strain calculation εep(t) = ε(t) − εvp(t) for T(t) exceeding the material creep temperature Tc ε(t) Relaxed stress calculation using the Prandtl hysteresis operators. εep(t) T1 T2 T3 Cyclic stress-strain curves σa εa Tinterp T1 T2 T3 Tinterp Low-cycle fatigue curvesPSW T d = 1/N Fatigue damage estimation (low-cycle fatigue Damage Operator Approach) Df Creep damage estimation by using the Larson-Miller parameter Dc σ(t), T(t) εep(t) σ(t) Df Dc D = Dc + Df Total damage Creep master curve σ LMP 1. 2. 3. 4. 5. 6. 7. material inputs DOA processes inputs from FEA Figure 11. Flowchart to illustrate the basics of the Damage Operator Approach by Nagode et al. 62 vol. 20/2018 TMF analysis of a steam turbine shaft (a) . Steam turbine domain nomenclature adopted in this paper. (b) . Sample contour plot of TMF damage distribution analyzed by the DOA in the main sealing part of the steam turbine. Figure 12. TMF prediction by the DOA. material viscoplasticity into account. The difference between these values and those presented in [12] is due to different approach to handling the stress-strain histories. The results presented in this paper have been evaluated by using only the signed von Mises stress from the FEA as an input to the DOA. The mechanical strain values have been obtained by using the Prandtl hysteresis operators, which are part of the DOA - Fig. 11. This approach gives a stress-strain response that more closely matches the physical be- havior of the material than the approach that had been used in [12]. In that previous paper, the signing was applied to both the stress and the strain histo- ries extracted from FEA. This naturally may lead to distorted hysteresis loops and may have resulted to, for instance, different stress means and amplitudes. If these stress and strain histories are used as a direct input to the function evaluating the TMF damage, different results may be obtained as well. 6. Conclusions This paper briefly presents the experimental program conducted in order to provide the necessary material and verification data for the adopted TMF approach. The TMF experimental setup was introduced and the experimentally determined stress-strain response was compared to the simulated one. The comparison shows good agreement of both experimental and numerical data. The DOA in the form allowing for the viscous correction of strain histories and stress relaxation has been implemented to the software codes and the results of the TMF prediction on the sample steam turbine shaft are presented in the text. Acknowledgements This project has received funding from the European Union’s Horizon 2020 research and innovation program under the Grant Agreement No. 653941. References [1] Share of renewables in energy consumption in the EU reached 17% in 2016. http://ec.europa.eu/eurostat/ documents/2995521/8612324/8-25012018-AP-EN.pdf/ 9d28caef-1961-4dd1-a901-af18f121fb2d. Accessed: 2018-07-30. [2] R. W. Neu, H. Sehitoglu. Thermomechanical Fatigue, Oxidation and Creep: Part I. Damage Mechanisms. Metall Trans A 20(9A):1755–1767, 1989. [3] R. W. Neu, H. Sehitoglu. Thermomechanical Fatigue, Oxidation and Creep: Part II. Life Prediction. Metall Trans A 20(9A):1769–1783, 1989. [4] K. N. Smith, P. Watson, T. H. Topper. A Stress-Strain Function for the Fatigue of Metals. J Mater, JMLSA 5(4):767–778, 1970. [5] M. Nihei, P. Heuler, C. Boller, T. Seeger. Evaluation of Mean Stress Effect on Fatigue Life by Use of Damage Parameters. Int J Fatigue 8(3):119–126, 1986. 63 http://ec.europa.eu/eurostat/documents/2995521/8612324/8-25012018-AP-EN.pdf/9d28caef-1961-4dd1-a901-af18f121fb2d http://ec.europa.eu/eurostat/documents/2995521/8612324/8-25012018-AP-EN.pdf/9d28caef-1961-4dd1-a901-af18f121fb2d http://ec.europa.eu/eurostat/documents/2995521/8612324/8-25012018-AP-EN.pdf/9d28caef-1961-4dd1-a901-af18f121fb2d M. Nesládek, J. Jurenka, M. Bartošák et al. Acta Polytechnica CTU Proceedings [6] M. Nagode, F. L angler, M. Hack. A Time-Dependent Damage Operator Approach to Thermo-mechanical Fatigue of Ni-resist D-5S. Int J Fatigue 33(5):692–699, 2011. doi:https://doi.org/10.1016/j.ijfatigue.2010.11.009. [7] M. Nagode. Continuous Damage Parameter Calculation under Thermomechanical Random Loading. MethodsX 1:81–89, 2014. https://doi.org/10.1016/j.mex.2014.07.004. [8] J. L. Chaboche. Time-Independent Constitutive Theories for Cyclic Plasticity. Int J Plasticity 2(2):149–188, 1986. [9] J. L. Chaboche. Constitutive Equations for Cyclic Plasticity and Cyclic Viscoplasticity. Int J Plasticity 5(3):247–302, 1989. [10] A. Palmgren. Die Lebensdauer von Kugellagern. Zeitschrift des Vereines Deutscher Ingenieure 68(14):339–341, 1924. [11] E. L. Robinson. Effect of Temperature Variation on the Long-Time Rupture Strength of Steels. Trans ASME 74:777–780, 1952. [12] M. Nesládek, M. Bartošák, J. Jurenka, et al. Thermo-Mechanical Fatigue Prediction of a Steam Turbine Shaft. MATEC Web Conf 165, 2018. https://doi.org/10.1051/matecconf/201816522016. 64 http://dx.doi.org/https://doi.org/10.1016/j.ijfatigue.2010.11.009 https://doi.org/10.1016/j.mex.2014.07.004 https://doi.org/10.1051/matecconf/201816522016 Acta Polytechnica CTU Proceedings 20:56–64, 2018 1 Introduction 2 Experimental work 2.1 Basic mechanical material tests 2.2 Verification tests 3 Numerical simulations 3.1 Constitutive material model 3.2 FEA of a steam turbine shaft 4 Thermo-mechanical fatigue prediction 5 Results and discussion 6 Conclusions Acknowledgements References