Temperature-Dependent Mechanical Properties of Graphene/Cu Nanocomposites with In-Plane Negative Poisson’s Ratios

Temperature-Dependent Mechanical Properties of Graphene/Cu Nanocomposites with In-Plane Negative Poisson’s Ratios

PDF

Yin Fan¹, Yang Xiang¹, Hui-Shen Shen²^,³

Research. Vol 2020 Article ID 5618021

Less

Research. Vol 2020 Article ID 5618021

• Research Article •

Temperature-Dependent Mechanical Properties of Graphene/Cu Nanocomposites with In-Plane Negative Poisson’s Ratios

Full

Yin Fan¹, Yang Xiang¹, Hui-Shen Shen²^,³

Affiliations

¹School of Engineering, Western Sydney University, Locked Bag 1797, Penrith, NSW 2751, Australia

²School of Aeronautics and Astronautics, Shanghai Jiao Tong University, Shanghai 200240, China

³School of Ocean and Civil Engineering, Shanghai Jiao Tong University, Shanghai 200240, China

Published: 2020-02-20 doi: 10.34133/2020/5618021

Outline

Abstract

Less

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 $1 1 0$ is selected as the matrix and an embedded-atom 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 Erhart-Albe 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 Lennard-Jones (L-J) potential. The simulation results showed that both graphene and copper matrix possess in-plane NPRs. The temperature-dependent mechanical properties of graphene/copper nanocomposites with in-plane NPRs are obtained for the first time.

Cite this Article

Yin Fan, Yang Xiang, Hui-Shen Shen. Temperature-Dependent Mechanical Properties of Graphene/Cu Nanocomposites with In-Plane Negative Poisson’s Ratios[J]. Research, 2020 , 2020 (2) : 5618021 - . DOI: 10.34133/2020/5618021

Full Text

Less

1. 　Introduction

Less

Graphene is an ideal reinforcement agent for high-performance composites due to its superb mechanical and physicochemical properties [1–6]. Recently, polymer- and metal-based composites [7–21] reinforced by graphene have gained tremendous attention from the scientific community. Compared with polymer-based 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 graphene-reinforced 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 single-layer 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 10-100 layers of carbon atoms. It is noted that the material properties of these multilayer graphene derivatives are far less superior to those of the single-layer 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 man-made 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 out-of-plane 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 out-of-plane direction of a structural element rather than in the in-plane direction [35–37]. The research on composites with in-plane 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 in-plane auxetic composites by using the matrix and the reinforcing phase that both possess the in-plane 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 face-centered cubic (FCC) lattice might have an auxetic phenomenon in the

1 1 ¯ 0

direction when subjected to uniaxial

1 1 0

loading. To date, the NPR for single-layer 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 C-C 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 graphene-reinforced copper composites are established and used for the prediction of the temperature-dependent NPRs and other material properties for the composites.

2. 　MD Models

Less

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

1 1 0

1 1 ¯ 0

, and

0 0 1

crystal orientations as shown in Figure 1(a). The simulation box containing 5632 copper atoms is of

40.89 Å × 40.89 Å × 39.76 Å

in dimensions. The boundary conditions are set as periodical in the three directions. An embedded-atom 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

(1) E tot = ∑ i F i ∑ j ≠ i ρ j a R i j + 12 ∑ i ∑ j ≠ i ϕ i j R i j,

where

F i ρ

is the embedding energy of atom

i

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 in-plane and out-of-plane 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. 　Single-Layer Graphene

Graphene is a two-dimensional 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

201.68 Å × 200.21 Å

is built, which is, respectively, referred to as the zig-zag direction and the armchair direction. A hybrid empirical bond order potential [54] for the C-C 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:

(2) E C = ∑ i > j f C r i j V R r i j − b i j + b j i 2 V Λ r i j,

where

V R

and

V Λ

are the pairwise attractive and repulsive contributions, respectively,

f C

is the cutoff function,

b

is the bond order, and

r

is the distance between the nearest two carbon atoms. The detailed definitions of these physical quantities can be found in the literature [54]. The C-C 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 in-plane properties, such as

E 11

E 22

G 12

ν 12

, and

ν 21

, the simulation tests in three directions are carried out. The force

f

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

M

atoms in the zigzag direction and

N

atoms in the armchair direction, then we can first obtain stretching rigidities (

C 11

and

C 22

) and shear rigidity (

C 66

) as:

(3) C 11 = E 11 h = M f L x L y Δ L x, (4) C 22 = E 22 h = N f L y L x Δ L y, (5) C 66 = G 12 h = M f L y L x Δ L x

(6) C 66 = G 12 h = N f L x L y Δ L y,

where

L x

and

L y

are the origin lengths of the graphene in the zig-zag 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

Δ L x

Δ L y

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

q

and the central deflection

W ¯

based on the classical plate theory:

(7) q L x 4 D 11 h e = A W 1 W ¯ h e + A W 3 W ¯ h e 3 + ⋯,

where

D 11

is the flexural rigidity of graphene in the zig-zag direction.

A W 1

and

A W 3

are parameters and defined as

(8) A W 1 = π 6 16 1 + 2 c 4 β 2 + c 3 β 4, (9) A W 3 = 3 π 6 64 1 g 13 + c 3 β 4 g 31 1 + 2 c 4 β 2 + c 3 β 4 + 1 + c 3 β 4 1 − ν 12 G ν 21 G, (10) g 13 = 1 + 18 c 4 β 2 + 81 c 3 β 2, (11) g 31 = 81 + 18 c 4 β 2 + c 3 β 4, (12) β = L x L y, (13) c 3 = C 22 C 11, (14) c 4 = ν 21 G + 2 C 66 C 11 1 − ν 12 G ν 21 G,

where superscript

G

represents graphene. Finally, the effective thickness of the graphene sheet (

h e

) can be derived from Equation (7):

(15) h e = 12 1 − ν 12 G ν 21 G q L x 4 / C 11 − A W 3 W ¯ 3 A W 1 W ¯ 1 / 2 .

The detailed derivation for the load-deflection 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]:

(16) α 11 = 1 L x 0 Δ L x Δ T = 1 L 0 L x − L x 0 T − T 0, (17) α 22 = 1 L y 0 Δ L y Δ T = 1 L 0 L y − L y 0 T − T 0,

where

α

is the TEC with unit K^-1.

T

and

T 0

are the simulation temperature and reference temperature, respectively.

L x 0

and

L y 0

are the lengths of graphene in the zig-zag direction and the armchair direction, respectively, corresponding to the reference temperature

T 0

2.3. 　Graphene-Reinforced Copper Nanocomposite

The MD model of the graphene-reinforced 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 Cu-Cu and C-C. The Lennard-Jones (L-J) potential which is theoretically a good molecular dynamics model for long and short distances for neutral atoms and molecule is adapted for the C-Cu interaction. Moreover, the C-Cu interaction has already been demonstrated to meet the prediction from the Lennard-Jones 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 L-J potential is used:

(18) E i j = 4 ε σ r 12 − σ r 6 .

In Equation (18),

r

is the distance between two atoms. Obviously, this potential depends on the depth

ε

and the equilibrium interatomic distance

σ

. These parameters are determined to be

ε = 0.01996 eV

and

σ = 3.225 Å

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 isothermal-isobaric (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

Less

3.1. 　Mechanical Properties and NPR of Single Crystal Copper

As a typical FCC metal, single crystal copper shows an in-plane auxetic property when the load is applied in the

1 1 0

direction. Table 1 shows the in-plane properties of the single crystal copper. From the table, in-plane NPRs (

ν 12

and

ν 21

) 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 (

E 11

and

E 22

) and shear modulus (

G 12

) 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 out-of-plane properties of the single crystal copper, including

E 33

G 13

, and

G 23

, are listed in Table 2. It is found that the out-of-plane elastic modulus is nearly half of the in-plane elastic modulus. The out-of-plane shear moduli

G 13

and

G 23

are almost four times of the in-plane shear modulus.

3.2. 　In-Plane Properties and Out-Of-Plane Bending Behavior of the Graphene Sheet

The variation of the stretching rigidities, the shear rigidity, the in-plane 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 C-C 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

2 × 10 − 5 nm

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 atomic-scale two-dimensional 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

q

3.75 × 10 − 3 GPa

. Then, the density of the graphene can be obtained, which is 4717 kg/m³. Finally, the effective moduli of graphene can also be obtained based on Equations (3)–(6). Table 4 presents the elastic moduli

E 11

and

E 22

as well as the shear modulus

G 12

for the graphene sheet.

3.3. 　Material Properties of Graphene/Cu Nanocomposites

In this section, five different graphene weight fractions for the graphene-reinforced 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]

(19) V G = w G w G + ρ G / ρ C − ρ G / ρ C w G,

where

V

and

w

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³ 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 in-plane 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 zig-zag direction of the graphene as well as the

1 1 0

crystal direction of the copper while the 22 direction is along the armchair direction of the graphene as well as the

1 1 ¯ 0

crystal direction of the copper. We found that

E 11

is slightly larger than

E 22

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 out-of-plane moduli,

E 33

G 13

, and

G 23

, with different

V G

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 in-plane Young’s moduli. It can be found that

E 33

is only improved by 7% with 13% volume fraction of graphene when compared with that of the pure copper. Unlike the in-plane shear modulus, the out-of-plane shear moduli,

G 13

and

G 23

, are actually weakened by adding graphene reinforcement. We observe that the out-of-plane 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 graphene-reinforced 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 in-plane auxetic. Unlike in [21] where out-of-plane NPRs were obtained for graphene/Cu composites, in the current study, the in-plane NPRs of the resulting composite are obtained. From Figure 4, it can be found that the in-plane NPRs (

ν 12

and

ν 21

) 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 ν₁₂ 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

α 11

and

α 22

for all three materials being very close, only their

α 11

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 graphene-reinforced 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 L-J potential in the present MD models of graphene-reinforced 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

Less

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 in-plane auxeticity are selected as reinforcements and matrix, and their detailed material properties are obtained through MD simulations. To determine the effective material properties of graphene-reinforced 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 in-plane 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

V G

is 5%.

Funding

Less

Australian Research Council(DP140104156)
National Natural Science Foundation of China(51779138)

References

Less

Reddy

C. D.

, Rajendran

, and Liew

K. M.

, “Equilibrium configuration and continuum elastic properties of finite sized graphene,” Nanotechnology, vol. 17, no. 3, pp. 864–870, 2006

Sun

, Li

Z.-J.

, and Zhang

X.-B.

, “Achieving of high density/utilization of active groups via synergic integration of C=N and C=O bonds for ultra-stable and high-rate lithium-ion batteries,” Research, vol. 2018, article 1936735, pp. 1–10, 2018

Geim

A. K.

, and Novoselov

K. S.

, “The rise of graphene,” Nature Materials, vol. 6, no. 3, pp. 183–191, 2007

Lee

, Wei

, Kysar

J. W.

, and Hone

, “Measurement of the elastic properties and intrinsic strength of monolayer graphene,” Science, vol. 321, no. 5887, pp. 385–388, 2008

Huang

, Yin

, Wu

, Qi

, He

, Zhang

, Yan

, Boey

, and Zhang

, “Graphene-based materials: synthesis, characterization, properties, and applications,” Small, vol. 7, no. 14, pp. 1876–1902, 2011

, Fang

, Zhang

, Tang

, Zhou

, and Yu

, “Chemically binding scaffolded anodes with 3D graphene architectures realizing fast and stable lithium storage,” Research, vol. 2019, article 8393085, pp. 1–9, 2019

Lin

, Xiang

, and Shen

H.-S.

, “Temperature dependent mechanical properties of graphene reinforced polymer nanocomposites - A molecular dynamics simulation,” Composites Part B: Engineering, vol. 111, pp. 261–269, 2017

Han

, Meng

, Araby

, Liu

, and 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

, Yu

, Jiang

, Zhou

, Deng

, and Xiao

, “A novel method for preparing graphene nanosheets/Al composites by accumulative extrusion-bonding process,” Carbon, vol. 152, pp. 932–945, 2019

, Tong

, Lin

, Chen

, Guo

, Xu

, and Zhou

, “Graphene-reinforced metal matrix nanocomposites–a review,” Materials Science and Technology, vol. 32, no. 9, pp. 930–953, 2016

Dixit

, Mahata

, Mahapatra

D. R.

, Kailas

S. V.

, and Chattopadhyay

, “Multi-layer graphene reinforced aluminum - Manufacturing of high strength composite by friction stir alloying,” Composites Part B: Engineering, vol. 136, pp. 63–71, 2018

Jagannadham

, “Thermal conductivity of copper-graphene composite films synthesized by electrochemical deposition with exfoliated graphene platelets,” Metallurgical and Materials Transactions B, vol. 43, no. 2, pp. 316–324, 2012

Hwang

, Yoon

, Jin

S. H.

, Lee

, Kim

T. S.

, Hong

S. H.

, and Jeon

, “Enhanced mechanical properties of graphene/copper nanocomposites using a molecular-level mixing process,” Advanced Materials, vol. 25, no. 46, pp. 6724–6729, 2013

Zhang

, Heim

F. M.

, Bartlett

J. L.

, Song

, Isheim

, and Li

, “Bioinspired, graphene-enabled Ni composites with high strength and toughness,” Science Advances, vol. 5, no. 5, article eaav5577, 2019

Naseer

, Ahmad

, Aslam

, Guan

B. H.

, Harun

W. S. W.

, Muhamad

, Raza

M. R.

, and German

R. M.

, “A review of processing techniques for graphene-reinforced metal matrix composites,” Materials and Manufacturing Processes, vol. 34, no. 9, pp. 957–985, 2019

Mokhalingam

, Kumar

, and Srivastava

, “Mechanical behaviour of graphene reinforced aluminum nano composites,” Materials Today: Proceedings, vol. 4, no. 2, pp. 3952–3958, 2017

Rezaei

, “Tensile mechanical characteristics and deformation mechanism of metal-graphene nanolayered composites,” Computational Materials Science, vol. 151, pp. 181–188, 2018

Montazeri

, and 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

Weng

, Ning

, Fu

, Hu

, Zhao

, Huang

, and Peng

, “Molecular dynamics study of strengthening mechanism of nanolaminated graphene/Cu composites under compression,” Scientific Reports, vol. 8, no. 1, p. 3089, 2018

Sharma

, Kumar

, and Chandra

, “Mechanical and thermal properties of graphene–carbon nanotube-reinforced metal matrix composites: a molecular dynamics study,” Journal of Composite Materials, vol. 51, no. 23, pp. 3299–3313, 2017

Zhang

, Lu

, Pei

, Li

, Wang

, and Tieu

, “The negative Poisson's ratio and strengthening mechanism of nanolayered graphene/Cu composites,” Carbon, vol. 143, pp. 125–137, 2019

Evans

K. E.

, Nkansah

M. A.

, Hutchinson

I. J.

, and Rogers

S. C.

, “Molecular network design,” Nature, vol. 353, no. 6340, p. 124, 1991

Fíla

, Zlámal

, Jiroušek

, Falta

, Koudelka

, Kytýř

, Doktor

, and Valach

, “Impact testing of polymer-filled auxetics using split Hopkinson pressure bar,” Advanced Engineering Materials, vol. 19, no. 10, article 1700076, 2017

Duc

N. D.

, Seung-Eock

, Tuan

N. D.

, Tran

, and Khoa

N. D.

, “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

Karathanasopoulos

, Reda

, and Ganghoffer

, “Designing two-dimensional metamaterials of controlled static and dynamic properties,” Computational Materials Science, vol. 138, pp. 323–332, 2017

Imbalzano

, Linforth

, Ngo

T. D.

, Lee

P. V. S.

, and Tran

, “Blast resistance of auxetic and honeycomb sandwich panels: comparisons and parametric designs,” Composite Structures, vol. 183, pp. 242–261, 2018

Jiang

, and Hu

, “Finite element modeling of multilayer orthogonal auxetic composites under low-velocity impact,” Materials, vol. 10, no. 8, p. 908, 2017

Alderson

K. L.

, Pickles

A. P.

, Neale

P. J.

, and Evans

K. E.

, “Auxetic polyethylene: the effect of a negative Poisson's ratio on hardness,” Acta Metallurgica et Materialia, vol. 42, no. 7, pp. 2261–2266, 1994

Choi

J. B.

, and Lakes

R. S.

, “Fracture toughness of re-entrant foam materials with a negative Poisson's ratio: experiment and analysis,” International Journal of Fracture, vol. 80, no. 1, pp. 73–83, 1996

Howell

, Prendergast

, and Hansen

, “Examination of acoustic behavior of negative Poisson's ratio materials,” Applied Acoustics, vol. 43, no. 2, pp. 141–148, 1994

Huang

, and Chen

, “Negative Poisson's ratio in modern functional materials,” Advanced Materials, vol. 28, no. 37, pp. 8079–8096, 2016

Saxena

K. K.

, Das

, and Calius

E. P.

, “Three decades of auxetics research − materials with negative Poisson's ratio: a review,” Advanced Engineering Materials, vol. 18, no. 11, pp. 1847–1870, 2016

Caddock

B. D.

, and Evans

K. E.

, “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

Herakovich

C. T.

, “Composite laminates with negative through-the-thickness Poisson's ratios,” Journal of Composite Materials, vol. 18, no. 5, pp. 447–455, 1984

, Shen

H.-S.

, and Wang

, “Nonlinear bending of sandwich beams with functionally graded negative Poisson's ratio honeycomb core,” Composite Structures, vol. 212, pp. 317–325, 2019

, Shen

H.-S.

, and 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

, Shen

H.-S.

, and Wang

, “Thermal post-buckling of sandwich beams with functionally graded negative Poisson's ratio honeycomb core,” International Journal of Mechanical Sciences, vol. 152, pp. 289–297, 2019

Nguyen

T. D.

, and Duc

N. D.

, “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

Gao

, Dong

, Li

, and Ren

, “Novel two-dimensional silicon dioxide with in-plane negative Poisson's ratio,” Nano Letters, vol. 17, no. 2, pp. 772–777, 2017

Alderson

, and Alderson

K. L.

, “Auxetic materials,” Proceedings of the Institution of Mechanical Engineers, Part G: Journal of Aerospace Engineering, vol. 221, no. 4, pp. 565–575, 2007

Lakes

, “Advances in negative Poisson's ratio materials,” Advanced Materials, vol. 5, no. 4, pp. 293–296, 1993

Milstein

, and Huang

, “Existence of a negative Poisson ratio in FCC crystals,” Physical Review B, vol. 19, no. 4, pp. 2030–2033, 1979

Zakharchenko

K. V.

, Katsnelson

M. I.

, and Fasolino

, “Finite temperature lattice properties of graphene beyond the quasiharmonic approximation,” Physical Review Letters, vol. 102, no. 4, article 046808, 2009

Jiang

J.-W.

, Chang

, Guo

, and Park

H. S.

, “Intrinsic negative Poisson’s ratio for single-layer graphene,” Nano Letters, vol. 16, no. 8, pp. 5286–5290, 2016

Fan

, Xiang

, and Shen

H.-S.

, “Temperature-dependent negative Poisson’s ratio of monolayer graphene: prediction from molecular dynamics simulations,” Nanotechnology Reviews, vol. 8, no. 1, pp. 415–421, 2019

Scarpa

, Adhikari

, and Phani

A. S.

, “Effective elastic mechanical properties of single layer graphene sheets,” Nanotechnology, vol. 20, no. 6, article 065709, 2009

Qin

, Qin

, and Hu

, “Origin of anisotropic negative Poisson’s ratio in graphene,” Nanoscale, vol. 10, no. 22, pp. 10365–10370, 2018

Korobeynikov

S. N.

, Alyokhin

V. V.

, and Babichev

A. V.

, “On the molecular mechanics of single layer graphene sheets,” International Journal of Engineering Science, vol. 133, pp. 109–131, 2018

Phliptom

, “Fast parallel algorithms for short-range molecular dynamics,” Journal of Computional Physics, vol. 117, no. 1, pp. 1–19, 1995

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

Lebedeva

I. V.

, Minkin

A. S.

, Popov

A. M.

, and Knizhnik

A. A.

, “Elastic constants of graphene: comparison of empirical potentials and DFT calculations,” Physica E: Low-dimensional Systems and Nanostructures, vol. 108, pp. 326–338, 2019

Tersoff

, “Empirical interatomic potential for carbon, with applications to amorphous carbon,” Physical Review Letters, vol. 61, no. 25, pp. 2879–2882, 1988

Brenner

D. W.

, “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

Erhart

, and Albe

, “Analytical potential for atomistic simulations of silicon, carbon, and silicon carbide,” Physical Review B, vol. 71, no. 3, article 035211, 2005

Foiles

S. M.

, Baskes

M. I.

, and Daw

M. S.

, “Embedded-atom-method 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

Daw

M. S.

, and Baskes

M. I.

, “Embedded-atom method: derivation and application to impurities, surfaces, and other defects in metals,” Physical Review B, vol. 29, no. 12, pp. 6443–6453, 1984

Shen

H.-S.

, “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

Islam

M. Z.

, Mahboob

, Lowe

R. L.

, and Bechtel

E. S.

, “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

Karadeniz

Z. H.

, and 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

Chang

S.-W.

, Nair

A. K.

, and Buehler

M. J.

, “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

Hong

, Li

, Zeng

X. C.

, and Zhang

, “Tuning thermal contact conductance at graphene-copper interface via surface nanoengineering,” Nanoscale, vol. 7, no. 14, pp. 6286–6294, 2015

Wang

, Wang

, Liu

, Crimp

M. A.

, Wang

, and Qu

, “Anisotropic thermal expansion coefficient of multilayer graphene reinforced copper matrix composites,” Journal of Alloys and Compounds, vol. 755, pp. 114–122, 2018

Zhang

, Xu

, Liu

, and Luo

S. N.

, “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

Montazeri

, and Panahi

, “MD-based estimates of enhanced load transfer in graphene/metal nanocomposites through Ni coating,” Applied Surface Science, vol. 457, pp. 1072–1080, 2018

Huang

S. P.

, Mainardi

D. S.

, and Balbuena

P. B.

, “Structure and dynamics of graphite-supported bimetallic nanoclusters,” Surface Science, vol. 545, no. 3, pp. 163–179, 2003

Rafiee

M. A.

, Rafiee

, Wang

, Song

, Yu

Z. Z.

, and Koratkar

, “Enhanced mechanical properties of nanocomposites at low graphene content,” ACS Nano, vol. 3, no. 12, pp. 3884–3890, 2009

Shen

H.-S.

A Two-Step Perturbation Method in Nonlinear Analysis of Beams, Plates and Shells, John Wiley & Sons Inc, 2013

Appendix

Less

Year 2020 volume 2020 Issue 2

PDF

Cite this Article

BibTeX

Article Info

doi: 10.34133/2020/5618021

Receive Date：2019-11-02
Online Date：2026-10-27
Published：2020-02-20

Article Data

Affiliations

History

Received：2019-11-02
Accepted：2020-01-19

Funding

Australian Research Council(DP140104156)

National Natural Science Foundation of China(51779138)

Affiliations

¹School of Engineering, Western Sydney University, Locked Bag 1797, Penrith, NSW 2751, Australia

²School of Aeronautics and Astronautics, Shanghai Jiao Tong University, Shanghai 200240, China

³School of Ocean and Civil Engineering, Shanghai Jiao Tong University, Shanghai 200240, China

References

Share

https://castjournals.cast.org.cn/joweb/research/EN/10.34133/2020/5618021

Share to

Scan QR to access full text

Cite this article

BibTeX

Citations

表12种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT

Table 1: In-plane material properties of single crystal copper under different temperatures (11 and 22 directions are, respectively, referred to as $1 1 0$ and $1 1 ¯ 0$ crystal orientations).

Temperature (K)	$E 11$ (GPa)	$E 22$ (GPa)	$G 12$ (GPa)	$ν 12$	$ν 21$	$α$ (×10^-6 K^-1)
300	106.33	106.71	18.648	-0.1976	-0.1949	16.3766
500	95.576	96.224	16.684	-0.2164	-0.2142	17.3444
700	86.733	87.289	14.544	-0.2352	-0.2413	18.3063
1000	70.642	70.674	11.615	-0.2857	-0.2905	20.6236

Table 2: Out-of-plane material properties of single crystal copper under different temperatures (11 and 22 directions are, respectively, referred to as $1 1 0$ and $1 1 ¯ 0$ crystal orientations).

Temperature (K)	$E 33$ (GPa)	$G 13$ (GPa)	$G 23$ (GPa)
300	54.306	68.526	69.585
500	49.792	63.267	63.553
700	43.934	59.295	59.722
1000	39.966	50.74	50.923

Table 3: In-plane properties of monolayer graphene under different temperatures (11 and 22 directions are, respectively, referred to as zig-zag and armchair directions).

Temperature (K)	$C 11$ (GPa nm)	$C 22$ (GPa nm)	$C 66$ (GPa nm)	$ν 12$	$ν 21$	$α 11$ (×10^-6 K^-1)	$α 22$ (×10^-6 K^-1)
300	357.27	339.82	128.91	-0.0409	-0.0392	-11.952	-11.948
500	353.24	330.04	131.71	-0.0412	-0.0389	-9.6269	-9.6267
700	349.71	327.29	137.62	-0.0436	-0.0411	-8.1185	-8.1177
1000	321.54	305.84	143.73	-0.0501	-0.0482	-6.6417	-6.6397

Table 4: Moduli of monolayer graphene under different temperatures (11 and 22 directions are, respectively, referred to as zig-zag and armchair directions).

Temperature (K)	$E 11$ (TPa)	$E 22$ (TPa)	$G 12$ (TPa)
300	2.2470	2.1372	0.8108
500	2.2216	2.0757	0.8284
700	2.1994	2.0584	0.8655
1000	2.0223	1.9235	0.9040

Table 5: The weight fractions and their corresponding volume fractions (the mass densities of graphene and copper are 4717 kg/m³ and 8960 kg/m³, respectively).

Weight fraction	Volume fraction
2.8%	5%
3.7%	7%
4.8%	9%
6.1%	11%
7.1%	13%

Table 6: Temperature-dependent properties of graphene-reinforced copper composites (11 and 22 directions are, respectively, referred to as zig-zag and armchair directions, and the 33 direction is referred to as the thickness direction).

Temperature (K)	$V G$	$E 11$ (GPa)	$E 22$ (GPa)	$E 33$ (GPa)	$G 12$ (GPa)	$G 13$ (GPa)	$G 23$ (GPa)
300	5%	207.55	196.69	61.454	66.389	33.617	32.327
	7%	249.62	235.67	68.840	80.214	26.881	26.895
	9%	275.44	256.48	69.651	91.706	25.517	24.949
	11%	307.06	286.97	71.640	108.42	23.616	22.802
	13%	319.77	304.99	73.993	120.34	20.001	19.445

500	5%	193.15	183.94	55.773	62.092	31.536	31.009
	7%	235.50	221.18	61.819	77.036	25.845	25.809
	9%	258.90	243.71	68.271	87.200	24.257	23.893
	11%	293.37	273.66	69.588	100.92	21.853	20.786
	13%	296.39	288.72	72.029	112.29	19.283	17.660

700	5%	180.50	171.58	47.192	58.314	28.928	28.200
	7%	219.61	206.17	53.817	72.686	24.115	23.587
	9%	242.91	229.60	62.517	82.040	23.627	22.989
	11%	271.55	254.96	67.880	97.670	17.003	16.290
	13%	284.59	275.48	69.542	110.79	18.057	17.396

1000	5%	155.35	148.86	40.818	53.053	24.109	23.701
	7%	197.96	178.74	49.213	65.355	22.997	22.394
	9%	220.43	199.26	60.216	76.439	21.923	21.839
	11%	248.34	210.12	64.923	93.537	16.211	15.926
	13%	253.87	223.40	66.792	99.681	16.731	15.988

Table 7: Temperature-dependent NPRs and TECs of graphene-reinforced copper composites (11 and 22 directions are, respectively, referred to as zig-zag and armchair directions, and the 33 direction is referred to as the thickness direction).

Temperature (K)	$V G$	$ν 12$	$ν 13$	$ν 23$	$α 11$ (×10^-6 K^-1)	$α 22$ (×10^-6 K^-1)
300	5%	-0.0649	0.6297	0.6512	1.4224	1.4194
	7%	-0.0594	0.5631	0.6034	1.0812	1.0813
	9%	-0.0486	0.5111	0.5088	0.9976	1.0021
	11%	-0.0364	0.4235	0.4598	0.8204	0.8204
	13%	-0.0380	0.4026	0.4348	0.7984	0.8014

500	5%	-0.0721	0.6298	0.6617	1.5037	1.5006
	7%	-0.0718	0.5726	0.6358	1.3848	1.3957
	9%	-0.0540	0.5494	0.5584	1.3174	1.3046
	11%	-0.0532	0.5345	0.5019	1.2199	1.2053
	13%	-0.0489	0.5025	0.4951	1.0232	1.0288

700	5%	-0.0751	0.6826	0.6742	1.6222	1.6278
	7%	-0.0736	0.6452	0.6294	1.4324	1.4286
	9%	-0.0525	0.5828	0.5763	1.3589	1.3489
	11%	-0.0581	0.5643	0.5011	1.2069	1.2066
	13%	-0.0271	0.5576	0.4317	1.1367	1.1422

1000	5%	-0.1118	0.6978	0.7063	1.7022	1.7331
	7%	-0.0994	0.6572	0.6453	1.7205	1.6899
	9%	-0.0875	0.6033	0.5791	1.5242	1.5416
	11%	-0.0676	0.5750	0.5854	1.3028	1.3016
	13%	-0.0415	0.5633	0.5955	1.3256	1.3083

Figure 1: The schematic diagrams of the (a) single crystal copper cell, (b) single layer graphene, and (c) graphene-reinforced copper composite.

Figure 2: The variation of in-plane Young’s moduli and shear modulus of composites with different graphene volume fractions at room temperature.

Figure 3: The variation of out-of-plane Young’s modulus and shear moduli of composites with different graphene volume fractions at room temperature.

Figure 4: The effect of temperature on NPRs of the composite ( $V G = 9 %$ ) and its component materials.

Figure 5: The effect of graphene volume fraction $V G$ on $ν 12$ of the composite under different temperature conditions.

Figure 6: The effect of temperature on TEC of the composite ( $V G = 9 %$ ) and its component materials.

Articles: Latest Articles; Most Read; Collections

Updates: Events; News; Multimedia

About: About Us

Contact

No. 86 Xueyuan South Road, Haidian District, Beijing

100081

010-62662699(Operation and Technical Support)

qkjq@cast.org.cn

Copyright © 2025 China Association for Science and Technology. All rights reserved. For all open access content, the relevant licensing terms apply.
Sponsored by the Office of the Leading Group for Cybersecurity and Informatization of CAST, and supported by Science and Technology Review Publishing House