- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Article Processing Charges ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Recently Accepted Articles ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents

Modelling and Simulation in Engineering

Volume 2012 (2012), Article ID 159806, 10 pages

http://dx.doi.org/10.1155/2012/159806

## A New Hyperbolic Shear Deformation Theory for Bending Analysis of Functionally Graded Plates

^{1}Departement of Civil Engineering, Ibn Khaldoun University of Tiaret, BP 78 Zaaroura, 14000 Tiaret, Algeria^{2}Laboratoire des Matériaux et Hydrologie, Université de Sidi Bel Abbes, BP 89 Cité Ben M’hidi, 22000 Sidi Bel Abbes, Algeria

Received 16 May 2012; Accepted 3 August 2012

Academic Editor: Guowei Wei

Copyright © 2012 Tahar Hassaine Daouadji 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.

#### Abstract

Theoretical formulation, Navier’s solutions of rectangular plates based on a new higher order shear deformation model are presented for the static response of functionally graded plates. This theory enforces traction-free boundary conditions at plate surfaces. Shear correction factors are not required because a correct representation of transverse shearing strain is given. Unlike any other theory, the number of unknown functions involved is only four, as against five in case of other shear deformation theories. The mechanical properties of the plate are assumed to vary continuously in the thickness direction by a simple power-law distribution in terms of the volume fractions of the constituents. Numerical illustrations concern flexural behavior of FG plates with metal-ceramic composition. Parametric studies are performed for varying ceramic volume fraction, volume fractions profiles, aspect ratios, and length to thickness ratios. Results are verified with available results in the literature. It can be concluded that the proposed theory is accurate and simple in solving the static bending behavior of functionally graded plates.

#### 1. Introduction

The concept of functionally graded materials (FGMs) was first introduced in 1984 by a group of material scientists in Japan, as ultrahigh temperature-resistant materials for aircraft, space vehicles, and other engineering applications. Functionally graded materials (FGMs) are new composite materials in which the microstructural details are spatially varied through nonuniform distribution of the reinforcement phase. This is achieved by using reinforcement with different properties, sizes, and shapes, as well as by interchanging the role of reinforcement and matrix phase in a continuous manner. The result is a microstructure that produces continuous or smooth change on thermal and mechanical properties at the macroscopic or continuum level (Koizumi, 1993 [1]; Hirai and Chen, 1999 [2]). Now, FGMs are developed for general use as structural components in extremely high-temperature environments. Therefore, it is important to study the wave propagation of functionally graded materials structures in terms of nondestructive evaluation and material characterization.

Several studies have been performed to analyze the mechanical or the thermal or the thermomechanical responses of FG plates and shells. A comprehensive review is done by Tanigawa (1995) [3]. Reddy (2000) [4] has analyzed the static behavior of functionally graded rectangular plates based on his third-order shear deformation plate theory. Cheng and Batra (2000) [5] have related the deflections of a simply supported FG polygonal plate given by the first-order shear deformation theory and third-order shear deformation theory to that of an equivalent homogeneous Kirchhoff plate. The static response of FG plate has been investigated by Zenkour (2006) [6] using a generalized shear deformation theory. In a recent study, Şimşek (2010) [7] has studied the dynamic deflections and the stresses of an FG simply supported beam subjected to a moving mass by using Euler Bernoulli, Timoshenko, and the parabolic shear deformation beam theory. Şimşek (2010) [8], Benchour et al. [9], and Abdelaziz et al. 2010 [10] studied the free vibration of FG beams having different boundary conditions using the classical, the first-order, and different higher-order shear deformation beam and plate theories. The nonlinear dynamic analysis of an FG beam with pinned-pinned supports due to a moving harmonic load has been examined by Şimşek (2010) [11] using Timoshenko beam theory.

The primary objective of this paper is to present a general formulation for functionally graded plates (FGPs) using a new higher-order shear deformation plate theory with only four unknown functions. The present theory satisfies equilibrium conditions at the top and bottom faces of the plate without using shear correction factors. The hyperbolic function in terms of thickness coordinate is used in the displacement field to account for shear deformation. Governing equations are derived from the principle of minimum total potential energy. Navier solution is used to obtain the closed-form solutions for simply supported FG plates. To illustrate the accuracy of the present theory, the obtained results are compared with three-dimensional elasticity solutions and results of the first-order, and the other higher-order theories (Table 1).

In this study, a new displacement models for an analysis of simply supported FGM plates are proposed. The plates are made of an isotropic material with material properties varying in the thickness direction only. Analytical solutions for bending deflections of FGM plates are obtained. The governing equations are derived from the principle of minimum total potential energy. Numerical examples are presented to illustrate the accuracy and efficiency of the present theory by comparing the obtained results with those computed using various other theories.

#### 2. Problem Formulation

Consider a plate of total thickness and composed of functionally graded material through the thickness. It is assumed that the material is isotropic, and grading is assumed to be only through the thickness. The plane is taken to be the undeformed mid plane of the plate with the -axis positive upward from the mid plane (Figure 1).

##### 2.1. Displacement Fields and Strains

The assumed displacement field is as follows: where and are the mid-plane displacements of the plate in the and directions, respectively; and are the bending and shear components of transverse displacement, respectively, while represents shape functions determining the distribution of the transverse shear strains and stresses along the thickness and is given as It should be noted that unlike the first-order shear deformation theory, this theory does not require shear correction factors. The kinematic relations can be obtained as follows: where

##### 2.2. Constitutive Relations

In FGM, material property gradation is considered through the thickness, and the expression given below represents the profile for the volume fraction where denotes a generic material property like modulus, and denote the property of the top and bottom faces of the plate, respectively, and is a parameter that dictates material variation profile through the thickness. Here, it is assumed that modules and vary according to the equation (5), and is assumed to be a constant. The linear constitutive relations are where

##### 2.3. Governing Equations

The governing equations of equilibrium can be derived by using the principle of virtual displacements. The principle of virtual work in the present case yields where is the top surface and is the applied transverse load.

Substituting and (6) into (8) and integrating through the thickness of the plate, (8) can be rewritten as whereThe governing equations of equilibrium can be derived from (9) by integrating the displacement gradients by parts and setting the coefficients , , , and zero separately. Thus one can obtain the equilibrium equations associated with the present shear deformation theory as follows: Using (6) in (10a) and (10b) the stress resultants of a sandwich plate made up of three layers can be related to the total strains by where where , , and so forth are the plate stiffness, defined bySubstituting from (12) into , we obtain the following equation:where , , and are the following differential operators:

##### 2.4. Exact Solution for a Simply Supported FGM Plate

Rectangular plates are generally classified in accordance with the type of support used. We are here concerned with the exact solution of (15a)–(15d) for a simply supported FG plate. The following boundary conditions are imposed at the side edges:To solve this problem, Navier assumed the transverse mechanical and temperature loads, , in the form of a double trigonometric series as where , , and represents the intensity of the load at the plate center.

Following the Navier solution procedure, we assume the following solution form for , , and that satisfies the following boundary conditions:
where *, **, *, and are arbitrary parameters to be determined subjected to the condition that the solution in (19) satisfies governing equations (15a)–(15d). One obtains the following operator equation:
where and is the symmetric matrix given by
in which

#### 3. Numerical Results and Discussions

The study has been focused on the static behavior of functionally graded plate based on the present new higher-order shear deformation model. Here, some representative results of the Navier solution obtained for a simply supported rectangular plate are presented.

A functionally graded material consisting of aluminum-alumina is considered. The following material properties are used in computing the numerical values (Bouazza et al. [14]).(i)Metal (aluminum, Al): GPa; .(ii)Ceramic (alumina, Al_{2}O_{3}): GPa; .Now, a functionally graded material consisting of aluminum and alumina is considered. Young's modulus for aluminum is 70 GPa while for alumina is 380 GPa. Note that, Poisson's ratio is selected constant for both and equal to 0.3. The various nondimensional parameters used are

It is clear that the deflection increases as the side-to-thickness ratio decreases. The same results were obtained in most literatures. In addition, the correlation between the present new higher-order shear deformation theory and different higher-order and first-order shear deformation theories is established by the author in his recent papers. It is found that this theory predicts the deflections and stresses more accurately when compared to the first- and third-order theories.

For the sake of completeness, results of the present theory are compared with those obtained using a new Navier-type three-dimensionally exact solution for small deflections in bending of linear elastic isotropic homogeneous rectangular plates. The center deflection and the distribution across the plate thickness of in-plane longitudinal stress and longitudinal tangential stress are compared with the results of the 3D solution and are shown in Tables 2 and 3. The present solution is realized for a quadratic plate, with the following fixed data: and three values for the plate thickness: , and . It is to be noted that the present results compare very well with the 3D solution. All deflections again compare well with the 3D solution and show good convergence with the average 3D solution.

In Table 4, the effect of volume fraction exponent on the dimensionless stresses and displacements of an FGM square plate () is given. This table shows comparison between results for plates subjected to uniform or sinusoidal distributed loads, respectively. As it is well known, the uniform load distribution always overpredicts the displacements and stresses magnitude. As the plate becomes more and more metallic, the difference increases for deflection and in-plane longitudinal stress while it decreases for in-plane normal stress . It is important to observe that the stresses for a fully ceramic plate are the same as that for a fully metal plate. This is because the plate for these two cases is fully homogeneous, and the stresses do not depend on the modulus of elasticity. Results in Table 4 should serve as benchmark results for future comparisons.

Tables 5 and 6 compare the deflections and stresses of different types of the FGM square plate () and FGM rectangular plate .The deflections decrease as the aspect ratio increases and this irrespective of the type of the FGM plate. All theories (SSDPT, PSDPT, and NHPSDT) give the same axial stress and for a fully ceramic plate . In general, the axial stress increases with the volume fraction exponent . The transverse shear stress for a FGM plate subjected to a distributed load. The results show that the transverse shear stresses may be indistinguishable. As the volume fraction exponent increases for FGM plates, the shear stress will increase, and the fully ceramic plates give the smallest shear stresses.

Figures 2 and 3 show the variation of the center deflection with the aspect and side-to-thickness ratios, respectively. The deflection is maximum for the metallic plate and minimum for the ceramic plate. The difference increases as the aspect ratio increases while it may be unchanged with the increase of side-to-thickness ratio. One of the main inferences from the analysis is that the response of FGM plates is intermediate to that of the ceramic and metal homogeneous plates (see also Table 4). It is to be noted that, in the case of thermal or combined loads and under certain conditions, the above response is not intermediate.

Figures 7 and 8 depict the through-the-thickness distributions of the shear stresses and , the in-plane longitudinal and normal stresses and , and the longitudinal tangential stress in the FGM plate under the uniform load. The volume fraction exponent of the FGM plate is taken as in these figures. Distinction between the curves in Figures 8 and 9 is obvious. As strain gradients increase, the inhomogeneities play a greater role in stress distribution calculations. The through-the-thickness distributions of the shear stresses and are not parabolic, and the stresses increase as the aspect ratio decreases. It is to be noted that the maximum value occurs at , not at the plate center as in the homogeneous case.

As exhibited in Figures 5 and 6, the in-plane longitudinal and normal stresses, and , are compressive throughout the plate up to and then they become tensile. The maximum compressive stresses occur at a point on the bottom surface and the maximum tensile stresses occur, of course, at a point on the top surface of the FGM plate. However, the tensile and compressive values of the longitudinal tangential stress, (cf. Figure 7), are maximum at a point on the bottom and top surfaces of the FGM plate, respectively. It is clear that the minimum value of zero for all in-plane stresses , and occurs at and this is irrespective of the aspect and side-to-thickness ratios.

Finally, the exact maximum deflections of simply supported FGM square plate are compared in Figure 4 for various ratios of module, (for a given thickness, ). This means that the deflections are computed for plates with different ceramic-metal mixtures. It is clear that the deflections decrease smoothly as the volume fraction exponent decreases and as the ratio of metal-to-ceramic modules increases.

#### 4. Conclusion

In this study, a new higher-order shear deformation model is proposed to analyze the static behavior of functionally graded plates. Unlike any other theory, the theory presented gives rise to only four governing equations resulting in considerably lower computational effort when compared with the other higher-order theories reported in the literature having more number of governing equations. Bending and stress analysis under transverse load were analyzed, and results were compared with previous other shear deformation theories. The developed theories give parabolic distribution of the transverse shear strains and satisfy the zero traction boundary conditions on the surfaces of the plate without using shear correction factors. The accuracy and efficiency of the present theories have been demonstrated for static behavior of functionally graded plates. All comparison studies demonstrated that the deflections and stresses obtained using the present new higher-order shear deformation theories (with four unknowns) and other higher shear deformation theories such as PSDPT and SSDPT (with five unknowns) are almost identical. The extension of the present theory is also envisaged for general boundary conditions and plates of a more general shape. In conclusion, it can be said that the proposed theory NHPSDT is accurate and simple in solving the static behaviors of FGM plates.

#### References

- M. Koizumi, “The concept of FGM,”
*Ceramic Transactions, Functionally Gradient Materials*, vol. 34, pp. 3–10, 1993. - T. Hirai and L. Chen, “Recent and prospective development of functionally graded materials in Japan,”
*Materials Science Forum*, vol. 308–311, pp. 509–514, 1999. - Y. Tanigawa, “Some basic thermoelastic problems for nonhomogeneous structural materials,”
*Applied Mechanics Reviews*, vol. 48, no. 6, pp. 287–300, 1995. View at Scopus - J. N. Reddy, “Analysis of functionally graded plates,”
*International Journal for Numerical Methods in Engineering*, vol. 47, no. 1–3, pp. 663–684, 2000. View at Scopus - Z. Q. Cheng and R. C. Batra, “Deflection relationships between the homogeneous Kirchhoff plate theory and different functionally graded plate theories,”
*Archives of Mechanics*, vol. 52, no. 1, pp. 143–158, 2000. View at Scopus - 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 · View at Google Scholar · View at Scopus - M. Şimşek, “Vibration analysis of a functionally graded beam under a moving mass by using different beam theories,”
*Composite Structures*, vol. 92, no. 4, pp. 904–917, 2010. View at Publisher · View at Google Scholar · View at Scopus - M. Şimşek, “Fundamental frequency analysis of functionally graded beams by using different higher-order beam theories,”
*Nuclear Engineering and Design*, vol. 240, no. 4, pp. 697–705, 2010. View at Publisher · View at Google Scholar · View at Scopus - A. Benachour, H. D. Tahar, H. A. Atmane, A. Tounsi, and M. S. Ahmed, “A four variable refined plate theory for free vibrations of functionally graded plates with arbitrary gradient,”
*Composites B*, vol. 42, no. 6, pp. 1386–1394, 2011. View at Publisher · View at Google Scholar · View at Scopus - H. H. Abdelaziz, H. A. Atmane, I. Mechab, L. Boumia, A. Tounsi, and A. B. E. Abbas, “Static analysis of functionally graded sandwich plates using an efficient and simple refined theory,”
*Chinese Journal of Aeronautics*, vol. 24, no. 4, pp. 434–448, 2011. View at Publisher · View at Google Scholar · View at Scopus - M. Şimşek, “Non-linear vibration analysis of a functionally graded Timoshenko beam under action of a moving harmonic load,”
*Composite Structures*, vol. 92, no. 10, pp. 2532–2546, 2010. View at Publisher · View at Google Scholar · View at Scopus - S. P. Timoshenko and S. Woinowsky-Krieger,
*Theory of Plates and Shells*, McGraw-Hill, New York, NY, USA, 1959. - H. Werner, “A three-dimensional solution for rectangular plate bending free of transversal normal stresses,”
*Communications in Numerical Methods in Engineering*, vol. 15, no. 4, pp. 295–302, 1999. View at Scopus - M. Bouazza, A. Tounsi, E. A. Adda Bedia, and M. Meguenni, “Stability analysis of fonctionnally graded plates subject to thermal load,”
*Advanced Structures Materials B*, vol. 15, pp. 669–680, 2011.