Analysis of crack initiation and propagation in Thermal Barrier Coatings using SEM-Based geometrical model with extended finite element method

Thermal Barrier Coatings (TBC) are widely used to protect the metallic components that operate at harsh conditions of elevated temperatures and oxidizing environments. Thermally grown oxide (TGO) causes cracks formation in the top coat (TC) that may lead to spallation failure of TBC. This work investigates effect of pores and TGO thickness on crack initiation and propagation due to thermal mismatch between TBC layers. Image processing is used to convert an SEM image, including pores, into a finite element (FE) model. An FE model using XFEM implemented in ABAQUS was developed to investigate crack initiation and propagation for various TGO thicknesses considering the effect of plastic deformation of BC, TGO and substrate. Results show that presence of pores changes the critical sites for crack initiation from the TC/TGO interface to be around the pores within the TC. Crack initiation temperatures and crack lengths were found to be affected with both TGO thickness and pores.     

Effect of material properties on residual stress distribution in thermal barrier coatings

Residual stress has a significant influence on the crack nucleation and propagation in thermal barrier coatings (TBC) system. In this work, the residual stress in the air plasma spraying (APS) TBC system during cooling process was numerically studied, and the influence of the material properties of each layer on the residual stress was investigated. The morphologies of the interface were described by a piecewise cosine function, and the amplitude for each segment gradually increases. The elasticity, plasticity and creep of top coat (TC), thermally grown oxide (TGO) layer and bond coat (BC) were considered and the elasticity and creep of the substrate layer were taken into account. The material properties of all layers vary with temperature. The results show that the material properties have complex influence on the residual stress during cooling. The effect of the material properties of TC and BC on the residual stress at the interface is relatively large, and that of TGO and substrate is relatively small. These results provide important insight into the failure mechanism of air plasma spraying thermal barrier coatings, and important guidance for the optimization of thermal barrier coating interfaces.     

Comprehensive dynamic failure mechanism of thermal barrier coatings based on a novel crack propagation and TGO growth coupling model

Comprehensive understanding of failure mechanism of thermal barrier coatings (TBCs) is essential to develop the next generation advanced TBCs with longer lifetime. In this study, a novel numerical model coupling crack propagation and thermally grown oxide (TGO) growth is developed. The residual stresses induced in the top coat (TC) and in the TGO are calculated during thermal cycling. The stresses in the TC are used to calculate strain energy release rates (SERRs) for in-plane cracking above the valley of undulation. The overall dynamic failure process, including successive crack propagation, coalescence and spalling, is examined using extended finite element method (XFEM). The results show that the tensile stress in the TC increases continuously with an increase in an undulation amplitude. The SERRs for TC cracks accumulate with cycling, resulting in the propagation of crack toward the TC/TGO interface. The TGO cracks nucleate at the peak of the TGO/bond coat (BC) interface and propagate toward the flank region of the TC/TGO interface. Both TC cracks and TGO cracks successively propagate and finally linkup leading to coating spallation. The propagation and coalescence behavior of cracks predicted by this model are in accordance with the experiment observations. Therefore, this study proposed coating optimization methods towards advanced TBCs with prolonged thermal cyclic lifetime.     

Influences of interface morphology and thermally grown oxide thickness on residual stress distribution in thermal barrier coating system

Calculation of residual stress with finite element method is a basic work in failure mechanism investigation in thermal barrier coating (TBC) system because the residual stress is main driving force for crack nucleation and propagation. In this work, a complicated cosine curve with gradually increasing amplitude was used to simulate interface morphologies between layers so as to study the residual stress behavior during the cooling process in air plasma spraying TBC system by finite element method. The substrate, thermally grown oxide (TGO) and top coat (TC) are considered to be elastic and bond coat (BC) elastic-perfectly plastic. The material properties are all temperature dependent. The stress result comparison between models with and without substrate shows the effect of substrate on the residual stress distribution around layers interfaces should not be ignored as the substrate influences the value of normal residual stress as well as the stress distribution along undulating interfaces. Then the model with substrate was used to study the residual stress evolution along interfaces during cooling down from the temperature of 1000 °C to room temperature. The influences of the thickness of TGO and the amplitude and wavelength of interface on the residual stress distributions near interfaces were considered. The results show that these influences are very complicated. Meanwhile, it's found that the hybrid roughness parameter containing information for height and spacing is more suitable to describe the interface complicacy. The results facilitate understanding the failure mechanism relevant to interface morphology and TGO thickness.     

Numerical study of residual stress and crack nucleation in thermal barrier coating system with plane model

The residual stresses could cause extensive damage to thermal barrier coatings and even failure. A finite element model of thermal barrier coating system had been designed to simulate the residual stresses and then to analyze the crack nucleation behavior. The distribution of normal and tangential stress components along top coat (TC) / thermally grown oxide (TGO) and TGO / bond coat (BC) interfaces are shown in this work. It is found that the maximum tensile stress along TC/TGO interface occurs in the peak region during heating-up, and that along TGO/BC interface is also located in the peak region, but during the process of cooling-down. A parameter correlating the normal stress component with corresponding tangential one was used to evaluate the interfacial cracks, indicating that cracks will initiate at the peak-off region of TC/TGO interface in the heating-up phase, but for TGO/BC interface, cracks will initiate at the peak position in the cooling-down phase.     

Finite element analysis of TGO thickness on stress distribution and evolution of 8YSZ thermal barrier coatings

According to the experimental research results of the thermally grown oxide (TGO) layered growth during the pre-oxidation process of 8 wt.% yttria-stabilized zirconia thermal barrier coating (TBC), a two-dimensional sinusoidal TC/bonding coat (BC) curve interface model of the longitudinal section of TBCs based on finite element simulation was constructed; the thickness and composition of the TGO layer relative to the TC/BC curve interfacial stress distribution and its evolution during the thermal cycling process were studied. The results show that when the TGO layer uses α-Al₂O₃ as the main oxide (black TGO), the thicker the black TGO layer, the more uniform the stress distribution of the TC/BC interface. When the TGO layer is dominated by spinel-structured Co and Cr oxides (gray TGO), the stress “band” of the TC/BC interface is destroyed; it shows the alternating phenomenon of tensile stress zone and compressive stress zone, and after the rapid random growth of TGO, the concentrated tensile stress increased by a large jump. Affected by the thickness of the prefabricated black TGO layer, there is a limit peak in the thickness of the black TGO layer, the normal stress at the TC/BC boundary is minimized, and the magnitude of the stress change is also minimized.     

Effects of TGO growth on the stress distribution and evolution of three-dimensional cylindrical thermal barrier coatings based on finite element simulations

Based on the ultrasonic C-scan results of 8YSZ coatings after thermal cycles, three-dimensional cylindrical numerical simulations of the physical geometry model of the thermal barrier coating (TBC) sinusoidal surfaces were conducted with finite elements to estimate the stress distribution and evolution law of the top coat (TC)/bond coat (BC) interface, including the centre and edge of the specimen affected by the dynamic growth of the thermally grown oxide (TGO). The results show that when a layer of TGO is grown on the TC/BC interface, compressive stress is uniformly distributed on the TGO interface, and the stress value decreases as a function of the TGO layer thickness. When the thickness of the TGO exceeds a certain value, the compressive stress of all parts of the interface gradually changes to tensile stress; meanwhile, the edges of the model affected by the crest and trough effects of the wave are reflected in the radial and circumferential directions, especially along the axial direction, with alternating concentrated tensile and compressive stresses. TGO growth imposes a minor influence on the magnitude and distributions of the radial and circumferential stresses at the BC interface. The linear elasticity, creep, fatigue, and stress accumulation effects of each layer of TBCs in each thermal cycle were fully considered in this model. The model not only interprets the crest and trough effects of the TC/BC surface interface during the growth of TGO, but also interprets the effects of the core and edge of the cylindrical model, further revealing the reason for which the core and edge of the TBC will most likely form cracks.     

Effect of TGO thickness,pores, and creep on the developed residual stresses in thermal barrier coatings under cyclic loading using SEM image-based finite element model

Thermal barrier coatings (TBC) allow the metallic internal components of gas turbine engines to operate at elevated temperatures near its melting points. Formation of thermally grown oxide (TGO) layers at the top coat (TC) and bond coat (BC) interface induces cracks in the TC that may lead to complete TBC failure due to spallation. An SEM image-based finite element (FE) model is developed using commercial finite element package ABAQUS to investigate the development of residual stresses resulting from cyclic loading of TBCs. The model includes thermo-mechanical material properties and considers the real interface between the coating layers. The model includes real pores based on an SEM image, taking advantage of image processing techniques. Effect of TC surface roughness and pores on the developed residual stresses during thermal cycling is investigated with respect to different TGO thicknesses. The analysis shows that presence of TC roughness causes stress concentration sites during heating that may force horizontal cracks to initiate and propagate with stress values that are indifferent to the TGO thickness. The pores are found to shift stress concentration regions from the TC/TGO interface to the vicinity of the pores during cooling, and that may cause horizontal cracks to start from within the TC with stresses that increase with TGO thickness. Moreover, the effect of creep for all layers on the generated residual stresses is studied. Considering creep gives lower stresses at the end of cooling, however, stress distribution remains the same with and without creep.     

Simulation of thermal barrier coating spallation induced by the initiation/growth/coalescence of internal crack and interfacial crack based on a real image model

In the furnace cycle test, the growth of oxide film leads to the propagation and coalescence of multiple cracks near the interface, which should be responsible for the spallation of thermal barrier coatings (TBCs). A TBC model with real interface morphology is created, and the near-interface large pore is retained. The purpose of this work is to clarify the mechanism of TBC spallation caused by successive initiation, propagation, and linkage of cracks near the interface during thermal cycle. The dynamic growth of thermally grown oxide (TGO) is carried out by applying a stress-free strain. The crack nucleation and arbitrary path propagation in YSZ and TGO are simulated by the extended finite element method (XFEM). The debonding along the YSZ/TGO/BC interface is evaluated using a surface-based cohesive behavior. The large-scale pore in YSZ near the interface can initiate a new crack. The ceramic crack can propagate to the YSZ/TGO interface, which will accelerate the interfacial damage and debonding. For the TGO/BC interface, the normal compressive stress and small shear stress at the valley hinder the further crack propagation. The growth of YSZ crack and the formation of through-TGO crack are the main causes of TBC delamination. The accelerated BC oxidation increases the lateral growth strain of TGO, which will promote crack propagation and coalescence. The optimization design proposed in this work can provide another option for developing TBC with high durability.     

Residual stress distribution along interfaces in thermal barrier coating system under thermal cycles

The residual interfacial stress plays an important role in crack initiating and propagating along the interface, which could result in delamination failure of the thermal barrier coatings (TBCs). In this study, the finite element model of air plasma spraying(APS) TBCs was established to assess the level and distribution of residual stress along top coat(TC)/thermally grown oxide (TGO) and bond coat (BC)/TGO interfaces under thermal cycles. Instead of using vertical stress S₂₂ in global coordinate system, the normal and tangential components in the local system along the interfaces, transformed from stress components S₁₁, S₂₂, and S₁₂ in the global one, were used to evaluate the way the cracks initiate and propagate along the interfaces. Firstly, the effect of the number of thermal cycles on residual stress was investigated. It was found that, for the TBCs model without TGO growth and crack, the impact of the number of thermal cycles on the stress is very insignificant and could be ignored. So the present study only chose to focus on the first thermal cycle. Then the influence of the TGO thickness and the interface amplitude on the normal and tangential residual stresses for both homogeneous and inhomogeneous temperature fields was explored. The results show that the TGO thickness, interface amplitude and temperature field affect the residual stress level and distribution, leading to different fracture mechanisms along TC/TGO and TGO/BC interfaces. Finally, the difference between the vertical stress in the global coordinate system and the normal stress in the local coordinate system was studied. Compared with vertical stress S₂₂, the stress components normal and tangential to the TC/TGO and TGO/BC interfaces are more appropriate to describing the stress distribution along the interfaces and predicting the propensity of crack initiating and propagating along the interfaces.     

An innovative model coupling TGO growth and crack propagation for the failure assessment of lamellar structured thermal barrier coatings

The failure of plasma-sprayed thermal barrier coating (TBC) is often caused by the coating spallation due to crack propagation. In this study, a new model with stacking lamellae is developed based on the cross-section micrograph to explore crack propagation behavior within the ceramic top coat (TC) during isothermal cycling. The dynamic growth process of thermally grown oxide (TGO) is simulated via material properties change step by step. The stress profiles in the lamellar model are first evaluated, and the pore and lamellar interface crack effects on the stress state are further explored. Then, the successive crack growth, linkage, and ultimate coating spallation process is simulated. The results show that the stress intensity in TC enhances with thermal cycling. Large stress concentration always occurs near the pore and lamellar interface crack, which can result in the incipient crack growth. Moreover, the lamellar interface crack also changes the stress distribution within the TC and at the TC/bond coat interface. The multiple crack propagation upon temperature cycling is explored, and the possible coalescence mechanism is proposed. The lamellar crack steadily propagates at the early stage. The crack length sharply increases before the occurrence of coating spallation. The simulated coat spalling path is in line with the experimental result. Therefore, the new lamellar model developed in this work is beneficial to further reveal coating failure mechanism and predict coating lifetime.     

Observation of Subcritical Spall Propagation of a Thermal Barrier Coating

Observations are reported of the room-temperature propagation of a spalling failure mode of a thermal barrier coating (TBC) from its bond coat after oxidation. The coating is a Y₂O₃-stabilized ZrO₂ coating formed by electron-beam deposition on a Ni-Co-Cr-Al-Y bond coat. The spall shape evolution and stress redistribution as the spall propagates are reported. The failure propagates primarily as an interface crack between the bond coat and the thermally grown aluminum oxide (TGO) formed on the underside of the TBC during oxidation. The observations are consistent with subcritical propagation of an interface crack between the TGO and bond coat assisted by the presence of moisture. An estimate of 9 J/m² is made of the fracture resistance in air of the interface.     

Model construction and effect of thermally grown oxide dynamic growth on distribution of thermal barrier coatings

A physical geometric model of the dynamic growth of thermally grown oxide (TGO) was established based on an analysis of the TGO growth of 8YSZ thermal barrier coatings during thermal cycling. Finite-element simulation was used to simulate the evolution law between the coating residual stress and thermal cycling, and the linear elasticity, creep effect, and stress accumulation in each thermal cycle were studied. The interface between the top coat (TC) and the bond coat (BC) was covered with a TGO layer that grew vertically and slowly in a layer-like manner. The stress in the TGO was distributed with a “layer” zonal gradient, and the TGO/BC boundaries were distributed uniformly with a large compressive stress, which decreased the TGO layer thickening. With the longitudinal rapid random TGO growth, the boundaries were subjected to a tensile stress, and a high tensile stress concentration area developed at the boundaries. The internal stress consisted of an alternating and mixed distribution of concentrated compressive and tensile stresses. The concentration area of the maximum equivalent stress was distributed in the one-layer TGO near the TC/TGO interface. When a microcrack formed at the TGO/BC boundaries, the crack was subjected to a tensile stress of different size, with a higher tensile stress at both ends, which facilitated crack expansion. Thus, the 8YSZ thermal barrier coating was prone to crack formation and expansion at the TGO/BC boundaries and in the TGO layer near the TC/TGO boundaries.     

Stresses and Microstructural Development of Thermal Barrier Coatings Using AIP/D-gun Two-Step Processing

After depositing a Ni–32Co–20Cr–8Al–0.5Y–1Si–0.03B (wt%) bond coat on a Ni-base superalloy using arc ion plating (AIP), a ceramic top coat with hollow spherical powder of ZrO₂–8 wt%Y₂O₃ (HSP–YSZ) was deposited using detonation gun (D-gun) spraying. Thermal exposure behaviors of thermal barrier coatings (TBCs) were investigated at 1100°C. The thermal growth oxides (TGO) layer thickened and became more undulated during thermal exposure. Yttrium aluminum garnet (YAG) was observed within TGO, which produced thickness imperfections and thus aided to build up out-of-plane stresses. As a result, radial cracks initiated at the defects around TGO imperfections and separation developed through crack nucleation, propagation, and coalescence at the weaker TGO/bond coat interface. With further thermal exposure, coalescence of interfacial separations created a connected crack. The TBC detached and final failure occurred at the TGO/bond coat interface, leading to spallation of TBC when cooling to ambient. The stress distributions in the TGO layer with different thermal exposure times were also measured using luminescence spectroscopy. The stresses were independent of time after a transient period from θ-Al₂O₃ to α-Al₂O₃ scale. It is suggested that the lifetime of AIP/D-gun TBCs with an HSP–YSZ coat is controlled by the initiation and linking of a sub-critical interfacial crack.     

Effect of segmented cracks on TGO growth and life of thick thermal barrier coating under isothermal oxidation conditions

This paper is devoted to a comparative study on the isothermal oxidation of thick thermal barrier coating (TTBC) with and without segmented cracks produced by atmospheric plasma spray (APS) process. Accordingly, the growth of thermally grown oxide (TGO) and its effect on the degradation of the coating were investigated. Thick top coat in both segmented crack and conventional thick TBC reduced the double layered TGO growth rate slightly. The segmented crack thick TBC demonstrated longer isothermal oxidation life in comparison with that of the conventional thick TBC at 1100 °C. The dominant failure mechanism was spallation due to lateral cracking within the TGO and/or within TBC near the TGO layer, called mixed failure. Stress, and consequently strain, induced on the TTBC due to progressive TGO growth, seems to be primarily responsible for the crack initiation and propagation leading to the coating failure. Increment of elastic energy stored within the top coat due to the increasing of TGO thickness, finally causes thick thermal barrier coating failure in high temperature isothermal oxidation.     

Vertical crack distribution effect on the TBC delamination induced by crack growth from the ceramic surface and near the interface

The propagation of vertical crack on the surface of thermal barrier coatings (TBCs) may affect the interface cracking and local spallation. This research aims to establish a TBC model incorporating multiple cracks to comprehensively understand the effects of vertical crack distribution on the coating failure. The continuous TGO growth and ceramic sintering are together introduced in this model. The influence of the vertical crack spacing and non-uniform distribution on the stress state, crack driving force, and dynamic propagation is examined. Moreover, the influence of coating thickness on the crack growth driving is also explored. The results show that large spacing will lead to early crack propagation. The uniform distribution of vertical cracks can delay the spallation. When the spacing is less than 4 times ceramic coat thickness, the cracking driving force will come in a steady-state stage with the increase of vertical crack length. Prefabrication of vertical cracks with spacing less than 0.72 mm on the coating surface can greatly decrease the strain energy. The results in this study will contribute to the construction of an advanced TBC system with long lifetime.     

Analysis of interfacial crack initiation mechanism of thermal barrier coatings in isothermal oxidation process based on interfacial stress state

In this paper, the interfacial stress state is used to analyze the interfacial crack initiation mechanism of the thermal barrier coatings (TBCs) during isothermal oxidation. The influence of thermal growth stress, initial residual stress, and creep behavior on the stress distribution is considered to have an accurate simulation result. A parameter that integrates the effects of interfacial normal and tangential stress is modified for evaluating interfacial crack initiation. It is found that, in the cooling stage, the interfacial cracks sprout at the top coat (TC)/thermally grown oxide (TGO) interface valley region and the TGO/bond coat (BC) interface peak region, which agrees with the experimental results. Furthermore, the influence of interfacial roughness on crack initiation is investigated. The result shows that different interfacial roughness affects the sprouting region of interfacial cracks and cracks within the TC layer.     

Effect of interface roughness and coating thickness on interfacial shear mechanical properties of EB-PVD yttria-partially stabilized zirconia thermal barrier coating systems

The effect of interface roughness and thickness of thermal barrier coating (TBC) on the interfacial shear mechanical properties of electron beam-physical vapor deposited (EB-PVD)-TBC was examined using as-sprayed and polished bond coats (BC) 200 μm and 500 μm TBC thickness systems, by using a barb test method. The residual compressive stress in the TBC layer from the interface to the top surface was measured, by using Raman spectroscopy. The interface toughness related to the interface roughness and the thickness of the TBC. The interface toughness was larger for the BC as-sprayed TBC system than for the BC polished TBC system. The delamination of the TBC propagated within the TBC layer adjacent to the interface for the BC as-sprayed TBC; for the BC polished TBC, this occurred at the interface between the TGO and the BC. Moreover, the interface toughness was larger in the 500 μm thickness TBC than in the 200 μm thickness TBC. The relation of interface toughness to interface roughness and thickness of the TBC was associated with the interface residual compressive stress and with the interface sliding friction during the delamination of TBC.     

Experimental and numerical life evaluation of TBCs with different BC and diffusion coating under cyclic thermal loading

This paper experimentally and numerically investigates the thermally grown oxide (TGO), lifetime, and stress values in thermal barrier coatings with different bond coat (BC), without top coat (TC), and diffusion coating under cyclic thermal loading. Scanning electron microscope (SEM) analysis shows that the atmospheric plasma spraying (APS) and high-velocity oxygen fuel (HVOF) fabricated sample has the highest and lowest TGO thickness and growth rate, respectively. The new coating with two BCs has a maximum lifetime of 102 cycles. After that, the lifetime of the coatings with HVOF-BC, diffusion coating, and APS-BC reach 84, 56, and 44 cycles, respectively. The diffusion coating does not have much effect on the TGO thickness; however, it delays the Al interdiffusion to the substrate. In the sample without the TC layer, oxygen contact with the BC layer has increased, leading to a rise in the BC oxidation rate. The numerical analysis of the stress values based on SEM images shows that the more intense TGO layer growth in APS coating caused an increase in TC layer stress values. Furthermore, the results show that the new coating with two-layer BC has the lowest stress values. The TC absence causes the loss of compressive stresses caused by TC on TGO and increases the tensile stress values in this layer.     

Stress states in plasma-sprayed thermal barrier coatings upon temperature cycling: Combined effects of creep,plastic deformation,and TGO growth

To ascertain material parameter effects on the stress states is beneficial to comprehend the crack growth behavior and delamination mechanism in thermal barrier coatings (TBCs). In this work, numerical models are established to explore the combined effects of material parameters including creep, plastic deformation, and thermally grown oxide (TGO) growth on the stress states upon temperature cycling. For all layers, thermal-physical properties reliant on temperature are incorporated into the model. The process of bond coat (BC) oxidation, namely TGO growth, is materialized by changing material properties with cycles. Based on the principle of a single variable, the residual stress states are explored using many different material combinations. The results indicate that the tensile stress in the ceramic top coat (TC) decreases with the increase in the TGO lateral strain distribution gradient. Increasing the BC yield strength or decreasing the TGO growth stress can reduce the tensile stress in TC if there is no creep in the model. When BC yield strength is relatively high (≥150 MPa), BC creep will strengthen the TC tensile stress. TGO creep can decrease the tensile stress in TC irrespective of TGO growth stress and BC creep. When TGO creep rate is higher than 10B^tgo, an exceedingly small tensile stress can always be achieved. This work could provide significant theory direction for material selection and composition control towards advanced TBCs with prolonged lifetime.