BENDING AND FREE VIBRATION ANALYSIS OF FUNCTIONALLY GRADED PLATES USING NEW EIGHT- UNKNOWN SHEAR DEFORMATION THEORY BY FINITE ELEMENT METHOD

In this paper, a new eight-unknown shear deformatio n theory is developed for bending and free vibration analysis of functionally graded plat es by finite element method. The theory based on full twelve-unknown higher order shear deformati on theory, simultaneously satisfy zeros transverse stresses at top and bottom surface of FG plates. A four-node rectangular element with sixteen degrees of freedom per node is used. Poisso n’s ratios, Young’s moduli and material densities vary continuously in thickness direction according to the volume fraction of constituents which is modeled as power law function s. Results are verified with available results in the literature. Parametric studies are performed for different power law index, side-tothickness ratios.


INTRODUCTION
Since it was invented by Japanese scientists in 1984 [1], functionally graded materials (FGMs) are increasingly and widely used in many fields, such as aerospace, marine, mechanical, and structural engineering due to its advantages compared to classical fiberreinforced laminated composites. The typical FGMs composed of ceramic and metal materials. The ceramic composition offers thermal barrier effects and protects the metal from corrosion and oxidation, and the metallic composition provides FGM toughness and strength.
For dynamic and static analysis of functionally graded plates and shells, many plate theories are developed. A review of shear deformation theories for isotropic and laminated plates was carried out by Ghugal and Shimpi [2] and Khandan et al. [3]. Focus on modeling of functionally graded plates and shells, Thai Huu-Tai and Kim Seung-Eock [4] reviewed various theoretical models to investigate their mechanical behavior. The classical plate theory (CPT) based on Kirchhoff assumptions and ignores the transverse shear deformation effect gives appropriate results for thin plates. First-order shear deformation theory (FSDT) takes into account the transverse shear deformation effect and needs a shear correction factor which is difficult to determine due to its dependence on many parameters. To overcome the weaknesses of FSDT, the higher-order shear deformation theories are proposed.
A comprehensive review of the various methods employed to study the static, dynamic and stability behavior of functionally graded plates can be found in work of Swaminathan et al. [5]. The review focuses on comparing the stress, vibration and buckling characteristics of FGM plates using different theories. Based on third order shear deformation theory with five displacement unknowns, Reddy [6] developed analytical and finite element solutions for static and dynamic analysis of functionally graded rectangular plates. El-Abbasi and Meguidin [7] used a new thick shell element to study the thermoelastic behavior of functionally graded plates and shells. They extended the four-nodded seven-parameter shell element to account for the varying elastic and thermal properties, as well as the temperature boundary conditions on both faces of FG plates and shells Oyekoya et al. [8] developed Mindlin-type element and Reissner-type element for the modelling of functionally graded plate subjected to buckling and free vibration. The Mindlintype element formulation is based on averaging of transverse shear distribution over plate thickness using Lagrangian interpolation. The Reissner-type element formulation is based on parabolic transverse shear distribution over plate thickness using Lagrangian and Hermitian interpolation. Talha and Singh [9] studied free vibration and static behavior of functionally graded plates using higher order shear deformation theory. A continuous isoparametric Lagrangian finite element with 13 degrees of freedom per node is employed for the modeling of functionally graded plates. Thai Huu-Tai and Choi Dong-Ho [10] presented finite element formulation of various four-unknown shear deformation theories for the bending and vibration analyses of functionally graded plates. To describe the primary variables, a four-node quadrilateral finite element is developed using Lagrangian and Hermitian interpolation functions. Three-dimensional graded finite element method based on Rayleigh-Ritz energy formulation has been applied to study the static response of the thick functionally graded plates [11].
In this paper, a new higher order displacement field based on twelve-unknown higher order shear deformation theory is developed to analyze the free vibration and buckling of functionally graded plates. The new eight-unknown higher order shear deformation theory is derived from the satisfaction of vanishing transverse shear stress at the top and bottom surfaces of the plate. The finite element model is developed for bending and free vibration analysis of power-law functionally graded plates. A C 1 continuous four-node quadrilateral plate element with sixteen degrees of freedom per node is employed. Lagrangian linear interpolation functions are used to describe the in-plane displacements and the rotation of normals about x, y axes; Hermitian cubic interpolation functions are given for the transverse displacement, rotation about z-axis, higherorder term of displacements and their first derivation.

KINEMATICS
The twelve-unknown higher order displacement field is given as follow [12]: where , , u v w denote the displacements of a point along the (x, y, z) coordinates. 0 0 0 , , u v w are corresponding displacements of a point on the midplane. x θ , y θ and z θ are the rotations of the line segment normal to the midplane about the y-axis, x-axis and z-axis , respectively. The functions 0 u * , 0 v * , 0 w * , x θ * , y θ * and z θ * are the higher order terms in the Taylor series expansion defined in the mid-plane.
For bending plates, the transverse shear stresses xz σ , yz σ must be vanished at the top and bottom surfaces. These conditions lead to the requirement that the corresponding transverse strains on these surfaces be zero. From , , , .
Thus, the displacement field (1) becomes: with: . 4 h c c h = = or in matrix notation as: where: displacement vector of any generic point within the plate; Following strain -displacement relation, the non-zero strains are given as: or: where: , , , , , , ; ; xz

CONSTITUTIVE EQUATION
Consider a rectangular FGM plate with the length a, width b, and thickness h. The x-, y-, and z-coordinates are taken along the length, width, and height of the plate, respectively, as shown in Fig. 1. The material properties of FGM plates are assumed to vary continuously through the thickness of the plate by a power law distribution as [6]: where V(z) represents the effective material property such as Young's modulus E, mass density ρ, and Poisson's ratio ν; subscripts m and c represent the metallic and ceramic constituents, respectively; and p is the volume fraction exponent. The stress-strain relationship for the FGM plate can be written as:  (12) in which:

FINITE ELEMENT FORMULATION
A C 1 continuous four-node quadrilateral plate bending element with sixteen degrees of freedom per node is used (Fig. 2). The Lagrangian linear interpolation functions ( ) ; ; ; , , ; T T  The displacement vector at any generic point can be written as: where: The strain vector is expressed by: The Hamilton's principle can be expressed as: and applying for each element: The strain energy of the FGM plate element is given by:

{ } { } { } [ ] [ ][ ]{ } { } [ ]{ }
The external work done on the plate element by distributed applied load may be written as: and {f} is mechanical load vector. The kinetic energy of the FGM plate can be expressed as:

NUMERICAL RESULTS
Matlab codes for finite element model have been built for numerical investigation. After checking convergence, a 10×10 mesh of four-node element has been used in the computation. The selective integration scheme based on Gauss-quadrature rules, with 3×3 for membrane, coupling, flexure and inertia terms and 2×2 for shear term. A rectangular FG plates with different boundary conditions as shown in Fig. 3 are considered (F-free, S-simply supported, and Cclamped). Material properties of the P-FG plate are given in Table 1. For convenience, the following dimensionless forms are used [13]: Table 1. Material properties used in the P-FG plate [13].

Example 1. Validation study
Dimensionless central deflections w of isotropic square plates (p = 0) with various values of thickness ratios a/h are presented in Table 2. The present results are compared with the solutions given by Thai, H.T., & Choi, D.H. [10] based on four-unknown shear deformation theories (zeros shape function -FSDT) and the analytical solutions reported by Zenkour [14] based on a mixed first-order shear deformation theory (MPT). It can be seen that the present solution is in close agreement with those solutions (errors <0.2 %).
Dimensionless fundamental frequencies ω of simply supported (SSSS) square FG plates (p = 0) with various values of thickness ratios a/h and power law index p are presented in Table 3. The comparison of the dimensionless fundamental frequencies of present results shows good agreement with analytical solutions of Thai H. T., & Kim S. E. [12] based on simple higher-order theory, and finite element results of Thai H. T., & Choi D.H. [9] based on four unknowns shear deformation theories.       Example 2. Effect of power law index p and side-to-thickness ratio a/h on the dimensionless central deflection w . In this example, the square FG plate with different boundary conditions under uniformly distributed load is considered. The calculated dimensionless central deflection with various power law index p = 0; 0.5; 1.0; 2; 5; 10 and a/h = 5; 10; 20; 50 are given in Table 4. Figures 4 and 5 show the variation of power law index p and side-to-thickness ratio a/h versus dimensionless central deflection. It is found that the dimensionless central deflection increases as power law index p increases, while dimensionless central deflection decreases as side-to-thickness ratio increase with all types of boundary conditions. Example 3. Effect of power law index p and side-to-thickness ratio a/h on the fundamental frequency ω Table 5 presents the dimensionless fundamental frequency for various power law index p = 0; 0.5; 1.0; 2; 5; 10 and a/h = 5; 10; 20; 50. Different boundary condition for each case is considered. The variation of dimensionless fundamental frequency versus power law index p and side-to-thickness ratio a/h is illustrated in Figures 6 and 7. It is observed that, for all types of boundary condition, dimensionless frequencies decreases as power law index and side-to-thickness ration increases. Effect of boundary conditions is clearly too, the dimensionless frequency of FG plate with boundary conditions CCCC is highest, and the lowest with SSSS boundary conditions.

CONCLUSIONS
In this study, the new eight-unknown shear deformation theory is used to analyze the bending and free vibration of rectangular fuctionally graded plates by finite element approach. The governing equations and boundary conditions are derived by employing the Hamilton's principle. Validation studies have been carried out to confirm the accuracy of the present formulation. The obtained result shows a good agreement with those available in the literature. Influence of power law index, side-to-thickness ratio on bending and vibration responses of FG plates have been investigated and discussed. The new eight unknowns shear deformation theory is accurate in predicting static and free vibration responses of FG plates.