Next Article in Journal
Investigation of Adhesive’s Material in Hermetic MEMS Package for Interfacial Crack between the Silver Epoxy and the Metal Lid during the Precondition Test
Next Article in Special Issue
Estimate of Coffin–Manson Curve Shift for the Porous Alloy AlSi9Cu3 Based on Numerical Simulations of a Porous Material Carried Out by Using the Taguchi Array
Previous Article in Journal
Electrodeposition and Corrosion Properties of Nickel–Graphene Oxide Composite Coatings
Previous Article in Special Issue
Investigation of the Relationship between Morphology and Thermal Conductivity of Powder Metallurgically Prepared Aluminium Foams
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mechanical Properties and Deformation Mechanisms of Graphene Foams with Bi-Modal Sheet Thickness by Coarse-Grained Molecular Dynamics Simulations

1
School of Mechanics and Civil Engineering, China University of Mining and Technology, Beijing 100083, China
2
LNM, Institute of Mechanics, Chinese Academy of Sciences, Beijing 100190, China
3
School of Engineering Science, University of Chinese Academy of Sciences, Beijing 100049, China
*
Author to whom correspondence should be addressed.
Materials 2021, 14(19), 5622; https://doi.org/10.3390/ma14195622
Submission received: 23 August 2021 / Revised: 23 September 2021 / Accepted: 24 September 2021 / Published: 27 September 2021
(This article belongs to the Special Issue New Advances in Characterization of Cellular Materials)

Abstract

:
Graphene foams (GrFs) have been widely used as structural and/or functional materials in many practical applications. They are always assembled by thin and thick graphene sheets with multiple thicknesses; however, the effect of this basic structural feature has been poorly understood by existing theoretical models. Here, we propose a coarse-grained bi-modal GrF model composed of a mixture of 1-layer flexible and 8-layer stiff sheets to study the mechanical properties and deformation mechanisms based on the mesoscopic model of graphene sheets (Model. Simul. Mater. Sci. Eng. 2011, 19, 54003). It is found that the modulus increases almost linearly with an increased proportion of 8-layer sheets, which is well explained by the mixture rule; the strength decreases first and reaches the minimum value at a critical proportion of stiff sheets ~30%, which is well explained by the analysis of structural connectivity and deformation energy of bi-modal GrFs. Furthermore, high-stress regions are mainly dispersed in thick sheets, while large-strain areas mainly locate in thin ones. Both of them have a highly uneven distribution in GrFs due to the intrinsic heterogeneity in both structures and the mechanical properties of sheets. Moreover, the elastic recovery ability of GrFs can be enhanced by adding more thick sheets. These results should be helpful for us to understand and further guide the design of advanced GrF-based materials.

1. Introduction

Graphene is the thinnest-known sheet, composed of a single layer of carbon atoms, which has been used as nanoscale building blocks to fabricate a good deal of macroscale bulk materials, such as one-dimensional graphene fibers [1,2,3,4], two-dimensional graphene papers [5] and three-dimensional (3D) graphene foams (GrFs) [6,7,8,9,10,11,12,13,14,15,16]. Among them, 3D GrFs have attracted much attention in recent years from material scientists to engineers due to their combined properties of excellent electrochemical stability [17], superior electrical conductivity [18], high energy absorption [19] and broad application prospects in advanced materials [20], stretchable electronic devices [21] and energy-storage components [22].
It is shown that the mechanical and physical properties of graphene sheets, as well as their macro-assembles of GrFs are greatly influenced by thickness of sheets. For example, for a single graphene sheet, the bending stiffness is proportional to the third power of thickness if there is no interlayer slippage, according to the theory of plates [23,24]. If the effect of interlayer slippage cannot be neglected, the stiffness of multilayer graphene sheets is still highly dependent on thickness based on the experimental study [25]. Roughly speaking, the thicker the graphene sheet, the greater the stiffness. The electrical conductivity of a single graphene sheet increases with an increased thickness according to the first principle calculation [26]. For their macro-assemblies of GrFs, Wang et al. [27,28] found that both elastic properties and Poisson’s ratio of GrFs increased with an increased thickness of sheets using coarse-grained molecular dynamics simulations (CGMD). Chen et al. [21] experimentally found that the thickness-dependent conductivity of GrFs reached the maximum value of ~10 S/cm at an optimal number of graphene layers of around 5, which was reproduced later by Liu et al. [26] using a combination of the first principle calculations and CGMD simulations.
Due to complicated preparation processes involving redox reactions and the growth/aggregation of graphene sheets by chemical vapor deposition (CVD) [10,29] or the self-assembly method [8,30,31], GrFs are always composed of sheets having multiple thicknesses rather than a single thickness. For example, Nieto et al. [10] made a freestanding GrF consisting of graphene sheets with ~3–10 layers, which was grown onto a nickel foam structure by a CVD method. Zhao et al. [32] made a GrF composed of sheets with 1–4 graphene layers by CVD and electrochemical strategies. More experimental reports about varied thickness of graphene sheets in GrFs can be found in [8,21,30,33,34].
However, to our best knowledge, all numerical models of GrFs published in references, no matter the full-atomic models [16,35,36,37,38], coarse-grained models [28,39,40,41,42] or finite element one [43], have a fairly idealized assumption that the thickness of constituent graphene sheets is the same throughout GrF systems. For example, all samples using full-atomic models [16,35,36,37,38] are composed of thinnest 1-layer graphene sheets, while other coarse-grained models [28,39,40,41,42] or finite element one [43] are composed of sheets with 1-layer or multi-layer single thickness. Although the GrFs composed of graphene sheets with multiple thicknesses have been widely observed in practical systems as mentioned above, the mechanics of this kind of GrFs have not been studied up to now, and the role of graphene sheets of different thicknesses in the system, as well as the synergy between them, remains elusive.
In this paper, a so-called bi-modal GrF model composed of 1-layer flexible sheets and 8-layer stiff sheets are proposed to systematically study these essential issues based on the coarse-grained model of graphene sheets with varied thickness proposed by Cranford et al. [44], as given in Section 2. Section 3 gives the tensile stress–strain curves of GrFs with varied proportions of stiff sheets, and the modulus increases almost linearly with an increased proportion of stiff sheets; there exists a critical proportion of stiff sheets ~30%, at which the strength reaches the minimum value. This critical phenomenon is explained by the analysis of structural connectivity by stiff sheets and the deformation energy of bi-modal GrFs. Furthermore, we find that high-stress states always emerge in stiff sheets, while flexible sheets experience larger deformation. The conclusion is given at the end.

2. Numerical Model of Bi-Modal GrFs and Methodology

We adopted the coarse-grained graphene model developed by Cranford et al. [44] and proven effective and computationally efficient in a series of studies including the mechanical deformations of a single graphene nanoribbon [44] and graphene macro-assemblies [28,45,46,47]. In this scheme, a square graphene sheet with the side length of 2.5 nm containing 264 carbon atoms Figure 1a-i is reduced to be a so-called coarse grain. So, a larger square graphene sheet with a side length of 75 nm can be represented by only 900 coarse grains connected by a set of linear springs and angle ones as shown in Figure 1a-ii,a-iii. The total energy of a graphene sheet is calculated by E total = i = 1 N B E B i + i = 1 N φ E φ i + i = 1 N θ E θ i + i = 1 N LJ E LJ i , where NB, Nφ, Nθ and NLJ are the number of bonds (linear springs), in-plane angles, out-of-plane angles and bead pairs, respectively; E B i , E φ i , E θ i and E LJ i are the bond energy, in-plane angle energy, out-of-plane angle energy and van der Waals energy, respectively, and are calculated by Equations (1)–(4) in simulations; the parameters of k B , k φ and k θ are the stiffness of the linear spring for bond stretching, the angle spring for in-plane shearing and the out-of-plane bending of sheet, respectively; the equilibrium distance of two bonded beads r 0 is 25 Å; the in-plane equilibrium angle φ 0 between two bonds connecting three neighbor beads is 90°; the out-of-plane equilibrium angle θ 0 between two bonds connecting three neighbor beads is 180°; ε is the depth of the potential well, and σ represents the zero-energy distance between two beads in the van der Waals potential function. r is the current distance between two beads. In practical GrF materials, a constituent graphene sheet contains about 1–10 graphene layers [10,21,48]. Here, a GrF composed of a mixture of 1-layer flexible sheets and 8-layer stiff sheets is used to study the mechanics of bi-modal GrFs. All mechanical parameters for the coarse-grained 1-layer and 8-layer sheets in our simulations are listed in Table 1, which was obtained based on the equivalent energy principle by Cranford et al. [44].
E B = k B ( r r 0 ) 2 / 2
E φ = k φ ( φ φ 0 ) 2 / 2
E θ = k θ ( θ θ 0 ) 2 / 2
E LJ = 4 ε ( ( σ / r ) 12 ( σ / r ) 6 )
Using the numerical synthesizing method described in detail at the end (Appendix A: Preparation of Numerical Model) of this paper, we obtain a series of well-equilibrated GrFs with a tunable proportion of 8-layer sheets, i.e., v = 0, 30%, 50%, 70% and 100%. Figure 1b shows the initial state of the GrF with the proportion of 8-layer sheet v = 50%, in which graphene sheets are randomly distributed in the foam system similar to that observed in a SEM experiment [10]. In order to mimic the strong connection between neighbor sheets due to strong physical crosslinks or chemical functional groups [10,49] in practical GrF systems, we adopt a crosslink model which has been widely used to study the large-deformation behavior and fracture mode of both carbon nanotube buckypapers [49,50] and graphene foams [27,46]. As shown in Figure 1b, a certain amount of crosslinks (green color) characterized by Equation (5) are added between neighbor sheets in all GrF samples to enhance inter-sheet connections. k C is the spring constant, and r is the current distance between two beads in different flakes with a equilibrium distance r 0 = 25 Å. Furthermore, we assume that the bonds in sheets as well as crosslinks between neighbor sheets would break when the local tensile strain surpasses the critical value of 12% according to the experimental and theoretical researches [51,52]. In addition, the equilibrium density of all GrF systems ranges from 46 to 303 mg/cm3, which is comparable to that obtained in experimental systems in the range of 1–1400 mg/cm3 [6,38,48].
E C = k C ( r r 0 ) 2 / 2
The four typical microstructures in bi-modal GrFs, i.e., edge–edge (E–E), edge–surface (E–S), point–surface (P–S) and surface–surface (S–S) for both flexible sheets and stiff sheets are shown in Figure 1c. The proportions of S–S contacts with a larger contact area slightly decrease while the proportions of E–S and E–E and P–S contacts with a smaller contact area slightly increase as an increased fraction of stiff sheets in GrFs as shown in Figure 1d. Accordingly, the areal density (the contact area divided by the system volume) decreases with the increased proportion of stiff sheets.
All systems are equilibrated by energy minimization and then relaxation at 300 K and 1 atmosphere under an isothermal–isobaric (NPT) ensemble with the criterion that the total energy fluctuation converges to less than 1%. In equilibration, the periodic boundary conditions are imposed in all three directions of the simulated samples. In the loading process, samples are uniaxially tensioned in the x-axis direction with a tensile strain rate 4.3 × 106  s 1 with a zero-pressure barostat in the other two unloading directions at room temperature. A time step of 10 fs is used. All simulations are performed in the large-scale atomic/molecular massively parallel simulator (LAMMPS) [53], and the results are visualized based on the Open Visualization Tool (OVITO) [54].

3. Results

3.1. Stress–Strain Curves in Uniaxial Tension

Figure 2a shows the uniaxial tensile stress–strain curves of the GrFs with varied proportions of 8-layer stiff sheets v = 0, 30%, 50%, 70% and 100%. The stress–strain curve of the GrF composed entirely of 1-layer sheets, i.e., v = 0, clearly shows four stages, as signified in the black curve: the linearly elastic stage, yielding stage, obvious hardening stage and the final fracture stage. As shown in Figure 2b-i, the snapshots in the front viewpoint at the tensile strain of 0, 0.25, 0.5 and 0.75 show that the system contracts obviously perpendicular to the tensile direction due to Poisson’s effect and breaks locally at the larger tensile strain of 0.5 and 0.75 as signified by yellow circles. For the GrFs composed of the mixture of 1-layer sheets and 8-layer ones at varied proportion v = 30%, 50% and 70% of 8-layer sheets, the stress generally increases with an increased v. Choosing the system with v = 50% as an example, in a wide range of tensile strain, the tensile stress is higher than that of the system containing a small amount of stiff sheets with v = 30%; the snapshots of the system at the tensile strain of 0, 0.25, 0.5 and 0.75 in Figure 2b-ii show that the system shrinks very little perpendicular to the tensile direction compared to the GrF composed entirely of 1-layer sheets. For the GrF composed entirely of 8-layer sheets, i.e., v = 100%, the tensile stress is much larger than that of all systems at any tensile strain; moreover, the system breaks seriously even at a smaller tensile strain of 0.25 as highlighted by yellow circles in Figure 2b-iii, and the system shrinks very little as the GrF with v = 50%. Qualitatively, it is seen that the fraction of stiff sheets greatly influences the macroscopic mechanical responses of the bi-modal GrFs.
Furthermore, the quantitative dependency between the modulus/strength of GrFs and the proportion of stiff sheets is numerically obtained and shown in Figure 2c. For each v, we conduct uniaxial tension simulations using five numerical samples to obtain the mean and variance of the modulus and strength of the system. The modulus increases almost linearly with the increased v because the stretching/shearing/bending modulus of 8-layer stiff sheets is much larger than that of 1-layer flexible sheets as given in Table 1, which is qualitatively consistent with the mixture rule [55]. Interestingly, we find that there exists a critical proportion of v~30%, at which the strength of the bi-modal GrFs reaches the minimum value, i.e., the strength of the GrF decreases first and reaches the minimum value at the critical v = ~30%, then increases monotonically with an increased v. In order to explain this critical phenomenon, we study the connectivity of stiff sheets in GrFs and find that the connectivity transforms from 0 to 1 as v > 30% as shown in Figure 2d, i.e., a complete bearing path by composed of entirely stiff sheets forms at this critical proportion, and the strength of the system would be instead determined by mechanical properties of stiff sheets. Comparing the GrFs with v = 0 and 30%, because flexible sheets have better deformation coordination ability, the number of broken bonds in the GrF with v = 0 is always smaller than that in the GrF with v = 30% as shown in Figure 2e, which is responsible for the larger strength of the GrF with v = 0. In other words, a small number of stiff sheets would deteriorate the bi-modal GrF because of their poor deformation coordination ability. For the system with v > 30%, the deformation of the system would be dominated by stiff sheets, because of the intrinsic high strength of stiff graphene sheets, the strength of the system increases monotonically with an increased proportion v.

3.2. The Distribution and Evolution of Deformation Energy

To further understand the deformation of bi-modal GrFs, we calculate the total stretching energy E B total = i = 1 N B E B i + j = 1 N C E C j , in-plane shearing energy E φ total = i = 1 N φ E φ i , out-of-plane bending energy E θ total = i = 1 N θ E θ i and van der Waals energy E LJ total = i = 1 N LJ E LJ i for the GrFs with varied v, where E B i , E φ i , E θ i , E LJ i and E C j are calculated by the Equations (1)–(5), respectively. Figure 3a shows that the total in-plane shearing energy of GrFs increases slightly with an increased tensile strain and saturates at a critical tensile strain of ~0.4; in particular, it is almost independent on v. This is consistent with our previous finding [27] that the shearing deformation of sheets is not the main deformation mode in GrFs and that the shearing energy of GrFs always keeps a small value regardless of the characteristic of the constituent graphene sheets. For the two other kinds of deformation modes, the stretching and bending of sheets, as shown in Figure 3b,c, on the one hand, the respective elastic energy of the two modes is larger than the corresponding shearing energy in Figure 3a, which indicates that the bending and stretching of sheets are still the main deformation modes in the bi-modal GrFs, as in in the GrFs composed of a single thickness of sheets [27]; on the other hand, both kinds of energy increase with an increased v at any given tensile strain; here, we note that the initial crosslink density, as well as the local connectivity of neighbor sheets, are almost the same for samples, so in the early stage of tension, because there is no bond breaking, the sheets in systems would experience almost similar stretching/shearing/bending deformations due to the same constrains imposed by crosslinks; however, due to the huge difference in stiffness and strength of 1-layer flexible sheets and 8-layer stiff ones, the immediate stress level, as well as the initial modulus of the GrF with more stiff sheets (larger v), show larger values as observed in Figure 2a,c. Moreover, the van der Waals adhesion energy Figure 3d exhibits a similar dependency on v as that of the stretching/bending energy; from the structural analysis in Figure 1d, we know that the GrF with more stiff sheets is apt to be a state with less inter-sheet contact area, which is responsible for the increase of the van der Waals energy as shown in Figure 3d. For the pure GrF or the GrF with 30%, the stretching/bending energy of the system increases much slightly with an increased tensile strain, as indicated by the slope bar; in sharp contrast, for the GrF containing more 8-layer sheets, e.g., v = 50%, 70% and 100%, the stretching/bending energy of the system increases obviously with an increased tensile strain, especially in the stage with the tensile strain <0.4; see the respective slope bar. This difference in the change of the stretching/bending energy at the critical value of v = ~30% indicates that the deformation mode of the system has shifted from being dominated by flexible sheets to being dominated by stiff sheets, which is consistent with the finding by the connectivity analysis as discussed in Figure 2d.

3.3. The Distribution and Evolution of Local Stress and Strain

In order to understand the role of flexible sheets and stiff ones in bi-modal GrFs in tension process, we study the distribution and evolution of local stress and strain in both 1-layer flexible and 8-layer stiff sheets. Figure 4a-i shows the distribution of local von-Mises stress in both flexible sheets and stiff ones in the GrF with v = 50%. In the initial state with the tensile strain ε = 0, the system is well equilibrated; then, with an increased tensile strain from 0 to 0.8, more and more high tensile stress regions as color coded by red emerge in the system. By separately displaying the stress distribution in flexible sheets and stiff ones in Figure 4a-ii,a-iii, we find that the high-stress regions mainly locate in stiff sheets, especially at the larger tensile strain of 0.4 and 0.8, which is mainly attributed to the poor deformation coordination ability of stiff sheets compared to the flexible ones. Furthermore, we count the number of the “high-stress beads”, defined as the magnitude of local von-Mises stress σvm being larger than the average value of the system, i.e., |σvm|> i = 1 N | σ v m i | / N , where N is the number of beads in the system. Figure 4b shows that the fraction of high-stress beads in stiff sheets is always larger than that in flexible ones; the curve related to 1-layer sheet is not monotonous, because the bonds in 1-layer sheets are broken ceaselessly with increasing strain. It should be noted that the aggregation of high-stress regions in stiff sheets also holds in the GrFs with the proportions of v = 30% and 70%, as shown in Figure S1 in Supplementary Information.
In order to compare the deformation behaviors of 1-layer flexible sheets and 8-layer stiff ones in the bi-modal GrFs with v = 50%, we calculate the average relative deformation for the three modes, bond stretching, in-plane shearing and out-of-plane bending of all 1-layer sheets and 8-layer ones by ε r ¯ = 1 / N r i = 1 N r | ( r i r 0 ) / r 0 | , ε φ ¯ = 1 / N φ i = 1 N φ | ( φ i φ 0 ) / φ 0 | and ε θ ¯ = 1 / N θ i = 1 N θ | ( θ i θ 0 ) / θ 0 | , respectively, where Nr, Nφ and Nθ are the number of bonds, the number of shearing angles and the number of bending angles in all flexible sheets or stiff ones. For the system with v = 50%, the part of flexible sheets shares the same value of Nr, Nφ and Nθ with the part of stiff ones. As shown in Figure 5, the average relative deformation for the three modes of flexible sheets (black line) is always larger than that of stiff sheets (red line) in the whole tension process. For example, in Figure 5a, the mean relative stretching length of bonds in flexible sheets is in the range of 0.002–0.0035 and is larger than that in stiff sheets at any tensile strain. For a more vivid comparison, the snapshots of the flexible part and stiff one at the tensile strain of 0.4 are color-coded by the local relative deformation, as shown in the inset in each figure. It is shown that high-strain regions mainly emerge in flexible sheets; in contrast to that, high-stress regions always emerge in stiff sheets in Figure 4, which is reasonable and easily understood. We also examine the deformation behaviors of flexible sheets and stiff ones in the bi-modal GrFs with proportions of v = 30% and 70% and obtain the same conclusion as shown in Figure S2 in Supplementary Information. Furthermore, from the snapshots in the three insets, we can see that the deformations of all sheets, especially flexible sheets, are very nonuniform, which is attributed to the high nonuniformity of foam structures and the heterogeneity of the mechanical properties of the two kinds of sheets.

3.4. Effect of the Proportion of Stiff Sheets on the Elastic Recovery of Bi-Modal GrFs

We also examine the elastic recovery of the GrFs with v = 0, 30%, 50%, 70% and 100% in the compression–uncompression process. The systems are first compressed to a maximum compressive strain of 0.8 at a loading rate of 4.3 × 106 s−1 in x direction and then uncompressed for a long time (20 ns) to fully relax. The stress–strain curves for the five GrFs in the compression–uncompression process are shown in Figure 6a, and none of them recover completely due to the microscopic plasticity unveiled in the works [47,56]. The residual strain for the five systems is 0.74, 0.39, 0.28, 0.17 and 0.07, respectively, which decreases with an increased v, as shown in the inset in Figure 6a. The corresponding snapshots of the compression–uncompression process of the three GrFs with v = 0, 50% and 100% are given in Figure 6b,d for intuitively showing the enhancement of the elastic recovery of GrFs with 8-layer sheets. This can be explained using our previous work [27] about the GrFs composed of mono-layer graphene sheets; at a given crosslink density, more elastic energy can be stored in the GrF composed of more stiff sheets as it is compressed. So, we can increase the proportion of stiff sheets in a bi-modal GrF to obtain a good elastic recovery ability in compression in practical applications.

4. Conclusions

We adopt coarse-grained molecular dynamics simulation to study the mechanical properties and deformation mechanisms of GrFs with bi-modal sheet thickness distribution. The modulus of GrFs increases almost linearly with the increased proportion of the 8-layer sheets, which can be well explained by the mixture rule. The strength of the bi-modal GrF decreases first and then increases with an increased proportion of 8-layer sheets; a minimum strength is obtained at a critical proportion of v ~ 30%. This is because, in the case of a smaller proportion of 8-layer sheets, v < 30%, first, stiff sheets do not form a complete bearing path in GrF, and the strength of the system is still mainly determined by the interconnected flexible sheets; second, due to the poor deformation coordination capacity of stiff sheets, when a certain number of stiff sheets are included in the GrF, more flexible sheets break in the tension process of the system. As the proportion of stiff sheets surpasses 30%, a complete bearing path form by stiff sheets, and the strength of the system is instead determined by the properties of stiff sheets. So, the strength of the system would increase with an increased proportion of stiff sheets. We also examine the distribution and evolution of the deformation energy of the stretching, shearing and bending of sheets and van der Waals adhesion energy between neighbor sheets in GrFs with varied proportion of 8-layer sheets. It is found that the stretching, bending and van der Waals energy would increase with the proportion of 8-layer sheets in GrFs. Furthermore, we examine the distribution of local von-Mises stress and the local stretching/shearing/bending strain of bi-modal GrFs. We find that a large number of high-stress regions non-uniformly distribute in 8-layer sheets and that only a small portion of them emerges in 1-layer sheets. In contrast, a large number of high-strain regions, no matter the local stretching, shearing or bending strain of sheets, non-uniformly locate in 1-layer sheets. Moreover, the elastic recovery of GrFs can be effectively enhanced by adding more 8-layer sheets in GrFs. These results would be helpful for us to understand and further guide the design of GrFs composed of both flexible and stiff graphene sheets.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ma14195622/s1, Figure S1: The distribution of the local von-Mises stress of the whole GrF system with the proportion of 8-layer sheet v = 30% and v = 70%. Figure S2: The comparison of the variation of bond length, in-plane shearing angle and out-of-plane bending angle in 1-layer soft sheets and 8-layer hard sheets with the proportion of 8-layer sheet v = 30% or v = 70%, respectively. The distribution of the bond strain, shear angle and bending angle in the whole GrF system with the proportion of 8-layer sheet v = 30% or v = 70% at the tensile strain of 0.4 for a direct comparison of the strain distribution in 1-layer soft sheets and 8-layer hard sheets.

Author Contributions

Conceptualization, C.W.; methodology, S.L. and M.L.; software, S.L. and M.L.; validation, S.L. and M.L.; formal analysis, C.W. and M.L.; investigation, M.L.; resources, C.W. and S.L.; data curation, M.L.; writing—original draft preparation, M.L.; writing—review and editing, C.W. and M.L.; visualization, M.L.; supervision, C.W. and S.L.; project administration, C.W. and S.L.; funding acquisition, C.W. and S.L. All authors have read and agreed to the published version of the manuscript.

Funding

S.L. acknowledges the support of NSFC through Grants #41472130 and the National Key Research and Development Program of China (Grant No. 2016YFC0600704). C.W. acknowledges the support of NSFC through Grants #11972348, #11790292, Strategic Priority Research Program of the Chinese Academy of Sciences (Grant No. XDB22040503).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing is not applicable.

Conflicts of Interest

The authors declare no competing interests.

Appendix A

Preparation of Numerical Model

Five steps are used to prepare the numerical bi-modal GrF with varied proportions of 8-layer stiff graphene sheets. The first step: 100 coarse-grained square graphene sheets with the side length of 75 nm are randomly placed in a much large cubic box to make sure there is no contact between them; the second step: an NPT assemble with 300 K by Langevin thermostat and 10 atmospheres by Berendsen barostat is used to make the system dense enough; under these conditions, the system reaches an equilibrium state in about 30 ns; the third step: a crosslink is added between the two neighbor beads if the distance between them is smaller than a critical value of 28 Å and if the total number of crosslinks in the system is 2751; the fourth step: for a given proportion of stiff sheets v, the same proportion of sheets are chosen randomly to set the mechanical properties of 8-layer stiff sheets given in Table 1, and other sheets are assigned the mechanical properties of 1-layer flexible sheets accordingly; the fifth step: an equilibrium simulation is conducted under the condition of 300 K and 1 atmosphere to obtain the final well-equilibrated state of the system with a total energy fluctuation convergence less than 1%.

References

  1. Sheng, L.; Wei, T.; Liang, Y.; Jiang, L.; Qu, L.; Fan, Z. Ultra-high toughness all graphene fibers derived from synergetic effect of interconnected graphene ribbons and graphene sheets. Carbon 2017, 120, 17–22. [Google Scholar] [CrossRef]
  2. Meng, F.; Lu, W.; Li, Q.; Byun, J.H.; Oh, Y.; Chou, T.W. Graphene-Based Fibers: A Review. Adv. Mater. 2015, 27, 5113–5131. [Google Scholar] [CrossRef] [PubMed]
  3. Xu, Z.; Zhang, Y.; Li, P.; Gao, C. Strong, conductive, lightweight, neat graphene aerogel fibers with aligned pores. ACS Nano 2012, 6, 7103–7113. [Google Scholar] [CrossRef]
  4. Xu, Z.; Sun, H.; Zhao, X.; Gao, C. Ultrastrong fibers assembled from giant graphene oxide sheets. Adv. Mater. 2013, 25, 188–193. [Google Scholar] [CrossRef]
  5. Zhang, L.; Alvarez, N.T.; Zhang, M.; Haase, M.; Malik, R.; Mast, D.; Shanov, V. Preparation and characterization of graphene paper for electromagnetic interference shielding. Carbon 2015, 82, 353–359. [Google Scholar] [CrossRef]
  6. Yavari, F.; Chen, Z.; Thomas, A.V.; Ren, W.; Cheng, H.M.; Koratkar, N. High sensitivity gas detection using a macroscopic three-dimensional graphene foam network. Sci. Rep. 2011, 1, 166. [Google Scholar] [CrossRef] [Green Version]
  7. Bi, H.; Xie, X.; Yin, K.; Zhou, Y.; Wan, S.; He, L.; Xu, F.; Banhart, F.; Sun, L.; Ruoff, R.S. Spongy graphene as a highly efficient and recyclable sorbent for oils and organic solvents. Adv. Funct. Mater. 2012, 22, 4421–4425. [Google Scholar] [CrossRef]
  8. Huang, X.; Qian, K.; Yang, J.; Zhang, J.; Li, L.; Yu, C.; Zhao, D. Functional nanoporous graphene foams with controlled pore sizes. Adv. Mater. 2012, 24, 4419–4423. [Google Scholar] [CrossRef]
  9. Hu, H.; Zhao, Z.; Wan, W.; Gogotsi, Y.; Qiu, J. Ultralight and highly compressible graphene aerogels. Adv. Mater. 2013, 25, 2219–2223. [Google Scholar] [CrossRef] [PubMed]
  10. Nieto, A.; Boesl, B.; Agarwal, A. Multi-scale intrinsic deformation mechanisms of 3D graphene foam. Carbon 2015, 85, 299–308. [Google Scholar] [CrossRef]
  11. Wu, Y.; Yi, N.; Huang, L.; Zhang, T.; Fang, S.; Chang, H.; Li, N.; Oh, J.; Lee, J.A.; Kozlov, M.; et al. Three-dimensionally bonded spongy graphene material with super compressive elasticity and near-zero Poisson’s ratio. Nat. Commun. 2015, 6, 6141. [Google Scholar] [CrossRef] [Green Version]
  12. Wang, J.; Shang, L.; Cheng, Y.; Ding, H.; Zhao, Y.; Gu, Z. Microfluidic Generation of Porous Particles Encapsulating Spongy Graphene for Oil Absorption. Small 2015, 11, 3890–3895. [Google Scholar] [CrossRef] [PubMed]
  13. Zhang, Q.; Zhang, F.; Medarametla, S.P.; Li, H.; Zhou, C.; Lin, D. 3D Printing of Graphene Aerogels. Small 2016, 12, 1702–1708. [Google Scholar] [CrossRef] [PubMed]
  14. Xu, X.; Zhang, Q.; Yu, Y.; Chen, W.; Hu, H.; Li, H. Naturally Dried Graphene Aerogels with Superelasticity and Tunable Poisson’s Ratio. Adv. Mater. 2016, 28, 9223–9230. [Google Scholar] [CrossRef] [PubMed]
  15. Gao, W.; Zhao, N.; Yao, W.; Xu, Z.; Bai, H.; Gao, C. Effect of flake size on the mechanical properties of graphene aerogels prepared by freeze casting. RSC Adv. 2017, 7, 33600–33605. [Google Scholar] [CrossRef] [Green Version]
  16. Jung, G.S.; Buehler, M.J. Multiscale Mechanics of Triply Periodic Minimal Surfaces of Three-Dimensional Graphene Foams. Nano Lett. 2018, 18, 4845–4853. [Google Scholar] [CrossRef]
  17. Shao, Q.; Tang, J.; Lin, Y.; Zhang, F.; Yuan, J.; Zhang, H.; Shinya, N.; Qin, L.-C. Synthesis and characterization of graphene hollow spheres for application in supercapacitors. J. Mater. Chem. A 2013, 1, 15423–15428. [Google Scholar] [CrossRef]
  18. Balázsi, K.; Furkó, M.; Liao, Z.; Fogarassy, Z.; Medved, D.; Zschech, E.; Dusza, J.; Balázsi, C. Graphene added multilayer ceramic sandwich (GMCS) composites: Structure, preparation and properties. J. Eur. Ceram. Soc. 2020, 40, 4792–4798. [Google Scholar] [CrossRef]
  19. Ma, Y.; Chang, H.; Zhang, M.; Chen, Y. Graphene-Based Materials for Lithium-Ion Hybrid Supercapacitors. Adv. Mater. 2015, 27, 5296–5308. [Google Scholar] [CrossRef]
  20. Wang, X.; Lu, M.; Wang, H.; Pei, Y.; Rao, H.; Du, X. Three-dimensional graphene aerogels-mesoporous silica frameworks for superior adsorption capability of phenols. Sep. Purif. Technol. 2015, 153, 7–13. [Google Scholar] [CrossRef]
  21. Chen, Z.; Ren, W.; Gao, L.; Liu, B.; Pei, S.; Cheng, H.-M. Three-dimensional flexible and conductive interconnected graphene networks grown by chemical vapour deposition. Nat. Mater. 2011, 10, 424–428. [Google Scholar] [CrossRef] [PubMed]
  22. Chen, S.; Bao, P.; Huang, X.; Sun, B.; Wang, G. Hierarchical 3D mesoporous silicon @ graphene nanoarchitectures for lithium ion batteries with superior performance. Nano Res. 2014, 7, 85–94. [Google Scholar] [CrossRef]
  23. Gibson, L.J.; Ashby, M.F. Cellular Solids: Structure and Properties; Cambridge University Press: Cambridge, UK, 1999; ISBN 9780521499118. [Google Scholar]
  24. Timoshenko, S.P.; Woinowsky-Krieger, S. Theory of Plates and Shells; McGraw-Hill College: London, UK, 1959; ISBN 9780070647794. [Google Scholar]
  25. Wang, G.; Dai, Z.; Xiao, J.; Feng, S.Z.; Weng, C.; Liu, L.; Xu, Z.; Huang, R.; Zhang, Z. Bending of Multilayer van der Waals Materials. Phys. Rev. Lett. 2019, 123, 116101. [Google Scholar] [CrossRef] [PubMed]
  26. Liu, F.; Wang, C.; Tang, Q. Conductivity Maximum in 3D Graphene Foams. Small 2018, 14, 1801458. [Google Scholar] [CrossRef] [PubMed]
  27. Wang, C.; Zhang, C.; Chen, S. Micro-mechanism and influencing factors of graphene foam elasticity. Carbon 2019, 148, 267–276. [Google Scholar] [CrossRef]
  28. Wang, C.; Zhang, C.; Chen, S. The microscopic deformation mechanism of 3D graphene foam materials under uniaxial compression. Carbon 2016, 109, 666–672. [Google Scholar] [CrossRef] [Green Version]
  29. Chae, S.J.; Güneş, F.; Kim, K.K.; Kim, E.S.; Han, G.H.; Kim, S.M.; Shin, H.; Yoon, S.M.; Choi, J.Y.; Park, M.H.; et al. Synthesis of large-area graphene layers on poly-nickel substrate by chemical vapor deposition: Wrinkle formation. Adv. Mater. 2009, 21, 2328–2333. [Google Scholar] [CrossRef]
  30. Qin, S.Y.; Liu, X.J.; Zhuo, R.X.; Zhang, X.Z. Microstructure-controllable graphene oxide hydrogel film based on a PH-responsive graphene oxide hydrogel. Macromol. Chem. Phys. 2012, 213, 2044–2051. [Google Scholar] [CrossRef]
  31. Xu, Y.; Sheng, K.; Li, C.; Shi, G. Self-assembled graphene hydrogel via a one-step hydrothermal process. ACS Nano 2010, 4, 4324–4330. [Google Scholar] [CrossRef] [PubMed]
  32. Zhao, Y.; Hu, C.; Hu, Y.; Cheng, H.; Shi, G.; Qu, L. A versatile, ultralight, nitrogen-doped graphene framework. Angew. Chem.-Int. Ed. 2012, 51, 11371–11375. [Google Scholar] [CrossRef]
  33. Yu, M.; Huang, Y.; Li, C.; Zeng, Y.; Wang, W.; Li, Y.; Fang, P.; Lu, X.; Tong, Y. Building three-dimensional graphene frameworks for energy storage and catalysis. Adv. Funct. Mater. 2015, 25, 324–330. [Google Scholar] [CrossRef]
  34. Lee, C.; Li, Q.; Kalb, W.; Liu, X.Z.; Berger, H.; Carpick, R.W.; Hone, J. Frictional characteristics of atomically thin sheets. Science 2010, 328, 76–80. [Google Scholar] [CrossRef] [Green Version]
  35. Baimova, J.A.; Rysaeva, L.K.; Liu, B.; Dmitriev, S.V.; Zhou, K. From flat graphene to bulk carbon nanostructures. Phys. Status Solidi 2015, 252, 1502–1507. [Google Scholar] [CrossRef]
  36. Patil, S.P.; Shendye, P.; Markert, B. Molecular Investigation of Mechanical Properties and Fracture Behavior of Graphene Aerogel. J. Phys. Chem. B 2020, 124, 6132–6139. [Google Scholar] [CrossRef]
  37. Qin, Z.; Jung, G.S.; Kang, M.J.; Buehler, M.J. The mechanics and design of a lightweight three-dimensional graphene assembly. Sci. Adv. 2017, 3, e1601536. [Google Scholar] [CrossRef] [Green Version]
  38. Zhang, X.; Zhong, L.; Mateos, A.; Kudo, A.; Vyatskikh, A.; Gao, H.; Greer, J.R.; Li, X. Theoretical strength and rubber-like behaviour in micro-sized pyrolytic carbon. Nat. Nanotechnol. 2019, 14, 762–769. [Google Scholar] [CrossRef] [Green Version]
  39. Liu, S.; Duan, K.; Li, L.; Wang, X.; Hu, Y. A multilayer coarse-grained molecular dynamics model for mechanical analysis of mesoscale graphene structures. Carbon 2021, 178, 528–539. [Google Scholar] [CrossRef]
  40. Shang, J.J.; Yang, Q.S.; Liu, X.; Wang, C. Compressive deformation mechanism of honeycomb-like graphene aerogels. Carbon 2018, 134, 398–410. [Google Scholar] [CrossRef] [Green Version]
  41. Shen, Z.; Ye, H.; Zhou, C.; Kröger, M.; Li, Y. Size of graphene sheets determines the structural and mechanical properties of 3D graphene foams. Nanotechnology 2018, 29, 104001. [Google Scholar] [CrossRef] [PubMed]
  42. Xia, W.; Vargas-Lara, F.; Keten, S.; Douglas, J.F. Structure and Dynamics of a Graphene Melt. ACS Nano 2018, 12, 5427–5435. [Google Scholar] [CrossRef] [PubMed]
  43. Mahdavi, S.M.; Adibnazari, S.; del Monte, F.; Gutiérrez, M.C. Compressive modulus and deformation mechanisms of 3DG foams: Experimental investigation and multiscale modeling. Nanotechnology 2021, 32, 485711. [Google Scholar] [CrossRef]
  44. Cranford, S.; Buehler, M.J. Twisted and coiled ultralong multilayer graphene ribbons. Model. Simul. Mater. Sci. Eng. 2011, 19, 54003. [Google Scholar] [CrossRef]
  45. Titov, A.V.; Král, P.; Pearson, R. Sandwiched Graphene−Membrane Superstructures. ACS Nano 2010, 4, 229–234. [Google Scholar] [CrossRef] [PubMed]
  46. Pan, D.; Wang, C.; Wang, T.C.; Yao, Y. Graphene Foam: Uniaxial Tension Behavior and Fracture Mode Based on a Mesoscopic Model. ACS Nano 2017, 11, 8988–8997. [Google Scholar] [CrossRef] [PubMed]
  47. Wang, C.; Pan, D.; Chen, S. Energy dissipative mechanism of graphene foam materials. Carbon 2018, 132, 641–650. [Google Scholar] [CrossRef] [Green Version]
  48. Bi, H.; Chen, I.W.; Lin, T.; Huang, F. A New Tubular Graphene Form of a Tetrahedrally Connected Cellular Structure. Adv. Mater. 2015, 27, 5943–5949. [Google Scholar] [CrossRef] [PubMed]
  49. Yang, T.; Wang, C.; Wu, Z. Crosslink-tuned large-deformation behavior and fracture mode in buckypapers. Carbon 2020, 159, 412–421. [Google Scholar] [CrossRef]
  50. Wang, C.; Wang, L.; Xu, Z. Enhanced mechanical properties of carbon nanotube networks by mobile and discrete binders. Carbon 2013, 64, 237–244. [Google Scholar] [CrossRef]
  51. Zhang, Y.Y.; Gu, Y.T. Mechanical properties of graphene: Effects of layer number, temperature and isotope. Comput. Mater. Sci. 2013, 71, 197–200. [Google Scholar] [CrossRef] [Green Version]
  52. Zhang, T.; Li, X.; Gao, H. Fracture of graphene: A review. Int. J. Fract. 2015, 196, 1–31. [Google Scholar] [CrossRef]
  53. Plimpton, S. Fast parallel algorithms for short-range molecular dynamics. J. Comput. Phys. 1995, 117, 1–19. [Google Scholar] [CrossRef] [Green Version]
  54. Stukowski, A. Visualization and analysis of atomistic simulation data with OVITO-the Open Visualization Tool. Model. Simul. Mater. Sci. Eng. 2010, 18, 015012. [Google Scholar] [CrossRef]
  55. Callister, W.D., Jr. Materials Science and Engineering: An Introduction; John Wiley & Sons Inc.: Hoboken, NJ, USA, 1999; ISBN 9780471320135. [Google Scholar]
  56. Nautiyal, P.; Boesl, B.; Agarwal, A. The mechanics of energy dissipation in a three-dimensional graphene foam with macroporous architecture. Carbon 2018, 132, 59–64. [Google Scholar] [CrossRef]
Figure 1. Coarse-grained model of bi-modal GrFs. (a-i) A square full-atomic graphene sheet with a side length of 2.5 nm; (a-ii) A square coarse-grained graphene sheet with a side length of 75 nm; (a-iii) three deformation modes of stretching, in-plane shearing and out-of-plane bending of graphene sheet; (b) The initial state of a well-equilibrated coarse-grained GrF consisting of the same number of 1-layer flexible (red) and 8-layer stiff (blue) coarse-grained graphene sheets; The crosslinks (green) are added between neighbor beads in different sheets. (c) Four typical microstructures of edge–edge (E–E), edge–surface (E–S), point–surface (P–S) and surface–surface (S–S) for two 1-layer sheets, a 1-layer sheet and a 8-layer sheet, and two 8-layer sheets; (d) The fraction of the four typical microstructures and the areal density (the contact area between sheets per unit volume) as an increased proportion of stiff sheets (the number of stiff sheets divided by the total number of sheets).
Figure 1. Coarse-grained model of bi-modal GrFs. (a-i) A square full-atomic graphene sheet with a side length of 2.5 nm; (a-ii) A square coarse-grained graphene sheet with a side length of 75 nm; (a-iii) three deformation modes of stretching, in-plane shearing and out-of-plane bending of graphene sheet; (b) The initial state of a well-equilibrated coarse-grained GrF consisting of the same number of 1-layer flexible (red) and 8-layer stiff (blue) coarse-grained graphene sheets; The crosslinks (green) are added between neighbor beads in different sheets. (c) Four typical microstructures of edge–edge (E–E), edge–surface (E–S), point–surface (P–S) and surface–surface (S–S) for two 1-layer sheets, a 1-layer sheet and a 8-layer sheet, and two 8-layer sheets; (d) The fraction of the four typical microstructures and the areal density (the contact area between sheets per unit volume) as an increased proportion of stiff sheets (the number of stiff sheets divided by the total number of sheets).
Materials 14 05622 g001
Figure 2. (a) The uniaxial tensile stress–strain curves of the GrFs with varied proportions of 8-layer stiff sheets v = 0, 30%, 50%, 70% and 100%; (b) The snapshots of the microstructural evolution of the GrFs with v = 0, 50% and 100% at four tensile strains of 0, 0.25, 0.5 and 0.75, respectively; (c) The modulus and strength of the GrFs as a function of the proportion of stiff sheets. (d) The connectivity of the GrFs with varied proportions of stiff sheets (The connectivity of the system is set to be 1 if the connected stiff sheets can run through the entire system in the loading direction; otherwise, it is 0); (e) The stress and the number of broken bonds in the GrFs with v = 0 and 30% as a function of tensile strain.
Figure 2. (a) The uniaxial tensile stress–strain curves of the GrFs with varied proportions of 8-layer stiff sheets v = 0, 30%, 50%, 70% and 100%; (b) The snapshots of the microstructural evolution of the GrFs with v = 0, 50% and 100% at four tensile strains of 0, 0.25, 0.5 and 0.75, respectively; (c) The modulus and strength of the GrFs as a function of the proportion of stiff sheets. (d) The connectivity of the GrFs with varied proportions of stiff sheets (The connectivity of the system is set to be 1 if the connected stiff sheets can run through the entire system in the loading direction; otherwise, it is 0); (e) The stress and the number of broken bonds in the GrFs with v = 0 and 30% as a function of tensile strain.
Materials 14 05622 g002
Figure 3. The distribution and evolution of three deformation energies, of stretching (a), in-plane shearing (b) and out-of-plane bending (c) of sheets, and the van der Waals potential energy (d) between neighboring sheets.
Figure 3. The distribution and evolution of three deformation energies, of stretching (a), in-plane shearing (b) and out-of-plane bending (c) of sheets, and the van der Waals potential energy (d) between neighboring sheets.
Materials 14 05622 g003
Figure 4. (a-i) The distribution of the local von-Mises stress of the whole GrF system with v = 50%; for a direct comparison of the stress distribution in 1-layer sheets and 8-layer sheets, the two parts are displayed separately in (a-ii) and (a-iii); (b) The comparison of the ratio of high stress beads in 1-layer sheets and 8-layer sheets in the GrF with v = 50%.
Figure 4. (a-i) The distribution of the local von-Mises stress of the whole GrF system with v = 50%; for a direct comparison of the stress distribution in 1-layer sheets and 8-layer sheets, the two parts are displayed separately in (a-ii) and (a-iii); (b) The comparison of the ratio of high stress beads in 1-layer sheets and 8-layer sheets in the GrF with v = 50%.
Materials 14 05622 g004
Figure 5. A comparison of the relative variation of (a) bond length, (b) in-plane shearing angle and (c) out-of-plane bending angle in 1-layer flexible sheets and 8-layer stiff sheets in the GrF with v = 50%, respectively. Inset: the distribution of the corresponding deformation in flexible sheets and stiff sheets at the tensile strain of 0.4.
Figure 5. A comparison of the relative variation of (a) bond length, (b) in-plane shearing angle and (c) out-of-plane bending angle in 1-layer flexible sheets and 8-layer stiff sheets in the GrF with v = 50%, respectively. Inset: the distribution of the corresponding deformation in flexible sheets and stiff sheets at the tensile strain of 0.4.
Materials 14 05622 g005
Figure 6. The effect of the proportion of 8-layer sheets on the elastic recovery of GrFs. (a) The stress-strain curves of GrFs with v = 0, 30%, 50%, 70% and 100% in the compression–uncompression process. Inset: the residual strain as a function of v; (bd) Four typical snapshots of the GrFs with v = 0, 50% and 100% at the compressive strain of 0, 0.4 and 0.8 and after recovery.
Figure 6. The effect of the proportion of 8-layer sheets on the elastic recovery of GrFs. (a) The stress-strain curves of GrFs with v = 0, 30%, 50%, 70% and 100% in the compression–uncompression process. Inset: the residual strain as a function of v; (bd) Four typical snapshots of the GrFs with v = 0, 50% and 100% at the compressive strain of 0, 0.4 and 0.8 and after recovery.
Materials 14 05622 g006
Table 1. The parameters of the main force field [44].
Table 1. The parameters of the main force field [44].
No. of Graphene Layers
Parameter18Units
k B 2351860kcal mol−1 Å−2
k C 18601860kcal mol−1 Å−2
r 0 2525Å
k φ 843567,480kcal mol−1 rad−2
φ 0 90°90°——
k θ 72.45466,543.5kcal mol−1 rad−2
θ 0 180°180°——
ε 473473kcal mol−1
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, S.; Lyu, M.; Wang, C. Mechanical Properties and Deformation Mechanisms of Graphene Foams with Bi-Modal Sheet Thickness by Coarse-Grained Molecular Dynamics Simulations. Materials 2021, 14, 5622. https://doi.org/10.3390/ma14195622

AMA Style

Liu S, Lyu M, Wang C. Mechanical Properties and Deformation Mechanisms of Graphene Foams with Bi-Modal Sheet Thickness by Coarse-Grained Molecular Dynamics Simulations. Materials. 2021; 14(19):5622. https://doi.org/10.3390/ma14195622

Chicago/Turabian Style

Liu, Shenggui, Mindong Lyu, and Chao Wang. 2021. "Mechanical Properties and Deformation Mechanisms of Graphene Foams with Bi-Modal Sheet Thickness by Coarse-Grained Molecular Dynamics Simulations" Materials 14, no. 19: 5622. https://doi.org/10.3390/ma14195622

APA Style

Liu, S., Lyu, M., & Wang, C. (2021). Mechanical Properties and Deformation Mechanisms of Graphene Foams with Bi-Modal Sheet Thickness by Coarse-Grained Molecular Dynamics Simulations. Materials, 14(19), 5622. https://doi.org/10.3390/ma14195622

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