Next Article in Journal
Effect of Low-Pressure Plasma Nitriding with Hollow Cathode Discharge on the Surface Microstructure of WC-Co Cermet
Next Article in Special Issue
In Situ High-Temperature Tensile Fracture Mechanism of PS-PVD EBCs
Previous Article in Journal
Preparation of High-Efficiency Flame-Retardant and Superhydrophobic Cotton Fabric by a Multi-Step Dipping
Previous Article in Special Issue
Rear Earth Oxide Multilayer Deposited by Plasma Spray-Physical Vapor Deposition for Envisaged Application as Thermal/Environmental Barrier Coating
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Experimental and Simulation Analysis of the Evolution of Residual Stress Due to Expansion via CMAS Infiltration in Thermal Barrier Coatings

1
Department of Mechanical Engineering, National Taiwan University of Science and Technology, Taipei 106335, Taiwan
2
State Key Laboratory for Strength and Vibration of Mechanical Structures, School of Aerospace Engineering, Xi’an Jiaotong University, Xi’an 710049, China
*
Authors to whom correspondence should be addressed.
Coatings 2021, 11(10), 1148; https://doi.org/10.3390/coatings11101148
Submission received: 24 August 2021 / Revised: 14 September 2021 / Accepted: 16 September 2021 / Published: 23 September 2021
(This article belongs to the Special Issue Defects, Stresses and Cracks in Thermal Barrier Coatings)

Abstract

:
The failure behavior of thermal barrier coatings (TBCs) involves multilayered systems infiltrated with calcium–magnesium–alumino-silicates (CMAS). The metastable tetragonal phase is mainly composed of 7YSZ (7 mol.% Y2O3-stabilized ZrO2), and it destabilizes into the Y-lean tetragonal phase, which may be induced by CMAS infiltration, and transforms into a monoclinic phase during cooling. The phase transformation leads to volume expansion around the CMAS-rich layer. Furthermore, it is shown that the spalling of the coating system emerges when the surface of the coating system is subjected to significant residual stress. In this study, a double-cantilever beam model is established to describe the macroscopic phenomenon of thermal buckling induced via CMAS. The result of the buckle height is used to demonstrate the consistency of the experiment and finite element simulation. The experimental parameters are imported into a multilayer cantilever beam model to analyze the interfacial stresses due to CMAS infiltration. The finite element results indicate that the phase transformation leads to damage in the coating system wherein the interfacial stresses due to phase transformation are 27% higher than those in the model without phase transformation.

1. Introduction

Thermal barrier coatings (TBCs) are widely applied in aerospace engineering to protect metallic substrates and other engine parts from high-temperature exposure. In the 7YSZ (7–8% yttria-stabilized zirconia) layer, two techniques—generation of lamellar and columnar structure—can be used to develop a TBC. Atmospheric plasma spraying (APS) coating is widely formed via a plasma jet, and the particle is spread in the lamellar microstructure. It is sensitive to thermal mismatch between ceramic and substrate. Electron beam–physical vapor deposition (EB-PVD) coating exhibits strain tolerance owing to its columnar structure. Recently, suspension plasma spraying (SPS) was used to produce a structure similar to that produced via EB-PVD. Hence, in this study, we focus on SPS, which possesses the advantage of EB-PVD and APS.
The main difference between pure zirconia and YSZ corresponds to the phase composition, and the metastable phase diagram of YSZ is shown in Figure 1a [1].
T 0 t m and T 0 c t denote the transition temperature of the tetragonal-to-monoclinic phase and the cubic-to-tetragonal phase, respectively. Undoped zirconia is characterized as an unstable material. When the temperature is under 1114 °C, zirconia exists in the monoclinic phase. When the temperature is between 1114 °C and 2369 °C, a phase transformation occurs, which leads to a stable phase, and the phase adapts to a cubic structure when the temperature exceeds 2369 °C. However, the phase transformation from the tetragonal-to-monoclinic phase during the cooling procedure [2] leads to volume expansion of 5%. This, in turn, results in a driving force that leads to microcrack propagation and interfacial delamination. Hence, many researchers focused on the stabilization of the tetragonal phase by adding oxides, such as yttria (Y2O3) and ceria (Ce2O), which exhibit a cation radius similar to that of zirconium. The main aim of the stabilization is to prevent 5% volume expansion due to the martensitic transformation during cooling. Herein, the t’ phase that can exhibit the best performance can be realized via yttria doping in the range of 7–8 wt.% (approximately 4.5 mol.%), as shown in Figure 1a by the hashed area. The t’ phase cannot transform into a monoclinic phase upon cooling, and it is often referred to as the non-transformable tetragonal phase [3].
Aircraft engines encounter particulate debris in atmosphere and volcanic ash at cruise altitude. Typically, volcanic ash consists of harmful mineral particulates, such as Al2O3, CaO, SiO2, MgO, and Fe2O3, which are often referred to as calcium–magnesium–alumino-silicates (CMAS). Thermochemical and mechanical properties of Eyjafjallajökull volcanic ash were studied by Webster et al. [4]. The thermochemical and mechanical reaction of CMAS can dramatically promote the destabilization of YSZ [5]. Phase transformation would occur, and the content of the zirconium and yttrium elements decreases with the evaporation of ZrO2 and Y2O3 during the supersonic atmospheric plasma spraying process [6]. The original YSZ is composed of 8 wt.% of Y2O3 as a phase stabilizer, and the effects of thermal exposure due to volcanic ash were examined on the original YSZ. After 10 h of thermal exposure, only 3.71 wt.% of Y2O3 was observed in the top region of YSZ. This indicated destabilization due to the molten volcanic ash that infiltrates the YSZ top region via chemical reaction: yttria is depleted by CMAS. Figure 1b introduces the process wherein YSZ is attacked by CMAS. Specifically, arrow (I) indicates that the non-transformable tetragonal phase is maintained in the stable tetragonal phase during heating. Further, arrow (II) indicates Y depletion due to the interaction between YSZ and CMAS at high temperatures. The melting CMAS dissolves the non-transformable tetragonal phase and extracts the yttrium from stabilized zirconia. This denotes that the grain boundary of zirconia becomes unstable because of the Y depletion region induced by CMAS infiltration [7]. Finally, the monoclinic phase is formed, given that arrow (III) denotes interaction with T 0 t m during cooling to room temperature, as shown in Figure 1b. The results of Raman spectra [8] indicated that the 7YSZ sample with CMAS exhibits characteristic peaks related to the monoclinic phase at approximately 600 °C during cooling.
The failure analysis in YSZ TBCs induced by CMAS infiltration has been investigated by several researchers. The degradation of EB-PVD TBCs due to the attack by CMAS deposits was examined by Peng et al. [9]. Based on the CMAS infiltration, the YSZ topcoat can be classified into two different layers: a CMAS-rich layer and a non-infiltration layer (original 7YSZ) [10]. During the thermal cycle, the CMAS-rich layer becomes denser and exhibits a higher thermal mismatch with the substrate. Based on thermodynamic laws and a mechanism-based constitutive theory, a couple of chemo-thermo-mechanical constitutive theories for TBCs under CMAS infiltration were reported by Xu et al. [11]. The buckling phenomenon in YSZ TBCs induced by CMAS was investigated by Shan et al. [12]. They discussed the significant volume expansion in the topcoat induced by CMAS infiltration. Based on the volume changes and phase transformation, the delamination and buckling of thin layers bonded to soft substrates were proposed by Reinoso et al. [13]. Li et al. discussed the temperature and stress fields in TBCs exposed to CMAS infiltration via experiment and simulation [14]. For the failure behavior of the coating system by the experiment and numerical method, the chemo-mechanical damage model on polycrystalline materials between mechanics and species diffusion was studied by Bai et al. [15]. Using the cohesive element method, Rezaei et al. [16,17] predicted the failure behavior and crack propagation of the coating system problem. Zhang et al. [18] conducted a numerical analysis and obtained the transient thermal stress in EB-PVD TBCs due to CMAS penetration. Dynamic CMAS infiltration failure analysis in TBCs was calculated by Tseng et al. [19]. Su et al. [20] numerically analyzed interfacial delamination in TBCs induced via CMAS infiltration. Naraparaju et al. [21] established a new model for estimating the CMAS infiltration depth in TBCs. The numerical study considering fluid–solid interactions between the CMAS and the topcoat was performed by Kabir et al. [22].
Generally, the failure due to expansion in TBCs induced by CMAS infiltration has not been clearly understood. Therefore, the objective of this study involves investigating the buckling mechanism and interfacial stress evolution of SPS TBCs induced via CMAS infiltration. In this study, the evolution of residual stress due to expansion via CMAS infiltration in TBCs is discussed using experiments and simulations. In the first part, a free-standing YSZ covered by CMAS is examined via experiments. Then, a simulation model of a double-cantilever beam is established with the experimental approach. Finally, the residual interfacial stress of TBCs is determined using the corresponding multilayer cantilever beam model.

2. Experiments

In the study, free-standing specimens (5 × 5 mm2) of SPS 7YSZ TBCs (approximately 450 μm thickness) were used. In the heating process of the experiment, the temperature of the dwelling time was 1400 °C, and the substrate melted during the onset temperature. The substrate was removed by applying hydrochloric acid to obtain a more precise characteristic, which was identical to that induced via CMAS infiltration. Furthermore, CMAS powder with a chemical composition of 33.2 CaO–6.4 MgO–11.9 Al2O3–48.5 SiO2 mol.% [23] was synthesized. The powder was ground using water for 6 h and placed inside an oven for 10 h. The mixing power was placed in a furnace and heated to 1350 °C for 4 h to generate CMAS glass. Hence, a fine frit powder of CMAS glass was obtained using a pestle and mortar. The CMAS powder was applied on the top surface of the free-standing SPS 7YSZ TBC specimens. The content of CMAS was controlled to 35 mg/cm2. The specimens were heated in the range of 25–1400 °C in a chamber furnace with a heating rate of 5 °C/min. Then, the specimens were heat-treated in a chamber furnace at 1400 °C for 4, 8, and 12 h, and they were cooled to room temperature with air cooling. Subsequently, cross sections of the CMAS-infiltrated SPS 7YSZ TBC specimens used in the aforementioned experiment were processed via standard metallographic techniques. These infiltrated specimens were observed via an optical microscope to determine the depth of the CMAS-rich layer and the buckle height. By scanning electron microscopy (ZEISS EVO10, Jena, Germany), the micromorphology of free-standing SPS 7YSZ TBCs infiltrated by CMAS at room temperature can be obtained, as shown in Figure 2. It can be identified with two parts. One is the CMAS-rich layer (inside the red frame), and the other is the non-infiltration layer.

Experimental Result

The infiltrated specimens were observed via an optical microscope to obtain the depth of the CMAS-rich layer and the buckle height. As shown in Figure 2, the depth of the CMAS-rich layer can be obtained from the average depth of 10 random reaction regions in the cross section infiltrated with CMAS, which is located inside the red frame. Two major experimental behaviors via different holding times were investigated. The former was the depth of the CMAS-rich layer, and the latter was the buckle height of YSZ infiltrated by CMAS. Note that each holding situation (4, 8, and 12 h) was performed with five different specimens for the experiment. Figure 3a shows the depth of the CMAS-rich layer with respect to different holding time values. The solid dots with bars and lines indicate the experimental data and the curve-fitting line. Based on previous studies [24], the depth R(Th) of the CMAS-rich layer can be defined by the parabolic law as follows:
R ( T h ) = 0.0286 T h
where Th denotes the holding time (h). Note that the depth R ( T h ) of the CMAS-rich layer follows the parabolic law by Wu et al. [25]. This study [24] also indicated that the infiltration behavior is initially faster and becomes slower as the holding time still operates.
As expected, the depth of the CMAS-rich layer was deeper when the holding time increased. The buckle height denotes the maximum distance between the bottom edge of the specimen and the horizontal ground. The buckle height of YSZ infiltrated with CMAS with respect to different holding times is shown in Figure 3b. The solid dots with bars and lines indicate the experimental data and the curve-fitting line. Notably, the buckle height of YSZ infiltrated with CMAS was more severely affected as the holding time increased.

3. Numerical Method

To investigate the infiltration process and interfacial stress evolution, two types of infiltrated coating models were used in this study, as shown in Figure 4. In Figure 4a, type I corresponds to a free-standing model, which is embedded by CMAS. The free-standing model was provided to confirm that the simulation and experimental results of the buckle height were in good agreement. The depth of the topcoat HTC was 0.2 mm. Notably, the depth R ( T h ) of the CMAS-rich layer was defined by the experimental results via Equation (1). The depth h ( T h ) of the infiltrated region was calculated by the simulation method, which will be introduced in Section 3.2. Furthermore, the depth R ( T h ) of the CMAS-rich layer included the depth of the topcoat. Type II corresponds to the TBC model and includes the topcoat, bond coat, and substrate, as shown in Figure 4b. The depths of topcoat HTC, bond coat HBC, and substrate HSUB were assumed as 0.2, 0.1, and 2 mm [26], respectively. Actually, the simulation result of interfacial stress showed the same value via different widths of the coating system, w (3.5 mm, 5 mm, 6.5 mm), which means that the width of the coating system does not generate the boundary effect. Hence, the width of the coating system, w, was set at 3.5 mm owing to saving computational efficiency. This model was established to examine the interfacial stress distribution of SPS TBCs due to CMAS infiltrating. The finite element models were constructed using the mesh with almost 45,000 elements and 414,000 elements for type I and type II models, respectively. All of the meshes were four-node quadrilateral elements (DC2D4 in the heat transfer analysis and CPE4T in the thermal structure analysis) for the topcoat, bond coat, and substrate.
To observe the buckle height behavior induced via CMAS infiltrating, the displacement was set to zero in the vertical direction at point A, while the bottom edge was fixed in the vertical direction to simulate the realistic service conditions, as shown in Figure 4. The left edge was assumed symmetric to simplify calculation processing, and the right edge was assumed with a traction-free condition. To confirm the agreement between the simulation and experimental results, heating was performed in the model from room temperature to holding temperature (25–1400 °C) at a heating rate of 5 °C/min. Subsequently, different holding times of 4, 6, 8, 10, and 12 h were introduced into the model. Specifically, the model was infiltrated with CMAS during this period. Finally, the environment in the model was cooled to room temperature at a cooling rate of 5 °C/min.

3.1. Material Properties

In the study, the material properties of TBCs, including Young’s modulus E, Poisson’s ratio ν, heat conductivity k, and coefficient of thermal expansion α, were used to characterize the elastic behavior of the TBCs, as listed in Table 1.
The effect of CMAS infiltration on the SPS coating is critical because it increases the driving force for interfacial delamination. However, it is difficult to measure the elastic and thermal properties of the compound of CMAS and the topcoat at service temperature via experiments. The structure of an infiltrated topcoat can be modeled as a laminate with the topcoat and CMAS phase. The elastic modulus of the compound can be determined by the Reuss rule-of-mixtures [27] as follows:
1 E S P S , C M A S = V S P S 1 E t Y S Z + V C M A S 1 E C M A S ,
where E S P S , C M A S denotes the effective elastic modulus of the SPS with CMAS. Additionally, E t Y S Z and E C M A S denote the elastic modulus of dense t’-YSZ and CMAS, respectively; V S P S denotes the volume fraction of the SPS coating, which is set at 90%; and V C M A S denotes ( 1 V S P S ) , which is set at 10%. In the case of the non-infiltrated topcoat, the effective Young’s modulus is defined as 40 GPa [28]. The results indicated that the effective Young’s modulus shows agreement between the experiment and the Reuss method. If the pores in the topcoat are filled with CMAS, the topcoat becomes a dense t’-YSZ with E t Y S Z = 200 GPa [29]. Given the stiffened mechanism due to the infiltration of CMAS, the effective Young’s modulus of the topcoat significantly increased from 40 to 175 GPa (non-infiltration to infiltration), as listed in Table 1. Notably, the effective Young’s modulus of the CMAS-rich layer is similar to the value evaluated via object-oriented finite elements [31]. Jackson et al. [27] also indicated that the effective modulus of the CMAS-rich layer increases to 180 GPa when the silicates crystallize and stiffen.
The thermal expansion coefficient of the infiltrated topcoat can be calculated via the Schapery rule-of-mixtures [32] as follows:
α S P S , C M A S = α C M A S E C M A S V C M A S + α t Y S Z E t Y S Z V t Y S Z E C M A S V C M A S + E t Y S Z V t Y S Z ,
Specifically, the coefficient of thermal expansion of CMAS ranged from 7–9 °C−1. Additionally, α C M A S was selected as 8.1 × 10 6 °C−1 in this analysis, and α S P S , C M A S denoted that the effective coefficient of thermal expansion from 11 × 10 6 °C−1 to 10.87 × 10 6 °C−1 (non-infiltration to infiltration), as listed in Table 1.

3.2. Expansion Induced by CMAS Infiltration

Based on our experimental results, the depth of the CMAS-rich layer can be modeled based on Equation (1). During CMAS infiltration, YSZ can be separated into a CMAS-rich layer and a non-infiltration region. Based on the element-change method USDFLD, the element properties of the topcoat are immediately transferred from the non-infiltration to the infiltration region based on Equation (1). Furthermore, experimental observations indicate that the volume of the infiltrated region is significantly affected by the separation of grain boundaries, dissolution of YSZ into CMAS, sintering of YSZ, and crystallization of molten CMAS. However, it is difficult to precisely determine the increase in the volume of the reaction region induced by CMAS infiltration at elevated temperatures. In previous studies related to the expansion model, Shaw and Duncombe [33] concluded a linear expansion at a constant rate at the beginning of the infiltration. Hence, the expansion decelerates and reaches a plateau. This denotes that equilibrium is realized with respect to expansion after the plateauing time. The depth h ( T h ) of the infiltrated region can be expressed in terms of the depth R ( T h ) of the CMAS-rich layer obtained via experiments and expansion due to CMAS infiltration, which can be related to the period of holding time. Specifically, the depth of the infiltrated region exhibits a square-root behavior as a function of the holding time as follows:
h ( T h ) = A T h + R ( T h ) ,
where A denotes the pre-exponential factor.
The thickening rate of the element row at the edge of right side, h , is as follows:
h = d ( h ( T h ) R ( T h ) ) d T h = 1 2 A t 1 / 2 ,
The growth strain rate of the edge of right side, ε , is defined as follows:
ε = h R ( T h ) = 1 2 R ( T h ) A t 1 / 2 .
The expansion in the reaction region is implemented via the subroutine UEXPAN in ABAQUS. It should be noted that expansion via CMAS infiltration provides the material behavior during the holding time for both the type I and type II models in the thermal structure analysis. The expansion strain increments of the reaction region are defined as isotropic swelling, which indicates that the element of the reaction region equally swells along the thickening and lateral directions with respect to a time increment. Conversely, the element transfers the properties by USDFLD and swells along the thickening and lateral directions when the YSZ corresponds to the reaction region based on Equation (6).

3.3. Expansion Induced via Phase Transformation

During the cooling stage, the phase transformation plays a crucial role in the thermomechanical failure induced via CMAS infiltration in 7YSZ TBCs. The experimental results indicated [6] that YSZ lacks Y3+ cations. This can lead to tetragonal–monoclinic transformation at the approximate range of 300–900 °C. This is dependent on the amount of yttrium oxide remaining in zirconia. Yashima et al. [3] examined the metastable–stable phase diagrams in the zirconia–yttrium system. The transition temperature T 0 t m is thermodynamically defined as the temperature where the free energies of monoclinic and tetragonal phases are equal. It is determined via experiments as follows:
T 0 t m = 1114 159.44 x ,
where x denotes the concentration of Y2O3 in mole percent. It is difficult to obtain the transition temperature T 0 t m , because the temperature is dependent on the microstructure of yttria-stabilized zirconia, heat treatment, and stress distribution. In the study, the range of transition temperature T 0 t m was assumed to be 600–500 °C, which indicates that the range of the concentration of Y2O3 was approximately 3.2–3.8%. It can be considered that the phase transformation started from 600 °C and was completed at 500 °C during the cooling process in the CMAS-rich layer. The phase transformation significantly affects the volume change, which corresponds to an increase of approximately 5%. In our simulation, the expansion was applied to implement the phase transformation process. Herein, the finite element model was considered for a 2-dimensional condition. Thus, the expansion value could be modified from volume to area. Hence, the value of the expansion area could be solved, and it corresponded to an increase of approximately 3%. Thus, the simulation forcibly satisfies the condition that the reaction region increases in the area by 3% when the temperature reaches the t-to-m martensitic transition temperature during the cooling stage. It should be noted that expansion via phase transformation provides the material behavior during the cooling time for both the type I and type II models in the thermal structure analysis. The increment in all the expansion strains can be implemented by UEXPAN and includes the thermal expansion strain, the expansion strain induced by CMAS infiltration, and the expansion strain due to transformation.

4. Numerical Result

Initially, it is necessary to clarify certain parameters between the simulation and the experiment. There are two major expansion factors for the CMAS-rich layer in the simulation and experimental processes. One of the factors is produced due to the separation of grain boundaries at high temperatures, and the other factor is associated with expansion induced by martensitic transformation during cooling. The latter is assumed to be 3% due to the modification from volumetric expansion (5%) to area expansion (3%) during the comparison between the simulation and experimental results of the buckle height. Based on Equation (6), the pre-exponential factor A must be determined to ensure agreement between the simulation and experimental results of the buckle height. Figure 5 shows that the buckle height was obtained via experimental and simulation results with different pre-factors. In Figure 5, solid dots and lines denote the experimental data and simulation data, respectively. In this case, the expansion rate was assumed to be 3%. A different pre-factor was selected in the simulation. The results show that increases in A indicate the occurrence of a more severe situation of grain boundary separation at high temperatures. This in turn generates a more critical buckle condition. Conversely, when A corresponds to zero, it indicates that the expansion does not occur during CMAS infiltration at high temperatures, and the buckle height value is restricted. Hence, it is concluded that pre-factor A can be defined as 0.0014 to obtain an identical buckle height value in the simulation and experiment.
Based on the aforementioned assumptions, the figure of the experiment and simulation of the buckle height induced by the infiltration of CMAS is listed in Table 2. The experimental and simulation results of YSZ induced by CMAS at room temperature for holding time values of 4 h, 8 h, and 12 h are shown in Table 2. The results suggest that the simulation solutions are in good agreement with the corresponding experimental results in terms of the buckle height of the YSZ specimens induced by CMAS infiltration.
The manner in which the two major expansion factors affect the buckle height is interesting. Figure 6 shows the buckle height of YSZ induced by CMAS infiltration with different expansion rates where the holding time is 12 h and pre-factor A is assumed to be 0.0014 in Equation (6). In this case, a thermal cycle, including heating, holding, and cooling process, is determined. During heating, the buckle height is maintained at zero because the specimen can be considered a free expansion with no CMAS infiltration. Further, the parabolic behavior of the buckle height is exhibited during the holding time. At high temperatures, CMAS gradually infiltrates YSZ specimens, and the CMAS separates the grain boundaries. Specifically, the dissolution and reprecipitation of zirconia grains occur, and the yttria-depleted zirconia phase is gradually formed. In the metastable zirconia–yttria phase diagram, arrow (II) denotes the CMAS infiltration during the holding time, as shown in Figure 1b. Additionally, the slope of the buckle height gradually decreases with respect to the depth of the CMAS-rich layer. The growth in the strain rate exhibits a parabolic behavior in the holding stage. The yttria-depleted zirconia phase (CMAS-rich layer) transforms from a tetragonal to a monoclinic phase during cooling. Arrow (III) denotes the interaction with the T0 curve during the cooling process, as shown in Figure 1b, and thus the high expansion of the CMAS-rich layer transforms to the monoclinic phase. If transformation is absent, the A–a line denotes the evolution of the buckle height during the cooling time, as shown in Figure 6. To simplify the expansion process induced via phase transformation, 600 °C and 500 °C are assumed to be the start and final times of the phase transformation. Points B and C indicate the value of the buckle height with an expansion rate of 1.5% and 3%, respectively, at a temperature of 500 °C. Hence, the B–b line and the C–c line can be considered the evolution of the buckle height with expansion rates of 1.5% and 3%, respectively, during cooling. Evidently, the buckling behavior is more harmful when the expansion rate increases, as shown in Figure 6.
The buckle height for different holding times was investigated with respect to three different expansion rates, as shown in Figure 7. In this case, the value of the buckle height was selected at room temperature, and points a, b, and c were selected, as shown in Figure 6. It should be noted that different holding times of 4, 6, 8, 10, and 12 h were implemented into the type I model. In this investigation, the expansion via phase transformation was more severe when the expansion ratio was 3%. On the contrary, the expansion ratio of 0% represents that phase transformation will not occur during the cooling step. The results indicate that the expansion rate significantly affects the value of the buckle height. When the expansion rate via phase transformation increases during the cooling stage, buckling is more dangerous. As expected, the value of the buckle height increases when the time of the free-standing YSZ infiltrated by CMAS is operated at high temperatures for a long time. This occurs because the depth of the CMAS-rich layer is deeper, and expansion due to the separation of grain boundaries and martensitic transformation is more significant.
In the second step, the residual interfacial stresses were determined at room temperature. Hence, the corresponding value, which includes expansion due to CMAS infiltration during the holding time and phase transformation during the cooling time, would transfer from the type I to the type II model. It includes the topcoat, bond coat, and substrate, as shown in Figure 4b. The maximum interfacial normal and shear stresses between topcoat and bond coat are shown in Figure 8a,b, respectively. This indicates that both interfacial normal and shear stresses are significantly affected by the phase transformation and the period of the holding time. As we know, the infiltrated region is deeper when the holding time is longer. In addition, the value of the buckle height increases when the TBCs infiltrated by the CMAS last for a long time. It can be observed that the interfacial stresses heightened owing to a larger expansion ratio and holding time. It can be speculated that increasing the expansion ratio via phase transformation and holding time would enhance the driving force for the crack evolution and interfacial delamination of the TBCs infiltrated by CMAS. This calculation indicates that the interfacial normal and shear stresses obtained via the phase transformation model are 27% higher than those obtained with the model without phase transformation. Additionally, the interfacial shear stress is more dominant than the interfacial normal stress. Previous studies [7] indicate that volume change can correspond to an increase in shear stress, which is as high as 10% and can fail the coating via spallation. Hence, the result of interfacial stress indicates that the TBCs exhibit a high probability of occurrence of interfacial delamination induced by CMAS infiltration during cooling. Hence, the shear-direction interfacial stress due to expansion by martensite transformation dominates the damage mechanism.

5. Conclusions

In this study, we examined the buckle behavior and residual stress evolution due to expansion via CMAS infiltration in SPS TBCs via experiments and simulations. First, the free-standing YSZ covered by CMAS was investigated via the experimental method. A simulation model of the double-cantilever beam via the experimental approach was established to investigate the buckle behavior induced by CMAS infiltration. The results indicate two major expansion factors of the CMAS-rich layer. The simulation method implements grain boundary separation at high temperatures and expansion induced via martensitic transformation during cooling. Finally, the interfacial stress is determined by the corresponding multilayer cantilever beam model. The result suggests that the buckle behavior is significantly affected by grain boundary separation and phase transformation expansion. The martensite phase transformation significantly affects the extent of damage to the coating system. Specifically, the interfacial normal and shear stresses are 27% with the model considering phase transformation compared to a model without phase transformation. Hence, the shear direction interfacial stress due to expansion via martensite transformation dominates the damage mechanism.

Author Contributions

Conceptualization, S.T. and D.L.; methodology, D.L.; software, S.T.; formal analysis, S.T. and D.L.; investigation, S.T.; resources, C.C. and X.F.; data curation, S.T.; writing—original draft preparation, S.T. and C.C.; writing—review and editing, S.T., C.C., and X.F.; visualization, S.T., C.C., and X.F.; project administration, C.C. and X.F. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the National Science and Technology Major Project of China (J2019-IV-0003-0070) and the National Science Foundation of China (1171101165).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author after obtaining permission from the authorized individual.

Acknowledgments

The authors are grateful to Peng Jiang for training and helpful discussion.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Evans, A.G.; Clarke, D.R.; Levi, C.G. The influence of oxides on the performance of advanced gas turbines. J. Eur. Ceram. Soc. 2008, 28, 1405–1419. [Google Scholar] [CrossRef]
  2. Schulz, U. Phase transformation in EB-PVD yttria partially stabilized zirconia thermal barrier coatings during annealing. J. Am. Ceram. Soc. 2000, 83, 904–910. [Google Scholar] [CrossRef]
  3. Yashima, M.; Kakihana, M.; Yoshimura, M. Metastable-stable phase diagrams in the zirconia-containing systems utilized in solid-oxide fuel cell application. Solid State Ion. 1996, 86, 1131–1149. [Google Scholar] [CrossRef]
  4. Webster, R.I.; Bansal, N.P.; Salem, J.A.; Opila, E.J.; Wiesner, V.L. Characterization of thermochemical and thermomechanical properties of Eyjafjallajökull volcanic ash glass. Coatings 2020, 10, 100. [Google Scholar] [CrossRef] [Green Version]
  5. Boissonnet, G.; Chalk, C.; Nicholls, J.; Bonnet, G.; Pedraza, F. Thermal insulation of CMAS (Calcium-Magnesium-Alumino-Silicates)-attacked plasma-sprayed thermal barrier coatings. J. Eur. Ceram. Soc. 2020, 40, 2042–2049. [Google Scholar] [CrossRef]
  6. Ma, G.; He, P.; Chen, S.; Kang, J.; Wang, H.; Liu, M.; Zhao, Q.; Li, G. Physicochemical properties of yttria-stabilized-zirconia in-flight particles during supersonic atmospheric plasma spray. Coatings 2019, 9, 431. [Google Scholar] [CrossRef] [Green Version]
  7. Aygun, A.; Vasiliev, A.L.; Padture, N.P.; Ma, X. Novel thermal barrier coatings that are resistant to high-temperature attack by glassy deposits. Acta Mater. 2007, 55, 6734–6745. [Google Scholar] [CrossRef]
  8. Garces, H.F.; Senturk, B.S.; Padture, N.P. In situ Raman spectroscopy studies of high-temperature degradation of thermal barrier coatings by molten silicate deposits. Scr. Mater. 2014, 76, 29–32. [Google Scholar] [CrossRef]
  9. Peng, H.; Wang, L.; Guo, L.; Miao, W.; Guo, H.; Gong, S. Degradation of EB-PVD thermal barrier coatings caused by CMAS deposits. Progress Natural Sci. Mater. Int. 2012, 22, 461–467. [Google Scholar] [CrossRef] [Green Version]
  10. Drexler, J.M.; Gledhill, A.D.; Shinoda, K.; Vasiliev, A.L.; Reddy, K.M.; Sampath, S.; Padture, N.P. Jet engine coatings for resisting volcanic ash damage. Adv. Mater. 2011, 23, 2419–2424. [Google Scholar] [CrossRef] [PubMed]
  11. Xu, G.N.; Yang, L.; Zhou, Y.C.; Pi, Z.P.; Zhu, W. A chemo-thermo-mechanically constitutive theory for thermal barrier coatings under CMAS infiltration and corrosion. J. Mech. Phys. Solids 2019, 133, 103710. [Google Scholar] [CrossRef]
  12. Shan, X.; Zou, Z.; Gu, L.; Yang, L.; Guo, F.; Zhao, X.; Xiao, P. Buckling failure in air-plasma sprayed thermal barrier coatings induced by molten silicate attack. Scr. Mater. 2016, 113, 71–74. [Google Scholar] [CrossRef]
  13. Reinoso, J.; Paggi, M.; Areias, P. A finite element framework for the interplay between delamination and buckling of rubber-like bi-material systems and stretchable electronics. J. Eur. Ceram. Soc. 2016, 36, 2371–2382. [Google Scholar] [CrossRef]
  14. Li, D.; Jiang, P.; Gao, R.; Sun, F.; Jin, X.; Fan, X. Experimental and numerical investigation on the thermal and mechanical behaviours of thermal barrier coatings exposed to CMAS corrosion. J. Adv. Ceram. 2021, 10, 551–564. [Google Scholar] [CrossRef]
  15. Bai, Y.; Santos, D.A.; Rezaei, S.; Stein, P.; Banerjee, S.; Xu, B.X. A chemo-mechanical damage model at large deformation: Numerical and experimental studies on polycrystalline energy materials. J. Solids Struct. 2021, 228, 111099. [Google Scholar] [CrossRef]
  16. Rezaei, S.; Wulfinghoff, S.; Reese, S. Prediction of fracture and damage in micro/nano coating systems using cohesive zone elements. Int. J. Solids Struct. 2017, 121, 62–74. [Google Scholar] [CrossRef]
  17. Rezaei, S.; Arghavani, M.; Wulfinghoff, S.; Kruppe, N.C.; Brögelmann, T.; Reese, S.; Bobzin, K. A novel approach for the prediction of deformation and fracture in hard coatings: Comparison of numerical modeling and nanoindentation tests. Mech. Mater. 2018, 117, 192–201. [Google Scholar] [CrossRef]
  18. Zhang, G.; Fan, X.; Xu, R.; Su, L.; Wang, T.J. Transient thermal stress due to the penetration of calcium-magnesium-alumino-silicate in EB-PVD thermal barrier coating system. Ceram. Int. 2018, 44, 12655–12663. [Google Scholar] [CrossRef]
  19. Tseng, S.; Chao, C.; Zhang, W.; Fan, X. Microcrack propagation induced by dynamic infiltration of calcium-magnesium-alumino-silicate in columnar structures for thermal barrier coatings. J. Chin. Inst. Eng. 2021, 44, 11–21. [Google Scholar] [CrossRef]
  20. Su, L.; Chen, X.; Wang, T.J. Numerical analysis of CMAS penetration induced interfacial delamination of transversely isotropic ceramic coat in thermal barrier coating system. Surf. Coat. Technol. 2015, 280, 100–109. [Google Scholar] [CrossRef] [Green Version]
  21. Naraparaju, R.; Chavez, J.J.G.; Niemeyer, P.; Hess, K.U.; Song, W.; Dingwell, D.B.; Lokachari, S.; Ramana, C.V.; Schulz, U. Estimation of CMAS infiltration depth in EB-PVD TBCs: A new constraint model supported with experimental approach. J. Eur. Ceram. Soc. 2019, 39, 2936–2945. [Google Scholar] [CrossRef]
  22. Kabir, M.R.; Sirigiri, A.K.; Naraparaju, R.; Schulz, U. Flow kinetics of molten silicates through thermal barrier coating: A numerical study. Coatings 2019, 9, 332. [Google Scholar] [CrossRef] [Green Version]
  23. Krämer, S.; Yang, J.; Levi, C.G.; Johnson, C.A. Thermochemical interaction of thermal barrier coatings with molten CaO–MgO–Al2O3–SiO2 (CMAS) deposits. J. Am. Ceram. Soc. 2006, 89, 3167–3175. [Google Scholar] [CrossRef]
  24. Mikulla, C.; Naraparaju, R.; Schulz, U.; Toma, F.L.; Barbosa, M.; Steinberg, L.; Leyens, C. Investigation of CMAS resistance of sacrificial suspension sprayed alumina topcoats on EB-PVD 7YSZ layers. J. Therm. Spray Technol. 2020, 29, 90–104. [Google Scholar] [CrossRef]
  25. Wu, D.; Zhang, H.; Shan, X.; Yang, F.; Guo, F.; Zhao, X.; Xiao, P.; Gong, S. A novel CMAS-resistant material based on thermodynamic equilibrium design: Apatite-type Gd10 (SiO4)6O3. J. Am. Ceram. Soc. 2020, 103, 3401–3415. [Google Scholar] [CrossRef]
  26. Li, B.; Fan, X.; Zhou, K.; Wang, T.J. Effect of oxide growth on the stress development in double-ceramic-layer thermal barrier coatings. Ceram. Int. 2017, 43, 14763–14774. [Google Scholar] [CrossRef]
  27. Jackson, R.W.; Zaleski, E.M.; Poerschke, D.L.; Hazel, B.T.; Begley, M.R.; Levi, C.G. Interaction of molten silicates with thermal barrier coatings under temperature gradients. Acta Mater. 2015, 89, 396–407. [Google Scholar] [CrossRef]
  28. Gupta, M.; Kumara, C.; Nylén, P. Bilayer Suspension Plasma-Sprayed Thermal Barrier Coatings with Enhanced Thermal Cyclic Lifetime: Experiments and Modeling. J. Therm. Spray Technol. 2017, 26, 1038–1051. [Google Scholar] [CrossRef] [Green Version]
  29. Levi, C.G.; Hutchinson, J.W.; Vidal-Sétif, M.H.; Johnson, C.A. Environmental degradation of thermal-barrier coatings by molten deposits. MRS Bull. 2012, 37, 932–941. [Google Scholar] [CrossRef] [Green Version]
  30. Cai, Z.; Jiang, J.; Wang, W.; Liu, Y.; Cao, Z. CMAS penetration-induced cracking behavior in the ceramic top coat of APS TBCs. Ceram. Int. 2019, 45, 14366–14375. [Google Scholar] [CrossRef]
  31. Sun, F.; Fan, X.; Zhang, T.; Jiang, P.; Yang, J. Numerical analysis of the influence of pore microstructure on thermal conductivity and Young’s modulus of thermal barrier coating. Ceram. Int. 2020, 46, 24326–24332. [Google Scholar] [CrossRef]
  32. Schapery, R.A. Thermal expansion coefficients of composite materials based on energy principles. J. Compos. Mater. 1968, 2, 380–404. [Google Scholar] [CrossRef]
  33. Shaw, T.M.; Duncombe, P.R. Forces between aluminum oxide grains in a silicate melt and their effect on grain boundary wetting. J. Am. Ceram. Soc. 1991, 74, 2495–2505. [Google Scholar] [CrossRef]
Figure 1. (a) ZrO2–Y2O3 phase diagram showing non-transformable tetragonal (t’) in the hashed area; (b) ZrO2–Y2O3 phase diagram under CMAS infiltration. The arrow denotes the interaction with CMAS. Reprinted with permission from [1] 2008 Elsevier.
Figure 1. (a) ZrO2–Y2O3 phase diagram showing non-transformable tetragonal (t’) in the hashed area; (b) ZrO2–Y2O3 phase diagram under CMAS infiltration. The arrow denotes the interaction with CMAS. Reprinted with permission from [1] 2008 Elsevier.
Coatings 11 01148 g001
Figure 2. SEM images of free-standing SPS 7YSZ TBCs infiltrated by CMAS for 4 h of holding time.
Figure 2. SEM images of free-standing SPS 7YSZ TBCs infiltrated by CMAS for 4 h of holding time.
Coatings 11 01148 g002
Figure 3. Experimental data with different holding times (4, 8, and 12 h): (a) depth of the CMAS-rich layer; (b) buckle height of YSZ induced via CMAS infiltration.
Figure 3. Experimental data with different holding times (4, 8, and 12 h): (a) depth of the CMAS-rich layer; (b) buckle height of YSZ induced via CMAS infiltration.
Coatings 11 01148 g003
Figure 4. Schematic view of two finite element models: (a) model of free-standing YSZ covered by dynamic CMAS infiltration; (b) model of TBCs covered by dynamic CMAS infiltration.
Figure 4. Schematic view of two finite element models: (a) model of free-standing YSZ covered by dynamic CMAS infiltration; (b) model of TBCs covered by dynamic CMAS infiltration.
Coatings 11 01148 g004
Figure 5. Buckle height of YSZ induced via CMAS infiltration at room temperature with different pre-factors and experiment data (the finite-element method is denoted by the spline line, and experimental method is denoted by the solid dot).
Figure 5. Buckle height of YSZ induced via CMAS infiltration at room temperature with different pre-factors and experiment data (the finite-element method is denoted by the spline line, and experimental method is denoted by the solid dot).
Coatings 11 01148 g005
Figure 6. Buckle height of YSZ induced by CMAS infiltration during a thermal cycle with different expansion rates.
Figure 6. Buckle height of YSZ induced by CMAS infiltration during a thermal cycle with different expansion rates.
Coatings 11 01148 g006
Figure 7. Buckle height of YSZ induced by CMAS infiltration for different holding times with three different expansion rates.
Figure 7. Buckle height of YSZ induced by CMAS infiltration for different holding times with three different expansion rates.
Coatings 11 01148 g007
Figure 8. Interfacial stress between topcoat and bond coat due to CMAS infiltration for different holding times with three different expansion rates: (a) max. interfacial normal stress; (b) max. shear normal stress.
Figure 8. Interfacial stress between topcoat and bond coat due to CMAS infiltration for different holding times with three different expansion rates: (a) max. interfacial normal stress; (b) max. shear normal stress.
Coatings 11 01148 g008
Table 1. Material properties of TBCs and the CMAS-rich layer [27,28,29,30].
Table 1. Material properties of TBCs and the CMAS-rich layer [27,28,29,30].
MaterialE (GPa)υk (W/mK)α (10−6 K−1)
CMAS-rich layer175.730.22.3310.87
Top coat400.21.0511
Bond coat1800.32016
Substrate2000.32014
Table 2. Comparison of thermal buckle performance obtained via (left) experimental and (right) finite element methods.
Table 2. Comparison of thermal buckle performance obtained via (left) experimental and (right) finite element methods.
Material Holding TimeExperimentFinite Element Method
4 h Coatings 11 01148 i001 Coatings 11 01148 i002
8 h Coatings 11 01148 i003 Coatings 11 01148 i004
12 h Coatings 11 01148 i005 Coatings 11 01148 i006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tseng, S.; Chao, C.; Li, D.; Fan, X. Experimental and Simulation Analysis of the Evolution of Residual Stress Due to Expansion via CMAS Infiltration in Thermal Barrier Coatings. Coatings 2021, 11, 1148. https://doi.org/10.3390/coatings11101148

AMA Style

Tseng S, Chao C, Li D, Fan X. Experimental and Simulation Analysis of the Evolution of Residual Stress Due to Expansion via CMAS Infiltration in Thermal Barrier Coatings. Coatings. 2021; 11(10):1148. https://doi.org/10.3390/coatings11101148

Chicago/Turabian Style

Tseng, Shaochen, Chingkong Chao, Dongxu Li, and Xueling Fan. 2021. "Experimental and Simulation Analysis of the Evolution of Residual Stress Due to Expansion via CMAS Infiltration in Thermal Barrier Coatings" Coatings 11, no. 10: 1148. https://doi.org/10.3390/coatings11101148

APA Style

Tseng, S., Chao, C., Li, D., & Fan, X. (2021). Experimental and Simulation Analysis of the Evolution of Residual Stress Due to Expansion via CMAS Infiltration in Thermal Barrier Coatings. Coatings, 11(10), 1148. https://doi.org/10.3390/coatings11101148

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop