Next Article in Journal
Cinnamate-Intercalated Layered Yttrium Hydroxide: UV Light-Responsive Switchable Material
Next Article in Special Issue
Design, Modeling, and Experimental Validation of an Active Microcatheter Driven by Shape Memory Effects
Previous Article in Journal
Low Power Compact 3D-Constructed AlScN Piezoelectric MEMS Mirrors for Various Scanning Strategies
Previous Article in Special Issue
Design Guidelines for Thin Diaphragm-Based Microsystems through Comprehensive Numerical and Analytical Studies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nonlinear Thermal/Mechanical Buckling of Orthotropic Annular/Circular Nanoplate with the Nonlocal Strain Gradient Model

by
Mostafa Sadeghian
*,
Arvydas Palevicius
and
Giedrius Janusas
Faculty of Mechanical Engineering and Design, Kaunas University of Technology, Studentu 56, 51424 Kaunas, Lithuania
*
Author to whom correspondence should be addressed.
Micromachines 2023, 14(9), 1790; https://doi.org/10.3390/mi14091790
Submission received: 14 July 2023 / Revised: 15 September 2023 / Accepted: 16 September 2023 / Published: 19 September 2023
(This article belongs to the Special Issue N/MEMS Intelligent Structures: Design, Manufacturing, and Control)

Abstract

:
This article presents the nonlinear investigation of the thermal and mechanical buckling of orthotropic annular/circular single-layer/bilayer nanoplate with the Pasternak and Winkler elastic foundations based on the nonlocal strain gradient theory. The stability equations of the graphene plate are derived using higher-order shear deformation theory (HSDT) and first-order shear deformation theory (FSDT) considering nonlinear von Karman strains. Furthermore, this paper analyses the nonlinear thermal and mechanical buckling of the orthotropic bilayer annular/circular nanoplate. HSDT provides an appropriate distribution for shear stress in the thickness direction, removes the limitation of the FSDT, and provides proper precision without using a shear correction coefficient. To solve the stability equations, the differential quadratic method (DQM) is employed. Additionally, for validation, the results are checked with available papers. The effects of strain gradient coefficient, nonlocal parameter, boundary conditions, elastic foundations, and geometric dimensions are studied on the results of the nondimensional buckling loads. Finally, an equation is proposed in which the thermal buckling results can be obtained from mechanical results (or vice versa).

1. Introduction

Nanoplates including graphene sheets can be named as one category of eminent materials for the next generation of nanodevices. Graphene sheets can be utilized in the manufacture of many nanodevices like sensors and memory devices [1] and there are other applications such as nanosheet resonators [2], mass sensors [3], and gas sensors [4]. Because performing highly accurate experiments at the nanoscale is very demanding and high-cost, some techniques such as continuum-based modelling and molecular dynamics (MD) have been employed recently that have the potential to consider the atomic size and length effects. MD modelling of nanoscale structures needs time-consuming procedures and complicated computations. Continuum modelling of nanostructures is computationally less expensive than MD modelling; therefore, continuum theories have been formulated and employed for the investigation of the mechanical properties of nanoscale structures. In other words, continuum modelling of nanostructures can assist in understanding the results of experimental measurements (or MD modelling) more properly, because it can decrease the volume and computational time. In order to make computations easier, continuum mechanics can be employed as an appropriate option for theoretical investigations. The atomic length scales and interatomic force coefficients can be incorporated into the constitutive equations. In the classical continuum model, the influences of size in nanostructures are not regraded [5,6].
Experimental investigations have proved that there are size-dependent effects on the mechanical properties at the nanostructures [7]. Therefore, modified continuum theories are required to consider small-scale influences and maintain the ease and computational efficiency of continuum theories. Consequently, new theories have been suggested to take into account nanoscale effects in nanostructures. Various nonclassical continuous models have been presented considering the material length scale coefficient, and multiple modifications to classical elasticity theories including Eringen’s nonlocal elasticity theory [8], the strain gradient theory [9], and the nonlocal strain gradient theory [10] have been proposed.
Narendar [11] perused the buckling of rectangular nanoplates on the refined basis of the two-variable plate theory considering nonlocal scale through the Navier technique. He observed that by increasing the size of the square nanoplate, the buckling load ratio decreases significantly. Jamalpoor and his colleagues [12] studied the vibration and buckling of nanoplates exposed to electric as well as magnetic potentials based on the nonlocal plate theory. They noticed that the buckling load increased by increasing the external magnetic potential. Li and Hu [13] used nonlocal continuum theory to study the buckling, vibration, and bending of composite nanobeams. They concluded that the critical buckling load can be increased by reducing the nonlocal coefficient when the nonlocal coefficient is larger than the material characteristic coefficient. Shaban and Alibeigloo [14] investigated the bending and vibration analysis of carbon nanotubes on elastic foundations using the nonlocal three-dimensional theory of elasticity. They observed that the static and vibrational behaviors of carbon nanotubes are significantly affected by elastic foundations. Rokni et al. [15] studied the free vibration of circular/annular plates considering variable thickness with various boundary conditions using three-dimensional elasticity. Shaban and Mazaheri [16] analyzed the electrostatic behavior of microsandwich panels with the functionally graded core using the nonlocal three-dimensional theory of elasticity. They observed that increasing the nonlocal coefficient led to an increase in the displacements. The authors of Ref. [17] perused the free vibration of functionally graded plates on the elastic foundation with elastically restrained edges using FSDT. In their paper, they studied the influences of coefficients including thickness-to-radius, material distribution, various combinations of constraints at edges on the frequency, foundation stiffness coefficients, mode shape, and modal stress.
It is noticed that in the former analyses, the small-scale effect in sort of the stress nonlocality led to a softening–stiffness effect, while the strain gradient size dependency resulted in a hardening–stiffness effect. Consequently, various nonclassical continuum theories of elasticity have been presented to prognosticate various size dependencies in the mechanical features of nanostructures. The theory of nonlocal strain gradients is one of the prominent models for analyzing size-dependent mechanical features. Contrary to classical elasticity theory, the nonlocal strain gradient theory can consider the influences of stiffness–hardening and stiffness–softening by assuming the strain and stress gradient. Therefore, Lim and his colleagues [10] studied a new size-dependent elasticity model named nonlocal strain gradient theory considered both softening and stiffening effects to illustrate the size dependency more accurately.
Some analyses have been carried out with respect to the nonlocal strain gradient theory [18]. For example, the authors of Ref. [19] investigated the buckling of the nanocrystal shell with the axial loads assuming thermal, magnetic, and electrical conditions with the nonlocal strain gradient theory. They deduced that the nonlocal parameter affects the critical buckling load more notably for nonlocal parameters ranging from 2 to 20 nm. Tanzadeh and Amoushahi [20] perused the buckling of rectangular nanoplates subjected to various uniaxial in-plane loads using the nonlocal strain gradient theory with the higher-order finite strip technique. They concluded that increasing the strain gradient coefficient resulted in an increase in buckling load. Cuong-Le et al. [21] studied the bending, buckling, and vibration of sandwich nanoplates based on the theory of nonlocal strain gradients. They noticed that the nonlocal and strain gradient coefficients result in stiffness reduction and stiffness hardening cases. Therefore, both of these coefficients play a significant effect in the buckling of sandwich nanoplate. Wang and his colleagues [22] perused the buckling of functionally graded nanotubes based on the nonlocal strain gradient theory with high-order theory via the generalized differential quadrature technique. Wang et al. [23] analyzed the buckling of porous functionally graded porous nanobeams considering hybrid effects using the nonlocal strain gradient principle with the aid of the generalized differential quadrature technique. Esen and Özmen [24] studied thermal vibration and buckling functionally graded porous rectangular nanoplates considering magnet electroelastic effects based on the nonlocal strain gradient theory and FDST. They revealed that the temperature increase influences the buckling characteristic of the nanoplate depending on the proportions of the constituents in the nanoplate. Tang and Qing [25] examined the buckling and vibration of functionally graded beams based on the theory of nonlocal strain gradient and using the Laplace transform method. They scrutinized that the effects of length-scale coefficients on the vibration frequencies and buckling loads increase with the rise in vibration and buckling order. Al-Furjan and his colleagues [26] conducted research on the dynamic buckling of carbon nanocones with the nonlocal strain gradient theory and FSDT. They demonstrated that by increasing the strain gradient coefficient, the dynamic instability region occurred at high frequencies. Magorzata Chwa and Muc [27] utilized the nonlocal strain theory and considered higher-order shear deformation theories to analyze buckling and free vibrations of rectangular nanoplates via the Rayleigh–Ritz method. Using nonlocal strain gradient elasticity theory and third-order shear deformation plate theory, Fan et al. [28] studied buckling and postbuckling porous functionally graded nanoplates. They noticed that by increasing the material characteristic gradient index, both nonlocal and strain gradient size influences are noticeable. Sadeghian et al. [29] studied the investigation of large deflection of the annular/circular nanoplate on the basis of nonlocal strain gradient theory. They concluded that both strain gradient and nonlocal coefficients have noticeable effects on decreasing or increasing the deflection of the nanoplate.
One of the important issues in evaluating the stability of various structures is the buckling analysis, which has drawn the attention of many researchers to perform research on it. The term buckling refers to the loss of stability, or in other words, when the structure changes from stable equilibrium to unstable equilibrium, the structure undergoes buckling. In fact, the behavior of the structure under load in which with a small increase in load a disproportionate increase in displacement occurs in the structure is called buckling, and the amount of force for which the structure buckles is called critical buckling load [30]. There are many papers that consider buckling analysis. For example, Zghal et al. [31] scrutinized the postbuckling of functionally graded nanotubes considering various mechanical loadings. In their article, they discussed the influences of the carbon-nanotubes (CNT) volume fractions, CNT distributions, gradient indexes, and boundary conditions. In another paper, Zghal and Dammak [32] studied the buckling of functionally graded structures subjected to various compression loads. They noticed that the evenly porous FGM structures have the minimum buckling load while the uneven ones possess the maximum one. Trabelsi et al. [33] studied thermal buckling of functionally graded shells. They demonstrated that the flexural rigidity of the FG cone exposed to thermal loads can be improved with the variation of the power-law parameter. Mehr et al. [34] studied functionally graded CNT-reinforced composite shells with thermal buckling. They noticed that increasing aspect ratio and CNT volume fraction cause higher critical buckling temperature. Van Do and Lee [35] investigated the thermal buckling of FGM plates based on the quasi-3D higher-order shear deformation theory. They noticed that the temperature distribution on the plate noticeably influences the thermal buckling responses of FGM plates.
DQM is one of the numerical methods that has been employed by many authors [36,37]. The initial purpose of this procedure is to employ Lagrange interpolation polynomials to field constants and to solve the equations at discrete grid nodes. Higher accuracies are gained by considering more grid nodes. With the aid of DQM, Han et al. [38] studied one-electrode microresonators on the basis of a generalized 1DOF principle. Duryodhana et al. [37] employed DQM to study the buckling (and free vibrations) of foamed composites. By means of DQM, Ren et al. [39] investigated the thermal and mechanical buckling behavior of heated rectangular plates considering the characteristics of the temperature-dependent material.
FSDT accounts for the shear deformation influences by the way of linear variation of in-plane displacements through the thickness. Because the FSDT violates the conditions of zero transverse shear stresses on the top and bottom surfaces, a shear correction factor (which can depend on many parameters) is needed to compensate for the error caused by a constant shear strain assumption through the thickness. The HSDTs account for the shear deformation influences and satisfy the zero transverse shear stresses on the top and bottom surfaces; therefore, a shear correction factor is not needed. HSDTs are modelled based on the assumption of higher-order variations of in-plane displacements or both in-plane and transverse displacements through the thickness. Some of HSDTs are computational costs because with each additional power of the thickness coordinate, an additional unknown can be added to the model [40]. Also, there is quite the difference between FSDT and HSDT behavior for thick plates, but for thin plates both of the theories predict approximately similar result [41].
Until now, the nonlocal strain gradient theory with HSDT and FSDT has not been carried out for studying the nonlinear thermal as well as mechanical buckling of the orthotropic annular/circular nanoplate via DQM. Furthermore, the buckling analysis of a bilayer annular/circular nanoplate with elastic foundations is also studied. The effects of the nonlocal parameter, strain gradient coefficient, boundary conditions, geometric dimensions, and elastic foundations on the results of the nondimensional buckling loads are investigated. From the results, it can be concluded that the effect of the nonlocal parameter on the buckling of the circular/annular nanoplate is more significant than the types of boundary conditions or elastic foundations. Also, it is noticed that a bilayer nanoplate cannot be assumed as an equivalent single-layer nanoplate (with the same thickness as the bilayer nanoplate) to gain the results of the buckling of the nanoplate. The results of the current article can be a benefit for the design as well as improvement of nanostructured devices, including circular gate transistors, microswitches, and so on.

2. Theory and Formulation

Figure 1a,b demonstrate an annular graphene plate with inner radius r i , outer radius r o , and constant thickness h on the Winkler–Pasternak elastic foundation (in which k p and k w define the Pasternak and Winkler elastic foundations) and under uniform extended buckling loads.
Taking into account the HSDT, the displacement field can be written in the r, θ, and z axes clarified by U, V, and W, respectively.
U ( r , θ , z ) = u 0 ( r ) z d w 0 ( r ) d r + g ( z ) ϕ ( r ) V ( r , θ , z ) = 0 W ( r , θ , z ) = w 0 ( r )
where u 0 and w 0 define the displacements of the midplane on the r and z axes, respectively. Furthermore, ϕ ( r ) is the rotation component around the θ axis. Also, the function g(z) can be mentioned as:
g ( z ) = f ( z ) + z y .
f (z) and y* can be assumed as various functions utilized in different papers which are listed in Table 1 (for instance, the Ambartsumian [42] model can be considered as 1 6 z 3 f ( z ) + h 2 8 y z ) .
The nonlinear strains, considering von Karman’s presumptions, can be written as Equations (3)–(7):
ε r = d U d r + 1 2 ( d W d r ) 2 = d u 0 d r z d 2 w 0 d r 2 + g z d ϕ d r + 1 2 d w 0 d r 2   ,
ε θ = U r = 1 r u 0 z d w 0 d r + g z ϕ ,
γ r θ = 0 ,
γ r z = d W d r + d U d z = ϕ d g z d z ,
γ θ z = 0 .
It is noted since the symmetric assumption is considered, Equations (5) and (7) are equal to zero.
The nonlocal form (NL) of the force and momentum resultants can be clarified as:
N r , N θ , Q r N L = h 2 h 2 σ r , σ θ , σ r z N L d z   ,
M r , M θ N L = h 2 h 2 σ r , σ θ N L z d z   ,
R r , R θ N L = h 2 h 2 σ r , σ θ N L f z d z   ,
R r z N L = h 2 h 2 σ r z N L f z d z   .
The potential energy of the system is the sum of the strain energy created due to the work of internal forces as well as the potential energy created by external forces.
Π = U + Ω
Π, U, and Ω are described as the potential energy of the entire system, the strain energy of the structure, and the potential energy of external forces, respectively. According to the concept of minimum potential energy, for a structure in equilibrium, the variation in the potential energy is equal to zero.
δ Π = δ U + δ Ω = 0
By zeroing δ Π , the coefficients of δ u 0 , δ w 0 , and δ ϕ must be equal to zero. The Euler-Lagrange equations (which are in non-local form, so they are marked using superscript NL) are obtained as
δ u 0 : N r N L + r d N r N L d r N θ N L = 0 ,
δ ϕ : y r d M r N L d r + M r N L M θ N L r Q r N L + R r N L + r d R r N L d r R θ N L r R r z N L = 0 ,
δ w 0 : r d 2 M r N L d r 2 + 2 d M r N L d r d M θ N L d r + d w 0 d r N r N L + d N r N L d r r + r N r N L d 2 w 0 d r 2 + k w w 0 + k p 2 w 0 r = 0 .
Also, Equation (16) can be rewritten as
δ w 0 : r d 2 M r N L d r 2 + 2 d M r N L d r d M θ N L d r + d w 0 d r N θ N L + r N r N L d 2 w 0 d r 2 + k w w 0 + k p 2 w 0 r = 0 .
The nonlocal strain gradient theory (assumed as a combination of both the strain gradient theory and the nonlocal stress field) was suggested by Lim et al. [10]:
( 1 μ 2 2 ) σ i j = C i j k l ( 1 l 2 2 ) ε k l , 2 = d 2 d r 2 + 1 r d d r
in which C i j k l , l , and μ define elastic, strain gradient (internal material length scale), and nonlocal parameters, respectively. In addition, the constitutive stress–strain equation in the very small scale is defined as [49]
( 1 μ 2 2 ) σ r σ θ σ r z = ( 1 l 2 2 ) Q 11 Q 12 0 Q 12 Q 22 0 0 0 G 13 ε r ε θ γ r z ,   Q 11 = E 1 1 ν 12 ν 21   ,   Q 22 = E 2 1 ν 12 ν 21 Q 12 = ν 12 E 2 1 ν 12 ν 21   ,   G 13 = E 1 2 ( 1 + v 12 ) .
It should be mentioned in Equation (19), E 1 and E 2 indicate the Young modulus along the 1 and 2 directions. Additionally, v 12 and v 21 are Poisson ratios in those directions. Also, G 13 signifies the shear modulus.
The nonlocal form can be mentioned as
1 μ 2 N r , N θ , Q r N L = h 2 h 2 1 μ 2 σ r , σ θ , σ r z N L d z .
The force and moment resultants in the local form can be written as
N r , N θ , Q r L = h 2 h 2 σ r , σ θ , σ r z L d z   ,
M r , M θ L = h 2 h 2 σ r , σ θ L z d z   ,
R r , R θ L = h 2 h 2 σ r , σ θ L f z d z   ,
R r z L = h 2 h 2 σ r z L f z d z   .
Moreover, the resultants in terms of displacements are gained as follows:
N r L = ( 1 l 2 2 ) { 1 1 ν 12 ν 21 E 1 h d u 0 d r + 1 2 d w 0 d r 2 + ν 12 E 2 h 1 r u 0 + E 1 d ϕ d r + ν 12 E 2 1 r ϕ h 2 h 2 f z d z } ,
N θ L = ( 1 l 2 2 ) { 1 1 ν 12 ν 21 ν 12 E 2 h d u 0 d r + 1 2 d w 0 d r 2 + E 2 h 1 r u 0 + ν 12 E 2 d ϕ d r + E 2 1 r ϕ h 2 h 2 f z d z } ,
M r L = ( 1 l 2 2 ) { 1 1 ν 12 ν 21 E 1 h 3 12 d 2 w 0 d r 2 + y d ϕ d r + ν 12 E 2 h 3 12 1 r d w 0 d r + y 1 r ϕ + E 1 d ϕ d r + ν 12 E 2 1 r ϕ h 2 h 2 z f z d z } ,
M θ L = ( 1 l 2 2 ) { E 2 1 ν 12 ν 21 ν 12 h 3 12 d 2 w 0 d r 2 + y d ϕ d r + h 3 12 1 r d w 0 d r + y 1 r ϕ + ν 12 d ϕ d r + 1 r ϕ h 2 h 2 z f z d z } ,
R r L = ( 1 l 2 2 ) { 1 1 ν 12 ν 21 E 1 d u 0 d r + 1 2 d w 0 d r 2 + ν 12 E 2 1 r u 0 h 2 h 2 f ( z ) d z + ν 12 E 2 1 r d w 0 d r + y 1 r ϕ + E 1 d 2 w 0 d r 2 + y d ϕ d r h 2 h 2 z f ( z ) d z + E 1 d ϕ d r + ν 12 E 2 1 r ϕ h 2 h 2 f ( z ) 2 d z } ,
R θ L = ( 1 l 2 2 ) { E 2 1 ν 12 ν 21 ν 12 d u 0 d r + 1 2 d w 0 d r 2 + 1 r u 0 h 2 h 2 f ( z ) d z + ν 12 d 2 w 0 d r 2 + y d ϕ d r + 1 r d w 0 d r + y 1 r ϕ h 2 h 2 z f ( z ) d z + ν 12 d ϕ d r + 1 r ϕ h 2 h 2 f ( z ) 2 d z } ,
Q r L = ( 1 l 2 2 ) { G 13 ϕ y h + G 13 ϕ h 2 h 2 f ( z ) + y d z } ,
R r z L = ( 1 l 2 2 ) { G 13 ϕ h 2 h 2 f z 2 d z + G 13 y ϕ h 2 h 2 f z d z } .
The equilibrium equations of a monolayer axisymmetric circular/annular nanoplate considering Pasternak and Winkler elastic foundations can be locally written as
δ u 0 : N r L + r d N r L d r N θ L = 0 ,
δ ϕ : y r d M r L d r + M r L M θ L r Q r L + R r L + r d R r L d r R θ L r R r z L = 0 ,
δ w 0 : r d 2 M r L d r 2 + 2 d M r L d r d M θ L d r + 1 μ 2 k w w 0 + k p 2 w 0 r + N θ L d w 0 d r + r N r L d 2 w 0 d r 2 + μ r 2 N r L d 2 w 0 d r 2 + 2 N θ L 1 r d w 0 d r = 0 ,
To derive the stability equations and determine the critical buckling force, the adjacent equilibrium method is used. Therefore, displacements are replaced by the following equations:
u = u 0 + u 1 ϕ = ϕ 0 + ϕ 1 w = w 0 + w 1 .
Therefore, for the resultants of force and moment
N r L = N 0 r L + N 1 r L ,   N θ L = N 0 θ L + N 1 θ L   M r L = M 0 r L + M 1 r L   ,   M θ L = M 0 θ L + M 1 θ L R r L = R 0 r L + R 1 r L   ,   R θ L = R 0 θ L + R 1 θ L   Q r L = Q 0 r L + Q 1 r L ,   R r z L = R 0 r z L + R 1 r z L .
By substituting Equation (37) into equilibrium Equations (33)–(35), Equations (38)–(40) are obtained:
N 0 r L N 0 θ L + r d N 0 r L d r + N 1 r L N 1 θ L + r d N 1 r L d r = 0 ,
y r d M 0 r L d r + M 0 r L M 0 θ L r Q 0 r L + R 0 r L R 0 θ L + r d R 0 r L d r r R 0 r z L + y r d M 1 r L d r + M 1 r L M 1 θ L r Q 1 r L + R 1 r L R 1 θ L + r d R 1 r L d r r R 1 r z L = 0 ,
2 d M 0 r L d r d M 0 θ L d r + r d 2 M 0 r L d r 2 + 1 μ 2 k w w 0 + k p 2 w 0 r + d N 1 r L d r r d w 0 d r + N 1 r L d w 0 d r + r N 1 r L d 2 w 0 d r 2 + N 0 r L d w 0 d r + d N 0 r L d r r d w 0 d r + r N 0 r L d 2 w 0 d r 2 + μ ( ( 2 N 1 r L ) d 2 w 0 d r 2 + ( 2 N 1 θ L ) ( 1 r d w 0 d r ) ) + μ ( ( 2 N 0 r L ) d 2 w 0 d r 2 + ( 2 N 0 θ L ) ( 1 r d w 0 d r ) ) + 2 d M 1 r L d r d M 1 θ L d r + r d 2 M 1 r L d r 2 + 1 μ 2 k w w 1 + k p 2 w 1 r + d N 1 r L d r r d w 1 d r + N 1 r L d w 1 d r + r N 1 r L d 2 w 1 d r 2 + N 0 r L d w 1 d r + d N 0 r L d r r d w 1 d r + r N 0 r L d 2 w 1 d r 2 + μ ( ( 2 N 1 r L ) d 2 w 1 d r 2 + ( 2 N 1 θ L ) ( 1 r d w 1 d r ) ) + μ ( ( 2 N 0 r L ) d 2 w 1 d r 2 + ( 2 N 0 θ L ) ( 1 r d w 1 d r ) ) = 0 .
In Equations (38)–(40), by omitting the terms that have only zero subscript and setting terms related to the preload state equal to zero and considering the uniform compressive load, we have
N 0 r L = N 0 θ L = N ¯ .
Therefore, the stability equations in terms of local stresses are written in the form of Equations (42)–(44):
N 1 r L N 1 θ L + r N 1 r L r = 0 ,
y M 1 r L M 1 θ L r Q 1 r L + r d M 1 r L d r + R 1 r L R 1 θ L + r d R 1 r L d r r R 1 r z L = 0 ,
2 d M 1 r L d r d M 1 θ L d r + r d 2 M 1 r L d r 2 + 1 μ 2 k w w 1 + k p 2 w 1 r + d N 1 r L d r r + N 1 r L d w 1 d r + r N 1 r L d 2 w 1 d r 2 N ¯ d w 1 d r + r d 2 w 1 d r 2 + μ 2 N 1 r L d 2 w 1 d r 2 + 2 N 1 θ L 1 r d w 1 d r = 0 .

2.1. Derivation of the Governing Equations for the Bilayer Annular/Circular Nanoplate

Numerous graphene layers can be employed to fix their weak buckling strengths. Therefore, graphene plates are set on top of each other by means of van der Waals bonds, generating graphene layers [42]. Van der Waals force is a general term used to explain the attraction of intermolecular forces among molecules (for more information please refer to Refs. [50,51]). Figure 2 shows a bilayer nanoplate (considering van der Waals forces between layers) on the Pasternak and Winkler elastic foundations.
Equilibrium equations considering the bilayer annular/circular nanoplate with the Pasternak–Winkler elastic foundation can be procured similar to the procedure used for the single-layer nanoplate. Furthermore, the annular/circular displacement fields of the axisymmetric bilayer nanoplate can be written as (i = 1 shows the top layer and i = 2 illustrates the bottom layer)
U i ( r , θ , z ) = u i ( r ) z d w i ( r ) d r + g z ϕ i ( r )   ,   i = 1 , 2 ,
V i ( r , θ , z ) = 0 ,   i = 1 , 2 ,
W i ( r , θ , z ) = w i ( r ) ,   i = 1 , 2 .
In addition, the strain equations are the same as those procured in the single-layer annular/circular nanoplate. It should be noted that considering the minimum potential energy method to achieve the equilibrium equations (as well as boundary conditions) of the top and bottom layers, the following equations are considered:
δ Ω 1 = r i r o 0 τ k o w 2 w 1 δ w 1 r d r d θ ,
δ Ω 2 = r i r o 0 τ k o w 2 w 1 k w w 2 + k p 2 w 2 δ w 2 r d r d θ ,
δ U = v σ 1 i j N L δ ε 1 i j d v + v σ 2 i j N L δ ε 2 i j d v i , j = r , θ ,
δ Π = δ U + δ Ω 1 + δ Ω 2 = 0 .
k o is the van der Waals stiffness constant between layers. Furthermore, the equilibrium equations in terms of local stresses for both layers can be mentioned as follows:
δ u i : N i r L + r d N i r L d r N i θ L = 0 , i = 1 , 2 ,
δ ϕ i : y r d M i r L d r + M i r L M i θ L r Q i r L + R i r L + r d R i r L d r R i θ L r R i r z L = 0 ,   i = 1 , 2 ,
δ w 1 : r d 2 M 1 r L d r 2 + 2 d M 1 r L d r d M 1 θ L d r + 1 μ 2 k 0 w 2 w 1 r + N 1 θ L d w 1 d r + r N 1 r L d 2 w 1 d r 2 N ¯ d w 1 d r + r d 2 w 1 d r 2 + μ 2 N 1 r L 2 w 1 r 2 + 2 N 1 θ L 1 r w 1 r = 0 ,
δ w 2 : r d 2 M 2 r L d r 2 + 2 d M 2 r L d r d M 2 θ L d r + 1 μ 2 k 0 w 2 w 1 k w w 2 + k p 2 w 2 r + N 2 θ L d w 2 d r + r N 2 r L d 2 w 2 d r 2 N ¯ d w 2 d r + r d 2 w 2 d r 2 + μ 2 N 2 r L 2 w 2 r 2 + 2 N 2 θ L 1 r w 2 r = 0 .
In addition, the stability equations can be presented for the thermal buckling of circular and annular plates. Moreover, the only difference between the equations of thermal buckling compared to mechanical buckling is in the strains. Consequently, the stress and moment resultants change. The nonlinear strains, taking into account von Karman’s assumptions, can be written as Equations (56)–(58):
ε r = d u d r z d 2 w d r 2 + g z d ϕ d r + 1 2 ( d w d r ) 2 α 11 Δ T ,
ε θ = u r z r d w d r + 1 r g z ϕ α 22 Δ T ,
γ r z = ϕ d g z d z .
The local stress resultants (considering thermal effects) according to displacements are as follows:
N r L = 1 1 ν 12 ν 21 E 1 h d u d r + 1 2 d w d r 2 + ν 12 E 2 h 1 r u + E 1 d ϕ d r + ν 12 E 2 1 r ϕ h 2 + h 2 f z d z h Δ T E 2 α 22 ν 12 + E 1 α 11 1 ν 12 ν 21 ,
N θ L = 1 1 ν 12 ν 21 ν 12 E 2 h + d u d r + 1 2 d w d r 2 + E 2 h 1 r u + ν 12 E 2 d ϕ d r + E 2 1 r ϕ h 2 + h 2 f z d z E 2 h Δ T α 11 ν 12 + E 1 α 22 1 ν 12 ν 21 .
In order to derive the stability equations from the equilibrium equations and determine the critical buckling temperature using the adjacent equilibrium technique, the equilibrium equations are obtained in the form of Equations (38)–(40). In those equations, by ignoring the terms that only have a zero subscript, zeroing the terms related to the preload state, and knowing that in thermal buckling
N 0 r L = N 0 θ L = N ¯ T = h Δ T E 2 α 22 ν 12 + E 1 α 11 1 ν 12 ν 21 .
Stability equations in terms of local stresses are written in the form of Equations (62)–(64):
N 1 r L N 1 θ L + r N 1 r L r = 0 ,
y M 1 r L M 1 θ L r Q 1 r L + r M 1 r L r + R 1 r L R 1 θ L + r R 1 r L r r R 1 r z L = 0 ,
2 M 1 r L r M 1 θ L r + r 2 M 1 r L r 2 + 1 μ 2 N 1 r L r r + N 1 r L w 1 r + r N 1 r L 2 w 1 r 2 N ¯ T w 1 r + r 2 w 1 r 2 + μ 2 N 1 r L 2 w 1 r 2 + 2 N 1 θ L 1 r w 1 r = 0 .

2.2. Equations of Thermal Stability for a Bilayer Annular/Circular Plate

The stability equations in terms of local stresses (considering the first layer (i = 1) and the second layer (i = 2)) can be written as follows:
δ u i : N i r L + r d N i r L d r N i θ L = 0 , i = 1 , 2 ,
δ ϕ i : y r d M i r L d r + M i r L M i θ L r Q i r L + R i r L + r d R i r L d r R i θ L r R i r z L = 0 ,   i = 1 , 2 ,
δ w 1 : r d 2 M 1 r L d r 2 + 2 d M 1 r L d r d M 1 θ L d r + 1 μ 2   k 0 w 2 w 1 r + N 1 θ L d w 1 d r + r N 1 r L d 2 w 1 d r 2 N ¯ T d w 1 d r + r d 2 w 1 d r 2 + μ 2 N 1 r L 2 w 1 r 2 + 2 N 1 θ L 1 r w 1 r = 0 ,
δ w 2 : r d 2 M 2 r L d r 2 + 2 d M 2 r L d r d M 2 θ L d r + 1 μ 2 k 0 w 2 w 1 k w w 2 + k p 2 w 2 r + N 2 θ L d w 2 d r + r N 2 r L d 2 w 2 d r 2 N ¯ T d w 2 d r + r d 2 w 2 d r 2 + μ 2 N 2 r L 2 w 2 r 2 + 2 N 2 θ L 1 r w 2 r = 0 .

2.3. Boundary Conditions

The boundary conditions for the annular/circular plate can be assumed as
Clamped (C):
u = w = ϕ = d w d r = 0 .
Simply supported (S):
u = w = M r = R r = 0 .
Free (F):
N r = M r = R r = Q r = 0   .

2.4. Nondimensional Assumptions

Because of the very small values that exist on the nanoscale, and for the ease of calculations, Equation (72) is introduced to convert the equations of the nanoplate nondimensional:
u = u 0 h ; w = w 0 r o ; ϕ = ϕ   ;   ψ = ψ ; N r = N r E 1 h ;   N θ = N θ E 1 h ;   Q r = Q r E 1 h   ; Q θ = Q θ E 1 h ; M r = M r E 1 h 2   ;   M θ = M θ E 1 h 2   ;   2 = d 2 d r 2 + 1 r d d r ; R r = R r E 1 h 2 ; R θ = R θ E 1 h 2 ; R r z = R r z E 1 h ;   r = r r o   ;   z = z h   ;   δ = h r o ; k w = k w r o E 1   ; k p = k p E 1 r o

3. Numerical Procedure

DQM is one of the most efficient and elegant techniques for the numerical solution of partial differential equations used by many authors (which are briefly explained in the Introduction). This method is obtained from the quadratic integration technique, which expresses the integral at one point in the direction of the domain and relies on all points along that direction. Moreover, the value of dependency can be obtained using weight constants:
a b f ( r ) d r = k = 1 n w k f k
where, w 1 , w 2 , , w n as well as f 1 , f 2 , , f n can be defined as weight constants and function values at discrete points, respectively. Belman et al. [52] (regarding quadratic integration) proposed that the derivative at one node of the function domain can be based on the function values at the entire nodes of the domain via weight constants.
d f d r r i = j = 1 N A i j f r j   ,   i = 1 , 2 , , N
in which A i j and N are the weight constant and the total number of nodes in the direction of r, respectively. For the first-order derivative, the weighting constants can be gained as
A i j 1 = P r i r i r j P r j ,
P r i = j = 1 N r i r j   ,   i j ,
A i i 1 = k = 1 N A i k 1 ,   i k .
For higher-order derivatives, the following equation can be defined as
d n f d r n r i = j = 1 N A i j ( n ) f ( r j ) ,   i = 1 , , N .
For the second- and higher-order derivatives, the weighting constants can be defined as follows:
A i j n = n A i j 1 A i i n 1 A i j n 1 r i r j   ,   i j ,
A i i n = j = 1 , i N A i j n   ,   i , j = 1 N .
In this study, the grid points are distributed according to the Chebyshev–Gauss–Lubato distribution as follows:
r i = r i + r o 2 cos i 1 N 1 π r o r i 2 ,   i = 1 N
where r i and r o are nodes at the start and end of the function.
It should be noted that for each node (considering a single-layer plate) there are three equations and four unknowns, whose unknowns are displacement components of the plate ( w 0 , φ , u 0 ) and critical buckling load in the case of mechanical/thermal axial symmetric analysis. In other words, there are 27 equations and 28 unknowns for 9 nodes. The desired unknown (mechanical buckling critical load or thermal buckling critical term) is obtained after applying the DQM. Therefore, by inserting Equations (25)–(32) (which are the resultants in terms of displacements) in the stability equations, the stability equations in terms of displacements are obtained. Then, by using Equation (74) (for first-order derivatives) or Equation (78) (for higher-order derivatives) and by using appropriate weighting constants (Equations (75)–(77) for the first-order derivative or Equations (79)–(80) for higher derivatives) the discretized form of equations will be obtained in terms of displacements. The matrix form of these equations can be written as follows:
C e f 3 N × 3 N u 0 1 φ 1 w 0 1 . . u 0 N φ N w 0 N = 0 .
C e f is the matrix form of coefficients of equations. For N nodes, the matrix form of equations includes 3*N algebraic equations and 3*N + 1 unknowns. Since solving these algebraic equations and unknowns are long and time-consuming, Maple software 2023 was used to solve it. So, in the computer program, the determinant of the matrix of coefficients should be equal to zero ( C e f = 0 ) to obtain the dimensionless critical load of mechanical buckling or the dimensionless term of thermal buckling. Then, an equation will be obtained that contains a Polynomial in terms of buckling load. By solving it, and in order to find the critical buckling load, the minimum value should be chosen.
Also, for a better understanding of applying DQM, an example of the discrete form of the equations can be seen in Appendix A.

4. Results and Discussion

In this part, various factors are studied to identify their effect on the nondimensional buckling load of the annular/circular nanoplate. For this purpose, HSDT, FSDT, and the nonlocal strain gradient theory are utilized with the aid of DQM. Moreover, to validate the procedure, the current results are compared with the results of references. Also, the following assumptions are considered (if not mentioned in the text):
k w = 1.13 GPa / nm   ;   k p = 1.13 Pa . m   , k 0 = 45 GPa / nm ,   r o = 5 nm , E 1 = 1765   ( GPa ) ,   E 2 = 1588   ( GPa ) , ν 12 = 0.3 ,   ν 21 = 0.27 ,   h = 0.34   ( nm ) .    
Figure 3 reveals the influence of the number of nodes employed in the DQM to achieve the results of this paper. As can be seen, after nine nodes the appropriate convergence is gained. So, the number of nine nodes is considered to obtain the numerical results.
To verify accuracy and validity, the present results (gained using both FSDT and HSDT) in the clamped boundary condition are compared with the reference, which shows good accuracy (Table 2). In addition, it can be concluded that by increasing the radius of the plate, the nondimensional buckling loads are decreased.
Table 3 compares the nondimensional buckling load gained by HSDT by considering various shape functions for a circular nanoplate at clamped as well as simply supported boundary conditions. Therefore, different functions are used to distribute the shear stress along the thickness, which can be summarized as:
g 1 z = h π sin π z h ,
g 2 z = 4 3 h 2 z 3 + z ,
g 3 z = h sinh z h z cosh 1 2 ,
g 4 z = z e 2 z h 2 ,
g 5 z = 5 3 h 2 z 3 + 5 4 z .
From Table 3, it can be deduced that using these shape functions results in obtaining almost the same values of dimensionless buckling loads for both clamped and simply supported boundary conditions.
Figure 4 illustrates the changes in the nondimensional buckling loads versus the nonlocal constant for the annular nanoplate at various boundary conditions. It can be seen that by increasing the nonlocal parameter the non-dimensional buckling load decreases. This may be persuaded by increasing the nonlocal parameter; the stiffness of the nanoplate is increased and, therefore, the nondimensional buckling of the plate is reduced. Also, by increasing the rigidity of the boundary condition, increasing the nonlocal parameter results in a significant decrease in the nondimensional buckling load. Moreover, as the nonlocal coefficient increases the effect of the boundary condition on the nondimensional load decreases.
From Figure 5 and Figure 6, it is distinguished that generally increasing the elastic foundation results in an increase in the nondimensional buckling load, which is due to the fact that assuming an elastic medium results in a stiffer structure. Also, it is noticed that in local analysis, the nondimensional buckling loads with the Winkler foundation are higher than values with the Pasternak foundation, but by increasing the nonlocal parameter the Pasternak foundation has more values of nondimensional buckling loads. Moreover, In the local analysis, by increasing the flexibility of the boundary condition, the nondimensional buckling load decreases, but it can be noticed in Figure 5 and Figure 6, by increasing the rigidity of the boundary condition and increasing the nonlocal coefficient, the nondimensional buckling load decreases. This can highlight that in nonlocal analyses, the influence of the nonlocal coefficient on the buckling analysis of nanoplate is more significant than the type of boundary conditions or elastic foundations.
Figure 7 and Figure 8 illustrate the nondimensional buckling load versus the nonlocal parameter obtained from HSDT and FSDT for two thickness-to-radius ratios including h/r = 0.05 and h/r = 0.06. As can be deduced by increasing the nonlocal parameter, the non-dimensional buckling load obtained from HSDT and FSDT differ significantly. Also, by comparing these two figures, it can be seen that in the higher thickness-to-radius ratios the difference between FSDT and HSDT increases more significantly.
Figure 9 shows the nondimensional buckling load versus the radius ratios (inner radius-to-outer radius ratios) for different boundary conditions and considering different values of strain gradients. It can be concluded that by increasing radius ratios the buckling loads decrease. Also, increasing the strain gradient results in a decrease in buckling loads which specifies that inclusion of the strain gradient coefficient makes the plate stiffer than that of the classical plate [54].

4.1. Thermal Analysis

In this section, to validate the research, the current results of thermal buckling of circular/annular plates are compared with the results of valid reference, and then the results of the relevant diagrams are presented. Figure 10 reveals the effect of the number of nodes employed in the DQM to achieve the results of thermal analysis. It is noted that T is the critical buckling temperature, which is considered as follows:
T = 10 3 α T Δ T .
Δ T is the difference between the critical buckling temperature and the reference temperature. The coefficient of thermal expansion is considered as
α T = 1.1 × 10 6 K 1   .
In order to validate the numerical results of the present article and the solution technique and compare them with the results of Ref. [55], the dimensionless thermal buckling parameter ( λ T ) is defined as follows:
λ T = 12 1 + ν α T   Δ T r o h 2 .
In this comparison, the dimensionless parameter of the thermal buckling of the isotropic graphene plate for the annular nanoplate assuming various boundary conditions, radius ratios, and thickness-to-radius ratios is given in Table 4. According to this table, it is clear that in smaller thickness-to-radius ratios, the results agree well with the results of Ref. [43], and as the h / r o and r i / r o ratios increase, a small amount of difference is observed. Moreover, the difference between the results is caused by the different types of theories.
To better observe the differences between the present study with the reference [55], the percentage errors are gathered in Table 5, which can be defined as follows:
E 1 ( % ) = E HSDT E Ref E Ref × 100 , E 2 ( % ) = E FSDT E Ref E Ref × 100
where E 1 ( % ) and E 2 ( % ) represent the HSDT and FSDT errors compared to the reference’s values, respectively. Also, E HSDT , E FSDT , and E Ref represent values of the dimensionless thermal buckling parameter ( λ T ) gained by the present study HSDT, FSDT, and the reference (obtained from Table 4), respectively.
Figure 11 shows the critical buckling temperature of the annular nanoplate versus the nonlocal parameter for different boundary conditions. As can be distinguished, the critical buckling temperature decreases with the rise of the nonlocal parameter (which has the same trend as the mechanical analysis).
Finally, by comparing the equations and checking the results for both thermal and mechanical analyses, it is found that with the following equation, thermal and mechanical buckling analyses can be related to each other:
T = N λ , λ = ( 1 ν 12 ν 21 ) α T ( 1 + E 2 E 1 ν 12 ) .
It is noted that in Equation (93), T is critical buckling temperature and N is mechanical load. In other words, the results of the thermal buckling can be gained using Equation (93) which can save the time of computations.

4.2. Bilayer Analysis

Figure 12, Figure 13, Figure 14 and Figure 15 demonstrate the variation in the critical buckling temperature in terms of nonlocal constants for the annular bilayer nanoplate with a thickness of h (0.34 nm) and the single-layer nanoplate considering a thickness of 2*h (0.68 nm) for different boundary conditions. In fact, this part studies the probability of replacing a single-layer nanoplate (considering double thickness) with a bilayer nanoplate. It can be concluded that a bilayer nanoplate cannot be considered a single-layer nanoplate with the same thickness (as the bilayer nanoplate) for obtaining the exact results of the buckling analysis. Furthermore, the results obtained from the equivalent single-layer nanoplate with double thickness overestimate the real results gained using a bilayer nanoplate. This can be justified because of the van der Waals interaction that exists between the layers of bilayer nanoplate, which is an approximately weak force. If the buckling loads are strong enough, then molecules break free of the van der Waals forces that hold them together. Thus, the results of a single-layer plate are not similar to those of a bilayer plate.

5. Conclusions

In this paper, the nonlinear mechanical and thermal buckling of the circular and annular (single layer/bilayer) nanoplate using the nonlocal strain gradient theory considering both HSDT and FSDT is carried out. DQM is employed to solve the governing equations of the annular/circular nanoplate. Also, the results were compared with other references that showed good agreement. Some of the results of this article can be summarized as follows:
  • In the higher thickness-to-radius ratios, the difference between FSDT and HSDT increases more significantly.
  • The influence of the nonlocal coefficient on the buckling analysis of circular/annular nanoplate is more significant than the types of boundary conditions or elastic foundations.
  • A bilayer nanoplate cannot be considered an equivalent single-layer nanoplate (with the same thickness as the bilayer nanoplate) to obtain the results of the buckling analysis of nanoplate. Furthermore, the results obtained from the equivalent single-layer nanoplate with double thickness overestimate the real results gained using a bilayer nanoplate.

Author Contributions

M.S.: Conceptualization, methodology, software, validation, formal analysis, investigation, and writing: preparation of original draughts. A.P.: Supervision, conceptualization, methodology, and formal analysis. G.J.: validation and investigation. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data sharing is unavailable due to privacy or ethical restrictions.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

As an example, in this part, the discrete forms of the dimensionless equations of stability for circular nanoplate on the Winkler elastic foundation considering FSDT using DQM are considered. The obtained FSDT stability equations in terms of displacements are as follows:
1 2 ( r ( 1 ν 12 ν 21 ) ) ( 2 r ( d u 0 ( r ) d r δ + r ( d w ( r ) d r ) 2 ν 12 α r ( d w ( r ) d r ) 2 2 α u ( r ) δ + 2 r 2 ( d 2 u 0 ( r ) d r 2 ) δ + 2 r 2 ( d w ( r ) d r ) ( d 2 w ( r ) d r 2 ) ) = 0 ,
1 12 1 r δ ( 1 ν 12 ν 21 ) ( δ 2 ( d 2 φ ( r ) d r 2 ) r 2 + δ 2 ( d φ ( r ) d r ) r α δ 2 φ ( r ) + 12 r 2 K s β φ ( r ) ν 12 ν 21 12 r 2 K s β φ ( r ) + 12 r 2 K s β ( d w ( r ) d r ) ν 12 ν 21 12 r 2 K s β ( d w ( r ) d r ) = 0 ,
1 r 3 δ ( K s β r 3 δ ( d φ ( r ) d r ) K s β r 3 δ ( d 2 w ( r ) d r 2 ) K s β r 2 δ φ ( r ) K s β r 2 δ ( d w ( r ) d r ) + N r 3 δ ( d 2 w ( r ) d r 2 ) + N r 2 δ ( d w ( r ) d r ) + K w w ( r ) r 3 μ N δ ( d 4 w ( r ) d r 4 ) r 3 μ N δ ( d w ( r ) d r ) + μ N δ ( d 2 w ( r ) d r 2 ) 2 μ N δ ( d 3 w ( r ) d r 3 ) r 2 μ K w ( d 2 w ( r ) d r 2 ) r 3 μ K w ( d w ( r ) d r ) r 2 ) = 0 .
Applying DQM, the discrete forms of Equations (A1)–(A3) are as follows:
δ ( 1 ν 12 ν 21 ) j = 1 9 A i j ( 1 ) u 0 ( r ) + δ 2 ( 1 ν 12 ν 21 ) j = 1 9 A i j ( 1 ) w 0 ( r j ) j = 1 9 A i j ( 1 ) w 0 ( r j ) δ ν 12 α ( 1 ν 12 ν 21 ) j = 1 9 A i j ( 1 ) w 0 ( r j ) j = 1 9 A i j ( 1 ) w 0 ( r j ) δ α ( 1 ν 12 ν 21 ) j = 1 9 w 0 ( r j ) + δ ( 1 ν 12 ν 21 ) j = 1 9 r j A i j ( 2 ) u 0 ( r ) + δ ( 1 ν 12 ν 21 ) j = 1 9 r j A i j ( 2 ) w 0 ( r j ) j = 1 9 A i j ( 1 ) w 0 ( r j ) = 0 ,
1 12 ( j = 1 9 r j C i j ( 2 ) φ ( r j ) ( 1 ν 12 ν 21 ) + ν 12 α j = 1 9 A i j ( 1 ) φ ( r j ) ( 1 ν 12 ν 21 ) ν 12 α j = 1 9 1 r j φ ( r j ) ( 1 ν 12 ν 21 ) ) δ + 1 12 ( j = 1 9 A i j ( 1 ) w ( r j ) ( 1 ν 12 ν 21 ) + ν 12 α j = 1 9 1 r j φ ( r j ) ( 1 ν 12 ν 21 ) ) δ 1 12 ( ν 12 α ( j = 1 9 A i j ( 1 ) φ ( r j ) ) ( 1 ν 12 ν 21 ) + α j = 1 9 1 r j φ ( r j ) ( 1 ν 12 ν 21 ) ) δ K s β ( j = 1 9 r j φ ( r j ) + j = 1 9 A i j ( 1 ) w ( r j ) ) δ = 0 ,
K s β ( j = 1 9 A i j ( 1 ) φ ( r j ) + j = 1 9 A i j ( 2 ) w 0 ( r j ) ) + K s β ( j = 1 9 φ ( r j ) + j = 1 9 1 r j A i j ( 1 ) w 0 ( r j ) ) N ( j = 1 9 A i j ( 2 ) w 0 ( r j ) + j = 1 9 1 r j A i j ( 1 ) w 0 ( r j ) ) K w j = 1 9 w 0 ( r j ) δ μ ( N ( j = 1 9 A i j ( 4 ) w 0 ( r j ) + 2 ( j = 1 9 1 r j 3 A i j ( 1 ) w 0 ( r j ) ) 2 ( j = 1 9 1 r j 2 A i j ( 2 ) w 0 ( r j ) ) + j = 1 9 1 r A i j ( 3 ) w 0 ( r j ) ) 1 δ K w ( j = 1 9 A i j ( 2 ) w 0 ( r j ) ) N ( j = 1 9 1 r j A i j ( 3 ) w 0 ( r j ) + j = 1 9 1 r 3 A i j ( 1 ) w 0 ( r j ) j = 1 9 1 r 2 A i j ( 2 ) w 0 ( r j ) + 1 δ K w ( j = 1 9 1 r A i j ( 1 ) w 0 ( r j ) ) = 0
where α = E 2 E 1 ; β = G 13 E 1 are considered; also, K s = 5 / 6 .

References

  1. Ji, Y.; Choe, M.; Cho, B.; Song, S.; Yoon, J.; Ko, H.C.; Lee, T. Organic nonvolatile memory devices with charge trapping multilayer graphene film. Nanotechnology 2012, 23, 105202. [Google Scholar] [CrossRef]
  2. Bunch, J.S.; van der Zande, A.M.; Verbridge, S.S.; Frank, I.W.; Tanenbaum, D.M.; Parpia, J.M.; Craighead, H.G.; McEuen, P.L. Electromechanical Resonators from Graphene Sheets. Science 2007, 315, 490–493. [Google Scholar] [CrossRef] [PubMed]
  3. Sakhaee-Pour, A.; Ahmadian, M.T.; Vafai, A. Applications of single-layered graphene sheets as mass sensors and atomistic dust detectors. Solid State Commun. 2008, 145, 168–172. [Google Scholar] [CrossRef]
  4. Lu, G.; Ocola, L.E.; Chen, J. Reduced graphene oxide for room-temperature gas sensors. Nanotechnology 2009, 20, 445502. [Google Scholar] [CrossRef] [PubMed]
  5. Krishna, S.; Sreedhar, I.; Patel, C.M. Molecular dynamics simulation of polyamide-based materials—A review. Comput. Mater. Sci. 2021, 200, 110853. [Google Scholar] [CrossRef]
  6. Roudbari, M.A.; Jorshari, T.D.; Lü, C.; Ansari, R.; Kouzani, A.Z.; Amabili, M. A review of size-dependent continuum mechanics models for micro- and nano-structures. Thin-Walled Struct. 2022, 170, 108562. [Google Scholar] [CrossRef]
  7. Li, L.; Tang, H.; Hu, Y. The effect of thickness on the mechanics of nanobeams. Int. J. Eng. Sci. 2018, 123, 81–91. [Google Scholar] [CrossRef]
  8. Eringen, A.C. Nonlocal polar elastic continua. Int. J. Eng. Sci. 1972, 10, 1–16. [Google Scholar] [CrossRef]
  9. Triantafyllidis, N.; Aifantis, E.C. A gradient approach to localization of deformation. I. Hyperelastic materials. J. Elast. 1986, 16, 225–237. [Google Scholar] [CrossRef]
  10. 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]
  11. Narendar, S. Buckling analysis of micro-/nano-scale plates based on two-variable refined plate theory incorporating nonlocal scale effects. Compos. Struct. 2011, 93, 3093–3103. [Google Scholar] [CrossRef]
  12. Jamalpoor, A.; Ahmadi-Savadkoohi, A.; Hosseini, M.; Hosseini-Hashemi, S. Free vibration and biaxial buckling analysis of double magneto-electro-elastic nanoplate-systems coupled by a visco- Pasternak medium via nonlocal elasticity theory. Eur. J. Mech.-A/Solids 2017, 63, 84–98. [Google Scholar] [CrossRef]
  13. Li, L.; Hu, Y. Buckling analysis of size-dependent nonlinear beams based on a nonlocal strain gradient theory. Int. J. Eng. Sci. 2015, 97, 84–94. [Google Scholar] [CrossRef]
  14. Shaban, M.; Alibeigloo, A. Three dimensional vibration and bending analysis of carbon nanotubes embedded in elastic medium based on theory of elasticity. Lat. Am. J. Solids Struct. 2014, 11, 2122–2140. [Google Scholar] [CrossRef]
  15. Taher, H.R.D.; Omidi, M.; Zadpoor, A.A.; Nikooyan, A.A. Free vibration of circular and annular plates with variable thickness and different combinations of boundary conditions. J. Sound Vib. 2006, 296, 1084–1092. [Google Scholar] [CrossRef]
  16. Shaban, M.; Mazaheri, H. Size-dependent electro-static analysis of smart micro-sandwich panels with functionally graded core. Acta Mech. 2021, 232, 111–133. [Google Scholar] [CrossRef]
  17. Shaban, M.; Alipour, M.M. Semi-analytical solution for free vibration of thick functionally graded plates rested on elastic foundation with elastically restrained edge. Acta Mech. Solida Sin. 2011, 24, 340–354. [Google Scholar] [CrossRef]
  18. Ghayesh, M.H.; Farajpour, A. Nonlinear mechanics of nanoscale tubes via nonlocal strain gradient theory. Int. J. Eng. Sci. 2018, 129, 84–95. [Google Scholar] [CrossRef]
  19. Gui, Y.; Wu, R. Buckling analysis of embedded thermo-magneto-electro-elastic nano cylindrical shell subjected to axial load with nonlocal strain gradient theory. Mech. Res. Commun. 2023, 128, 104043. [Google Scholar] [CrossRef]
  20. Tanzadeh, H.; Amoushahi, H. Buckling analysis of orthotropic nanoplates based on nonlocal strain gradient theory using the higher-order finite strip method (H-FSM). Eur. J. Mech.-A/Solids 2022, 95, 104622. [Google Scholar] [CrossRef]
  21. Cuong-Le, T.; Nguyen, K.D.; Hoang-Le, M.; Sang-To, T.; Phan-Vu, P.; Wahab, M.A. Nonlocal strain gradient IGA numerical solution for static bending, free vibration and buckling of sigmoid FG sandwich nanoplate. Phys. B Condens. Matter 2022, 631, 413726. [Google Scholar] [CrossRef]
  22. Wang, P.; Gao, Z.; Pan, F.; Moradi, Z.; Mahmoudi, T.; Khadimallah, M.A. A couple of GDQM and iteration techniques for the linear and nonlinear buckling of bi-directional functionally graded nanotubes based on the nonlocal strain gradient theory and high-order beam theory. Eng. Anal. Bound. Elem. 2022, 143, 124–136. [Google Scholar] [CrossRef]
  23. Wang, S.; Kang, W.; Yang, W.; Zhang, Z.; Li, Q.; Liu, M.; Wang, X. Hygrothermal effects on buckling behaviors of porous bi-directional functionally graded micro-/nanobeams using two-phase local/nonlocal strain gradient theory. Eur. J. Mech.-A/Solids 2022, 94, 104554. [Google Scholar] [CrossRef]
  24. Esen, I.; Özmen, R. Thermal vibration and buckling of magneto-electro-elastic functionally graded porous nanoplates using nonlocal strain gradient elasticity. Compos. Struct. 2022, 296, 115878. [Google Scholar] [CrossRef]
  25. Tang, Y.; Qing, H. Elastic buckling and free vibration analysis of functionally graded Timoshenko beam with nonlocal strain gradient integral model. Appl. Math. Model. 2021, 96, 657–677. [Google Scholar] [CrossRef]
  26. Al-Furjan, M.S.H.; Farrokhian, A.; Keshtegar, B.; Kolahchi, R.; Trung, N.-T. Higher order nonlocal viscoelastic strain gradient theory for dynamic buckling analysis of carbon nanocones. Aerosp. Sci. Technol. 2020, 107, 106259. [Google Scholar] [CrossRef]
  27. Chwał, M.; Muc, A. Buckling and Free Vibrations of Nanoplates—Comparison of Nonlocal Strain and Stress Approaches. Appl. Sci. 2019, 9, 1409. [Google Scholar] [CrossRef]
  28. Fan, F.; Safaei, B.; Sahmani, S. Buckling and postbuckling response of nonlocal strain gradient porous functionally graded micro/nano-plates via NURBS-based isogeometric analysis. Thin-Walled Struct. 2021, 159, 107231. [Google Scholar] [CrossRef]
  29. Sadeghian, M.; Palevicius, A.; Janusas, G. Nonlocal Strain Gradient Model for the Nonlinear Static Analysis of a Circular/Annular Nanoplate. Micromachines 2023, 14, 1052. [Google Scholar] [CrossRef]
  30. Timoshenko, S.; Woinowsky-Krieger, S. Theory of Plates and Shells; McGraw-Hill: New York, NY, USA, 1959; Volume 2. [Google Scholar]
  31. Zghal, S.; Trabelsi, S.; Dammak, F. Post-buckling behavior of functionally graded and carbon-nanotubes based structures with different mechanical loadings. Mech. Based Des. Struct. Mach. 2022, 50, 2997–3039. [Google Scholar] [CrossRef]
  32. Zghal, S.; Dammak, F. Buckling responses of porous structural components with gradient power-based and sigmoid material variations under different types of compression loads. Compos. Struct. 2021, 273, 114313. [Google Scholar] [CrossRef]
  33. Trabelsi, S.; Zghal, S.; Dammak, F. A Finite Element Procedure for Thermal Buckling Analysis of Functionally Graded Shell Structures. In Design and Modeling of Mechanical Systems—IV; Springer International Publishing: Cham, Switzerland, 2020; pp. 409–416. [Google Scholar]
  34. Mehar, K.; Kumar Panda, S.; Devarajan, Y.; Choubey, G. Numerical buckling analysis of graded CNT-reinforced composite sandwich shell structure under thermal loading. Compos. Struct. 2019, 216, 406–414. [Google Scholar] [CrossRef]
  35. Do, V.N.V.; Lee, C.-H. Quasi-3D higher-order shear deformation theory for thermal buckling analysis of FGM plates based on a meshless method. Aerosp. Sci. Technol. 2018, 82–83, 450–465. [Google Scholar] [CrossRef]
  36. Szekrényes, A. Application of differential quadrature method to delaminated first-order shear deformable composite plates. Thin-Walled Struct. 2021, 166, 108028. [Google Scholar] [CrossRef]
  37. Duryodhana, D.; Waddar, S.; Bonthu, D.; Pitchaimani, J.; Powar, S.; Doddamani, M. Buckling and free vibrations behaviour through differential quadrature method for foamed composites. Results Eng. 2023, 17, 100894. [Google Scholar] [CrossRef]
  38. Han, J.; Li, L.; Jin, G.; Ma, W.; Feng, J.; Jia, H.; Chang, D. Qualitative Identification of the Static Pull-In and Fundamental Frequency of One-Electrode MEMS Resonators. Micromachines 2018, 9, 614. [Google Scholar] [CrossRef]
  39. Ren, Y.; Huo, R.; Zhou, D. Thermo-mechanical buckling analysis of non-uniformly heated rectangular plates with temperature-dependent material properties. Thin-Walled Struct. 2023, 186, 110653. [Google Scholar] [CrossRef]
  40. Thai, H.-T.; Kim, S.-E. A simple higher-order shear deformation theory for bending and free vibration analysis of functionally graded plates. Compos. Struct. 2013, 96, 165–173. [Google Scholar] [CrossRef]
  41. Sreehari, V.M. A detailed comparison study of first order and higher order shear deformation theories in the analysis of laminated composite plate. IOP Conf. Ser. Mater. Sci. Eng. 2017, 225, 012018. [Google Scholar] [CrossRef]
  42. Ambartsumian, S. On the theory of bending plates. Q. Appl. Math. 1958, 5, 69–77. [Google Scholar]
  43. Reddy, J.N. A Simple Higher-Order Theory for Laminated Composite Plates; American Society of Mechanical Engineers (ASME): New York, NY, USA, 1984. [Google Scholar]
  44. Reissner, E. On Tranverse Bending of Plates, Including the Effect of Transverse Shear Deformation; National Academy of Sciences, Engineering and Medicine: Washington, DC, USA, 1974. [Google Scholar]
  45. Touratier, M. An efficient standard plate theory. Int. J. Eng. Sci. 1991, 29, 901–916. [Google Scholar] [CrossRef]
  46. Soldatos, K. A transverse shear deformation theory for homogeneous monoclinic plates. Acta Mech. 1992, 94, 195–220. [Google Scholar] [CrossRef]
  47. Aydogdu, M. A new shear deformation theory for laminated composite plates. Compos. Struct. 2009, 89, 94–101. [Google Scholar] [CrossRef]
  48. Mantari, J.; Oktem, A.; Soares, C.G. A new higher order shear deformation theory for sandwich and composite laminated plates. Compos. Part B Eng. 2012, 43, 1489–1499. [Google Scholar] [CrossRef]
  49. Li, Q.; Wu, D.; Chen, X.; Liu, L.; Yu, Y.; Gao, W. Nonlinear vibration and dynamic buckling analyses of sandwich functionally graded porous plate with graphene platelet reinforcement resting on Winkler–Pasternak elastic foundation. Int. J. Mech. Sci. 2018, 148, 596–610. [Google Scholar] [CrossRef]
  50. Farajpour, A.; Arab Solghar, A.; Shahidi, A. Postbuckling analysis of multi-layered graphene sheets under non-uniform biaxial compression. Phys. E Low-Dimens. Syst. Nanostructures 2013, 47, 197–206. [Google Scholar] [CrossRef]
  51. Pradhan, S.C.; Phadikar, J.K. 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]
  52. Bellman, R.; Casti, J. Differential quadrature and long-term integration. J. Math. Anal. Appl. 1971, 34, 235–238. [Google Scholar] [CrossRef]
  53. Farajpour, A.; Mohammadi, M.; Shahidi, A.R.; Mahzoon, M. Axisymmetric buckling of the circular graphene sheets with the nonlocal continuum plate model. Phys. E Low-Dimens. Syst. Nanostructures 2011, 43, 1820–1825. [Google Scholar] [CrossRef]
  54. Xu, X.-J.; Wang, X.-C.; Zheng, M.-L.; Ma, Z. Bending and buckling of nonlocal strain gradient elastic beams. Compos. Struct. 2017, 160, 366–377. [Google Scholar] [CrossRef]
  55. Sepahi, O.; Forouzan, M.R.; Malekzadeh, P. Thermal buckling and postbuckling analysis of functionally graded annular plates with temperature-dependent material properties. Mater. Des. 2011, 32, 4030–4041. [Google Scholar] [CrossRef]
Figure 1. (a) The annular nanoplate with Pasternak–Winkler elastic foundation; (b) the graphene plate with uniform buckling loads.
Figure 1. (a) The annular nanoplate with Pasternak–Winkler elastic foundation; (b) the graphene plate with uniform buckling loads.
Micromachines 14 01790 g001
Figure 2. A schematic of the bilayer nanoplate considering the Pasternak and Winkler elastic foundations.
Figure 2. A schematic of the bilayer nanoplate considering the Pasternak and Winkler elastic foundations.
Micromachines 14 01790 g002
Figure 3. The effect of the number of nodes on the dimensionless mechanical buckling load.
Figure 3. The effect of the number of nodes on the dimensionless mechanical buckling load.
Micromachines 14 01790 g003
Figure 4. The nondimensional buckling load of the annular nanoplate with respect to the nonlocal parameter in diverse boundary conditions.
Figure 4. The nondimensional buckling load of the annular nanoplate with respect to the nonlocal parameter in diverse boundary conditions.
Micromachines 14 01790 g004
Figure 5. Nondimensional buckling loads of the annular nanoplate with respect to the nonlocal parameter for different values of elastic foundations at various boundary conditions (including C-C and S–S).
Figure 5. Nondimensional buckling loads of the annular nanoplate with respect to the nonlocal parameter for different values of elastic foundations at various boundary conditions (including C-C and S–S).
Micromachines 14 01790 g005
Figure 6. The nondimensional buckling load of the annular nanoplate with respect to the nonlocal parameter for different values of elastic foundations at various boundary conditions (including C–S and S–C).
Figure 6. The nondimensional buckling load of the annular nanoplate with respect to the nonlocal parameter for different values of elastic foundations at various boundary conditions (including C–S and S–C).
Micromachines 14 01790 g006
Figure 7. The nondimensional buckling load of the annular nanoplate with respect to the nonlocal parameter for FSDT and HSDT considering h/r = 0.05 for the simply supported boundary condition.
Figure 7. The nondimensional buckling load of the annular nanoplate with respect to the nonlocal parameter for FSDT and HSDT considering h/r = 0.05 for the simply supported boundary condition.
Micromachines 14 01790 g007
Figure 8. The nondimensional buckling load of the annular nanoplate with respect to the nonlocal parameter for FSDT and HSDT considering h/r = 0.06 for the simply supported boundary condition.
Figure 8. The nondimensional buckling load of the annular nanoplate with respect to the nonlocal parameter for FSDT and HSDT considering h/r = 0.06 for the simply supported boundary condition.
Micromachines 14 01790 g008
Figure 9. The nondimensional buckling load of the annular nanoplate versus the radius ratios.
Figure 9. The nondimensional buckling load of the annular nanoplate versus the radius ratios.
Micromachines 14 01790 g009
Figure 10. The critical buckling temperature for the annular nanoplate with respect to the number of nodes.
Figure 10. The critical buckling temperature for the annular nanoplate with respect to the number of nodes.
Micromachines 14 01790 g010
Figure 11. The critical buckling temperature for the annular nanoplate versus the nonlocal parameter.
Figure 11. The critical buckling temperature for the annular nanoplate versus the nonlocal parameter.
Micromachines 14 01790 g011
Figure 12. Critical buckling temperature for single-layer and bilayer annular nanoplate with respect to the nonlocal parameter for the C–C condition.
Figure 12. Critical buckling temperature for single-layer and bilayer annular nanoplate with respect to the nonlocal parameter for the C–C condition.
Micromachines 14 01790 g012
Figure 13. Critical buckling temperature for single-layer and bilayer annular nanoplate versus the nonlocal parameter for S–S condition.
Figure 13. Critical buckling temperature for single-layer and bilayer annular nanoplate versus the nonlocal parameter for S–S condition.
Micromachines 14 01790 g013
Figure 14. Critical buckling temperature for single-layer and bilayer annular nanoplate versus the nonlocal parameter for S–C condition.
Figure 14. Critical buckling temperature for single-layer and bilayer annular nanoplate versus the nonlocal parameter for S–C condition.
Micromachines 14 01790 g014
Figure 15. Critical buckling temperature for single-layer and bilayer annular nanoplate versus the nonlocal parameter for C–S condition.
Figure 15. Critical buckling temperature for single-layer and bilayer annular nanoplate versus the nonlocal parameter for C–S condition.
Micromachines 14 01790 g015
Table 1. Recommended functions used in HSDT by various authors.
Table 1. Recommended functions used in HSDT by various authors.
Model g ( z ) Function
Ambartsumian [42] 1 6 z 3 + h 2 8 z
Reddy [43] 4 3 h 2 z 3 + z
Reissner [44] 5 3 h 2 z 3 + 5 4 z
Touratier [45] h π sin π z h
Soldatos [46] h sinh z h z cosh 1 2
Aydogdu [47] z e 2 z h 2
Mantari [48] h π sin π z h e m cos π z h + m π h z         ,         m 0
Table 2. Comparison of the nondimensional buckling loads of the circular nanoplate gained by the present study with a reference.
Table 2. Comparison of the nondimensional buckling loads of the circular nanoplate gained by the present study with a reference.
RadiusReferenceThe Percentage of the Nondimensional Buckling Loads
μ = 0μ = 0.25μ = 1μ = 2.25μ = 4
4Ref. [53]0.9430 0.7671 0.4918 0.3077 0.2019
Present FSDT results0.93530.75540.47940.30080.1943
Present HSDT results0.92110.74350.34330.24320.0913
6Ref. [53]0.41910.3803 0.2977 0.2186 0.1593
Present FSDT results0.42220.38220.29750.21730.1577
Present HSDT results0.415930.37680.29380.21070.1446
8Ref. [53]0.2358 0.2230 0.1918 0.1555 0.1229
Present FSDT results0.23780.22550.19330.15620.1231
Present HSDT results0.23520.22230.19090.15450.1218
10Ref. [53]0.1509 0.1455 0.1316 0.1134 0.0951
Present FSDT results0.15320.14760.13320.11450.0957
Present HSDT results0.15090.14550.13140.11320.0947
Table 3. Comparison of the dimensionless buckling load gained by considering various functions of HSDT for a circular nanoplate under both simply supported and clamped boundary conditions.
Table 3. Comparison of the dimensionless buckling load gained by considering various functions of HSDT for a circular nanoplate under both simply supported and clamped boundary conditions.
ConditionFunctionμ = 0μ = 1μ = 2μ = 3μ = 4
Cg1(z)0.012960.007200.004850.004000.00355
g2(z)0.012950.007090.004800.003960.00353
g3(z)0.012960.007170.004840.003990.00355
g4(z)0.012970.007220.004860.004000.00356
g5(z)0.012920.007060.004780.003950.00352
Sg1(z)0.010960.007750.005210.004260.00377
g2(z)0.010960.007670.005170.004230.00375
g3(z)0.010960.0077320.005200.004260.00377
g4(z)0.010960.0077630.005220.004270.00377
g5(z)0.010950.007600.005130.004210.00373
Table 4. Comparison of the dimensionless thermal buckling parameter ( λ T ) gained by the present study (HSDT and FSDT) with the reference for the annular nanoplate under clamped boundary conditions.
Table 4. Comparison of the dimensionless thermal buckling parameter ( λ T ) gained by the present study (HSDT and FSDT) with the reference for the annular nanoplate under clamped boundary conditions.
r i r o = 0.2 r i r o = 0.4 r i r o = 0.6
h r o Ref. [55]Present Study (HSDT)Present Study (FSDT)Ref. [55]Present Study (HSDT)Present Study (FSDT)Ref. [55]Present Study (HSDT)Present Study (FSDT)
S–S0.0518.54718.49618.51528.64128.6428.64260.19860.20260.198
0.117.72517.65017.68026.92126.92226.92153.25153.28353.252
0.1516.51016.42716.46824.47224.48224.47244.66344.75544.663
0.215.06914.98715.03521.70921.73721.70936.43636.61836.436
S–C0.0540.71740.68740.63861.16761.10761.166122.304121.893122.302
0.137.44537.42837.38653.96053.88953.96195.93395.16095.934
0.1533.02233.02332.97945.10845.00145.10870.61268.08570.613
0.228.33728.34028.30636.68935.86736.68951.60651.53951.606
C–S0.0527.96726.948727.93349.95648.89849.956109.571107.756109.574
0.125.91124.35625.94844.22642.68944.22886.18883.33286.189
0.1523.08521.29423.13337.13635.2437.13863.63560.69663.636
0.220.03318.18420.07830.34928.4630.35046.67446.29946.674
Table 5. The differences in percentage errors between the present study with the Ref. [55].
Table 5. The differences in percentage errors between the present study with the Ref. [55].
r i r o = 0.2 r i r o = 0.4 r i r o = 0.6
h / r o E 1 ( % ) E 2 ( % ) E 1 ( % ) E 2 ( % ) E 1 ( % ) E 2 ( % )
S–S0.050.27490.17250.00340.00340.00660
0.10.42310.25380.003700.06000.0018
0.150.50270.25430.040800.20590
0.20.54410.22560.128900.49950
S–C0.050.07360.19400.09800.00160.33600.0016
0.10.04530.15750.13150.00180.80570.0010
0.150.00300.13020.237203.57870.0014
0.20.01050.10932.240400.12980
C–S0.053.64100.12152.117801.65640.0027
0.16.00130.14273.47530.00453.31360.0011
0.157.75820.20795.10550.00534.61850.0015
0.29.22970.22466.22420.00320.80340
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Sadeghian, M.; Palevicius, A.; Janusas, G. Nonlinear Thermal/Mechanical Buckling of Orthotropic Annular/Circular Nanoplate with the Nonlocal Strain Gradient Model. Micromachines 2023, 14, 1790. https://doi.org/10.3390/mi14091790

AMA Style

Sadeghian M, Palevicius A, Janusas G. Nonlinear Thermal/Mechanical Buckling of Orthotropic Annular/Circular Nanoplate with the Nonlocal Strain Gradient Model. Micromachines. 2023; 14(9):1790. https://doi.org/10.3390/mi14091790

Chicago/Turabian Style

Sadeghian, Mostafa, Arvydas Palevicius, and Giedrius Janusas. 2023. "Nonlinear Thermal/Mechanical Buckling of Orthotropic Annular/Circular Nanoplate with the Nonlocal Strain Gradient Model" Micromachines 14, no. 9: 1790. https://doi.org/10.3390/mi14091790

APA Style

Sadeghian, M., Palevicius, A., & Janusas, G. (2023). Nonlinear Thermal/Mechanical Buckling of Orthotropic Annular/Circular Nanoplate with the Nonlocal Strain Gradient Model. Micromachines, 14(9), 1790. https://doi.org/10.3390/mi14091790

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