Document Type: Research Paper
Authors
Department of Solid Mechanics, Faculty of Mechanical Engineering, University of Kashan, Kashan, Iran
Abstract
Keywords
Nonlinear MagnetoNonlocal Vibration Analysis of Coupled Piezoelectric MicroPlates Reinforced with Agglomerated CNTs
S. Amir^{*}, A.R. Vossough, H. Vossough, E. Arshid
Department of Solid Mechanics, Faculty of Mechanical Engineering, University of Kashan, Kashan, Iran
KEYWORDS 

ABSTRACT 
Coupled system Agglomerated CNTs DQM Piezoelectric Magnetic field


The aim of this article is to analyze nonlinear electromagneto vibration of a doublepiezoelectric composite microplatesystem (DPCMPS) pursuant to the nonlocal piezoelasticity theory. The two microplates are assumed to be connected by an enclosing elastic medium, which is simulated by the Pasternak foundation. Both of piezoelectric composite microplates are made of polyvinylidene fluoride (PVDF) reinforced by agglomerated carbon nanotubes (CNTs). The MoriTanaka model is employed to compute the mechanical properties of composite. Applying nonlinear straindisplacement relations and contemplating charge equation for coupling between electrical and mechanical fields, the motion equations are derived in consonance to the energy method and Hamilton's principle. These equations can't be solved analytically as a result of their nonlinear terms. Hence, the differential quadrature method (DQM) is employed to solve the governing differential equations for the case when all four ends are clamped supported and free electrical boundary conditions. The frequency ratio of DPCMPS is inspected for three typical vibrational states, namely, outofphase, inphase and the case when one microplate is fixed in the DPCMPS. A detailed parametric study is conducted to scrutinize the influences of the small scale coefficient, stiffness of the internal elastic medium, the volume fraction of the CNTs, agglomeration and magnetic field. The results reveal that with increasing volume fraction of the CNTs, the frequency of the structure increases. This study might be beneficial for the design and smart control of nano/micro devices such as MEMS and NEMS. 
Nanocomposites hold the promise of advances that exceed those achieved in recent decades in composite materials. The nanostructure generated by a nanophase in the polymer matrix represents a radical alternative to the structure of conventional polymer composites. These complex hybrid materials integrate the predominant surfaces of nanoparticles and the polymeric structure into a novel nanostructure, which produces critical fabrication and interface implementations leading to extraordinary properties [1]. PVDF is an ideal piezoelectric matrix as a result of characteristics including flexibility in thermoplastic conversion techniques, excellent dimensional stability, abrasion and corrosion resistance, high strength, and capability of maintaining its mechanical properties at elevated temperature. Consequently, it has found multiple applications in nanocomposites in a wide range of industries including oil and gas, petrochemical, wire and cable, electronics, automotive, and construction. Boron nitride nanotubes (BNNTs) applied as the matrix reinforcers, apart from having high mechanical, electrical, and chemical properties, present more resistant to oxidation than other conventional nano reinforces such as carbon nanotubes (CNTs). Hence, they are applied for hightemperature applications [26]. Both PVDF and BNNT are smart materials since they have piezoelectric properties.
Piezoelectricity is a classical discipline traced to the original work of Jacques and Pierre Curie around 1880. This phenomenon describes the relations between mechanical strains on a solid and its resulting electrical behavior resulting from changes in the electric polarization. One can generate an electrical output from a solid resulting from mechanical strains or can create a mechanical distortion resulting from the application of an electrical perturbation [7]. Piezoelectric materials have been applied to manufacture various sensors, conductors, actuators, etc. in fact, they have become one of the smart materials nowadays [8].
Regarding research development into the application of smart nanocomposite (which PVDF and BNNT are as matrix and reinforcer, respectively), Barzoki et al. [8] explored electrothermomechanical torsional buckling of a piezoelectric polymeric cylindrical shell reinforced by DWBNNTs with an elastic core. They concluded that the higher the infill core, the higher is dimensionless critical torsional buckling load. In another research, Barzoki et al. [9] investigated nonlinear buckling response of embedded piezoelectric cylindrical shell reinforced with BNNT under electrothermomechanical loadings applying harmonic differential quadrature method (HDQM). They discovered that the critical buckling load increases when the piezoelectric effect is considered. Ghorbanpour et al. [10] demonstrated nonlinear vibration and stability of a smart composite micro tube made of PVDF reinforced by BNNTs embedded in an elastic medium under electro thermal loadings is inspected. They concluded that the stability of the system is strongly dependent on the imposed electric potential and the volume percent of BNNTs reinforcement.
In recent years, small scale effect in micro and nano applications of the beam, plate, and shell type structures has been utilized on the basis of nonlocal elasticity theory which was initiated in the papers of Eringen [1113]. He regarded the stress state at a given point as a function of the strain states of all points in the body, while the local continuum mechanics assumes that the stress state at a given point depends uniquely on the strain state at the same point. Shen et al. [14] explored the nonlocal plate model for nonlinear vibration of singlelayer graphene sheets (SLGS) in thermal environments. Their results revealed that with properly selected small scale parameters and material properties, the nonlocal plate model could provide a remarkably accurate prediction of the graphene sheet behavior under nonlinear vibration in the thermal environment. Pradhan and Kumar [15] reported vibration analysis of orthotropic graphene sheets applying nonlocal elasticity theory and DQM. Their results indicated that the nonlocal effect increases as size of graphene sheet are decreased. Amir [16] examined Orthotropic patterns of viscoPasternak foundation in nonlocal vibration of orthotropic graphene sheet under thermomagnetic fields based on new firstorder shear deformation theory. The results indicate that the stability of singlelayer graphene sheet is strongly dependent on applied magnetic field. Ghorbanpour et al. [17] investigated Pasternak foundation effect on the axial and torsional wave propagation in the embedded doublewalled carbon nanotubes (DWCNTs) applying nonlocal elasticity cylindrical shell theory. They concluded that the frequencies are dependent on small scale coefficient and shear modulus of the elastic medium.
The mentioned studies above on the nanostructures are on the basis of the nonlocal elasticity theory, which is not proper for direct application in the piezoelectric materials. Recently, Eringen’s nonlocal elasticity theory was extended by Zhou et al. [1820] for the piezoelectric materials. In the nonlocal piezoelectric materials, the stress state and the electric displacement at a given point are, respectively, as a function of the strain state and electric potential of all points in the body. Ke et al. [21] employed a nonlocal piezoelasticity model to nonlinear vibration analyze of the piezoelectric nanobeams. They applied DQM to study the effects of nonlocal parameter, temperature change and the external electric voltage on the nonlinear frequency of the piezoelectric nanobeams. Sobhy and Zenkour [22] discussed the magnetic field effect on thermomechanical buckling and vibration of viscoelastic sandwich nanobeams with CNT reinforced face sheets on a viscoelastic substrate. Moreover, Differential quadrature method for vibration analysis of electrorheological sandwich plate with CNT reinforced nanocomposite face sheets subjected to electric field examined by Ghorbanpour Arani et al. [23]. Furthermore, in smallscales, nonlinear dynamic buckling analysis of embedded micro cylindrical shells reinforced with agglomerated CNTs applying strain gradient theory was researched by Tohidi et al. [24].
With respect to developmental works on mechanical behavior analysis of nano and micro plates, it should be noted that none of the researches mentioned above have considered a coupled doubleplate system. Herein, Murmu and Adhikari [25] analyzed vibration of the nonlocal double nanoplate system (NDNPS). Their study highlighted that the smallscale effects considerably influence the transverse vibration of NDNPS. Also, they elucidated that the boost of the stiffness of coupling springs in the NDNPS reduces the small scale effects during the asynchronous modes of vibration. Moreover, buckling behavior of the NDNPS was investigated by Murmu et al. [26], who indicated that the nonlocal effects in the coupled system are higher within creasing values of the nonlocal parameter for the case of synchronous buckling modes than in the asynchronous buckling modes. Furthermore, their analytical results indicated that the increase of the stiffness of the coupling springs in the doubleGSsystem reduces the nonlocal effects during the asynchronous modes of buckling. Exact solution for nonlocal vibration of doubleorthotropic nanoplates embedded in elastic medium was reported by Pouresmaeeli et al. [27], who manifested that the frequency of double orthotropic nanoplates is always smaller than that of double isotropic nanoplates. The three papers [2527] have contemplated the Winkler model for simulation of elastic medium between two nanoplates. In this simplified model, a proportional interaction between pressure and deflection of SLGS is assumed, which is demonstrated in the form of discrete and independent vertical springs. Whereas, Pasternak suggested taking into account not only the normal stresses but also the transverse shear deformation and continuity among the spring elements, and its subsequent applications for developing the model for buckling analysis, which proved to be more accurate than the Winkler model. Recently, analysis of the coupled system of double layered graphene sheets (CSDLGSs) embedded in a viscoPasternak foundation is manifested by Ghorbanpour Arani et al. [28] who indicated that the frequency ratio of the CSDLGSs is more than the SLGS. To the best of our knowledge, none of the works in the literature have contemplated the nonlinear terms in the governing equations for a coupled system. This study aims to deliberated nonlinear terms for vibration analysis of a DPCMPS in which two microplates are connected by an enclosing Pasternak foundation.
None of the aforementioned studies [2528] have considered smart coupled structures, while these structures may be applied in mechanical behavior control of coupled micro and nano structures. Recently, buckling analysis and smart control of SLGS using elastically coupled PVDF nanoplate using the nonlocal piezoelasticity were examined by Ghorbanpour et al. [29], who indicated that the imposed external voltage is an effective controlling parameter for buckling of the SLGS. Moreover, their results revealed that the effect of external voltage becomes more prominent at higher nonlocal parameter and shear modulus. But paper [29] is linear analysis, and just one of two plates is smart.
However, to date, no study was reported in the literature on the vibration of an elastically coupled DPCMPS. Motivated by these considerations, in order to ameliorate the optimum design of smart microstructure, the authors aim to investigate the electromagneto nonlinear nonlocal vibration of an elastically coupled DPCMPS. Herein, the two PVDF microplates reinforced by agglomerated CNTs are coupled by an enclosing Pasternak foundation. deliberated the nonlinear straindisplacement relations and charge equation, the nonlinear governing equations are derived applying energy method and Hamilton's principle. Hence, the DQM is presented to solve the nonlinear governing equations and estimate the frequency. In the present study, the influences of nonlocal parameters, temperature gradient, elastic medium constants, agglomeration and volume fraction of CNTs and magnetic field in polymer have been taken into account.
Pursuant to the theory of nonlocal piezoelasticity, the stress tensor and the electric displacement at a reference point depend not only on the strain components and electric field components in the same position but on all other points of the body as well. The nonlocal constitutive behavior for the piezoelectric material can be computed as follows [21]:
(1) 

(2) 
where and are, respectively, the nonlocal stress tensor and local stress tensor, and are the components of the nonlocal and local electric displacement. is the nonlocal modulus. is the Euclidean distance, and is defined that l is the external characteristic length, e_{0} denotes a constant appropriate to each material, and a is an internal characteristic length of the material. Consequently, e_{0}a is a constant parameter that is acquired with molecular dynamics, experimental results, experimental studies, and molecular structure mechanics. In order to gain a constitutive equation of the nonlocal elasticity can be written as [30]:
(3) 
where the parameter μ=(e_{0}a)^{2} denotes the small scale effect on the response of structures in nano/ micro size and is the Laplacian operator in the above equation. Similarly, Eq. (2) can be written as [20]:
(4) 
Based on the classical plate theory (CPT), which satisfies Kirchhoff assumption, the displacement field is expressed as [31]:
(5) 
where (u, v, w) denote the total displacements of a point along with the (x, y, z) coordinates and (u_{0}, v_{0}, w_{0}) are the displacements of points on the midplane. The vonKármán nonlinear strains associated with the above displacement field can be expressed in the following form [32]:
(6) 
On the basis of the CPT, shear strains ε_{xz}, ε_{yz} are contemplated negligible. Hence, the strain equations in terms of the midplane displacements are derived by substituting the Eq. (5) into the Eq. (6) as follows:
(7) 
The strain components ε_{xx}, ε_{yy,} and γ_{xy} at an arbitrary point of the sheet are related to the middle surface strains and curvatures tensor as follows:
(8) 
where ( , , ) are components of the membrane strains (middle surface strains) tensor and ( , , ) are components of the bending strain (curvature) tensor.
An elastically coupled DPCMPS having the length l, the width b, and the thickness h, assuming that [32], is illustrated in Fig. 1.
The origin of the Cartesian coordinate system is contemplated at one corner of the middle surface of the microplate. The x, y, and z axes are taken in conformity with the length, width, and thickness of the microplates, respectively. The two microplates are created of PVDF and reinforced by CNTs in xdirection so that both microplates are identical.
Fig. 1. Schematic of doublesmart composite microplatesystem. 
The DPCMPS is subjected to uniform temperature change and polarized in xdirection. The two microplates are coupled by an elastic medium, which is simulated by the Pasternak foundation. As is well known this foundation model is characterized by two parameters: the Winkler constant k_{w} and shear constant k_{g}.
In a piezoelectric material, the application of an electric field to it will cause a strain proportional to the mechanical field strength, and vice versa. Pursuant to a piezoelectric microplate under electrothermal loads, constitutive equations can be represented as [33]:
(9) 

(10) 
where e_{ij}, ϵ_{ij} (i, j=1…,6), α_{k} (k=x, y), and ΔT are piezoelectric constants, dielectric constants, thermal expansion coefficients, and temperature gradient, respectively. C_{ij} is a component of stiffness tensor. Electric field tensor E can be written in term of electric potential ϕ as [34]:
(11) 
In this section, the effective modulus of the composite shell reinforced by CNTs is developed. Different methods are available to estimate the overall properties of a composite [1]. As a result of its simplicity and accuracy, even at high volume fractions of the inclusions, the MoriTanaka method [1] is applied in this section. To begin with, the CNTs are assumed to be aligned and straight with the dispersion of uniform in the polymer. The matrix is assumed to be elastic and isotropic, with Young’s modulus E_{m} and the Poisson’s ratio ν_{m}. The constitutive relations for a layer of the composite with the principal axes parallel to the r, θ and z directions are [1]:
(12) 
where σ_{ij}, ε_{ij}, γ_{ij}, k, m, n, l, and p are the stress components, the strain components, and the stiffness coefficients, respectively. In consonance to the MoriTanaka method, the stiffness coefficients are shown by [1]:
(13) 

where C_{m} and C_{r} are the volume fractions of the matrix and the CNTs respectively, and k_{r} ، l_{r} ، n_{r} ، p_{r} , m_{r} are the Hills elastic modulus for the CNTs [1]. The experimental results reveal that most of the CNTs are bent and centralized in one area of the polymer. These regions with concentrated CNTs are assumed in this section to have spherical shapes and are deliberated as ‘‘inclusions’’ with different elastic properties from the surrounding material. The total volume V_{r} of CNTs can be divided into the following two parts [1]:
(14) 
where and are the volumes of CNTs dispersed in the inclusions, concentrated regions, and in the matrix, respectively. Introduce two parameters ξ and ζ describe the agglomeration of CNTs:
(15) 

(16) 
Despite that, the average volume fraction C_{r} of CNTs in the composite is:
(17) 
Assume that all the orientations of the CNTs are completely random. Hence, the effective bulk modulus (K) and effective shear modulus (G) may be written as:
(18) 

(19) 
where:
(20) 

(21) 

(22) 

(23) 
where χ_{r}, β_{r}, δ_{r}, η_{r} may be calculated as:
(24) 

(25) 

(26) 

(27) 
where K_{m} and G_{m} are the bulk and shear moduli of the matrix which can be written as:
(28) 

(29) 
Furthermore, β and α can be obtained from:
(30) 

(31) 

(32) 
Finally, the elastic modulus (E) and Poisson’s ratio (υ) can be computed as:
(33) 

(34) 
The governing differential equations of motion are derived applying Hamilton's principle, which is given as [35]:
(35) 
where δU is the virtual strain energy which is obtained by the following relation [36]:
(36) 
δK is the virtual kinetic energy and is defined as following [37]:
(37) 
Also, δV is the virtual work operated by externally applied forces and is acquired by:
(38) 
in which q_{m} can be written as:
(39) 
where η is the magnetic permeability; is the gradient operator; H_{x} is the magnetic field. Furthermore, q_{e} can be written as:
(40) 
where k_{w} and k_{g} are spring and shear constants of an elastic medium, respectively.
The motion equations can be derived by applying Eq. (35) as follows:
(41) 
where:
(42) 
where (m_{0}, m_{2}) are mass moments of inertia, and ρ_{0} denotes the density of the material. Meanwhile, the force resultants (N_{xx}, N_{yy}, N_{xy}) and the moment resultants (M_{xx}, M_{yy}, M_{xy}) of the plate can be defined as:
(43) 
Charge equation for coupling electrical and mechanical fields is:
(44) 
In this study, transverse vibration is investigated (i.e., u_{0}=v_{0}=0).
As can be observed, the coupled governing equations contain nonlinear terms and should be solved applying a numerical method such as DQM. In this method, the differential equations are changed into a firstorder algebraic equation by employing appropriate weighting coefficients. Weighting coefficients are not correlated to any special problem and only depend on the grid spacing. For the implementation of the DQ approximation, consider a function f(ζ, η) which has the field on a rectangular domain ( and ) with n_{ζ } n_{η} grid points along x and y axes. According to DQ method, the r^{th} derivative of a function f(x, y) can be defined as [35]:
(45) 
where are weighting coefficients and defined as:
(46) 
where M(ζ_{i}) is Lagrangian operators which can be presented as:
(47) 
The weighting coefficients for the second, third, and fourth derivatives are defined as:
(48) 
In a similar method, the weighting coefficients for ydirection can be acquired. The coordinates of grid points are chosen as:
(49) 
In order to carry out the eigenvalue analysis, the domain and boundary points are separated, and in vector forms, they are denoted as {d} and {b}, respectively. Hence, the discretized form of the motion equations together with the boundary conditions can be expressed in matrix form as:
(50) 
in which [M], [K_{L}], and [K_{NL}] are the mass matrix, linear stiffness matrix, and nonlinear stiffness matrix. This nonlinear equation can now be solved by applying a direct iterative process as follows:
(51) 
where ε_{0} is a small value number, and in the present analysis, it is taken to be 0.1%.
Mechanical, thermal, and electrical properties of PVDF matrix and CNT reinforcement are chosen from Ref. [9]. The final converged solution applying the numerical procedure above is illustrated as the influences of the elastic medium, nonlocal parameter, volume percent of CNT, CNTs agglomeration and temperature change on the frequency of the structure.
Since no reference to such a work is found todate in the literature, its validation is not possible. Despite that, the present work could be partially validated based on a simplified analysis suggested by Shen et al. [14] on thermal nonlinear vibration of the SLGS for which the coupled plate and volume percent of CNTs in the polymer were ignored. For this purpose, a SLGS with C_{r}=0, T=300 K, l=9.496 nm, b=4.877 nm, h=0.145 nm, ρ_{0}=5624 kg/m^{3}, k_{w}=k_{g}=0 and e_{0}a=0.67 nm is considered. Table 1 reveals the result of validation exercise by showing nonlineartolinear frequency for dissimilar dimensionless amplitude (w/h) and temperature. As it can be observed, the results acquired are in good agreement with those expressed in [14].
Fig. 2 manifests the effects of the nonlocal parameter on the frequency versus the orientation angle of CNTs.
As can be observed, the frequency of the system decreases with considering size effects. This is due to the fact that the contemplating of nonlocal parameter decreases the interaction force between microplate atoms, and that leads to a softer structure.
The effect of volume percent of CNTs on the frequency versus the orientation angle of CNTs is shown in Fig. 3. It is clear that the frequency increases with increasing the volume percent of CNTs, and the influence of volume percent of CNTs on the frequency become more prominent at the middle angle. It is because with increasing the volume percent of CNTs, the stiffness of structure increases.
Fig. 4 manifests agglomeration effects on the frequency versus orientation angle of CNTs. As can be seen, contemplating agglomeration effects leads to lower frequency since the stability of the system decreases.
Fig. 5 illustrates the influence of the thermal gradient (ΔT) on the frequency versus the orientation angle of CNTs. It is evident that an increase in temperature change does not considerable effect on the frequency.
Fig. 6 indicates the effect of temperature variations on the results. It is observed by increasing the temperature difference, the frequency enhances as a result to change in the mechanical properties of the structure.
Fig. 2. The effects of the nonlocal parameter on the frequency versus the orientation angle of CNTs. 

Fig. 3. The effects of volume percent of CNTs on the frequency versus the orientation angle of CNTs. 



Fig. 4. The effects of agglomeration of CNTs on the frequency versus the orientation angle of CNTs. 
Table 1. Comparing dimensionless nonlinear frequency obtained in the present study and those of Shen et al. [14].
Ω_{NL}/Ω_{L} 
Ref. 
Temperature 

w/h=2 
w/h=1.5 
w/h=1 
w/h=0.5 

1.2933 
1.1742 
1.0802 
1.0208 
Present work 
T=300 K 
1.2900 
1.1720 
1.0798 
1.0205 
Shen et al. [14] 

0.256 
0.188 
0.037 
0.029 
Difference (%) 

1.4492 
1.2738 
1.1292 
1.0339 
Present work 
T=400 K 
1.4485 
1.2719 
1.1289 
1.0337 
Shen et al. [14] 

0.048 
0.149 
0.027 
0.019 
Difference (%) 

1.8498 
1.5372 
1.2668 
1.0731 
Present work 
T=500 K 
1.8477 
1.5355 
1.2663 
1.0728 
Shen et al. [14] 

0.114 
0.111 
0.039 
0.028 
Difference (%) 
Fig. 5. The effects of thermal gradient on the frequency versus the orientation angle of CNTs. 
Fig. 6. The effects of temperature variations on the frequency. 
The effect of elastic medium on the frequency of the structure is portrayed in Fig. 7. It can be found that considering elastic medium leads to higher frequency. Moreover, deliberating Pasternak medium predicts a higher frequency with respect to Winkler medium. It is due to the Pasternak medium, the normal and shear constant is taken into account.
The effect of foundation constants is discussed in Fig. 8. It can be observed that adding the spring to the structure, leads the frequency to increase and adding the shear layer, has a similar effect on the results as well.
Fig. 7. The effects of the elastic medium on the frequency versus the orientation angle of CNTs. 
Fig. 8. The effects of foundation constants on the frequency. 
The effect of the magnetic field on the frequency of structure is illustrated in Fig. 9. It can be found that increasing the magnetic field, the frequency increases. It is due to the increasing magnetic field leads to higher stiffness.
The vibration response of piezoelectric nano/micro composites has applications in designing many NEMS/MEMS devices such as hydraulic sensors and actuators. In the present study, electromagneto nonlinear vibration of a doublepiezoelectric composite microplate made of PVDF reinforced by CNTs is inspected deliberating agglomeration effects. The internal elastic medium between two microplates is simulated as Pasternak foundation. Considering charge equation, the nonlinear motion equations are derived based on nonlocal piezoelasticity theory. The DQM is applied in order to acquire the nonlinear frequency ratio of the DPCMPS so that the effects of the small scale coefficient, stiffness of the internal elastic medium, the volume fraction, and orientation angle of the CNTs reinforcement, temperature change and agglomeration are discussed.
Fig. 9. The effects of the magnetic field on the frequency versus the orientation angle of CNTs. 
The results of this study are validated by Shen et al. [14]. The results reveal that along with increasing geometrical aspect ratio, the effect of coupling elastic medium between two piezoelectric composite microplates decreases. Furthermore, the effects of small scale parameters and volume percent lead to a higher frequency. It is noteworthy to mention that the frequency of structure considering the agglomeration of CNTs becomes lower.
The authors would like to thank the reviewers for their valuable comments and suggestions to improve the clarity of this study.
The authors are thankful to the University of Kashan for supporting this work by Grant No. 891255/2.
[1] Pissis P. Thermoset nanocomposites for engineering applications: iSmithers Rapra Publishing; 2007.
[2] Merhari L. Hybrid nanocomposites for nanotechnology: Springer; 2009.
[3] Schwartz M. Encyclopedia of Smart Materials, 2 Volume Set. Encyclopedia of Smart Materials, 2 Volume Set, by Mel Schwartz, pp 1176 ISBN 0471177806 WileyVCH, March 2002 2002: 1176.
[4] Topolov VY, Bowen CR. Electromechanical properties in composites based on ferroelectrics: Springer Science & Business Media; 2008.
[5] Yang J. The mechanics of piezoelectric structures: World Scientific; 2006.
[6] Brockmann TH. Theory of Adaptive Fiber Composites: from piezoelectric material behavior to dynamics of rotating structures: Springer Science & Business Media; 2009.
[7] Arnau A. Piezoelectric transducers and applications: Springer; 2004.
[8] Barzoki AM, Arani AG, Kolahchi R, Mozdianfard M. Electrothermomechanical torsional buckling of a piezoelectric polymeric cylindrical shell reinforced by DWBNNTs with an elastic core. Applied Mathematical Modelling 2012; 36(7): 298395.
[9] Barzoki AM, Arani AG, Kolahchi R, Mozdianfard M, Loghman A. Nonlinear buckling response of embedded piezoelectric cylindrical shell reinforced with BNNT under electro–thermomechanical loadings using HDQM. Composites Part B: Engineering 2013; 44(1): 7227.
[10] Arani AG, Shajari A, Amir S, Loghman A. Electrothermomechanical nonlinear nonlocal vibration and instability of embedded microtube reinforced by BNNT, conveying fluid. Physica E: Lowdimensional Systems and Nanostructures 2012; 45: 10921.
[11] Eringen AC. Nonlocal polar elastic continua. International journal of engineering science 1972; 10(1): 116.
[12] Eringen AC. On differential equations of nonlocal elasticity and solutions of screw dislocation and surface waves. Journal of applied physics 1983; 54(9): 470310.
[13] Eringen AC. Nonlocal continuum field theories: Springer Science & Business Media; 2002.
[14] Shen L, Shen HS, Zhang CL. Nonlocal plate model for nonlinear vibration of single layer graphene sheets in thermal environments. Computational Materials Science 2010; 48(3): 6805.
[15] Pradhan S, Kumar A. Vibration analysis of orthotropic graphene sheets embedded in Pasternak elastic medium using nonlocal elasticity theory and differential quadrature method. Computational Materials Science 2010; 50(1): 23945.
[16] Amir S. Orthotropic patterns of viscoPasternak foundation in nonlocal vibration of orthotropic graphene sheet under thermomagnetic fields based on new firstorder shear deformation theory. Proceedings of the Institution of Mechanical Engineers, Part L: Journal of Materials: Design and Applications 2019; 233(2): 197208.
[17] Arani AG, Barzoki AM, Kolahchi R, Loghman A. Pasternak foundation effect on the axial and torsional waves propagation in embedded DWCNTs using nonlocal elasticity cylindrical shell theory. Journal of mechanical science and technology 2011; 25(9): 2385.
[18] Zhou ZG, Wang B. The scattering of harmonic elastic antiplane shear waves by a Griffith crack in a piezoelectric material plane by using the nonlocal theory. International Journal of Engineering Science 2002; 40(3): 30317.
[19] Zhou ZG, Wu LZ, Du SY. Nonlocal theory solution for a Mode I crack in piezoelectric materials. European Journal of MechanicsA/Solids 2006; 25(5): 793807.
[20] Zhou ZG, Du SY, Wu LZ. Investigation of antiplane shear behavior of a Griffith permeable crack in functionally graded piezoelectric materials by use of the nonlocal theory. Composite structures 2007; 78(4): 57583.
[21] Ke LL, Wang YS, Wang ZD. Nonlinear vibration of the piezoelectric nanobeams based on the nonlocal theory. Composite Structures 2012; 94(6): 203847.
[22] Sobhy M, Zenkour AM. Magnetic field effect on thermomechanical buckling and vibration of viscoelastic sandwich nanobeams with CNT reinforced face sheets on a viscoelastic substrate. Composites Part B: Engineering 2018; 154: 492506.
[23] Arani AG, Jamali S, Zarei HB. Differential quadrature method for vibration analysis of electrorheological sandwich plate with CNT reinforced nanocomposite facesheets subjected to electric field. Composite Structures 2017; 180: 21120.
[24] Tohidi H, HosseiniHashemi S, Maghsoudpour A. Nonlinear sizedependent dynamic buckling analysis of embedded micro cylindrical shells reinforced with agglomerated CNTs using strain gradient theory. Microsystem Technologies 2017; 23(12): 572744.
[25] Murmu T, Adhikari S. Nonlocal vibration of bonded doublenanoplatesystems. Composites Part B: Engineering 2011; 42(7): 190111.
[26] Murmu T, Sienz J, Adhikari S, Arnold C. Nonlocal buckling behavior of bonded doublenanoplatesystems. Journal of Applied Physics 2011; 110(8): 084316.
[27] Pouresmaeeli S, Fazelzadeh S, Ghavanloo E. Exact solution for nonlocal vibration of doubleorthotropic nanoplates embedded in elastic medium. Composites Part B: Engineering 2012; 43(8): 338490.
[28] Arani AG, Shiravand A, Rahi M, Kolahchi R. Nonlocal vibration of coupled DLGS systems embedded on ViscoPasternak foundation. Physica B: Condensed Matter 2012; 407(21): 412331.
[29] Arani AG, Kolahchi R, Vossough H. Buckling analysis and smart control of SLGS using elastically coupled PVDF nanoplate based on the nonlocal Mindlin plate theory. Physica B: Condensed Matter 2012; 407(22): 445865.
[30] Amir S, Bidgoli EMR, Arshid E. Sizedependent vibration analysis of a threelayered porous rectangular nano plate with piezoelectromagnetic face sheets subjected to pre loads based on SSDT. Mechanics of Advanced Materials and Structures 2018: 115.
[31] Reddy JN. Mechanics of laminated composite plates and shells: theory and analysis: CRC press; 2004.
[32] Vinson JR. Plate and panel structures of isotropic, composite and piezoelectric materials, including sandwich construction: Springer Science & Business Media; 2006.
[33] Tan P, Tong L. Microelectromechanics models for piezoelectricfiberreinforced composite materials. Composites science and technology 2001; 61(5): 75969.
[34] Arshid E, Kiani A, Amir S. Magnetoelectroelastic vibration of moderately thick FG annular plates subjected to multi physical loads in thermal environment using GDQ method by considering neutral surface. Proceedings of the Institution of Mechanical Engineers, Part L: Journal of Materials: Design and Applications 2019: 1464420719832626.
[35] Arshid E, Khorshidvand AR. Free vibration analysis of saturated porous FG circular plates integrated with piezoelectric actuators via differential quadrature method. ThinWalled Structures 2018; 125: 22033.
[36] Arshid E, Khorshidvand AR, Khorsandijou SM. The effect of porosity on free vibration of SPFG circular plates resting on viscoPasternak elastic foundation based on CPT, FSDT and TSDT. Structural Engineering and Mechanics 2019; 70(1): 97112.
[37] Arshid E, Khorshidvand AR. Flexural Vibrations analysis of saturated porous circular plates using differential quadrature method. Iranian Journal of Mechanical Engineers 2017; 19: 78100.
[38] Arani AG, Kolahchi R, Barzoki AAM, Mozdianfard MR, Farahani SMN. Elastic foundation effect on nonlinear thermovibration of embedded doublelayered orthotropic graphene sheets using differential quadrature method. Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science 2013; 227(4): 86279.