Research Article  Open Access
Yin Fan, Yang Xiang, HuiShen Shen, "TemperatureDependent Mechanical Properties of Graphene/Cu Nanocomposites with InPlane Negative Poisson’s Ratios", Research, vol. 2020, Article ID 5618021, 12 pages, 2020. https://doi.org/10.34133/2020/5618021
TemperatureDependent Mechanical Properties of Graphene/Cu Nanocomposites with InPlane Negative Poisson’s Ratios
Abstract
Negative Poisson’s ratio (NPR), also known as “auxetic”, is a highly desired property in a wide range of future industry applications. By employing molecular dynamics (MD) simulation, metal matrix nanocomposites reinforced by graphene sheets are studied in this paper. In the simulation, single crystal copper with crystal orientation is selected as the matrix and an embeddedatom method (EAM) potential is used to describe the interaction of copper atoms. An aligned graphene sheet is selected as reinforcement, and a hybrid potential, namely, the ErhartAlbe potential, is used for the interaction between a pair of carbon atoms. The interaction between the carbon atom and copper atom is approximated by the LennardJones (LJ) potential. The simulation results showed that both graphene and copper matrix possess inplane NPRs. The temperaturedependent mechanical properties of graphene/copper nanocomposites with inplane NPRs are obtained for the first time.
1. Introduction
Graphene is an ideal reinforcement agent for highperformance composites due to its superb mechanical and physicochemical properties [1–6]. Recently, polymer and metalbased composites [7–21] reinforced by graphene have gained tremendous attention from the scientific community. Compared with polymerbased composites, composites with metal matrix can provide higher stiffness and strength and service under a higher temperature, which makes them have the potential for more critical applications in industry, especially in the aerospace and automobile industry. Several researchers have reported in recent years [9–21] on experimental and numerical studies of graphenereinforced metal composites. The numerical studies [16–21] mainly based on the method of molecular dynamics (MD) simulations indicated that the material properties of metal composites can be significantly improved by applying singlelayer graphene reinforcement. In most experimental reports, the graphene used to reinforce metal matrix composites (MMCs) is mainly multilayered graphene (MLG) or graphene nanoplatelets (GNPs) which contain 10100 layers of carbon atoms. It is noted that the material properties of these multilayer graphene derivatives are far less superior to those of the singlelayer graphene counterpart [14]. It was also observed [15] that increasing the number of graphene layers may degrade the mechanical and other properties of these MMCs.
Negative Poisson’s ratio (NPR), also known as “auxetic”, is a phenomenon associated with a material that experiences transverse expansion when a tensile load is applied in the longitudinal direction. On the other hand, if the material is under a compressive load, it contracts transversely. The NPR materials are a novel class of materials, which have advantages that conventional materials may not possess [22]. Materials with NPR have triggered significant interest in the scientific community because of their unique physical properties [23–27] including high shear modulus, high indentation resistance [28], large plane strain fracture toughness [29], and excellent shock and sound absorption abilities [30, 31], and all of which are highly desirable properties for engineering applications. Specifically, the NPR materials can be used in microelectromechanical systems (MEMS) and fastening devices due to their auxetic characteristics, and they can be also added into engineering structures (e.g., sandwich or laminated panels in the aerospace industry) as energy absorption components. The natural materials found with NPR are polymers with microporous characteristics. Such polymer consists of a randomly or periodically distributed 3D array of particles interconnected by fibrils in the microscopic scale. For instance, microporous polytetrafluoroethylene (PTFE) film has a network of individual fibrils which can be extended by several microns in both directions when under uniaxial tension [32, 33]. On the other hand, auxetic materials could be achieved via designing and tailoring their microscopic compositions and structures. Such manmade materials consist of a large number of unit cells, which are arranged into certain periodical truss patterns such as reentrant type, chiral type, or rotating units [32]. Recently, it has been reported that laminated composites reinforced with microfibers can gain outofplane NPR [34] through design and optimization of stacking sequences. However, there are still some limitations in engineering applications for the above NPR materials. For example, the majority of materials fabricated only have NPR in certain directions where the stiffness is low in those directions. Hence, such materials may only be used as the core for sandwiches which need to be enhanced by composite face sheets. In addition, NPR is usually achieved in the outofplane direction of a structural element rather than in the inplane direction [35–37]. The research on composites with inplane NPR is still in the conceptual stage [38, 39] due to the unavailability of such type of NPR materials and their associated material properties for structural applications. Nevertheless, as hinted in the literature [40], there is a possible approach to obtain inplane auxetic composites by using the matrix and the reinforcing phase that both possess the inplane auxetic property.
In most polycrystalline metal materials, the value of Poisson’s ratio is in the vicinity of 1/3. However, it is different for metal materials with single crystals where NPR may exist in certain crystallographic directions [41]. Early research [42] suggested that some single crystal metals including copper with a facecentered cubic (FCC) lattice might have an auxetic phenomenon in the direction when subjected to uniaxial loading. To date, the NPR for singlelayer graphene has been observed by researchers using different approaches, including the Monte Carlo method [43], molecular dynamics and statics [44, 45], molecular structural mechanics [46], first principle theory [47], and membrane theory [48].
In this study, the classical MD method using the LAMMPS package [49] is employed to construct the single crystal copper and the software package OVITO [50] is utilized to generate graphic presentations from the MD simulation results. For graphene, discrepancies exist in the MD studies on Poisson’s ratio due to the selection of different potentials for the CC bond. From existing literature [51], the absolute NPR (0.158) can be obtained by MD simulations with the Tersoff potential [52] but cannot be achieved with the potentials based on a reactive empirical bond order (REBO) potential [53]. Hence, a hybrid potential [54] based on both REBO and Tersoff potentials is used to address this issue. Finally, the MD models of graphenereinforced copper composites are established and used for the prediction of the temperaturedependent NPRs and other material properties for the composites.
2. MD Models
2.1. Single Crystal Copper
Copper is an ideal substrate that supports growing graphene layers. In this study, we create a MD model of single crystal copper along the , , and crystal orientations as shown in Figure 1(a). The simulation box containing 5632 copper atoms is of in dimensions. The boundary conditions are set as periodical in the three directions. An embeddedatom method (EAM) potential [55] developed by Dawn and Baskes [56] is employed to create the interaction between each pair of two copper atoms. The mathematical form of the general total energy for a pair of copper atoms is where is the embedding energy of atom which is a function of the atomic electron density , is the pair potential interaction, and and are the element types of atoms I and J. In the case of this study, and are the same. The lattice constant and the atomic weight of the FCC copper are selected as 3.615 Å and 63.55, respectively.
In the simulation, the tensile and shear behaviors controlled by strain are applied to this MD model to study its inplane and outofplane mechanical properties. The temperature varying from 300 K to 1000 K is also considered in the simulation to observe the effect of temperature on the mechanical properties of the single crystal copper.
2.2. SingleLayer Graphene
Graphene is a twodimensional material. However, we cannot ignore its effective thickness in the mechanical analysis of Young’s modulus and shear modulus. The MD model of the monolayer graphene with is built, which is, respectively, referred to as the zigzag direction and the armchair direction. A hybrid empirical bond order potential [54] for the CC interaction based on REBO [53] and Tersoff [52] potentials is employed herein. The cohesive energy of this potential is written as a sum over individual bond energies: where and are the pairwise attractive and repulsive contributions, respectively, is the cutoff function, is the bond order, and is the distance between the nearest two carbon atoms. The detailed definitions of these physical quantities can be found in the literature [54]. The CC bond length is well known as 1.42 Å, and the atom weight is 12.011. For the sake of convenience, this hybrid potential is named as the Erhart and Albe potential (EAP).
To determine the inplane properties, such as , , , , and , the simulation tests in three directions are carried out. The force is only applied on the carbon atoms closest to the boundary for each simulation, as shown in Figure 1(b). If we assume that there are atoms in the zigzag direction and atoms in the armchair direction, then we can first obtain stretching rigidities ( and ) and shear rigidity () as: or where and are the origin lengths of the graphene in the zigzag direction and the armchair direction, respectively, and denotes the deformation caused by the force. It is worth noting that Equations (5) and (6) are valid when or is very small.
There is divergence on the effective thickness of graphene which affects the values of its moduli. Usually, the graphite thickness is used to estimate the effective thickness of graphene. However, this method is not accurate for a mechanical concept. In this study, we treat the graphene sheet as a thin plate in the framework of mechanics and the effective thickness of the graphene sheet will be obtained after the central deflection of graphene is measured under a uniform load. Shen [57] has derived the relationship between the load and the central deflection based on the classical plate theory: where is the flexural rigidity of graphene in the zigzag direction. and are parameters and defined as where superscript represents graphene. Finally, the effective thickness of the graphene sheet () can be derived from Equation (7):
The detailed derivation for the loaddeflection relationship in Equation (7) is presented in the appendix. In the simulation, the boundary atoms are fixed and a small transverse force is applied to the rest of the atoms.
Once the effective thickness is determined by MD simulations, the effective Young’s modulus and shear modulus of the graphene sheet can be obtained uniquely from Equations (3) to (6). The TECs at different temperatures can also be obtained from the relaxation process. For the computation of the TEC, the standard relationship is used [58, 59]: where is the TEC with unit K^{1}. and are the simulation temperature and reference temperature, respectively. and are the lengths of graphene in the zigzag direction and the armchair direction, respectively, corresponding to the reference temperature .
2.3. GrapheneReinforced Copper Nanocomposite
The MD model of the graphenereinforced copper nanocomposite is illustrated in Figure 1(c). The EAM and EAP potentials are also used in the simulation of composite materials for the interactions of CuCu and CC. The LennardJones (LJ) potential which is theoretically a good molecular dynamics model for long and short distances for neutral atoms and molecule is adapted for the CCu interaction. Moreover, the CCu interaction has already been demonstrated to meet the prediction from the LennardJones potential in previous investigations [60, 61]. This potential has also been widely used in existing literature [62–64] on the MD study of carbon/copper molecular structures. To estimate the interaction between copper and carbon atoms, the LJ potential is used:
In Equation (18), is the distance between two atoms. Obviously, this potential depends on the depth and the equilibrium interatomic distance . These parameters are determined to be and based on existing studies [65]. On the other hand, the common divisor of the two atoms’ lattice parameters is considered in the geometric design of the MD models for the composites to reduce the effect of mismatch, which will lead to internal stress in the MD models during the phases of relaxation and loading.
In the simulation, five different weight fractions (2.8%, 3.7%, 4.8%, 6.1%, and 7.1%) of the graphene sheet are taken into account. The temperature effect is also considered in the simulation. The time step is set as 1 fs, and the thermal relaxation consumes 40000 steps within the context of an isothermalisobaric (NPT) ensemble. The tensile deformation of the simulation box is controlled by a strain rate of 10^{6}/fs, and the whole load process takes 100000 steps while it takes only 10000 steps for the shear deformation.
3. Results and Discussion
3.1. Mechanical Properties and NPR of Single Crystal Copper
As a typical FCC metal, single crystal copper shows an inplane auxetic property when the load is applied in the direction. Table 1 shows the inplane properties of the single crystal copper. From the table, inplane NPRs ( and ) are obtained through linear fitting MD simulation results and it is found that auxeticity is more significant with the increase in temperature. However, the elastic moduli ( and ) and shear modulus () are decreased as temperature rises. We also find that the shear modulus is much smaller than the elastic moduli. The thermal expansion coefficient is also obtained and listed in Table 1. The outofplane properties of the single crystal copper, including , , and , are listed in Table 2. It is found that the outofplane elastic modulus is nearly half of the inplane elastic modulus. The outofplane shear moduli and are almost four times of the inplane shear modulus.


3.2. InPlane Properties and OutOfPlane Bending Behavior of the Graphene Sheet
The variation of the stretching rigidities, the shear rigidity, the inplane Poisson’s ratios, and the thermal expansion coefficients of the graphene from 300 K to 1000 K is shown in Table 3. It can be seen that the stretching rigidities are decreased slightly as temperature increases. However, the trend of the shear rigidity is opposite to that of the stretching rigidities, which is in agreement with the MD results of Lin et al. [7]. Like the single crystal copper, the NPR of the graphene sheet is increased when temperature rises. However, the thermal expansion coefficients of the graphene sheet are negative, which means that graphene contracts when temperature increases. The effect of different loading directions on NPR has been investigated in our previous work [45] where the changes of NPRs with an increased load in the zigzag and armchair directions were also studied. We found that the gap between the NPRs in the zigzag and armchair directions is remarkable when the load is small but the gap is gradually reduced with the increase in the load. It was believed that the geometric configuration in different chiral directions plays a key role, because the distance between neighboring carbon atoms is increased once deformation of the graphene is occurred, which causes the reduction of the CC bonds. With the increase in the interatomic distance, the effect of geometric configuration becomes weak, which results in the reduction of the gap between the NPRs in those two directions.

The bending behavior is then performed to determine the effective thickness of graphene. As mentioned above, the small force is set to be 10^{4} nN and 1 million steps are taken for the balance of the whole system. The center deflection of graphene is calculated as the average value among transverse displacements of the nearest carbon atoms to the center point. In our study, the strain load of 10^{6}/fs is applied on the MD model and the averaged lattice parameters for the time step are converged to within 0.0001 nm. The strain load 10^{6}/fs with unit fs^{1} means that we apply about as an incremental displacement per femtosecond on the MD model of graphene. The equilibrium state depends on the difference of lattice parameters within a time step. If the difference is lower than 0.0001 nm, the whole model can be regarded to be in an equilibrium state. The graphene is an atomicscale twodimensional hexagonal lattice made by carbon atoms. In a dimensionally fixed graphene, the geometric center point must be in a hexagon cell. In our computation, the transverse displacements of the six carbon atoms of the hexagon cell in the center area of graphene are averaged to obtain the deflection of the graphene for each time point. Applying Equation (15), we obtain the effective thickness of the graphene being 0.159 nm with the deflection of 5.784 Å at room temperature when the equivalent applied pressure is . Then, the density of the graphene can be obtained, which is 4717 kg/m^{3}. Finally, the effective moduli of graphene can also be obtained based on Equations (3)–(6). Table 4 presents the elastic moduli and as well as the shear modulus for the graphene sheet.

3.3. Material Properties of Graphene/Cu Nanocomposites
In this section, five different graphene weight fractions for the graphenereinforced copper composites will be taken into consideration. The calculation of the graphene weight fractions can be achieved based on the number of carbon and copper atoms and their respective atom weight. Then, the volume fractions can be directly obtained by [66] where and denote the volume fraction and the weight fraction, respectively. is the mass density. The superscripts or subscripts G and C represent graphene and copper, respectively. For example, one of our MD models for composites contains 112064 atoms in all, of which 14720 atoms are carbon and the rest are copper atoms. Accordingly, the obtained graphene weight fraction is 2.8%. Taking 8960 kg/m^{3} as the mass density of copper, the volume fraction of graphene is then calculated to be 5%. The considered five weight fractions and their corresponding volume fractions are listed in Table 5.

Figure 2 illustrates the effect of graphene volume fraction on the inplane Young’s moduli and shear modulus of the composites at room temperature (300 K). It is worth noting that the 11 direction is along the zigzag direction of the graphene as well as the crystal direction of the copper while the 22 direction is along the armchair direction of the graphene as well as the crystal direction of the copper. We found that is slightly larger than and they are both improved as the volume fraction of graphene is increased. We also notice that the rising trend of the moduli becomes slower when the graphene volume is higher than 7%. The outofplane moduli, , , and , with different are depicted in Figure 3. Although the Young’s modulus in the thickness direction is enhanced by increasing the volume fraction of graphene, the incremental rate is relatively smaller when compared with that of the inplane Young’s moduli. It can be found that is only improved by 7% with 13% volume fraction of graphene when compared with that of the pure copper. Unlike the inplane shear modulus, the outofplane shear moduli, and , are actually weakened by adding graphene reinforcement. We observe that the outofplane shear properties are degraded by as much as half when the composite contains only 5% volume fraction of graphene. The detailed moduli of the composites with different volume fractions of graphene under various temperatures are presented in Table 6.

In order to study the auxetic behavior of the graphenereinforced copper composites, the NPRs of the composite with 9% graphene volume fraction and its component materials are all shown in Figure 4. As mentioned above, the graphene sheet and single crystal copper are both inplane auxetic. Unlike in [21] where outofplane NPRs were obtained for graphene/Cu composites, in the current study, the inplane NPRs of the resulting composite are obtained. From Figure 4, it can be found that the inplane NPRs ( and ) of the composite are between those of graphene and copper but closer to the graphene NPRs. With increase in temperature, the NPRs of the composite are decreased and the NPRs can reach about 0.09 at 1000 K. The effect of graphene volume fraction on ν_{12} of the composite is shown in Figure 5. Unlike the effect of temperature, the increase in graphene content may weaken auxetic characteristics of the composite. The TECs of the copper and graphene as well as the composite varying from 300 K to 1000 K are shown in Figure 6. Due to and for all three materials being very close, only their is compared in Figure 6. Obviously, the TEC of graphene is negative but the TEC of copper is positive. The TEC of the composite is still positive but lower than that of the copper. In other words, the addition of graphene with negative TECs decreases the TECs of the copper based composite. Hence, internal stress may exist in the graphenereinforced copper composites once temperature is varied. It is observed that the TEC trend of the composite is slightly different from that of graphene or copper due to the coupling effect. Table 7 shows detailed NPRs and TECs of composites from 300 K to 1000 K. It is worth noting that the interaction between the carbon atom and copper atom is approximated by the LJ potential in the present MD models of graphenereinforced composites. Hence, the influence of graphene reinforcement on the mechanical properties and Poison’s ratio of the composites mainly depends on the mechanical loading transfer at the interface of the graphene sheet and copper matrix.

4. Conclusion
This study not only presents novel MD results for MMCs reinforced by graphene sheets but also provides a new approach for the design of auxetic materials and structures. The graphene sheet and single crystal copper both showing inplane auxeticity are selected as reinforcements and matrix, and their detailed material properties are obtained through MD simulations. To determine the effective material properties of graphenereinforced metal matrix composites (GRMMCs), the MD simulations are carried out under different graphene volume fractions and temperature ranges from 300 K to 1000 K. The results showed that the inplane NPRs of GRMMCs are decreased as the temperature rises or the graphene volume fraction increases and the maximum NPRs may reach 0.11 at 1000 K when the volume fraction of graphene is 5%.
Appendix
A graphene sheet with length , width , and effective thickness under a distributed transverse pressure is considered. When the effective thickness is thin enough, the classical thin plate theory is adapted and the displacement field of an arbitrary point (, ) on the graphene can be expressed as where , , and are midplane displacements in a Cartesian coordinate (, , and ) and is the time. Then, the corresponding strains in the von Kármán type can be written as
The stressstrain relationship in a plane stress state can be expressed as where
By employing the Hamilton principle, we have where where is the area of the plate. Using the variational principle, the von Kármántype nonlinear equations for the thin plate can be derived as [57] where is the stress function defined by , , and and the linear operators are defined by and contains the geometric nonlinearity terms in the von Kármán sense and can be given as in Equations (A.9)–(A.11); , , and () are reduced stiffness matrices which are defined as where , , and are the stiffness coefficients defined by
The solutions of Equations (A.7) and (A.8) can be expanded in the perturbation forms as
By substituting Equations (A.15)–(A.17) into Equations (A.7) and (A.8) and using a twostep perturbation approach [67], the solution of load can be written as
Equation (A.18) is the loaddeflection relationship that we used in this study to obtain the effective thickness of the graphene. Substituting into Equation (A.18) and multiplying both sides of Equation (A.18) by , one has
From Equation (A.19), the final expression of can be obtained as
Conflicts of Interest
The authors declare that there is no conflict of interest.
Acknowledgments
This study was supported by the National Natural Science Foundation of China (Grant 51779138) and the Australian Research Council (Grant DP140104156). The authors are very grateful for these financial supports.
References
 C. D. Reddy, S. Rajendran, and K. M. Liew, “Equilibrium configuration and continuum elastic properties of finite sized graphene,” Nanotechnology, vol. 17, no. 3, pp. 864–870, 2006. View at: Publisher Site  Google Scholar
 T. Sun, Z.J. Li, and X.B. Zhang, “Achieving of high density/utilization of active groups via synergic integration of C=N and C=O bonds for ultrastable and highrate lithiumion batteries,” Research, vol. 2018, article 1936735, pp. 1–10, 2018. View at: Publisher Site  Google Scholar
 A. K. Geim and K. S. Novoselov, “The rise of graphene,” Nature Materials, vol. 6, no. 3, pp. 183–191, 2007. View at: Publisher Site  Google Scholar
 C. Lee, X. Wei, J. W. Kysar, and J. Hone, “Measurement of the elastic properties and intrinsic strength of monolayer graphene,” Science, vol. 321, no. 5887, pp. 385–388, 2008. View at: Publisher Site  Google Scholar
 X. Huang, Z. Yin, S. Wu et al., “Graphenebased materials: synthesis, characterization, properties, and applications,” Small, vol. 7, no. 14, pp. 1876–1902, 2011. View at: Publisher Site  Google Scholar
 P. Wu, Z. Fang, A. Zhang et al., “Chemically binding scaffolded anodes with 3D graphene architectures realizing fast and stable lithium storage,” Research, vol. 2019, article 8393085, pp. 1–9, 2019. View at: Publisher Site  Google Scholar
 F. Lin, Y. Xiang, and H.S. Shen, “Temperature dependent mechanical properties of graphene reinforced polymer nanocomposites  A molecular dynamics simulation,” Composites Part B: Engineering, vol. 111, pp. 261–269, 2017. View at: Publisher Site  Google Scholar
 S. Han, Q. Meng, S. Araby, T. Liu, and M. Demiral, “Mechanical and electrical properties of graphene and carbon nanotube reinforced epoxy adhesives: experimental and numerical analysis,” Composites Part A: Applied Science and Manufacturing, vol. 120, pp. 116–126, 2019. View at: Publisher Site  Google Scholar
 G. Wu, Z. Yu, L. Jiang et al., “A novel method for preparing graphene nanosheets/Al composites by accumulative extrusionbonding process,” Carbon, vol. 152, pp. 932–945, 2019. View at: Publisher Site  Google Scholar
 Z. Hu, G. Tong, D. Lin et al., “Graphenereinforced metal matrix nanocomposites–a review,” Materials Science and Technology, vol. 32, no. 9, pp. 930–953, 2016. View at: Publisher Site  Google Scholar
 S. Dixit, A. Mahata, D. R. Mahapatra, S. V. Kailas, and K. Chattopadhyay, “Multilayer graphene reinforced aluminum  Manufacturing of high strength composite by friction stir alloying,” Composites Part B: Engineering, vol. 136, pp. 63–71, 2018. View at: Publisher Site  Google Scholar
 K. Jagannadham, “Thermal conductivity of coppergraphene composite films synthesized by electrochemical deposition with exfoliated graphene platelets,” Metallurgical and Materials Transactions B, vol. 43, no. 2, pp. 316–324, 2012. View at: Publisher Site  Google Scholar
 J. Hwang, T. Yoon, S. H. Jin et al., “Enhanced mechanical properties of graphene/copper nanocomposites using a molecularlevel mixing process,” Advanced Materials, vol. 25, no. 46, pp. 6724–6729, 2013. View at: Publisher Site  Google Scholar
 Y. Zhang, F. M. Heim, J. L. Bartlett, N. Song, D. Isheim, and X. Li, “Bioinspired, grapheneenabled Ni composites with high strength and toughness,” Science Advances, vol. 5, no. 5, article eaav5577, 2019. View at: Publisher Site  Google Scholar
 A. Naseer, F. Ahmad, M. Aslam et al., “A review of processing techniques for graphenereinforced metal matrix composites,” Materials and Manufacturing Processes, vol. 34, no. 9, pp. 957–985, 2019. View at: Publisher Site  Google Scholar
 A. Mokhalingam, D. Kumar, and A. Srivastava, “Mechanical behaviour of graphene reinforced aluminum nano composites,” Materials Today: Proceedings, vol. 4, no. 2, pp. 3952–3958, 2017. View at: Publisher Site  Google Scholar
 R. Rezaei, “Tensile mechanical characteristics and deformation mechanism of metalgraphene nanolayered composites,” Computational Materials Science, vol. 151, pp. 181–188, 2018. View at: Publisher Site  Google Scholar
 A. Montazeri and A. Mobarghei, “Nanotribological behavior analysis of graphene/metal nanocomposites via MD simulations: new concepts and underlying mechanisms,” Journal of Physics and Chemistry of Solids, vol. 115, pp. 49–58, 2018. View at: Publisher Site  Google Scholar
 S. Weng, H. Ning, T. Fu et al., “Molecular dynamics study of strengthening mechanism of nanolaminated graphene/Cu composites under compression,” Scientific Reports, vol. 8, no. 1, p. 3089, 2018. View at: Publisher Site  Google Scholar
 S. Sharma, P. Kumar, and R. Chandra, “Mechanical and thermal properties of graphene–carbon nanotubereinforced metal matrix composites: a molecular dynamics study,” Journal of Composite Materials, vol. 51, no. 23, pp. 3299–3313, 2017. View at: Publisher Site  Google Scholar
 C. Zhang, C. Lu, L. Pei, J. Li, R. Wang, and K. Tieu, “The negative Poisson's ratio and strengthening mechanism of nanolayered graphene/Cu composites,” Carbon, vol. 143, pp. 125–137, 2019. View at: Publisher Site  Google Scholar
 K. E. Evans, M. A. Nkansah, I. J. Hutchinson, and S. C. Rogers, “Molecular network design,” Nature, vol. 353, no. 6340, p. 124, 1991. View at: Publisher Site  Google Scholar
 T. Fíla, P. Zlámal, O. Jiroušek et al., “Impact testing of polymerfilled auxetics using split Hopkinson pressure bar,” Advanced Engineering Materials, vol. 19, no. 10, article 1700076, 2017. View at: Publisher Site  Google Scholar
 N. D. Duc, K. SeungEock, N. D. Tuan, P. Tran, and N. D. Khoa, “New approach to study nonlinear dynamic response and vibration of sandwich composite cylindrical panels with auxetic honeycomb core layer,” Aerospace Science and Technology, vol. 70, pp. 396–404, 2017. View at: Publisher Site  Google Scholar
 N. Karathanasopoulos, H. Reda, and J. Ganghoffer, “Designing twodimensional metamaterials of controlled static and dynamic properties,” Computational Materials Science, vol. 138, pp. 323–332, 2017. View at: Publisher Site  Google Scholar
 G. Imbalzano, S. Linforth, T. D. Ngo, P. V. S. Lee, and P. Tran, “Blast resistance of auxetic and honeycomb sandwich panels: comparisons and parametric designs,” Composite Structures, vol. 183, pp. 242–261, 2018. View at: Publisher Site  Google Scholar
 L. Jiang and H. Hu, “Finite element modeling of multilayer orthogonal auxetic composites under lowvelocity impact,” Materials, vol. 10, no. 8, p. 908, 2017. View at: Publisher Site  Google Scholar
 K. L. Alderson, A. P. Pickles, P. J. Neale, and K. E. Evans, “Auxetic polyethylene: the effect of a negative Poisson's ratio on hardness,” Acta Metallurgica et Materialia, vol. 42, no. 7, pp. 2261–2266, 1994. View at: Publisher Site  Google Scholar
 J. B. Choi and R. S. Lakes, “Fracture toughness of reentrant foam materials with a negative Poisson's ratio: experiment and analysis,” International Journal of Fracture, vol. 80, no. 1, pp. 73–83, 1996. View at: Publisher Site  Google Scholar
 B. Howell, P. Prendergast, and L. Hansen, “Examination of acoustic behavior of negative Poisson's ratio materials,” Applied Acoustics, vol. 43, no. 2, pp. 141–148, 1994. View at: Publisher Site  Google Scholar
 C. Huang and L. Chen, “Negative Poisson's ratio in modern functional materials,” Advanced Materials, vol. 28, no. 37, pp. 8079–8096, 2016. View at: Publisher Site  Google Scholar
 K. K. Saxena, R. Das, and E. P. Calius, “Three decades of auxetics research − materials with negative Poisson's ratio: a review,” Advanced Engineering Materials, vol. 18, no. 11, pp. 1847–1870, 2016. View at: Publisher Site  Google Scholar
 B. D. Caddock and K. E. Evans, “Microporous materials with negative Poisson's ratios. I. Microstructure and mechanical properties,” Journal of Physics D: Applied Physics, vol. 22, no. 12, pp. 1877–1882, 1989. View at: Publisher Site  Google Scholar
 C. T. Herakovich, “Composite laminates with negative throughthethickness Poisson's ratios,” Journal of Composite Materials, vol. 18, no. 5, pp. 447–455, 1984. View at: Publisher Site  Google Scholar
 C. Li, H.S. Shen, and H. Wang, “Nonlinear bending of sandwich beams with functionally graded negative Poisson's ratio honeycomb core,” Composite Structures, vol. 212, pp. 317–325, 2019. View at: Publisher Site  Google Scholar
 C. Li, H.S. Shen, and H. Wang, “Nonlinear vibration of sandwich beams with functionally graded negative Poisson's ratio honeycomb core,” International Journal of Structural Stability and Dynamics, vol. 19, no. 3, article 1950034, 2019. View at: Publisher Site  Google Scholar
 C. Li, H.S. Shen, and H. Wang, “Thermal postbuckling of sandwich beams with functionally graded negative Poisson's ratio honeycomb core,” International Journal of Mechanical Sciences, vol. 152, pp. 289–297, 2019. View at: Publisher Site  Google Scholar
 T. D. Nguyen and N. D. Duc, “Evaluation of elastic properties and thermal expansion coefficient of composites reinforced by randomly distributed spherical particles with negative Poisson's ratio,” Composite Structures, vol. 153, pp. 569–577, 2016. View at: Publisher Site  Google Scholar
 Z. Gao, X. Dong, N. Li, and J. Ren, “Novel twodimensional silicon dioxide with inplane negative Poisson's ratio,” Nano Letters, vol. 17, no. 2, pp. 772–777, 2017. View at: Publisher Site  Google Scholar
 A. Alderson and K. L. Alderson, “Auxetic materials,” Proceedings of the Institution of Mechanical Engineers, Part G: Journal of Aerospace Engineering, vol. 221, no. 4, pp. 565–575, 2007. View at: Publisher Site  Google Scholar
 R. Lakes, “Advances in negative Poisson's ratio materials,” Advanced Materials, vol. 5, no. 4, pp. 293–296, 1993. View at: Publisher Site  Google Scholar
 F. Milstein and K. Huang, “Existence of a negative Poisson ratio in FCC crystals,” Physical Review B, vol. 19, no. 4, pp. 2030–2033, 1979. View at: Publisher Site  Google Scholar
 K. V. Zakharchenko, M. I. Katsnelson, and A. Fasolino, “Finite temperature lattice properties of graphene beyond the quasiharmonic approximation,” Physical Review Letters, vol. 102, no. 4, article 046808, 2009. View at: Publisher Site  Google Scholar
 J.W. Jiang, T. Chang, X. Guo, and H. S. Park, “Intrinsic negative Poisson’s ratio for singlelayer graphene,” Nano Letters, vol. 16, no. 8, pp. 5286–5290, 2016. View at: Publisher Site  Google Scholar
 Y. Fan, Y. Xiang, and H.S. Shen, “Temperaturedependent negative Poisson’s ratio of monolayer graphene: prediction from molecular dynamics simulations,” Nanotechnology Reviews, vol. 8, no. 1, pp. 415–421, 2019. View at: Publisher Site  Google Scholar
 F. Scarpa, S. Adhikari, and A. S. Phani, “Effective elastic mechanical properties of single layer graphene sheets,” Nanotechnology, vol. 20, no. 6, article 065709, 2009. View at: Publisher Site  Google Scholar
 Z. Qin, G. Qin, and M. Hu, “Origin of anisotropic negative Poisson’s ratio in graphene,” Nanoscale, vol. 10, no. 22, pp. 10365–10370, 2018. View at: Publisher Site  Google Scholar
 S. N. Korobeynikov, V. V. Alyokhin, and A. V. Babichev, “On the molecular mechanics of single layer graphene sheets,” International Journal of Engineering Science, vol. 133, pp. 109–131, 2018. View at: Publisher Site  Google Scholar
 S. Phliptom, “Fast parallel algorithms for shortrange molecular dynamics,” Journal of Computional Physics, vol. 117, no. 1, pp. 1–19, 1995. View at: Publisher Site  Google Scholar
 A. Stukowski, “Visualization and analysis of atomistic simulation data with OVITO  the Open Visualization Tool,” Modelling and Simulation in Materials Science and Engineering, vol. 18, no. 1, article 015012, 2009. View at: Publisher Site  Google Scholar
 I. V. Lebedeva, A. S. Minkin, A. M. Popov, and A. A. Knizhnik, “Elastic constants of graphene: comparison of empirical potentials and DFT calculations,” Physica E: Lowdimensional Systems and Nanostructures, vol. 108, pp. 326–338, 2019. View at: Publisher Site  Google Scholar
 J. Tersoff, “Empirical interatomic potential for carbon, with applications to amorphous carbon,” Physical Review Letters, vol. 61, no. 25, pp. 2879–2882, 1988. View at: Publisher Site  Google Scholar
 D. W. Brenner, “Empirical potential for hydrocarbons for use in simulating the chemical vapor deposition of diamond films,” Physical Review B, vol. 42, no. 15, pp. 9458–9471, 1990. View at: Publisher Site  Google Scholar
 P. Erhart and K. Albe, “Analytical potential for atomistic simulations of silicon, carbon, and silicon carbide,” Physical Review B, vol. 71, no. 3, article 035211, 2005. View at: Publisher Site  Google Scholar
 S. M. Foiles, M. I. Baskes, and M. S. Daw, “Embeddedatommethod functions for the fcc metals Cu, Ag, Au, Ni, Pd, Pt, and their alloys,” Physical Review B, vol. 33, no. 12, pp. 7983–7991, 1986. View at: Publisher Site  Google Scholar
 M. S. Daw and M. I. Baskes, “Embeddedatom method: derivation and application to impurities, surfaces, and other defects in metals,” Physical Review B, vol. 29, no. 12, pp. 6443–6453, 1984. View at: Publisher Site  Google Scholar
 H.S. Shen, “Nonlinear analysis of composite laminated thin plates subjected to lateral pressure and thermal loading and resting on elastic foundations,” Composite Structures, vol. 49, no. 2, pp. 115–128, 2000. View at: Publisher Site  Google Scholar
 M. Z. Islam, M. Mahboob, R. L. Lowe, and E. S. Bechtel, “Characterization of the thermal expansion properties of graphene using molecular dynamics simulations,” Journal of Physics D: Applied Physics, vol. 46, no. 43, article 435302, 2013. View at: Publisher Site  Google Scholar
 Z. H. Karadeniz and D. Kumlutas, “A numerical study on the coefficients of thermal expansion of fiber reinforced composite materials,” Composite Structures, vol. 78, no. 1, pp. 1–10, 2007. View at: Publisher Site  Google Scholar
 S.W. Chang, A. K. Nair, and M. J. Buehler, “Geometry and temperature effects of the interfacial thermal conductance in copper– and nickel–graphene nanocomposites,” Journal of Physics: Condensed Matter, vol. 24, no. 24, article 245301, 2012. View at: Publisher Site  Google Scholar
 Y. Hong, L. Li, X. C. Zeng, and J. Zhang, “Tuning thermal contact conductance at graphenecopper interface via surface nanoengineering,” Nanoscale, vol. 7, no. 14, pp. 6286–6294, 2015. View at: Publisher Site  Google Scholar
 X. Wang, X. Wang, M. Liu, M. A. Crimp, Y. Wang, and Z. Qu, “Anisotropic thermal expansion coefficient of multilayer graphene reinforced copper matrix composites,” Journal of Alloys and Compounds, vol. 755, pp. 114–122, 2018. View at: Publisher Site  Google Scholar
 S. Zhang, Y. Xu, X. Liu, and S. N. Luo, “Competing roles of interfaces and matrix grain size in the deformation and failure of polycrystalline Cu–graphene nanolayered composites under shear loading,” Physical Chemistry Chemical Physics, vol. 20, no. 36, pp. 23694–23701, 2018. View at: Publisher Site  Google Scholar
 A. Montazeri and B. Panahi, “MDbased estimates of enhanced load transfer in graphene/metal nanocomposites through Ni coating,” Applied Surface Science, vol. 457, pp. 1072–1080, 2018. View at: Publisher Site  Google Scholar
 S. P. Huang, D. S. Mainardi, and P. B. Balbuena, “Structure and dynamics of graphitesupported bimetallic nanoclusters,” Surface Science, vol. 545, no. 3, pp. 163–179, 2003. View at: Publisher Site  Google Scholar
 M. A. Rafiee, J. Rafiee, Z. Wang, H. Song, Z. Z. Yu, and N. Koratkar, “Enhanced mechanical properties of nanocomposites at low graphene content,” ACS Nano, vol. 3, no. 12, pp. 3884–3890, 2009. View at: Publisher Site  Google Scholar
 H.S. Shen, A TwoStep Perturbation Method in Nonlinear Analysis of Beams, Plates and Shells, John Wiley & Sons Inc, 2013.
Copyright
Copyright © 2020 Yin Fan et al. Exclusive Licensee Science and Technology Review Publishing House. Distributed under a Creative Commons Attribution License (CC BY 4.0).