# A Molecular Dynamics Study on Wrinkles in Graphene with Simply Supported Boundary under In-Plane Shear

Authors: Jianzhang Huang; Qiang Han
Journal: Journal of Nanomaterials
 (2017)

---

## Abstract

The formation and evolution mechanisms of wrinkling in a rectangular single layer graphene sheet (SLGS) with simply supported boundary subjected to in-plane shear displacements are investigated using molecular dynamics simulations. Through investigating the out-of-plane displacements of the key point atom, we clarify the wrinkling growth and propagation process. Our results show that the boundary condition plays important roles in the wrinkling deformation. And the dependence of wrinkling parameters on the applied shear displacements is captured. Based on the elasticity theory, the formation mechanism of graphene wrinkling is revealed from the viewpoint of elastic energy. The effects of aspect ratio of graphene, temperature, and loading velocity on graphene wrinkling parameters and patterns are also investigated.


---

## Body

## 1. Introduction

Graphene, a research hotspot of nanoengineering, is an important constituent part of nanoelectromechanical systems due to its remarkable mechanical, thermal, and electrical properties [1]. Its fantastic physical properties make graphene an ideal candidate component for nanodevices, such as field-effect transistor, nanosensors, and supercapacitor [2–5]. As the real 2-dimensional (2D) material with atomic thickness, graphene is found to wrinkle easily because of its relatively low bending rigidity [6, 7]. Wrinkles can be induced on graphene by several reasons such as cooling-heating treatment, surface morphology growth of substrate, surface structural of substrate, topological defects of graphene, mechanical strain, and the boundary conditions [8–12]. Wrinkling, a kind of structural instability, certainly affects the graphenes remarkable physical property greatly, that is, mechanical, electrical, and thermal properties [13–16]. Dong et al. reported that the wrinkled graphene significantly increased the rate capability and cycling ability of MoO3C-reduced graphene oxide hybrid nanobelts, because the wrinkling morphology of graphene provides continuous electron pathways and stress buffering effects [17]. It is found that the wrinkle size and its top curvature improve the field enhancement factor of the graphene and decrease the electron affinities and ionization potentials [15]. Qin et al. investigated the mechanical properties of wrinkled graphene generated by topological defects and found that the interlayer shear mechanical properties are significantly enhanced due to the geometrical locking effect [18]. Meanwhile, it is possible to control the configuration and morphology of graphene wrinkles to satisfy the demand of nanoengineering applications [9, 19–21]. Therefore, understanding the deformation mechanism and controlling method of graphene wrinkling have fatal research significance.The edge boundary condition greatly influences the wrinkling morphology in graphene [10, 11]. But to the best of our knowledge, up to now it is not clear of the effects of simply supported boundary condition on graphene wrinkling. To fill the knowledge gap with regard to the simply supported boundary on graphene wrinkling, we investigate the wrinkling mechanism of shear-loaded zigzag rectangular SLGS with simply supported boundary condition investigated based on molecular dynamics (MD) simulations in this paper. The wrinkling deformation mechanism is elaborated through chasing the development of the atomic out-of-plane displacement development of the key atom and wrinkling parameters. The formation mechanism of wrinkling is also discussed by elasticity theory. We also investigate the effects of aspect ratio, temperature, and loading velocity on graphene wrinkling in this paper.


## 2. Simulation Model and Method

Figure1 shows a zigzag rectangular graphene sheet model (12.053 × 4.052 nm2) with specific edge boundary conditions subjected to in-plane shear. The shear displacement loading is applied on its top layer of atoms, and all the degrees of freedom of the bottom layer atoms are fixed. The left and right edges of graphene are simply supported boundary, of which only the degree of freedom in z-direction is fixed. The boundary condition can be expressed by(1)w=0,ifx=-lorl,where w is the out-of-plane displacement and the length of graphene is 2l. According to the elasticity theory, the wrinkling stability problem of graphene can be represented by a governing equation. The critical wrinkling strains and wrinkling modes can be obtained through the solution of the governing equation. The boundary condition affects the stress distribution in the graphene and has decisive influence on the wrinkling stability and modes. Meanwhile, the boundary condition plays a crucial role in the solution of governing equation of wrinkling problem. Therefore, the wrinkling stability is very sensitive to the boundary condition.Figure 1
Graphene model with given edge boundary condition.A key pointA (0.1418 nm; 0.0614 nm) and section line I (y=0.0614nm) are set upon the graphene surface as seen in Figure 1. The C-C bond interaction is described using AIREBO potential [22] and the length of C-C bond is 0.141 nm with the cut-off distance of 0.2 nm. We employ the MD simulation using the large-scale atomic/molecular massively parallel simulator (LAMMPS) with a time step of 1 fs for integration. The temperature is maintained at 0.01 K using the Nose-Hoover thermostat. After the equilibration, the shear strain is loaded in the way of applying shear displacements d on the loading edge 0.01 nm at every 10000 steps and repeated 200 times.


## 3. Results and Discussion

Self-similar wrinkles with distinct direction angle are generated on the graphene surface by the external mechanical strain. We define the global wrinkling parametersηχ¯ to describe the wrinkling characteristics. This ηχ¯ is used to reflect the mean value of the wrinkling characteristics value defined as the following expression:(2)ηχ¯=1n∑i=1nηχi.As shown in Figure2, ηχ represents an arbitrary specific wrinkling parameter, such as wrinkling amplitude ηA, wavelength ηλ, the ratio of wrinkling amplitude to wavelength ηA/ηλ, and the wrinkling direction angle ηα as well. The wrinkling direction angle ηα is defined as the angle between wrinkling direction and x-axis direction. The parameter ηA/ηλ may reflect the morphology of the graphene wrinkles. The wrinkling number n represents the total number of wrinkling crests and troughs. For example, n=7 (4, 3) indicates that there are 7 wrinkles including 4 wrinkling crests and 3 troughs.Figure 2
Wrinkling parameters of cross-line I.As shown in Figure3, based on the development of atomic out-of-plane displacement of the key point A, the wrinkling developing evolution can be divided into five stages, that is (I) prophase stage, (II) prometaphase stage, (III) metaphase stage, (IV) anaphase stage, and (V) telophase stage. Graphene structure is in a stable stage at first after the shear displacements load. As the d increases to 0.054 nm, wrinkles are formed on the graphene surface and wrinkling evolution begins. Prophase stage occurs in 0.054nm≤d≤0.106nm, during which graphene wrinkles rapidly and wrinkling numbers are 4 (2, 2) as shown in Figure 3. When d reaches 0.106 nm, a new wrinkling crest is formed from the right simply supported edge increasing the wrinkling numbers to 5 (3, 2) as shown in Figure 4(a) and evolution goes into prometaphase stage, which happens in 0.106nm<d≤0.353nm. After this, the evolution process is pushed into the metaphase stage (0.353nm<d≤0.782nm). In this stage, a wrinkling trough is generated from the right edge as shown in Figure 4(b) and n=6 (3, 3). The anaphase stage starts with d=1.151nm and terminates by d=1.151nm. During anaphase stage, wrinkling numbers rise to 7 (4, 3) owing to the formation of the new wrinkling crest from the right edge as shown in Figure 4(c). When d reaches 1.151 nm, the evolution process is driven into the last stage, telophase. During the previous four stages, the reciprocating change of key point atomic displacements means that wrinkles develop unstably. Unlike the previous four stages, in telophase stage, wrinkling angle decreases and wrinkling number keeps 7 (4, 3) as shown in Figure 4(d) and wrinkling grows stably till the graphene breaks, corresponding to d=1.316nm.Figure 3
Out-of-plane displacements of key pointA and wrinkling propagation.Figure 4
Wrinkling patterns at differentd (unit: Å).
(a)
d
=
0.21
n
m
(b)
d
=
0.44
n
m
(c)
d
=
0.88
n
m
(d)
d
=
1.24
n
mTable1 lists the computing wrinkling parameters results of cross-line I. It is seen that the wrinkling amplitude ηA¯ and ratio of amplitude to wavelength ηA/ηλ¯ grow with the increasing d, but the wrinkling wavelength ηλ¯ shows a reduction during wrinkling evolution process, indicating that the wrinkles are becoming slender. It is found that the wrinkling number n increases from 5 to 7 as the shear load goes on. Meanwhile, the amplitude and wavelength keep on progressing. Comparing with the opening researches, the wrinkling numbers in the graphene with simply supported edge boundary are more than the one with free boundary. As d increases, ηα¯ has a reduction and the maximum of out-of-plane displacements wE grows. But at the eve fracture of graphene (d=1.24nm), wE reaches 0.367 nm. The variations of the wrinkling parameters over d are shown in Figure 5. It is observed that two turning points occur on the curve owing to the generation of new wrinkles, respectively, which has influence on the mean value of wrinkling parameters. It shows that ηA¯,  ηA/ηλ¯, and wE have increasing trend over d, while their increasing rates slow down at larger d. The wrinkling wavelength ηλ¯ decreases rapidly at lower d and then slowly reduces as d increases. Due to the generation of the new wrinkle, the wrinkling angle ηα¯ exhibits a slight increment.Table 1
Wrinkling parameters overd.
d (nm)
η
A
¯ (nm)
η
λ
¯ (nm)
η
A
/
η
λ
¯
n (crest, trough)
η
α
¯ (∘)
w
E
0.21
0.142
2.409
0.059
5(3,2)
44.493
0.172
0.44
0.212
1.874
0.113
6(3,3)
42.954
0.233
0.88
0.278
1.554
0.179
7(4,3)
42.895
0.305
1.24
0.338
1.523
0.222
7(4,3)
42.097
0.367Figure 5
Curves of wrinkling parameters overd.
(a)
Amplitude
(b)
Wavelength
(c)
Wrinkling ratio
(d)
Direction angle
(e)
w
EThe graphene initially undergoes in-plane deformation once subjected to shear loading. As the loading displacement increases, owing to the finite thickness of graphene, the graphene displays out-of-plane deformation to release the strain energy and then forms the wrinkles to reach a low-energy configuration. Based on the elasticity theory [23], the strain energy per unit area of wrinkled graphene can be represented by(3)U¯=Us¯+Ub¯,where Us¯ and Ub¯ are stretching and bending terms of the energy, respectively. They can be expressed as(4)Ub¯=12Et3121-ν2π4A24λ4,Us¯=Et2γ24+Et2π2A2γ16H2,where γ=d/H and H and t are the width and thickness of graphene. A and λ represent the wrinkling amplitude and wavelength. E and ν are Young’s modulus and Poisson’s ratio of graphene.Use the relation of(5)A=21-νπλ.The strain energy per unit area U¯ can be expressed in terms of the wrinkling amplitude A as(6)U¯=Etγ21-νγt2121+νA2+π2A216H2.It is known, from (6), the energy increases as the wrinkling amplitude rises. After a new wrinkle is produced, a reduction of the wrinkling amplitude occurs, as seen in Figure 5(a). The reduction of wrinkling amplitude means that the energy of the new wrinkling configuration is lower than the ones before the new wrinkle occurs, which implies that a more stable wrinkling configuration is formed. As the loading continues, the energy keeps elevating and pushing the progress of the wrinkling. Until the occurrence of second new wrinkle, the energy decreases and forms a stable configuration. It can be observed that the formation of wrinkling is a low-energy process.


### 3.1. Aspect Ratio Effects

We define the aspect ratio as the ratio of length to width of the graphene sheet with a fixed width of 4.025 nm. The effect on graphene wrinkling of chosen aspect ratio in ranges of 1 to 13 is studied in this paper. The wrinkling parameters, patterns, and atoms displacements in cross-line under samed=0.88nm are shown in Table 2 and Figures 6 and 7. The parameter nηλ¯/L stands for the extent of wrinkling in graphene, where L represents the length of graphene sheet. It is observed that the wrinkling amplitude ηA¯ and ratio ηA/ηλ¯ decrease with the increasing aspect ratio. On the contrary, the wrinkling wavelength ηλ¯, wrinkling number n, and extent parameter nηλ¯/L increase as the aspect ratio increases. The wrinkles become stubby and flat with a wider coverage extent, as shown in Figures 6 and 7.Table 2
Wrinkling parameters in different aspect ratio(d=0.88nm).
Aspect ratio
η
A
¯ (nm)
η
λ
¯ (nm)
η
A
/
η
λ
¯
n (crest, trough)
n
η
λ
¯
/
L
wE
1
0.321
1.543
0.208
2(1,1)
0.762
0.328
2
0.305
1.576
0.193
4(2,2)
0.783
0.316
3
0.283
1.570
0.181
7(4,3)
0.904
0.308
5
0.292
1.644
0.178
11(5,6)
0.893
0.323
8
0.286
1.661
0.172
18(9,9)
0.922
0.323
13
0.278
1.676
0.166
30(15,15)
0.955
0.317Figure 6
Wrinkling patterns in different aspect ratio(d=0.88nm), unit: Å. (a) Aspect ratio 1. (b) Aspect ratio 2. (c) Aspect ratio 3. (d) Aspect ratio 5. (e) Aspect ratio 8. (f) Aspect ratio 13.
(a)
(b)
(c)
(d)
(e)
(f)Figure 7
Atoms displacements in cross-line I with different aspect ratio(d=0.88nm).


### 3.2. Temperature Effects

We investigate the temperature effects on wrinkling parameters and patterns of rectangular graphene with aspect ratio of 3 under temperature ranging from 0.01 K to 600 K. The wrinkling parameters and patterns atd=0.47nm are listed in Table 3 and Figure 8. As seen in Table 3, it shows that temperature has great influence on graphene wrinkling process and wrinkling parameters. The existence of temperature intensifies the thermal vibration of graphene atoms, which significantly influences the wrinkling deformation process. As a consequence, severe variation of the wrinkling parameters occurs as the temperature rises. Generally, it displays an increasing trend of ηA¯ and ηA/ηλ¯ as temperature increases. ηλ¯ and nηλ¯/L are enlarged with the increasing temperature, when temperature is less than 400 K, after which ηλ¯ and nηλ¯/L reduce as temperature increases. Obvious discrepancy between wrinkling patterns under different temperatures is observed from Figure 8, indicating the significant temperature effect on wrinkling pattern, which affects the wrinkling parameters.Table 3
Wrinkling parameters under different temperatures(d=0.47nm).
Temperature (K)
η
A
¯ (nm)
η
λ
¯ (nm)
η
A
/
η
λ
¯
n (crest, trough)
nηλ¯/L
w
E
0.01
0.219
1.859
0.118
6(3,3)
0.918
0.276
100
0.226
1.823
0.124
6(3,3)
0.900
0.242
100 to 0.01
0.217
1.869
0.116
6(3,3)
0.930
0.237
200
0.218
1.836
0.119
6(3,3)
0.906
0.250
200 to 0.01
0.217
1.867
0.117
6(3,3)
0.929
0.237
300
0.234
1.866
0.126
6(3,3)
0.921
0.275
300 to 0.01
0.218
1.875
0.116
6(3,3)
0.934
0.237
400
0.237
1.918
0.124
6(3,3)
0.947
0.271
400 to 0.01
0.217
1.871
0.116
6(3,3)
0.931
0.237
500
0.237
1.881
0.126
6(3,3)
0.928
0.281
500 to 0.01
0.217
1.871
0.116
6(3,3)
0.931
0.237
600
0.236
1.860
0.126
6(3,3)
0.918
0.313
600 to 0.01
0.217
1.868
0.117
6(3,3)
0.930
0.237Figure 8
Wrinkling patterns under different temperatures(d=0.47nm), unit: Å.
(a)
300 K
(b)
600 K
(c)
300 K anneals to 0.01 K
(d)
600 K anneals to 0.01 KWe perform an annealing treatment simulation to wrinkling graphene and investigate the effect of annealing on graphene wrinkling. As seen in Figure9, after graphene wrinkles, the temperature decreases to 0.01 K gradually. After the temperature is maintained in 0.01 K, the wrinkling parameters and configuration are obtained and shown in Table 3 and Figure 8. The wrinkling parameters under different temperature situation become close after annealing treatment, indicating that the annealing treatment eliminates the effect of temperature on wrinkling parameters. But, as for wrinkling configuration, the wrinkling pattern is still the same as the one before annealing treatment. It means that the annealing treatment cannot recover the wrinkling pattern to low temperature.Figure 9
Annealing treatment simulation, the system temperature decreases to 0.01 K gradually.Based on the elasticity theory, the graphene structure wrinkles when it is subjected to critical strain. With the existence of temperature, the atoms of graphene exhibit thermal motion, which produces nonnegligible atomic out-of-plane displacement compared to the thickness of graphene. Therefore ripples are formed upon the graphene [6]. The out-of-plane of the wrinkling bifurcation point significantly influences the wrinkling configuration.


### 3.3. Loading Velocity Effects

The wrinkling of graphene sheet with aspect ratio of 3 under different loading velocity is investigated in this section. The loading velocity ranges from6×10-5 to 2.2×10-4nm/ps and the other simulation settings remain unchanged. The wrinkling parameters and patterns of graphene wrinkles under d=0.88nm are shown in Table 4 and Figure 10. It shows that as the loading velocity increases the wrinkling parameters change slightly, indicating that wrinkling parameters have low dependency upon the loading velocity. Based on the results, in all loading velocity cases graphene generates the same number wrinkles, but the wrinkling patterns are different as seen in Table 4 and Figure 10. It illustrates that loading velocity has impact on wrinkling patterns but slightly affects wrinkling parameters. The loading velocity not only affects the mechanical property performance of materials but also affects the deformation response of the materials [24, 25]. Under different loading velocity, the relaxation time for the rearrangement of carbon atoms and bonds is different. On the eve of the wrinkling of graphene, the rearrangement situation of the carbon atoms has great influence on the wrinkling configuration after graphene wrinkles.Table 4
Wrinkling parameters under different loading velocity(d=0.88nm).
Velocity(10-5nm/ps)
η
A
¯ (nm)
η
λ
¯ (nm)
η
A
/
η
λ
¯
n (crest, trough)
n
η
λ
¯
/
L
w
E
6
0.280
1.580
0.178
7(3,4)
0.916
0.305
8
0.281
1.573
0.179
7(3,4)
0.912
0.306
10
0.283
1.570
0.181
7(4,3)
0.910
0.308
12
0.279
1.582
0.176
7(4,3)
0.917
0.303
14
0.280
1.575
0.178
7(3,4)
0.913
0.305
18
0.280
1.576
0.178
7(4,3)
0.914
0.304
22
0.280
1.569
0.179
7(3,4)
0.910
0.305Figure 10
Wrinkling patterns under different loading velocity(d=0.88nm), unit: Å. (a) 8×10-5nm/ps. (b) 1.8×10-4nm/ps.
(a)
(b)


## 3.1. Aspect Ratio Effects

We define the aspect ratio as the ratio of length to width of the graphene sheet with a fixed width of 4.025 nm. The effect on graphene wrinkling of chosen aspect ratio in ranges of 1 to 13 is studied in this paper. The wrinkling parameters, patterns, and atoms displacements in cross-line under samed=0.88nm are shown in Table 2 and Figures 6 and 7. The parameter nηλ¯/L stands for the extent of wrinkling in graphene, where L represents the length of graphene sheet. It is observed that the wrinkling amplitude ηA¯ and ratio ηA/ηλ¯ decrease with the increasing aspect ratio. On the contrary, the wrinkling wavelength ηλ¯, wrinkling number n, and extent parameter nηλ¯/L increase as the aspect ratio increases. The wrinkles become stubby and flat with a wider coverage extent, as shown in Figures 6 and 7.Table 2
Wrinkling parameters in different aspect ratio(d=0.88nm).
Aspect ratio
η
A
¯ (nm)
η
λ
¯ (nm)
η
A
/
η
λ
¯
n (crest, trough)
n
η
λ
¯
/
L
wE
1
0.321
1.543
0.208
2(1,1)
0.762
0.328
2
0.305
1.576
0.193
4(2,2)
0.783
0.316
3
0.283
1.570
0.181
7(4,3)
0.904
0.308
5
0.292
1.644
0.178
11(5,6)
0.893
0.323
8
0.286
1.661
0.172
18(9,9)
0.922
0.323
13
0.278
1.676
0.166
30(15,15)
0.955
0.317Figure 6
Wrinkling patterns in different aspect ratio(d=0.88nm), unit: Å. (a) Aspect ratio 1. (b) Aspect ratio 2. (c) Aspect ratio 3. (d) Aspect ratio 5. (e) Aspect ratio 8. (f) Aspect ratio 13.
(a)
(b)
(c)
(d)
(e)
(f)Figure 7
Atoms displacements in cross-line I with different aspect ratio(d=0.88nm).


## 3.2. Temperature Effects

We investigate the temperature effects on wrinkling parameters and patterns of rectangular graphene with aspect ratio of 3 under temperature ranging from 0.01 K to 600 K. The wrinkling parameters and patterns atd=0.47nm are listed in Table 3 and Figure 8. As seen in Table 3, it shows that temperature has great influence on graphene wrinkling process and wrinkling parameters. The existence of temperature intensifies the thermal vibration of graphene atoms, which significantly influences the wrinkling deformation process. As a consequence, severe variation of the wrinkling parameters occurs as the temperature rises. Generally, it displays an increasing trend of ηA¯ and ηA/ηλ¯ as temperature increases. ηλ¯ and nηλ¯/L are enlarged with the increasing temperature, when temperature is less than 400 K, after which ηλ¯ and nηλ¯/L reduce as temperature increases. Obvious discrepancy between wrinkling patterns under different temperatures is observed from Figure 8, indicating the significant temperature effect on wrinkling pattern, which affects the wrinkling parameters.Table 3
Wrinkling parameters under different temperatures(d=0.47nm).
Temperature (K)
η
A
¯ (nm)
η
λ
¯ (nm)
η
A
/
η
λ
¯
n (crest, trough)
nηλ¯/L
w
E
0.01
0.219
1.859
0.118
6(3,3)
0.918
0.276
100
0.226
1.823
0.124
6(3,3)
0.900
0.242
100 to 0.01
0.217
1.869
0.116
6(3,3)
0.930
0.237
200
0.218
1.836
0.119
6(3,3)
0.906
0.250
200 to 0.01
0.217
1.867
0.117
6(3,3)
0.929
0.237
300
0.234
1.866
0.126
6(3,3)
0.921
0.275
300 to 0.01
0.218
1.875
0.116
6(3,3)
0.934
0.237
400
0.237
1.918
0.124
6(3,3)
0.947
0.271
400 to 0.01
0.217
1.871
0.116
6(3,3)
0.931
0.237
500
0.237
1.881
0.126
6(3,3)
0.928
0.281
500 to 0.01
0.217
1.871
0.116
6(3,3)
0.931
0.237
600
0.236
1.860
0.126
6(3,3)
0.918
0.313
600 to 0.01
0.217
1.868
0.117
6(3,3)
0.930
0.237Figure 8
Wrinkling patterns under different temperatures(d=0.47nm), unit: Å.
(a)
300 K
(b)
600 K
(c)
300 K anneals to 0.01 K
(d)
600 K anneals to 0.01 KWe perform an annealing treatment simulation to wrinkling graphene and investigate the effect of annealing on graphene wrinkling. As seen in Figure9, after graphene wrinkles, the temperature decreases to 0.01 K gradually. After the temperature is maintained in 0.01 K, the wrinkling parameters and configuration are obtained and shown in Table 3 and Figure 8. The wrinkling parameters under different temperature situation become close after annealing treatment, indicating that the annealing treatment eliminates the effect of temperature on wrinkling parameters. But, as for wrinkling configuration, the wrinkling pattern is still the same as the one before annealing treatment. It means that the annealing treatment cannot recover the wrinkling pattern to low temperature.Figure 9
Annealing treatment simulation, the system temperature decreases to 0.01 K gradually.Based on the elasticity theory, the graphene structure wrinkles when it is subjected to critical strain. With the existence of temperature, the atoms of graphene exhibit thermal motion, which produces nonnegligible atomic out-of-plane displacement compared to the thickness of graphene. Therefore ripples are formed upon the graphene [6]. The out-of-plane of the wrinkling bifurcation point significantly influences the wrinkling configuration.


## 3.3. Loading Velocity Effects

The wrinkling of graphene sheet with aspect ratio of 3 under different loading velocity is investigated in this section. The loading velocity ranges from6×10-5 to 2.2×10-4nm/ps and the other simulation settings remain unchanged. The wrinkling parameters and patterns of graphene wrinkles under d=0.88nm are shown in Table 4 and Figure 10. It shows that as the loading velocity increases the wrinkling parameters change slightly, indicating that wrinkling parameters have low dependency upon the loading velocity. Based on the results, in all loading velocity cases graphene generates the same number wrinkles, but the wrinkling patterns are different as seen in Table 4 and Figure 10. It illustrates that loading velocity has impact on wrinkling patterns but slightly affects wrinkling parameters. The loading velocity not only affects the mechanical property performance of materials but also affects the deformation response of the materials [24, 25]. Under different loading velocity, the relaxation time for the rearrangement of carbon atoms and bonds is different. On the eve of the wrinkling of graphene, the rearrangement situation of the carbon atoms has great influence on the wrinkling configuration after graphene wrinkles.Table 4
Wrinkling parameters under different loading velocity(d=0.88nm).
Velocity(10-5nm/ps)
η
A
¯ (nm)
η
λ
¯ (nm)
η
A
/
η
λ
¯
n (crest, trough)
n
η
λ
¯
/
L
w
E
6
0.280
1.580
0.178
7(3,4)
0.916
0.305
8
0.281
1.573
0.179
7(3,4)
0.912
0.306
10
0.283
1.570
0.181
7(4,3)
0.910
0.308
12
0.279
1.582
0.176
7(4,3)
0.917
0.303
14
0.280
1.575
0.178
7(3,4)
0.913
0.305
18
0.280
1.576
0.178
7(4,3)
0.914
0.304
22
0.280
1.569
0.179
7(3,4)
0.910
0.305Figure 10
Wrinkling patterns under different loading velocity(d=0.88nm), unit: Å. (a) 8×10-5nm/ps. (b) 1.8×10-4nm/ps.
(a)
(b)


## 4. Conclusions

In summary, in the present paper the wrinkling of shear-loaded rectangular SLGS with simply supported boundary condition is investigated through MD simulations. It reveals that, due to the graphene shape and edge boundary condition, the wrinkles are produced on the graphene surface, and the new wrinkles are generated from the right edge driving the wrinkles to propagate to the left. Generation of new wrinkles changes the direction angles of the wrinkles group. The wrinkling revolution process is divided into five stages: (I) prophase stage, (II) prometaphase stage, (III) metaphase stage, (IV) anaphase stage, and (V) telophase stage. The wrinkling amplitude, maximum of out-of-plane displacements, and ratio of amplitude and wavelength rise with the increasing shear displacements, but the wrinkling wavelength reduces in general. The wrinkling formation is the energy-releasing process during the deformation of graphene under shear loading. The aspect ratio of graphene and temperature greatly affects the graphene wrinkling. The larger aspect ratio makes the wrinkles flatter and wider extent coverage. Wrinkling deformation of graphene is significantly influenced by temperature due to the intensive atomic thermal vibration. And it is found that the annealing process could eliminate the temperature effect on wrinkling parameter but could not recover graphene wrinkling configuration to low temperature. The loading velocity has little effect on wrinkling parameters but affects the wrinkling patterns. Because of the rearrangement of carbon atoms, the loading velocity has influence on the wrinkling configuration when graphene wrinkles.



---

*Source: 1326790-2017-03-07.xml*