Document Type: Research Paper
Authors
^{1} Department of Mechanical Engineering, Arak University, Arak, 3815688349, Iran
^{2} Institute of Nanosciences & Nanotechnolgy, Arak University, Arak, 3815688349, Iran
Abstract
Keywords

Mechanics of Advanced Composite Structures 6 (2019) 117 – 129


Semnan University 
Mechanics of Advanced Composite Structures journal homepage: http://MACS.journals.semnan.ac.ir 
Analytical Approach for Thermoelectromechanical Vibration of Piezoelectric Nanoplates Resting on Elastic Foundations based on Nonlocal Theory
K. Khorshidi ^{a,b}^{*}, M. Karimi ^{a}
^{a }Department of Mechanical Engineering, Arak University, Arak, 3815688349, Iran
^{b }Institute of Nanosciences & Nanotechnolgy, Arak University, Arak, 3815688349, Iran
Paper INFO 

ABSTRACT 
Paper history: Received 20180729 Received in revised form 20181229 Accepted 20190217 
In the present work, thermoelectro vibration of the piezoelectric nanoplates resting on the elastic foundations using nonlocal elasticity theory are considered. Inplane and transverse displacements of the nanoplate have been approximated by six different modified shear deformation plate theories considering transverse shear deformation effects and rotary inertia. Moreover, two new distributions of transverse shear stress along the thickness of the nanoplate were introduced for the first time. The equations of motion were derived by implementing Hamilton’s principle and solved using analytical method for various boundary conditions including SSSS, CSSS, CSCS, CCSS and CCCC. Based on a comparison with the previously published results, the accuracy of the results was confirmed. Finally, the effects of different parameters such as boundary conditions, variations of the thickness to length ratio, aspect ratio, increasing temperature, external voltage, foundation coefficients and length scale on the natural frequency of the plate were shown and discussed in details. 



Keywords: Analytical method Elastic foundation Piezoelectric nanoplates Nonlocal theory 


© 2019 Published by Semnan University Press. All rights reserved. 
Piezoelectric materials such as ZnO, ZnS, PZT, etc. are widely used as actuators and sensors due to their electromechanical coupling effects. In recent years, piezoelectric nanostructures such as nanohelices, nanowires, nanorings, etc., have been in the spotlight thanks for desires in having small systems in comparison with the macroscaled system and are used as essential components in different industries as nanogenerators, chemical sensors, lightemitting diodes, etc. [13]. Structures with the small dimensions between a few nanometers and 100 nanometers (nanostructures) do not treat the same as structures with macro scale. In fact, the behavior of the structure is dependent on the length scale and, conventional continuum model cannot be used on such scale. Unlike the classical theories which do not consider size effect, other developed theories such as strain gradient [4], couple stress [5] and nonlocal theory [6] are sensitive to length scale in material behavior. Nonlocal theory of Eringen as one of the most applied theories could capture the nonlocal effects and is beneficial for nanostructures. In local classical theories, there is a onetoone relation between stress and strain of each point i.e. the state of stress in every point is a function of the strain at that point; while in nonlocal theory, the stress of a reference point is dependent on strains of the whole domain.
Nanostructures held attentions due to their unique properties. Recently, nanostructures such as nanobeams and nanoplates are widely used in nanoelectromechanical (NEM) devices. Dynamic behavior of nanoplates based on classic plate theory (CPT) has been discussed using Eringen model [79]. Pradhan and Phadikar [10] utilized the nonlocal theory based on firstorder shear deformation theory (FSDT) so as to predict the vibrational behavior of nanoplates. Aghababaei and Reddy [11] employed thirdorder shear deformation plate theory (TSDT) for bending and free vibration of nanoplates. Khorshidi and Asgari [12] studied the free vibration analysis of functionally graded rectangular nanoplates based on nonlocal exponential shear deformation theory. Alibeigloo [13] used the threedimensional theory of elasticity for analyzing the free vibration of nanoplates employing nonlocal continuum mechanics. Extensive studies on piezoelectric structures can be found in open literature [1416]. Khorshidi and Pagoli [17] presented an analytical solution for vibrating circular plates coupled with piezoelectric layers. Dynamic response of piezoelectric nanobeam according to Timoshenko beam theory and nonlocal theory has been evaluated by Ke and Wang [18]. Hung and Yu illustrated the efficacy of the surface piezoelectricity on the electromechanical response of piezoelectric ring [19]. Jandaghian and Rahmani [20] studied the sizedependent vibration analysis of functionally graded piezoelectric plates subjected to thermoeletromechanical loading with Kirchhoff plate theory. Nikoei and Hassani [21] employed an isogeometric analysis approach for free vibration of laminated composite plates integrated with piezoelectric using ReissnerMindlin theory.
The elastic foundation was used to model a rather soft material in contact with the plate surface. Elastic mediums could be simulated with various models. First of all, Winkler presented the simplest foundation model because He did not consider the interaction between lateral springs. For improving interaction between springs, different foundation models were proposed later such as Pasternak [22], Vlasov [23] and FilonenkoBorodich [24] models. Pasternak is a useful twoparameter foundation model which captures the interaction of adjacent springs by joining ends of springs to a shear layer. Many studies are dedicated to the dynamic behavior of plates resting Pasternak foundation [25, 26]. Moradi et al. [27, 28] analyzed the nanocomposite plates and sandwich plates reinforced by wavy carbon nanotubes resting on elastic foundation. They used a meshfree method and firstorder shear deformation theory in their analysis.
Classical plate theory was proposed by Kirchhoff [29, 30]. He assumed that straight lines normal to the midplane before deformation, remain straight and normal to the midplane after deformation. Since Kirchhoff theory ignores the transverse shear deformation effects, it is not proper for moderately thick plates. Therefore, it is devoted to the thin plates. Firstorder shear deformation theory [31, 32] incorporates the shear deformation effects with a constant transverse shear deformation distribution along the thickness of the plate. Thus, it violates the stressfree conditions at the bottom and top of the plate and needs a shear correction factor to compensate this error. In order to get more accurate results and avoid using the shear correction factor, higherorder shear deformation theories (HSDT) have been developed. Reddy [33] employed a parabolic transverse shear stress distribution along the thickness of the plate. His model did not need shear correction factor because of satisfying free stress conditions at the bottom and top of the plate. Various distributions of transverse shear stress through the thickness of the plate can be found in the literature [3437]. For example, Sayyad and Ghugal [38, 39] presented exponential and trigonometric shear deformation theories for bending and free vibration analysis of moderately thick plates. Different methods such as analytical, numerical and semianalytical methods have been utilized for solving the equations of motion according to literature. Ke et al. [40] investigated the vibration of the Mindlin piezoelectric nanoplate for different boundary conditions using the differential quadrature method. Khorshidi et al. [41] presented the Navier solution for vibrating piezoelectric nanoplate. They considered only fully simplysupported boundary condition in their analytical solution.
In the present study, the vibration analysis of piezoelectric nanoplates under external voltage and temperature resting on a Pasternak foundation for five different boundary conditions has been illustrated. Dynamic analysis was carried out for six types of different modified shear deformation theory and nonlocal elasticity theory. Governing equations were reduced to a single equation using a simple method and analytical solution for this single equation was obtained for five different boundary conditions. In the result section, efficacy of different variables such as foundation coefficients, nonlocal parameter, aspect ratio, temperature rising and external voltage on fundamental frequency was discussed in detail.
As previously mentioned, the stress at each point in the nonlocal theory is related to the strain field in whole domain of the body. This theory has been developed by Zenkur [42] for capturing thermal effects. Based on this theory, stress and electric displacement at a reference point are expressed as:
(1) 

(2) 


(3) 
(4) 
In the above equations, and (i, j, k= 1,2, 3) denote the components of stress field, electric displacement, strain, displacement and electric field respectively. _{ }are elastic constants; and are piezoelectric and dielectric constants; and are thermal moduli and pyroelectric constants respectively. represents the Kernel function in nonlocal theory such that is Euclidean distance and τ=e_{0}a/L. where e_{0}, a and L denote the material constant, internal and external characteristic length respectively. Since the integral form of nonlocal elasticity theory is complicated, Eringen [6] obtained a differential form of nonlocal theory in the following equations.

(5) 
(6) 
where represent the Laplace operator. Under the assumption of plane stress conditions, following constitutive relations would be derived by expanding (5) and (6):
(7) 


(8) 
where

(9) 
Generally, the modified shear deformation theories are displacementbased ones which are recognized to be strong methods for upgrading the accuracy of the results [43]. Against the classical plate theory, modified shear deformation theories consider both rotary inertia and shear deformation effects. Fig. 1 not only shows the distribution of transverse shear stress along the thickness of the plate, but also verifies that these theories satisfy the tangential traction free boundary condition. Thus, a shear correction factor is not required in such theories.
According to former studies, modified shear deformation theories have been developed using hyperbolic, exponential, polynomial and trigonometric functions along the thickness of the plate. Based on the modified shear deformation theory, the displacement field could be expressed in the following form:
Fig. 1. Distribution of transverse shear stress through thickness of plate

(10) 
where u_{1}, u_{2} and u_{3} represent displacements of an arbitrary point along x, yand z axis, respectively; w is the outplane displacement of the midplane in the nanoplate (on the zdirection); ξ and ψ are rotation functions in xoz and yoz plane; u and v denote inplane displacements of midplane surface of the plate along xaxis and yaxis, respectively. Based on this theory, displacements u_{1} and u_{2} include two parts; First part is the same as classical plate theory; second part which is considered to capture shear deformations. Depending on choosing f_{i}(z) (i=1,6), these shear deformation theories are varied according to Table 1. In fact, the function f_{i}(z) specifies the distribution of transverse shear stress along thickness coordinate. The last two functions f(z) in Table 1 are suggested for the first time in this article. It is worth mentioning that the inplane displacement components could be neglected (u and v=0) in the analysis of transverse vibration due to the homogeneity of the structure [12].
Consider a rectangular peizoelectric nanoplate with length L_{1 }along x axis (0<x<L_{1}), width L_{2} along y axis (0<y< L_{2}) and thickness h in zdirection (h/2<z<h/2) as shown in Fig. 2. Nanoplate is subjected to external voltage V_{0}, temperature rising ∆T and is rested on an elastic foundation. Wang [3] approximated the electric potential as a combination of cosine and linear variation in order to satisfy the Maxwell equation:

(11) 
where γ=π/h and is electric potential in midplane of nanoplate. By neglecting the inplane displacements and substituting displacement field (10) in (4) and electric potential (11) in (3), linear strains and electric field could be found as below:

(12) 
(13) 
Strain energy (U), kinetic energy (T) and work done by external forces (W_{f}) in piezoelectric nanoplate could be obtained as follow:
(14) 


(15) 

(16) 
Table 1. Different shear deformation theories
f_{i}(z) 
Theory 

exponential [38] 
f_{1}(z) 

trigonometric [39] 
f_{2}(z) 

hyperbolic [35] 
f_{3}(z) 

parabolic [36] 
f_{4}(z) 

1^{st} suggestion 
f_{5}(z) 

2^{nd} suggestion 
f_{6}(z) 
Fig. 2. Geometry of piezoelectric nanoplate subjected to the thermoelectro loads resting on elastic foundation 
where A is the midplane of nanoplate (at z=0); k_{s} and k_{p} are spring stiffness and Pasternak’s modulus respectively. (F_{Px}, F_{Py}), (F_{Tx}, F_{Ty}) and (F_{Ex}, F_{Ey}) are normal forces created by axial force P, temperature rising ΔT and external electric voltage V_{0} given as:
(17) 

(18) 

(19) 
Hence, Hamilton’s principle would be transformed to:
(20)
By incorporating the Eqs. (14) through (19) into Eq. (20) and integrating by parts, the following equations would be obtained [41]:
(21) 

(22) 

(23) 

(24) 
where the operator L_{i} is defined by
(25) 

(26) 


(27) 
(28) 


(29) 

(30) 

(31) 

(32) 
(33) 

(34) 
Moreover, A_{i} and B_{i} are given as:
(35) 


(36) 

(37) 
(38) 

(39) 

(40) 

(41) 

(42) 

(43) 
Then, and can be obtained in terms of w according to Cramer method in solving systems of linear equations into Eqs. (21) through (23).
(44) 
where the operator is defined as:
(45) 

(46) 

(47) 

(48) 
Finally, by taking the operator from both sides of Eq. (24) and using Eqs. (44) through (48), transverse displacement equation of piezoelectric nanoplate can be obtained as follows:
(49) 
Analytical solution for Eq. (49) can be applied to different boundary conditions according to Table 2. The solution that satisfies appropriate boundary conditions can be expressed in the following form:
(50) 
where ω is the natural frequency associated with (mth, nth) mode. Eigenfunctions X_{m}(x) and Y_{n}(y) which are listed in Table 2 are chosen so that, satisfy at least geometric boundary conditions. By substituting Eq. (50) into Eq. (49) and using orthogonal properties of trigonometric functions, partial differential equation (49) would be converted to an algebraic equation and the solution will be found.
In this section, numerical results for vibration analysis of piezoelectric nanoplate subjected to thermoelectro loads resting on the elastic foundation under different boundary conditions including, SSSS, CSSS, CSCS, CSCS, CCSS, CCCC are illustrated. For all calculations, length of nanoplate (L_{1}) is considered 50nm and thickness (h) is taken as 5nm, otherwise they are specified. The results obtained from the modified shear deformation theories are so similar to each other and it is not possible to show this slight difference between them through the figures; consequently, it is necessary to choose one theory for interpreting the results.
Accordingly, all figures were plotted based on f_{6}(z). Material properties of PZT4 used in this section are given as [40]:
Moreover, the following parameters were used for illustrating results:
Table 2. Admissible functions for various boundary conditions [44]

Boundary conditions 

functions X(x) and Y(y) 



SSSS 

CSSS 

CSCS 

CCSS 

CCCC 
For verification of this model, results obtained from simply supported square nanoplate have been compared with corresponding ones in the open literature. Table 3 indicates an excellent agreement between the present work and the former studies results. According to this table, there is a bit different among the results obtained from various shear deformation theories. These differences are due to the fact that, function f(z) have different expansions through the thickness in various theories. It is worth to mention that every extra power in the expansion of function f(z) through the thickness of the structure includes additional unknown variables in those theories. Additionally, physical interpretation of these unknown variables are difficult [46]. Thus, it is better to use such distributions that are simpler with acceptable accuracy. Although two new proposed theories are simpler than other modified shear deformation theory, they are nearly identical in accuracy. The expansion of the functions f(z) through the thickness of the structure for various distribution mentioned in Table 1 are given as follow:
….

Fundamental frequencies of the simply supported piezoelectric nanoplate are presented in Table 4 based on different theories for various values of aspect ratios and thickness ratios. According to this table, the results obtained based on f_{5}(z) are equal to the results obtained by f_{4}(z) up to seven decimal places due to the fact that both of these distributions are thirdorder polynomial along the thickness of the plate. Moreover, it is seen that the fundamental frequencies of the structure based on various f_{i}(z) are a little different and can be sorted in the order of f_{6}(z)> f_{1}(z)> f_{2}(z)> f_{3}(z)> f_{4}(z)= f_{5}(z). In fact, by decreasing the thickness of the structure, the difference among these theories declines since shear deformation effects could be neglected in the thin plates. The results in this table have been extracted for the values of µ=0.1 and ΔT=V_{0}=K_{w}=K_{p}=0.
Fig. 3 shows the variation of the first four dimensionless frequencies of simplysupported nanoplate with the nonlocal parameter. This figure depicts that nonlocal parameter is more highlighted in higher modes because wavelength gets smaller by increasing the number of modes. Thus, the nonlocal parameter is more significant for smaller wavelengths. If nonlocal model views as atoms connected to each other by springs, in the case of local elasticity, the stiffness of these springs takes an infinite value. Thus, frequencies and stiffness of structure decreased with the increase in the nonlocal parameter, as Fig. 4 emphasizes on this matter.
Fig. 3. Variation of first four dimensionless frequencies of SSSS nanoplate with nonlocal parameter
Table 3. Comparison of fundamental dimensionless frequency of SSSS square piezoelectric nanoplate with varying nonlocal parameter
e_{0}a/L_{1} 
exponential 
trigonometric 
hyperbolic 
parabolic 
1^{st} suggestion 
2^{nd} suggestion 
[40] 
[45] 
0 
0.60590 
0.60582 
0.60580 
0.60580 
0.60580 
0.60593 
0.6068 
0.6290 
0.1 
0.55371 
0.55364 
0.55361 
0.55361 
0.55361 
0.55374 
0.5545 
0.5748 
0.2 
0.45292 
0.45286 
0.45285 
0.45285 
0.45285 
0.45295 
0.4536 
0.4702 
0.3 
0.36362 
0.36357 
0.36356 
0.36356 
0.36356 
0.36364 
0.3641 
0.3775 
0.4 
0.29712 
0.29709 
0.29708 
0.29707 
0.29707 
0.29714 
0.2976 
0.3085 
0.5 
0.24871 
0.24868 
0.24867 
0.24867 
0.24867 
0.24872 
0.2491 
0.2582 
Table 4. Comparison fundamental frequency (GHz) of SSSS piezoelectric nanoplate with various theories 

S 
f_{1}(z) 
f_{2}(z) 
f_{3}(z) 
f_{4}(z) 
f_{5}(z) 
f_{6}(z) 

0.01 
0.5 
0.3993743 
0.3993740 
0.3993739 
0.3993739 
0.3993739 
0.3993745 
1 
0.6188323 
0.6188315 
0.6188312 
0.6188312 
0.6188312 
0.6188327 

2.5 
1.8720008 
1.8719913 
1.8719888 
1.8719887 
1.8719887 
1.8720048 

0.05 
0.5 
1.9852550 
1.9852150 
1.9852033 
1.9852030 
1.9852030 
1.9852717 
1 
3.0644812 
3.0643795 
3.0643506 
3.0643500 
3.0643500 
3.0645238 

2.5 
9.0365220 
9.0354308 
9.0351154 
9.0351089 
9.0351089 
9.0369784 

0.1 
0.5 
3.9011997 
3.9008939 
3.9008005 
3.9007993 
3.9007993 
3.9013267 
1 
5.9565367 
5.9557790 
5.9555485 
5.9555453 
5.9555453 
5.9568515 

2.5 
16.471277 
16.464190 
16.461682 
16.461695 
16.461695 
16.474162 

0.2 
0.5 
7.3288709 
7.3267856 
7.3260518 
7.3260551 
7.3260551 
7.3297205 
1 
10.823112 
10.818284 
10.816484 
10.816503 
10.816503 
10.825062 

2.5 
25.831823 
25.796501 
25.778487 
25.779217 
25.779217 
25.845245 
Fig. 4. Variation of dimensionless fundamental frequency with nonlocal parameter for various boundary conditions
Variation of the first four dimensionless frequencies of simplysupported nanoplate with thickness ratio is displayed in Fig. 5. As the thickness of the plate increases, frequency rises because increasing the thickness would improve the strain energy and rigidity of the structure as shown in this figure. Moreover, this figure reveals that nonlocal effect is insignificant for thin plates. As boundary condition gets stronger support, the rigidity of structure rises and frequency of vibration increases. Thus, frequency parameter is lowest in SSSS and highest in CCCC just as Fig. 6 approves this. The effect of aspect ratio on fundamental frequency parameter for SSSS, CSCS and CCCC nanoplates is plotted in Figs. 79. For a constant length of nanoplate, the width of nanoplate gets smaller by increasing the aspect ratio. On the other hand, it is obvious that the dynamic behavior of nanoplate is considerably dependent on the dimensions of nanoplate. Hence, as these figures depict, the nonlocal effect is more notable for higher aspect ratios. This behavior was observed for other two boundary conditions.
Fig. 5. Variation of first four dimensionless frequencies of SSSS nanoplate with length to thickness ratio
Fig. 6. Variation of dimensionless fundamental frequency with thickness ratio for various boundary conditions
Fig. 7. Variation of fundamental dimensionless frequencies of SSSS nanoplate with aspect ratio
Fig. 8. Variation of fundamental dimensionless frequencies of CSCS nanoplate with aspect ratio
Fig. 9. Variation of fundamental dimensionless frequencies of CCCC nanoplate with aspect ratio
Figs. 1011 show the effect of external voltage on dimensionless fundamental frequency for a piezoelectric nanoplate. It can be seen from the Fig. 9 that fundamental frequency is quite dependent on the external voltage. Exerting negative and positive voltage creates the compressive and tensile forces, respectively. Compressive force decreases or weakens the stiffness and tensile force, increases or improves the stiffness so fundamental frequencies decrease by increasing the external voltage. Additionally, it can be understood from Fig. 10 that decreasing the frequency is more prominent in larger values of the nonlocal parameters.
It is observed from Figs. 1112 that changing temperature would not have much effect on fundamental frequency for all sets of boundary condition and increasing the temperature, causes a slight reduction in the stiffness and natural frequency of piezoelectric nanoplate. It is seen from Figs. 1013 that external voltage has the most effect on the fundamental frequency; while changing the temperature has not much influence on fundamental frequency. This is due to the fact that the coefficient of electrical load is much more than the coefficient of thermal load according to the Eqs. (18) and(19), i.e. .
Fig. 10. Variation of fundamental frequency of SSSS nanoplate with external voltage(L_{1}=50nm, L_{2}=25nm h=10nm)
Fig. 11. Variation of dimensionless fundamental frequency of nanoplate with external voltage for various boundary conditions
Fig. 12. Variation of dimensionless fundamental frequency of SSSS nanoplate with increasing temperature (S=0.2)
Fig. 13. Variation of dimensionless fundamental frequency of nanoplate with increasing temperature (S=0.2, μ=0)
Fig. 14 shows the effects of Winkler and shearing layer coefficient on fundamental frequency. As it could be seen, foundation enhances the frequency of the structure by increasing the strain energy and rigidity of nanoplate. Furthermore, it is observed that fundamental frequency is considerably dependent on shearing layer coefficient than Winkler parameter for all sets boundary conditions. In order to show the effect of foundation parameters on the frequency parameter, the variation of the dimensionless fundamental frequency of the structure versus foundation parameters is plotted in Figs. 1516. In addition, it is seen that nonlocal effects are more prominent in larger values of foundation parameters.
The free vibration analysis of piezoelectric nanoplates subjected to electricalthermal loads resting on the elastic foundation using nonlocal elasticity theory based on the various modified shear deformation theories was studied. Two new distributions of shear stress along thickness were introduced for the first time in this article. Governing equations were derived using Hamilton’s principle. The system of governing equations was converted to a single partial differential equation using a simple approach which is the equation of transverse vibration of nanoplate. The transverse vibration equation was solved for five different boundary conditions including SSSS, CSSS, CSCS, CCSS, CCCC, and effects of different parameters such as thickness to length ratio, aspect ratio, increasing temperature, external voltage, foundation coefficients and length scale on natural frequencies were illustrated in detail. The numerical results show that:
Fig. 14. Variation of dimensionless fundamental frequency of nanoplate with increasing foundation parameter (S=0.2, μ=0.1)
Fig. 15. Variation of dimensionless fundamental frequency of SSSS nanoplate with increasing shearing layer coefficient
Fig. 16. Variation of dimensionless fundamental frequency of SSSS nanoplate with Winkler coefficient
References
[1] Jiang S, Zhou D, Gong S, Lu W. Study of piezoelectric ceramic materials for hightemperature and highfrequency applications. Sensors and Actuators A: Physical. 1998; 69(1): 14.
[2] Zhang S, Xia R, Lebrun L, Anderson D, Shrout TR. Piezoelectric materials for high power, high temperature applications. Materials letters. 2005; 59(27): 34715.
[3] Wang ZL. Piezoelectric nanostructures: From growth phenomena to electric nanogenerators. MRS bulletin. 2007;32(2):10916.
[4] Lam DC, Yang F, Chong A, Wang J, Tong P. Experiments and theory in strain gradient elasticity. Journal of the Mechanics and Physics of Solids. 2003; 51(8): 1477508.
[5] Yang F, Chong A, Lam DCC, Tong P. Couple stress based strain gradient theory for elasticity. International Journal of Solids and Structures. 2002; 39(10): 273143.
[6] Eringen AC. On differential equations of nonlocal elasticity and solutions of screw dislocation and surface waves. Journal of applied physics. 1983; 54(9): 470310.
[7] He X, Kitipornchai S, Liew K. Resonance analysis of multilayered graphene sheets used as nanoscale resonators. Nanotechnology. 2005; 16(10): 2086.
[8] Kitipornchai S, He X, Liew K. Continuum model for the vibration of multilayered graphene sheets. Physical Review B. 2005; 72(7): 075443.
[9] Behfar K, Naghdabadi R. Nanoscale vibrational analysis of a multilayered graphene sheet embedded in an elastic medium. Composites science and technology. 2005; 65(78): 115964.
[10] Pradhan S, Phadikar J. Nonlocal elasticity theory for vibration of nanoplates. Journal of Sound and Vibration. 2009; 325(12): 20623.
[11] Aghababaei R, Reddy J. Nonlocal thirdorder shear deformation plate theory with application to bending and vibration of plates. Journal of Sound and Vibration. 2009; 326(12): 27789.
[12] Khorshidi K, Asgari T, Fallah A. Free vibrations analysis of functionally graded rectangular nanoplates based on nonlocal exponential shear deformation theory. Mechanics of Advanced Composite Structures. 2015;2(2):7993.
[13] Alibeigloo A. Free vibration analysis of nanoplate using threedimensional theory of elasticity. Acta Mechanica. 2011; 222(12): 149.
[14] Ding H, Chen W. Three dimensional problems of piezoelasticity: Nova Science Publishers; 2001.
[15] Tiersten HF. Linear Piezoelectric Plate Vibrations: Elements of the Linear Theory of Piezoelectricity and the Vibrations Piezoelectric Plates: Springer; 2013.
[16] Yang J. An introduction to the theory of piezoelectricity: Springer Science & Business Media; 2004.
[17] Khorshidi K, Pagoli M. Analytical solution for sound radiation of vibrating circular plates coupled with piezoelectric layers. Mechanics of Advanced Composite Structures. 2016; 3(2): 8998.
[18] Ke LL, Wang YS. Thermoelectricmechanical vibration of piezoelectric nanobeams based on the nonlocal theory. Smart Materials and Structures. 2012; 21(2): 025018.
[19] Huang GY, Yu SW. Effect of surface piezoelectricity on the electromechanical behavior of a piezoelectric ring. Physica status solidi (b). 2006; 243(4).
[20] Jandaghian AA, Rahmani O. Sizedependent free vibration analysis of functionally graded piezoelectric plate subjected to thermoelectromechanical loading. Journal of Intelligent Material Systems and Structures. 2017; 28(20): 303953.
[21] Nikoei S, Hassani B. Static and free vibration isogeometric analysis of laminated composite plates integrated with piezoelectric using Reissner–Mindlin theory. Modares Mechanical Engineering. 2018; 17(11): 18191.
[22] Pasternak P. On a new method of analysis of an elastic foundation by means of two constants [Gosudarstvennoe Izdatelstvo Literaturi po Stroitelstvu I Arkhitekture]. Moscow: USSR. 1954.
[23] Vlasov V. Structural mechanics of thinwalled threedimensional systems. Stroizdat, Moscow. 1949.
[24] FilonenkoBorodich M. Some approximate theories of elastic foundation. Uchenyie Zapiski Moskovkogo Gosudarstuennogo Universiteta Mekhanika, Moscow. 1940; 46: 318.
[25] RahbarRanji A, Shahbaztabar A. Free vibration analysis of nonhomogeneous orthotropic plates resting on Pasternak elastic foundation by RayleighRitz method. Journal of Central South University. 2016; 23(2): 41320.
[26] Zhong R, Wang Q, Tang J, Shuai C, Liang Q. Vibration characteristics of functionally graded carbon nanotube reinforced composite rectangular plates on Pasternak foundation with arbitrary boundary conditions and internal line supports. Curved and Layered Structures. 2018; 5(1): 1034.
[27] Moradidastjerdi R, MomeniKhabisi H, Baghbani R. Meshfree Dynamic Analyses of FGM Sandwich Plates Resting on A Pasternak Elastic Foundation. Mechanics of Advanced Composite Structures. 2017; 4(2): 15368.
[28] Moradi Dastjerdi R, Payganeh G, Rajabizadeh Mirakabad S, Jafari MofradTaheri M. Static and free vibration analyses of functionally graded nanocomposite plates reinforced by wavy carbon nanotubes resting on a Pasternak elastic foundation. Mechanics of Advanced Composite Structures. 2016 ;3(2): 12335.
[29] Kirchoff G. Uber das Gleichgewicht und die Bewegung einer elastischen Scheibe. Journal fur die reine und angewandte Mathematik (Crelle's Journal). 1850; 40: 5188.
[30] Kirchhoff G. Ueber die Schwingungen einer kreisförmigen elastischen Scheibe. Annalen der Physik. 1850; 157(10): 25864.
[31] Reissner E. The effect of transverse shear deformation on the bending of elastic plates. J appl Mech. 1945: A69A77.
[32] Mindlin RD. Influence of rotatory inertia and shear on flexural motions of isotropic, elastic plates. J appl Mech. 1951; 18: 318.
[33] Reddy JN. A simple higherorder theory for laminated composite plates. Journal of applied mechanics. 1984; 51(4): 74552.
[34] Touratier M. An efficient standard plate theory. International journal of engineering science. 1991; 29(8): 90116.
[35] Karama M, Afaq K, Mistou S. A new theory for laminated composite plates. Proceedings of the Institution of Mechanical Engineers, Part L: Journal of Materials: Design and Applications. 2009; 223(2): 5362.
[36] Reissner, E. On transverse bending of plates, including the effects of transverse shear deformation. International Journal of Solids and Structures;1975; 11; 569–573.
[37] Panc V. Theories of elastic plates: Springer Science & Business Media; 1975.
[38] Sayyad AS, Ghugal YM. Bending and free vibration analysis of thick isotropic plates by using exponential shear deformation theory. Applied and Computational mechanics. 2012; 6(1): 65–82.
[39] Ghugal YM, Sayyad AS. Free vibration of thick orthotropic plates using trigonometric shear deformation theory. Latin American Journal of Solids and Structures. 2011; 8(3): 22943.
[40] Ke LL, Liu C, Wang YS. Free vibration of nonlocal piezoelectric nanoplates under various boundary conditions. Physica E: Lowdimensional Systems and Nanostructures. 2015; 66: 93106.
[41] Khorshidi K, Siahpush A, Fallah A. ElectroMechanical free vibrations analysis of composite rectangular piezoelectric nanoplate using modified shear deformation theories. Journal of Science and Technology of Composites. 2017;4(2):15160.
[42] Zenkour A, Sobhy M. Nonlocal elasticity theory for thermal buckling of nanoplates lying on Winkler–Pasternak elastic substrate medium. Physica E: Lowdimensional Systems and Nanostructures. 2013; 53: 2519.
[43] Liu D, Li X. An overall view of laminate theories based on displacement hypothesis. Journal of composite materials. 1996;30(14):153961.
[44] Sobhy M. Buckling and free vibration of exponentially graded sandwich plates resting on elastic foundations under various boundary conditions. Composite Structures. 2013; 99:7687.
[45] Liu C, Ke LL, Wang YS, Yang J, Kitipornchai S. Thermoelectromechanical vibration of piezoelectric nanoplates based on the nonlocal theory. Composite Structures. 2013; 106: 16774.
[46] Reddy J, Robbins D. Theories and computational models for composite laminates. Applied mechanics reviews. 1994;47(6):14769.