## Mathematical Modeling and Optimization of Functionally Graded Structures

View this Special IssueResearch Article | Open Access

# Thermoelastic Characteristics in Thermal Barrier Coatings with a Graded Layer between the Top and Bond Coats

**Academic Editor:**Abdelouahed Tounsi

#### Abstract

A graded layer was introduced at the interface between the top and bond coats to reduce the risk of failure in a thermal barrier coating (TBC) system, and the thermoelastic behavior was investigated through mathematical approaches. Two types of TBC model with and without the graded layer, subject to a symmetric temperature distribution in the longitudinal direction, were taken into consideration to evaluate thermoelastic behaviors such as temperature distribution, displacement, and thermal stress. Thermoelastic theory was applied to derive two governing partial differential equations, and a finite volume method was developed to obtain approximations because of the complexity. The TBC with the graded layer shows improved durability in thermoelastic characteristics through mathematical approaches, in agreement with the experimental results. The results will be useful in discovering technologies for enhancing the thermomechanical properties of TBCs.

#### 1. Introduction

The deposition process for thermal barrier coatings (TBCs) exerts a critical influence on the determination of the thermomechanical properties, such as elastic modulus, thermal conductivity, and coefficient of thermal expansion. Moreover, failure phenomena usually occur at the interface between the top and bond coats because of the mismatch of mechanical and thermal properties as the TBC system cools from a high operating temperature to ambient temperature. To improve the thermal durability of TBCs for protecting hot components such as combustors and turbine parts, functionally graded materials (FGMs) and composite materials (CMs) are widely employed [1–9]. FGMs consist of two or more distinct material phases whose volume fractions continuously vary with space variables, thus yielding continuously varying mechanical and thermal properties in the FGM system. CMs are products composed of two or more distinct constituent materials involving different physical or chemical properties such that their constituents remain still distinguishable at the macroscopic or microscopic scale within the structure. The original purpose of these materials is, simultaneously, to resist high temperature and ensure low thermal conductivity through TBC techniques.

CMs make use of improved characteristics such as strength, fatigue life, stiffness, temperature-dependent behavior, thermal insulation, thermal conductivity, and weight. For example, TBCs are applied in building aircraft to obtain light and strong materials in surfacing while insulating against hazardous physical factors such as heat and lightning. In aircraft, aerospace, automotive, and shipping applications, fiber-reinforced and resin-matrix CMs contribute to the improvement through their higher strength and stiffness benefits [10–12]. However, owing to different macroscopic or microscopic scales within the structure, CMs are exposed to the risk of delamination during service. On the other hand, FGMs represent improved bonding strength, toughness, wear and corrosion resistance, and reduced residual and thermal stresses because of the graded material properties. These outstanding advantages of FGMs over conventional composites and monolithic materials have led to extensive study for potential applications as structural elements, such as FGM beams [13, 14], plates [15, 16], shells [17], and cylinders [18, 19], and have provoked the interest of researchers toward the development of potential new structural and functional applications. Therefore, FGMs and CMs are being explored to determine their characteristics from various points of view.

The microstructural characterizations of TBCs are greatly influenced by the size, shape, and density of the feedstock powder [20–23]. In addition, the deposition process parameter is an essential element in determining its microstructure. Cernuschi et al. [24] and Jung et al. [25] analyzed the microstructure features of TBCs manufactured using the air plasma spray process, and the effects of thermal annealing on the microstructure of TBCs deposited by electron beam physical vapor deposition have been investigated by Kulkarni et al. [26] and Schulz and Schmücker [27]. Moreover, the study of microstructural behavior has been extended to functionally graded (FG) TBC materials. The microstructure and phase stability of FG TBCs based on LaMgAl_{11}O_{19} (LaMa)/YSZ have been investigated by Chen et al. [28], and the results indicated that all of the LaMa and LaMa-containing intermediate composite coatings suffer irreversible phase transformations induced by the recrystallization of amorphous LaMa coating. Schulz et al. [29] tried to optimize the thermal, wear, and corrosion properties of FGMs and found that TBCs deposited onto Cu substrates by pulsed laser deposition showed improved spallation behavior by a graded lamella microstructure with improved interface fracture toughness. Shaw [30] studied the thermal residual stresses in multilayered and compositionally graded plates and coatings.

In the present paper, the thermoelastic characteristics of circular disk TBC specimens with and without an FG layer between the top and bond coats are analyzed through mathematical approaches. Two partial differential equations are derived based on thermoelastic theory, and the thermoelastic characteristics, such as temperature distribution profiles, displacement, and stresses, are determined through mathematical approaches. Because of the complexity of the governing equations, a finite volume approach is adopted to analyze the thermoelastic characteristics.

#### 2. Mathematical Modeling

##### 2.1. Thermoelastic Formulation

According to the assumption that the temperature profiles pressure the circular disk to the longitudinal direction only, all quantities are independent of circumferential -direction. Equilibrium equations in polar coordinates thus are Let be the displacement to the radial direction, and let be the displacement to the longitudinal direction. The strain-displacement relations are as follows: In addition, the strain-stress relations due to the action of thermal loading can be expressed as where .

The combination of (1)–(3) leads to the following governing equations:

where .

##### 2.2. Temperature Distribution Formulation

Since the circular disk is subject to a uniform temperature loading in the longitudinal -direction only (see Figure 1), the differential equation for the temperature distribution profile is The general solution of (5) is where and are integration constants for the layer temperature distribution. The integration constants will be determined based on the following boundary conditions: but only two boundary values and are known, and so heat flux at each layer point is considered. The heat flux at layer point is expressed as where is the heat flux into layer, is the conductivity, and is the length of layer. At each layer point, the integral constants for the temperature distribution profile are determined uniquely by solving the following linear system: The is a known value solved by layer linear system and is a given initial value, and thus the number of equations will determine ’s unknown coefficients. Here, is the number of layers. The temperature distribution profiles obtained through this process will be applied to (4a) and (4b) for the thermoelastic characteristics.

##### 2.3. FG Layer Formulation

The elastic modulus , Poisson’s ratio , thermal expansion coefficient , and thermal conductivity are assumed to vary exponentially with the variable ; that is,

The inner area of the circular disk consists of 100% of the top coat material, whereas the outer area of the disk has 100% of the bond coat material; that is, the properties in (10a), (10b), (10c), and (10d) have values equal to those of the materials in the top and bond coats, respectively, at the inner and outer areas of the disk. Thus, the constants in (4a) and (4b) can be determined as

The subscripts and on a variable represent the properties of the constituent materials in the top and bond coats, respectively. The mechanical and thermal properties obtained through this process will be used to analyze the thermoelastic characteristics of graded layer between the top coat and bond coat.

##### 2.4. Discretization of the Governing Equation Based on a Finite Volume Method

The governing equations for the microstructural characterization of TBCs are too complicated to be solved analytically, and a finite volume method is adopted for numerical approximations. The model is divided up into control volumes and integrates the field equations over each control volume (see Figure 2). In formulating a numerical approximation for the present problem, the finite surface mesh is denoted by , and the discretizations for the governing equations are developed based on the following relations at the adjacent locations: In the previous equations a subscript 1/2 implies the value of the displacement at the boundary of the control surface. According to the above relations at the adjacent locations, (4a) is discretized as follows:Similar process is adjusted to (4b) and arrives at

The finite volume model is developed based on the following boundary conditions:

#### 3. Results and Discussion

The temperature distribution profiles and thermoelastic characteristics were investigated based on the mechanical and thermal properties shown in Table 1 using the processes in Sections 2.2 and 2.3. The mechanical and thermal properties of the FG layer are obtained based on (11a), (11b), (11c), and (11d). The temperature distribution profiles of the TBCs with and without the FG layer are shown in Figure 3. At the top coat of the TBC without the FG layer (hereafter CM), the decrease in temperature is linear, while the TBC with the FG layer (hereafter FGM) exhibits an exponential decrease. Over the entire domain, the circle with the FG layer is under a loading from a lower temperature, which implies an improved durability for the FGM. The rate of temperature decrease becomes slower as the metal concentration increases along the longitudinal direction and demonstrates that the mathematical approach is reliable and reasonable.

| |||||||||||||||||||||||||||||||||||

*Properties of the top coat prepared using an air plasma spray with METCO 204 NS-G powder. |

The longitudinal and radial displacement profiles are presented in Figures 4 and 5, respectively. Unlike the displacement profile of the CM at the top and bond coats, the longitudinal displacement of the FGM fluctuates at the top coat, the FG layer, and the bond coat (see Figure 4). At the substrate, the two models exhibit similar behavior, although a lower displacement develops along the -direction in the FGM. The radial displacement profiles of each layer in the FGM are displayed in Figure 5(a). At the top coat, a larger fluctuation appears near the center, while the magnitude of the radial displacement at the FG layer is almost constant for the interval , which implies that the durability is increased through the FG layer. A small radial displacement develops at the bond coat because of the relaxation of delamination. The radial behaviors at the bond coat in both the CM and the FGM are compared in Figure 5(b). The magnitude of the radial displacement of the FGM is almost negligible, whereas the bond coat of the CM undulates with a large magnitude in comparison with the FGM.

**(a)**

**(b)**

Figure 6 represents the longitudinal and circumferential stresses. The tensile longitudinal stress appearing at the interface between the bond coat and the substrate of the CM converts into a compressive stress in the FGM, and the FGM experiences larger compressive longitudinal stress at the top coat, the FG layer, and the bond coat (see Figure 6(a)), indicating that the FGM will improve the thermal stability of the TBC system during service. A tensile longitudinal stress develops at the substrate in the FGM. In the circumferential stress, the magnitude is larger at the top coat, the FG layer, and the bond coat in the FGM; on the other hand, a smaller magnitude appeared at the substrate (see Figure 6(b)). Both models are influenced by the compressive circumferential stress over the entire domain.

**(a)**

**(b)**

The radial stress profiles are displayed in Figure 7. At the top coat, the largest undulation occurs near the center and the magnitude decreases as the radius increases (see Figure 7(a)). A fluctuation with the smallest magnitude develops at the bond coat because of the alleviation of stress at the FG layer. A comparison of the radial stress at the bond coat between the CM and the FGM is shown in Figure 7(b). For the interval , a large fluctuation is produced in the CM, while the influence of the radial stress is negligible in the FGM, showing the superiority of the FGM.

**(a)**

**(b)**

Overall, as compared with the models, the FG layer in TBCs provides buffer zone interactions between the top coat and the bond coat and absorbs critical shocks. This phenomenon is shown through the thermoelastic behavior at the bond coat of the FGM: (i) the magnitude of the radial displacement is trivial in comparison with that of the CM, (ii) the tensile longitudinal stress at the interface between the bond coat and the substrate in the CM converts into a compressive stress, and (iii) a negligible fluctuation is developed along the radius in comparison with that of the CM. The results demonstrate that the thermal durability of the TBC can be improved by introducing an FG layer between the top coat and the bond coat, and the mathematical analysis presented is a reasonable approach. The results and analysis performed here contribute further understanding to the behavior of TBCs.

#### 4. Conclusion

The thermoelastic behavior of the FGM (circular disk TBC with FG layer) was investigated and compared with that of the CM (circular disk TBC without FG layer). The FGM is subject to a lower temperature loading over the entire domain, which reduces the magnitude of the longitudinal displacement. The FG layer in the FGM converts the tensile longitudinal stress appearing at the interface between the bond coat and the substrate of the CM into a compressive stress and increases the magnitude of the compressive stress at the top coat, the FG layer, and the bond coat. Moreover, an insignificant radial stress is developed at the bond coat in the FGM, contributing to the improved performance of the FG layer. Consequently, the thermoelastic characteristics obtained for temperature, displacements, and stresses demonstrate that the FG layer will decrease the risk of failure in the TBC system and the thermal durability of the TBC system can be improved during the lifetime cycle.

#### Acknowledgments

This research was supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology (2010-0024178). This work was supported by a National Research Foundation of Korea (NRF) grant funded by the Korean Government (MEST) (2012-0009450) and by a grant from the Fundamental R&D Program for Core Technology of Materials funded by the Korean Ministry of Knowledge Economy (10041233).

#### References

- X. Cao, F. Jin, I. Jeon, and T. J. Lu, “Propagation of Love waves in a functionally graded piezoelectric material (FGPM) layered composite system,”
*International Journal of Solids and Structures*, vol. 46, no. 22-23, pp. 4123–4132, 2009. View at: Publisher Site | Google Scholar | Zentralblatt MATH - M. L. Pines and H. A. Bruck, “Pressureless sintering of particle-reinforced metal-ceramic composites for functionally graded materials: part I. Porosity reduction models,”
*Acta Materialia*, vol. 54, no. 6, pp. 1457–1465, 2006. View at: Publisher Site | Google Scholar - E. Ercenk, U. Sen, and S. Yilmaz, “Structural characterization of plasma sprayed basalt-SiC glass-ceramic coatings,”
*Ceramics International*, vol. 37, no. 3, pp. 883–889, 2011. View at: Publisher Site | Google Scholar - S. Das, A. K. Mukhopadhyay, S. Datta, G. C. Das, and D. Basu, “Hard glass-ceramic coating by microwave processing,”
*Journal of the European Ceramic Society*, vol. 28, no. 4, pp. 729–738, 2008. View at: Publisher Site | Google Scholar - L. Zhao, W. Q. Chen, and C. F. Lü, “Symplectic elasticity for bi-directional functionally graded materials,”
*Mechanics of Materials*, vol. 54, pp. 32–42, 2012. View at: Publisher Site | Google Scholar - A. J. Goupee and S. S. Vel, “Transient multiscale thermoelastic analysis of functionally graded materials,”
*Composite Structures*, vol. 92, no. 6, pp. 1372–1390, 2010. View at: Publisher Site | Google Scholar - G. J. Nie, Z. Zhong, and S. Chen, “Analytical solution for a functionally graded beam with arbitrary graded material properties,”
*Composites B*, vol. 44, no. 1, pp. 274–282, 2013. View at: Publisher Site | Google Scholar - R. F. Gibson, “Modal vibration response measurements for characterization of composite materials and structures,”
*Composites Science and Technology*, vol. 60, no. 15, pp. 2769–2780, 2000. View at: Publisher Site | Google Scholar - K. Y. Maalawi, “Use of material grading for enhanced buckling design of thin-walled composite rings/long cylinders under external pressure,”
*Composite Structures*, vol. 93, no. 2, pp. 351–359, 2011. View at: Publisher Site | Google Scholar - I. M. Daniel and O. Ishai,
*Engineering Mechanics of Composite Materials*, Oxford University Press, New York, NY, USA, 2nd edition, 2006. - X.-H. Zhao and W. F. Chen, “The effective elastic moduli of concrete and composite materials,”
*Composites B*, vol. 29, no. 1, pp. 31–40, 1998. View at: Publisher Site | Google Scholar - K. B. Shin and S. H. Hahn, “Evaluation of the structural integrity of hybrid railway carriage structures including the ageing effects of composite materials,”
*Composite Structures*, vol. 68, no. 2, pp. 129–137, 2005. View at: Publisher Site | Google Scholar - Y. Li, H. Zhang, N. Zhang, and Y. Dai, “Stress analysis of functionally gradient beam using effective principal axes,”
*International Journal of Mechanics and Materials in Design*, vol. 2, no. 3-4, pp. 157–164, 2005. View at: Publisher Site | Google Scholar - Z. Zhong and T. Yu, “Analytical solution of a cantilever functionally graded beam,”
*Composites Science and Technology*, vol. 67, no. 3-4, pp. 481–488, 2007. View at: Publisher Site | Google Scholar - X. L. Chen and K. M. Liew, “Buckling of rectangular functionally graded material plates subjected to nonlinearly distributed in-plane edge loads,”
*Smart Materials and Structures*, vol. 13, no. 6, pp. 1430–1437, 2004. View at: Publisher Site | Google Scholar - E. Feldman and J. Aboudi, “Buckling analysis of functionally graded plates subjected to uniaxial loading,”
*Composite Structures*, vol. 38, no. 1–4, pp. 29–36, 1997. View at: Publisher Site | Google Scholar - K. D. Kim, Gilson Rescober Lomboy, and S. C. Han, “Geometrically non-linear analysis of Functionally Graded Material (FGM) plates and shells using a four-node quasi-conforming shell element,”
*Journal of Composite Materials*, vol. 42, no. 5, pp. 485–511, 2008. View at: Publisher Site | Google Scholar - Y. Obata and N. Noda, “Steady thermal stresses in a hollow circular cylinder and a hollow sphere of a functionally gradient material,”
*Journal of Thermal Stresses*, vol. 17, no. 3, pp. 471–487, 1994. View at: Publisher Site | Google Scholar - K. M. Liew, S. Kitipornchai, X. Z. Zhang, and C. W. Lim, “Analysis of the thermal stress behaviour of functionally graded hollow circular cylinders,”
*International Journal of Solids and Structures*, vol. 40, no. 10, pp. 2355–2380, 2003. View at: Publisher Site | Google Scholar | Zentralblatt MATH - X. Wang, R. T. Wu, and A. Atkinson, “Characterisation of residual stress and interface degradation in TBCs by photo-luminescence piezo-spectroscopy,”
*Surface and Coatings Technology*, vol. 204, no. 15, pp. 2472–2482, 2010. View at: Publisher Site | Google Scholar - A. C. Karaoglanli, E. Altuncu, I. Ozdemir, A. Turk, and F. Ustel, “Structure and durability evaluation of YSZ+Al
_{2}O_{3}composite TBCs with APS and HVOF bond coats under thermal cycling conditions,”*Surface and Coatings Technology*, vol. 205, no. 2, pp. S369–S373, 2011. View at: Publisher Site | Google Scholar - Y. H. Sohn, E. Y. Lee, B. A. Nagaraj, R. R. Biederman, and R. D. Sisson Jr., “Microstructural characterization of thermal barrier coatings on high pressure turbine blades,”
*Surface and Coatings Technology*, vol. 146-147, pp. 132–139, 2001. View at: Publisher Site | Google Scholar - C. Batista, A. Portinha, R. M. Ribeiro, V. Teixeira, M. F. Costa, and C. R. Oliveira, “Morphological and microstructural characterization of laser-glazed plasma-sprayed thermal barrier coatings,”
*Surface and Coatings Technology*, vol. 200, no. 9, pp. 2929–2937, 2006. View at: Publisher Site | Google Scholar - F. Cernuschi, P. Bison, and A. Moscatelli, “Microstructural characterization of porous thermal barrier coatings by laser flash technique,”
*Acta Materialia*, vol. 57, no. 12, pp. 3460–3471, 2009. View at: Publisher Site | Google Scholar - S.-I. Jung, J. H. Kim, J. H. Lee, Y. G. Jung, U. Paik, and K. S. Lee, “Microstructure and mechanical properties of zirconia-based thermal barrier coatings with starting powder morphology,”
*Surface and Coatings Technology*, vol. 204, no. 6-7, pp. 802–806, 2009. View at: Publisher Site | Google Scholar - A. Kulkarni, A. Goland, H. Herman et al., “Advanced neutron and X-ray techniques for insights into the microstructure of EB-PVD thermal barrier coatings,”
*Materials Science and Engineering A*, vol. 426, no. 1-2, pp. 43–52, 2006. View at: Publisher Site | Google Scholar - U. Schulz and M. Schmücker, “Microstructure of ZrO
_{2}thermal barrier coatings applied by EB-PVD,”*Materials Science and Engineering A*, vol. 276, no. 1-2, pp. 1–8, 2000. View at: Publisher Site | Google Scholar - X. Chen, L. Gu, B. Zou, Y. Wang, and X. Cao, “New functioanlly graded thermal barrier coating system based on LaMgAl11O19/YSZ prepared by air plasma spraying,”
*Surface and Coatings Technology*, vol. 206, no. 8-9, pp. 2265–2274, 2012. View at: Publisher Site | Google Scholar - U. Schulz, M. Peters, F. W. Bach, and G. Tegeder, “Graded coatings for thermal, wear and corrosion barriers,”
*Materials Science and Engineering A*, vol. 362, no. 1-2, pp. 61–80, 2003. View at: Publisher Site | Google Scholar - L. L. Shaw, “Thermal residual stresses in plates and coatings composed of multi-layered and functionally graded materials,”
*Composites B*, vol. 29, no. 3, pp. 199–210, 1998. View at: Google Scholar

#### Copyright

Copyright © 2013 Seokchan Kim et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.