International audienceWe characterize the macroscopic e ffective behavior of a graphene sheet modeled by a hexagonal lattice of elastic bars, using Gamma-convergence

Annie Raoult

Hervé Le Dret

Comptes Rendus Mathématique

C. R. Acad. Sci. Paris, Ser. I 349 (2011) 111–114Contents lists available at ScienceDirectC. R. Acad. Sci. Paris, Ser. Iwww.sciencedirect.comMathematical Problems in Mechanics/Calculus of VariationsHomogenization of hexagonal latticesHomogénéisation d’un réseau hexagonalHervé Le Dret a, Annie Raoult ba UPMC Univ Paris 06, UMR 7598 LJLL, Paris, F-75005 Franceb Laboratoire MAP5, Université Paris Descartes, UMR CNRS 8145, 45, rue des Saints Pères, 75270 Paris cedex 06, Francea r t i c l e i n f o a b s t r a c tArticle history:Received 8 November 2010Accepted 14 December 2010Available online 23 December 2010Presented by Philippe G. CiarletWe characterize the macroscopic effective behavior of a graphene sheet modeled by ahexagonal lattice of elastic bars, using Γ -convergence.© 2010 Académie des sciences. Published by Elsevier Masson SAS. All rights reserved.r é s u m éNous identifions le comportement macroscopique d’une feuille de graphène modélisée parun réseau hexagonal de barres élastiques. Nous utilisons pour cela les techniques de laΓ -convergence.© 2010 Académie des sciences. Published by Elsevier Masson SAS. All rights reserved.Version française abrégéeNous considérons une feuille de graphène modélisée par un réseau hexagonal de barres élastiques, voir [3]. Nous cher-chons à en caractériser le comportement mécanique macroscopique. Il s’agit d’un réseau complexe et à notre connaissance,les travaux récents sur le passage à la limite discret vers continu, voir [1,2], ne s’appliquent pas.Le principe de notre analyse consiste à associer aux inconnues discrètes d’une part des fonctions affines par morceauxdécrivant le comportement des atomes qui sont situés sur les sommets d’une triangulation, et d’autre part, des fonctionsconstantes par morceaux permettant de prendre en compte les atomes restants. L’énergie (1) du système discret se réécritsous la forme (2), qui permet d’utiliser les techniques du calcul des variations et de l’homogénéisation. Il ne s’agit toutefoispas d’un problème d’homogénéisation classique, car il fait intervenir deux champs indépendants, des densités qui s’annulentsur la moitié des triangles et des énergies valant +∞ en dehors d’espaces de dimension finie.Nous effectuons une analyse de Γ -convergence quand la longueur des barres tend vers 0. L’énergie limite obtenue prendla forme d’une énergie élastique non linéaire membranaire, voir Proposition 3.3 de la version anglaise, dont la densitéhomogénéisée est donnée par les formules (4) et (6). Pour la démonstration, nous adaptons certains arguments de l’homo-généisation périodique [7,8], au contexte du passage du discret au continu. En particulier, nous introduisons une versiondiscrète de la méthode de slicing de De Giorgi. La densité homogénéisée satisfait le principe d’indifférence matérielle. Lesrotations d’angle multiple entier de π3 appartiennent à son groupe de symétrie matérielle.Pour conclure, nous faisons quelques remarques sur la règle de Cauchy–Born dans le contexte du modèle de feuille degraphène considéré ici. Nous montrons qu’une feuille de graphène n’obéit pas à cette règle au niveau discret. Cependant, lemodèle limite satisfait une version affaiblie de la règle de Cauchy–Born.E-mail addresses: herve.le_dret@upmc.fr (H. Le Dret), annie.raoult@parisdescartes.fr (A. Raoult).1631-073X/$ – see front matter © 2010 Académie des sciences. Published by Elsevier Masson SAS. All rights reserved.doi:10.1016/j.crma.2010.12.012112 H. Le Dret, A. Raoult / C. R. Acad. Sci. Paris, Ser. I 349 (2011) 111–114Fig. 1. •: type 1 nodes, ◦: type 2 nodes, bar types shown on rightmost triangle.1. IntroductionWe consider a graphene sheet modeled by a hexagonal network of elastic bars, see [3]. We are interested in derivingan equivalent continuum mechanics model for the deformations of the sheet by means of a homogenization procedurewhen the rest lengths of the bars go to 0, using Γ -convergence techniques in order to obtain rigorous convergence results.There is a comprehensive body of work on the homogenization of discrete networks, see for instance [1,2]. However, to ourknowledge, the present problem remained to be solved.Let us describe the reference configuration of the sheet. Let (e1, e2) be an orthonormal basis in R2, and s =√3e1,t = √3( 12 e1 +√32 e2) and p = 13 (s + t). In our description, the network is comprised of two types of nodes, the type 1 nodesthat sit at points is + jt with (i, j) ∈ Z2, and the type 2 nodes that sit at points is + jt + p again with (i, j) ∈ Z2, see Fig. 1.The actual graphene sheet is a scaled version of the basic hexagonal network by a small factor ε, furthermore cropped to fitinside a given open set ω of R2. Thus, when at rest, the nodes are located at ε(is + jt) or ε(is + jt + p) depending on theirtype and for those integers i and j that are such that the corresponding points fall within ω. This description using twotypes of nodes is consistent with the standard description of complex lattices, which are superpositions of a finite numberof shifted simple Bravais lattices, see [4,5].For simplicity, we will ignore external forces. It is fairly straightforward to add corresponding terms in the energy. Weimpose a boundary condition of place on part of ∂ω.The internal energy of the sheet is assumed to only derive from elastic bars that join type 1 nodes to their nearestneighboring type 2 nodes. There are thus three types of bars: type 1 bars parallel to s–p, type 2 bars parallel to t–p, andtype 3 bars parallel to p. Under loading and/or boundary conditions, the sheet deforms in R3 and stores an elastic energythat is assumed to be of the formEε =∑all barsk((deformed bar length) − ε)2, (1)where k is a stiffness coefficient. This specific form is not crucial for the later convergence result, but we keep it to fix ideas.Deformed bar lengths are computed using the relative displacements of nodes. The deformation of type α nodes is givenby a function χεα : Z2 → R3, which is actually only defined on the subset of Z2 corresponding to nodes in ω. A typicaldeformed bar length thus assumes the form |χε1 (i, j) − χε2 (i′, j′)| where | · | denotes the Euclidean norm in R3 and theinteger pairs (i, j) and (i′, j′) are related to one another in a way that is simple to work out depending on the type of thebar. The deformed configuration of the sheet minimizes the total energy among all possible configurations.2. Continuous formulationIn order to derive a limit continuous model, we replace the discrete unknowns χεα by continuous counterparts whilekeeping the same energy, see [1–3]. First, we define a piecewise affine function ϕε by declaring that ϕε(ε(is+ jt)) = χε1 (i, j).Since type 1 nodes are located at the vertices of a triangulation T ε , this determines ϕε . Type 2 nodes are taken into accountvia a piecewise constant deviation vector defined by γ ε(x) = χε2 (i, j) − χε1 (i, j) in the triangles that contain a type 2 node,which we call full triangles, and γ ε(x) = 0 in those that do not contain a type 2 node, which we call empty triangles. It isthen a simple matter to express the deformed lengths 1, 2 and 3 of each type of bar with these new unknowns as1 =∣∣ε∂sϕε(x) − γ ε(x)∣∣, 2 = ∣∣ε∂tϕε(x) − γ ε(x)∣∣ and 3 = ∣∣γ ε(x)∣∣,where ∂sϕ = Dϕ(s), ∂tϕ = Dϕ(t) are the directional derivatives in directions s and t , and x is any point in a full trianglecontaining these particular three bars. This asymmetric treatment of nodes and bars of different types is in a sense artificialfrom the point of view of mechanics. It is however, as already noticed before, in keeping with the standard description ofcomplex lattices.Let Y = T f ∪ Te be the parallelogram obtained from the reference full triangle T f whose vertices are (0,0), (1,0) and(0,1) in the oblique coordinate system based on s and t , and the reference empty triangle Te whose vertices are (1,0),H. Le Dret, A. Raoult / C. R. Acad. Sci. Paris, Ser. I 349 (2011) 111–114 113(1,1) and (0,1), and L be the Z-lattice generated by s and t . From now on, all computations will be performed in thisoblique coordinate system. Let L(R2;R3) be the space of linear mappings from R2 into R3. Because all the interveningquantities are piecewise constant, we can rewrite the elastic energy as an integral Eε = Iε(ϕε,γ ε) whereIε(ψ, δ) =∫ωW ε(ε−1x, Dψ(x), δ(x))dx, (2)and W ε : R2 × L(R2;R3) × R3 → R is defined byW ε(y, g, τ ) ={2k[(|g(s) − ε−1τ | − 1)2 + (|g(t) − ε−1τ | − 1)2 + (ε−1|τ | − 1)2] if y ∈ T f + L,0 if y ∈ Te + L.We denote by A(ε) (resp. C(ε)) the space of functions that are piecewise affine continuous (resp. piecewise constant andzero in empty triangles) on T ε and satisfy the boundary condition of place. We extend the sheet energy functional to thespace H = L2(ω;R3) × L2(ω;R3) by letting Iε(ψ, δ) = +∞ whenever ψ /∈ A(ε) or δ /∈ C(ε).It is clear that we have rephrased the equilibrium of the sheet as a problem in the calculus of variations: Find(ϕε,γ ε) ∈ H such thatIε(ϕε,γ ε) = inf(ψ,δ)∈H Iε(ψ, δ). (3)Our objective now is to let ε → 0 and find a limit problem that describes the asymptotic behavior of the continuoussheet deformation ϕε and deviation vector γ ε . This is a periodic nonlinear variational homogenization problem, see [7,8],with several differences compared with the classical case: the energy functionals are +∞ outside of finite-dimensionalsubspaces of H that depend on ε and the densities are not coercive on the unit cell Y .3. Gamma-limit of the energiesFor our purposes here, we are interested in computing (Γ - limε→0 Iε) in the strong topology of H . It should be notedthat the specific form of the energy used earlier plays no role in the ensuing analysis, and the Γ -convergence result holdstrue for more general energies, defined for instance on W 1,p(ω;R3) with p < +∞. We start with an a priori bound. LetH1Γ (ω;R3) denote the space of H1 deformations that satisfy the boundary condition of place.Proposition 3.1. For all (ψ, δ) ∈ H and given any sequence (ψε, δε) → (ψ, δ) in H such that Iε(ψε, δε)  M for some M independentof ε, there exists C independent of ε such that∥∥ψε∥∥H1(ω;R3)  C and∥∥δε∥∥L2(ω;R3) → 0.Consequently, for any Γ -convergent subsequence, we have (Γ - limε→0 Iε)(ψ, δ) = +∞ if ψ /∈ H1Γ (ω;R3) or δ 	= 0.Note that, even though W ε vanishes on empty triangles, the fact that (ψε, δε) belongs to A(ε) × C(ε) enables us torecover uniform coercivity.Let us now introduce the functionW0(y, g) = infτ∈R3W ε(y, g, τ ), (4)which is defined on R2 × L(R2;R3), no longer depends on ε and vanishes for y ∈ Te + L, and the functionalIε0(ψ) =∫ωW0(ε−1x, Dψ(x))dx if ψ ∈ A(ε), Iε0(ψ) = +∞ if ψ ∈ L2(ω;R3) \ A(ε). (5)The following proposition gives the connection between the functionals Iε and Iε0 .Proposition 3.2. For any Γ -convergent subsequence, we have, for all ψ ∈ H1Γ (ω;R3),(Γ - limε→0 Iε)(ψ,0) =(Γ - limε→0 Iε0)(ψ).We define a homogenized energy density on L(R2;R3) byWhom(g) = infk∈N{1k2(infθ∈A(kY )∫kYW0(y, g + Dθ(y)) dy)}, (6)where A(kY ) denotes the set of piecewise affine functions on the mesh defined on kY that vanish on ∂(kY ).114 H. Le Dret, A. Raoult / C. R. Acad. Sci. Paris, Ser. I 349 (2011) 111–114Proposition 3.3. For any Γ -convergent subsequence, we have, for all ψ ∈ H1Γ (ω;R3),(Γ - limε→0 Iε0)(ψ) =∫ωWhom(Dψ(x))dx.The proof of Proposition 3.3 takes arguments from [7,8], adapted to the discrete to continuous context. In particular,it uses a discrete version of De Giorgi’s slicing method. The uniqueness of the Γ -limit and the compactness of functionalsequences with respect to Γ -convergence imply that the whole sequence Γ converges to the above limit. In addition, thebounds given in Proposition 3.1 show that the usual De Giorgi argument applies, i.e., the limit points of the sequence ofminimizers in the weak topology of H1 minimize the Γ -limit functional. We have thus obtained an asymptotic descriptionof the mechanical behavior of the graphene sheet.It is easy to see that the homogenized energy density is frame-indifferent and that its material symmetry group containsall rotations of angle an integer multiple of π3 , see also [9]. Finally, the homogenized density vanishes on deformations thatcompress the sheet, which is normal behavior for a membrane model, see Section 4.4. Remarks on the Cauchy–Born ruleThe Cauchy–Born rule postulates that if a homogeneous deformation is imposed on the boundary of a monoatomiccrystal, all the atoms move according to the same homogeneous deformation, see [4,6]. In the case of a graphene sheet asmodeled here, the Cauchy–Born rule fails. Indeed, assume that g is a linear mapping from R2 to R3 such that the radius ofthe circumcircle to the triangle g(εT f ) is strictly less than ε. Then, there are two points on the normal to g(εT f ) passingthrough its circumcenter that are at distance ε from the three vertices. Choosing τ corresponding to one of these pointsyields W ε(y, g, τ ) = 0, whereas W ε(y, g, g(εp)) > 0 for such a g . This corresponds to an energy-minimizing, out-of-planeand nonhomogeneous deformation of type 2 nodes. Even if deformations are restricted to the plane, it is possible to fold ahexagonal network on itself with zero energy, whereas a homogeneous compression has a strictly positive energy cost.If we look now at the limit problem, then by the abstract properties of Γ -convergence, we know that the homogenizeddensity Whom is quasiconvex, that is to say∫ωWhom(g + Dθ(x)) dx ∫ωWhom(g)dxfor all g ∈ L(R2;R3) and θ Lipschitz vanishing on the boundary of ω. Hence, the homogeneous deformation φ(x) = g(x) inω is an energy minimizer for the boundary condition φ(x) = g(x) on ∂ω. This is an expression of the Cauchy–Born rule interms of energy minimization. Hence, it can be argued that a slightly weaker version of the Cauchy–Born rule is recoveredin the homogenization limit, in the sense that homogeneous deformations minimize the limit energy for homogeneousboundary data, but that they may not be the only minimizers.References[1] R. Alicandro, M. Cicalese, A general integral representation result for continuum limits of discrete energies with superlinear growth, SIAM J. Math.Anal. 36 (1) (2004) 1–37.[2] A. Braides, M.S. Gelli, From Discrete to Continuum: A Variational Approach, Lecture Notes, SISSA, Trieste, 2000.[3] D. Caillerie, A. Mourad, A. Raoult, Discrete homogenization in graphene sheet modeling, J. Elasticity 84 (2006) 33–68.[4] W. E, P. Ming, Cauchy–Born rule and the stability of crystalline solids: Static problems, Arch. Rational Mech. Anal. 183 (2007) 241–297.[5] P. Engel, Geometric Crystallography, Kluwer Academic Publishers, 1986.[6] G. Friesecke, F. Theil, Validity and failure of the Cauchy–Born hypothesis in a two-dimensional mass-spring lattice, J. Nonlinear Sci. 12 (2002) 445–478.[7] P. Marcellini, Periodic solutions and homogenization of nonlinear variational problems, Ann. Mat. Pura Appl. IV (CXVII) (1978) 139–152.[8] S. Müller, Homogenization of nonconvex integral functionals and cellular elastic materials, Arch. Rational Mech. Anal. 99 (3) (1987) 189–212.[9] A. Raoult, D. Caillerie, A. Mourad, Elastic lattices: Equilibrium, invariant laws and homogenization, Ann. Univ. Ferrara 54 (2008) 297–318.

Homogenization of hexagonal lattices

Le Dret, Hervé

Raoult, Annie

Hal-Diderot

HAL Descartes

Numérisation de Documents Anciens Mathématiques

https://comptes-rendus.academie-sciences.fr/mathematique/item/10.1016/j.crma.2010.12.012.pdf

Homogenization of hexagonal lattices

Abstract

Similar works

Full text

Available Versions

Comptes Rendus Mathématique

Hal-Diderot

HAL Descartes

Numérisation de Documents Anciens Mathématiques