Nonlinear Electromechanical Stability of a Functionally Graded Circular Plate Integrated With Functionally Graded Piezoelectric Layers 1

This research develops nonlinear electromechanical stability of a circular functionally graded plate integrated with functionally graded piezoelectric layers under compressive radial force. Geometric nonlinearity is considered in the strain-displacement relation using Von-Karman relation. The structure is loaded under mechanical and electrical loads. Distribution of electric potential is considered along the radial and thickness direction. The top and bottom of both piezoelectric layers is short-circuited. The effect of various values of non homogenous index for both functionally graded (FG) and functionally graded piezoelectric (FGP) layers can be considered on the responses of the system. Furthermore, a comprehensive study for evaluation of geometric parameters can be performed on the critical loads of the structure.

For example a combination of ceramic and metal can be used as functionally graded materials.The material properties are changed gradually from metal to ceramic.This change may be described by a function along the thickness or other dimension of the structure.Composition of functionally graded materials with piezoelectric elements proposes new intelligent materials that can be studied in this paper.The piezoelectric effect has been presented scientifically by Pierre and Jacques Curie in 1880.Piezoelectric structures are very applicable in the industrial systems as sensor or actuator in various geometries such as plates, cylinders and shells.In order to control the distribution of the displacement or electric potential in a piezoelectric structure, functionally graded piezoelectric material (FGPM) can be used.An investigation on the literature can justifies necessity of this research.Wu et al. (2002) have used GDQR (generalized differential quadrature rule) for free vibration analysis of solid circular plates.O¨ zakc et al. (2003) have focused on the buckling load optimization of variable thickness circular and annular plates using finite element approach.Ma and Wang (2003) have presented deflection bending of a functionally graded circular plate using the classical nonlinear von Karman plate theory.The plate has been subjected to various types of loading such as mechanical and thermal loadings.As the results of that study, nonlinear bending and critical buckling temperature and thermal post-buckling behavior of the FGM plates were discussed.The stability of parametric vibrations of circular plate subjected to in-plane forces was analyzed using the Liapunov method by Tylikowski and Frischmuth (2003).Zhou et al. (2003) employed Chebyshev-Ritz method in three-dimensional free vibration analysis of circular and annular plates.They used a linear analysis with small strain assumption.Based on the geometric properties of circular and annular plates, the vibration was divided into three distinct categories: axisymmetric vibration, torsional vibration and circumferential vibration.Li et al. (2004) used a large deflection bending analysis of an axisymmetric simply supported circular plate.The incremental load technique was developed for solving the bending problem of a thin circular plate with large deflection.They have found that the employed technique has capability in solution of engineering problems.A circular plate containing piezoelectric layers as actuator under static and dynamic mechanical and electrical loads and using Kirchhoff plate model have been sudied by Sekouri et al. (2004).Experiments using a thin circular aluminum plate structure with distributed piezoelectric actuators were also conducted to verify the analysis and the computer simulations.Kang et al. (2005) presented a closed form solution for finding the natural frequencies and mode shapes of a circular orthotropic plate using Rayleigh-Ritz method.They have found that the obtained results have capability in designing of circular plates such as wood disk as an orthotropic composite material.Nosier and Fallah (2009) studied axisymmetric and asymmetric behavior of functionally graded circular plates under transverse mechanical loading using the first-order shear deformation plate theory with von Karman non-linearity.By introducing a stress function and a potential function, the problem were uncoupled to form equations describing the interior and edge-zone problems of FG plates.A perturbation technique, in conjunction with Fourier series method to model the problem asymmetries, is used to obtain the solution for various clamped and simply supported boundary conditions.Camier et al. (2009) studied Large-amplitude and geometrically nonlinear vibrations of free-edge circular plates with geometric imperfections.Free vibration analysis of circular thin plate with three types of boundary conditions have been studied by Yalcin et al. (2009).The solution procedure has been performed using differential transform method.The obtained results using this semi-numerical-analytical solution technique have been compared with results of Bessel function solution.Axisymmetric bending and buckling of perfect functionally graded solid circular plates have been studied by Saidi et al. (2009) based on the unconstrained third-order shear deformation plate theory (UTST).The obtained results have been compared with those results extracted using lower order shear deformation theories.Sahraee and Saidi (2009) investigated axisymmetric bending and stretching of functionally graded (FG) circular plates subjected to uniform transverse loading based on fourth-order shear deformation plate theory (FOST).Gradation of used material has been considered along the thickness direction based on a power law function.Vivio and Vullo (2010) introduced a new analytical method for evaluation of elastic stresses and deformations in the solid and annular circular plates with variable thickness subjected to transverse loading.Arefi (2013) and Arefi and Nahas (2014) presented nonlinear analysis of the functionally graded piezoelectric cylinder and sphere, respectively.Three dimensional analysis of a functionally graded piezoelectric shell under multi fields has been studied by Arefi (2014).
An investigation on the literature indicates that there is no published work to study the nonlinear electromechanical stability of a functionally graded circular plate integrated with functionally graded piezoelectric materials under mechanical and electrical loads.Some useful information about linear and nonlinear analysis of functionally graded piezoelectric materials can be considered in the literature [Asemi et al. 2014, Arefi and Rahimi 2011, Arefi and Rahimi 2012 (a, b, c, d, e, f), Rahimi et al, 2011, Khoshgoftar et al, 2009].

FORMULATION
Fundamental equations for mechanical stability of the functionally graded plates integrated with two functionally graded piezoelectric layers at top and bottom of plate is developed in the present section.Classical plate theory (CPT) is used for simulation of deformation components of the plate. ,,rz  is used for components of coordinate system and , uw is used for symmetric components of plate deformation.Based on the above assumptions, the deformation components can be given as follows: [Arefi and Rahimi, 2012] The strain components can be derived from nonlinear strain-displacement relation (Von-Karman) [Arefi and Rahimi, 2012 Using above relation, three in-plane components of strain can be obtained as follows: The simplified notation can be shown as follows: After determination of strain components, the constitutive equations can be separately derived for both functionally graded and functionally graded piezoelectric layers.These equations for functionally graded layer are: [Arefi and Rahimi, 2012;Gaur and Rana, 2014] And for functionally graded piezoelectric layers are: [Arefi and Rahimi, 2012] In order to complete the constitutive relations, in this step, the electric potential distribution must be defined.A three dimensional distribution of electric potential can be represented as multiplication of two functions one through radial direction () r  and another through thickness direc- tion ()  fz.
where, () fz must satisfy electric potential boundary conditions along the thickness (transverse) direction.The short-circuited boundary conditions is considered for both top and bottom piezoelectric layers.By employing a second order approximation for electric potential distribution along the z direction, we'll have: [Ebrahimi and Rastgo, 2008] Using the above equation, the electric filed components can be derived as follows: ) By substitution of electric field components from Eq. ( 9) into constitutive equations of piezoelectric layers (Eq.( 6)), we will have: In order to attain the final governing equations, the resultant of force and moments per unit width must be evaluated.
As mentioned above, the integral must be evaluated along the thickness direction.Since the plate is included two different materials, above integral must be decomposed into two integral Substituting the stress relations in terms of strain components and electric potential presents the resultant of force and moments as follows: where, mentioned coefficients , ii EA may be considered in Appendix A.
Defined components of the resultant of force and moments must satisfy equilibrium equations as follows: Substitution of resultant of force and moments in equilibrium equations yields: Electric displacement equations along the three directions are [Gaur and Rana, 2014]: Discharge equation that implies divergence of electric displacement vanishes through the piezoelectric section yields third equation as follows: [Ebrahimi and Rastgo, 2008 In this step, we can collect three essential equations for evaluation of the results of the problem.These equations include two mechanical equations and one electrical equation.
By considering the displacement fields and electric potential as follows: sin( )sin( ), sin( )sin( ), sin( )sin( ) We can evaluate the results of the problem in terms of different geometrical and material properties.Before final evaluation of the results, the distribution of material properties must be defined.

RESULTS AND DISCUSSION
Before solution of the problem, it is appropriate to define the material properties for the FG and FGP layers.For FG layer, it is assumed that the bottom of the plate is steel and top of that is ceramic.Therefore the distribution of the material properties for FG layer is (Ebrahimi and Rastgo 2008): where, () e h is thickness of elastic solid section of the plate and n is the non-homogenous index of ceramic-metal section of the plate.Figure 1 show the distribution of modulus of elasticity along the thickness direction of the plate in terms of different values of non homogenous index.
The distribution of the mechanical and electrical properties for the two FGP layers can be supposed as a power function along the thickness direction as follows (Khoshgoftar et al 2009;Arefi and Rahimi, 2010): where, i E represents the value of the all mechanical and electrical components at e zh  and e h is thickness of the piezoelectric section.Other numerical parameters are considered as:  ), it is observed that for increasing the value of m from 1.5 till 2, the buckling load increases and then with increasing the m from 2 till 3, the buckling load decreases.As another results of this research, the effect of stiffness of material can be considered on the results of the problem.This investigation can be performed by employing a dimensionless parameter such as ( c m E E ).This study can be performed for three values of (

CONCLUSION
Electromechanical stability of a functionally graded circular plate integrated with two functionally graded piezoelectric layers under radial compressive load has been studied in this paper.All mechanical and electrical properties can be varied along the thickness direction.The effect of different geometrical and material parameters has been considered on the buckling load of the paper.The obtained results in this paper can direct engineers for production of electromechanical structures in technical application.This analysis shows that emploting a functionally graded material offers various options for optimized design.Evaluation of stability and buckling load is important for application of the piezoelectric structures in different conditions.The present results is applicable for engineer in fabrication of piezoelectric structures as electromechanical elements (sensor or actuator).Some important results are expressed as follows: 1. Investigation on the effect of ratio of thickness

Figure 1 :
Figure 1: Distribution of variable modulus of elasticity of FGM along the thickness direction.

Figure 2 :
Figure 2: Buckling load of a FGP circular plate for of different values of piezoelectric thickness ( e p h h) in terms of non homogeneous index.

Figure 3 :
Figure 3: Buckling load of a FGP circular plate for different values of outer radius ( b m a  ) in terms of non homogeneous index.

Figure 4 :
Figure 4: Buckling load of a FGP circular plate for different values of modulus elasticity of ceramic ( c ) ]