Next Article in Journal
Low-Temperature Vapor-Solid Growth of ZnO Nanowhiskers for Electron Field Emission
Next Article in Special Issue
Environmental Barrier Coatings Made by Different Thermal Spray Technologies
Previous Article in Journal
Smart, Photocatalytic and Self-Cleaning Asphalt Mixtures: A Literature Review
Previous Article in Special Issue
Structural Stabilization of Mullite Films Exposed to Oxygen Potential Gradients at High Temperatures
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Crack Initiation Criteria in EBC under Thermal Stress

1
Japan Fine Ceramics Center, Nagoya 456-8587, Japan
2
National Institute for Materials Science, Ibaraki 305-0047, Japan
3
Institute of Industrial Science, The University of Tokyo, Tokyo 153-8505, Japan
4
Institute of Biomaterials and Bioengineering, Tokyo Medical and Dental University, Tokyo 101-0062, Japan
5
Faculty of Art and Regional Design, Saga University, Saga 840-8502, Japan
*
Author to whom correspondence should be addressed.
Coatings 2019, 9(11), 697; https://doi.org/10.3390/coatings9110697
Submission received: 12 September 2019 / Revised: 14 October 2019 / Accepted: 21 October 2019 / Published: 24 October 2019
(This article belongs to the Special Issue Environmental Barrier Coatings)

Abstract

:
For design of multi-layered environmental barrier coatings (EBCs), it is essential to assure mechanical reliability against interface crack initiation and propagation induced by thermal stress owing to a misfit of the coefficients of thermal expansion between the coating layers and SiC/SiC substrate. We conducted finite element method (FEM) analyses to evaluate energy release rate (ERR) for interface cracks and performed experiment to obtain interface fracture toughness to assess mechanical reliability of an EBC with a function of thermal barrier (T/EBC; SiC/SiAlON/mullite/Yb-silicate gradient composition layer/Yb2SiO5 with porous segment structure) on an SiC/SiC substrate under thermal stress due to cooling in fabrication process. Our FEM analysis revealed that a thinner SiAlON layer and a thicker mullite layer are most suitable to reduce ERRs for crack initiation at the SiC/SiAlON, SiAlON/mullite and mullite/Yb2Si2O7 interfaces. Interface fracture tests of the T/EBC with layer thicknesses within the proposed range exhibited fracture at the SiC/SiAlON and SiAlON/mullite interfaces. We also estimated the approximate fracture toughness for the SiC/SiAlON and SiAlON/mullite interfaces and lower limit of fracture toughness for the mullite/Yb2Si2O7 interface. Comparison between ERR and fracture toughness indicates that the fabricated T/EBC possesses sufficient mechanical reliability against interface crack initiation and propagation.

1. Introduction

Silicon carbide (SiC) fiber reinforced SiC matrix (SiC/SiC) composites are one-third lighter and have approximate 100–200 K higher heat resistance than current heat-resistant super alloys (Nickel-based super alloys) [1,2,3]. Application of SiC/SiC to advanced hot-section components in airplane engines is expected for improving fuel consumption and curbing emission of carbon dioxide [2,3]. SiC/SiC composites react with oxygen to form silica, which hinders degradation of the composites at high temperature over 1373–1473 K [4,5]. In the condition of high temperature and humidity (i.e., combustion gas environment), however, the silica reacts with water vapor to form silicic acid (Si(OH4)) gas [4,5], leading to wall thinning of the composites. Thus, application of the SiC/SiC composites to the hot-section components inevitably requires environmental barrier coatings (EBCs).
To achieve superior environmental shielding performance and thermomechanical durability, EBCs consist of several layers with different shielding functions [6,7]. For using EBCs at about 1673 K, which is the durable temperature of the heat-resistant SiC fiber, it is of immense importance to select a material with the following properties as the water vapor shielding layer; significant wall-thinning resistance in the combustion gas environment, and coefficients of thermal expansion (CTEs) being close to that of the SiC/SiC. Klemm conducted burner-rig tests for various heat-resistant materials under the condition simulating the actual combustion gas environment to quantitatively evaluate relationship between CTEs and wall-thinning rate [8]. The result indicated that ytterbium silicate (Yb2SiO5; YbMS and Yb2Si2O7; YbDS) has prominent wall-thinning resistance. He reported that the resistance of YbMS is particularly outstanding and that the CTE of YbMS is slightly larger than that of the SiC/SiC substrate, which is almost equal to that of YbDS. Yb-silicate is unlikely to fracture by phase transition during thermal cycle because phase transition causing the large volume change does not appear in the material throughout a wide temperature range from room temperature to high temperature. On the basis of these features, Yb-silicate is considered to be an attractive material as the water vapor shielding layer. By the way, YbMS is expected to be applied as the thermal barrier coating (TBC) for SiC/SiC substrate because the thermal conductivity of YbMS [9] is smaller than that of Y2O3-ZrO2 (YSZ), which is used as TBC for Nickel-based super alloys.
Under the combustion gas environment, EBCs are exposed to a steep gradient of oxygen chemical potential (dµO). It is known that, in an oxide layer under the environment of high temperature and dµO, the oxide ion migrates from high oxygen partial pressure (PO₂) to low PO₂ while the cation migrates in the opposite direction following the Gibbs-Duhem equation. This suggests that oxygen permeation occurs because of migration of cations as well as oxide ions. In addition, the cation migration may promote phase decomposition and fracture of the multi-layer component. For improving the phase stability of EBCs, therefore, it is of primary importance to prevent the migration of cations. Kitaoka et al. conducted oxygen permeability tests for representative heat-resistant oxides (Al6Si2O13 (mullite) and Yb-silicate) under high temperature [10,11,12]. They clarified that the oxygen shielding of mullite was superior to that of Yb-silicate and that the mullite was decomposed due to the motion of Al ions from the low-PO₂ to high-PO₂ side producing deficiency of Al ions in the layer on the low-PO₂ side. Kitaoka et al. demonstrated that the phase of mullite is stabilized by depositing a SiAlON layer on the low-PO₂ side of the mullite layer because it feeds Al ions [13].
On the basis of the abovementioned studies, we propose an EBC with a function of thermal barrier as shown in Figure 1, which we call T/EBC hereafter. The structure consists of layers as follows (in order from the side of SiC/SiC substrate); the SiC layer for smoothing the surface of the substrate, the SiAlON layer for bonding with the substrate and ensuring structural stability of the mullite, the mullite dense layer for oxygen shielding, the Yb-silicate dense layer for water vapor shielding, and the layer of YbMS with porous segment structure for thermal shock resistance. Here, for decreasing thermal stress that occurs by the difference in CTEs of the coating layers and the substrate, the Yb-silicate dense layer is formed with gradient composition to have YbDS and YbMS on the substrate and top sides, respectively. Hereafter it is referred to as the Yb-silicate gradient composition layer. The topmost YbMS layer with a porous segment structure (topcoat) is employed for thermal shock resistance so as to relax strain due to rapid temperature change in the thickness direction associated with start and stop of the engine.
T/EBC is fabricated by depositing these layers on a SiC/SiC substrate at high temperature and cooling it to room temperature. During the cooling process, the thermal stress may cause interface crack initiation/propagation and delamination along interfaces between the layers with the function of environmental barrier. This leads to deterioration of the mechanical reliability of T/EBC. Thus, it is essential to design T/EBC in such a way that crack initiation and propagation are prevented during the cooling process.
In general, fracture of brittle materials follows the Griffith theory; i.e., a crack initiates or propagates when the energy release rate (ERR) exceeds fracture toughness for the corresponding fracture mode. While fracture toughness is a material property that should be evaluated in experiments, ERR can be calculated by numerical simulations through, for example, finite element method (FEM) calculations. Unlike conventional TBCs having simple structures, the mechanical behavior of T/EBC is expected to be complicated due to its complex structure consisting of multiple layers on a substrate. Numerical simulation can be an efficient tool to evaluate ERR even in such a complex structure.
In designing T/EBC structure with sufficient mechanical reliability, we must optimize a wide range of parameters, including conditions of multi-layer deposition processes and structural dimensions such as layer thicknesses. Layer thicknesses should affect much the mechanical state and play a crucial role in the mechanical stability of T/EBC under thermomechanical loading, which the component should undergo in operation. It is of importance to understand how such structural parameters influence the mechanical reliability. The aim of this study is to propose the condition of layer thicknesses for preventing interface crack initiation and propagation during the cooling process in the fabrication of the T/EBC. We perform FEM analysis for a T/EBC model with varying thicknesses of the coating layers to investigate the dependence of the layer thicknesses on ERR of crack initiation at interfaces, where thermal stress concentrates significantly. The analysis result facilitates the determination of layer thicknesses to prevent interface crack initiation. Then, the T/EBC is fabricated within the proposed thicknesses of the coating layers. We calculate ERR for interface crack initiation and propagation by conducting FEM analysis to a model of the fabricated T/EBC. Comparing the ERRs and fracture toughness obtained by interface fracture tests, we assess the validity of the proposed layer thicknesses.

2. Experiment and Simulation Procedures

2.1. Fabrication Procedure of T/EBC

In this study, the substrate was SiC/SiC (50 mm square and 3 mm in thickness) whose matrix was formed by chemical vapor infiltration and subsequent melt infiltration. Before T/EBC deposition, a SiC layer with an adequate thickness was deposited by chemical vapor deposition in order to smoothen the surface to be coated.
Each constituent layer of T/EBC shown in Figure 1 is a complex oxide or an oxynitride containing multiple cations. When these compounds are heated to a high temperature, the composition of vapor containing cation differs substantially from that of compounds in solid phase (incongruent vaporization). Thus, deposition methods accompanying the melting process of raw materials, such as plasma spraying, bring about the incongruent vaporization of the materials during deposition, resulting in significantly deviated compositions of coating layers.
In this study, therefore, dual electron beam-physical vapor deposition (EB-PVD) was employed to prepare layers of compounds which were evaporated incongruently. Two target raw materials located in the deposition chamber were evaporated using two electron guns regulated individually, which enabled independent control of vapor pressures of gases generated from the targets. Compound layers with arbitrary compositions were deposited on a substrate placed at a proper distance from the targets. During deposition, small amount of source gas was introduced in the coating chamber, and the coated surface was heated up to about 1473 K by irradiation of a direct diode laser (wavelength of 915 nm) to promote full crystallization and densification of each layer. After the deposition finished, the specimen was cooled from the deposition temperature to room temperature at the average cooling rate of 20 K/min.
Table 1 lists the target materials (Ingots A and B) and source gases used in the deposition of the layers. During deposition of the Yb-silicate gradient composition layer, the ratio of two electron beam powers irradiated to two targets was gradually changed with coating time to control the composition toward the thickness direction. The YbMS topcoat was deposited on a rotating sample for out-of-surface to form the porous segmented structure owing to the shadowing effect.
The cross-sectional microstructure of T/EBCs was observed using a scanning electron microscope (SEM; SU8000 SEM, HITACHI, Tokyo, Japan) with energy dispersive spectroscopy (EDS).

2.2. ERR Calculation by FEM Analysis

For ensuring the mechanical reliability of T/EBC in the fabrication procedure, we need to determine its layer thicknesses with which ERR for interface crack initiation due to the thermal stress, G init th , can be reduced so that interface crack is not likely to initiate. In general, it is difficult to fabricate T/EBC without any initial interface crack such as microcrack. Thus, under the proposed layer thicknesses for preventing interface crack initiation, evaluating the behavior of interface crack propagation is also required. This means that we need to calculate ERR for interface crack propagation due to the thermal stress, G prop th , and compare with it and interface fracture toughness.
In this study, in order to investigate the effect of layer thicknesses on G init th , we conducted two-step FEM analyses using a T/EBC model with varying layer thickness; (1) Thermal stress analysis: FEM analysis for evaluating the thermal stress state after the cooling process and (2) Interface crack introduction analysis: FEM analysis for calculating ERR by introducing interface crack in the thermal stress state that was obtained in Step (1). Then, for obtaining G prop th , we conducted the two-step FEM analyses to the T/EBC model with the optimized layer thickness to prevent interface crack initiation. ABAQUS (ver.6.14.6, Dassault Systemes, Vélizy-Villacoublay, France) was used in the calculations.

2.2.1. FEM Model of T/EBC

The state of thermal stress in T/EBC is affected by the thicknesses of the coating layers. In particular, the difference in CTEs of the SiAlON and mullite layers is large compared to those at other interfaces (CTEs are explained in Section 2.2.2), and thus the thicknesses of the SiAlON and mullite layers are expected to significantly affect the state of thermal stress after cooling. We examined the simulation models mimicking T/EBC with various thicknesses of the SiAlON and mullite layers as shown in Figure 2. Note that the FEM model used in this analysis is a model of the vicinity of interface edge of T/EBC deposited sample, which is indicated by a blue rectangle in Figure 2.
In this study, the Yb-silicate gradient composition layer was modeled with five layers with different ratios of the YbDS and YbMS as follows; YbDS:YbMS = 100:0, 80:20, 50:50, 20:80, and 0:100. For computational efficiency and simplicity, the segmented YbMS layer was regarded as a homogenized material with an anisotropic (in-plane isotropic) mechanical property equivalent to the segmented layer [14]; hereafter it is referred to as the YbMS anisotropic homogeneous layer (Figure 2). The SiC/SiC substrate was treated as an in-plane isotropic material owing to the orientation of SiC fibers. T/EBC structure was represented by a 2-D simulation model, where displacements of the left and bottom edges of the model along the x and y directions, respectively, (shown in Figure 2) were constrained. A plane strain state was assumed to the z direction. In this analysis, the thicknesses of the SiAlON and mullite layers were changed (5, 15 and 25 μm), while the thicknesses of other layers were fixed as below; the SiC, Yb-silicate gradient composition and topcoat YbMS layers were 25, 100 (20 μm times five layers) and 200 μm, respectively. The thickness and width of the SiC/SiC substrate were 3000 and 4000 μm, respectively. The regions near all interfaces, where strong stress concentration is expected, were divided into a finer mesh (0.5 μm × 0.5 μm).

2.2.2. Thermal Stress Analysis

Since the experimental temperature is decreased gradually during the cooling process in fabrication of the T/EBC as explained in Section 2.1, the temperature dependence of the elastic properties (Young’ modulus and Poisson’s ratio) and the thermal property (CTE) must be taken into account for the thermal stress analysis. In addition, it is necessary to take account of the effect of creep on the mechanical state for the conditions of the experimental temperature and the cooling rate in the present deposition process. Thus, in this study, a series of the thermal stress analyses was carried out to evaluate the thermal stress distribution in T/EBC, where the creep effect was considered. Thermal stress is accumulated in T/EBC throughout the cooling process; i.e., stress is assumed to be null at the high temperature in the deposition process, and we evaluate the thermal stress state in T/EBC after cooling to room temperature from the high temperature. The thermal stress is dominated by stress in the x-axis direction in Figure 2 mainly, and thus the loading parallel to the interface, in other words, mode II-rich thermal loading is applied to T/EBC.
The temperature condition for this analysis was determined from the experimental condition of the cooling process after deposition of T/EBC (cooling from a high temperature in the deposition process to room temperature at cooling rate of 20 K/min). A uniform temperature distribution in T/EBC was assumed during the cooling process. Although the surface temperature of the coated layers is up to 1473 K in the deposition process as described in Section 2.1, in this analysis we set the deposition temperature (the initial temperature) to be 1673 K, which is an assumed operation temperature of T/EBC and provides an estimation on the safer side. In addition, the cooling rate of 20 K/min can cause a considerable effect of creep during the cooling process. In general, it is known that creep emerges significantly at temperatures higher than half the melting point [15]. Thus, in this analysis, the stress was set to be null at the initial state of 1673 K, and the development of stress distribution was calculated at a cooling rate of 20 K/min while creep was considered within the temperature range of 1673–1173 K. Creep was assumed to be negligible for SiC and SiC/SiC substrate. Since creep was not taken into account for all materials (i.e., the analysis is independent of time) below 1173 K, the system was cooled immediately from 1173 to 303 K.
As we consider a wide temperature range, we adopted temperature-dependent material properties; Young’s modulus, E, Poisson’s ratio, ν, CTE (reference temperature: 1673 K), α’, and creep properties (creep coefficient, C, and creep index, n) listed in Table 2, Table 3, Table 4, Table 5, Table 6, Table 7 and Table 8, respectively. The elastic properties (i.e., E and ν) of the YbMS anisotropic homogeneous layer and the SiC/SiC substrate were regarded as anisotropic from the microstructure and the fiber orientation, respectively. Also, E and ν of the mullite layer were regarded as anisotropic on the basis of the experimental result of the in-plane and out-of-plane loading tests on the deposited mullite layer (along the x and y axes in Figure 2, respectively). The other layers were dealt with as isotropic materials. CTE and the creep properties of each layer and substrate were assumed to be isotropic. The temperature-dependent material properties were calculated from experimental measurements and literature (see Appendix A).

2.2.3. Interface Crack Introduction Analysis

ERR owing to the thermal stress, G th , is defined as the reduction in strain energy, Π , stored in a structure per increased crack area. When the crack area, A, increases to A + d A , ERR is evaluated as follows:
G th = d Π d A .
Thus, it is necessary to obtain the strain energies stored in T/EBC with different crack lengths to evaluate G init th ( = G th | A 0 ) and G prop th . We carried out FEM analyses for T/EBC models with various crack length under the stress condition obtained from the thermal stress analysis, and the strain energy in T/EBC was evaluated as a function of crack length.
The geometry, the boundary conditions and the mesh sizes of the model in this analysis were set in the same way as in the thermal stress analysis. The initial stress condition was derived from the stress distribution of the uncracked T/EBC model obtained by the thermal stress analysis. The temperature was kept constant at 303 K (temperature after the cooling process) during the analyses, and thus we used the material properties E and ν at 303 K shown in Table 2, Table 3, Table 4, Table 5 and Table 6. Interface cracks were introduced from the right edge of the simulation model shown in Figure 2. According to a preliminary FEM analysis, the SiC/SiAlON, SiAlON/mullite and mullite/YbDS interfaces were chosen as the objective interfaces to evaluate G init th and G prop th (the detail of the preliminary analysis is found in Appendix B).

2.2.4. Evaluation of ERR for Interface Crack Initiation and Propagation

We evaluated G init th and G prop th from the strain energy in T/EBC with different crack lengths obtained by the FEM analyses in Section 2.2.3 in the following ways.

Evaluation of G init th

G init th is defined as ERR in the limit when A approaches 0, as follows:
G init th = lim A 0 d Π d A .
Note that G init th is equal to the first-order coefficient of the Taylor expansion of Π ( A ) around A = 0 with the inverted sign. Therefore, we approximated Π ( A ) as a polynomial from the FEM analyses for the T/EBC models with several different crack lengths, and G init th was evaluated as the first-order coefficient of the polynomial with the inverted sign. Two cases are possible according to the sign of G init th : for G init th > 0 , crack initiation occurs if the interface fracture toughness, Γ , is smaller than G init th ; for G init th < 0 , no crack initiation is expected. From the result of the preliminary analysis (Appendix B), we examined crack length a = 2.5, 5.0 and 7.5 μm to evaluate G init th at the objective interfaces.

Evaluation of G prop th

ERR when a crack propagates from A to A + Δ A is given by numerical differentiation as follows:
G prop th ( A + Δ A 2 ) = d Π d A | A + Δ A / 2 Π ( A + Δ A ) Π ( A ) Δ A .
From the result of the preliminary analysis (Appendix B), the a for evaluation of G prop th was determined as 2.5, 5.0, 7.5, 10, 50, 100, 250, ..., 3500 μm in increments of 250 μm between 250 and 3500 μm.

2.3. Interface Fracture Test

Interfacial toughness was measured by an interface fracture test. The test method was a modification of the one designed for EBCs on SiC/SiC composites having a weak inter-laminar strength [16]. Figure 3 shows schematic of the test set-up.
A block of specimen with a height, L, of 4 mm and a width, w, of 3 mm was prepared using a diamond blade saw in such a way that one of the fiber directions in the woven fabric of the SiC/SiC substrate was oriented to the longitudinal side of the specimen. The surface was as cut without polishing to avoid the coating delamination during the processing. The cut surface was smooth enough for us to distinguish the interface in the cross-sectional observation by conventional optical microscopy.
The specimen was set in a steel jig like a horizontal die, the upper part of which was adjustable. The specimen was put in the die and slid horizontally from the back side by a punch as indicated by the dashed arrow 1 in Figure 3 so that only the coating protruded from the die while the substrate was buried. Then the upper part was adjusted in the vertical direction (the dashed arrow 2 in Figure 3) to fix the specimen. The jig with the fixed specimen was placed under a universal testing machine (load cell capacity: 500 N, EZ-LX, Shimadzu Corporation, Kyoto, Japan) having a steel pushing plate with a thickness of 0.5 mm and a width of 5 mm.
The edge of the coating was compressively loaded by the plate with a constant crosshead speed of 0.024 mm/min till the coating was delaminated from the substrate; thus, the mode II-rich loading at the interface was achieved with little loading to the substrate in the inter-laminar shear direction. The shear stress ought not to concentrate at only an interface between SiC/SiC substrate and T/EBC (SiC/SiAlON interface); i.e., it occurs along some interfaces in T/EBC. The pushing load, P, and the crosshead displacement, u, were measured during the test. The coating near the loading point was observed in the cross-sectional direction of the specimen by optical microscopy to determine when a crack was initiated at the loaded edge and the corresponding load, P init . The number of the tests was six.
The energy release rate, G init , associated with the initiation of an interface crack at the loaded edge under an applied load of P init is given by
G init = ( P init ) 2 2 w 2 h coat { 1 E coat 1 E sub ( 1 D + ( h coat / 2 + h sub δ ) 2 I ( h coat ) 2 ) } .
Definition of parameters and details of deriving Equation (4) are shown in Appendix C. The G init is defined as fracture toughness for the interface crack initiation, Γ init . Note that, Γ init is a nominal fracture toughness value of interface crack initiation including the effect of thermal residual stress.
After the tests, four of the specimens were mounted in resin and their cross-sections of the delaminated coating and substrate were observed; the cross-section was polished to mirror finish. SEM observation (VE–7800, Keyence corporation, Osaka, Japan) and energy dispersive X-ray (EDX) analysis (EX74120 attached to field emission SEM (JSM6500F, JEOL Ltd., Tokyo, Japan)) of the cross section were done to identify the fractured interface. The remaining two specimens were served for the analyses of fracture surfaces; optical and SEM observations and EDX analysis were done.

3. Results of T/EBC Fabrication

3.1. Layer Thickness Condition for Preventing the Objective Interface Crack Initiation

Figure 4a–c show the relationships between the SiAlON layer thickness, hSA, and G init th at the SiC/SiAlON, SiAlON/mullite and mullite/YbDS interfaces, respectively, with various mullite layer thicknesses, hMu. From Figure 4a, the ERR for interface crack initiation at the SiC/SiAlON interface, G init ,   SC / SA th , is nearly zero for any combination of layer thicknesses. Thus, a crack is unlikely to initiate at the SiC/SiAlON interface during the cooling process after deposition if thicknesses of the SiAlON and mullite layers are both within the range of 5–25 μm.
The ERR for interface crack initiation at the SiAlON/mullite interface, G init ,   SA / Mu th , depends on the thicknesses of hSA and hMu as shown in Figure 4b, which is found to be expressed as the following equation within the examined range of hSA and hMu (5–25 µm):
G init , SA / Mu th ( h SA , h Mu ) = k = 0 2 l = 0 2 c k l ( h SA ) k ( h Mu ) l
with the coefficients, c k l , listed in Table 9. Figure 5 shows the contour chart of G init ,   SA / Mu th as a function of hSA and hMu shown in Equation (5). The level of G init ,   SA / Mu th is indicated by the depth of color (violet); a deeper-colored region corresponds to a lower G init ,   SA / Mu th and thus a better combination of layer thicknesses. This result suggests that one of the SiAlON or mullite layers must be thick and the other be thin in order to decrease G init ,   SA / Mu th and hinder crack initiation at the SiAlON/mullite interface.
As shown in Figure 4c, it is found that a thin mullite layer (hMu = 5 µm) results in a rise in the ERR for interface crack initiation at the mullite/YbDS interface, G init ,   Mu / YD th , with increasing hSA, while a thicker mullite layer (hMu = 15 and 25 µm) reduces G init ,   Mu / YD th to almost zero regardless of hSA. Therefore, a thick mullite layer is necessary to decrease G init ,   Mu / YD th , i.e., to suppress crack initiation at the mullite/YbDS interface.
From the above results, it is found that a thin SiAlON layer and a thick mullite layer are required for reducing G init th and preventing crack initiation at the SiC/SiAlON, SiAlON/mullite and mullite/YbDS interfaces, where strong stress concentration is expected during the cooling process in fabrication of T/EBC as shown in Figure 1.

3.2. Phase Structure of T/EBC

In Section 3.1, we found that thinning the SiAlON layer together with thickening the mullite layer was effective in suppressing interface crack formation during cooling where a large thermal stress is expected, i.e., the SiC/SiAlON, SiAlON/mullite and mullite/YbDS interfaces. In this study, therefore, the thickness of the SiAlON layer and the mullite layer were determined to be 5 and 20 µm, respectively, and T/EBC was prepared by the dual EB-PVD process described in Section 2.1. We confirmed no delamination at all the interfaces of the fabricated T/EBC. Cross-sectional SEM images and elemental mapping images of T/EBC were shown in Figure 6a,b, respectively, where we found the thicknesses of constituent layers to be; SiAlON: about 5 μm, mullite: about 20 μm, Yb-silicate gradient composition layer: about 100 μm, and porous segmented YbMS layer: about 200 μm. That is, T/EBC was successfully fabricated nearly as designed.

4. Behavior of Interface Crack Initiation and Propagation in T/EBC

To investigate the behavior of crack initiation and propagation during the cooling process in fabrication of T/EBC shown in Section 3.2, calculated ERRs ( G init th , G prop th ) must be compared with interface fracture toughness obtained by experiment. Thus, we carried out the interface fracture toughness tests (Section 2.3) to obtain Γ , which is to be compared with G init th and G prop th with hSA = 5 µm and hMu = 20 µm.

4.1. G i n i t t h and G p r o p t h for Objective Interface Cracks in Fabricated T/EBC

Table 10 lists G init th of the objective interface cracks for the T/EBC model with hSA = 5 µm and hMu = 20 µm. From this result, the SiAlON/mullite interface is found to possess the highest G init th of the objective interfaces.
Figure 7 shows the relationships between G prop th and a at the objective interfaces in the T/EBC with hSA = 5 µm and hMu = 20 µm. In the G prop th a relationships at all the objective interfaces, we observe three stages with increasing a: (A) G prop th shows a relatively rapid increase; (B) G prop th decreases; (C) G prop th increases again.
The reason for the decrease in G prop th at Stage B is explained by the distribution of out-of-plane thermal stress in the vicinity of the interface edge. Figure 8a shows the σy distributions in the coating layers along the y axis, which are obtained at the positions rx distant from the right interface edge in the x direction (rx = 0, 5, 10, 20, 30 and 40 μm). The highest level of σy was observed at the interface edge (rx = 0 μm), and the σy distribution reduces immediately to a negligible level within rx = 30 μm. Thus, the effect of σy on the coating layers is significant but extremely localized in the vicinity of interface edge, while that is negligible inside the simulation model. Note that the effect of σx on G prop th is marginal for sufficiently small a because we find σx ~ 0 at the interface edge (rx = 0 μm) directly from the balance of stress component in the x direction. Figure 8b shows the σx distributions in the coating layers along the y axis at various rx. These results indicate that σy has a dominating effect on G prop th for very short cracks. The mechanism of the second rise in G prop th at Stage C is attributed to the increase in the level of σx inside the simulation model as shown in Figure 8b. Therefore, the three-stage behavior in the G prop th a relationship is explained as follows: in the beginning, G prop th rises as a crack propagates owing to the strong σy near the interface edge (Stage A); then G prop th decreases because of a steep drop in σy inside the model (Stage B); for a sufficiently long crack, G prop th is mainly affected by σx and rises again with increasing σx (Stage C).

4.2. Interface Fracture Toughness of T/EBC

Figure 9 shows a typical Pu curve during the tests. By the cross-sectional optical microscopy, the onset of interface fracture near the loaded edge was identified at Point A in Pu curve before the maximum load. The load at Point A is defined as P init . The crack then propagated down through the interface to the midway till the load reached the maximum; finally it grew unstably to the complete delamination at the maximum load. Buckling or compressive fracture of the coating was not observed during the test. The measured P init values for respective tests are listed in Table 11.
Examples of the cross-sectional SEM images of the specimens after the tests were shown in Figure 10. At the loaded edge where the interface crack started to propagate, fracture occurred below or above the SiAlON layer, i.e., either at the interface between the SiC and SiAlON layers (Figure 10a) or at the interface between the SiAlON and mullite layers (Figure 10b). The crack path was along either of the two interfaces, and it was switched to each other by crack kinking across the SiAlON layer. However, crack propagation within a layer, such as a crack passing through the inside of layer parallel to the coating plane, was not observed.
Figure 11a shows a fracture surface near the loaded edge after the tests observed by optical microscopy. The fracture surface was classified into three according to their colors: the largest black area denoted by Region A, the gray area by Region B, and the small white area by Region C as illustrated in Figure 11a. EDX maps of a selected area in the fracture surface (indicated by a dashed rectangle in Figure 11a) are shown in Figure 11b. Mapping of element distributions shows the concentration of Si in Region A. In Region B the existence of Si, Al and N was detected. Region C was characterized by Yb concentration. Apparently Al also looks rich in Region C, but it is due to misdetection because the Mα line of Yb is very close to Kα line of Al which was used for the mapping of Al distribution. These results suggest that Region A with the largest area corresponds to the fracture surface created by the interface fracture between the SiC and SiAlON layers observed in Figure 10a; Region B showing the second largest area corresponds to the surface created by the SiAlON/mullite interface fracture (Figure 10b). We could not observe the fractured interface corresponding to Region C in the cross-sectional SEM observation of the T/EBC and substrate (Figure 10) because of the limited area of Region C. To determine the exact location of the fracture surface in Region C, we cut after the tests one of the substrate in the cross-section including Region C. The result is shown in Figure 12 suggesting that the fracture occurred in the Yb-silicate gradient composition layer.
To calculate the line fractions of Regions A, B and C along the loaded edge, the fracture surfaces near the edge were analyzed using image processing software (Image J 1.48v, National Institute of Health, Bethesda, MD, USA) (Table 12). As shown in Table 12, fracture surface was dominantly composed of Regions A and B; the ratio of Region C was relatively small. Thus the crack initiation at the edge was supposed to occur primarily either at the interface between the SiC and SiAlON layers or at the interface between the SiAlON and mullite layers.
The energy release rate for the edge crack initiation along the SiC/SiAlON interface under P init (in Table 11) was calculated from Equations (A14)–(A18) and (4), regarding the SiC/SiC plate and the SiC layer as “substrate” and the rest of the layers as “coating” in Equation (A14); the mean value of 6.4 J/m2 was obtained. When the SiAlON layer is counted as part of the substrate in addition to the SiC/SiC and SiC, the mean energy release rate for the crack along the SiAlON/mullite interface becomes 6.5 J/m2. These values of energy release rates are very close to each other because the thickness of the SiAlON layer lying between the two interfaces is quite small compared to the total thickness of the system. Under the assumption that the edge crack was initiated simultaneously through the edge (i.e., in the direction of depth) at P init regardless of the fractured interfaces, we can estimate that the SiC/SiAlON interface fracture and the SiAlON/mullite interface fracture both occur with significant proportions at almost the same energy release rates. This suggests that the interface toughness of these two interfaces are very close; both the interfaces should have an interface toughness of ~6.4 J/m2. The mean energy release rate at P init when the edge crack was initiated along the mullite/YbDS interface, which was supposed to have large thermal stress according to the FEM analysis, was calculated to be 6.8 J/m2. Since no fracture occurred at this interface at P init , we can expect that the toughness of the interface should be larger than 6.8 J/m2. The abovementioned results indicate that the approximate Γ init for the SiC/SiAlON and SiAlON/mullite interfaces are 6.4 J/m2 and the lower limit of Γ init for the mullite/YbDS interface is 6.8 J/m2.
Table 13 shows comparison between G init th of the SiC/SiAlON, SiAlON/mullite and mullite/YbDS interface cracks and Γ init of the corresponding interfaces. G init th s of every objective interface crack are smaller than Γ init . This result suggests that the effect of thermal residual stress on the intrinsic fracture toughness for interface crack initiation is sufficiently small. Therefore, in this study, we regard the nominal value as the fracture toughness of interface crack initiation. Table 13 also suggests that the SiC/SiAlON, SiAlON/mullite and mullite/YbDS interface cracks are not likely to initiate by cooling in the T/EBC fabrication process.
Figure 13a–c show comparisons between G prop th and fracture toughness for interface crack propagation, Γ prop , at the SiC/SiAlON, SiAlON/mullite and mullite/YbDS interfaces. Here, we assumed that fracture surfaces formed by interface crack initiation and propagation were the same and regarded Γ prop as equal to Γ init . As shown in Figure 13a,c, G prop th of the SiC/SiAlON and mullite/YbDS interface cracks are smaller than Γ prop of the corresponding interfaces. The results suggest that these interface cracks are not likely to propagate regardless of initial interface crack length. Figure 13b suggests that, if the T/EBC has an initial length of the SiAlON/mullite interface crack of about 1200 μm, the crack should propagate instantly because G prop th exceeds Γ prop at the crack length of 1200 µm. However, delamination along the SiAlON/mullite interface was not observed in the experiment (see Section 3.2). These results suggest that there was no initial crack exceeding the length threshold at the SiAlON/mullite interface.
The results shown in this section reveal that the T/EBC with the proposed layer thicknesses can be fabricated without delamination along interfaces by cooling in the fabrication process. This is confirmed by the result that no delamination along interface is observed as described in Section 3.2.

4.3. Future Prospects

In this study, we mainly focused on the mechanical reliability of T/EBC during the cooling process in fabrication. From the practical viewpoint, however, it is also essential to evaluate the mechanical reliability under an operating condition. In general, an EBC is exposed to a thermal cycle condition with a high humidity in operation, which can induce a microstructural change and chemical transformation of the coating layers through a reaction with heated oxygen and water vapor. Consequently, changes in the material properties of each layer due to microstructural change and chemical transformation are crucial to the mechanical state in T/EBC. Thus, it is necessary to evaluate a time dependent ERR by the FEM analysis where changes in material properties over time are incorporated. Furthermore, the interface fracture toughness changes presumably over time owing to microstructural change and chemical transformation, which should also be taken into account for a reliable design of T/EBC. Nevertheless, it is beyond the scope of this paper and will be considered in our future work.

5. Conclusions

To assess the mechanical reliability of T/EBC (SiC/SiAlON/mullite/Yb-silicate gradient composition layer/YbMS with porous segment structure) against interface crack initiation and propagation due to thermal stress induced in fabricated process, we carried out FEM analysis to evaluate ERR for interface cracks and performed experiment to obtain interface fracture toughness. Our FEM calculations revealed that G init th of the objective interfaces decreases by making the SiAlON layer thinner and the mullite layer thicker. We fabricated the T/EBC with layer thicknesses within the proposed range (hSA = 5 µm and hMu = 20 µm) and confirmed no delamination along the interfaces. In the interface fracture test for the fabricated T/EBC, fracture surfaces were found at the SiC/SiAlON and SiAlON/mullite interfaces. We estimated the approximate fracture toughness for the SiC/SiAlON and SiAlON/mullite interfaces and minimum limit of fracture toughness for the mullite/YbDS interface. Comparison between the fracture toughness by the experiment and calculated G init th and G prop th indicated that the fabricated T/EBC possesses sufficient mechanical reliability against interface cracks.
Of course, the present study considers only the stress state just after the fabrication process and does not assure the mechanical reliability in operation under thermal cycle with a high humidity, which is regarded as our future work.

Author Contributions

Conceptualization, E.K. and Y.U.; Fabrication of T/EBC, N.Y. and T.Y.; Measurement of material properties, N.Y., T.Y., H.K. and T.A.; FEM analysis, E.K. and A.K.; Interface fracture test, H.K.; Writing–original draft preparation for fabrication process and result, N.Y. and T.Y.; Writing–original draft preparation for FEM analysis and result, E.K. and A.K.; Writing–original draft preparation for interface fracture test procedure and result, H.K.; Writing–review and editing, S.K. and Y.U.; Supervision, H.K., T.A., S.K. and Y.U.; Project administration, S.K.; Funding acquisition, H.K., T.A., S.K. and Y.U.

Funding

This work was supported by Council for Science, Technology and Innovation (CSTI), Cabinet Office, Japan, through Cross-ministerial Strategic Innovation Promotion Program (SIP), “Structural Materials for Innovation” (Funding agency: JST, Japan Science and Technology Agency).

Acknowledgments

The authors would like to thank Kaori Sakata, Shoko Taniguchi and Takahumi Kawano for technical assistance.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Procedure for Calculating the Temperature-Dependent Material Properties

Appendix A.1. Young’s Modulus, Poisson’s Ratio and Shear Modulus

In this study, Young’s moduli E at 303–1673 K of the SiC and SiAlON layers were considered to be identical to that of their bulk, Ebulk. E at the same temperature range of the mullite, YbDS and YbMS layers was estimated from the measured values of coating specimens deposited by the dual EB-PVD process with the same conditions described in Section 2.1, Efilm. The outline of the estimation method of E of each layer is shown below;
Step 1: Derivation of an approximation formula of relation between temperature, T, and Ebulk of measured values or literature values.
Step 2: Estimation of Ebulk at the temperature range from 303 to 1673 K based on the approximation formula derived in Step 1.
Step 3: Measurement of Efilm at 303 K of a film/substrate layered specimen.
Step 4: Estimation of temperature dependence of the E, E ( T ) , of the coating layers using the temperature dependence of Ebulk and the ratio of E film ( 303 ) to E bulk ( 303 ) shown in Equation (A1);
E ( T ) = E bulk ( T ) E film ( 303 ) E bulk ( 303 ) .
Es of the SiC and SiAlON layers were estimated through Steps 1 and 2. The values of E of the mullite, YbDS and YbMS layers were estimated through Steps 1–4. Table A1 and Table A2 summarize the estimation procedure for each layer at Step 1 and Steps 2–3, respectively. For the Yb-silicate gradient composition layer, E was obtained from Young’s moduli of the YbMS and YbDS layers by the following Equation [17]:
E β : 1 β = ( E YM ) β ( E YD ) 1 β ,   β   =   0 ,   0.2 ,   0.5 ,   0.8 ,   and   1 .
Here, β and 1 − β denote volume content of YbMS and YbDS, respectively.
Table A1. Estimation procedure of E of SiC, SiAlON, mullite, YbDS and YbMS layers (Step 1).
Table A1. Estimation procedure of E of SiC, SiAlON, mullite, YbDS and YbMS layers (Step 1).
Coating LayerStep 1
Specimen (Bulk)Measurement MethodTemperature Range T [K]
SCWafer 1JIS R1602303–1573
SASintered body 2Resonant method 5303–973
MuSintered body 3JIS R1602303–1273
YbDSSintered body 4JIS R1602303–1273
YbMSReference [18]Reference [18]300–1600 [18]
1 CVD-SiC wafer, ADMAP Inc.. 2 Raw powder: BSI3–001B, AG Materials Inc.. Sintering condition: 2023 K for two hours, applied pressure 40MPa, PN2 = 0.6MPa. 3 Raw powder: KM-101, KCM Corporation. Sintering condition: 1573 K for fifty hours in air 2023 K for five hours in air. 4 Raw powder: self-build by ultrasonic spray pyrolysis. Sintering condition: 1673 K for two hours in air Jet mill grinding 1773 K for five hours in air. 5 Resonant method with cantilever bending and torsion mechanism. Dimensions of a specimen are 2 mm × 10 mm × 60 mm.
Table A2. Detail of estimation procedure of E of SiC, SiAlON, mullite, YbDS and YbMS layers (Steps 2 and 3).
Table A2. Detail of estimation procedure of E of SiC, SiAlON, mullite, YbDS and YbMS layers (Steps 2 and 3).
Coating LayerStep 2Step 3
Approximation FormulaSpecimen (Film/Substrate)Measurement Method
SCLinear
SALinear
Mu E = E 0 F T exp ( T 0 T )   1Mu/MuNanoindentation test 2
YbDS E = E 0 F T exp ( T 0 T ) 1YbDS/YbDSNanoindentation test 2
YbMS E = E 0 F T exp ( T 0 T ) 1YbMS/YbDSNanoindentation test 2
1 Reference [19], Optimization of E 0 ,   F ,   and   T 0 . 2 Oliver and Pharr method.
The values of ν of the SiC, SiAlON and YbDS layers were considered to be identical to those of their bulk. The outline of the estimation procedure of ν of each layer at Step 1 is shown in Table A3.
Step 1: Measurement of ν of the bulk, ν bulk in some temperature range.
Step 2: Obtainment of ν in the temperature range from 303 to 1673 K by the assumption that ν bulk outside the temperature range in Step 1 was identical to that measured at adjacent temperature.
Table A3. Estimation procedure of ν of SiC, SiAlON and YbDS layers (Step 1).
Table A3. Estimation procedure of ν of SiC, SiAlON and YbDS layers (Step 1).
Coating LayerStep 1
Specimen (Bulk)Measuring MethodTemperature Range T [K]
SCDescription in Table A1JIS R1602 and ASTM C848Description in Table A1
SADescription in Table A1Description in Table A1Description in Table A1
YbDSDescription in Table A1JIS R1602 and ASTM C848Description in Table A1
ν i j of the anisotropic mullite layer in the temperature range from 303 to 1673 K was estimated using ν i j of the bulk specimen described in Table A1 and Efilm of the mullite layer derived from the procedure shown in Table A2. ν of the YbMS used in this estimation was literature value [9] and temperature-independent. For the Yb-silicate compositional gradient layers, ν was obtained from Young’s moduli and bulk moduli, B, of YbMS and YbDS by the following equations:
ν β : 1 β = 1 2 ( 1 E β : 1 β 3 B β : 1 β ) ,     B β : 1 β = ( B YM ) β ( B YD ) 1 β .
The YbMS anisotropic homogeneous layer is characterized by anisotropic Young’s modulus, E i ani , Poisson’s ratio, ν i j ani , and shear modulus, G i j ani (i, j = x, y, z; ij). A method of calculations for these parameters will be discussed in Reference [14]. Shear modulus of mullite layer was estimated by Equation (A4).
G i j = E i E j E i + E j + 2 E j ν i j .

Appendix A.2. CTE (Reference Temperature: 1673 K)

α′ was calculated from CTE at the reference temperature of 303 K (α; evaluated as a function of temperature from experimental measurements) as follows:
α ( T ) = α ( T ) × ( T 303 ) α ( 1673 ) × ( 1673 303 ) ( T 1673 ) × ( α ( 1673 ) × ( 1673 303 ) + 1 ) .
α of the SiC and SiAlON layers was considered to be identical to that of their bulk. The values of α of the sintered bodies and the wafer in the temperature range from 573 to 1673 K were measured based on JIS R1618. Then, α of the sintered bodies and the wafer at 323–1673 K was estimated by extrapolation of the relationship between α and T. α of the mullite, YbDS and YbMS layers at given temperature was estimated using measured values of film/substrate layered specimens deposited with the same conditions as in Section 2.1. The outline of the calculation procedure of α of each layer is described below;
Step 1: Derivation of relationship between strain, ε , and T for each film/substrate layered specimen (described in Table A2), based on digital image correlation method.
Step 2: Approximation of the ε T relationship derived in Step 1 by Equation (A6).
ε = N 1 ( T 303 ) + N 2 ( T 303 ) 2 .
Step 3: Calculation of α in the temperature range from 323 to 1673 K using the approximation formula derived in Step 2.
The value of α of the Yb-silicate gradient composition layers was obtained by the law of mixture as follows [20]:
α β : 1 β = α YM β B YM + α YD ( 1 β ) B YD β B YM + ( 1 β ) B YD .
In this study, α of the YbMS anisotropic homogeneous layer was regarded as identical to that of the YbMS dense layer.

Appendix A.3. Creep Properties

In this study, we assumed that the creep property of each layer follows the strain-hardening law expressed as
ε ˙ c = C σ n ,
where ε ˙ c denotes creep strain rate.
As described in Section 2.2.2, we assumed that creep was negligible for the SiC layer. Creep properties of the SiAlON, mullite, YbDS and YbMS layers were estimated by the method described below;
Step 1: Estimation of n at a given temperature and the average of n within the temperature range from the graph of the true stress and the strain rate derived from literature values and/or measured values. Noted that the average of n was used in the FEM analysis.
Step 2: Estimation of C at a given temperature using the average of n, true stress and strain rate derived in Step 1.
Step 3: Estimation of C in the temperature range from 1173 to 1673 K using C * and Q/R in Equation (A9) based on C versus inverse of T plot derived in Step 2.
C = C * exp ( Q R T ) .
The estimation method of the creep properties of each layer at Step 1 is shown in Table A4. For the Yb-silicate gradient composition layers, C and n were obtained as follows:
C β : 1 β = ( C YM ) β ( C YD ) 1 β ,
n β : 1 β = ( n YM ) β ( n YD ) 1 β .
The creep property of the YbMS anisotropic homogeneous layer was regarded as identical to that of the YbMS dense layer.
Table A4. Method of creep property estimation of the SiAlON, mullite, YbDS and YbMS layers (Step 1).
Table A4. Method of creep property estimation of the SiAlON, mullite, YbDS and YbMS layers (Step 1).
Coating LayerStep 1
SpecimenMeasurement MethodTemperature Range T [K]
SAReference [21]Reference [21]1723–1873
MuReference [22]Reference [22]1638–1713
YbDSDescription in Table A1Uniaxial compressive creep test 21623–1723
YbMSSintered body 1Uniaxial compressive creep test 21623–1723
1 Raw powder: self-build by ultrasonic spray pyrolysis. Sintering condition: 1873 K for two hours in air Jet mill grinding 1823 K for 0.5 hours applied pressure 50 MPa in Ar, 1873 K for five hours in air. 2 Dimensions of a specimen are 3 mm × 2 mm × 2 mm.

Appendix B. Preliminary FEM Analyses: Determination of Objective Interfaces and Crack Lengths

As was explained in Introduction, an interface crack in the T/EBC is likely to occur at interfaces with strong thermal stress concentration due to the difference in the CTEs from the SiC/SiC substrate. Here we determined the objective interfaces and crack lengths from the stress distribution obtained by a preliminary FEM analysis. The thermal stress analysis was carried out in the way shown in Section 2.2.2 for the uncracked T/EBC model with the SiAlON and mullite layers both being 25 μm thick.
Figure A1 shows the distribution of σx near the interface edge in the constituent layers of the T/EBC after the cooling process. Thermal stress is found to be concentrated particularly at the SiC/SiAlON, SiAlON/mullite and mullite/YbDS interfaces. Thus, we focused on crack initiation and propagation at these objective interfaces in this study.
From the balance of the stress component in the x direction, we deduce σx ~ 0 at the interface edge and thus a dominating effect of σy on G init th . Therefore, we determined the crack lengths examined to evaluate G init th , based on the σy distribution in the vicinity of the interface edge. Figure A2 shows the σy distribution in the coating layers along the y axis at the positions rx distant from the interface edge (rx = 0, 5, 10, 60, 100, 300, 500 and 1000 μm). The highest level of σy was observed at the interface edge (rx = 0 μm) while the σy distribution reduces to a negligible level inside the model. On the basis of the results, we examined the crack lengths of a = 2.5, 5.0 and 7.5 μm to evaluate G init th . For evaluation of G prop th , the crack length a was varied as a = 2.5, 5.0, 7.5, 10, 50, 100, 250, ..., 3500 μm in increments of 250 μm between 250 and 3500 μm.
Figure A1. Thermal stress distribution near right interface edge of T/EBC deposited on SiC/SiC substrate without interface cracks.
Figure A1. Thermal stress distribution near right interface edge of T/EBC deposited on SiC/SiC substrate without interface cracks.
Coatings 09 00697 g0a1
Figure A2. σy distributions in coating layers along the y axis at the positions rx distant from right SiC/substrate interface edge in x direction.
Figure A2. σy distributions in coating layers along the y axis at the positions rx distant from right SiC/substrate interface edge in x direction.
Coatings 09 00697 g0a2

Appendix C. Derivation Process of Energy Release Rate under Applying Mechanical Loading

Figure A3 shows the loads and moments applied to the specimen. Pb and Mb ( b = 1 ,   2 ) are the load and moment from surrounding jigs, respectively, and P is the applied load by the pushing plate. Equilibria of forces and moments are given by the following equations, respectively [16,23,24]:
P P 1 P 2 = 0 ,
and
P ( h coat 2 + h sub δ ) M 1 M 2 = 0 ,
where hsub and hcoat are the thicknesses of the substrate and coating, respectively, and δ is the distance from the bottom of the substrate to the neutral axis of the specimen. When L is small compared to hsub, large P2 and M2 can be generated, but in this case these are negligible. The specimen is a multi-layer composite of ten layers including the substrate with different Young’s moduli, Ep, and thicknesses, hp ( p = 1 ,   2 ,   , 10 ) , so δ is given in a quite complicated form. Fortunately, in the loading condition for this test, the load applied to the coating is almost purely compressive parallel to the coating plane and the contribution of moment to the strain energy seems relatively small, thus the Young’s modulus of the coating, Ecoat, and the substrate, Esub, may be simply approximated with the Voigt model in the rule of mixture, respectively,
E sub = 1 h sub E s h s , E coat = 1 h coat E t h t , s + t = 10 ,
where s, t, hsub and hcoat are varied according to the fractured interface. Thus δ can be obtained as the neutral axis position of a simple two-layer composite (single equivalent coating layer with Ecoat and substrate with Esub):
δ = Λ h coat , Λ = 1 + 2 Σ η + 2 Σ η 2 2 η ( 1 + Σ η ) ,
where
Σ = E coat E sub , η = h coat h sub , E coat = E coat 1 ( ν coat ) 2 , E sub = E sub 1 ( ν sub ) 2 .
ν coat and ν sub are Poisson’s ratios of coating and substrate, respectively.
Energy release rate for the edge crack under an applied load of P is denoted by G and expressed as [24]
G = 1 2 w 2 E coat ( P 2 h coat ) + 1 2 w 2 E sub { ( P 1 ) 2 D h coat ( M 1 ) 2 I ( h coat ) 3 + ( P 2 ) 2 h sub + 12 ( M 2 ) 2 ( h sub ) 3 } .
Here, D is the dimensionless cross-sectional area and I is the moment of inertia of area written as
D = 1 η + Σ , I = Σ { ( Λ 1 η ) 2 ( Λ 1 η ) + 1 3 } + Λ η ( Λ 1 η ) + 1 η 3 .
I is again calculated assuming that the material system is a two-layer composite with a coating and substrate. Substituting Equations (A12), (A13) and P = P init into (A17) and ignoring P2 and M2, Equation (4) is derived.
Figure A3. Loads and moments applied to specimen.
Figure A3. Loads and moments applied to specimen.
Coatings 09 00697 g0a3

References

  1. Pollock, T.M.; Tin, S. Nickel-based superalloys for advanced turbine engines: Chemistry, microstructure and properties. J. Propul. Power 2006, 22, 361–374. [Google Scholar] [CrossRef]
  2. Spitsberg, I.; Steibel, J. Thermal and environmental barrier coatings for SiC/SiC CMCs in aircraft engine applications. Int. J. Appl. Ceram. Technol. 2004, 1, 291–301. [Google Scholar] [CrossRef]
  3. Padture, N.P. Advanced structural ceramics in aerospace propulsion. Nat. Mater. 2016, 15, 804–809. [Google Scholar] [CrossRef] [PubMed]
  4. Opila, E.J. Oxidation and volatilization of silica formers in water vapor. J. Am. Ceram. Soc. 2003, 86, 1238–1248. [Google Scholar] [CrossRef]
  5. Robinson, R.C.; Smialek, J.L. SiC recession caused by SiO2 Scale volatility under combustion conditions: I, experimental results and empirical model. J. Am. Ceram. Soc. 1999, 82, 1817–1825. [Google Scholar] [CrossRef]
  6. Poerschke, D.L.; Van Sluytman, J.S.; Wong, K.B.; Levi, C.G. Thermochemical compatibility of ytterbia-(hafnia/silica) multilayers for environmental barrier coatings. Acta Mater. 2013, 61, 6743–6755. [Google Scholar] [CrossRef]
  7. Richards, B.T.; Young, K.A.; de Francqueville, F.; Sehr, S.; Begley, M.R.; Wadley, H.N.G. Response of ytterbium disilicate-silicon environmental barrier coatings to thermal cycling in water vapor. Acta Mater. 2016, 106, 1–14. [Google Scholar] [CrossRef]
  8. Klemm, H. Silicon nitride for high-temperature applications. J. Am. Ceram. Soc. 2010, 93, 1501–1522. [Google Scholar] [CrossRef]
  9. Lu, M.-H.; Xiang, H.-M.; Feng, Z.-H.; Wang, X.-Y.; Zhou, Y.-C. Mechanical and thermal properties of Yb2SiO5: A promising material for T/EBCs applications. J. Am. Ceram. Soc. 2016, 99, 1404–1411. [Google Scholar] [CrossRef]
  10. Kitaoka, S.; Matsudaira, T.; Yokoe, D.; Kato, T.; Takata, M. Oxygen permeation mechanism in polycrystalline mullite at high temperatures. J. Am. Ceram. Soc. 2017, 100, 3217–3226. [Google Scholar] [CrossRef]
  11. Wada, M.; Matsudaira, T.; Kawashima, N.; Kitaoka, S.; Takata, M. Mass transfer in polycrystalline ytterbium disilicate under oxygen potential gradients at high temperatures. Acta Mater. 2017, 135, 372–381. [Google Scholar] [CrossRef]
  12. Kitaoka, S.; Matsudaira, T.; Wada, M.; Yokoi, T.; Yamaguchi, N.; Kawashima, N.; Ogawa, T.; Yokoe, D.; Kato, T. Mass transfer in Yb silicates under oxygen potential gradients at high temperatures. AMTC Lett. 2019, 6, 150–151. [Google Scholar]
  13. Kitaoka, S.; Matsudaira, T.; Kawashima, N.; Yokoe, D.; Takata, M. Structural stabilization of mullite films exposed to oxygen potential gradients at high temperatures. Coatings 2019, 9, 630. [Google Scholar] [CrossRef]
  14. Kawai, E.; Kubo, A.; Umeno, Y. Homogenization method for calculating energy release rate in environmental barrier coatings with columnar layer for ceramics. manuscript in preparation.
  15. The Society of Materials Science. Strength and Fracture of Materials, 4th ed.; The Society of Materials Science: Kyoto, Japan, 2011; p. 144. (In Japanese) [Google Scholar]
  16. Aoki, Y.; Inoue, J.; Kagawa, Y.; Igashira, K. A simple method for measurement of shear delamination toughness in environmental barrier coatings. Surf. Coat. Technol. 2017, 321, 213–218. [Google Scholar] [CrossRef]
  17. Kagawa, Y.; Hatta, H. Tailoring Ceramic Composites, 1st ed.; Agne Shofusha: Tokyo, Japan, 1990; pp. 96–97. (In Japanese) [Google Scholar]
  18. Tian, Z.; Zheng, L.; Wang, J.; Wan, P.; Li, J.; Wang, J. Theoretical and experimental determination of the major thermo-mechanical properties of RE2SiO5 (RE = Tb, Dy, Ho, Er, Tm, Yb, Lu, and Y) for environmental and thermal barrier coating applications. J. Eur. Ceram. Soc. 2016, 36, 189–202. [Google Scholar] [CrossRef]
  19. Wachtman, J.B., Jr.; Tefft, W.E.; Lam, D.G., Jr.; Apstein, C.S. Exponential temperature dependence of Young’s modulus for several oxides. Phys. Rev. 1961, 122, 1754–1759. [Google Scholar] [CrossRef]
  20. Turner, P.S. Thermal-expansion stresses in reinforced plastics. J. Res. Natl. Bur. Stand. 1946, 37, 239–250. [Google Scholar] [CrossRef]
  21. Chihara, K.; Hiratsuka, D.; Tatami, J.; Wakai, F.; Komeya, K. High-temperature deformation of α-SiAlON nanoceramics without additives. Scr. Mater. 2007, 56, 871–874. [Google Scholar] [CrossRef]
  22. Okamoto, Y.; Fukudome, H.; Hayashi, K.; Nishikawa, T. Creep deformation of polycrystalline mullite. J. Eur. Ceram. Soc. 1990, 6, 161–168. [Google Scholar] [CrossRef]
  23. Suo, Z.; Hutchinson, J.W. Interface crack between two elastic layers. Int. J. Fract. 1990, 43, 1–18. [Google Scholar] [CrossRef]
  24. Hutchinson, J.W.; Suo, Z. Mixed mode cracking in layered materials. In Advances in Applied Mechanics; Hutchinson, J.W., Wu, T.Y., Eds.; Elsevier: Amsterdam, The Netherlands, 1991; Volume 29, pp. 63–191. [Google Scholar]
Figure 1. Schematic of thermal/environmental barrier coating (T/EBC) consisting of SiAlON (SA; bonding layer), mullite (Mu; oxygen shielding), layer of gradient composition of Yb2SiO5 (YbMS) and Yb2Si2O7 (YbDS) (Yb-silicate gradient composition layer; water vapor shielding) and top YbMS with porous segment structure (thermal shock resistance).
Figure 1. Schematic of thermal/environmental barrier coating (T/EBC) consisting of SiAlON (SA; bonding layer), mullite (Mu; oxygen shielding), layer of gradient composition of Yb2SiO5 (YbMS) and Yb2Si2O7 (YbDS) (Yb-silicate gradient composition layer; water vapor shielding) and top YbMS with porous segment structure (thermal shock resistance).
Coatings 09 00697 g001
Figure 2. Schematic illustration of finite element method (FEM) analyses model and boundary condition.
Figure 2. Schematic illustration of finite element method (FEM) analyses model and boundary condition.
Coatings 09 00697 g002
Figure 3. Schematics of experimental set-up of interface fracture test: (a) Cross-sectional side view; (b) Front view.
Figure 3. Schematics of experimental set-up of interface fracture test: (a) Cross-sectional side view; (b) Front view.
Coatings 09 00697 g003
Figure 4. Relationship between SiAlON layer thickness and ERR for interface crack initiation due to thermal stress at (a) SiC/SiAlON, (b) SiAlON/mullite and (c) mullite/YbDS interfaces with various mullite layer thicknesses.
Figure 4. Relationship between SiAlON layer thickness and ERR for interface crack initiation due to thermal stress at (a) SiC/SiAlON, (b) SiAlON/mullite and (c) mullite/YbDS interfaces with various mullite layer thicknesses.
Coatings 09 00697 g004
Figure 5. Contour chart of G init , SA / Mu th ( h SA , h Mu ) interpolated with Equation (5).
Figure 5. Contour chart of G init , SA / Mu th ( h SA , h Mu ) interpolated with Equation (5).
Coatings 09 00697 g005
Figure 6. (a) Cross-sectional SEM images of T/EBC and (b) elemental mapping images focused on SiC/SiAlON/mullite layers.
Figure 6. (a) Cross-sectional SEM images of T/EBC and (b) elemental mapping images focused on SiC/SiAlON/mullite layers.
Coatings 09 00697 g006
Figure 7. Relationships between ERR for interface crack propagation due to the thermal stress and crack length at objective interfaces in T/EBC with hSA = 5 µm and hMu = 20 µm; (a) SiC/SiAlON interface: (a–1) 0 ≤ a ≤ 3500 μm and (a–2) 0 ≤ a ≤ 100 μm, (b) SiAlON/mullite interface: (b–1) 0 ≤ a ≤ 3500 μm and (b–2) 0 ≤ a ≤ 200 μm and (c) mullite/YbDS interface: (c–1) 0 ≤ a ≤ 3500 μm and (c–2) 0 ≤ a ≤ 100 μm.
Figure 7. Relationships between ERR for interface crack propagation due to the thermal stress and crack length at objective interfaces in T/EBC with hSA = 5 µm and hMu = 20 µm; (a) SiC/SiAlON interface: (a–1) 0 ≤ a ≤ 3500 μm and (a–2) 0 ≤ a ≤ 100 μm, (b) SiAlON/mullite interface: (b–1) 0 ≤ a ≤ 3500 μm and (b–2) 0 ≤ a ≤ 200 μm and (c) mullite/YbDS interface: (c–1) 0 ≤ a ≤ 3500 μm and (c–2) 0 ≤ a ≤ 100 μm.
Coatings 09 00697 g007aCoatings 09 00697 g007b
Figure 8. (a) σy distributions, and (b) σx distributions in coating layers along the y axis at the positions rx distant from right SiC/substrate interface edge in the x direction.
Figure 8. (a) σy distributions, and (b) σx distributions in coating layers along the y axis at the positions rx distant from right SiC/substrate interface edge in the x direction.
Coatings 09 00697 g008aCoatings 09 00697 g008b
Figure 9. Typical load-crosshead displacement curve during test (Specimen No. 2).
Figure 9. Typical load-crosshead displacement curve during test (Specimen No. 2).
Coatings 09 00697 g009
Figure 10. Examples of cross-sectional SEM images of specimens after interface fracture test; Fracture at (a) SiC/SiAlON interface and (b) SiAlON/mullite interface.
Figure 10. Examples of cross-sectional SEM images of specimens after interface fracture test; Fracture at (a) SiC/SiAlON interface and (b) SiAlON/mullite interface.
Coatings 09 00697 g010
Figure 11. (a) Top view of fracture surface of specimen after interface fracture test and (b) EDX maps of selected area in fracture surface indicated by a dashed rectangle in (a).
Figure 11. (a) Top view of fracture surface of specimen after interface fracture test and (b) EDX maps of selected area in fracture surface indicated by a dashed rectangle in (a).
Coatings 09 00697 g011
Figure 12. Cross-sectional SEM observation including Region C.
Figure 12. Cross-sectional SEM observation including Region C.
Coatings 09 00697 g012
Figure 13. Comparison between ERR for crack propagation due to thermal stress and fracture toughness for crack propagation at (a) SiC/SiAlON interface, (b) SiAlON/mullite interface and (c) mullite/YbDS interface.
Figure 13. Comparison between ERR for crack propagation due to thermal stress and fracture toughness for crack propagation at (a) SiC/SiAlON interface, (b) SiAlON/mullite interface and (c) mullite/YbDS interface.
Coatings 09 00697 g013
Table 1. Raw material targets and source gases for preparation of T/EBC.
Table 1. Raw material targets and source gases for preparation of T/EBC.
Coating LayerRaw Material TargetsSource Gas
Ingot AIngot B
SAAl2O3SiNH3
MuAl2O3SiO2O2
Yb-silicate gradient compositionYb2O3SiO2O2
YbMS with porous segment structureYb2O3SiO2O2
Table 2. Young’s moduli (E) of isotropic elastic layers.
Table 2. Young’s moduli (E) of isotropic elastic layers.
Temp.
T [K]
E [GPa]
SCSAYbDS:YbMS
100:080:2050:5020:800:100
30347623212812111110296.0
37347423112812111110295.9
47347222912712011010195.8
57346922712611911010195.5
67346722612511810910095.0
77346422412411710810094.5
87346222312311610799.093.8
97345922112211510698.193.0
107345722012011410597.292.1
117345421811911310496.291.2
127345121711811210395.190.2
137344921511611010294.089.2
147344621311510910192.988.1
157344421211310899.391.786.9
167344121011210698.090.485.7
Table 3. Poisson’s ratios (ν) of isotropic elastic layers.
Table 3. Poisson’s ratios (ν) of isotropic elastic layers.
Temp.
T [K]
ν
SCSAYbDS:YbMS
100:080:2050:5020:800:100
3030.210.340.260.250.250.240.23 1
3730.210.350.250.250.240.23
4730.210.350.250.250.240.23
5730.210.350.250.250.240.23
6730.210.350.250.250.240.23
7730.210.350.250.250.240.23
8730.210.360.250.250.240.23
9730.210.390.250.250.240.23
10730.210.390.240.240.240.23
11730.210.390.240.240.240.23
12730.210.390.230.230.230.23
13730.210.390.230.230.230.23
14730.200.390.230.230.230.23
15730.190.390.230.230.230.23
16730.190.390.230.230.230.23
1 No temperature dependence is assumed.
Table 4. Elastic properties (E, ν, and G) of SiC/SiC substrate.
Table 4. Elastic properties (E, ν, and G) of SiC/SiC substrate.
Temp.
T [K]
E [GPa]νG [GPa]
ExEyνxyνxzνyzGxyGxz
30315059.90.25 10.1 10.08 137.568.1
37314758.936.867.0
47314457.536.065.4
57314056.135.163.8
67313754.734.262.2
77313353.333.360.6
87313051.932.459.0
97312650.531.557.4
107312349.130.755.7
117311947.729.854.1
127311646.328.952.5
137311244.828.050.9
147310943.427.149.3
157310542.026.347.7
167310140.625.446.1
1 No temperature dependence is assumed.
Table 5. Elastic properties (E, ν, and G) of mullite layer.
Table 5. Elastic properties (E, ν, and G) of mullite layer.
Temp.
T [K]
E [GPa]νG [GPa]
ExEyνxyνxzνyzGxyGxz
3031032140.195 10.3 10.405 155.039.6
37310221354.839.4
47310221154.339.1
57310120953.838.7
67310020753.338.4
77398.920552.838.0
87398.020452.337.7
97397.020251.837.3
107396.120051.336.9
117395.119850.836.6
127394.119650.336.2
137393.219449.835.8
147392.219249.335.5
157391.219048.835.1
167390.318848.234.7
1 No temperature dependence is assumed.
Table 6. Elastic properties (E, ν, and G) of YbMS anisotropic homogenous layer.
Table 6. Elastic properties (E, ν, and G) of YbMS anisotropic homogenous layer.
Temp.
T [K]
E [GPa]νG [GPa]
ExEyνxyνxzνyzGxyGxz
3030.135960.003 10.24 10.23 10.1350.055
3730.13595.90.1350.054
4730.13595.80.1350.054
5730.13595.50.1340.054
6730.13495.00.1340.054
7730.13394.50.1330.054
8730.13293.80.1320.053
9730.13193.00.1310.053
10730.13092.10.1300.052
11730.12991.20.1280.052
12730.12790.20.1270.051
13730.12689.20.1250.051
14730.12488.10.1240.050
15730.12286.90.1220.049
16730.12185.70.1210.049
1 No temperature dependence is assumed.
Table 7. Coefficients of thermal expansion (CTEs) at reference temperature of 1673 K (α’) of substrate and coating layers.
Table 7. Coefficients of thermal expansion (CTEs) at reference temperature of 1673 K (α’) of substrate and coating layers.
Temp.
T [K]
α’ [×10−6/K]
Sub.SCSAMuYbDS:YbMSAni.
100:080:2050:5020:800:100
3235.104.993.616.254.624.754.965.205.395.39
3735.165.053.686.344.724.855.075.345.545.54
4735.285.173.816.504.905.055.315.615.845.84
5735.405.283.926.675.095.265.545.886.146.14
6735.515.374.026.845.275.465.786.156.446.44
7735.635.484.107.015.465.666.016.426.746.74
8735.755.534.157.185.655.876.256.707.047.04
9735.865.634.257.355.836.076.486.977.347.34
10735.985.734.357.526.026.286.727.247.647.64
11736.105.834.457.696.206.486.967.517.947.94
12736.215.934.557.866.396.697.197.788.248.24
13736.336.034.658.036.576.897.438.058.548.54
14736.456.134.758.206.767.097.668.338.848.84
15736.566.234.858.376.957.307.898.609.149.14
16936.706.354.978.577.177.548.188.929.509.50
Table 8. Creep properties (C and n) of coating layers.
Table 8. Creep properties (C and n) of coating layers.
Temp.
T [K]
C1
SAMuYbDS:YbMSAni.
100:080:2050:5020:800:100
11732.62 × 10−152.66 × 10−202.26 × 10−121.95 × 10−121.57 × 10−121.27 × 10−121.10 × 10−121.10 × 10−12
12232.42 × 10−149.45 × 10−191.25 × 10−111.10 × 10−119.08 × 10−127.50 × 10−126.60 × 10−126.60 × 10−12
12731.88 × 10−132.53 × 10−176.05 × 10−115.41 × 10−114.57 × 10−113.86 × 10−113.44 × 10−113.44 × 10−11
13231.25 × 10−125.29 × 10−162.60 × 10−102.35 × 10−102.03 × 10−101.75 × 10−101.59 × 10−101.59 × 10−10
13737.24 × 10−128.85 × 10−151.00 × 10−99.22 × 10−108.10 × 10−107.13 × 10−106.54 × 10−106.54 × 10−10
14233.71 × 10−111.22 × 10−133.53 × 10−93.28 × 10−92.93 × 10−92.63 × 10−92.44 × 10−92.44 × 10−9
14731.70 × 10−101.40 × 10−121.14 × 10−81.07 × 10−89.73 × 10−98.86 × 10−98.32 × 10−98.32 × 10−9
15237.04 × 10−101.37 × 10−113.40 × 10−83.23 × 10−82.99 × 10−82.76 × 10−82.62 × 10−82.62 × 10−8
15732.67 × 10−91.16 × 10−109.48 × 10−89.08 × 10−88.53 × 10−88.00 × 10−87.67 × 10−87.67 × 10−8
16239.31 × 10−98.61 × 10−102.48 × 10−72.40 × 10−72.28 × 10−72.17 × 10−72.10 × 10−72.10 × 10−7
16733.01 × 10−85.67 × 10−96.12 × 10−75.98 × 10−75.76 × 10−75.55 × 10−75.42 × 10−75.42 × 10−7
n2
SAMuYbDS:YbMSAni.
100:080:2050:5020:800:100
1.091.200.9600.9761.001.021.041.04
1 The unit of C depends on the corresponding creep index n. 2 No temperature dependence is assumed.
Table 9. Coefficients for polynomial approximation of G init , SA / Mu th ( h SA , h Mu ) in Equation (5), c k l 1.
Table 9. Coefficients for polynomial approximation of G init , SA / Mu th ( h SA , h Mu ) in Equation (5), c k l 1.
ckll
012
k02.61−2.83 × 10−18.06 × 10−3
1−2.82 × 10−13.60 × 10−2−1.03 × 10−3
27.48 × 10−3−9.68 × 10−42.77 × 10−5
1 Unit of c k l : μ m ( k + l ) .
Table 10. G init th for T/EBC model with hSA = 5 µm and hMu = 20 µm.
Table 10. G init th for T/EBC model with hSA = 5 µm and hMu = 20 µm.
Interface G init th   [ J / m 2 ]
SiC/SiAlON0
SiAlON/mullite0.324
Mullite/YbDS−0.069
Table 11. Crack initiation load ( P init ) obtained by interface fracture test.
Table 11. Crack initiation load ( P init ) obtained by interface fracture test.
Specimen Number P init   [ N ]
183
270
371
464
576
676
Average73
Table 12. Line fractions of Regions A, B and C at loaded edge.
Table 12. Line fractions of Regions A, B and C at loaded edge.
RegionLine Fraction
A0.31
B0.56
C0.13
Table 13. Comparison between G init th of target interface cracks and fracture toughness for interface crack initiation ( Γ init ) of corresponding interfaces.
Table 13. Comparison between G init th of target interface cracks and fracture toughness for interface crack initiation ( Γ init ) of corresponding interfaces.
Interface G init th   [ J / m 2 ] Γ init   [ J / m 2 ]
SiC/SiAlON06.4 1
SiAlON/mullite0.3246.4 1
Mullite/YbDS−0.0696.8 2
1 Approximate value. 2 Lower limit value.

Share and Cite

MDPI and ACS Style

Kawai, E.; Kakisawa, H.; Kubo, A.; Yamaguchi, N.; Yokoi, T.; Akatsu, T.; Kitaoka, S.; Umeno, Y. Crack Initiation Criteria in EBC under Thermal Stress. Coatings 2019, 9, 697. https://doi.org/10.3390/coatings9110697

AMA Style

Kawai E, Kakisawa H, Kubo A, Yamaguchi N, Yokoi T, Akatsu T, Kitaoka S, Umeno Y. Crack Initiation Criteria in EBC under Thermal Stress. Coatings. 2019; 9(11):697. https://doi.org/10.3390/coatings9110697

Chicago/Turabian Style

Kawai, Emi, Hideki Kakisawa, Atsushi Kubo, Norio Yamaguchi, Taishi Yokoi, Takashi Akatsu, Satoshi Kitaoka, and Yoshitaka Umeno. 2019. "Crack Initiation Criteria in EBC under Thermal Stress" Coatings 9, no. 11: 697. https://doi.org/10.3390/coatings9110697

APA Style

Kawai, E., Kakisawa, H., Kubo, A., Yamaguchi, N., Yokoi, T., Akatsu, T., Kitaoka, S., & Umeno, Y. (2019). Crack Initiation Criteria in EBC under Thermal Stress. Coatings, 9(11), 697. https://doi.org/10.3390/coatings9110697

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