Research Article  Open Access
Bending Analysis of Functionally Graded Plates in the Context of Different Theories of Thermoelasticity
Abstract
The quasistatic bending response is presented for a simply supported functionally graded rectangular plate subjected to a throughthethickness 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 highspeed 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 ceramicmetal 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. Throughthethickness 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 twodimensional 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 coordinate system and are functions of the spatial coordinates, are the displacements along the axes and respectively, and and are the rotations about the y and xaxes.
The strain components will be
where
3. Theories of Thermoelasticity
In addition, the stressstraintemperature 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 stresstemperature 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 stresstemperature 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 hightemperature 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 (CT), Lord and Shulman's theory (LS), and Green and Lindsay's theory (GL) 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 inplane shear stress are given at and 1/2, respectively. It can be found from Table 1 that the results obtained by the LS model agree well with those obtained by CT whereas the GL 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 GL theory is accurate to predict temperature, displacement, and stresses, so some results have been plotted in Figures 1–8. The throughthethickness variation of the temperature and stresses is plotted in Figures 1, 2, 3, and 4 for different values of the time parameter The throughthethickness variation of the longitudinal stress transverse shear stress and inplane 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 inplane 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 throughthethickness 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 GL 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 GL model give an accurate prediction comparing with those obtained by the other two models; the results of LS model agree well with those obtained by CT model;the temperature may be independent of the parameters used in the different theories; for higher values of the relaxation times, LS and GL models may be failed to get accurate solution comparing with the CT 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.
References
 A. M. Zenkour, “A comprehensive analysis of functionally graded sandwich plates—part 1: deflection and stresses,” International Journal of Solids and Structures, vol. 42, no. 1819, pp. 5224–5242, 2005. View at: Publisher Site  Google Scholar
 A. M. Zenkour, “A comprehensive analysis of functionally graded sandwich plates—part 2: buckling and free vibration,” International Journal of Solids and Structures, vol. 42, no. 1819, pp. 5243–5258, 2005. View at: Publisher Site  Google Scholar
 A. M. Zenkour, “Generalized shear deformation theory for bending analysis of functionally graded plates,” Applied Mathematical Modelling, vol. 30, no. 1, pp. 67–84, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 A. M. Zenkour, “Benchmark trigonometric and 3D elasticity solutions for an exponentially graded thick rectangular plate,” Archive of Applied Mechanics, vol. 77, no. 4, pp. 197–214, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 R. C. Wetherhold, S. Seelman, and J. Wang, “Use of functionally graded materials to eliminate or control thermal deformation,” Composites Science and Technology, vol. 56, no. 9, pp. 1099–1104, 1996. View at: Publisher Site  Google Scholar
 S. Suresh and A. Mortensen, “Functionally graded metals and metalceramic composites—part 2: thermal mechanical behavior,” International Materials Reviews, vol. 42, no. 3, pp. 85–116, 1997. View at: Google Scholar
 G. V. Praveen and J. N. Reddy, “Nonlinear transient thermoelastic analysis of functionally graded ceramicmetal plates,” International Journal of Solids and Structures, vol. 35, no. 33, pp. 4457–4476, 1998. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 H. W. Lord and Y. Shulman, “A generalized dynamical theory of thermoelasticity,” Journal of the Mechanics and Physics of Solids, vol. 15, no. 5, pp. 299–309, 1967. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 A. E. Green and K. A. Lindsay, “Thermolasticity,” Journal of Elasticity 2, vol. 2, pp. 1–7, 1972. View at: Publisher Site  Google Scholar
 J. Ignaczak, “Uniqueness in generalized thermoelasticity,” Journal of Thermal Stresses, vol. 2, no. 2, pp. 171–175, 1979. View at: Publisher Site  Google Scholar
 R. S. Dhaliwal and A. Singh, Dynamic Coupled Thermoelasticity, Hindustan Publishing, New Delhi, India, 1980.
 R. Dhaliwal and H. Sherief, “Generalized thermoelasticity for anisotropic media,” Quarterly of Applied Mathematics, vol. 33, no. 1, pp. 1–8, 1980. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 D. S. Chandrasekharaiah, “Thermoelasticity with second sound,” Applied Mechanics Reviews, vol. 39, pp. 355–376, 1986. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 R. S. Dhaliwal and H. S. Saxena, “Generalized magnetothermoelastic waves in an infinite elastic solid with a cylindrical cavity,” Journal of Thermal Stresses, vol. 14, no. 4, pp. 353–369, 1991. View at: Publisher Site  Google Scholar  MathSciNet
 A. E. Green and P. M. Naghdi, “Thermoelasticity without energy dissipation,” Journal of Elasticity, vol. 31, no. 3, pp. 189–208, 1993. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 H. Sherief and R. Dhaliwal, “A uniqueness theorem and a variational principle for generalized thermoelasticity,” Journal of Thermal Stresses, vol. 3, no. 2, pp. 223–231, 1980. View at: Publisher Site  Google Scholar
 H. Sherief and R. Dhaliwal, “Generalized thermoelasticity for anisotropic media,” Quarterly of Applied Mathematics, vol. 33, no. 1, pp. 1–8, 1980/81. View at: Google Scholar  MathSciNet
 K. Daneshjoo and M. Ramezani, “Coupled thermoelasticity in laminated composite plates based on GreenLindsay model,” Composite Structures, vol. 55, no. 4, pp. 387–392, 2002. View at: Publisher Site  Google Scholar
 A. E. Green and N. Laws, “On the entropy production inequality,” Archive for Rational Mechanics and Analysis, vol. 45, no. 1, pp. 47–53, 1972. View at: Google Scholar  MathSciNet
 R. B. Hetnarski and J. Ignaczak, “Nonclassical dynamical thermoelasticity,” International Journal of Solids and Structures, vol. 37, no. 12, pp. 215–224, 2000. View at: Publisher Site  Google Scholar  MathSciNet
 A. Bagri and M. R. Eslami, “Generalized coupled thermoelasticity of functionally graded annular disk considering the LordShulman theory,” Composite Structures, vol. 83, no. 2, pp. 168–179, 2008. View at: Publisher Site  Google Scholar
 A. S. ElKaramany and M. A. Ezzat, “Thermal shock problem in generalized thermoviscoelasticity under four theories,” International Journal of Engineering Science, vol. 42, no. 7, pp. 649–671, 2004. View at: Publisher Site  Google Scholar  MathSciNet
Copyright
Copyright © 2009 A. M. Zenkour 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.