Nanoscale mass nanosensor based on the vibration analysis of embedded magneto-electro-elastic nanoplate made of FGMs via nonlocal Mindlin plate theory

In the current paper, the sensitivity performance of functionally graded magneto-electro-elastic (FG-MEE) nanoplate with attached nanoparticles as a nanosensor is analyzed based on nonlocal Mindlin plate assumption. Power law distribution model is employed to display how the material properties of FG-MEE nanoplate vary across the thickness direction. It is supposed that FG-MEE nanoplate is under initial external electric and magnetic potentials. Boundary condition of each edge of FG-MEE nanoplate is assumed to be simply supported. Furthermore, a Pasternak substrate is applied for modelling the total reaction pressure between nanoplate and foundation. Partial differential equations and corresponding boundary conditions are first achieved using Hamilton’s variational principle and then analytically solved to determine the frequency shift utilizing Navier’s approach. Numerical examples are performed to elucidate the dependency of the sensitivity performance of FG-MEE nanosensor on the volume fraction exponent, nonlocal parameter, total attached mass and location of the nanoparticle, aspect ratio, mode number, initial external electric voltage, initial external magnetic potential, and Pasternak medium coefficients. It is clearly indicated that these factors have highly significant impacts on the variations of frequency shift.


Introduction
In the recent years, with the rapid growth in the use of the composite materials in a wide range of engineering structures such as aircraft structures (Baker and Baker 2004) and sensor and actuators (Tressler et al. 1999;Akdogan et al. 2005;Li et al. 2008) because of their excellent mechanical properties, the importance of recognizing the composite material behavior increases day to day. Magneto-electro-elastic (MEE) composite materials are an important class of smart composite materials combining piezoelectric and piezomagnetic phases. The capacity of converting energy among magnetism, electricity, or elasticity into another form has caused that the MEE materials become suitable for smart applications. Up to now, various analytical or numerical works have been carried out to investigate MEE plate, beam, and shell which contain studies on the static behavior (Wu and Tsai 2007;Wu et al. 2010), linear and nonlinear vibration (Ramirez et al. 2006;Razavi and Shooshtari 2015;Shirbani et al. 2017) and on the post-buckling behavior (Zhou et al. 2003). For the first time, the excellent behavior of MEE composite materials experimentally was reported by Van Den Boomgaard et al. (1974). A discrete layer model is applied to predict the natural frequencies of orthotropic functionally graded (FG) MEE plates combining various percentages of barium titanate (BaTiO 3 ) and cobalt iron oxide (CoFe 2 O 4 ) by Ramirez et al. (2006). Nonlinear forced vibration of a MEE plate is investigated via first-order shear deformation theory by Shooshtari and Razavi (2016). Liu et al. (2016) developed a MEE 3D model for studying the components of the MEE field. Surface effects include surface piezoelectricity, surface elasticity and surface piezomagneticity investigated on the bending, buckling and free vibration of the MEE Euler-Bernoulli beams by Xu et al. (2016).
During the recent years, the investigation of the mechanical behavior of MEE structures at nano-and micro-scales have attracted the interest of some researchers. It have proven that the classical approaches cannot predict satisfactorily the mechanical behavior of nano-and micro-structures (Pradhan and Phadikar 2009; Ansari et al. 2010;Hosseini et al. 2016a, b, c). Experiments and atomistic simulations (Chen et al. 2006;Stan et al. 2007) demonstrated that the influence of small size becomes significant when the dimensions of the structure become small and reach to the micro/nano length scale. Since the traditional classical elasticity and plasticity theories have no material internal length scale, these theories are not capable of taking into account the influence of small size scale. Therefore, these theories modified and included one or several independent internal length scale parameters. Some of the most popular modified theories include nonlocal elasticity theory suggested by Eringen (1972Eringen ( , 1983, modified couple stress theory presented (MCST) by Yang et al. (2002) and modified strain gradient theory (MSGT) suggested by Lam et al. (2003).
Our literature review here will focus on the analysis of MEE nanostructures based on the nonlocal elasticity of Eringen. According to the suggested model of Eringen, the elastic strain is described by a Fredholm type integral equations in which the stress is result of a convolution between the local response to an elastic strain and a smoothing kernel dependent on a nonlocal parameter. But some researchers showed the constitutive boundary conditions and paradoxes of nonlocal integral model (Romano and Romano et al. 2017a, b;Apuzzo et al. 2017). They demonstrated that Eringen's nonlocal integral model admits no solution. On the basis of the nonlocal elasticity theory of Eringen, Ansari et al. (2015c) proposed a third-order shear deformable beam model to study forced vibration behavior of magnetoelectro-thermo-elastic nanobeams by taking into account the effects of small scale. Recently, by considering the von Kármán's assumption, the influences of the external electric and magnetic potential on the nonlinear natural frequencies of the MEE nanoplates are shown by Farajpour et al. (2016). Wang et al. (2016) presented a novel twodimensional linear elastic theory for demonstrating surface, nonlocal effects, poling directions, piezoelectric phase materials on the MEE nanoplate. Numerical results illustrated that driving frequency and output power density of magnetic energy harvesters have obvious size-dependent behavior because of surface effects. On the basis of the nonlocal Mindlin plate approach, Li et al. (2014) analytically carried out elastic instability and free vibration behavior of MEE nanoplate resting on Pasternak medium. It was shown that buckling load and natural frequency decrease (increase) linearly with electric potential (magnetic potential). Recently, the generalized differential quadrature method in conjunction with pseudo arc-length continuation are applied to calculate the nonlinear buckling and postbuckling of MEE Mindlin nanoplate under thermal load by Ansari and Gholami (2017). A mathematical model based on Euler-Bernoulli beam assumption is proposed by   Ansari et al. 2015b;Li et al. 2016;Ma et al. 2017) and MEE nanoplates Wu et al. 2015;Kiani et al. 2017) by taking into account the nonlocal parameter. Functionally graded materials (FGMs) are a special subgroup of novel composite materials which contain the heterogeneous property and have a continuous variation of material properties from one surface to another. Due to the numerous profits containing high mechanical and thermal resistance, FGMs are vastly used in sensors and thermogenerators (Müller et al. 2003;Rahimi et al. 2011), biomedical fields (Pompe et al. 2003;Parthasarathy et al. 2011), aerospace (Kumar et al. 2013). With respect to the development of science and wide application of FGM nanostructures, the study of the size effect on the mechanical behavior of FGM nanostructures is of considerable practical interest. On the basis of nonlocal elasticity theory in conjunction with Gurtin-Murdoch elasticity theory, Hosseini and Jamalpoor (2015) reported the important role of the surface effects, i.e., surface elasticity, residual stresses and surface density on the free vibration of a double-FGM viscoelastic nanoplates-system under thermal load. They considered that the material traits of the plate follow power law repartition in the thickness direction.
Numerical results expressed that an increase in the amount of the power law index makes a reduction of the natural frequencies. A modified power-law model in conjunction with the refined four-variable plate assumption are implemented by  to study the influence of FGMs on the vibration frequency of smart piezoelectric plates under various boundary conditions. Wave propagation phenomena of a FG-MEE nanorod is investigated via nonlocal continuum mechanics by Narendar (2016). Furthermore, there are a number of literatures regarding the effects of FGMs on the nanostructures (Natarajan et al. 2012;Ansari et al. 2015a;Barretta et al. 2016;Thang et al. 2017;Mechab et al. 2016;Jamalpoor and Kiani 2017;Hosseini et al. 2017).
It should be pointed that sometimes one or more concentrated particles are attached to the structures. Thus, we need the mass detection technology for detecting the position of attached masses with high accuracy and sensitivity. In order to find the exact point of attached masses, sensors play a significant role. The main idea for determining the location of added masses such as bacterium/ virus, biomolecules, and buckyballs is to measure the resonant frequency shift of the nanosensor caused by variations in total mass of the system (Zhou et al. 2014;Asemi et al. 2015;Karličić et al. 2015;Shi et al. 2015;Jalali et al. 2015;Sadeghzadeh 2016). It is that piezoelectric (ZnO-BaTiO 3 ) (Guo et al. 2012;Alluri et al. 2015), magnetic/ piezoelectric (Yeh et al. 2016), and magnetic (CoFe 2 O 4 )  (Kim et al. 2009;Taei et al. 2016) can be suitable candidates for sensor devices because of having excellent properties like dimensional stability, abrasion, and corrosion resistance.
The above-reviewed papers explicitly illustrated that while there exist a considerable number of studies that deal with the static and dynamic responses of MEE nanoplate based on different (Kirchhoff or Mindlin) plate assumptions, nanoscale mass detection based on vibration FG-MEE Mindlin nanoplate seem to be nonexistent. Accordingly, in this work, a first-order shear deformable (Mindlin) plate theory in conjunction with the nonlocal elasticity theory is applied to take into account the size effect on the sensitivity of the sensor. In this regard, the key novelties of the presented study are summarized as follows: • Development of a nonlocal Mindlin plate theory for free vibration of FG-MEE nanoplate resting on the Pasternak medium. • FG-MEE Mindlin nanoplate with attached concentrated nanoparticles is presented as a mass nanosensor according to the vibration analysis.
Partial differential equations of motion of the system with multiple added masses are derived by employing Hamilton's principle, then by applying the analytical Navier type solution, the frequency shift of the system is presented in the explicit closed-form.

Geometrical configuration
Consider a simply supported thick rectangular (l x Â l y ) nanoplate with uniform thickness h consisting of BaTiO 3 -CoFe 2 O 4 MEE, which is exposed to electric U x; y; z; t ð Þand magnetic potentials W x; y; z; t ð Þ, as shown in Fig. 1. As can be seen from Fig. 1, it is considered that nanoplate with attached nanoparticles is resting on the Pasternak-type foundation with the transverse and shear coefficients k w and k G , respectively. Also, the location of any rigid nanoparticle is assumed to be at (x e ; y e Þ; e ¼ 1; 2; . . .. Furthermore, the material properties of MEE nanoplate are supposed to be graded in the thickness direction based on the power-law model in terms of the volume fractions of the constituents. It should be noted that in FG-MEE nanoplate, the materials properties vary continuously from a fully BaTiO 3 phase at the top surface (z ¼ h 2 ) to a fully CoFe 2 O 4 phase at the bottom surface (z ¼ À h 2 ).

FGMs
On the basis of the power law distribution, the effective material characteristics of MEE nanoplate vary continuously across the plate thickness and can be proposed as below: where, PM 1 and PM 2 are the properties of CoFe 2 O 4 phase and BaTiO 3 phase, respectively. Also, V 1 and V 2 refer to the volume fractions of CoFe 2 O 4 rich and BaTiO 3 rich, respectively, which can be written as in which K is the volume fraction exponent and always has a positive value.

Displacement field
In order to obtain more accurate estimation of mechanical behavior of thick plate, we apply Mindlin plate hypothesis (first-order shear deformation theory) which incorporates in which u 0 and v 0 denote the in plane deflections of the mid-surface; h x and h y refer to the rotation of the middle plane in the x and y directions, respectively; and w 0 is the lateral displacement.

Strains
The nonzero components of the linear strain-displacement relationships can be expressed as

Nonlocal elasticity theory
It is considered via Eringen's nonlocal elastic theory (Eringen 1972(Eringen , 1983) that the stress, electric displacement, and magnetic induction at a reference spot of MEE nanostructures are supposed to be a function of the strain, electric and magnetic components at each spot in material domain of body. According to the main idea of this nonlocal theory, the basic relations for FG-MEE substance without considering body force can be written as follows (Ansari et al. 2015b;Farajpour et al. 2016) where the terms r nl Ã T (the superscript T signifies transpose) represent, respectively, the nonlocal forms of the stress tensor, electric displacement and magnetic induction. Also, r l ij , D l i and B l i demonstrate, respectively, the local forms of the stress tensor, electric displacement and magnetic induction.
indicates the Kernel function, and the Euclidean distance between reference point x and any points x 0 in the body is illustrated with the symbol x À x 0 j j. Furthermore, in Eq. (5), E m and H n are electric field and magnetic field vectors, respectively. Also, symbol V is the volume of MEE q 33 0 699.7 g 11 (10 -9 C 2 /Nm 2 ) 11.2 0.08 g 22 11.2 0.08 g 33 12.6 0.093 substance. Nonlocal parameter that takes into account the small scale effects into the constitutive equations is demonstrated by e 0 a b , which e 0 is a constant parameter determined experimentally for matching theoretical model with empirical results; a is the internal characteristic length (e.g. lattice parameter, granular size); and b refers to an external property of length (e.g. crack length, wavelength). However, the strain-driven integral constitutive nonlocal relation (5) proposed by Eringen (1983) is inapplicable to nonlocal elastic problems on bounded structural domains Barretta 2016, 2017a, b;Romano et al. 2017a, b;Apuzzo et al. 2017). In this paper we assume the following phenomenological differential constitutive law (consequent but not equivalent to Eq. 5) to investigate size-dependent behavior of Mindlin plates.
in which r 2 ¼ o 2 ox 2 þ o 2 oy 2 denotes the two-dimensional Laplacian operator. By not considering the size effect (e 0 a ¼ 0), the traditional classical constitutive relations of FG-MEE structures can be acquired from nonlocal basic relations. ½ C, e, ½ q, g, ½ d and ½ l denote, respectively, the elastic, piezoelectric, piezomagnetic, dielectric, magnetoelectric and magnetic constant matrices, which here are supposed to be varied continuously as a power law function of the plate thickness and are illustrated as below By supposing the plane stress state, the reduced material constants C ij z ð Þ; e ij z ð Þ; q ij z ð Þ; g ij z ð Þ; d ij z ð Þ and l ij z ð Þ can be denoted as follows It should be noted that if the electric and magnetic fields are supposed to be the negative gradients of scalar electric U x; y; z; t ð Þ and magnetic potentials Wðx; y; z; tÞ, respectively, Maxwell equations in the quasi-static approximation are satisfied as However, in this work, the boundary conditions of magnetic and electric fields at the bottom and top surfaces of the FG-MEE nanoplate are taken to be as where the terms / 0 and W 0 (À/ 0 and ÀW 0 Þ indicate, respectively, the external electric potential and magnetic potential exerting to the top (bottom) surface of FG-MEE nanoplate. With respect to the electric-magnetic boundary conditions, the explicit distribution of electric and magnetic potentials are expressed as a combination of cosine and linear variation in the following form (Farajpour et al. 2016;Ansari and Gholami 2017) in which / x; y; t ð Þ and w x; y; t ð Þ demonstrate the spatial variations related to the electric potential and magnetic potential along the x and y directions, respectively.

Governing equations and boundary conditions
In order to determine the governing differential equations of the motion of present model and related boundary conditions, in this section, Hamilton's principle can be applied in the below form as where dT; dU and dV, respectively, denote the virtual kinetic energy, the virtual strain energy and the virtual work exerted by external force q. Here, q refers to the transverse load per unit area of FG-MEE nanoplate due to the added nanoparticle with mass m e for eth mass in position (x e ; y e ), Pasternak elastic medium, and external loads which can be expressed as (Shen et al. 2012;Hosseini et al. 2016a) illustrate the normal loads induced by initial external electric potential and initial external magnetic potential along x and y coordinate directions, respectively. Also, the Dirac delta function in the two two-dimensional state is shown with symbol d x À x e ð Þd y À y e ð Þ, which define as d x À x e ð Þd y À y e ð Þ¼0 forðx; yÞ 6 ¼ ðx e ; y e Þ Z 1 À1 Z 1 À1 d x À x e ð Þd y À y e ð Þdxdy ¼ 1:

Virtual kinetic energy
The virtual kinetic energy on the basis of the first shear order deformation theory can be written as follows where I 0 ; I 1 and I 2 denote the mass inertia terms which can be expressed as I 0 ; I 1 ; I 2 f g¼ R h 2 À h 2 fqðzÞ; qðzÞz; qðzÞz 2 gdz.

Virtual strain energy
The variation of strain energy using the first shear order deformation theory can be expressed as in which parameters N ij ; M ij and Q ii can be defined by the following integral expressions The term k s is the shear correction factor and here equals to 5=6.

Virtual work
Variation of work applied by external loads denotes by

The equilibrium equations in terms of the displacements
Finally, by implementing relations (15-18) into Hamilton's variational principle and integrating by parts, and factoring the coefficients of du 0 ; dv 0 ; dw 0 ; dh x ; dh y ; d/ and dw, the following equilibrium equations of motion are derived as and the related natural and necessary boundary conditions with normal vectors n x and n y are obtained along borders x and y as below The following nonlocal expressions in terms of displacements can be written by integrating Eq. (6) over the area of nanoplate cross-section where parameters N E xxi ; N M xxi ; N E yyi ; N M yyi ; A ij ; B ij ; D ij ; E ij ; F ij ; G ij ; H ij ; J i ; L i ; Q i ; X i and P i are represented in integral forms, which are provided in appendix A. By substituting Eqs. (21-1-21-14) into (19-1-19-7), the neoclassical governing equations of FG-METE Mindlin nanoplate in terms of displacements can be obtained as:

Analytical solution for vibration behavior of FG-MEE nanoplates
In the following section, the simply supported boundary conditions are considered on each the four ends of rectangular nanoplate. Thus, the deflections and moments at each end of FG-MEE nanoplate are assumed to be zero. Also, the short circuit type of edge supports which is used for the actuator case is selected for the electric and magnetic potentials (Ansari and Gholami 2016). On the basis of the short circuit type, the electric potential and magnetic potential are supposed to be zero at the ends of FG-MEE nanoplate. Overall, the related boundary conditions can be expressed as  Ansari and Gholami 2016).
Another comparison study is also carried out through Fig. 2 in which the fundamental frequency shift calculated according to the present model for different values of nonlocal parameter and different added particles at location x 1 ¼ 0:5l x ; y 1 ¼ 0:5l y without considering FGMs, Pasternak medium, and piezomagnetic effect, are compared to those of piezoelectric Kirchhoff nanoplate model of Asemi et al. (2015). Figure 2 denotes an excellent agreement between the outcomes of current paper and the results of Ref. (Asemi et al. 2015).

Benchmark results
In this subsection, the benchmark results are presented in the tabular and graphical form to better understanding and explaining the significant parameters such as the position of single or multiple attached masses on the natural frequency shift changes. To achieve this goal, we consider the three different configurations (single mass, two masses, and four masses) for the number of added particles. It is assumed that the total mass of added masses is identical for each state. However, the details of these three cases are given in Table 2. Furthermore, it is supposed that FG-MEE nanoplate to be composed of barium titanate (BaTiO 3 ) and cobalt iron oxide (CoFe 2 O 4 ) with effective material characteristics given in Table 3 (Aboudi 2001). Moreover, the utilized geometrical characteristics of nanoplate are assumed to be as h ¼ 1 nm; l x ¼ l y ¼ 10 nm, and half wave numbers are assumed to be m ¼ n ¼ 1.
Initially, to study the influences of length to thickness ratio and material property gradient index on the changes of frequency shift, Table 4 is provided for the three figurations of the attached nanoparticles. Here, we took It is clear that an increase in the volume fraction of barium titanate part (BaTiO 3 )/cobalt iron oxide component (CoFe 2 O 4 ) leads to decreasing/increasing the bending stiffness of FG-MEE nanoplate, and consequently, frequency shift or sensibility efficiency of nanosensor decrease/increases. Furthermore, another remarkable point derived from this table is that the sensitivity property of the nanosensor reduces by increasing the number of added masses. Also, it can be found that frequency shift significantly decreases as the length to thickness ratio increases. Figure 3 gives a comparison of variations of frequency shift against to the total attached mass for different values of the nonlocal parameter, The most important distinctions are as follows. It is observed that by taking into account the large values of the nonlocal parameter, the sensitivity performance of nanosensor reduces, particularly for large values of total attached mass (mainly due to the softening effect of the nonlocal parameter on the nanostructures). Furthermore, it is demonstrated that increasing the total attached mass value from 10 À21 to 10 À19 kg leads to a notable increase in the frequency shift of FG-MEE nanoplate. Moreover, it can be resulted from this figure that the influence of various figurations of attached nanoparticles on the frequency shift become lower for large values of the total attached mass and the nonlocal parameter. It is also found that for large values of nonlocal parameter, FG-MEE nanosensor is insensitive to the variation of the total attached mass.
Depicted in Fig. 4 is the influence of the aspect ratio on the sensing property of FG-MEE nanosensor for different values of the total attached mass, As can be seen, an increase in the aspect ratio causes to monotonically decreasing the frequency shift of FG-MEE nanoplate and the ability of nanosensor to sense attached particles reduces. Also, it is clear that the case of single mass has the maximum value of frequency shift.
To study the effects of the transverse and shear stiffness parameters of Pasternak elastic substrate on the sensitivity performance of FG-MEE nanosensor, Figure 5 is plotted. Here, we did not consider the external electric and external magnetic potentials, and the nonlocal parameter is supposed to be 0.2. In particular, it is easily deduced that by taking into account stiffer polymer matrix, sensitivity property of nanosensor increases. Furthermore, it is noticed that under the maximum value of the total attached mass, the influence of Pasternak foundation on the frequency shift increases. Figure 6 reveals the influence of the applied external electric voltage on the frequency shift curves versus the total attached mass for the three cases of mass, interesting result is perhaps the fact that a rise in the value of the electric voltage form À0:2 to 0:05 V causes a decrease in the sensitivity behavior of FG-MEE nanosensor. This may be explained by the fact that imposing the positive (negative) applied external electric potential, the compressive (tensile) in-plane loads tend to generate in FG-MEE nanosensor, and consequently, frequency shift reduces (increases). Furthermore, it is found from figure that when FG-MEE nanosensor is subjected to the positive external electric potential, the influence of the total attached mass on the sensitivity property remarkably become low. The frequency shift curves of FG-MEE nanosensor with respect to the total attached mass for different values of external magnetic potential are displayed in Fig. 7. Here, we took It is illustrated that contrary to the influence of external electric voltage on the variation of frequency shift, by increasing external magnetic potential from À0:01 to 0:02 A, the sensitivity property of FG-MEE nanoplate according to the mass sensor considerably increases. Also, it should be noted that when FG-MEE nanosensor is subjected to the negative external magnetic potential, the effect of the total attached mass on the sensitivity property significantly become low. Figure 8a-d show the effect of the higher modes on the difference between the natural frequency of FG-MEE nanoplate with and without added nanoparticle corresponding to the different positions of single attached nanoparticle. The results are calculated for It can be seen that the sensitivity property of nanosensor increases when the attached nanoparticle is closer to the center of FG-MEE nanoplate. This implies that the sensitivity property of nanosensor vanishes when added nanoparticle is located at the edges of FG-MEE nanoplate. Moreover, it is denoted that by increasing mode number the frequency shift remarkably increases.

Concluding remarks
The main contribution of the current study was to provide a nanosize mass detection model basis on the vibrating FG-MEE Mindlin nanoplate resting on Pasternak foundation. Also, external electric voltage and external magnetic potential were exerted to the FG-MEE nanosensor. Eringen's differential constitutive law was employed to take into account small-scale effect into the constitutive equations. By using the power law distribution, the effective material properties of rectangular FG-MEE nanoplate change gradually in the direction of nanoplate thickness. Hamilton's principle was utilized in conjunction with nonlocal Mindlin plate theory to derive the size-dependent governing partial differential equations of motion. Afterward, outlined discussions were depicted to demonstrate how to change the sensitivity property of FG-MEE nanosensor by varying the position of attached nanoparticle, the number of nanoparticles, the values of the nonlocal parameter, initial external electric and magnetic potentials, material property gradient index, the dimension of the FG-MEE nanoplate, and Pasternak substrate coefficients. Some noteworthy points derived from numerical results can be listed as K 63 ¼ À J 2 r mp ð Þ 2 À L 2 R 2 r np ð Þ 2 ; Þmp ð Þ; K 65 ¼ À G 22 þ L 2 ð Þ R np ð Þ; K 66 ¼ À 1 r Q 1 mp ð Þ 2 À R 2 r X 1 np ð Þ 2 À P 1 r; K 67 ¼ K 76 ¼ À 1 r Q 2 mp ð Þ 2 À R 2 r X 2 np ð Þ 2 À P 2 r; K 71 ¼ À F 11 mp ð Þ; K 72 ¼ À H 11 R np ð Þ; Þmp ð Þ; K 75 ¼ À H 22 þ k s L 3 ð Þ R np ð Þ; K 7 6 ¼ À 1 r Q 3 mp ð Þ 2 À R 2 r X 3 np ð Þ 2 À P 3 r: