#### Abstract

The quasistatic bending response is presented for a simply supported functionally graded rectangular plate subjected to a through-the-thickness temperature field under the effect of various theories of generalized thermoelasticity, namely, classical dynamical coupled theory, Lord and Shulman's theory with one relaxation time, and Green and Lindsay's theory with two relaxation times. The generalized shear deformation theory obtained by the first author is used. Material properties of the plate are assumed to be graded in the thickness direction according to a simple exponential law distribution in terms of the volume fractions of the constituents. The numerical illustrations concern quasistatic bending response of functionally graded square plates with two constituent materials are studied using the different theories of generalized thermoelasticity

#### 1. Introduction

In the recent years, functionally graded materials (FGMs) have gained considerable attention in many engineering applications. FGMs are considered as a potential structural material for future high-speed spacecraft and power generation industries. FGMs are new materials, microscopically inhomogeneous, in which the mechanical properties vary smoothly and continuously from one surface to the other [1ā7].

The effect of thermal loading on the displacement and stress fields for FGM plates and shells has been studied by a number of authors. For example, Wetherhold et al. [5] have considered the use of functionally graded materials to eliminate or control thermal deformation in beams and plates. Suresh and Mortensen [6] have discussed the large deformation of graded multilayered composites under mechanical and thermal loads. Praveen and Reddy [7] have investigated the response of functionally graded ceramic-metal plates using a finite element that accounts for the transverse shear strains, rotary inertia, and moderately large rotations in the von Karman sense.

The theory of thermoelasticity that includes the effect of temperature change has been well established. According to this theory, the temperature field is coupled with the elastic strain field. This theory covers a wide range of extensions of classical dynamical coupled thermoelasticity. Lord and Shulman [8] and Green and Lindsay [9] have extended the coupled theory of thermoelasticity by introducing the thermal relaxation times in the constitutive equations. Additional thermoelasticity theories have been presented and investigated by other researchers [10ā22].

Lord and Shulman [8] have considered isotropic solids and introduced one relaxation time parameter into the Fourier heat conduction equation; that is, both the heat flux and its time derivative are considered in the heat conduction equation. The heat equation associated with this theory is thus hyperbolic. A direct consequence is that the paradox of infinite speed of propagation inherent in both the uncoupled and coupled theories of classical thermoelasticity is eliminated and the heat wave feature can be modeled by the generalized thermoelasticity. The Green and Lindsay's theory [9] does not violate the Fourier's law of heat conduction when the body under consideration has a center of symmetry. In this theory, both the equations of motion and heat conduction are hyperbolic but the equation of motion is modified and differs from that of classical coupled thermoelasticity theory.

In this paper, a generalized nonclassical dynamic coupled thermoelasticity analysis is carried out on a functionally graded material plate. Through-the-thickness temperature distribution varying according to exponential law is considered and then temperature and stress behavior are presented for the mentioned plate. The governing equations of FGM plate for two-dimensional generalized thermoelastic problems are derived within the framework of the classical coupled theory, Lord and Shulman's theory, and Green and Lindsay's theory. The material properties of the functionally graded plate are assumed to vary continuously through the thickness according to an exponential law distribution of the volume fraction of the constituents. A generalized shear deformation theory is presented to obtain the governing equations. An exact solution for the coupled governing equations under simply supported boundary conditions is obtained. Numerical results are provided to show the influence of the material properties, and a temperature field on the displacement and stresses.

#### 2. Formulation of the Problem

We consider a solid rectangular plate of length width and thickness made of functionally graded materials. The material properties of the FGM plate are assumed to be function of the volume fraction of the constituent materials. Using the rectangular Cartesian coordinates we take the functional graded between the physical properties and for ceramic and metal FGM plate where and are the corresponding properties of ceramic (top surface) and metal (bottom surface), respectively.

The displacements of a material point located at in the plate may be written as
where and represent the differentiation with respect to and are the displacements corresponding to the co-ordinate system and are functions of the spatial co-ordinates, are the displacements along the axes and respectively, and and are the rotations about the *y*- and *x-*axes.

The strain components will be

where

#### 3. Theories of Thermoelasticity

In addition, the stress-strain-temperature relations for the linear thermoelastic materials are given, according to the generalized theories of thermoelasticity, by where and are the absolute temperature, reference temperature, strain tensor components, components of tensor of stress-temperature moduli, components of tensor of elastic moduli, and the first relaxation time of Green and Lindsay's theory, respectively.

In details, we can rewrite the stress components in the deferent theories of generalized thermoelasticity as follows:
where
and the material properties and are functions of In the absence of body forces and internal heat generation, the heat conduction equation will be in the form
where and are additional relaxation times. A comma followed by index *j* denotes partial differentiation with respect to the position of a material particle. A superimposed dot indicates partial derivative with respect to time In addition to the elastic coefficients and the material properties and are also functions of

It is clear that, by setting in (3.2) and in (3.4), we get the field equations for the conventional coupled theory of thermoelasticity; whereas when and the equations reduce to the Lord and Shulman's theory and when and and are nonvanishing, the equations reduce to the Green and Lindsay's theory.

Note that the stress-temperature modulus is given in terms of Young's modulus Poisson's ratio , and the thermal expansion coefficient by the relation

Generally, this study assumes that Youngās modulus *E*, Poissonās ratio material density thermal expansion coefficient specific heat capacity and thermal conductivity coefficient of the FGM change continuously through the thickness direction of the plate and obey the gradation relation given in (2.1).

#### 4. Solution of the Problem

To solve the problem, we obtain the stress and moment resultants for the FGM plate by integrating the stress components given in (3.2) over the thickness and written as

Here, the coefficients and are defined by

By using Hamilton's principle, the governing equations can be obtained in the form The edges of the plate are assumed to be simply supported and maintained at the reference temperature. That is, For the present problem, the solution for the change in temperature is sought in the form where and This temperature identically satisfies the boundary conditions given in (4.4) at the edges of the plate. The function will be obtained from the solution of the heat equation (3.4). In addition, we assume the following solution form for that satisfies the boundary conditions:

where and are arbitrary parameters. Finally, we get the stress components
Substitution for *T* from (4.5) into (3.4), with aid of (4.6), gives the following partial differential equation with variable coefficients:
Inserting into the above ordinary differential equation the material properties and from (2.1), we obtain
where

The exact closed form solution of (4.9) is given by where and are the integration constants, and the functions and are the modified Bessel's functions of the first and second kinds, respectively, in which

Note that the term represents the characteristic time of heat conduction through length The function is given by where Note that the constants and are given from the temperature boundary conditions at the lower and upper surfaces of the plate

#### 5. Numerical Results and Discussion

We present exact results for a simply supported FG square plate subjected to a transient thermal load. Since it is common in high-temperature applications to employ a ceramic top layer as a thermal barrier to a metallic structure, we choose the constituent materials of the FG plate to be Aluminum (Al) and Silicon (SiC) having the following material properties:

The dimensions of the simply supported FG plate are ām. In addition, the values of different parameters are used as follows: Numerical results are presented in terms of the nondimensional variables defined as The results of the classical coupled theory (C-T), Lord and Shulman's theory (L-S), and Green and Lindsay's theory (G-L) for various values of time parameter are listed in Table 1. The relaxation times for these theories are chosen to be The dimensionless temperature , displacement longitudinal stress transverse shear stress and in-plane shear stress are given at and 1/2, respectively. It can be found from Table 1 that the results obtained by the L-S model agree well with those obtained by C-T whereas the G-L model gives an accurate prediction of the results that slightly differ from the above two models. The temperature given in the context of all theories may be unchanged.

It is well known that G-L theory is accurate to predict temperature, displacement, and stresses, so some results have been plotted in Figures 1ā8. The through-the-thickness variation of the temperature and stresses is plotted in Figures 1, 2, 3, and 4 for different values of the time parameter The through-the-thickness variation of the longitudinal stress transverse shear stress and in-plane shear stress changes significantly as a function of time. For example, at the magnitude of the temperature and longitudinal stress is maximum at a point on the top surface of the plate. However at the maximum values of the transverse shear stress occur at *z* = 0.298 for lower values of The magnitude of the in-plane stress is maximum at a point on the top surface of the plate and it increases as increases. Since material properties and the temperature change vary through the thickness, then the through-the-thickness variation of the longitudinal stress is nonlinear and the maximum values of the transverse shear stress dose not occur at the center of the plate.

Now, we discuss the effect of relaxation times in G-L model. Let the first relaxation time and the second relaxation time be double of it. Figures 5, 6, 7, and 8 present the displacement and stresses versus the time parameter using different values of the relaxation time It is clear that all results are very sensitive to the variation of the relaxation time. The results may be inaccurate for higher values of

#### 6. Conclusion

In this paper, the numerical illustrations concern quasistatic bending response of FG square plates are studied in the context of the generalized thermoelasticity theories. A refined shear deformation theory is used for this purpose. Material properties of the plate are assumed to be graded in the thickness direction according to a simple exponential law distribution in terms of the volume fractions of the constituents. An exact solution for the present problem is obtained. Numerical results are provided to show the influence of the material properties, and a temperature field on the displacement and stresses.

From these results, we can conclude that

the results of G-L model give an accurate prediction comparing with those obtained by the other two models; the results of L-S model agree well with those obtained by C-T model;the temperature may be independent of the parameters used in the different theories; for higher values of the relaxation times, L-S and G-L models may be failed to get accurate solution comparing with the C-T model.#### Acknowledgment

The investigators would like to express their appreciation to the Deanship of Scientific Research at King AbdulAziz University for their financial support of this study, Grant no. 180/428.