Nonlocal Elasticity Theory for Transient Analysis of Higher-Order Shear Deformable Nanoscale Plates

The small scale effect on the transient analysis of nanoscale plates is studied.The elastic theory of the nano-scale plate is reformulated using Eringen’s nonlocal differential constitutive relations and higher-order shear deformation theory (HSDT). The equations of motion of the nonlocal theories are derived for the nano-scale plates. The Eringen’s nonlocal elasticity of Eringen has ability to capture the small scale effects and the higher-order shear deformation theory has ability to capture the quadratic variation of shear strain and consequently shear stress through the plate thickness. The solutions of transient dynamic analysis of nano-scale plate are presented using these theories to illustrate the effect of nonlocal theory on dynamic response of the nano-scale plates. On the basis of those numerical results, the relations between nonlocal and local theory are investigated and discussed, as are the nonlocal parameter, aspect ratio, side-to-thickness ratio, nano-scale plate size, and time step effects on the dynamic response. In order to validate the present solutions, the reference solutions are employed and examined. The results of nano-scale plates using the nonlocal theory can be used as a benchmark test for the transient analysis.


Introduction
Due to their superior properties, nanoscale plates such as graphene [1], the two-dimensional (2D) counterpart of threedimensional (3D) graphite, have attracted attention within the fields of solid-state physics, materials science, and nanoelectronics.Experimentation with nano-scale size specimens is both difficult and expensive.Development of appropriate mathematical models for nanostructures, therefore, is essential.Approaches to the modeling of nanostructures are classified into three main categories: atomistic [2,3], continuum [4][5][6], and hybrid atomistic-continuum mechanics [7][8][9].The continuum mechanics approach is less computationally expensive and, furthermore, generates results that tend to be in good agreement with those of the atomistic and hybrid approaches [10,11].
As small size analysis using local theory over-predicts results, correct prediction of micro-/nanostructures requires consideration of small-scale effect.Peddieson et al. [12] applied nonlocal elasticity to the formulation of a nonlocal version of the Euler-Bernoulli beam model and concluded that nonlocal continuum mechanics could potentially play a useful role in nanotechnology applications.One of the wellknown continuum mechanics theories that include small scale effects with good accuracy is Eringen's nonlocal theory [13,14].This, compared with classical continuum mechanics theory, can predict the behavior of the large nanosized structures, without a large number of equations.
Nonlocal elasticity theory has been widely applied to the bending, vibration, and buckling behavior of one-dimensional (1D) nanostructures including nanobeams, nanorods, and carbon nanotubes.Contrastingly, there has been no such work related to the transient analysis of nano-scale plate based on higher-order shear deformation theory (HSDT).However, some researchers have applied nonlocal elasticity theory to the study of various applications of micro-and nanostructures [15,16].Understanding the importance of employing nonlocal elasticity theory to small scale structures, 2 Journal of Nanomaterials a number of researchers have reported on static, dynamic and stability analysis of 1D and 2D micro-/nanostructures [17][18][19][20][21][22][23][24][25].Kiani [26] studied the free vibration of embedded single-walled nanotubes accounting for nonlocal effect using a meshless method.Nonlocal continuum-based modeling of a nanoplate subjected to a moving nanoparticle was presented by Kiani [27,28].Jomehzadeh and Saidi [29] investigated the decoupling of the nonlocal elasticity equations for three dimensional vibration analysis of nanoplates.
These studies [17][18][19][20][21][22][23][24][25] were based on classical and firstorder theories of plates.The analysis of 2D nano-scale plates with accurate stress fields requires further study.The thirdorder shear deformation theory (TSDT) of Reddy [30] is based on a displacement field that includes the cubic term in the thickness coordinate; hence, the transverse shear strain and stress are represented as quadratic through the plate thickness and vanish on the bounding planes of the plate.Some studies, in fact, have incorporated the TSDT to obtain more accurate results [31][32][33].It seems appropriate and potentially fruitful, therefore, to explore TSDT's extension to include size effects.The present study, accordingly, applies the nonlocal TSDT to the transient dynamic response of nanoscale plates.
To avoid the resonant behavior of such structures, a vibration analysis of their design is very important.Also, most of these structures, whether they are used in the civil, marine or aerospace field, are subjected to dynamic loads; therefore, there is a need to assess their transient response.However, transient analysis of nano-scale plates has not received adequate consideration.Work that exists on the vibration analysis of laminated composite and FGM plates can be found in [34,35].
In the present work, nonlocal elasticity theory is applied to a transient analysis of nano-scale plates.Navier's method is used to solve the governing equations for simply supported boundary conditions.The effects of (i) the nonlocal parameter, (ii) size of the nano-scale plates, (iii) aspect ratios, and (iv) side-to-thickness ratios on the nondimensional dynamic responses are investigated.The present work will be helpful to the design of nano-electro-mechanical system or microelectro-mechanical systems devices incorporating nano-scale plates.

Formulation
In classical local elasticity theory, stress at a point depends only on the strain at that point, whereas, in nonlocal elasticity theory, the stress at a point is a function of the strains at all points on the continuum.In other words, in nonlocal elasticity theory, stress at a point is determined by both stress at that point and spatial derivatives of it.According to Eringen [13,14], the nonlocal constitutive behavior of a Hookean solid is represented by the differential constitutive relation where   is the local stress tensor,   is the nonlocal stress tensor, and the nonlocal parameter  is defined by in which  0 is the material constant set by the experiment and  is the internal characteristic length.The relations between stress resultants in local theory and in nonlocal theory are defined by integrating plate thickness into (1) as where L = 1 − ∇ 2 and where  and  take the symbols , , the superscript  denotes the quantities in local third-order shear deformation theory, and ℎ is the thickness of the plate.
In general, differential operator ∇ in (1) is the 3D Laplace operator.For 2D problems, the operator ∇ may be reduced to the 2D one.Thus, the linear differential operator L becomes wherein it is clear that operator L is independent of the  direction.
The third-order shear deformation theory (TSDT) extends the first-order theory by assuming that shear strain and consequently shear stress are not constant through the plate thickness [30].The displacement field of the third-order theory of plates is given by where  0  and  0 are the in plane displacements and transverse displacement of point on the midplane (i.e.,  = 0), respectively,  1 = 4/3ℎ 2 ,   is the transverse displacement of the mid-plane of the plate, and   denotes the slope of the transverse normal on the mid-plane.
The variational statements facilitate the direct derivation of the equations of motion in terms of the displacements.The governing equations of the third-order nonlocal plate theory can be derived, using the dynamic version of the principle of virtual displacement (Hamilton's principle), as where  is the virtual strain energy,  is the virtual work done by external forces, and  is the virtual kinetic energy.
It can be seen from ( 1) that nonlocal behavior enters into the problem through the constitutive relations.The principal of virtual work, because it is independent of such relations, can be applied to derive the equilibrium equations of nonlocal nano-scale plates.On the basis of the principle of virtual displacements, then the governing equations are obtained [32] where ) and  2 = 3 1 , in which  is the mass density.From these, by setting L = 1, the equations of motion of the conventional third-order plate theory are obtained.

Navier Solutions of Nonlocal Third-Order Shear Deformation Theory
Here, analytical solutions to the free and forced vibration of simply supported plates are presented using the nonlocal third-order plate theory in order to illustrate the small scale effects on deflections and natural frequencies of the plates.For the set of simply supported boundary conditions, the Navier solution can be obtained [32].According to Navier solution theory, the generalized displacements at middle of the plane ( = 0) are expanded in double Fourier series as where Λ 1 = cos  sin  ⋅     , Λ 2 = sin  cos  ⋅     , and Λ 3 = sin  sin ⋅    , in which,  = /,  = /,  = √ −1,  is the time and   is the natural frequency.Substituting ( 9) into ( 8), the matrix form obtained is where {Δ} = {  ,   ,   ,   ,   }, [K] is the stiffness matrix, [M] is the mass matrix, dot-superscript convention indicates the differentiation with respect to the time variable , and {} is the force vector.For this study, the Newmark- time integration method is used.In order to investigate the transient response of a laminated composite and nano-scale plate, a parametric study is carried out according to different dynamic loading conditions.The load acting on the plate is distributed over its entire surface and varied with time according to  (, , ) =    () , (11) where   = ∑ ∞ =1 ∑ ∞ =1   sin  sin , the   is the coefficient for some typical load [32], () = {1, 0 ≤ }: rectangular step loading, and 1, 0 ≤  ≤  1 : rectangular pulse loading, 0,  ≥  1 (12) in which  1 = 1000 s.

Numerical Results and Discussion
To validate this present study, several numerical examples are solved to test the performance in free and forced vibration analyses.Examples include composite materials for confirmation of crucial features or for comparison with previous published analysis results.Specially, a mistake in a reference is corrected for the purpose of an accurate comparison.

Validation.
In order to validate the results of the forced vibration analysis, a transient analysis of the simply supported cross-ply (0 ∘ /90 ∘ ) plate shown in Figure 1 is carried out under sinusoidal step loading intensity  0 = 10 with Δ = 5 s.The rectangular plate, shown in Figure 1, is analyzed for the natural vibration of laminated composite plate, and results are presented in Figure 2. The geometric and material properties applied are as follows: The center deflection versus time of a cross-ply plate (0 ∘ /90 ∘ ) is compared with the results by Reddy [32]. Figure 2 shows the present results and the results by Reddy [32].It can be seen that the difference between the two results is almost negligible.
The results are presented in the nondimensional form using the equation: The results are presented in the non-dimensional form using the equation: The results based on HSDT plate theory with various values of nonlocal parameter  are compared with those reported by Aghababaei and Reddy [36].It can be seen that when the aspect ratios are 1 and 2, the two sets of results are identical.Table 1 also shows that the material property ( = 30×10 6 ) presented by Aghababaei and Reddy [36] is not accurate, and that, for accurate comparison by (15), it should be modified.

Forced Vibration of Nanoscale Plate under
Step Load.The non-dimensional geometric and material data on the nanoscale plate are provided in Table 2.
The results of a forced vibration analysis of laminated composite plate under a distributed rectangular step load with Δ = 5 s and a variable nonlocal parameter are plotted in Figure 3.The nonlocal parameters are  = 0, 0.25, 1.0, 2.25, and 4, because  0  in (2) should be smaller than 2.0 nm for carbon nanotubes (Q.Wang and C. M. Wang [37]).The results show the profound scale effect for the higher values of nonlocal parameter.In other words, the local elasticity theory under-predicts the dynamic response.Thus also, if local elasticity results are used for experimental determination of Young's modulus, the value will be underpredicted.

Forced Vibration of Nanoscale Plate under Pulse Load.
A forced vibration analysis of center deflection under the rectangular pulse loading is carried out using the same material.The plot of the center deflection versus time is shown in Figure 4.As can be seen, the maximum deflection increases until the applied loading is removed after  1 = 1000 s, due to the nonlocal parameters increase due to stiffness degradation.After removal of the load, however, a change of the maximum deflection occurs.This phenomenon is expected to vary with the point at which the load is removed.Accordingly, the amplitude also changes after the load removal.
Figure 5 shows the variation of dynamic response with respect to the length of the nano-scale plate, which varies from 5 nm to 20 nm.The nonlocal variable is assumed to be 1.As the size of the nano-scale plate increases, the amplitude and frequency of the dynamic response decrease.The dynamic response due to the change of the aspect ratio (/ = 0.5, 0.75, 1.0) of the nano-scale plate is shown in Figure 6.The nonlocal variable is again assumed to be 1.The amplitude and the period of the dynamic response decrease, as the aspect ratio decreases.When the aspect ratio changes from 1 to 0.75, the maximum deflection is reduced by 45%: when the aspect ratio changes from 0.75 to 0.5, the reduction is 64%.
Figure 7 shows the dynamic response of a nano-scale plate of variable side-to-thickness ratio under a suddenly applied pulse load.Typically in the case of isotropic plates or laminated composite plates, the side-to-thickness ratio is more than 20, and the effect of shear deformation is very small.Similarly, in the case of the nano-scale plate shown in Figure 7 with a side-to-thickness ratio of 20 or greater, there is no effect on the maximum deflection.But when the side-to-thickness ratio decreases, the frequency of the dynamic response is increased; though there is no difference in its amplitude.Notably, the frequency with a the side-tothickness ratio of 50 is twice that with a the side-to-thickness ratio of 100, due to the doubling of the thickness of the nanoscale plate.
Next, the effect of the time step on the accuracy of the solutions is investigated using a nano-scale plate under distributed pulse loading.Figure 8 shows the nondimensionalized center transverse deflection for three different time steps.A larger time step reduces the amplitude and increases the period.For all of the time steps below 10 s, the effect is not noticeable on the graph.In all the subsequent, cases then the time step = 5 s is used.The difference between the maximum deflection of time step 5 s and that of 50 s is about 3%.It is expected therefore that, in order to obtain an accurate dynamic response, an appropriate choice of the loading time interval is required.Figure 9 shows the dynamic response of the nano-scale plate under a suddenly applied pulse load and a variable time step.The change of the amplitude occurred according to the time in which the load is removed.
In order to study the effect of nonlocal parameter on the maximum deflection, the predicted results of the present nonlocal plate model have been illustrated in Figure 10.As it is obvious from Figure 10, the predicted results of nonlocal plate model generally increase with the nonlocal parameter, irrespective of the aspect ratio (/).As the nonlocal parameter increases, the ratios ( (/=1.0)/ (/=0.5)and  (/=0.75)/ (/=0.5)) of the proposed nonlocal plate decrease.In other words, the effect of nonlocal parameters on maximum deflection is more significant when the aspect ratio is 0.5 (/ = 0.5).
Figure 11 plots the absolute value of maximum deflection of the nano-scale plate under suddenly applied step and pulse load with variable  1 .The maximum deflection when changing from 710 to 940 shows a negative value, and when changing from 950 to 1180, it shows a positive value.It would be expected that this absolute value of maximum deflection curve maintains a constant shape and that the maximum value of each deflection curves does not exceed the maximum deflection when subjected to a step load.

Conclusions
In this paper, equations of motion under nonlocal thirdorder shear deformation theory (TSDT) as well as analytical solutions for bending and free vibration are derived to elucidate the effect of nonlocal parameters.In order to illustrate the effects of nonlocal theories on plate response, numerical results for simply supported rectangular plates are presented.It is found that nonlocal theory can be applied to the analysis of nano-scale plates where the small size effects are significant.
Using the Eringen's nonlocal elasticity theory, the dynamic responses of nano-scale plates are studied.The higher-order shear deformation theory (HSDT) is applied to capture the quadratic variation of shear strain and consequently shear stress as function of plate thickness.It is determined that the effect of nonlocal constitutive relations is to increase the magnitude of deflections and decrease the amplitude of frequencies.Additionally, it is found that the difference between nonlocal theory and local theory is significant for high values of nonlocal parameter.
From the present work, the following conclusions are drawn.
(1) The dynamic response computed by nonlocal elasticity is always larger compared with the local elasticity value.The stiffness of a nano-scale plate by nonlocal elasticity is small compared with that by local elasticity.Correspondingly, the amplitude and period are large.
(2) As the size of a nano-scale plate decreases, the effect of nonlocal elasticity becomes more significant and dynamic response increased.With aspect ratio decreases of / = 1.0 to 0.75 and 0.75 to 0.5, the amplitude and frequency of the dynamic response decrease from 1.0 to 0.75, the maximum deflection is reduced by 45%, and from 0.75 to 0.5, there is a 64% reduction.
(3) In the case of a nano-scale plate of side-to-thickness ratio larger than 20, the effect of shear deformation is very small.When the side-to-thickness ratio decreases, there is no difference in the amplitude of the dynamic response, but its frequency increases.
(4) The effect of a larger time step is to reduce the amplitude and increase the period.It is expected therefore that, in order to obtain accurate dynamic response, an appropriate choice of loading time interval is required.
The dynamic response of the nano-scale plate can serve as a benchmark for future guidelines to nano-scale structural design under nonlocal elasticity theory.Especially, the maximum deflection of a nano-scale plate under a suddenly applied step and pulse load with variable  1 will be contributed for the purpose of a future benchmark test applicable to the transient analysis of nano-scale structures.Further, it will be necessary to include, in transient analysis, the nonlocal elasticity theory for nano-scale shell and other boundary conditions.

4 Figure 3 : 4 Figure 4 :
Figure 3: Dynamic response of nano-scale plate under suddenly applied step load with variable nonlocal parameter.

Figure 10 :
Figure 10: Maximum deflection of nano-scale plate of various aspect ratio with nonlocal parameter.

1 Figure 11 :
Figure 11: Absolute value of maximum deflection of nano-scale plate under suddenly applied step and pulse load with variable  1 ( = 1).

Table 1 :
Nondimensional frequency of simply supported nanoscale plate with various nonlocal parameter ().
) 4.2.Free Vibration of Nanoscale Plate.The analytical free vibration responses are numerically evaluated here for an

Table 2 :
Nondimensional material and geometric properties.
isotropic plate to discuss the effects of nonlocal parameter  on the plate vibration response.Table1shows the nondimensional frequencies of simply supported nano-scale plates with various values of nonlocal parameter .The nanoscale plate is made of the following material properties: