UDK 539.42:669.24:519.61/.64 Original scientific article/Izvirni znanstveni članek ISSN 1580-2949 MTAEC9, 49(5)773(2015) THERMAL FATIGUE OF SINGLE-CRYSTAL SUPERALLOYS: EXPERIMENTS, CRACK-INITIATION AND CRACK-PROPAGATION CRITERIA TOPLOTNO UTRUJANJE MONOKRISTALNIH SUPERZLITIN: PREIZKUSI, MERILA ZA NASTANEK IN NAPREDOVANJE RAZPOKE Leonid Getsov1, Artem Semenov2, Sergey Semenov2, Alexander Rybnikov1, Elena Tikhomirova3 1NPO CKTI, Saint-Petersburg, Russia 2Saint-Petersburg State Polytechnical University, Russia 3Klimov Company, Saint-Petersburg, Russia guetsov@yahoo.com, semenov.artem@googlemail.com Prejem rokopisa - received: 2014-10-03; sprejem za objavo - accepted for publication: 2014-11-06 doi:10.17222/mit.2014.251 Thermal-fatigue tests were performed on corset samples made from five single-crystal nickel-based superalloys with various crystallographic orientations ((001),(011),(111)) at different temperatures and cycle durations. The proposed criteria for the thermal-fatigue crack initiation and crack propagation are discussed and compared with experimental results. The sensitivity of durability to the crystallographic orientation and to thermal-cycle parameters is analyzed experimentally and by means of the corresponding finite-element simulations. Keywords: nickel-based superalloy, single crystal, thermal fatigue, crack initiation, crack propagation Preizkusi toplotnega utrujanja so bili izvršeni na vzorcih korzetne oblike. Vzorci so bili izdelani iz petih monokristalnih nikljevih superzlitin z različno kristalografsko orientacijo ((001),(011),(111)), preizkusi pa so bili narejeni pri različnih temperaturah trajanja ciklov. Predlagana merila za nastanek razpok in njihovo napredovanje pri toplotnem utrujanju so opisana in primerjana z eksperimentalnimi rezultati. Odvisnost zdržljivosti od kristalografske orientacije in parametrov toplotnih ciklov je eksperimentalno analizirana in ugotovljena z ustrezno simulacijo s končnimi elementi. Ključne besede: superzlitina na osnovi niklja, monokristal, toplotno utrujanje, nastanek razpok, napredovanje razpoke 1 EXPERIMENTAL DATA The durability of the single-crystal blades of modern gas turbines is mainly determined by their ability to resist the thermal-cycling stress. In Russia, the testing of the thermal fatigue of corset samples (Figure 1) in vacuum conditions, with and without lateral thermal-bar- Figure 1: Specimen for thermal-fatigue tests Slika 1: Vzorec za preizkus toplotnega utrujanja rier coatings is widespread. During the test, the shear-band formation, the conditions of crack nucleation and the rate of crack propagation are investigated. Experimental studies were performed on five single-crystal alloys (ZhS32, ZhS36 and others) developed by the VIAM research centre, with different alloying and, most importantly, different carbon contents. The chemical compositions of the alloys are given in Table 1. The tests of the thermal-fatigue resistance were carried out for different crystallographic orientations (CGOs) of the specimens, various minimum (Tmin) and maximum temperatures (r^ax) of the cycle, with and without a dwell time at Tmax.1,2 The test results were characterized with the number of cycles to the macrocrack initiation and to the failure of the samples and the in-plane and out-of-plane residual strains (£i„, £out) in the central parts of the specimens. The changing of Tmax and CGO has a significant influence on the thermal-cycling durability (Table 2). The results of a comparative analysis of the thermal-fatigue durabilities of the five considered alloys are given in Figure 2, for a thermal cycle of 150^900 °C (Figures 2a and 2c) and for a cycle of 500^1050 °C (Figures 2b and 2d). Table 1: Chemical compositions of Ni-based single-crystal alloys (w/%) Tabela 1: Kemijska sestava monokristalnih zlitin na osnovi Ni (w/%) Cr Orientation AT/°C T ■ /°C min Tmax/°C Number of cycles to failure [001] 750 150 900 951 200 950 450 250 1000 63 500 450 950 2635 500 1000 1220 550 1050 356 [011] 750 150 900 560 250 1000 95 Re Al Ti Nb Re La Ru Ze Si Alloy 1 0.120.18 4.3-5.6 3.5-4.5 5.6-6.3 1.4-1.! Alloy 0.05 4.5 2.0 5.9 1.0 1.0 2.0 Alloy 3 <0.02 2.85 6.2 6.0 3.9 <0.2 Alloy 4 <0.02 4.2 4.05 6.0 0.8 <0.2 Alloy 5 <0.02 3.2-3.8 4.6-5.4 4.2-4. 2.2-2.i 5.7-6.3 0.250.55 8.1-8.6 0.01-0.2 <0.02 <0.2 Figure 2: Comparison of thermal-fatigue resistance values for alloys 1-5 (number of cycles (a, b) and transversal residual strains at failure (c, d) for cycles of 150«900 °C (a, c) and 500«1050 °C (b, d)). Slika 2: Primerjava odpornosti proti toplotnemu utrujanju zlitin 1-5 (število ciklov (a, b) in prečne preostale napetosti pri porušitvi (c, d) pri ciklih od 150«900 °C (a, c) in 500«1050 °C (b, d)) Table 2: Effects of cycle temperature and CGO (alloy 2) Tabela 2: Vpliv temperature cikliranja in CGO (zlitina 2) Orientation [001] [011] Ar/°c 750 500 750 150 200 250 450 500 550 150 250 rm,x/°C 900 950 1000 950 1000 1050 900 1000 Number of cycles to failure 951 450 63 2635 1220 356 560 95 Table 3: Effects of dwell time at rmax and of aging on the thermal-fatigue resistance for alloy 3 at thermal cycle of 500»1050 °C Tabela 3: Vpliv zdržljivosti pri Tmax in staranja na odpornost proti toplotnemu utrujanju zlitine 3 pri toplotnem cikliranju 500»1050 °C Condition and test mode Without a dwell time at Tm^X Dwell time of 2 min at Tmax Aging at 850 °C, 500 h Number of cycles to macrocrack 670 27 85 Number of cycles to failure 746 91 242 16 12 13 eouf/% -1 17 £in, £out - residual strains in 2 directions transversal to the sample axis The obtained results indicate a relatively weak influence of the chemical composition of the single-crystal alloys on the thermal-fatigue durability (the differences are less than 2-5 times). On the other hand, the long-term strengths and creep-resistance values of the single-crystal alloys with different compositions differ quite significantly.3 None of the alloys among 1-5 demonstrates the maximum durability throughout the whole thermal-cycling-parameter range. A significant decrease in the durability is observed after the dwell time at the maximum temperature of the cycle (Table 3). The preliminary aging at 850 °C also leads to a considerable decrease in the durability and residual strains 2 CRITERION FOR THE TMF-CRACK INITIATION The prediction of the thermo-mechanical-fatigue (TMF) failures of the single-crystal materials is per- 4 formed on the basis of the deformation criterion.4-6 The macrocrack-initiation criterion is the condition for achieving the critical value of the total damage initiated by different mechanisms: D1(A8 pq) + D2(Ae ;)+03(8 pq)+D,(8 ;) = 1 (1) Criterion (1) is based on the linear damage summation of the cyclic plastic strain: ,p ( A8 pq, )k D1 t! ) cyclic creep strain: (A8 D2 =L 1=1 C,(T^) one-sided accumulated plastic strain: 8 pq D = max 8 p( T) (2) (3) (4) (5) and one-sided accumulated creep strain: 8 eq D4 = max —— 4 0 < ^ < ^ max 8 cr( T ) where Ci, C2, k, m, 8 p!, 8 T are the material parameters, depending on the temperature and crystallographic orientation. Usually the material parameters k =2, m = 5/4, Ci = (8rp)k, C2 = (3/48rc)m are used.4 Different norms of the strain tensor are considered as the equivalent strains for (1)-(5), including the maximum shear strain in the slip system with normal nim} to the slip plane and slip direction l(011): 8eq = n{111} ■ e ■ l(011) (6) the maximum principal strain (the maximum eigenvalue of the strain tensor): 8eq = 81 = max arg{det(E->ll) = 0} (7) the von Mises strain intensity: 8eq [(8^ - 8, )2 + (8, - 8 , )2 + (8 , - 8, )2]+ Y+ 7J, + 7^ ) (8) the Hill strain intensity: 8eq [(- 8y )2 + (8y - 8 , )2 + (8 , - )2]+ K1*( YI + YJ, + YL ) (9) Figure 3: Damage-field distribution after the 10th cycle for sample 43 from alloy 3 with near to (001) orientation Slika 3: Razporeditev polja poškodbe po 10. ciklu za vzorec 43 iz zlitine 3 z orientacijo blizu (001) and the maximum shear strain: 1 = 2 ( 81 - 8 3 ) (10) / \ 2 4 8 r ( 001) -1 v 8 r ( 011) ^ Equivalent strains (6) and (9) correspond to the crystallographic failure mode, while equivalent strains (7), (8) and (10) correspond to the non-crystallographic failure mode. In (9): K* = 9 3 RESULTS OF THE FINITE-ELEMENT SIMULATION In order to verify criterion (1) a stress-strain analysis of the single-crystal corset samples (Figure 1) was carried out using the finite-element (FE) program PANTOCRATOR7 with an application of physical models of elastoplasticity. These material models take into account the fact that inelastic deformation occurs in accordance with the crystal slip system due to a slip mechanism and, therefore, the deformation is strongly sensitive to the crystallographic orientation. The elasto-plastic and visco-elasto-plastic material models8,9 with nonlinear kinematic and isotropic hardening, also including the self-hardening of each system and the latent hardening10 between the slip systems, are used in the FE simulations. The application of the visco-elastic models leads to unrealistically overestimated levels of the stress. The obtained results for the inhomogeneous stress, strain and damage fields allowed us to find the location of the crack initiation. The damage field is determined with criterion (1) on the basis of an analysis of the strain-field evolution using the experimental data on creep- and elastoplastic-deformation curves. The typical damage field distribution after the 10th cycle (20 °C^Tmax = 900 °C^Tmin = 150 °C) is presented in Figure 3 for the sample from alloy 3 with the orientation near to (001). The best prediction (in comparison with Figure 4: Influence of the crystal orientation on the cyclic stressstrain curve. Results of finite-element simulations of thermal cycles with Tmin = 150 °C and Tmax = 900 °C (the central point of a specimen). Slika 4: Vpliv orientacije kristala na ciklično krivuljo napetost -raztezek. Rezultati simulacije s končnimi elementi za toplotne cikle s Tmin = 150 °C in Tmax = 900 °C (sredina vzorca). eq -1 Figure 5: Comparison of deformed states (for clarity increased 10 times) of the central sections of the specimens with orientations (001), <011) and (111) for the cooling phase (T = 150 °C) after the 10th cycle. Colors show the damage-field distribution. Slika 5: Primerjava deformacijskih stanj (za boljšo prepoznavnost povečano 10-krat) v srednjem prerezu vzorcev z orientacijo (001), (011) in (111) pri fazi ohlajanja (T = 150 °C) po 10. ciklu. Barve prikazujejo področje razporeditve poškodbe. the experiment) of the number of cycles for the crack initiation is given for strain measure (6) in this loading case. The results of the FE simulations show that the cry-stallographic orientation has a significant influence on the stress-strain state of the samples (Figures 4 and 5) which was also confirmed with the experiments.10 The results of the verification of criterion (1) for alloys 2 and 3 are given in Figure 6. The computed number of the cycles to crack initiation demonstrates a satisfactory accuracy. The thermal-fatigue failure process for the corset specimens shows the presence of two distinct modes (cry- Figure 7: Diagram of thermal-fatigue failure modes (alloy 1) Slika 7: Diagram vrste poškodb pri toplotnem utrujanju (zlitina 1) stallographic and non-crystallographic) of the failure. Their occurrence depends on the CGO and on temperatures Tmin and Tmax of the cycle. The experimental results and the Arrhenius approximation of the boundary between the failure modes are given in Figure 7. The difference in the failure mode is taken into consideration for the damage computation based on equation (1) by choosing the equivalent strain measure (6)-(10) and analysing the crack propagation by choosing the direction of the crack growth. The results of the experiments and computations showed that the dwell time has a significant effect on the thermal-fatigue durability. Figure 8 shows the computed effect of the influence of the dwell time at Tmax on the thermal-fatigue durability for the specimen from alloy 3 and for the thermal cycle of Tmn = 500 °C ^ Tmax = 1050 °C. A comparison of the numbers of the cycles to crack initiation computed with equation (1) and obtained with the experimental results is given in Table 4. Figure 6: Comparison of experimental data with the results of computations of thermal-fatigue durability for alloys 2 and 3 Slika 6: Primerjava eksperimentalnih podatkov z rezultati izračunov odpornosti proti termičnemu utrujanju za zlitini 2 in 3 Figure 8: Influence of the dwell time at Tmax on the thermal-fatigue durability (alloy 3) Slika 8: Vpliv časa zadržanja pri Tmax na zdržljivost pri toplotnem utrujanju (zlitina 3) Table 4: Comparison of the numbers of the cycles to crack initiation obtained with the computation and with the experiment for various dwell times at rmax (alloy 3) Tabela 4: Primerjava izra~unanega {tevila ciklov do nastanka razpoke in eksperimentalnih rezultatov pri različnih časih zadržanja na T max (zlitina 3) Computation with (1) Experiment Dwell time 0 s 480 670 120 s 45 27 Alloys 1-5 have significantly different creep characteristics and creep-rupture strengths, therefore it can be expected that the thermal-fatigue durability in the case of a dwell time at Tmax varies significantly for different alloys, contrary to the cases of no dwell time at Tmax. 4 CRITERION FOR THE TMF-CRACK PROPAGATION The direction of a crack propagation in a single crystal is defined by the crystallographic structure (typically it is plane (111) for the considered alloys1112 in the crystallographic mode of failure) and, correspondingly, it is defined by the stress state during the non-crystallogra-phic failure mode (a high temperature). The criterion for the thermal-fatigue-crack propaga-tion,13 based on a linear sum of the contributions of the fatigue and creep, is used in order to determine the rate of crack propagation and the number of cycles to failure. The criterion is based on using stress-intensity factor A^eff (for the description of the fatigue) and the C*-inte-gral (for the description of the creep during the dwell time at the maximum temperature of the cycle). The mathematical formulation of the criterion represents a generalization of the Paris power-type equation by taking into account two mechanisms of cracking:13 dn dN = B()m + A(C *(r))qdr (11) where A, B, m, q are the material constants, defined separately from the fatigue experiments as dn/dN(A^eff) and creep experiments as dn/dt(C*) crack-growth data. The crack-propagation process in the corset sample (Figure 1) from alloy 2 is simulated cycle by cycle using Figure 9: Views of the specimen from alloy 2 after the failure: a) left half, b) right half Slika 9: Videz vzorca zlitine 2 po porušitvi: a) leva polovica, b) desna polovica Figure 10: Plastic-strain-intensity-field distribution in the specimen from alloy 2 with a crack of a length of 0.5 mm Slika 10: Razporeditev polja intenzitete plastične deformacije v vzorcu iz zlitine 2 z dolžino razpoke 0,5 mm the FE program ABAQUS.14 The specimen has its CGO near to (001). The deviation from the axial orientation is 5.47 ° and from the azimuthal orientation it is 41.97 ° (the Euler angles for rotations ZX'Z' are: (p = 354.7 9 = 41.7 ^ = 89.8 °). The CGO was determined with Laue diffraction patterns. The location of the crack initiation is defined with an FE computation in section 3 and it corresponds with the experimentally obtained result. The crack propagates in plane (111) (Figure 9), except in the rupture stage. The crack-propagation process is simulated only in plane (111), without out-of-plane deviations, using a step-by-step technique with FE remeshing for the environment of the crack front at every step. The initial crack front has an elliptic form. The AJen integral is used to calculate A^eff. The fracture zone is established by comparing the calculated J integral with the Jic value. The details of the simulation are described in.15 In the experiment the sample failed after 560 thermal cycles (Tmin = 150 °C Tmax = 900 °C); the registered crack-initiation stage lasted for 435 cycles and the crack-propagation stage took 125 cycles. A comparison of the computed number of cycles for the crack-propagation stage based on the FE simulation (Figure 10) with the experimental results is given in Table 5. A simulated crack-front evolution is shown in Figure 11a. It demonstrates a good correlation with the data of a fractography analysis (Figure 11b). Table 5: Comparison of the estimated number of cycles for the crack-propagation stage with the experimental results (alloy 2) Tabela 5: Primerjava določenega števila ciklov za fazo napredovanja razpok z eksperimentalnimi rezultati (zlitina 2) FE computation with (11) Experiment Number of cycles 199 125 Figure 11: Evolution of the thermal-fatigue crack-front propagation on a fracture surface: a) results of finite-element simulations, b) data of a fractography analysis (specimen from alloy 2) Slika 11: Razvoj napredovanja fronte toplotne utrujenostne razpoke na površini preloma: a) rezultati simulacije s končnimi elementi, b) podatki iz analize preloma (vzorec iz zlitine 2) 5 CONCLUSIONS 1. Thermal-fatigue tests were performed on corset samples from five Russian single-crystal nickel-based superalloys with various crystallographic orientations «001), <011), <111» under different temperatures and cycle durations with the aims to determine the characteristics of the thermal-fatigue resistance, compare the alloys and determine the thermal-fatigue-failure criterion. 2. The crack-initiation criterion for the single-crystal alloys under thermal cyclic loadings is proposed and discussed. A satisfactory correlation is observed between the proposed deformation criterion (1) and the obtained experimental results. 3. The finite-element simulations of the single-crystal corset specimens under the thermal cyclic loading were performed using different material models. The obtained results indicate an applicability of the proposed deformation criteria for failure for the five considered single-crystal alloys at temperatures of up to 1100 °C. 4. A high sensitivity of the thermal-fatigue durability to the crystallographic orientation and to the thermal-cycle parameters is observed experimentally and also in the corresponding finite-element simulations of the corset specimens. 5. The thermal-fatigue failures of the corset specimens show the presence of crystallographic and non-crys-tallographic modes of failure. The occurrence of these modes depends on the crystallographic orientation and the temperature parameters of the cycle. 6. The proposed crack-propagation criterion (11) under the thermal cyclic loadings demonstrates a satisfactory accuracy in comparison with the experimentally obtained number of cycles to failure and the crack-front evolution indicated by the data of a fracto-graphy analysis. 7. To obtain reliable analysis results for the thermal-fatigue strength and durability of the single-crystal blades of GTE, further detailed studies of the material characteristics, including the crack-resistance parameters, carried out over a wide temperature range, are required. Acknowledgements The present study is supported by Russian Science Foundation, grant No.15-19-00091. 6 REFERENCES ^L. B. Getsov, N. I. Dobina, A. I. Rybnikov, A. S. Semenov, A. A. Staroselsky, N. V. Tumanov, Strength of Materials, 40 (2008) 5, 538-551, doi:10.1007/s11223-008-9076-1 2 L. B. Getsov, Materials and strength of gas turbine components 1-2, Publ. House, Rybinsk 2011 3 N. V. Petrushin, I. L. Svetlov, O. G. Ospenkova, Heat-resistant casting nickel alloys, All materials, Encyclopedic Handbook 6, 2012, 15-19 4 A. S. Semenov, L. B. Getsov, Strength of Materials, 46 (2014) 1, 38-48, doi:10.1007/s11223-014-9513-2 5L. B. Getsov, A. S. Semenov, E. A. Tikhomirova, A. I. Rybnikov, Mater. Tehnol., 48 (2014) 2, 101-106 6L. Getsov, A. Semenov, A. Staroselsky, Mater. Tehnol., 42 (2008) 1, 3-12 7 A. S. Semenov, Proc. of V Int. Conf. Scientific and technical problems of forecasting the reliability and durability of the structures and methods for their solution, St. Petersburg, 2003, 466-480 8G. Cailletaud, International Journal of Plasticity, 8 (1992) 1, 55-73, doi:10.1016/0749-6419(92)90038-E 9 J. Besson, G. Cailletaud, J. L. Chaboche, S. Forest, Non-linear mechanics of materials, 1s' ed., Springer, Netherlands 2010, doi:10.1007/978-90-481-3356-7 10 U. F. Kocks, T. J. Brown, Acta Metall., 14 (1966) 2, 87-98, doi:10.1016/0001-6160(66)90290-2 11 P. A. S. Reed, I. Sinclair, X. D. Wu, Metallurgical and Materials Transactions A, 31 (2000) 1, 109-123, doi:10.1007/s11661-000-0058-6 12 J. Telesman, L. J. Ghosn, Accelerated fatigue crack growth behavior of PWA 1480 single crystal alloy and its dependence on the deformation mode, In: S. Reichman, D. N. Duhl, G. Maurer, S. Antolovich, C. Lund (Eds.), Superalloys 1988, The Metallurgical Society, 1988, 615-624 13 A. S. Semenov, S. G. Semenov, A. A. Nasarenko, L. B. Getsov, Mater. Tehnol., 46 (2012) 3, 197-203 14 ABAQUS Analysis User's Manual, Version 6.10 15 S. G. Semenov, A. S. Semenov, L. B. Getsov, B. E. Melnikov, Proc. Int. XLI Summer School - Conference Advanced Problems in Mechanics, St. Petersburg, 2013, 74-81