Symplectic Method for the Thin Piezoelectric Plates

: The symplectic method for a thin piezoelectric plate problem is developed. The Hamiltonian canonical equation of thin piezoelectric plate is given by using the variational principle. By applying the separation of variables method, we can obtain symplectic orthogonal eigensolutions. As an application, the problem of a thin piezoelectric plate with full edges simply supported under a uniformly distributed load is discussed, and analytical solutions of the deﬂection and potential of a piezoelectric thin plate are obtained. A numerical example shows that the solutions converge very rapidly. The advantage of this method is that it does not need to assume the predetermined function in advance, so it has better universality. It may also be applied to the problem of thin piezoelectric plate buckling and vibrating.


Introduction
Recent years have seen an increase in the use of piezoelectric materials in structural shape control, vibration and noise control, shape detection, and other fields due to their good electromechanical coupling characteristics [1][2][3][4][5]. At present, there are many models to study the piezoelectric effect. Dunn and Wienecke [6] applied Green's function to transversely isotropic piezoelectric solids. Hill [7] applied the boundary element method to the problems of 3D piezoelectricity. Ding et al. [8] obtained some exact solutions for the 3D piezoelectric problem. The above research are based on the Lagrange system.
The symplectic method for solving the elastic theory was proposed by Zhong [9,10]. The characteristic of this method is that it does not assume the trial function in advance but introduces the problem into the Hamiltonian system and uses the separated variable method to solve the elastic problem under specific boundary conditions. This method differs from the traditional semi-inverse method, it does not require a trial function, so it is more universal. The symplectic method can easily be adapted to various thin plate problems, thick plate problems, or a combination of vibration and buckling. Therefore, many researchers use the symplectic method to study problems in mechanics and engineering science. The 2D and 3D symplectic structures of transversely isotropic piezoelectric media have been derived by Gu [11][12][13][14]. Xu et al. [15] studied the problem of interfacial cracks in non-ideal interfacial magnetoelectroelastic bimaterials. Some exact solutions for the bending and vibration of rectangular plates were also obtained by Li and Zhong [16][17][18]. In a study by Zhang and Deng et al. [19], non-local Timoshenko beams with gradational materials were studied for free vibration. Xu and Deng et al. [20] analyzed the free vibration and buckling of natural-fiber-reinforced composite plates with partial or internal cracks supported by corners. Leung and Mao [21] presented the symplectic Galerkin method for non-linear vibrations of beams and plates. Jia et al. [22] calculated the free vibration of orthotropic and isotropic cylindrical shells whose thicknesses were both uniform and stepped under general boundary conditions. The above reflects the wide application of the symplectic method in various fields. Therefore, in the present study we have further developed the symplectic method for thin piezoelectric plates in the Hamiltonian system. The basic equation of the thin piezoelectric plate is transformed into a Hamiltonian canonical equation. By applying the separation of variables method, we can obtain symplectic orthogonal eigensolutions. In this paper, the analytical solutions for a thin piezoelectric plate with full edges simply supported are obtained. Figure 1 shows the coordinate system of a thin piezoelectric plate subjected to a uniformly distributed load q and applied voltage V, where 0 ≤ x ≤ a, 0 ≤ y ≤ b and h is the thickness of the plate.

Hamilton Canonical Equations
vibration of orthotropic and isotropic cylindrical shells whose thicknesses w form and stepped under general boundary conditions. The above reflects the wide application of the symplectic method in v Therefore, in the present study we have further developed the symplectic m piezoelectric plates in the Hamiltonian system. The basic equation of the thin plate is transformed into a Hamiltonian canonical equation. By applying the variables method, we can obtain symplectic orthogonal eigensolutions. In t analytical solutions for a thin piezoelectric plate with full edges simply sup tained.  The constitutive equations of a transverse-isotropic piezoelectric plate [23] ( ) 11    The constitutive equations of a transverse-isotropic piezoelectric plate are as follows [23] σ xx = C 11 ε xx + C 12 ε yy + C 13 ε zz − e 31 E z , σ yy = C 12 ε xx + C 11 ε yy + C 13 ε zz − e 31 E z , σ zz = C 13 ε xx + C 13 ε yy + C 33 ε zz − e 31 E z , σ xy = 2C 66 ε xy ,

Hamilton Canonical Equations
The strain-displacement relations are given as below: in which C 66 = (C 11 − C 12 )/2. σ ij , ε ij , C ij are the stress, strain, and elastic constants, respectively. D i , E i , e ij , and κ ij , respectively represent the electric displacements, the electric field intensity, and the piezoelectric and the dielectric constants. According to Kirchhoff's hypothesis [24], i.e., ∂w/∂z = 0, ε xz = 0, ε yz = 0, the displacement components u, v and w in the x, y and z directions are described as follows: Substituting Equation (5) into the Equation (2), we have: Given that the thin piezoelectric plate has a large in-plane size to thickness ratio and it was subjected to the potential across its thickness, the electric field component in the length direction can be ignored [25,26]. Therefore, it is assumed that the electric field exists only along the z direction and the electric field intensity vector can be expressed as follows: where φ is electric potential and D x and D y are approximately assumed to be zero. Thus, the electric boundary conditions are as follows: Hence, the stress and electric potential displacement of the thin piezoelectric plate can be deduced from Equation (1): According to Equations (4) and (7)-(9), the expressions for the electric potential and electric field can be obtained as: According to the principle of virtual work, the expressions of the axial force and bending moments can be obtained as [26]: where M ij , Q i are the bending moments and shear forces, respectively. N ij are the internal forces (axial forces) caused by the electrical coupling effect. The bending moments, shear forces, and axial forces are related to stresses by: where By Equations (11) and (12), the bending equation of the piezoelectric thin plate is as follows: The plate functional is [27]: where Set the derivative of w over y be θ ∂w ∂y where the dot indicates ∂/∂y. In combination with Equation (12) and introducing the Lagrange multiplier by V y , Equation (14) is modified as follows: Set T = −V y − e 31 Vθ, by Equation (17), we have: (17) can be expressed in the following matrix form: where It is obvious that , is the symplectic matrix and I 2 is the 2 × 2 unit matrix. Therefore, H is a Hamiltonian operator matrix, and Equation (18) is the Hamilton canonical equations for the thin piezoelectric plate.
In order to the separate variables in Equation (19), we let: where X(x) = w(x), θ(x), T(x), M yy (x) T . Substituting Equation (20) into Equation (19) gives where µ is the eigenvalue and X(x) is the corresponding eigenvector.
the eigenvalue equation is derived by extending the determinant which has roots: . The eigenvector X(x) of H satisfies: where A i , B i , F i , G i (i = 1, 2, 3, 4) are undetermined constants, and these constants are usually not independent. Combined with the boundary conditions, the above constants can be determined, then the analytical solutions can be obtained.

Symplectic Solutions of a Thin Piezoelectric Plate with Full Edges Simply Supported
In order to illustrate the application of the symplectic method in the thin piezoelectric plate problem, consider a thin piezoelectric plate whose full edges are simply supported as is shown in Figure 2.

Symplectic Solutions of a Thin Piezoelectric Plate with Full Edges Sim ported
In order to illustrate the application of the symplectic method in the thin plate problem, consider a thin piezoelectric plate whose full edges are simp as is shown in Figure 2. The full edges simply supported satisfies the boundary conditions:

Analytical Solutions of Thin Piezoelectric Plate
By Equations (32) and (33), we have the solution of the state vector Z, which is: The constants C 1 , C 2 , C 3 , C 4 can be solved by the boundary conditions (27) as follows: Finally, the bending deflection of the thin piezoelectric plate with full edges simply supported under a uniformly distributed load is: where By Equation (10), the electric potential can be expressed as: where h n (y) = 1 − cosh(ξy) cosh(ξb/2) .
When an applied voltage is not considered, Equations (44) and (46) can be rewritten as: Especially, it is worth mentioning that when the piezoelectric material factor is not considered, Equation (48) is exactly the same as the classical plate bending solution [24].

Numerical Results
In this section, the numerical results are discussed in detail to validate the effectiveness of the proposed method.
The dimensionless bending deflections w (qa 4 /D 11 ) at the midpoint (a/2, 0) of the classical thin plate are calculated as the special case of the present results, shown in Table 1. A comparison with the known results of ref. [24] displays excellent agreement. PZT-4 material is selected as an example and the material properties are C 11 = 139 GPa, C 12 = 77.8 GPa, e 31 = −6.98 N/Vm, κ 33 = 5.47 × 10 −9 F/m [28]. We assume that the length of the thin plate is a = 1 m, the width is b = 1 m, the thickness is h = 20 mm, the uniformly distributed load is q = −2 MPa, V = 20 V. Table 2 presents w of thin piezoelectric plates with full edges simply supported n is the expanded term. The results show that the solutions converge very rapidly. The bending deflection diagram of a thin piezoelectric plate with full edges simply supported under a uniformly distributed load is shown in Figure 3.
In this section, the numerical results are discussed in detail to validate the effectiveness of the proposed method.
The dimensionless bending deflections 4 11 ( / ) w qa D at the midpoint ( ) / 2,0 a of the classical thin plate are calculated as the special case of the present results, shown in Table 1. A comparison with the known results of ref. [24] displays excellent agreement.   Table 2 presents w of thin piezoelectric plates with full edges simply supported n is the expanded term. The results show that the solutions converge very rapidly. The bending deflection diagram of a thin piezoelectric plate with full edges simply supported under a uniformly distributed load is shown in Figure 3.

Conclusions
The symplectic method for the thin piezoelectric plate problem is established in this paper. The analytical solutions are obtained for the bending deflection and the electric potential. The accuracy, validity, and reliability of the proposed symplectic method are verified by comparison studies with respect to the published literature. Compared with

Conclusions
The symplectic method for the thin piezoelectric plate problem is established in this paper. The analytical solutions are obtained for the bending deflection and the electric potential. The accuracy, validity, and reliability of the proposed symplectic method are verified by comparison studies with respect to the published literature. Compared with the traditional methods, the implementation of the symplectic method is simple. It can also be applied to other problems with piezoelectric materials.