Abstract

The bending-driven failure test is a reliable and efficient method for evaluating the quality and load-bearing capacity of thermal barrier coatings (TBCs). This study utilizes the discrete element method (DEM) to examine the damage evolution behavior and mechanical properties of TBCs with two different coating thicknesses under four-point bending (4PB) conditions at the microscale. The results reveal that during the bending process, both thin and thick coatings experience tensile instability fractures and the formation of transverse cracks perpendicular to the interface, with crack spacing ranging from approximately one to two times the coating thickness. Thicker coatings exhibit larger crack spacing and a significantly higher delamination damage evolution rate at the interface compared to thinner coatings, displaying more pronounced delamination characteristics. While thick coatings demonstrate stronger deformation resistance, their higher bending modulus and load-bearing capacity lead to the accumulation of more cracks under equivalent strain conditions, increasing the risk of crack propagation and failure. Additionally, the pores in the coating's microstructure promote crack branching and deflection, resulting in an expanded fractured area and a negative impact on the TBC system's lifespan. This study also analyzes variations in the load–displacement curve, particle contact states, strain energy, and acoustic emission counts. By integrating experimental results, it explores the relationship between different load stages and coating damage evolution. These findings provide a theoretical basis for identifying coating failure modes in 4PB tests and offer valuable insights for the design and optimization of TBC performance.

1 Introduction

Coating the surface of high-temperature alloy turbine blades with thermal barrier coatings (TBCs) can reduce the metal substrate temperature by 100–300 °C, significantly extending the blade's service life while enhancing engine efficiency and fuel economy [1]. In operational environments, TBCs are subjected to mechanical loads such as tension and bending, which can induce surface microcracks or delamination. These defects accelerate coating detachment and degradation. The fracture and delamination of the coating compromise its protective function, exposing the substrate material to high temperatures or corrosive media, ultimately increasing the risk of substrate failure [2]. Studying the damage behavior and fracture mechanisms of TBCs is essential for enhancing the service life and durability of components, as well as for gaining a deeper understanding of the damage patterns in coating systems.

In recent years, numerous researchers have focused on predicting the lifetime and understanding the failure mechanisms of TBCs. Studies have revealed that defects often make the interface a weak point in the mechanical performance of TBCs. The primary fracture modes observed include surface cracking, internal interface cracking within the coating, and cracking at the coating–substrate interface, particularly under high temperatures, erosion, or other extreme conditions. These fracture phenomena are largely attributed to residual stresses induced by oxidation or thermal cycling, as well as stress concentrations resulting from external loads. Previous studies [36] have discussed the effects of geometric and physical factors on coating failure caused by residual stresses. For example, Yu et al. investigated residual stresses in TBCs during the cooling process using numerical simulations and analyzed how the material properties of each layer influence these stresses [4]. Bhattacharyya et al. analyzed residual stress distribution in TBCs via plate bending analysis, comparing two material systems [5]. In addition, Refs. [710] provide studies on fracture modes induced by external loads. For instance, Li et al. analyzed the crack propagation behavior of TBCs by simulating particle erosion [7]. Zhang et al. conducted tensile tests on TBCs and studied the microscopic morphology of their cross-sectional fractures [8]. Šulák et al. conducted low cycle fatigue (LCF) tests in a symmetrical push-pull cycle under strain control at 900 °C. Cyclic hardening/softening curves, cyclic stress–strain curves, and fatigue life curves of both the TBCs-coated and uncoated materials were assessed [10].

Three-point bending (3PB) and four-point bending (4PB) tests are commonly used methods for studying the fracture behavior of brittle coatings, employed to analyze and evaluate the mechanical properties of the coating/substrate [1114], damage characteristics [1518], and influencing factors [1921]. For example, Zhao et al. improved the 4PB specimen and combined it with finite element analysis to accurately evaluate the fracture toughness of the TBCs from a fracture mechanics perspective. Through experimental analysis and theoretical calculations, they derived the bond strength required for the initiation of interface edge cracks and the fracture toughness of the system [11]. Chen et al. employed a sandwiched four-point bend specimen to evaluate the crack growth resistance in plasma-sprayed TBCs, discovering a rising crack growth resistance curve and obtaining the steady-state strain energy release rate [14]. Wang et al. employed the acoustic emission (AE) technique to characterize the failure behavior of TBCs under three-point bending tests. The results indicated that the main failure mechanisms of TBCs under three-point bending conditions include the debonding of the metal coating from the substrate and the propagation of the horizontal crack along the bond coat (BC)/substrate interface under the action of the flexural moment [17]. Eriksson et al. conducted tensile bending fatigue tests on TBCs and observed that under low-load conditions, the TBCs exhibited extensive delamination damage, while under high-load conditions, the TBCs demonstrated relatively better integrity. At the same time, they used finite element (FE) analysis to investigate potential transverse crack structures [18]. Frommherz et al. measured the delamination toughness and its sintering behavior using an improved 4PB apparatus. The research results indicate that the increase in coating stiffness is related to changes in its microstructure. The specific spray structure of TBCs has a great influence on the sinter stability which results in a characteristic temperature dependency [19]. In addition, Martins et al. conducted tests on coatings using an improved 4PB test to study the effect of the morphology curvature of the bond coat on the fracture mechanical properties of TBCs [20]. Despite the extensive research on TBCs in the aforementioned studies, the relationship between the performance of these complex coatings and their microstructure remains unclear. The discrete element method (DEM), a widely used numerical simulation technique for discontinuous media materials, can effectively represent the discrete characteristics and interactions of bulk materials. This method, by simulating particle collisions and frictional interactions, can capture the behavior of materials under complex loading conditions, making it particularly suitable for performance prediction and microstructural analysis of coating materials.

To investigate the effects of localized damage and cracks on the overall mechanical behavior of coatings, this study employs DEM and conducts a four-point bending (4PB) numerical simulation by constructing a rectangular DEM model of the coating system. The failure modes of TBCs under bending loads were analyzed at the microscale, and the simulation results were validated against existing experimental data. Additionally, the study examines two coatings with different thicknesses to evaluate the influence of coating geometry on structural mechanical properties and crack propagation behavior. The TBCs analyzed consist of a ceramic top layer (TC) made of lanthanum zirconate (LZO) and a bond coat (BC) composed of NiCrAlY, both prepared using atmospheric plasma spraying (APS), with a high-temperature alloy GH4169 (Inconel 718, Zhejiang Flysun Special Steel Co., Ltd, China) serving as the substrate.

2 Numerical Simulation of the Four-Point Bending Test for the Thermal Barrier Coatings

2.1 Calibration of the Coating's Microscopic Parameters.

The DEM model consists of a collection of discrete circular particles, where the contact forces between particles follow the force-displacement law, and the particles move based on Newton's second law. As a DEM, the bonded particle model includes a contact bond model (CBM) and a parallel bond model (PBM) [22]. Different constitutive models give rise to distinct mechanical properties in the medium. In this study, the ceramic layer material, LZO, is brittle, so the parallel bond model is used, whereas the bond coat material, NiCrAlY, and the high-temperature alloy substrate, GH4169 (Inconel 718), are plastic, and thus, the contact bond model is applied.

Since the DEM model reflects the macroscopic mechanical properties of the coating based on the mechanical properties of particles and their interaction behaviors, a systematic calibration of the model's microscopic parameters is required prior to numerical simulations to ensure accurate prediction of the macroscopic mechanical performance of TBCs. The symbols and definitions of the parameters involved in the calibration are detailed in Table 1. During the calibration process, initial values are assigned to various microscopic parameters. Numerical simulations of typical mechanical behaviors, including uniaxial compression, three-point bending, fracture toughness, and uniaxial tension, are conducted for each functional layer of the coating, as shown in Fig. 1. Subsequently, based on the obtained stress–strain curves and load–displacement curves, the equivalent macroscopic mechanical parameters of the coating are extracted, following the calculation methods provided in Eqs. (1)(5).

Fig. 1
Four calibration tests were conducted: (a) uniaxial compression, (b) three-point bending, (c) fracture toughness, and (d) uniaxial tension
Fig. 1
Four calibration tests were conducted: (a) uniaxial compression, (b) three-point bending, (c) fracture toughness, and (d) uniaxial tension
Close modal
Table 1

Micro and macro parameters used in parallel bonded DEM model

Microscopic parametersMacro parameters
Parallel bond effective modulus E¯*Elastic modulus E
Linear effective modulus E*Poisson's ratio ν
Contact stiffness ratio k*Compressive strength σbc
Parallel bond stiffness ratio k¯*Bending strength σbb
Parallel bond tensile strength σ¯cFracture toughness KIC
Parallel bond cohesion c¯Tensile strength σb
Friction factor μ
Parallel bond friction angle ϕ¯
Distribution factor of moment β¯
Quantity ratio Rf
Strength ratio Rσ
Contact bond shear strength Sσ
Contact bond tensile strength Tσ
Modulus ratio RE
Microscopic parametersMacro parameters
Parallel bond effective modulus E¯*Elastic modulus E
Linear effective modulus E*Poisson's ratio ν
Contact stiffness ratio k*Compressive strength σbc
Parallel bond stiffness ratio k¯*Bending strength σbb
Parallel bond tensile strength σ¯cFracture toughness KIC
Parallel bond cohesion c¯Tensile strength σb
Friction factor μ
Parallel bond friction angle ϕ¯
Distribution factor of moment β¯
Quantity ratio Rf
Strength ratio Rσ
Contact bond shear strength Sσ
Contact bond tensile strength Tσ
Modulus ratio RE

On this basis, the correlation between microscopic and macroscopic parameters is analyzed using Eq. (6), taking the relationship between the microscopic and macroscopic parameters of the TC layer in Table 2 as an example. Key microscopic parameters that are highly correlated with the target macroscopic properties are identified and fitted using multivariate linear regression to establish multivariate linear regression equations between microscopic parameters and macroscopic mechanical properties (e.g., Eqs. (7)(11) represent the regression models for the TC layer). The final microscopic parameters obtained from the regression analysis are subsequently incorporated into the DEM numerical model for validation simulations. The resulting macroscopic parameters from the simulations are compared with experimental data for each material to evaluate the rationality of the calibrated parameters and the accuracy of the model predictions.

Table 2

Pearson correlation coefficient of micro- and macro-parameters of TC

EνσbcσbbKIC
E¯*0.826−0.164−0.0270.0280.075
E*/E¯*0.1740.0760.1260.1100.041
k¯*−0.2980.913−0.063−0.114−0.070
σ¯c0.0110.0330.6760.6040.724
c¯/σ¯c−0.004−0.0580.3620.028−0.001
μ0.003−0.00060.1650.0780.062
ϕ¯−0.0020.0160.0470.0020.118
β¯0.024−0.072−0.285−0.103−0.242
Rf−0.3460.233−0.363−0.625−0.526
Rσ0.009−0.0640.2580.3780.152
RE0.084−0.017−0.0930.0120.073
EνσbcσbbKIC
E¯*0.826−0.164−0.0270.0280.075
E*/E¯*0.1740.0760.1260.1100.041
k¯*−0.2980.913−0.063−0.114−0.070
σ¯c0.0110.0330.6760.6040.724
c¯/σ¯c−0.004−0.0580.3620.028−0.001
μ0.003−0.00060.1650.0780.062
ϕ¯−0.0020.0160.0470.0020.118
β¯0.024−0.072−0.285−0.103−0.242
Rf−0.3460.233−0.363−0.625−0.526
Rσ0.009−0.0640.2580.3780.152
RE0.084−0.017−0.0930.0120.073
In summary, the calibration of microscopic parameters is essentially a dynamic optimization process aimed at continuously adjusting the microscopic parameters to achieve a high degree of agreement between numerical simulation results and experimental measurements. The final microscopic parameters for each material layer are listed in Table 3, and a comparison of the macroscopic mechanical parameters obtained from numerical simulations and experimental results is provided in Table 4.
(1)
(2)
(3)
(4)
(5)
(6)
(7)
(8)
(9)
(10)
(11)
where ΔεAi is the change in axial strain in the experiment, Δσi and ΔεTi are the corresponding changes in axial stress and lateral strain, respectively. Lk is the span of the support points in the bending and fracture toughness tests, Pmax is the maximum axial loading force in the test, b is the width of the test sample's cross section, h is the height of the test sample, and hl is the notch height in the fracture toughness test. r is the Pearson correlation coefficient, n is the total number of corresponding micro- or macro-parameters, Xi represents the individual micro-parameters, Yi represents the individual macro parameters, X¯ is the average value of the micro-parameters in each group, and Y¯ is the average value of the macro parameters in each group.
Table 3

Micro-parameters of numerical specimen

Hierarchical structureMicroscopic parameter
TCE¯*(GPa)E*/E¯*k¯*σ¯c(MPa)c¯/σ¯cμϕ¯()β¯RfRσRE
403.6150.16194.9871.50.1150.70.50.6580.5
BCE*(GPa)k*μTσ(Pa)Sσ(Pa)
482.029109.429111719.122850
SubstrateE*(GPa)k*μTσ(Pa)Sσ(Pa)
350313124,924950
Hierarchical structureMicroscopic parameter
TCE¯*(GPa)E*/E¯*k¯*σ¯c(MPa)c¯/σ¯cμϕ¯()β¯RfRσRE
403.6150.16194.9871.50.1150.70.50.6580.5
BCE*(GPa)k*μTσ(Pa)Sσ(Pa)
482.029109.429111719.122850
SubstrateE*(GPa)k*μTσ(Pa)Sσ(Pa)
350313124,924950
Table 4

Macroscopic parameter simulation of coating system-experimental comparison

Hierarchical structureMacroscopic parameterSimulated valueExperimental value [2326]Error
TCE(GPa)161.1591757.91%
ν0.2640.2816.05%
σbc(MPa)263.4502505.38%
σbb(MPa)50.475569.87%
KIC(MPam1/2)1.9801.6007.61%
BCE(GPa)188.2102005.90%
ν0.2780.3007.33%
σb(MPa)413.5994160.58%
SubstrateE(GPa)237.2072207.82%
ν0.3010.3102.90%
σb(MPa)1212.20011505.41%
Hierarchical structureMacroscopic parameterSimulated valueExperimental value [2326]Error
TCE(GPa)161.1591757.91%
ν0.2640.2816.05%
σbc(MPa)263.4502505.38%
σbb(MPa)50.475569.87%
KIC(MPam1/2)1.9801.6007.61%
BCE(GPa)188.2102005.90%
ν0.2780.3007.33%
σb(MPa)413.5994160.58%
SubstrateE(GPa)237.2072207.82%
ν0.3010.3102.90%
σb(MPa)1212.20011505.41%

2.2 The Establishment of Discrete Element Method Model.

To systematically investigate the effect of coating thickness on the mechanical properties and damage behavior of the TBCs, two different coating thicknesses were designed in this study. The first group (G1) has a TC layer thickness of 300 µm, a BC layer thickness of 150 µm, and a substrate thickness of 750 µm; the second group (G2) has a TC layer thickness of 600 µm, a BC layer thickness of 150 µm, and a substrate thickness of 750 µm. The length of both sets of samples is 12 mm. The Particle Flow Code (PFC) software was used to build DEM models for simulating the 4PB test. The model consists of a collection of particles bonded together, with the particle positions randomly distributed in an irregular pattern. The G1 model contains 83,587 particles, while the G2 model contains 93,890 particles. In the simulation, the upper end of the model is fixed by a support disk, while the lower end is moved upward at a speed of 16 µm/s by a loading disk to replicate the 4PB test process. The schematic of the loading and support arrangement of the model is shown in Fig. 2.

Fig. 2
Numerical simulation of 4PB test
Fig. 2
Numerical simulation of 4PB test
Close modal

2.3 Mechanisms of Crack Formation and Damage in Thermal Barrier Coatings.

Figure 3 shows the crack propagation and fracture process of G1 during the 4PB test. Multiple cracks formed during the bending process, with the main crack located near the loading disk. The enlarged view in Fig. 3 further reveals the details of a typical crack. Initially, when the stress on the TC layer surface exceeds its fracture limit, transverse cracks perpendicular to the interface appear at the interface. Due to the sufficiently high fracture toughness of the BC layer, these cracks are blocked at the TC/BC interface (Fig. 3(a)), corresponding to point A(A) on the curve in Fig. 9.

Fig. 3
Damage stages of TBCs during 4PB: (a) transverse crack initiation, (b) crack deflection and microcracking, (c) crack propagation, (d) interface crack growth, (e) TC/BC delamination, and (f) BC/substrate delamination (The lines in the TC layer denote cracks.)
Fig. 3
Damage stages of TBCs during 4PB: (a) transverse crack initiation, (b) crack deflection and microcracking, (c) crack propagation, (d) interface crack growth, (e) TC/BC delamination, and (f) BC/substrate delamination (The lines in the TC layer denote cracks.)
Close modal

Then, as the bending continues, new transverse cracks initiate on the surface of the TC layer. At the same time, the deflection tendency of the transverse cracks gradually increases, leading to the formation of interface microcracks at the TC/BC interface (Fig. 3(b)), corresponding to point B(B) on the curve in Fig. 9. However, such transverse cracks do not propagate before reaching steady-state expansion conditions. In such a way, multiple transverse cracks rapidly propagate (some transverse cracks initiate from the TC/BC interface and extend upward to the surface of the TC layer), causing further damage (Fig. 3(c)), corresponding to point C(C) on the curve in Fig. 9.

When the transverse cracks in the TC layer approach saturation, interface cracks at the TC/BC interface begin to expand rapidly. Additionally, transverse cracks initiate from this interface and propagate through the BC layer, eventually being blocked at the BC/substrate interface (Fig. 3(d)), corresponding to point D(D) on the curve in Fig. 9. Then, steady-state expansion of the interface cracks at the TC/BC interface leads to the convergence of these cracks, resulting in delamination at the TC/BC interface (Fig. 3(e)), corresponding to point E(E) on the curve in Fig. 9. Finally, when the transverse cracks in the BC layer reach the BC/substrate interface, the cracks deflect into this interface and begin propagating along it, eventually causing delamination between the BC layer and the substrate (Fig. 3(f)), corresponding to point F(F) on the curve in Fig. 9.

For G2, the cracking process is similar to that of G1. However, compared to G1, G2 cracks at lower deflection values and strain rates (Fig. 11). By comparing the fracture conditions of two coating thicknesses under the same strain conditions (Fig. 4), it can be observed that the thick coating exhibits a wider TC/BC interface delamination zone and a narrower BC/substrate interface delamination zone compared to the thin coating. Furthermore, the results show that as the TC layer thickness increases, the delamination location shifts toward the TC/BC interface, which is similar to the experimental results of Jiang et al. [27].

Fig. 4
Delamination behavior of TBCs under 0.6% strain: (a) G1 and (b) G2 (cracks were concealed)
Fig. 4
Delamination behavior of TBCs under 0.6% strain: (a) G1 and (b) G2 (cracks were concealed)
Close modal
The microstructure of the TBCs produced by atmospheric plasma spraying exhibits various microstructural features, including pores, lamellar boundaries, oxide inclusions, and different phases. These may lead to coating failure [28]. However, the quantity and spatial distribution of these pores exhibit randomness. Due to its applicability to brittle materials such as ceramics and geotechnical materials, the Weibull distribution is widely used to describe the size and strength distribution characteristics of material defects [29,30]. In the numerical simulation of the TC layer, this study employs the Weibull distribution to statistically model the size and distribution characteristics of large pores and, in conjunction with Eqs. (12) and (13), constructs a coating model containing pore defects. Additionally, the study investigates the influence of pores on the crack propagation behavior of the coating.
(12)
(13)
where xi and yi represent the corresponding pore's x and y coordinates, Rxi and Ryi are random numbers uniformly distributed between 0 and 1, and β and α are the size and shape parameters, respectively.

As previously mentioned, when there are no pore defects in the coating, cracks propagate along a direction perpendicular to the induced maximum tensile stress (Fig. 5(a)). However, due to the presence of microcracks and pores ahead of the main crack, in the TC layer (Fig. 5(b)), cracks nucleate at the valley regions of the top coat surface and around the preexisting pores. When the crack extends to an existing pore, phenomena such as crack penetration through the pore, crack path deflection, and crack bifurcation (indicated by the arrows) occur. The crack propagation behavior is consistent with the experimental observations reported in Ref. [31] (Fig. 5(c)). The crack branching mechanism is a complex phenomenon, influenced by various factors, including the crack tip velocity, stress intensity factor, and their rates [32].

Fig. 5
Comparison of crack propagation in TBCs: (a) defect-free, (b) porous defect, and (c) porous defect from Ref. [31]
Fig. 5
Comparison of crack propagation in TBCs: (a) defect-free, (b) porous defect, and (c) porous defect from Ref. [31]
Close modal

2.4 Thermal Barrier Coatings Force Chain Distribution Characteristics Analysis.

During the 4PB test, the overall particle contact state of the TBCs transitions from tensile stress in the center region of the coating to compressive stress in the shoulder region (near the end constraints). The force chain distribution characteristics of G1 are shown in Fig. 6. Initially, tensile contact primarily occurs on the upper surface of the TC layer and the lower surface of the substrate, while the shoulder region and BC layer are mainly in a compressed state (Fig. 6(a)).

Fig. 6
Evolution of particle contact states in TBCs during the 4PB test: (a) initial contact state, (b) initiation of transverse crack in the TC layer, (c) transverse crack in the TC layer approaching saturation, (d) initiation of transverse crack in the BC layer, (e) delamination at the TC/BC interface, and (f) delamination at the BC/substrate interface
Fig. 6
Evolution of particle contact states in TBCs during the 4PB test: (a) initial contact state, (b) initiation of transverse crack in the TC layer, (c) transverse crack in the TC layer approaching saturation, (d) initiation of transverse crack in the BC layer, (e) delamination at the TC/BC interface, and (f) delamination at the BC/substrate interface
Close modal

When the transverse cracks in the TC layer nucleate, the contact state at the crack changes to compression, and tensile contact moves outward from the crack, concentrating near the BC/substrate interface (Fig. 6(b)), corresponding to Fig. 3(a). As the transverse cracks in the TC layer approach saturation, the contact state in the BC layer changes to tensile, with tensile contact primarily concentrated in the BC layer (Fig. 6(c)), corresponding to Figs. 3(b) and 3(c). When transverse cracks nucleate in the BC layer, tensile contact in both the TC and BC layers decreases, with tensile contact concentrated near the BC/substrate interface (Fig. 6(d)), corresponding to Fig. 3(d). When delamination occurs at the TC/BC interface, the contact state of the TC layer changes to compression, and tensile contact in the BC layer decreases (Fig. 6(e)), corresponding to Fig. 3(e). When delamination occurs at the BC/substrate interface, a significant amount of tensile contact appears in the fractured delaminated part of the TC layer, while tensile contact is primarily concentrated near the substrate at the BC/substrate interface (Fig. 6(f)), corresponding to Fig. 3(f).

Additionally, it is worth noting that the study found that before crack initiation, the particles at the crack initiation location were all in tensile contact, and crack formation was primarily caused by tensile instability fracture. Tensile failure dominates the damage mechanism. At the same time, the shoulder region (end constraints) remains in a compressed contact state, with no significant damage observed. Therefore, it can be inferred that the coating exhibits strong damage resistance under compressive conditions. For G2, similar evolutionary behavior and the same damage pattern are observed under the 4PB load. This is because, in both coatings, the damage mechanism is dominated by tensile failure.

3 Theoretical Analysis

3.1 Strain and Stress at the Coating Surface Cracking.

The cracks in TBCs largely depend on the instantaneous stresses induced by strain. The stresses in TBCs include stresses caused by external loads and initial residual stresses. The former is induced by thermal or mechanical loads. The latter arises from the temperature gradient and deformation differences between the coating and substrate during the cooling process from the molten state. The deflection of the TBCs can be recorded (in this simulation, the maximum vertical displacement of particles between the two loading disks is recorded), and then, the corresponding strain values can be evaluated using Eqs. (14)(16) based on the geometric two-layer relations [33].
(14)
where ε is the strain corresponding to the deflection δ, h0 is the distance between the top/bond coating interface and the neutral axis of the composite beam, and ρ is the radius of curvature of the neutral axis.
(15)
where Es and Ec are the elastic moduli of the substrate and TC, respectively, hs and hc are the thickness of the substrate and the ceramic layer, respectively. It is noteworthy that the BC layer has similar mechanical properties to the substrate, and therefore can be simplified as a single layer with the substrate [34]. In this study, as shown in Table 4, Es = 161.159 GPa and Ec = 237.207 GPa.
(16)
where L is the difference between the outer and inner spans in the 4PB load.
The stress σ of elastic deformation of TBCs under bending can be expressed as:
(17)
where E* is the equivalent elastic modulus and σr is the initial residual stress. The distance between the coating and the neutral axis of the composite beam is
(18)
The equivalent elastic modulus E* is calculated using the following equation:
(19)
The initial residual stress σr is calculated using the following equation [35]:
(20)
(21)
where R represents the initial curvature of the TBCs during the cooling process, caused by residual stress due to thermal expansion mismatch. v is the Poisson's number, and the indices “c” and “s” refer to coating and substrate, respectively.

Figure 7 presents the stress–strain curve for the elastic deformation phase in this simulation. For G1, the yield strain during the elastic deformation phase is 0.05%, the yield strength is 187 MPa, and the bending modulus is 0.37 MPa. In contrast, G2 exhibits a relatively smaller yield strain of 0.04% during the elastic deformation phase, with a corresponding yield strength of 185 MPa and a bending modulus of 0.46 MPa. The results indicate that the thicker coating has a higher bending modulus, demonstrating greater resistance to deformation, whereas the thinner coating may experience a substrate effect, where part of the stress is borne by the substrate, leading to relatively lower stress under the same strain conditions.

Fig. 7
Stress–strain curve of the TBCs during the elastic deformation stage
Fig. 7
Stress–strain curve of the TBCs during the elastic deformation stage
Close modal

3.2 Mechanical Damage Mechanics Analysis of the Bending Test.

Figure 8 is a schematic diagram of the load–displacement analysis of the TBCs under bending load. As shown in the figure, when the load exceeds the elastic limit load Fe of the coating (where cracks begin to form in the TBCs), the curve deflects, and the tangent slope starts to decrease. If the coating remains undamaged after the load exceeds the elastic limit load, the displacement corresponding to load F1 is OH (point N in the figure). When damage occurs, the actual displacement u1 is OJ (point S in the figure), so the damage displacement ud1 is HJ, i.e.,
(22)
where ke is the slope of the elastic (linear) step.
Fig. 8
Load–displacement analysis schematic diagram
Fig. 8
Load–displacement analysis schematic diagram
Close modal
After the load reaches point S and continues to increase ΔF, the corresponding actual displacement u2=u1+Δu is OK (see point T in the load–displacement curve). When there is no damage, the displacement is OI, corresponding to point Q in the figure. Therefore, the damage displacement at point T, ud2, is IK, i.e.,
(23)
Consequently, the damage displacement associated with ΔF increased by Δud=ud2ud1 [36]. From Eqs. (22) and (23), we have
(24)

From Eq. (24), it can be seen that the displacement change is decomposed into elastic response (ΔF/ke) and damage displacement (Δud), thereby quantitatively distinguishing between elastic deformation and nonlinear damage effects. By calculating the variation of damage displacement with load, the accumulation pattern of material damage, the evolution trend of nonlinear behavior, and the characteristics of stiffness degradation and energy dissipation can be revealed.

4 Results and Discussion

4.1 Load and Displacement.

Figure 9 shows the normalized load–displacement curves obtained from the G1 and G2 models. To facilitate comparison with the data from Ref. [37] and make the curves more intuitive, the load–displacement curves were normalized according to Eqs. (25) and (26). All the curves in the figure show an initial linear range, which then changed to a nonlinear portion, indicative of an elastic-plastic response. In the OA(OA) stage, as the load applied by the loading disk increases, the TBCs exhibits linear deformation, and the accumulated strain energy provides the energy required for subsequent crack initiation. When the load reaches point A(A), it reaches the critical value, and the initial nonlinear “toe” region can be observed. In AB(AB) stage, the load decreases with increasing displacement. This may be due to cracks in the TC layer being saturated and then blocked at the TC/BC interface, resulting in the instantaneous release of accumulated strain energy [38]. Subsequently, in BE(BE) stage, as the loading displacement further increases, the load continues to increase. Finally, in EF(EF) stage, delamination occurs at the TC/BC interface and the BC/substrate interface, leading to the release of strain energy again. As displacement increases, the load growth slows down.

At the same time, Fig. 9 shows that the curve trends for samples with different coating thicknesses are generally consistent, but at the same displacement, the load for the thin coating is lower than that for the thick coating. This may be because the thick coating has a larger equivalent bending stiffness. Additionally, by comparing the numerical simulation results with the 4PB experimental results for yttria-stabilized zirconia (YSZ) coatings from Ref. [37], it is found that the curve trends are consistent, and the normalized results for load and displacement parameters roughly match. This validates that the DEM model can qualitatively investigate the relationship between load and displacement in the TBCs during the bending process
(25)
(26)
where Fnormal is the normalized load on the loading disk, F is the load applied to the loading disk, and Fmax is the maximum load measured during the experiment. unormal is the normalized displacement of the loading disk, u is the displacement of the loading disk, and umax is the maximum displacement during the loading disk test.

4.2 Crack Density.

Figure 10 presents the fracture modes of the TBCs under 4PB loading, and the fracture characteristics are consistent with those in Ref. [39] (Fig. 10(a)). Under a strain rate of 0.6%, the crack density of the thin coating G1 (Fig. 10(b)) is significantly higher than that of the thick coating G2 (Fig. 10(c)). Additionally, it can be observed that the spacing between transverse cracks is approximately one–two times the coating thickness. These values align well with the shear lag model used to determine crack spacing in thin films on ductile substrates [40]. According to this model, the typical crack spacing lc should satisfy 3<(σYSlc)/(σbChc)<23, where hc is the coating thickness, σYS is the yield strength of the substrate, and σbC is the tensile strength of the coating. Using σYS=718.5MPa from Ref. [41] and σbC=413.6MPa as calibrated in Table 4 (with the tensile strength of the TC layer being much smaller than that of the BC layer, it is neglected), the relationship 1<lc/hc<2 is obtained, which is in good agreement with the observed values from this simulation.

Fig. 9
Normalized load–vertical displacement curve
Fig. 9
Normalized load–vertical displacement curve
Close modal
Fig. 10
Cross-sectional crack images in the TBCs under 4PB tests: (a) Ref. [39], (b) G1, and (c) G2
Fig. 10
Cross-sectional crack images in the TBCs under 4PB tests: (a) Ref. [39], (b) G1, and (c) G2
Close modal

Figure 11 shows the relationship between the crack density and strain in the TBCs, where crack density is defined as the ratio of the number of transverse cracks to their maximum distribution distance. From the figure, it can be seen that the crack density in G1 is generally higher than that in G2, but during the initial stage of crack nucleation, the crack density in G2 is slightly higher than in G1. The growth trend of crack density is in good agreement with the 4PB results for YSZ TBCs presented in Ref. [42], and the values are similar. Moreover, for cracks in the BC layer, G2 produces cracks earlier than G1 at lower strain rates. Combining the comparative analysis in Fig. 3, it is found that once the strain reaches 0.5%, further increasing the strain does not significantly increase the number of cracks in the crack density but instead more easily leads to delamination of the coating. This result is consistent with the experimental observations in Ref. [35], indicating that delamination of the coating is the primary failure mode at high strain stages.

Fig. 11
Crack density-strain curve of the TBCs under 4PB tests (the arrow indicates the corresponding point in Fig. 3)
Fig. 11
Crack density-strain curve of the TBCs under 4PB tests (the arrow indicates the corresponding point in Fig. 3)
Close modal

4.3 Variation of Strain Energy and Acoustic Emission Counts.

The AE characteristics of TBCs are related closely to the formation of microcracks. In PFC, the formation of a microcrack is accompanied by a release of strain energy, which causes an AE event. The AE characteristics can be simulated from the statistics of microcrack numbers over a certain interval during loading [43]. Figure 12 shows the relationship between strain energy and AE counts as a function of strain in the TBCs. The results indicate that as the bending strain of the specimen increases, the strain energy of the material continuously accumulates, with part of the strain energy being converted into the energy required for crack propagation. Crack propagation is accompanied by local damage, leading to the formation of crack surfaces and changes in the internal stress state of the material, thereby affecting the release and transfer of strain energy. During this process, the rate of crack propagation is closely related to the applied load and the accumulation rate of the material's strain energy.

Fig. 12
Graph of strain energy and AE counts variation
Fig. 12
Graph of strain energy and AE counts variation
Close modal

Furthermore, the experiment further revealed the significant impact of coating thickness on crack formation and strain energy variation. The thick coating G2, due to its larger volume and complex structural characteristics, exhibits significantly higher strain energy and AE counts than the thin coating G1 under the same strain conditions. This is primarily because the thick coating can withstand larger stress concentration areas during loading, thereby accumulating more elastic energy and generating more crack sources.

The energy release rate is an important parameter for measuring the energy required for crack propagation. In the 4PB experiment, the strain energy release rate is closely related to the crack propagation rate. When the energy release rate exceeds the critical energy required for crack propagation, the crack will rapidly propagate. Otherwise, the crack will stagnate. The critical strain energy release rate (Gss) can be calculated using Eqs. (27)(31) [44]. The calculated critical strain energy release rate for G1 is 129.6 J/m2, while for G2, it is 338.8 J/m2. The thick coating G2 has a higher critical strain energy release rate and releases more energy during crack propagation, while still retaining a higher strain energy storage, indicating its stronger energy-bearing capacity. This further indicates that the thick coating exhibits superior damage resistance under loading conditions.
(27)
(28)
(29)
(30)
(31)
where P is the load when fracture occurs, b represents the width of the specimen (set as 1 in this study). I is the inertia moment, and the indices “c” and “s” refer to coating and substrate, respectively.

5 Conclusion

Through simulation experiments and theoretical analysis, the damage evolution behavior and mechanical properties of TBCs with two different thicknesses under 4PB tests were clarified. The study confirmed that the DEM model can effectively simulate the mechanical response of the TBCs, providing a reliable tool for coating performance prediction and design optimization. The main conclusions are summarized as follows:

  1. Under continuous bending, the damage process exhibits six stages: The transverse cracks first initiate and propagate in the TC layer. The cracks then multiply and become saturated; local interface cracks propagate between the TC and BC layers. Transverse cracks initiate and propagate in the BC layer. Local interface cracks propagate between the BC layer and the substrate; Finally, when some adjacent interface cracks connect with each other, interface delamination occurs.

  2. During the bending process, the overall particle contact state of the TBCs transitions from the tensile region at the center of the coating to the compressive region at the shoulder area (near the end constraints). Both coating thicknesses experience tensile instability fracture during the bending process, with tensile failure being the dominant mechanism of coating damage.

  3. This study systematically analyzes the relationship between load and displacement variation based on experimental load–displacement curves, exploring the underlying mechanisms and the relationship between different loading stages and the damage evolution of TBCs. The study also indicates that coating thickness has a significant impact on mechanical properties: thicker coatings have higher load-bearing capacity, but they also carry a greater risk of crack accumulation.

  4. By observing the fracture modes of the TBCs under 4PB, it was found that the crack density in the thin coating is significantly higher than that in the thick coating, with the spacing between transverse cracks being approximately one–two times the coating thickness. Furthermore, the thick coating generates cracks earlier than the thin coating at lower strain rates. The average crack spacing in the thick coating is larger than that in the thin coating, and its damage evolution rate is faster, exhibiting more pronounced interface delamination. Based on the comparative analysis of damage evolution, it can be observed that when the strain reaches 0.5%, further increasing the strain does not significantly increase the crack density but rather leads to easier delamination of the coating. These results indicate that coating delamination is the primary failure mode at high strain stages.

  5. By analyzing the variations of strain energy and AE counts with strain in the TBCs, it was found that as the bending strain increases, the strain energy of the material continuously accumulates and is partially converted into crack propagation energy. Crack propagation is accompanied by local damage, which alters the internal stress state and affects the release of strain energy. Compared to the thin coating G1, the thick coating G2 exhibits higher strain energy and AE counts under the same strain conditions. The strain energy release rate is closely related to the crack propagation rate, and when the release rate exceeds the critical value, cracks propagate rapidly. The thick coating G2 has a higher critical release rate, allowing it to release more energy during propagation while maintaining higher strain energy storage, demonstrating stronger energy-bearing capacity.

This study employs the DEM to analyze, at the microscale, the damage evolution mechanism of TBCs under 4PB conditions, providing new theoretical insights into the failure modes of TBCs. The study reveals the significant impact of coating thickness on damage behavior and mechanical properties, particularly in terms of differences in crack propagation and interface delamination damage evolution. This research provides a theoretical basis for optimizing coating design, enhancing coating load-bearing capacity, and predicting coating failure, offering important academic and engineering significance for the application and performance enhancement of TBCs.

Acknowledgment

Yafeng Li greatly appreciates the support from the Natural Science Foundation of Tianjin (21JCYBJC01400) and the China Scholarship Council (CSC) Scholarship.

Conflict of Interest

There are no conflicts of interest.

Data Availability Statement

The datasets generated and supporting the findings of this article are obtainable from the corresponding author upon reasonable request.

References

1.
Hegab
,
A.
,
Dahuwa
,
K.
,
Islam
,
R.
,
Cairns
,
A.
,
Khurana
,
A.
,
Shrestha
,
S.
, and
Francis
,
R.
,
2021
, “
Plasma Electrolytic Oxidation Thermal Barrier Coating for Reduced Heat Losses in IC Engines
,”
Appl. Therm. Eng.
,
196
, p.
117316
.
2.
Qiao
,
D.
,
Yan
,
W. G.
,
Zeng
,
W.
,
Man
,
J. X.
,
Xue
,
B. R.
, and
Bian
,
X. D.
,
2024
, “
Mechanical Loading Effect on Stress States and Failure Behavior in Thermal Barrier Coatings
,”
Crystals
,
14
(
1
), p.
2
.
3.
Jiang
,
J. S.
,
Zou
,
Z. H.
,
Wang
,
W. Z.
,
Zhao
,
X. F.
,
Liu
,
Y. Z.
, and
Cao
,
Z. M.
,
2018
, “
Effect of Internal Oxidation on the Interfacial Morphology and Residual Stress in Air Plasma Sprayed Thermal Barrier Coatings
,”
Surf. Coat. Technol.
,
334
, pp.
215
226
.
4.
Yu
,
Q. M.
, and
He
,
Q.
,
2018
, “
Effect of Material Properties on Residual Stress Distribution in Thermal Barrier Coatings
,”
Ceram. Int.
,
44
(
3
), pp.
3371
3380
.
5.
Bhattacharyya
,
A.
, and
Maurice
,
D.
,
2019
, “
Residual Stresses in Functionally Graded Thermal Barrier Coatings
,”
Mech. Mater.
,
129
, pp.
50
56
.
6.
Yu
,
Q. M.
, and
Cen
,
L.
,
2017
, “
Residual Stress Distribution Along Interfaces in Thermal Barrier Coating System Under Thermal Cycles
,”
Ceram. Int.
,
43
(
3
), pp.
3089
3100
.
7.
Li
,
Y. F.
,
Gao
,
H. F.
,
Wang
,
L.
,
Sun
,
Y. L.
, and
Zhang
,
J.
,
2024
, “
Discrete Element Modeling of Effect of Interfacial Roughness and Pre-Crack on Crack Propagation in Thermal Barrier Coatings
,”
Comput. Part. Mech.
,
12
(
1
) pp.
1
9
.
8.
Zhang
,
X. F.
,
Zhou
,
K. S.
,
Liu
,
M.
,
Deng
,
C. M.
,
Deng
,
C. G.
,
Mao
,
J.
, and
Deng
,
Z. Q.
,
2018
, “
Mechanisms Governing the Thermal Shock and Tensile Fracture of PS-PVD 7YSZ TBC
,”
Ceram. Int.
,
44
(
4
), pp.
3973
3980
.
9.
Jing
,
F. L.
,
Yang
,
J. J.
,
Tang
,
S. B.
,
Wen
,
Q.
,
Zhang
,
T.
,
Wu
,
J.
, and
Fan
,
X. L.
,
2022
, “
Quantitative Characterization of the Interfacial Damage in EB-PVD Thermal Barrier Coating
,”
Coatings
,
12
(
7
), p.
984
.
10.
Sulák
,
I.
,
Obrtlík
,
K.
,
Celko
,
L.
,
Chráska
,
T.
,
Jech
,
D.
, and
Gejdos
,
P.
,
2018
, “
Low Cycle Fatigue Performance of Ni-Based Superalloy Coated With Complex Thermal Barrier Coating
,”
Mater. Charact.
,
139
, pp.
347
354
.
11.
Zhao
,
P. F.
,
Sun
,
C. A.
,
Zhu
,
X. Y.
,
Shang
,
F. L.
, and
Li
,
C. J.
,
2010
, “
Fracture Toughness Measurements of Plasma-Sprayed Thermal Barrier Coatings Using a Modified Four-Point Bending Method
,”
Surf. Coat. Technol.
,
204
(
24
), pp.
4066
4074
.
12.
Zhao
,
P. F.
,
Li
,
X. D.
,
Shang
,
F. L.
, and
Li
,
C. J.
,
2011
, “
Interlamellar Cracking of Thermal Barrier Coatings With TGOs by Non-Standard Four-Point Bending Tests
,”
Mater. Sci. Eng. A
,
528
(
25–26
), pp.
7641
7647
.
13.
Chen
,
Y.
,
Zhang
,
X.
,
Zhao
,
X. F.
,
Markocsan
,
N.
,
Nylén
,
P.
, and
Xiao
,
P.
,
2019
, “
Measurements of Elastic Modulus and Fracture Toughness of an Air Plasma Sprayed Thermal Barrier Coating Using Micro-Cantilever Bending
,”
Surf. Coat. Technol.
,
374
, pp.
12
20
.
14.
Chen
,
Z. X.
,
Qian
,
L. H.
, and
Zhu
,
S. J.
,
2010
, “
Determination and Analysis of Crack Growth Resistance in Plasma-Sprayed Thermal Barrier Coatings
,”
Eng. Fract. Mech.
,
77
(
11
), pp.
2136
2144
.
15.
Jing
,
J. N.
,
Dong
,
L. H.
,
Wang
,
H. D.
, and
Jin
,
G.
,
2018
, “
Application of Acoustic Emission Technology for Quantitative Characterization of Plasma-Sprayed Coatings Subjected to Bending Fatigue Tests
,”
J. Therm. Spray Technol.
,
27
(
7
), pp.
1090
1102
.
16.
Long
,
H.
,
Liang
,
L. H.
, and
Wei
,
Y. G.
,
2019
, “
Failure Characterization of Solid Structures Based on an Equivalence of Cohesive Zone Model
,”
Int. J. Solids Struct.
,
163
, pp.
194
210
.
17.
Wang
,
L.
,
Ni
,
J. X.
,
Shao
,
F.
,
Yang
,
J. S.
,
Zhong
,
X. H.
,
Zhao
,
H. Y.
,
Liu
,
C. G.
,
Tao
,
S. Y.
,
Wang
,
Y.
, and
Li
,
D. Y.
,
2017
, “
Failure Behavior of Plasma-Sprayed Yttria-Stabilized Zirconia Thermal Barrier Coatings Under Three-Point Bending Test via Acoustic Emission Technique
,”
J. Therm. Spray Technol.
,
26
(
1–2
), pp.
116
131
.
18.
Eriksson
,
R.
,
Chen
,
Z.
, and
Jonnalagadda
,
K. P.
,
2017
, “
Bending Fatigue of Thermal Barrier Coatings
,”
ASME J. Eng. Gas Turbines Power
,
139
(
12
), p.
122101
.
19.
Frommherz
,
M.
,
Scholz
,
A.
,
Oechsner
,
M.
,
Bakan
,
E.
, and
Vaßen
,
R.
,
2016
, “
Gadolinium Zirconate/YSZ Thermal Barrier Coatings: Mixed-Mode Interfacial Fracture Toughness and Sintering Behavior
,”
Surf. Coat. Technol.
,
286
, pp.
119
128
.
20.
Martins
,
J. P.
,
Yu
,
H. Y.
,
Chen
,
Y.
,
Brewster
,
G.
,
McIntyre
,
R.
, and
Xiao
,
P.
,
2021
, “
Effect of Bond Coat Topography on the Fracture Mechanics and Lifetime of Air-Plasma Sprayed Thermal Barrier Coatings
,”
Surf. Coat. Technol.
,
421
, p.
127447
.
21.
Vaßen
,
R.
,
Bakan
,
E.
, and
Schwartz-Lückge
,
S.
,
2021
, “
Influence of Substrate Removal Method on the Properties of Free-Standing YSZ Coatings
,”
Coatings
,
11
(
4
), p.
449
.
22.
Wang
,
X.
, and
Tian
,
L. G.
,
2018
, “
Mechanical and Crack Evolution Characteristics of Coal-Rock Under Different Fracture-Hole Conditions: a Numerical Study Based on Particle Flow Code
,”
Environ. Earth Sci.
,
77
(
8
), pp.
1
10
.
23.
Ji
,
Z. Z.
,
Li
,
Y. F.
,
Wang
,
L.
,
Dube
,
T.
, and
Zhang
,
J.
,
2023
, “
Discrete Element Modeling of Delamination Behavior in Thermal Barrier Coating
,”
Comput. Part. Mech.
,
10
(
5
), pp.
1019
1029
.
24.
Wang
,
L.
,
Wang
,
Y.
,
Sun
,
X. G.
,
Pan
,
Z. Y.
,
He
,
J. Q.
, and
Li
,
C. G.
,
2011
, “
Influence of Pores on the Surface Microcompression Mechanical Response of Thermal Barrier Coatings Fabricated by Atmospheric Plasma Spray-Finite Element Simulation
,”
Appl. Surf. Sci.
,
257
(
6
), pp.
2238
2249
.
25.
Zhou
,
Y. C.
, and
Hashida
,
T.
,
2001
, “
Coupled Effects of Temperature Gradient and Oxidation on Thermal Stress in Thermal Barrier Coating System
,”
Int. J. Solids Struct.
,
38
(
24–25
), pp.
4235
4264
.
26.
Qi
,
H. Y.
,
Yang
,
X. G.
, and
Wang
,
Y. M.
,
2009
, “
Interfacial Fracture Toughness of APS Bond Coat/Substrate Under High Temperature
,”
Int. J. Fract.
,
157
(
1–2
), pp.
71
80
.
27.
Jiang
,
P.
,
Fan
,
X. L.
,
Sun
,
Y. L.
,
Li
,
D. J.
,
Li
,
B.
, and
Wang
,
T. J.
,
2017
, “
Competition Mechanism of Interfacial Cracks in Thermal Barrier Coating System
,”
Mater. Des.
,
132
, pp.
559
566
.
28.
Di Girolamo
,
G.
,
Marra
,
F.
,
Schioppa
,
M.
,
Blasi
,
C.
,
Pulci
,
G.
, and
Valente
,
T.
,
2015
, “
Evolution of Microstructural and Mechanical Properties of Lanthanum Zirconate Thermal Barrier Coatings at High Temperature
,”
Surf. Coat. Technol.
,
268
, pp.
298
302
.
29.
Nakamura
,
S.
,
Tanaka
,
S.
,
Furushima
,
R.
,
Sato
,
K.
, and
Uematsu
,
K.
,
2009
, “
Estimation of Weibull Modulus From Coarser Defect Distribution in dry-Pressed Alumina Ceramics
,”
J. Ceram. Soc. Jpn.
,
117
(
1366
), pp.
742
747
.
30.
Zhang
,
Z. H.
,
Kameda
,
J.
,
Sakurai
,
S.
, and
Sato
,
M.
,
2005
, “
Through-Thickness Dependence of In-Plane Cracking Behavior in Plasma-Sprayed Thermal Barrier Coatings
,”
Metall. Mater. Trans. A
,
36A
(
7
), pp.
1841
1854
.
31.
Amer
,
M.
,
Curry
,
N.
,
Arshad
,
M.
,
Hayat
,
Q.
,
Janik
,
V.
,
Nottingham
,
J.
, and
Bai
,
M. W.
,
2023
, “
Unraveling the Cracking Mechanisms of Air Plasma-Sprayed Thermal Barrier Coatings: An In-Situ SEM Investigation
,”
Coatings
,
13
(
9
), p.
1493
.
32.
Sun
,
Y. N.
,
Edwards
,
M. G.
,
Chen
,
B.
, and
Li
,
C. F.
,
2021
, “
A State-of-the-Art Review of Crack Branching
,”
Eng. Fract. Mech.
,
257
, p.
108036
.
33.
Chiu
,
C.-C.
, and
Case
,
E. D.
,
1991
, “
Elastic Modulus Determination of Coating Layers as Applied to Layered Ceramic Composites
,”
Mater. Sci. Eng. A
,
132
, pp.
39
47
.
34.
Ray
,
A. K.
, and
Steinbrech
,
R. W.
,
1999
, “
Crack Propagation Studies of Thermal Barrier Coatings Under Bending
,”
J. Eur. Ceram. Soc.
,
19
(
12
), pp.
2097
2109
.
35.
Andritschky
,
M.
,
Alpuim
,
P.
,
Stöver
,
D.
, and
Funke
,
C.
,
1999
, “
Study of the Mechanics of the Delamination of Ceramic Functional Coatings
,”
Mater. Sci. Eng. A
,
271
(
1–2
), pp.
62
69
.
36.
Liu
,
H. Y.
,
Wei
,
Y. G.
,
Liang
,
L. H.
,
Liu
,
X. H.
,
Wang
,
Y. B.
, and
Ma
,
H. S.
,
2018
, “
Damage Characterization Model of Ceramic Coating Systems Based on Energy Analysis and Bending Tests
,”
Ceram. Int.
,
44
(
5
), pp.
4807
4813
.
37.
Kucuk
,
A.
,
Berndt
,
C. C.
,
Senturk
,
U.
,
Lima
,
R. S.
, and
Lima
,
C. R. C.
,
2000
, “
Influence of Plasma Spray Parameters on Mechanical Properties of Yttria Stabilized Zirconia Coatings. I: Four Point Bend Test
,”
Mater. Sci. Eng. A
,
284
(
1–2
), pp.
29
40
.
38.
Zhu
,
Q.
,
He
,
W.
,
Zhu
,
J. G.
,
Zhou
,
Y. X.
, and
Chen
,
L.
,
2018
, “
Investigation on Interfacial Fracture Toughness of Plasma-Sprayed TBCs Using a Three-Point Bending Method
,”
Surf. Coat. Technol.
,
353
, pp.
75
83
.
39.
Jiang
,
P.
,
Fan
,
X. L.
,
Sun
,
Y. L.
,
Li
,
D. J.
, and
Wang
,
T. J.
,
2018
, “
Bending-Driven Failure Mechanism and Modelling of Double-Ceramic-Layer Thermal Barrier Coating System
,”
Int. J. Solids Struct.
,
130
, pp.
11
20
.
40.
Hu
,
M.
, and
Evans
,
A. J. A. M.
,
1989
, “
The Cracking and Decohesion of Thin Films on Ductile Substrates
,”
Acta Metall.
,
37
(
3
), pp.
917
925
.
41.
Kong
,
Y. H.
,
Liu
,
R. Y.
,
Chen
,
G. S.
,
Xie
,
L. X.
, and
Zhu
,
S. G.
,
2012
, “
Effect of Different Solution Treatments on Room Temperature Tensile Properties of GH4169 Alloy
,”
Rare Met. Mater. Eng.
,
41
, pp.
616
620
.
42.
Zhang
,
X. C.
,
Watanabe
,
M.
, and
Kuroda
,
S.
,
2013
, “
Effects of Residual Stress on the Mechanical Properties of Plasma-Sprayed Thermal Barrier Coatings
,”
Eng. Fract. Mech.
,
110
, pp.
314
327
.
43.
Yoon
,
J. S.
,
Zang
,
A. N.
, and
Stephansson
,
O.
,
2014
, “
Numerical Investigation on Optimized Stimulation of Intact and Naturally Fractured Deep Geothermal Reservoirs Using Hydro-Mechanical Coupled Discrete Particles Joints Model
,”
Geothermics
,
52
, pp.
165
184
.
44.
Li
,
H.
,
Khor
,
K. A.
, and
Cheang
,
P.
,
2002
, “
Young's Modulus and Fracture Toughness Determination of High Velocity Oxy-Fuel-Sprayed Bioceramic Coatings
,”
Surf. Coat. Technol.
,
155
(
1
), pp.
21
32
.