Journal of Nanomaterials

Volume 2013 (2013), Article ID 416189, 6 pages

http://dx.doi.org/10.1155/2013/416189

## Thermally Induced Asymmetric Buckling of Circular Monolayer Graphene

Department of Mechanical Engineering, Kun Shan University, Tainan 71003, Taiwan

Received 25 June 2013; Revised 4 August 2013; Accepted 7 August 2013

Academic Editor: Sheng-Rui Jian

Copyright © 2013 Haw-Long Lee and Win-Jin Chang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### 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.

#### 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 radius and thickness subjected to temperature change 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] where is the displacement along the thickness of the graphene. and are the radial and circumferential coordinates. is the Laplacian operator in polar coordinates. is the nonlocal parameter revealing the nanoscale size effect. and are the thermal buckling load and flexural bending rigidity of graphene, respectively. They are where , , and are Young’s modulus, Poisson’s ratio, and thermal expansion coefficient of graphene, respectively.

The corresponding boundary conditions are The nondimensional variables are introduced and defined as follows: where , , , and 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: where

##### 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: where is the number of nodal diametrical lines.

Substituting the relation of (9) into (5), we have The general solution of the earlier equation should be expressed as where and are the Bessel functions of the first and second kinds of order . are arbitrary constants depending on the boundary conditions.

Substituting the boundary conditions given by (6) into (11), we can obtain the characteristic functions as where are the roots of the Bessel function of the first kind of order for mode .

Once the parameter is given, the following nondimensional buckling temperature can be obtained from (8): In addition, the nondimensional buckling mode shape can be expressed as There are nodal diametrical lines, which are straight lines with no curve and no displacement. In addition, there are nodal circles, which are no displacement. The radius of nodal circle can be solved by

##### 2.3. Axisymmetric Buckling Mode

For the axisymmetric buckling mode (i.e., ), using (10), the differential equation can be expressed as The general solution can be written as where and are the Bessel functions of the first and second kinds of order zero and are the arbitrary constants.

Applying the boundary conditions, the characteristic function, , can be obtained. The result is the same as that derived by Farajpour et al. [21].

Using (14)-(15) and setting , 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 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. Figure 2 shows the nondimensional critical buckling temperature as a function of nonlocal parameters for asymmetric monolayer graphene with clamped boundary conditions with for the first five modes. The parameter value of 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 3 depicts the critical buckling temperature ratio for asymmetric (i.e., ) and axisymmetric (i.e., ) 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., ) buckling mode shape of circular monolayer graphene with for and , respectively. In practical applications, the lower buckling modes are easier to reach than the higher-order modes. Therefore, the parameters and are selected in the analysis. The nodal circles can be seen for the figures. There are circular nodes for different values of . According to the calculation, the buckling temperatures are and for and , respectively.

Figure 5 shows the asymmetric buckling mode shape and nodal lines of circular monolayer graphene with for different values of and . The nodal lines are straight diametrical lines and circles . There are nodal circulars and 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 for and , for and , for and , and for and are obtained from this calculation.

Figure 6 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 are listed in Table 1. There are two nodal circles for the case of . The radii of the two nodal circles increase with increasing the value of 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 and nonlocal parameter, both circles may be out of the range.

#### 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 were nodal circulars and 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.

#### Acknowledgment

The authors wish to thank the National Science Council of the Republic of China in Taiwan for providing financial support for this study under Projects NSC 101–2221-E-168-006 and NSC 101–2221-E-168 −009.

#### References

- C. Stampfer, E. Schurtenberger, F. Molitor, J. Güttinger, T. Ihn, and K. Ensslin, “Tunable graphene single electron transistor,”
*Nano Letters*, vol. 8, no. 8, pp. 2378–2383, 2008. View at Publisher · View at Google Scholar · View at Scopus - S. Park, J. An, J. W. Suk, and R. S. Ruoff, “Graphene-based actuators,”
*Small*, vol. 6, no. 2, pp. 210–212, 2010. View at Publisher · View at Google Scholar · View at Scopus - T. Mueller, F. Xia, and P. Avouris, “Graphene photodetectors for high-speed optical communications,”
*Nature Photonics*, vol. 4, no. 5, pp. 297–301, 2010. View at Publisher · View at Google Scholar · View at Scopus - W. Bao, F. Miao, Z. Chen et al., “Controlled ripple texturing of suspended graphene and ultrathin graphite membranes,”
*Nature Nanotechnology*, vol. 4, no. 9, pp. 562–566, 2009. View at Publisher · View at Google Scholar · View at Scopus - F. Guinea, M. I. Katsnelson, and A. K. Geim, “Energy gaps and a zero-field quantum hall effect in graphene by strain engineering,”
*Nature Physics*, vol. 6, no. 1, pp. 30–33, 2010. View at Publisher · View at Google Scholar · View at Scopus - V. M. Pereira, A. H. Castro Neto, and N. M. R. Peres, “Tight-binding approach to uniaxial strain in graphene,”
*Physical Review B*, vol. 80, no. 4, Article ID 045401, 2009. View at Publisher · View at Google Scholar · View at Scopus - E. V. Castro, M. P. López-Sancho, and M. A. H. Vozmediano, “Effect of pressure on the magnetism of bilayer graphene,”
*Physical Review B*, vol. 84, no. 7, Article ID 075432, 2011. View at Publisher · View at Google Scholar · View at Scopus - H. S. Shen, L. Shen, and C. L. Zhang, “Nonlocal plate model for nonlinear bending of single-layer graphene sheets subjected to transverse loads in thermal environments,”
*Applied Physics A*, vol. 103, no. 1, pp. 103–112, 2011. View at Publisher · View at Google Scholar · View at Scopus - D. C. C. Lam, F. Yang, A. C. M. Chong, J. Wang, and P. Tong, “Experiments and theory in strain gradient elasticity,”
*Journal of the Mechanics and Physics of Solids*, vol. 51, no. 8, pp. 1477–1508, 2003. View at Publisher · View at Google Scholar · View at Scopus - F. Yang, A. C. M. Chong, D. C. C. Lam, and P. Tong, “Couple stress based strain gradient theory for elasticity,”
*International Journal of Solids and Structures*, vol. 39, no. 10, pp. 2731–2743, 2002. View at Publisher · View at Google Scholar · View at Scopus - A. C. Eringen,
*Nonlocal Continuum Field Theories*, Springer, New York, NY, USA, 2002. - H. L. Lee and W. J. Chang, “Vibration analysis of fluid-conveying double-walled carbon nanotubes based on nonlocal elastic theory,”
*Journal of Physics Condensed Matter*, vol. 21, no. 11, Article ID 115302, 2009. View at Publisher · View at Google Scholar · View at Scopus - W. J. Chang and H. L. Lee, “Buckling temperature of a single-walled carbon nanotube using nonlocal Timoshenko beam model,”
*Journal of Computational and Theoretical Nanoscience*, vol. 7, no. 11, pp. 2367–2371, 2010. View at Publisher · View at Google Scholar · View at Scopus - H. L. Lee, J. C. Hsu, and W. J. Chang, “Frequency shift of carbon-nanotube-based mass sensor using nonlocal elasticity theory,”
*Nanoscale Research Letters*, vol. 5, no. 11, pp. 1774–1778, 2010. View at Publisher · View at Google Scholar · View at Scopus - H. L. Lee and W. J. Chang, “Dynamic modelling of a single-walled carbon nanotube for nanoparticle delivery,”
*Proceedings of the Royal Society A*, vol. 467, no. 2127, pp. 860–868, 2011. View at Publisher · View at Google Scholar · View at Scopus - H. L. Lee and W. J. Chang, “Vibration analysis of viscoelastic carbon nanotubes,”
*Micro & Nano Letters*, vol. 7, no. 12, pp. 1308–1316, 2012. View at Google Scholar - A. C. Eringen, “Linear theory of nonlocal elasticity and dispersion of plane waves,”
*International Journal of Engineering Science*, vol. 10, no. 5, pp. 425–435, 1972. View at Google Scholar · View at Scopus - O. Frank, G. Tsoukleri, J. Parthenios et al., “Compression behavior of single-layer graphenes,”
*ACS Nano*, vol. 4, no. 6, pp. 3131–3138, 2010. View at Publisher · View at Google Scholar · View at Scopus - M. Neek-Amal and F. M. Peeters, “Buckled circular monolayer graphene: a graphene nano-bowl,”
*Journal of Physics Condensed Matter*, vol. 23, no. 4, Article ID 045002, 2011. View at Publisher · View at Google Scholar · View at Scopus - T. Natsuki, J. Shi, and Q. Ni, “Buckling instability of circular double-layered graphene sheets,”
*Journal of Physics Condensed Matter*, vol. 24, no. 13, Article ID 135004, 2012. View at Publisher · View at Google Scholar · View at Scopus - A. Farajpour, M. Mohammadi, A. R. Shahidi, and M. Mahzoon, “Axisymmetric buckling of the circular graphene sheets with the nonlocal continuum plate model,”
*Physica E*, vol. 43, no. 10, pp. 1820–1825, 2011. View at Publisher · View at Google Scholar · View at Scopus