# Thermally Induced Asymmetric Buckling of Circular Monolayer Graphene

Authors: Haw-Long Lee; Win-Jin Chang
Journal: Journal of Nanomaterials
 (2013)

---

## Abstract

The asymmetric buckling behaviors of circular monolayer graphene with clamped boundary condition subjected to temperature change are analytically studied based on the nonlocal elasticity theory, including the small length effect. The axisymmetrical and asymmetric critical buckling temperatures and mode shape of different order modes are obtained. According to the analysis, the asymmetric critical buckling temperature of monolayer graphene is larger than the axisymmetric one. The axisymmetrical and asymmetric critical buckling temperatures decrease with increasing nonlocal parameter. In addition, nodal diametrical lines and nodal circles can be found from the modal shapes. In order to avoid destruction of the sensors due to buckling of the structure, they can be placed at the nodal diametrical lines or nodal circles.


---

## Body

## 1. Introduction

Graphene has recently attracted significant research interest due to its potential applications in nanobiological devices and nanoelectromechanical systems such as single-electron transistor [1], actuators [2], and photodetectors [3]. Because graphene is only one atom layer thick, it is easily amenable to external influences, including thermal loading and mechanical deformation [4–8]. The external mechanical loading can affect its electronic structure and physical properties. For example, there were significant changes in graphene’s band structure for uniaxial strains of up to 15% [5]. The magnetic behavior of multilayer graphene can be affected by mechanical transverse deformation [7].In the last years, continuum theory has received the attention of many researchers for analysis of nanostructures. This is because available experiments on nanoscale materials are difficult to perform, and molecular dynamic simulations are highly time-consuming. Accordingly, various size-dependent continuum theories such as strain gradient theory [9], couple stress theory [10], and nonlocal elasticity theory [11] have been proposed. Among these theories, nonlocal elasticity theory including the small-scale effect has often been used to study the buckling and vibration responses of nanostructures [12–16]. The theory was initiated by Eringen [17] in the early 1970s and stated that the stresses at a point depend on strains at all the points of the body.This paper focuses on the buckling behavior of graphene. In recent years, some researchers have studied the buckling of graphene by using different methods. For example, Frank et al. [18] experimentally studied the compression buckling strain of a graphene sheet. Neek-Amal and Peeters [19] investigated the stability of circular monolayer graphene subjected to a radial load using molecular dynamics simulations. Natsuki et al. [20] studied the buckling properties of circular double-layered graphene sheets using plate theory and found that buckling stability is significantly affected by the buckling mode shapes.In this paper, thermally induced asymmetric buckling of circular monolayer graphene is studied based on the nonlocal elasticity theory by using the analytical method. According to the analysis, the axisymmetrical and asymmetric buckling temperatures and mode shape for the graphene with various modes are obtained. In addition, nodal diametrical lines and nodal circles in the buckling mode shape are also investigated.


## 2. Physical Model and Mathematical Formulation

### 2.1. Governing Equation and Boundary Conditions

Consider that a circular graphene sheet is clamped on its edge, which is modeled as a clamped circular plate with the radiusR and thickness h subjected to temperature change T as depicted in Figure 1. Based on the nonlocal elasticity theory, the nonlocal constitutive relation of circular plate under the thermal load is given by [21]
(1)D∇4w(r,θ)+[1-(e0a)2∇2]NT∇2w(r,θ)=0,
where w is the displacement along the thickness of the graphene. r and θ are the radial and circumferential coordinates. ∇2=((∂2/∂r2)+(1/r)(∂/∂r)+(1/r2)(∂2/∂θ2)) is the Laplacian operator in polar coordinates. e0a is the nonlocal parameter revealing the nanoscale size effect. NT and D are the thermal buckling load and flexural bending rigidity of graphene, respectively. They are
(2)NT=EαT1-μh,D=Eh312(1-μ2),
where E,  μ, and α are Young’s modulus, Poisson’s ratio, and thermal expansion coefficient of graphene, respectively.Figure 1
A clamped circular graphene sheet with the radiusR under temperature change T.The corresponding boundary conditions are(3)w(R)=0,∂w(R)∂r=0.
The nondimensional variables are introduced and defined as follows:
(4)ξ=rR,w-=wR,ε=e0aR,T*=12αTR2(1+μ)h2,
where ξ, w-, ε, and T* denote the nondimensional radius, outplane displacement, nonlocal parameter, and buckling temperature.Using the dimensionless variables given by (4), the governing equation and associated boundary conditions can be simplified to the following dimensionless form:
(5)∇-2(∇-2w-+λ2w-)=0,(6)w-(1)=0,∂w-(1)∂ξ=0,
where
(7)∇-2=(∂2∂ξ2+∂ξ∂ξ+∂2ξ2∂θ2),(8)λ2=T*1-ε2T*.


### 2.2. Asymmetric Buckling Mode

It is well known that the normal vibration of an elastic linear system is harmonic. The displacement of monolayer graphene for a harmonic vibration can therefore be separated as follows:(9)w-(ξ,θ)=f(ξ)cosnθ,n=0,1,2,3,…,
where n is the number of nodal diametrical lines.Substituting the relation of (9) into (5), we have
(10)(d2dξ2+dξdξ-n2ξ2)[d2dξ2+dξdξ+(λ2-n2ξ2)]f(ξ)=0.
The general solution of the earlier equation should be expressed as
(11)f(ξ)=A1ξn+A2ξ-n+A3Jn(λξ)+A4Yn(λξ),n=1,2,3,
where Jn and  Yn are the Bessel functions of the first and second kinds of order n.  A1,…,A4 are arbitrary constants depending on the boundary conditions.Substituting the boundary conditions given by (6) into (11), we can obtain the characteristic functions as
(12)Jn+1(λmn)=0,
where λmn are the roots of the Bessel function of the first kind of order n+1 for mode m.Once the parameterλmn is given, the following nondimensional buckling temperature can be obtained from (8):
(13)Tmn*=λmn21+ε2λmn2.
In addition, the nondimensional buckling mode shape can be expressed as
(14)w-mn(ξ,θ)=AmnJn(λmnξ)cosnθ.
There are n nodal diametrical lines, which are straight lines with no curve and no displacement. In addition, there are m nodal circles, which are no displacement. The radius of nodal circle ρ can be solved by
(15)Jn(λmnρ)=0.


### 2.3. Axisymmetric Buckling Mode

For the axisymmetric buckling mode (i.e.,n=0), using (10), the differential equation can be expressed as
(16)(d2dξ2+dξdξ)(d2dξ2+dξdξ+λ2)f(ξ)=0.
The general solution can be written as
(17)fs(ξ)=B1+B2lnξ+B3J0(λξ)+B4Y0(λξ),
where J0 and Y0 are the Bessel functions of the first and second kinds of order zero and B1,…,B4 are the arbitrary constants.Applying the boundary conditions, the characteristic function,J1(λm0)=0, can be obtained. The result is the same as that derived by Farajpour et al. [21].Using (14)-(15) and setting n=0, the nondimensional buckling temperature, buckling shape, and radii of nodal circles of axisymmetric monolayer graphene can also be determined. It is noted that there are no nodal diametrical lines, but there are m nodal circles for axisymmetric buckling mode.


## 2.1. Governing Equation and Boundary Conditions

Consider that a circular graphene sheet is clamped on its edge, which is modeled as a clamped circular plate with the radiusR and thickness h subjected to temperature change T as depicted in Figure 1. Based on the nonlocal elasticity theory, the nonlocal constitutive relation of circular plate under the thermal load is given by [21]
(1)D∇4w(r,θ)+[1-(e0a)2∇2]NT∇2w(r,θ)=0,
where w is the displacement along the thickness of the graphene. r and θ are the radial and circumferential coordinates. ∇2=((∂2/∂r2)+(1/r)(∂/∂r)+(1/r2)(∂2/∂θ2)) is the Laplacian operator in polar coordinates. e0a is the nonlocal parameter revealing the nanoscale size effect. NT and D are the thermal buckling load and flexural bending rigidity of graphene, respectively. They are
(2)NT=EαT1-μh,D=Eh312(1-μ2),
where E,  μ, and α are Young’s modulus, Poisson’s ratio, and thermal expansion coefficient of graphene, respectively.Figure 1
A clamped circular graphene sheet with the radiusR under temperature change T.The corresponding boundary conditions are(3)w(R)=0,∂w(R)∂r=0.
The nondimensional variables are introduced and defined as follows:
(4)ξ=rR,w-=wR,ε=e0aR,T*=12αTR2(1+μ)h2,
where ξ, w-, ε, and T* denote the nondimensional radius, outplane displacement, nonlocal parameter, and buckling temperature.Using the dimensionless variables given by (4), the governing equation and associated boundary conditions can be simplified to the following dimensionless form:
(5)∇-2(∇-2w-+λ2w-)=0,(6)w-(1)=0,∂w-(1)∂ξ=0,
where
(7)∇-2=(∂2∂ξ2+∂ξ∂ξ+∂2ξ2∂θ2),(8)λ2=T*1-ε2T*.


## 2.2. Asymmetric Buckling Mode

It is well known that the normal vibration of an elastic linear system is harmonic. The displacement of monolayer graphene for a harmonic vibration can therefore be separated as follows:(9)w-(ξ,θ)=f(ξ)cosnθ,n=0,1,2,3,…,
where n is the number of nodal diametrical lines.Substituting the relation of (9) into (5), we have
(10)(d2dξ2+dξdξ-n2ξ2)[d2dξ2+dξdξ+(λ2-n2ξ2)]f(ξ)=0.
The general solution of the earlier equation should be expressed as
(11)f(ξ)=A1ξn+A2ξ-n+A3Jn(λξ)+A4Yn(λξ),n=1,2,3,
where Jn and  Yn are the Bessel functions of the first and second kinds of order n.  A1,…,A4 are arbitrary constants depending on the boundary conditions.Substituting the boundary conditions given by (6) into (11), we can obtain the characteristic functions as
(12)Jn+1(λmn)=0,
where λmn are the roots of the Bessel function of the first kind of order n+1 for mode m.Once the parameterλmn is given, the following nondimensional buckling temperature can be obtained from (8):
(13)Tmn*=λmn21+ε2λmn2.
In addition, the nondimensional buckling mode shape can be expressed as
(14)w-mn(ξ,θ)=AmnJn(λmnξ)cosnθ.
There are n nodal diametrical lines, which are straight lines with no curve and no displacement. In addition, there are m nodal circles, which are no displacement. The radius of nodal circle ρ can be solved by
(15)Jn(λmnρ)=0.


## 2.3. Axisymmetric Buckling Mode

For the axisymmetric buckling mode (i.e.,n=0), using (10), the differential equation can be expressed as
(16)(d2dξ2+dξdξ)(d2dξ2+dξdξ+λ2)f(ξ)=0.
The general solution can be written as
(17)fs(ξ)=B1+B2lnξ+B3J0(λξ)+B4Y0(λξ),
where J0 and Y0 are the Bessel functions of the first and second kinds of order zero and B1,…,B4 are the arbitrary constants.Applying the boundary conditions, the characteristic function,J1(λm0)=0, can be obtained. The result is the same as that derived by Farajpour et al. [21].Using (14)-(15) and setting n=0, the nondimensional buckling temperature, buckling shape, and radii of nodal circles of axisymmetric monolayer graphene can also be determined. It is noted that there are no nodal diametrical lines, but there are m nodal circles for axisymmetric buckling mode.


## 3. Results and Discussion

To study the effect of small length scale on the buckling behaviors of circular monolayer graphene, critical buckling temperature and mode shape are analyzed for different mode numbers and nonlocal parameters. Figure2 shows the nondimensional critical buckling temperature as a function of nonlocal parameters for asymmetric monolayer graphene with clamped boundary conditions with n=1 for the first five modes. The parameter value of e0a/R=0 implies that the nonlocal effect is neglected. It can be seen that the effect of nonlocal parameter on the critical buckling temperature is significant, especially at higher-order modes. The critical buckling temperature decreases with increasing nonlocal parameter. This is because the internal interaction force increases as the nonlocal parameter increases.Figure 2
The critical buckling temperature as a function of nonlocal parameter for asymmetric monolayer graphene withn=1.Figure3 depicts the critical buckling temperature ratio for asymmetric (i.e., n=1) and axisymmetric (i.e., n=0) monolayer graphene at different modes. The critical buckling temperature of asymmetric monolayer graphene is larger than that of axisymmetric ones. The discrepancy is larger for a lower-order mode. Figures 4(a) and 4(b) illustrate the axisymmetric (i.e., n=0) buckling mode shape of circular monolayer graphene with e0a/R=0.05 for m=1 and m=2, respectively. In practical applications, the lower buckling modes are easier to reach than the higher-order modes. Therefore, the parameters m=1 and m=2 are selected in the analysis. The nodal circles r=const can be seen for the figures. There are m circular nodes for different values of m. According to the calculation, the buckling temperatures are T10*=14.16 and T20*=43.86 for m=1 and m=2, respectively.Figure 3
The critical buckling temperature ratio for asymmetric and axisymmetric monolayer graphene.Axisymmetric buckling mode shape of circular monolayer graphene withe0a/R=0.05 for (a) m=1 and (b) m=2.
(a)
(b)Figure5 shows the asymmetric buckling mode shape and nodal lines of circular monolayer graphene with e0a/R=0.05 for different values of m and n. The nodal lines are straight diametrical lines θ=const and circles r=const. There are m nodal circulars and n straight diametrical lines for the case analyzed in the figure. The sensors used to measure modal data can be placed at the nodal diametrical lines and nodal circles in order to avoid destruction of the sensors due to structure buckling. The buckling temperatures of T11*=24.73 for m=1 and n=1, T21*=60.19 for m=2 and n=1, T12*=36.95 for m=1 and n=2, and T22*=76.9 for n=2 and m=2 are obtained from this calculation.Asymmetric buckling mode shape of circular monolayer graphene withe0a/R=0.05 for (a) m=1, n=1, (b) m=2, n=1, (c) m=1, n=2, and (d) m=2, n=2.
(a)
(b)
(c)
(d)Figure6 depicts the nondimensional radius of nodal circle of mode 1 as a function of nonlocal parameter. It can be seen that the radius of nodal circle obviously increases with increasing the nonlocal parameter. The nondimensional radii of nodal circles as a function of nonlocal parameter with m=2 are listed in Table 1. There are two nodal circles for the case of m=2. The radii of the two nodal circles increase with increasing the value of n and nonlocal parameter. The large one can easily exceed the range of the graphene when the nonlocal parameter is too large. For the case of larger n and nonlocal parameter, both circles may be out of the range.Table 1
Nondimensional radii of nodal circles as a function of nonlocal parameter withm=2.
e
0
a
/
R
n
=
0
n
=
1
n
=
2
S
L
S
L
S
L
0
0.343
0.787
0.455
0.833
0.526
0.862
0.01
0.344
0.789
0.457
0.836
0.529
0.866
0.02
0.346
0.795
0.462
0.845
0.536
0.879
0.03
0.350
0.804
0.470
0.860
0.548
0.899
0.04
0.356
0.817
0.480
0.879
0.565
0.926
0.05
0.363
0.834
0.494
0.904
0.585
0.960
0.06
0.372
0.854
0.510
0.934
0.610
0.999
0.07
0.382
0.877
0.528
0.967
0.637
0.08
0.393
0.902
0.549
0.668
0.09
0.405
0.931
0.571
0.700
0.1
0.419
0.595
0.735
0.11
0.433
0.620
0.772
0.12
0.448
0.647
0.810
0.13
0.464
0.675
0.850
0.14
0.480
0.704
0.891
0.15
0.498
0.733
0.933
0.16
0.515
0.764
0.17
0.534
0.795
0.18
0.552
0.826
0.19
0.571
0.859
0.2
0.591
0.891
S: radius of smaller nodal circle.
L: radius of lager nodal circle.Figure 6
Nondimensional radius of nodal circle as a function of nonlocal parameter.


## 4. Conclusions

The effect of the small length scale on the critical buckling temperature and mode shape of circular monolayer graphene with clamped boundary condition was analyzed using the nonlocal elasticity theory. The results showed that the critical buckling temperature decreased with increasing nonlocal parameter. The critical buckling temperature of asymmetric monolayer graphene was larger than that of axisymmetric one. It can be observed that there werem nodal circulars and n straight diametrical lines for the asymmetric buckling of circular monolayer graphene. In addition, the radii of nodal circles obviously increased with increasing the nonlocal parameter.



---

*Source: 416189-2013-09-05.xml*