Next Article in Journal
Theoretical Investigation of a Simple Design of Triple-Band Terahertz Metamaterial Absorber for High-Q Sensing
Next Article in Special Issue
Size-Dependent Free Vibrations of FG Polymer Composite Curved Nanobeams Reinforced with Graphene Nanoplatelets Resting on Pasternak Foundations
Previous Article in Journal
Fertility Detection of Hatching Eggs Based on a Convolutional Neural Network
Previous Article in Special Issue
Modified Nonlocal Strain Gradient Elasticity for Nano-Rods and Application to Carbon Nanotubes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Buckling and Free Vibrations of Nanoplates—Comparison of Nonlocal Strain and Stress Approaches

Institute of Machine Design, Cracow University of Technology, ul. Warszawska 24, 31-155 Kraków, Poland
*
Author to whom correspondence should be addressed.
Appl. Sci. 2019, 9(7), 1409; https://doi.org/10.3390/app9071409
Submission received: 27 February 2019 / Revised: 26 March 2019 / Accepted: 1 April 2019 / Published: 3 April 2019
(This article belongs to the Special Issue Recent Advances in Non-Local Modelling of Nano-Structures)

Abstract

:
The buckling and free vibrations of rectangular nanoplates are considered in the present paper. The refined continuum transverse shear deformation theory (third and first order) is introduced to formulate the fundamental equations of the nanoplate. Besides, the analysis involve the nonlocal strain and stress theories of elasticity to take into account the small-scale effects encountered in nanostructures/nanocomposites. Hamilton’s principle is used to establish the governing equations of the nanoplate. The Rayleigh-Ritz method is proposed to solve eigenvalue problems dealing with the buckling and free vibration analysis of the nanoplates considered. Some examples are presented to investigate and illustrate the effects of various formulations.

Graphical Abstract

1. Introduction

Concerning nano-objects, it is observed that classical (local) continuum theories presenting scale-free relations are inadequate for the mechanical analysis because of the absence of any material length scale parameters, see e.g., works by Peddieson et al. [1], Lu et al. [2], Barretta et al. [3]. At the nanometer scale, the size effects often may become prominent. Some non-classical continuum approaches have been applied to describe the mechanical behavior of small-sized structures, such as couple-stress theory, e.g., papers by Chen and Li [4], Wang [5], Ma et al. [6], Kim et al. [7], nonlocal elasticity theory, e.g., works by Barretta et al. [8], Hu et al. [9], Arefi et al. [10], Arash and Wang [11], strain gradient theories, e.g., papers by Askes and Aifantis [12], Akgoz and Civalek [13], Lim et al. [14], Barretta et al. [15].
A continuum is described as classical if the displacement field completely defines the deformation process. So, every point of a classical continuum (or Cauchy-type) has three displacement degrees of freedom. In the classical theory of elasticity, the symmetric tensor of stress and strain exists. The classical continuum mechanics is widely used for many problems. However, discrepancies between theory and experiment were observed around holes and notches when the stress gradient is present [16]. The classical elasticity also fails in the vibration analysis, wave propagation etc., where the influence of the material microstructure is dominant. The same problem is connected with the mechanical analysis of granulated materials, thermosetting polymers, and nanomaterials.
Firstly, the non-classical continuum was discussed at the end of the nineteenth century, e.g., by Voigt (1887) [17]. At the beginning of the twenty century, the Cosserat brothers (1909) [18] formulated the couple-stress continuum theory to overcome the insufficiency of the classical approach. The Cosserat continuum model besides the ordinary (force-stresses) considers the couple-stresses. Hence, the displacement degrees of freedom and independent, rotational degrees of freedom are assigned to every particle of the Cosserat continuum. The stress tensor in Cosserat medium is non-symmetrical. The Cosserat theory in the primary form did not involve the broad public. It was excessively complex to be applied because of its generality and geometrical non-linearity [19]. In the following years, the activities were arising on an extension of the couple-stress theory, e.g., by Toupin [20] and Mindlin and Tiersten [21]. In the 1970s and 1980s, Eringen and Eringen and Edelen [22,23,24] derived stress gradient theory, and then Aifantis [25,26] elaborated strain gradient theory. These gradient theories have a reduced amount of higher order terms [7,27], are much simpler than those proposed in the 1960s. Mindlin and co-workers [21] made some assumptions to formulate a much simpler version of the Cosserat theory. They significantly enriched the Cosserat theory by introducing the linearised form of the couple-stress theory for perfectly elastic, centrosymmetric-isotropic materials. They introduced an internal length parameter which brings equations with and without couple-stresses. The parameter is small in comparison with body dimensions, but its influence may become important as dimensions of a body or phenomenon are in the order of the length parameter. Mindlin’s theory of elasticity incorporates kinematic quantities at macro-scale and micro-scale so it is a multi-scale description [28]. Mindlin and co-workers adopted the couple-stress theory to solve problems of wave propagation, vibration, stress-concentration and nuclei of strain. In the following years, the gradient theories of elasticity were developed having the additional higher order spatial derivatives of stress, strain or acceleration. New generations of gradient theories are focused on minimisation of the number of additional constitutive parameters [12].
Eringen and co-workers [22,23,24] introduced the stress gradient theory from the integral nonlocal theories. They formulated that the stress state at a given point x depends not only on the local stress at that point but is the weighted average of the local stress of all points in the neighbourhood of x. A nonlocal stress σ i j N L o c is defined as [24]:
σ i j N L o c ( x ) = α ( | x x | , τ ) σ i j ( x ) d V ( x )
where σ i j N L o c is nonlocal stress tensor at point x, σ i j is classical (local) stress tensor at x′ within domain V, α is nonlocal modulus that depends on the distance x x and on the dimensionless length scale τ expressed as τ = e 0 θ / L with e 0 θ nonlocal parameter, e 0 nonlocal material constant, θ and L internal and external characteristic lengths, respectively.
For particular nonlocal problems originally studied in [24], instead of the nonlocal integral Equation (1), Eringen formulated an equivalent differential law as:
σ i j N L o c e 0 θ 2   2 σ i j N L o c = σ i j
where 2 is the Laplace operator.
The problem of the equivalence between the original Eringen integral convolution Equation (1) and the differential form Equation (2) and associated mechanical paradoxes are discussed in the literature—see [1,29,30] for the detailed description.
In the 1984 Aifantis [31] formulated a simple model of strain gradient plasticity by introducing only one additional constitutive constant. Firstly, the strain gradient theory was applied to dislocation description. Next, in the 1990s, Aifantis and co-workers [12,25,26,32,33] presented another simple model for strain gradient elasticity having one additional constitutive constant based on the strain gradient plasticity. They proposed gradient elasticity theory for finite deformations and infinitesimal deformations and extended the linear elastic constitutive relations with the Laplacian of the strains as [12]:
σ i j = C i j k l ( ε k l ( e 0 θ ) 2 2 ε k l )
Recent achievements in mechanics of generalized continuum inform that deformation measures describing the generalized continuum, such as Cosserat theories, are kinematically redundant. The above original concept and the detailed discussion were first introduced in the paper by Romano et al. [34] and next referenced in [35,36,37].
Applications of nonlocal continuum mechanics in case of carbon nanostructures such as carbon nanotubes and graphene were studied by [38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55] among others. The nonlocal vibrations and buckling of graphene were anlyzed by Pradhan [39], Pradhan and co-workers [40,41] and Shakaee-Pour [42]. Aghababaei and Reddy [43] used nonlocal shear deformation theory to investigate the bending and vibration of plates. Fazelzadeh and Ghavanloo [45] applied the nonlocal anisotropic elastic shell model for vibrations of nanostructures. Nami and Janghorban [47] studied the resonance of functionally graded micro/nano plates involving nonlocal stress and strain gradient elasticity. Ansari et al. [48] applied the first order shear deformation theory and the modified couple-stress theory to analyze the vibration of microplates. Barretta and Marotti de Sciarra [49] presented the nonlocal elastostatic problem and pointed out the nonlocality effect. Sladek et al. [51] applied the nonlocal theories for a large deformation of piezoelectric nanoplates. Zenkour [53] analyzed the thermoelastic vibrations of nanoplates based on the nonlocal theory of elasticity. Currently, a monographic elaboration of nonlocal elasticity approaches is presented by Ghavanloo et al. [56]. The problem of mechanical properties of nanostructures and nanocomposites in local and nonlocal regime was also previously considered by authors [57,58,59,60,61,62,63,64,65,66].
The main aim of the current paper is to investigate the effect of various formulations on the description of buckling and free vibrations of a rectangular nanoplate. The nonlocal stress and strain theories are involved, and the influence of small-scale parameter on the mechanical behavior of the nanoplate is studied and discussed.

2. Formulation of the Problem

Let us consider a nanoplate of length Lx, width Ly and thickness t as shown in Figure 1. A Cartesian coordinate system (x, y, z) is used to label the material points of the nanoplate in the unstressed reference configuration. The transverse normal stress σzz is negligible in comparison with in-plane stresses σxx, σyy and transverse out of plane stresses σxz, σyz. In addition to the vibration analysis, the elastic buckling of nanoplate under uniform uniaxial compressive loadings as well as biaxial loadings ( N x 0 , N y 0 ) are considered in this work.

2.1. Kinematic Relations

The 3-D displacements U 1 , U 2 and U 3 at any point of the nanoplate in the x, y and z directions (see also Figure 1), respectively, are expressed in the following form:
U 1 ( x , y , z ) = u ( x , y ) z w x + F z ψ 1 ( x , y ) U 2 ( x , y , z ) = v ( x , y ) z w y + F z ψ 2 ( x , y ) U 3 ( x , y , z ) = w x , y
where u, v, w are the displacements of a generic point on the reference mid-surface, ψ1, ψ2 are the rotations of normal to the mid-surface about the y- and x-axes, respectively. The explicit form of the function F(z) is presented in Table 1 and Figure 2.
The distributions of the transverse shear corrections functions F(z) used in the literature [67,68] are demonstrated in Figure 2 for third order transverse shear deformation theory (TSDT). The differences in the values of the functions are negligibly small. In HSDT there is no need to use shear correction factors since the form of the functions leads to the parabolic distributions of transverse shear stresses.
Based on the above expressions, the total, linear 3-D strain tensor can be written as follows:
ε i j x , y , z = 1 2 U i x j + U j x i ,   x 1 = x ,   x 2 = y , x 3 = z ,   i , j = 1 , 2 , 3

2.2. Constitutive Nonlocal Equations and Variational Formulation

The brief overview of nonlocal theories has been presented in the Introduction. We have focused on the classical (local) and modified couple-stress formulations as well as the stress and strain gradient theories. Now, on the base of the local continuum, the constitutive nonlocal equations for the nanoplate are presented. Generally, both strain and stress gradients can influence the constitutive response of materials.
The local (classical) constitutive relations for a shear deformable nanoplate (orthotropic) are in the following form:
  • local stress theory
    { σ 11 σ 22 σ 12 σ 23 σ 13 } = [ C 11 C 12 0 0 0 C 12 C 22 0 0 0 0 0 C 66 0 0 0 0 0 C 44 0 0 0 0 0 C 55 ] × { ε 11 ε 22 ε 12 ε 23 ε 13 }
  • local strain theory
    { σ 11 σ 22 σ 12 σ 23 σ 13 } = [ C 11 C 12 0 0 0 C 12 C 22 0 0 0 0 0 C 66 0 0 0 0 0 C 44 0 0 0 0 0 C 55 ] × ( 1 ( e 0 θ ) 2   2 ) { ε 11 ε 22 ε 12 ε 23 ε 13 }
    2 represents the classical 2-D Laplacian operator.
In the Cartesian system of coordinates the nonlocal functional of the strain energy of orthotropic plates becomes:
U L = 1 2 Ω t / 2 t / 2 σ N L o c ε d x d y d z
where Ω denotes the space occupied by the nanoplate, { ε } is the local strain tensor and { σ N L o c } is the nonlocal stress tensor. The nonlocal stress tensor { σ N L o c } is defined as follows:
  • the nonlocal stress theory
    { σ } = [ C ] { ε } = ( 1 ( e 0 θ ) 2   2 ) { σ N L o c } =   { σ N L o c }
  • the nonlocal strain theory
    { σ N L o c } = [ C ] ( 1 ( e 0 θ ) 2   2 ) { ε }
    where { σ } is the local stress tensor, [ C ] is the stiffness matrix and is an operator.
The components of the strain tensor are functions of the kinematic variables (4), (5) so that inserting the definitions (9) or (10) in Equation (8) one can find that the functional U L becomes a function of the kinematic variables u, v, w, ψ1, ψ2 in the right-hand side of (4) only.
Besides, the work done by the exerted loads on the nanoplate is as follows (see Figure 1):
W = 1 2 Ω N x 0 w x 2 + N y 0 w y 2 d x d y
Finally, the functional of the total potential energy in the Lagrange form can be written in the following form:
Π L = U L + W
The kinetic energy T of the analyzed nanoplate can be expressed as:
T = 1 2 Ω t / 2 t / 2 ρ ( U ˜ i τ ) 2 d x d y d z
where ρ and τ denotes the nanoplate density and time, respectively, and i = 1, 2, 3.
When the nanoplate is oscillating in a normal mode, the motion is harmonic so that the solutions may be taken in the form of a harmonic function of time τ with a natural circular frequency ωmn (m, n denote the modes of vibration), such as
{ U ˜ i } ( x , y , z , τ ) = { U i } ( x , y , z ) ( c o s ω m n τ )
where { U i } and { U ˜ i } are the components of displacements dependent and independent on time τ , respectively.
The Hamilton principle is applied to derive governing equations which is defined in a classical way as energy functional:
H = U L + W + T
The Hamilton’s principle states that
δ τ 1 τ 2 H d τ = 0
where δ denotes the variation and τ 1 and τ 2 denotes the prescribed time range.
For the simplicity of the further considerations, the Hamilton function is replaced by the more convenient form in the case of the nonlocal stress theory, i.e.,:
H ˜ = H ,   where   = 1 e 0 2 θ 2   2
using the definition (9) and assuming that 1   is equal to the unity operator.
As it may be seen for CPT and FSDT/HSDT the governing relations of motion derived from the first variations of the Hamilton functional H ˜ are reduced to the addition of the linear operator acting on the right side of equilibrium equations. The stress resultants defined for the local stress components are written as:
N i j = t / 2 t / 2 σ i j z d z ,   Q i 3 = 5 6 σ i 3 t ,   M i j = t / 2 t / 2 σ i j z 2 d z ,   i , j = x , y
where Nij, Mij and Qi3 are the resultant forces, moments and shear forces.
The sense of the transformation (17) is obvious by analysis of the relations written in Table 2.

3. The Rayleigh-Ritz Method

The Rayleigh quotient is involved in the min-max theorem to obtain exact values of all eigenvalues (see the analysis presented by Muc et al. [69,70]. It can also be applied in eigenvalue algorithms to get an eigenvalue approximation from an eigenvector approximation. Exactly, this is the base for Rayleigh quotient iteration. A careful selection for the trial function that fulfills boundary conditions and depends on the set of variational parameters must be given in advance. Therefore, the Rayleigh-Ritz variational principle is a powerful tool method for the approximate solution of eigenvalue problems where a trial function (or functions) is introduced.
For simply-supported boundary conditions the displacement trial functions u(x,y), v(x,y), w(x,y), ψ1(x,y) and ψ2(x,y) (Equation (4)) are expressed in terms of Fourier functions that satisfy the geometric and stresses boundary conditions. Thus the essential boundary conditions are satisfied, and the Fourier polynomials provide a complete and orthogonal set leading to a relatively fast convergence. The trial functions are defined in the whole space Ω occupied by the nanoplate.
Assuming the buckling (free vibrations) displacements of the form:
u ( x , y ) = m , n = 1 U m n cos ( α m x ) s i n ( β n y ) v ( x , y ) = m , n = 1 V m n sin ( α m x ) cos ( β n y ) w ( x , y ) = m , n = 1 W m n sin ( α m x ) cos ( β n y ) ψ 1 ( x , y ) = m , n = 1 Ψ 1 m n cos ( α m x ) sin ( β n y ) ψ 2 ( x , y ) = m , n = 1 Ψ 2 m n sin ( α m x ) cos ( β n y ) where   α m = m π L x   ,   β n = m π L y   ,   m ,   n natural   wavenumbers .
The explicit form of the Rayleigh-Ritz quotient is given below:
  • the nonlocal stress theory
    λ = Ω t 2 t 2 d x d y d z ( σ 11 ε 11 + σ 12 ε 12 + σ 22 ε 22 + σ 13 ε 13 + σ 23 ε 23 ) ( W + T ) ,   { σ } = [ C ] { ε }
  • the nonlocal strain theory
    λ = Ω t 2 t 2 d x d y d z ( σ 11 ε 11 + σ 12 ε 12 + σ 22 ε 22 + σ 13 ε 13 + σ 23 ε 23 ) ( W + T ) , { σ } = [ C ] ( 1 ( e 0 θ ) 2   2 ) { ε }
    where λ is the Rayleigh quotient.
The Equations (20) and (21) are written in the Lagrange form since the strains are functions of the kinematical variables u, v, w, ψ1, ψ2 reduced to three for CPT.

4. Numerical Results and Discussion

We consider the structural setup appearing in Figure 1. We likewise investigate how different constitutive formulations and the introduced kinematical hypotheses for the nanoplate (here a graphene sheet) affect the solutions. The properties of the graphene sheets are assumed as follows: the thickness t = 0.34 nm, the length of the nanoplate Lx = Ly = 10 nm, Young’s modulus E = 1.06TPa and Poisson’s ratio ν = 0.25. However there are many articles that deal with the quantitative estimation of characteristics of nanoscale materials using local/nonlocal continuum mechanics with well-established parameters fitted to experimental results—see, e.g., papers by Shi et al. [71,72,73].
For the simplicity of writing the results in a more compact manner let us assume that:
C 11 = C 22 = E 1 ν 2 , C 12 = ν C 11 , C 44 = C 55 = C 66 = 1 ν 2 C 11
The nanoplates are subjected to uni-axial uniform compression loading (see Figure 1).
The numerical analysis of buckling loads and natural frequencies is conducted with the use of the symbolic package Mathematica. The program allows us to integrate symbolically overall geometrical variables x, y and z appearing in the functionals (20) and (21). Since for CPT and FSDT the equilibrium equations are presented in the explicit form in Table 2, the appropriate buckling loads and natural frequencies can be evaluated almost automatically inserting the Fourier expansions (19) to Equations (20) and (21). Let us also note that in plate equations (Equation (4)) the terms including in-plane displacements, the normal deflection and two rotations are uncoupled. Therefore, the governing relations are a set of three equations only—see Table 2.

4.1. Buckling Loads–Uniform, Uni-Axial Compression along the x-axis ( N y 0 = 0 )

4.1.1. Nonlocal Stress Theory

Taking into account the uniform, uni-axial compression loading along the x-axis (assuming N y 0 = 0 ) the following relation for the nonlocal stress theory are considered:
  • CPT:
    N x 0 = D ( α m 2 + β n 2 ) 2 α m 2 [ 1 + e 0 2 θ 2 ( α m 2 + β n 2 ) ]
  • FSDT and HSDT:
    The determinant of the coefficient matrix derived from Equation (20) must be equal to zero, i.e.,:
    det K 33 K 33 K 33 K 43 K 44 K 45 K 53 K 54 K 55 = 0
The determinant is cut off to three coupled displacement terms (3–w, 4–ψ1, 5–ψ2) and in this way, the coefficients vary in the range 3–5.
K 33 F S D T = ( α m 2 + β n 2 ) N x 0 t C 44 α m 2 [ 1 + e 0 2 θ 2 ( α m 2 + β n 2 ) ] ,   K 34 F S D T = α m ,   K 35 F S D T = β n K 43 F S D T = K 34 F S D T , K 44 F S D T =   1 D t C 44 α m 2 1 + 1 ν 2 , K 45 F S D T = D t C 44 α m β n 1 + ν 2 , K 53 F S D T = K 35 F S D T ,   K 54 F S D T = K 45 F S D T ,   K 55 F S D T =   1 D t C 44 β n 2 1 + 1 ν 2 , K 33 H S D T = D ( α m 2 + β n 2 ) 2 N x 0 α m 2 [ 1 + e 0 2 θ 2 ( α m 2 + β n 2 ) ] K 34 H S D T = ξ 1 α m α m 2 + β n 2 , K 35 H S D T = ξ 1 β n α m 2 + β n 2 , K 43 H S D T = K 34 H S D T , K 44 H S D T = ξ 2 α m 2 + 1 ν 2 β n 2 + ξ 3 , K 45 H S D T = ξ 2 1 ν 2 α m β n , K 53 H S D T = K 35 H S D T , K 54 H S D T = K 45 H S D T , K 55 H S D T = ξ 2 1 ν 2 α m 2 + β n 2 + ξ 3
D = E t 3 12 1 ν 2 , ξ 1 = E 1 ν 2 t / 2 t / 2 z F ( z ) d z , ξ 2 = E 1 ν 2 t / 2 t / 2 F 2 ( z ) d z , ξ 3 = E 2 1 + ν t / 2 t / 2 d F ( z ) d z 2 d z
Computing the values of the parameters ξ1, ξ2, ξ3 one can verify easily that for the third order polynomial and exponential approximations of the functions F(z) (see Equation (4), Table 1 and Figure 2) the differences of their values are negligibly small.
In order to compare various effects of constitutive relation formulations (the nonlocal stress and strain theory) as well as of kinematical relations (Equation (4)) the buckling loads will be further referred to the value:
N ¯ = D α m 2 + β n 2 2 α m 2
It corresponds to the buckling loads evaluated with the use of the classical plate theory. In addition, the results of the analysis carried out in Refs [74,75,76] show that for the Lx/Ly > 1 the buckling mode n = 1 and m ≥ 1. The presented plots are just restricted to the abovementioned limits.
Figure 3 demonstrates the classical results—the formulations of 2-D kinematical relations have a significant effect on the values of buckling loads which increase with the growth of the t/Lx ratio.
The decrease of buckling loads for more accurate kinematical formulations in the thickness (z-direction) is generally accepted. However it is essential to point out that the percentage error between classical and higher order plate theories is also a function of mechanical properties of nanoplates, particularly in view of the values of Gxz/E and Gyz/E ratios, where Gxz and Gyz are the transverse (out-of-plane) shear moduli—see, e.g., Muc [76].
As it may be seen from Equations (23)–(26) the nonlocal stress theory formulation is strictly connected with the existence of the multiplication factor 1 / [ 1 + e 0 2 θ 2 ( α m 2 + β n 2 ) ] . The influence of the constitutive relations of nonlocal stress theory on the buckling loads is plotted in Figure 4.

4.1.2. Nonlocal Strain Theory

Taking into account as the above, uniform uni-axial compression loading along the x-axis, the following relation for the nonlocal strain theory is considered:
  • CPT:
    N x 0 = D ( α m 2 + β n 2 ) 2 [ 1 + e 0 2 θ 2 ( α m 2 + β n 2 ) ] α m 2
  • FSDT and HSDT:
    The fundamental relations for FSDT or HSDT can be derived in a similar manner as in the previous section. Since they are more complicated as for the nonlocal stress theory, we do not demonstrate them in the explicit form. Instead of, we show the effects in the form of graphs plotted in Figure 5. In general, the nonlocal strain theory leads to the growth of buckling loads in comparison with the nonlocal stress theory (Figure 4). The increase factor is equal to 1 + e 0 2 θ 2 ( α m 2 + β n 2 ) . Those effects are identical to the reported in the literature by Lim et al. [14] for beams.

4.2. Free Vibrations

Using the Rayleigh-Ritz method, it is possible to derive the natural frequencies easily for nanoplates. For the classical plate theory the frequencies are described by the simple equations having the following form:
  • the nonlocal stress theory:
    ω m n 2 ρ t = D ( α m 2 + β n 2 ) 2 1 + e 0 2 θ 2 ( α m 2 + β n 2 )
  • the nonlocal strain theory:
    ω m n 2 ρ t = D ( α m 2 + β n 2 ) 2 [ 1 + e 0 2 θ 2 ( α m 2 + β n 2 ) ]
In the dimensionless form of free vibrations, the influence of the constitutive equations is identical to those plotted in Figure 4 and Figure 5 for CPT.
Using the transverse shear deformation theory, the fundamental flexural frequencies are derived in a slightly different manner presented—see work [57] for the detailed description. In this case, the eigenproblem is transformed to (Table 2):
det K ω 2 M = 0   where   M = ρ t ρ t 3 / 12 ρ t 3 / 12
The fundamental frequency is the lowest value obtained from Equation (31).

5. Conclusions and Future Works

The present work is devoted to the characterization of the formulation of basic equations on the eigenfrequencies and buckling loads for the nanoplates considered as a 2-D continuum. The analysis allows us to formulate the following conclusions:
  • The description of the displacements with the use of the third order polynomial or exponential approximations leads to the identical results.
  • The first order shear deformation theory seems to be a good approximation; the use of the higher order shear deformation theory can be necessary as transverse shear effects are significantly enhanced by the low values of the Gxz/E and Gyz/E ratios.
  • In the strain or stress nonlocal problems, the evaluation of eigenvalues is reduced to the operation of division/multiplication by the identical factor; therefore two different physical formulations have different effects on the searched values understood in the sense of weakening or enhancement the buckling/vibration effects in comparison with the classical global approach.
The proposed Rayleigh-Ritz method can be easily applied to the analysis of different variants of boundary conditions—the proposed method of solution is used, e.g., by Shi et al. [77].
In view of the opposite effects on the nonlocal strain and stress formulations, it seems to be reasonable to investigate also Mindlin’s nonlocal theories.
In our opinion, there is an open question about the limits of the continuum theory application to the description of nanostructures. It is necessary to extend eigenproblems analysis of nanostructures/nanocomposites to the area of molecular dynamic problems—see, e.g., works by Muc and Barski [78] and Barski et al. [79].

Author Contributions

A.M. and M.C.—conceptualization and methodology, formal analysis, original draft preparation, review and editing, visualization.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Peddieson, J.; Buchanan, G.R.; McNitt, R.P. Application of nonlocal continuum models to nanotechnology. Int. J. Eng. Sci. 2003, 41, 305–312. [Google Scholar] [CrossRef]
  2. Lu, P.; Zhang, P.Q.; Lee, H.P.; Wang, C.M.; Reddy, J.N. Non-local elastic plate theories. Proc. R. Soc. A 2007, 463, 3225–3240. [Google Scholar] [CrossRef]
  3. Barretta, R.; Feo, L.; Luciano, R.; Marotti de Sciarra, F. Application of an enhanced version of the Eringen differential model to nanotechnology. Compos. B 2016, 96, 274–280. [Google Scholar] [CrossRef]
  4. Chen, W.J.; Li, X.P. Size-dependent free vibration analysis of composite laminated Timoshenko beam based on new modified couple stress theory. Arch. Appl. Mech. 2013, 83, 431–444. [Google Scholar] [CrossRef]
  5. Wang, Q. Wave propagation in carbon nanotubes via nonlocal continuum mechanics. J. Appl. Phys. 2005, 98, 124301. [Google Scholar] [CrossRef]
  6. Ma, H.M.; Gao, X.L.; Reddy, J.N. A non-classical Mindlin plate model based on a modified couple stress theory. Acta Mech. 2011, 220, 217–235. [Google Scholar] [CrossRef]
  7. Kim, J.; Żur, K.K.; Reddy, J.N. Bending, free vibration, and buckling of modified couples stress-based functionally graded porous micro-plates. Compos. Struct. 2019, 209, 879–888. [Google Scholar] [CrossRef]
  8. Barretta, R.; Feo, L.; Luciano, R.; Marotti de Sciarra, F. An Eringen-like model for Timoshenko nanobeams. Compos. Struct. 2016, 139, 104–110. [Google Scholar] [CrossRef]
  9. Hu, Y.G.; Liew, K.M.; Wang, Q.; He, X.Q.; Yakobson, B.I. Nonlocal shell model for elastic wave propagation in single- and double-walled carbon nanotubes. J. Mech. Phys. Sol. 2008, 56, 3475–3485. [Google Scholar] [CrossRef]
  10. Arefi, M.; Bidgoli, E.; Dimitri, R.; Bacciocchi, M.; Tornabene, F. Nonlocal bending analysis of curved nanobeams reinforced by graphene nanoplatelets. Compos. B 2019, 166, 1–12. [Google Scholar] [CrossRef]
  11. Arash, B.; Wang, Q. A review on the application of nonlocal elastic models in modeling of carbon nanotubes and graphenes. Comput. Mater. Sci. 2012, 51, 303–313. [Google Scholar] [CrossRef]
  12. Askes, H.; Aifantis, E.C. Gradient elasticity in statics and dynamics: An overview of formulations, length scale identification procedures, finite element implementations and new results. Int. J. Sol. Struct. 2011, 48, 1962–1990. [Google Scholar] [CrossRef]
  13. Akgoz, B.; Civalek, O. Bending analysis of FG microbeams resting on Winkler elastic foundation via strain gradient elasticity. Compos. Struct. 2015, 134, 294–301. [Google Scholar] [CrossRef]
  14. Lim, C.W.; Zhang, G.; Reddy, J.N. A higher-order nonlocal elasticity and strain gradient theory and its applications in wave propagation. J. Mech. Phys. Solids 2015, 78, 298–313. [Google Scholar] [CrossRef]
  15. Barretta, R.; Canadija, M.; Marotti de Sciarra, F. Modified nonlocal strain gradient elasticity for nano-rods and application to carbon nanotubes. Appl. Sci. 2019, 9, 514. [Google Scholar] [CrossRef]
  16. Nowacki, W. Couple stresses in thermoelasticity. Rozprawy Inżynierskie 1968, 16, 441–471. (In Polish) [Google Scholar]
  17. Voigt, W. Theoretische Studien über die Elastizitatsverhaltnisse der Kristalle. Abh Koniglichen Gesellschaft Wiss Gottingen 1887, 43. [Google Scholar]
  18. Cosserat, E.; Cosserat, F. Théorie des corps déformables; Herman et Fils: Paris, France, 1909. [Google Scholar]
  19. Altenbach, J.; Altenbach, H.; Eremeyev, V.A. On generalized Cosserat-type theories of plates and shells: A short review and bibliography. Arch. Appl. Mech. 2010, 80, 73–92. [Google Scholar] [CrossRef]
  20. Toupin, R.A. Theory of elasticity with couple stresses. Arch. Ration Mech. Anal. 1964, 17, 85–112. [Google Scholar] [CrossRef]
  21. Mindlin, R.D.; Tiersten, H.F. Effects of couple-stesses in linear elasticity. Arch. Ration. Mech. Anal. 1962, 11, 415–448. [Google Scholar] [CrossRef]
  22. Eringen, A.C. Nonlocal polar elastic continua. Int. J. Eng. Sci. 1972, 10, 1–16. [Google Scholar] [CrossRef]
  23. Eringen, A.C.; Edelen, D.G.B. On nonlocal elasticity. Int. J. Eng. Sci. 1972, 10, 233–248. [Google Scholar] [CrossRef]
  24. Eringen, A.C. On differential equations of nonlocal elasticity and solutions of screw dislocation and surface waves. J. Appl. Phys. 1983, 54, 4703–4710. [Google Scholar] [CrossRef]
  25. Aifantis, E. On the role of gradients in the localization of deformation and fracture. Int. J. Eng. Sci. 1992, 30, 1279–1299. [Google Scholar] [CrossRef]
  26. Aifantis, E. Gradient deformation models at nano, micro, and macro scales. J. Eng. Mater Technol. 1999, 121, 189–202. [Google Scholar] [CrossRef]
  27. Yang, F.; Chong, A.C.M.; Lam, D.C.C.; Tong, P. Couple stress based strain gradient theory for elasticity. Int. J. Solids Struct. 2002, 39, 2731–2743. [Google Scholar] [CrossRef]
  28. Mindlin, R.D. Micro-structure in linear elasticity. Arch. Rat. Mech. Anal. 1964, 16, 52–78. [Google Scholar] [CrossRef]
  29. Fernandez-Saez, J.; Zaera, R.; Loya, J.A.; Reddy, J.N. Bending of Euler-Bernoulli beams using Eringen’s integral formulation: A paradox resolved. Int. J. Eng. Sci. 2016, 99, 107–116. [Google Scholar] [CrossRef]
  30. Romano, G.; Barretta, R.; Diaco, M.; Marotti de Sciarra, F. Constitutive boundary conditions and paradoxes in nonlocal elastic nanobeams. Int. J. Mech. Sci. 2017, 121, 151–156. [Google Scholar] [CrossRef]
  31. Aifantis, E. On the microstructural origin of certain inelastic models. J. Eng. Mater. Technol. 1984, 106, 326–330. [Google Scholar] [CrossRef]
  32. Altan, S.; Aifantis, E. On the structure of the mode III crack-tip in gradient elasticity. Scripta Metall. Mater. 1992, 26, 319–324. [Google Scholar] [CrossRef]
  33. Aifantis, E. Strain gradient interpretation of size effects. Int. J. Fract. 1999, 95, 299–314. [Google Scholar] [CrossRef]
  34. Romano, G.; Barretta, R.; Diaco, M. Micromorphic continua: Non-redundant formulations. Contin. Mech. Thermodyn. 2016, 28, 1659–1670. [Google Scholar] [CrossRef]
  35. Barbagallo, G.; Madeo, A.; d’Agostino, M.V.; Abreu, R.; Ghiba, I.-D.; Neff, P. Transparent anisotropy for the relaxed micromorphic model: Macroscopic consistency conditions and long wave length asymptotics. Int. J. Sol. Struct. 2017, 120, 7–30. [Google Scholar] [CrossRef]
  36. Neff, P.; Madeo, A.; Barbagallo, G.; d’Agostino, M.V.; Abreu, R.; Ghiba, I.-D. Real wave propagation in the isotropic-relaxed micromorphic model. Proc. R. Soc. A 2017, 473, 20160790. [Google Scholar] [CrossRef]
  37. Sourki, R.; Hosseini, S.A. Coupling effects of nonlocal and modified couple stress theories incorporating surface energy on analytical transverse vibration of a weakened nanobeam. Eur. Phys. J. Plus 2017, 132, 184. [Google Scholar] [CrossRef]
  38. Reddy, J.N. Nonlocal theories for bending, buckling and vibration of beams. Int. J. Eng. Sci. 2007, 45, 288–307. [Google Scholar] [CrossRef]
  39. Pradhan, S. Buckling of single layer graphene sheet based on nonlocal elasticity and higher order shear deformation theory. Phys. Lett. A 2009, 373, 4182–4188. [Google Scholar] [CrossRef]
  40. Pradhan, S.; Murmu, T. Small scale effect on the buckling of single-layered graphene sheets under biaxial compression via nonlocal continuum mechanics. Comput. Mater. Sci. 2009, 47, 268–274. [Google Scholar] [CrossRef]
  41. Pradhan, S.; Phadikar, J. Small scale effect on vibration of embedded multilayered graphene sheets based on nonlocal continuum models. Phys. Lett. A 2009, 373, 1062–1069. [Google Scholar] [CrossRef]
  42. Sakhaee-Pour, A. Elastic buckling of single-layered graphene sheet. Comput. Mater. Sci. 2009, 45, 266–270. [Google Scholar] [CrossRef]
  43. Aghababaei, R.; Reddy, J.N. Nonlocal third-order shear deformation plate theory with application to bending and vibration of plates. J. Sound Vib. 2009, 326, 277–289. [Google Scholar] [CrossRef]
  44. Civalek, O.; Demir, C. Bending analysis of microtubules using nonlocal Euler–Bernoulli beam theory. Appl. Math. Model 2011, 35, 2053–2067. [Google Scholar] [CrossRef]
  45. Fazelzadeh, S.A.; Ghavanloo, E. Nonlocal anisotropic elastic shell model for vibrations of single-walled carbon nanotubes with arbitrary chirality. Compos. Struct. 2012, 94, 1016–1022. [Google Scholar] [CrossRef]
  46. Shen, Z.-B.; Tang, H.-L.; Li, D.-K.; Tang, G.-J. Vibration of single-layered graphene sheet-based nanomechanical sensor via nonlocal Kirchhoff plate theory. Comput. Mater. Sci. 2012, 61, 200–205. [Google Scholar] [CrossRef]
  47. Nami, M.R.; Janghorban, M. Resonance behavior of FG rectangular micro/nano plate based on nonlocal elasticity theory and strain gradient theory with one gradient constant. Compos. Struct. 2014, 111, 349–353. [Google Scholar] [CrossRef]
  48. Ansari, R.; Faghih Shojaei, M.; Mohammadi, V.; Gholami, R.; Darabi, M.A. Nonlinear vibrations of functionally graded Mindlin microplates based on the modified couple stress theory. Compos. Struct. 2014, 114, 124–134. [Google Scholar] [CrossRef]
  49. Barretta, R.; Marotti de Sciarra, F. Analogies between nonlocal and local Bernoulli–Euler nanobeams. Arch. Appl. Mech. 2015, 85, 89–99. [Google Scholar] [CrossRef]
  50. Sarrami-Foroushani, S.; Azhari, M. Nonlocal buckling and vibration analysis of thick rectangular nanoplates using finite strip method based on refined plate theory. Acta Mech. 2016, 227, 721–742. [Google Scholar] [CrossRef]
  51. Sladek, J.; Sladek, V.; Hrcek, S.; Pan, E. The nonlocal and gradient theories for a large deformation of piezoelectric nanoplates. Compos. Struct. 2017, 172, 119–129. [Google Scholar] [CrossRef]
  52. De Domenico, D.; Askes, H. Stress gradient, strain gradient and inertia gradient beam theories for the simulation of flexural wave dispersion in carbon nanotubes. Compos. B 2018, 153, 285–294. [Google Scholar] [CrossRef]
  53. Zenkour, A.M. A novel mixed nonlocal elasticity theory for thermoelastic vibration of nanoplates. Compos. Struct. 2018, 185, 821–833. [Google Scholar] [CrossRef]
  54. Zhu, J.; Lv, Z.; Liu, H. Thermo-electro-mechanical vibration analysis of nonlocal piezoelectric nanoplates involving material uncertainties. Compos. Struct. 2019, 208, 771–783. [Google Scholar] [CrossRef]
  55. Apuzzo, A.; Barretta, R.; Faghidian, S.A.; Luciano, R.; Marotti de Sciarra, F. Nonlocal strain gradient exact solutions for functionally graded inflected nano-beams. Compos. B 2019, 164, 667–674. [Google Scholar] [CrossRef]
  56. Ghavanloo, E.; Rafii-Tabar, H.; Fazelzadeh, S.A. Computational continuum mechanics of nanoscopic structures. Nonlocal elasticity approaches; Springer: New York, NY, USA, 2019. [Google Scholar]
  57. Muc, A. Design and identification methods of effective mechanical properties for carbon nanotubes. Mater. Design 2010, 31, 1671–1675. [Google Scholar] [CrossRef]
  58. Muc, A. Modelling of carbon nanotubes behaviour with the use of a thin shell theory. J. Theor. Appl. Mech. 2011, 49, 531–540. [Google Scholar]
  59. Muc, A. Natural frequencies of rectangular laminated plates-introduction to optimal design in aeroelastic problems. Aerospace 2018, 5, 95. [Google Scholar] [CrossRef]
  60. Muc, A.; Chwał, M. Vibration control of defects in carbon nanotubes. Solid Mech. Appl. 2011, 30, 239–246. [Google Scholar]
  61. Muc, A.; Banaś, A.; Chwał, M. Free vibrations of carbon nanotubes with defects. Mech. Mech. Eng. 2013, 17, 157–166. [Google Scholar]
  62. Chwał, M. Free vibrations analysis of carbon nanotubes. Adv. Mater. Res. 2014, 849, 94–99. [Google Scholar]
  63. Chwał, M. Nonlocal analysis of natural vibrations of carbon nanotubes. J. Mater. Eng. Perform. 2018, 27, 6087–6096. [Google Scholar]
  64. Muc, A.; Jamróz, M. Homogenization models for carbon nanotubes. Mech. Compos. Mater. 2004, 40, 101–106. [Google Scholar] [CrossRef]
  65. Chwał, M.; Muc, A. Transversely isotropic properties of carbon nanotube/polymer composites. Compos. B 2016, 88, 295–300. [Google Scholar]
  66. Chwał, M. Deformations and tensile fracture of carbon nanotubes based on the numerical homogenization. Acta Phys. Pol. A 2017, 131, 440–442. [Google Scholar]
  67. Reddy, J.N. A simple higher-order theory for laminated composite plates. J. Appl. Mech. 1984, 51, 745–752. [Google Scholar] [CrossRef]
  68. Hosseini-Hashemi, Sh.; Fadaee, M.; Es’haghi, M. A novel approach for in-plane/out-of-plane frequency analysis of functionally graded circular/annular plates. Int. J. Mech. Sci. 2010, 52, 1025–1035. [Google Scholar] [CrossRef]
  69. Muc, A.; Kędziora, P.; Stawiarski, A. Buckling enhancement of laminated composite structures partially covered by piezoelectric actuators. Eur. J. Mech. A 2019, 73, 112–125. [Google Scholar] [CrossRef]
  70. Muc, A.; Chwał, M.; Barski, M. Remarks on experimental and theoretical investigations of buckling loads for laminated plated and shell structures. Compos. Struct. 2018, 203, 861–874. [Google Scholar] [CrossRef]
  71. Shi, J.-X.; Natsuki, T.; Lei, X.-W.; Ni, Q.-Q. Equivalent Young’s modulus and thickness of graphene sheets for the continuum mechanical models. Appl. Phys. Lett. 2014, 104, 223101. [Google Scholar] [CrossRef]
  72. Shi, J.-X.; Ni, Q.-Q.; Lei, X.-W.; Natsuki, T. Study on wave propagation characteristics of double-layer graphene sheets via nonlocal Mindlin-Reissner plate theory. Int. J. Mech. Sci. 2014, 84, 25–30. [Google Scholar] [CrossRef]
  73. Shi, J.-X.; Ni, Q.-Q.; Lei, X.-W.; Natsuki, T. Nonlocal vibration analysis of nanomechanical systems resonators using circular double-layer graphene sheets. Appl. Phys. A 2014, 115, 213–219. [Google Scholar] [CrossRef]
  74. Muc, A. Choice of design variables in the stacking sequence optimization for laminated structures. Mech. Compos. Mater. 2016, 52, 211–224. [Google Scholar] [CrossRef]
  75. Muc, A.; Chwał, M. Analytical discrete stacking sequence optimization of rectangular composite plates subjected to buckling and FPF constraints. J. Theor. Appl. Mech. 2016, 54, 423–436. [Google Scholar] [CrossRef]
  76. Muc, A. Transverse shear effects in stability problems of laminated shallow shells. Compos. Struct. 1989, 12, 171–180. [Google Scholar] [CrossRef]
  77. Shi, J.W.; Nakatani, A.; Kitagawa, H. Vibration analysis of fully clamped arbitrarily laminated plate. Compos. Struct. 2004, 63, 115–122. [Google Scholar] [CrossRef]
  78. Muc, A.; Barski, M. Design of particulate-reinforced composite materials. Materials 2018, 11, 234. [Google Scholar] [CrossRef]
  79. Barski, M.; Chwał, M.; Kędziora, P. Molecular dynamics in simulation of magneto-rheological fluids behavior. Key Eng. Mater. 2013, 542, 11–27. [Google Scholar] [CrossRef]
Figure 1. An analyzed continuum nanoplate: (a) Geometry and buckling loading, (b) Example of deformed shape.
Figure 1. An analyzed continuum nanoplate: (a) Geometry and buckling loading, (b) Example of deformed shape.
Applsci 09 01409 g001
Figure 2. Variations of the function F(z) through the plate thickness t.
Figure 2. Variations of the function F(z) through the plate thickness t.
Applsci 09 01409 g002
Figure 3. The effects of the accuracy of 2-D kinematical formulations on the buckling loads—a square nanoplate, m = n = 1.
Figure 3. The effects of the accuracy of 2-D kinematical formulations on the buckling loads—a square nanoplate, m = n = 1.
Applsci 09 01409 g003
Figure 4. The influence of the constitutive relations (nonlocal stress theory) on the buckling loads.
Figure 4. The influence of the constitutive relations (nonlocal stress theory) on the buckling loads.
Applsci 09 01409 g004
Figure 5. The influence of the constitutive relations (nonlocal strain theory) on the buckling loads.
Figure 5. The influence of the constitutive relations (nonlocal strain theory) on the buckling loads.
Applsci 09 01409 g005
Table 1. Function F(z) for variants of plate theories.
Table 1. Function F(z) for variants of plate theories.
Classical Plate Theory (CPT)First Order Shear
Deformation Theory
(FSDT)
Higher Order Shear
Deformation Theory
(HSDT)–Reddy [67]
Higher Order Shear
Deformation Theory
(HSDT)–Hosseini-Hashemi et al. [68]
F(z)0z; the second term in Equation (10) is eliminated z 1 4 3 z 2 t 2 z exp 2 z 2 t 2
Table 2. The governing equations of motion for nanoplates – the nonlocal stress formulation.
Table 2. The governing equations of motion for nanoplates – the nonlocal stress formulation.
CPTFSDT/HSDT
Equations of equilibrium N x x x + N x y y = ρ t u ¨
N x y x + N y y y = ρ t v ¨
2 M x x x 2 + 2 2 M x y x y + 2 M y y y 2 = ( ρ t w ¨ + N x 0 2 w x 2 + N y 0 2 w y 2 )
N x x x + N x y y = ρ t u ¨
N x y x + N y y y = ρ t v ¨
Q x 3 x + Q y 3 y = ( ρ t w ¨ + N x 0 2 w x 2 + N y 0 2 w y 2 )
M x x x + M x y y Q x 3 = ρ t 3 ψ ¨ 1 / 12
M x y x + M y y y Q y 3 = ρ t 3 ψ ¨ 2 / 12
Kinematical variablesu, v, wu, v, w, ψ1, ψ2
Stress resultants N x x , N x y , N y y , M x x , M x y , M y y N x x , N x y , N y y , M x x , M x y , M y y , Q x 3 , Q y 3

Share and Cite

MDPI and ACS Style

Chwał, M.; Muc, A. Buckling and Free Vibrations of Nanoplates—Comparison of Nonlocal Strain and Stress Approaches. Appl. Sci. 2019, 9, 1409. https://doi.org/10.3390/app9071409

AMA Style

Chwał M, Muc A. Buckling and Free Vibrations of Nanoplates—Comparison of Nonlocal Strain and Stress Approaches. Applied Sciences. 2019; 9(7):1409. https://doi.org/10.3390/app9071409

Chicago/Turabian Style

Chwał, Małgorzata, and Aleksander Muc. 2019. "Buckling and Free Vibrations of Nanoplates—Comparison of Nonlocal Strain and Stress Approaches" Applied Sciences 9, no. 7: 1409. https://doi.org/10.3390/app9071409

APA Style

Chwał, M., & Muc, A. (2019). Buckling and Free Vibrations of Nanoplates—Comparison of Nonlocal Strain and Stress Approaches. Applied Sciences, 9(7), 1409. https://doi.org/10.3390/app9071409

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop