The Third-Order Shear Deformation Theory for Modeling the Static Bending and Dynamic Responses of Piezoelectric Bidirectional Functionally Graded Plates

)is work is the first exploration of the static bending and dynamic response analyses of piezoelectric bidirectional functionally graded plates by combining the third-order shear deformation theory of Reddy and the finite element approach, which can numerically model mechanical relations of the structure. )e present approach and mechanical model are confirmed through the verification examples. )e geometrical and material study is conducted to evaluate the effects of the feedback coefficients, volume fraction parameter, and constraint conditions on the static and dynamic behaviors of piezoelectric bidirectional functionally graded structures, and this work presents a wide variety of static and dynamic behaviors of the plate with many interesting results. )ere are manymeanings that have not beenmentioned by any work, especially the working performance of the structure is better than that when the feedback parameter of the piezoelectric component is added, that is, the piezoelectric layer increases the working efficiency. Numerical investigations are the important basis for calculating and designing related materials and structures in technical practice.


Introduction
In recent decades, due to the science and technology are developed so rapidly, many intelligent materials appear to serve the complicated requirements in engineering practice such as functionally graded materials (FGMs), piezoelectrics, electroactive polymers, magnetostrictive, and shape memory alloys, in which functionally graded materials are one of the most important and popular structures, which are used and considered widely by scientists worldwide. Traditionally, FGMs are a remarkable type of composite structure, where the material characteristics change gradually and smoothly from one surface to the other ones [1][2][3][4], where they can be dealt with and designed for special duties and applications. A huge number of methods, which are founded by the particulate process, layer process, and melting process are employed to manufacture the FG structures [5]. In recent years, scientists have used two-dimensional functionally graded materials to significantly lessen the thermal stress phenomena in machine parts that are subjected to intent thermal loads [6][7][8][9]. Moreover, as mentioned above, hightech structures require complex technical conditions; therefore, one structure may need to integrate many kinds of smart materials, for example, piezoelectric bidirectional functionally graded structures [10][11][12]. is is also the mixed material plate investigated in this study.
Before focusing on the main purpose of this paper, some reviews of piezoelectric bidirectional functionally graded structures and related ones are discussed. For the analysis of the static bending phenomena, om et al. [10] introduced an analysis of two-directional FG plates by using the FE method (FEM) and a developed third-order shear theory (TSDT). In their investigation, the plate kinematics were expressed by using the TSDT without any shear correction coefficients. e third-order shear deformation theory also relied on the scrupulous kinematic of displacement components, which was presented to govern over other aforementioned approaches and established based on elastic formulations. e material characteristics were simulated to be changed smoothly in two directions, and their effective characteristics were calculated through the law of commixture. Tam et al. [11] firstly presented the material optimization of the acclimatization minimization of twodirectional FG structures subjected to static loads. In this approach, the distribution of the two-directional FG mechanical characteristics was represented by using the 2D nonuniform rational B-spline (NURBS) functions. In addition, 1-D NURBS fundamental functions were employed to present the accurate structural parameter and presume the unrevealed solution in the FE analysis (FEA), where design parameters were the ceramic capacity fraction coefficients at driving points by the x and z directions established by the isogeometric analysis (IGA) approach. e optimization issue was solved by a tool called the differential evolution (DE) algorithm. For vibration analyses, Jinqiang et al. [12] considered active oscillation control of the piezoelectric FG plate by employing the piezoelectricity material as an actuator, in which the Kirchhoff multilayered plate theory was employed, and the equation of motion for the piezoelectric FG plate was derived by using Rayleigh-Ritz approach and Hamilton's principle. Wankhade and Bajoria [13] studied the vibration behaviors of piezoelectric laminated composite plates with the coupled electromechanical actuation, in which the efficient FE formulation with the higher-order shear deformation (HSDT) plate theory was used to combine the influences of shear strain components in the computational model. An isoperimetric 8-noded rectangular plate element was employed based on the linear electric potential distribution by the thickness direction to derive the mechanical relations. Nam et al. [14] investigated free and forced oscillation behaviors of clever functionally graded plates, which were stiffened by adding graphene platelets (GPLs). en, the active control of FG plates with piezoelectricity layers was analyzed. To contribute numeral solutions to fundamental issues, a computational strategy with a C 0 -HSDT polygonal FE formulation (PFEM) was developed, which was suitable for simulating both thin and thick structures. And, for dynamic problems, some important works can be counted as follows. Yas et al. [15] used a semianalytical method called the state space-based differential quadrature method (SSDQM) to investigate the 3D free oscillation of a multidirectional FG piezoelectric (FGP) annular plate structure, which was rested on Pasternak's two-parameter elastic foundations. Wu et al. [16] introduced a review of 3D analytical approaches for the investigation of multilayered and FG piezoelectric plate and shell structures. e literature considered the 3D analysis of multilayered and FG piezoelectric plates with 191 references. Yongbo et al. [17] studied analytically the deflection field of piezoelectric 2D-FGM structures, where the material characteristics varied exponentially both along with the longitudinal and horizontal directions. e double equations were introduced by variation rule and presenting segmentation of variables was employed. en, in the complex space, which was consisted of the initial and dual variables, the problem could be solved by using the symplectic development. Rasool et al. [18] studied the abridged dynamic displacements of a developed piezoelectric sandwich structure subjected to the mix of thermo-electro-mechanical loads. In the considered piezoelectric sandwich structure, two active piezoelectric components were added to the surfaces of a submissive advanced plate with two thin GRNC layers and a thick porous composite core. In the GRNC sheets, the graphene component was considered to be changed by the functional dismemberment rules. By employing developed Halpin-Tsai equations, the thermal-dependent mechanical characteristics of nanocomposite layers have been defined. Reddy's third-order shear deformation plate theory (called TSDT) was used to establish the dual reigning equations, and MLS shape functions were employed in the meshless method. Wang and Huang [19] investigated analytically the dynamic responses of piezoelectric sensor components placed in an elastic foundation subjected to the electromechanical loads persuaded by piezoelectricity layers. A sensor/actuator mechanical model considering the deformation in two longitudinal and transverse directions of the piezoelectric layers was modified. e dynamic loading transfer between the piezoelectricity layers and the elastic medium was examined by using a Fourier transform approach and integral expressions in terms of the interfacial shear and normal stress components. Quang and his co-workers [20] examined the static and dynamic response of thin and moderate thickness plates based on a cell-based smoothed discrete shear gap method (CS-DSG3) using the FSDT. In which, the CS-DSG3 wasdeveloped for investigation and active oscillation control of the functionally graded plates with piezoelectric sensor and actuator components. Besides, mechanical behaviors of the rod, plate, and shell structures taking into account the piezoelectricity were presented abundantly in documents [25][26][27].
To deal with the structures made of functionally graded materials resting on elastic foundations, most recently, some works used the HSDT to simulate the static bending, the buckling, and free vibration of these structures, readers can find relevant studies in [21][22][23].
To the best of our best knowledge, there seems to be no papers dealing with the static bending and dynamic behavior of piezoelectric two-directional FG material (2D-FGM) plates; however, this is an interesting issue with a great help in applying this structure in engineering practice, which can be used to make aircraft wings, turbines, etc. erefore, in this paper, both the static and dynamic problems of this mixed structure are carried out by using the FEM based on the TSDT of Reddy. is is a computational theory that does not need any shear correction coefficients but still ensures the accuracy of the simulation. is paper also shows that when the feedback coefficients of the piezoelectric effect are added, the performance of the plate increases significantly, in the other words, the vibration reduction is better, and 2 Advances in Materials Science and Engineering therefore, increasing the efficiency of use against structures with variable mechanical properties are integrated with piezoelectricity layers. e rest of this work is structured into the following four sections. Section 2 briefly introduces piezoelectric bidirectional FG plates. Section 3 presents a clear FE formulation for piezoelectricity 2D-FGM plate structures. Section 4 shows the accuracy problems and numerical explorations of both static bending and dynamic response analyses of the plate. Some new explorations are summed up in Section 5.

Piezoelectric Bidirectional Functionally
Graded Plates is paper focuses on the static and dynamic behaviors of the piezoelectric 2D-FGM plate with the core layer made from three different materials, which are represented as 1, 2, and 3, respectively. Material characteristics including Poisson's ratio, Young's modulus, and mass density are ] i , E i , and ρ i (i � 1-3). Two outer surfaces are the piezoelectricity actuator and sensor layers in Figure 1.
e mechanical characteristics of the core layer are contingent on the volume fraction exponents of three different materials as follows [10]: where V i (i � 1-3) are the volume fraction of materials, and [10] in which, n and q ≥ 0 are the volume fraction indices.

Finite Element Formulation of Piezoelectricity 2D-FGM Plates
Herein, the TSDT of Reddy is employed; thus, the displacement components at one point in the structure are expressed as [24] u(x, y, z) with u 0 , v 0 , and w 0 are the displacement components at (x, y, z � 0) and φ x an d φ y are the transversal normal rotations in xz and yz surfaces. e strain parts are calculated as Equation (4) is then shortened in the matrix expression as Advances in Materials Science and Engineering e linear piezoelectric fastening between the electric field and the elastic field can be defined by the converse piezoelectric and direct equations and the constitutive and electric displacement relations [21,22,[24][25][26][27][28][29][30].
For the core 2D-FGM layer, For the two piezoelectric layers, For the core 2D-FGM layer, For the two surface piezoelectric layers, So, the stiffness coefficients D b and D s vary in both x-and y-directions, and this is completely different from the 1D-FGM plate.
−E v is the component of the electricity, where it is calculated by the transformation of the potential as [23]    Advances in Materials Science and Engineering In this work, it is considered that the electric field is applied only to the structure in the depth direction, and the symbol E z is the electric component, which is calculated as -D is the displacement vector due to the electricity; -e represents the stress coefficients related to the piezoelectricity, which has the following specific expression: -p is the component (tensor) of the dielectricity at unchanged mechanistic strain, which has the description as follows: e internal load can be presented as in which is work employs the plate element of four-node, and each one contains five DOFs: Advances in Materials Science and Engineering 5 where N j is the interpolation functions (Lagrangian type) for the four-node plate element and ξ and η are the values in the natural coordinate system. According to the previous interpolation method, the strain fields are shortened as en, the electric field vectors (in the sensor and actuator layers) are defined by two degrees of freedom ϕ a and ϕ s as where B a ϕ � 0, 0, where h s and h a are the depths of the sensor and actuator layers. e whole strain energy of the piezoelectric FG plate can be computed as where F is the facade load, and expression (22) Advances in Materials Science and Engineering 7 where the stiffness matrices are defined as where Q b and Q a are the electric potentials, which act on the sensor and actuator layers, respectively. e expression of the kinetic energy of the plate element is expressed as where ρ, ρ b , and ρ a are, respectively, the densities of the core, sensor, and actuator layers, N is Lagrangian shape function matrices, and L and M e are calculated as en, the weak form type Galerkin of the controlling equations of the plate can be formed by employing Hamilton's principle, where it can be expressed as in which Substituting the expressions of the kinetic function into equation (28), then group them individually into the displacement and potential variables, and the resulting equation is expressed as From the second and third equations in (29), ϕ ae and ϕ se can be obtained; then, adding them into the first expression, the forced oscillation equation of the plate element in the matrix form is obtained as When the external load acts on the structure, the structure is going to be bent; thus, the voltage will be varied; then, the voltage of the sensor will be passed to the controlling center to restrain the acting voltage return to the actuator one to derive an electric load that proceeds opposing to the way, where the structure is bending, getting back the structure t1o its balance before the deformation phenomenon. e operation is expressed as presented in Figure 2. Now, consider G v as a velocity feedback coefficient and G d as a displacement reaction coefficient of the controlling station, and these reaction coefficients deputize the relation in the middle of the input voltage ϕ a and the output voltage in the piezoelectric sensor layer ϕ s : Herein, by neglecting the aspect of the concerned Q a and Q s , from the third-expression of equation (29), the voltage vector of the piezoelectric sensor component is calculated as follows: is can be explained that, as the system is vibrated, the voltage engendered in the piezoelectric sensor component due to the strain phenomenon of the structure is about to engender an initial indication to the controlling station, which is going to affect directly on the piezoelectric actuator layer. erefore, the stress and strain components due to the piezoelectricity will generate a voltage indication applied back to it, which controls the strain and stress in the desired way.
By adding equations (31) and (32) into the second and third equations in (29), we obtain the following expression: Finally, by adding equation (34) into equation (30), the forced oscillation equation of plate element is obtained as in which and C d e is the damping matrix caused by the influence of the feedback parameter, which is defined as Neglecting the damping coefficient generated by the structure, the forced oscillation equation of the element is expressed as where C e is the Rayleigh structural dampness matrix, where it is computed by the mass and stiffness matrices: By operating the process of coming together with the element vectors and matrices and deleting the fixed DOFs, the forced oscillation equation of the piezoelectric structure is obtained as follows: . As the volume fraction exponent n � 2, G d � 5, and the deflection of the plate at y � b/2 is presented in Figure 3 with different mesh sizes, and compared with those of Nguyen-Quang et al. [20], one can see that the mesh size is increased (the number of elements increases) and the deflection curve of the plate converges to the result shown in document [20]. It can be seen that, with the mesh size 8 × 8, the numerical data of this work meets a good agreement. erefore, for all following explorations, this mesh size will be employed. To be more convenient, the nondimensional deflection is normalized as follows:

Computed Results of Static
(1) Effect of the Volume Fraction Indices n and q. By changing the values of volume fraction indices n and q in the range of 0 to 10, the nondimensional deflection w * is plotted in Figure 4. When n varies from 0 to 2, if the value of q increases, the nondimensional deflection w * will increase. However, when the value of n is higher than 2, the increase of q will guide to a decrease in the nondimensional deflection w * . is is because changing both n Advances in Materials Science and Engineering and q results in a change in the substance characteristics; therefore, the mechanical behaviors of the structure are also changed.
(2) Effect of the Feedback Coefficient G d . Figure 5 presents the nondimensional defection of the plate at y � b/2 with the value of the feedback coefficient G d varies in a range of 0 to 10 for the SSSS plate (n � 10 and q � 0.5). When increasing G d , the maximum deflection of the plate decreases, and the deflection line is no longer symmetrical through the x � a/2. is is due to the fact that the material also varies in the Ox direction of the plate. Figures 5-8 present the deflection of the plate at y � b/2 for the CFFF plate. As the reaction coefficient G d increases, the nondimensional maximum deviation of the structure degrades. is can be illuminated that, due to the increase of G d , the stiffness matrix is increased by an amount as presented in equation (36), so the structure becomes stiffer and causes the deflection to decrease. e characteristics of the piezoelectric layers are presented in this section. e load is changed over time by the following formula:

Numerical Results and
with the load amplitude P 0 � 10 4 N/m 2 , and the function F (t) is expressed as where t 1 � 2 ms. e comparative results of dynamic behaviors of the center point of the piezoelectric FG plate are shown in Figure 9 corresponding to the case of G v � 0 and G d � 0 and compared with those of the CS-DSG3 [20]. Readers can see that the present theory and mechanical model are verified.   (1) Effect of the Feedback Coefficient G v . Now, to examine the influence of G V on the proficiency and to eliminate the vibration of the structure, the plate with n � 2 and q � 1 is considered. Nondimensional deflexion w * of the center location of the outermost edge of the plate (x � a, y � b/2) is presented in Figures 10 and 11. e results indicate that As the value of the response coefficient G v increases, nondimensional deflection w * will progress to 0 faster. is  means that the proficiency to eliminate the vibration of the plate is better. At the same time, the maximum value of the nondimensional deflection w * also decreases, and this shows that the plate becomes stronger. is is because when G v is increased, the stiffness matrix of the structure is increased according to equation (37), causing the vibration energy of the plate to be dissipated; thus, the amplitude of the oscillation of the plate decreases.

Numerical Results and
As the value of the feedback coefficient G V becomes greater, the culmination value of nondimensional deflection w * decreases gradually.
(2) Influence of the Depth of the Piezoelectric Component. e changing the piezoelectric layer thickness h s will affect the overall stiffness and mass of the plate; therefore, it will significantly affect the response of the piezoelectric 2D-FGM   plate. Figures 12-14 show the change of the nondimensional deflection w * over time as the depth of the piezoelectric layer h s is increased, h s /h � 0.1-0.5, respectively. Some discussions are given as follows: (1) As the depth of the piezoelectric layer h s increases, maximum deflection w * decreases. At the same time, the shape of the deflection curve over time is also different. is is because both stiffness, mass, and structural damping matrices of the structure are changed, so both the maximum value and the responsive shape of the structure also change.
(2) When h s /h � 0.1 and 0.2, changing the value of G d from 0-5, it can be seen that the maximum value of the deflection w * does not change much, but when G d increases, nondimensional deflection w * will reach 0 faster. In addition, Figures 12-14 show that, at the moment of the free oscillation, w * will

Conclusions
is work carries out static bending and dynamic behavior analysis of piezoelectric 2D-FG plates by using the TSDT of Reddy and the FEM. e paper confirms the aptness of the present approach through comparison with the previously published data. rough the parameter study, some conclusions are highlighted as follows: (1) As the value of the feedback coefficient G d becomes greater, the culmination deflexion of the structure decreases (2) As the value of the feedback coefficient G v gets larger, the culmination deflexion of the structure decreases in the forced oscillation problem, and the vibration of the plate reduces faster, which is the advantage of the plate when the piezoelectric layer is integrated (3) When the depth of the piezoelectricity component is boosted, the culmination deflection of the structure decreases, and the shape that responds to the deflection of the plate over time also changes.

Data Availability
e data used to support the findings of this study are included within the article.