High contrast homogenisation in nonlinear elasticity under small loads

Abstract

We study the homogenisation of geometrically nonlinear elastic composites with high contrast. The composites we analyse consist of a perforated matrix material, which we call the “stiff” material, and a “soft” material that fills the remaining pores. We assume that the pores are of size $0\mathrm{<}\mathit{\epsilon }\ll 1$ and are periodically distributed with period ε. We also assume that the stiffness of the soft material degenerates with rate ${\mathit{\epsilon }}^{2\mathit{\gamma }}$ , $\mathit{\gamma }\mathrm{>}0$ , so that the contrast between the two materials becomes infinite as $\mathit{\epsilon }\downarrow 0$ . We study the homogenisation limit $\mathit{\epsilon }\downarrow 0$ in a low energy regime, where the displacement of the stiff component is infinitesimally small. We derive an effective two-scale model, which, depending on the scaling of the energy, is either a quadratic functional or a partially quadratic functional that still allows for large strains in the soft inclusions. In the latter case, averaging out the small scale-term justifies a single-scale model for high-contrast materials, which features a non-linear and non-monotone effect describing a coupling between microscopic and the effective macroscopic displacements.

1. Introduction

We consider a geometrically nonlinear elastic composite material that consists of a “stiff” matrix material and periodically distributed pores filled by a “soft” material: for $\mathit{\epsilon }\mathrm{>}0$ and a fixed scaling parameter $\mathit{\gamma }\mathrm{>}0$ we consider the energy functional of non-linear elasticity $$\begin{array}{}\text{(1)}& {\mathcal{I}}_{\mathit{\epsilon }}\mathrm{(}\mathit{u}\mathrm{)}:={{\displaystyle \int }}_{\mathrm{\Omega }}\mathrm{(}{\mathit{\epsilon }}^{2\mathit{\gamma }}{\mathit{W}}^0\mathrm{(}\nabla \mathit{u}\mathrm{)}{\mathit{\chi }}_{\mathit{\epsilon }}+{\mathit{W}}^1\mathrm{(}\nabla \mathit{u}\mathrm{)}\mathrm{(}1-{\mathit{\chi }}_{\mathit{\epsilon }}\mathrm{)}\mathrm{)}\mathit{d}\mathit{x}-{{\displaystyle \int }}_{\mathrm{\Omega }}{\mathit{f}}_{\mathit{\epsilon }}·\mathit{u}\mathit{d}\mathit{x}\mathrm{,}\mathit{u}\in {\mathit{H}}^1\mathrm{(}\mathrm{\Omega }\mathrm{)}.\end{array}$$ Here Ω denotes a Lipschitz domain in ${\mathbb{R}}^{\mathit{d}}$ (the reference domain of the elastic body) and $\mathit{u}:\mathrm{\Omega }\to {\mathbb{R}}^{\mathit{d}}$ is a deformation satisfying clamped boundary conditions: $\mathit{u}\mathrm{(}\mathit{x}\mathrm{)}=\mathit{x}$ on $\partial \mathrm{\Omega }$ . We denote by ${\mathit{f}}_{\mathit{\epsilon }}:\mathrm{\Omega }\to {\mathbb{R}}^{\mathit{d}}$ the density of the applied body forces, ${\mathit{W}}^0$ and ${\mathit{W}}^1$ are frame-indifferent, non-degenerate energy densities (see Section 2 below for the precise assumptions), and ${\mathit{\chi }}_{\mathit{\epsilon }}$ denotes the indicator function of the pores, i.e. of the domain occupied by the “soft” material component. As will be made precise in Section 2, we assume that the pores are of size ε and are periodically distributed in the interior of Ω with period ε. As can be seen from (1), in the homogenisation limit $\mathit{\epsilon }\downarrow 0$ the stiffness of the “soft” material degenerates with rate ${\mathit{\epsilon }}^{2\mathit{\gamma }}$ ( $\mathit{\gamma }\mathrm{>}0$ ), while the stiffness of the “stiff” material remains unchanged. Hence, the contrast between the soft material (occupying the pores) and the stiff material (occupying the perforated matrix) becomes infinite in the limit $\mathit{\epsilon }\downarrow 0$ . We therefore refer to the corresponding limit procedure as high-contrast homogenisation. Our goal is to identify the effective behaviour of the minimisation problem associated with ${\mathcal{I}}_{\mathit{\epsilon }}$ by studying its limit under a proper rescaling.

Summary and discussion of our result. To illustrate our result, here in the introduction we restrict ourselves to the special case $\mathit{\gamma }=1$ . If we assume that the density of the body forces is small in magnitude, in the sense that ${\mathit{f}}_{\mathit{\epsilon }}={\mathit{\epsilon }}^{\mathit{\alpha }}\mathit{f}$ for some $\mathit{\alpha }⩾1$ , and has vanishing first moment, i.e. ${\int }_{\mathrm{\Omega }}\mathit{f}\mathrm{(}\mathit{x}\mathrm{)}·\mathit{x}\mathit{d}\mathit{x}=0$ , then (1) can be expressed as $$\begin{array}{rcl}{\mathit{I}}_{\mathit{\epsilon }}^{\mathit{\alpha }}\mathrm{(}\mathit{\phi }\mathrm{)}& :=& {\displaystyle \frac{1}{{\mathit{\epsilon }}^{2\mathit{\alpha }}}}{\mathcal{I}}_{\mathit{\epsilon }}\mathrm{(}\mathit{u}\mathrm{)}\\ \text{(2)}& & =& {{\displaystyle \int }}_{\mathrm{\Omega }}\mathrm{\left(}{\displaystyle \frac{1}{{\mathit{\epsilon }}^{2\mathrm{(}\mathit{\alpha }-1\mathrm{)}}}}{\mathit{W}}^0\mathrm{\right(}\mathit{I}+{\mathit{\epsilon }}^{\mathit{\alpha }}\nabla \mathit{\phi }\mathrm{)}{\mathit{\chi }}_{\mathit{\epsilon }}+{\displaystyle \frac{1}{{\mathit{\epsilon }}^{2\mathit{\alpha }}}}{\mathit{W}}^1\mathrm{(}\mathit{I}+{\mathit{\epsilon }}^{\mathit{\alpha }}\nabla \mathit{\phi }\mathrm{)}\mathrm{(}1-{\mathit{\chi }}_{\mathit{\epsilon }}\mathrm{)}\mathrm{)}\mathit{d}\mathit{x}-{{\displaystyle \int }}_{\mathrm{\Omega }}\mathit{f}·\mathit{\phi }\mathit{d}\mathit{x}\mathrm{,}\end{array}$$ where $$\begin{array}{}\text{(3)}& \mathit{\phi }\mathrm{(}\mathit{x}\mathrm{)}={\displaystyle \frac{\mathit{u}\mathrm{(}\mathit{x}\mathrm{)}-\mathit{x}}{{\mathit{\epsilon }}^{\mathit{\alpha }}}}\mathrm{,}\mathit{x}\in \mathrm{\Omega }\mathrm{,}\mathit{\phi }\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega };{\mathbb{R}}^{\mathit{d}}\mathrm{)}\mathrm{,}\end{array}$$ denotes the scaled displacement, and I stands for the identity matrix in ${\mathbb{R}}^{\mathit{d}\times \mathit{d}}$ . In this paper we analyse the asymptotics of the minimisation problem associated with ${\mathit{I}}_{\mathit{\epsilon }}^{\mathit{\alpha }}$ in the limit $\mathit{\epsilon }\downarrow 0$ by appealing to the concept of Γ-convergence. The latter goes back to De Giorgi (e.g. see [20] for a standard reference). In a metric setting it is defined as follows:

Definition 1 (Γ-convergence).

Let $\mathrm{(}\mathit{X}\mathrm{,}\mathit{d}\mathrm{)}$ denote a metric space. A sequence of functionals ${\mathit{I}}_{\mathit{\epsilon }}:\mathit{X}\to [-\infty \mathrm{,}\infty ]$ Γ-converges to a functional $\mathit{I}:\mathit{X}\to [-\infty \mathrm{,}\infty ]$ , if

(lower bound). For every ${\mathit{x}}_0\in \mathit{X}$ and every ${\mathit{x}}_{\mathit{\epsilon }}\to {\mathit{x}}_0$ in X we have ${lim\hspace{0.17em}inf}_{\mathit{\epsilon }\to 0}{\mathit{I}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{x}}_{\mathit{\epsilon }}\mathrm{)}⩾{\mathit{I}}_0\mathrm{(}{\mathit{x}}_0\mathrm{)}$ .

(recovery sequence). For every ${\mathit{x}}_0\in \mathit{X}$ there exists a sequence ${\mathit{x}}_{\mathit{\epsilon }}\to {\mathit{x}}_0$ in X such that ${lim}_{\mathit{\epsilon }\to 0}{\mathit{I}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{x}}_{\mathit{\epsilon }}\mathrm{)}={\mathit{I}}_0\mathrm{(}{\mathit{x}}_0\mathrm{)}$ .

In that case we call ${\mathit{I}}_0$ the Γ-limit of the sequence ${\mathit{I}}_{\mathit{\epsilon }}$ .

A fundamental property of Γ-convergence is the following fact: If a sequence of functionals Γ-converges and the functionals are equicoercive, then the associated sequence of minima (resp. minimisers) converge (up to a subsequence) to the minimum (resp. a minimiser) of the Γ-limit, and any minimiser of the Γ-limit can be obtained as a limit of a minimizing sequence of the original functionals. Thanks to this property, Γ-convergence is especially useful for the study of the asymptotics of parametrised minimisation problems. The Γ-limit, if it exists, is unique; yet, the question whether a sequence of functionals Γ-converges or not, and the form of the Γ-limit depend on the topology of X. In particular, a sequence of functionals is more likely to Γ-converge in a stronger topology, while it is more likely to be equicoercive in a weaker topology. Therefore, it is natural to consider the strongest notion of convergence on X for which the functionals remain equicoercive. In the situation we are interested in, namely the asymptotics of the functionals ${\mathit{I}}_{\mathit{\epsilon }}^{\mathit{\alpha }}$ defined on $\mathit{X}={\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ , it turns out that due to the presence of high-contrast, a natural and appropriate notion of convergence on ${\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ is a variant of two-scale convergence that we discuss next. Let us first recall the standard notion of two-scale convergence from [32] and [1]:

Definition 2.

We say that a sequence ${\mathit{f}}_{\mathit{\epsilon }}\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}$ weakly two-scale converges to $\mathit{f}\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\times \mathit{Y}\mathrm{)}$ if the sequence ${\mathit{f}}_{\mathit{\epsilon }}$ is bounded and $$\begin{array}{c}\underset{\mathit{\epsilon }\downarrow 0}{lim}{{\displaystyle \int }}_{\mathrm{\Omega }}{\mathit{f}}_{\mathit{\epsilon }}\mathrm{(}\mathit{x}\mathrm{)}\mathit{\phi }\mathrm{(}\mathit{x}\mathrm{,}\mathit{x}\mathrm{/}\mathit{\epsilon }\mathrm{)}\mathit{d}\mathit{x}={\iint }_{\mathrm{\Omega }\times \mathit{Y}}\mathit{f}\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}\mathit{\phi }\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}\mathit{d}\mathit{x}\mathit{d}\mathit{y}\end{array}$$ for all $\mathit{\phi }\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{,}{\mathit{C}}_{\mathrm{\#}}\mathrm{(}\mathit{Y}\mathrm{)}\mathrm{)}$ , where ${\mathit{C}}_{\mathrm{\#}}\mathrm{(}\mathit{Y}\mathrm{)}$ denotes the Banach space of continuous, $\mathit{Y}:={\mathrm{(}0\mathrm{,}1\mathrm{)}}^{\mathit{d}}$ -periodic functions on ${\mathbb{R}}^{\mathit{d}}$ . We say that a sequence ${\mathit{f}}_{\mathit{\epsilon }}\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}$ strongly two-scale converges to $\mathit{f}\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\times \mathit{Y}\mathrm{)}$ if the sequence ${\mathit{f}}_{\mathit{\epsilon }}$ weakly two-scale converges to f and one has $\Vert {\mathit{f}}_{\mathit{\epsilon }}{\Vert }_{{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}}\to \Vert \mathit{f}{\Vert }_{{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\times \mathit{Y}\mathrm{)}}$ as $\mathit{\epsilon }\downarrow 0$ . For vector-valued functions two-scale convergence is defined component-wise.

In the study of the asymptotics of ${\mathit{I}}_{\mathit{\epsilon }}^{\mathit{\alpha }}$ we work with a variant of this definition that is tailor-made to capture the effects of high-contrast. For $\mathit{\gamma }=1$ it can be summarised as follows (for details and the general case see Section 3.1): Given a sequence of displacements ${\mathit{\phi }}_{\mathit{\epsilon }}$ in ${\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ we consider the unique decomposition ${\mathit{\phi }}_{\mathit{\epsilon }}={\mathit{g}}_{\mathit{\epsilon }}^0+{\mathit{g}}_{\mathit{\epsilon }}^1$ into a function ${\mathit{g}}_{\mathit{\epsilon }}^0\in {\mathit{H}}_0^1\mathrm{(}{\mathrm{\Omega }}_{\mathit{\epsilon }}^0\mathrm{)}$ and a contribution ${\mathit{g}}_{\mathit{\epsilon }}^1\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ that is harmonic in ${\mathrm{\Omega }}_{\mathit{\epsilon }}^0$ . We then write ${\mathit{\phi }}_{\mathit{\epsilon }}\stackrel{2}{\mathrm{\rightharpoonup }}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}$ and say that ${\mathit{\phi }}_{\mathit{\epsilon }}$ converges to a pair $\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}$ with ${\mathit{g}}^0\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega };{\mathit{H}}_0^1\mathrm{(}{\mathit{Y}}_0\mathrm{)}\mathrm{)}$ and ${\mathit{g}}^1\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ , if ${\mathit{g}}_{\mathit{\epsilon }}^1\mathrm{\rightharpoonup }{\mathit{g}}^1$ weakly in ${\mathit{H}}^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ , and ${\mathit{g}}_{\mathit{\epsilon }}^0\stackrel{2}{\mathrm{\rightharpoonup }}{\mathit{g}}^0$ , $\mathit{\epsilon }\nabla {\mathit{g}}_{\mathit{\epsilon }}^0\stackrel{2}{\mathrm{\rightharpoonup }}{\nabla }_{\mathit{y}}{\mathit{g}}^0$ weakly two-scale. The component ${\mathit{g}}^1$ of a limit pair $\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}$ describes the (scaled) macroscopic displacement of the body, and ${\mathit{g}}^0$ is a two-scale function describing the (scaled) microscopic displacement on the pores relative to the deformed matrix material. As a main result (see Theorem 1 and Theorem 3) we prove (in fact in a slightly more general situation) that ${\mathit{I}}_{\mathit{\epsilon }}^{\mathit{\alpha }}$ Γ-converges with respect to the type of convergence introduced above. It turns out that two different regimes emerge for $\mathit{\alpha }\mathrm{>}1$ and $\mathit{\alpha }=1$ .

In the small strain regime ( $\mathit{\alpha }\mathrm{>}1$ ), the strain ${\mathit{\epsilon }}^{\mathit{\alpha }}\nabla \mathit{\phi }$ becomes infinitesimally small in the entire domain Ω, and the limit behaviour is expressed by a linearised, two-scale energy ${\mathit{I}}_{\text{small}}:{\mathit{L}}^2\mathrm{(}\mathrm{\Omega };{\mathit{H}}_0^1\mathrm{(}{\mathit{Y}}_0\mathrm{)}\mathrm{)}\times {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega };{\mathbb{R}}^{\mathit{d}}\mathrm{)}\to \mathbb{R}$ , $$\begin{array}{}\text{(4)}& {\mathit{I}}_{\text{small}}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}& :=& {{\displaystyle \int }}_{\mathrm{\Omega }\times \mathit{Y}}{\mathit{Q}}^0\mathrm{(}{\nabla }_{\mathit{y}}{\mathit{g}}^0\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}\mathrm{)}\mathit{d}\mathit{y}+{\mathit{Q}}_{hom}^1\mathrm{(}\nabla {\mathit{g}}^1\mathrm{(}\mathit{x}\mathrm{)}\mathrm{)}\mathit{d}\mathit{x}\text{(5)}& & & -{{\displaystyle \int }}_{\mathrm{\Omega }}\mathrm{(}{{\displaystyle \int }}_{{\mathit{Y}}_0}{\mathit{g}}^0\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}\mathit{d}\mathit{y}+{\mathit{g}}^1\mathrm{(}\mathit{x}\mathrm{)}\mathrm{)}·\mathit{f}\mathrm{(}\mathit{x}\mathrm{)}\mathit{d}\mathit{x}.\end{array}$$ Here ${\mathit{Q}}^0$ and ${\mathit{Q}}^1$ are the quadratic forms of the quadratic expansions of ${\mathit{W}}^0$ and ${\mathit{W}}^1$ at the identity, and ${\mathit{Q}}_{hom}^1$ denotes the homogenised energy density obtained from ${\mathit{Q}}^1$ , see (12) and (23) for details. The functional ${\mathit{I}}_{\text{small}}$ is the two-scale Γ-limit of the sequence ${\mathit{I}}_{\mathit{\epsilon }}^{\mathit{\alpha }}$ in the sense that (cf. Theorem 1):

(lower bound). For every $\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}$ and every ${\mathit{\phi }}_{\mathit{\epsilon }}\stackrel{2}{\mathrm{\rightharpoonup }}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}$ we have ${lim\hspace{0.17em}inf}_{\mathit{\epsilon }\downarrow 0}{\mathit{I}}_{\mathit{\epsilon }}^{\mathit{\alpha }}\mathrm{(}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}⩾{\mathit{I}}_{\text{small}}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}$ .

(recovery sequence). For every $\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}$ there exists a sequence ${\mathit{\phi }}_{\mathit{\epsilon }}\stackrel{2}{\mathrm{\rightharpoonup }}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}$ such that ${lim}_{\mathit{\epsilon }\downarrow 0}{\mathit{I}}_{\mathit{\epsilon }}^{\mathit{\alpha }}\mathrm{(}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}⩾{\mathit{I}}_{\text{small}}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}$ .

In addition, we prove that the functionals ${\mathrm{(}{\mathit{I}}_{\mathit{\epsilon }}^{\mathit{\alpha }}\mathrm{)}}_{\mathit{\epsilon }\mathrm{>}0}$ are equicoercive and deduce convergence of the associated minimisation problems (see Theorem 1 and Proposition 1). In Theorem 2 we establish a two-scale expansion showing that if ${\mathit{\phi }}^{\mathit{\epsilon }}$ is an (almost) minimiser of ${\mathit{I}}_{\mathit{\epsilon }}^{\mathit{\alpha }}$ , then $$\begin{array}{}\text{(6)}& {\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{(}\mathit{x}\mathrm{)}\approx {\mathit{g}}_{\ast }^1\mathrm{(}\mathit{x}\mathrm{)}+\mathit{\epsilon }\mathit{\psi }\mathrm{(}\mathit{x}\mathrm{,}\mathit{x}\mathrm{/}\mathit{\epsilon }\mathrm{)}+{\mathit{g}}_{\ast }^0\mathrm{(}\mathit{x}\mathrm{,}\mathit{x}\mathrm{/}\mathit{\epsilon }\mathrm{)}\mathrm{,}\end{array}$$ where $\mathrm{(}{\mathit{g}}_{\ast }^0\mathrm{,}{\mathit{g}}_{\ast }^1\mathrm{)}$ is a minimiser of ${\mathit{I}}_{\text{small}}$ , and ψ denotes a corrector function that only depends on ${\mathit{g}}_{\ast }^1$ . Finally, we illustrate that by averaging out the fast variable y the limit ${\mathit{I}}_{\text{small}}$ can be further simplified. In fact, Proposition 2 shows that ${\mathit{I}}_{\mathit{\epsilon }}^{\mathit{\alpha }}$ Γ-converges (with respect to weak convergence in ${\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}$ ) to the functional ${\overline{\mathit{I}}}_{\text{small}}:{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}\to \mathbb{R}\cup \{+\infty \}$ given by $$\begin{array}{rcl}{\overline{\mathit{I}}}_{\text{small}}\mathrm{(}\mathit{\phi }\mathrm{)}& :=& min\left\{{{\displaystyle \int }}_{\mathrm{\Omega }}{\mathit{Q}}_{hom}^1\mathrm{\right(}\nabla {\mathit{g}}^1\mathrm{)}+\bar{\mathit{Q}}\mathrm{(}\bar{\mathit{G}}\mathrm{)}:{\mathit{g}}^1\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}\mathrm{,}\bar{\mathit{G}}\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}\text{ with }{\mathit{g}}^1+\bar{\mathit{G}}=\mathit{\phi }\}\\ \text{(7)}& & & -{{\displaystyle \int }}_{\mathrm{\Omega }}\mathit{f}·\mathit{\phi }\mathit{d}\mathit{x}\mathrm{,}\end{array}$$ with a (positive definite) quadratic form $\bar{\mathit{Q}}:{\mathbb{R}}^{\mathit{d}}\to [0\mathrm{,}\infty \mathrm{)}$ defined by $$\begin{array}{c}\bar{\mathit{Q}}\mathrm{(}\bar{\mathit{G}}\mathrm{)}:=min\left\{{{\displaystyle \int }}_{{\mathit{Y}}^0}{\mathit{Q}}^0\mathrm{\right(}\nabla {\mathit{g}}^0\mathrm{(}\mathit{y}\mathrm{)}\mathrm{)}\mathit{d}\mathit{y}:{\mathit{g}}^0\in {\mathit{H}}_0^1\mathrm{(}{\mathit{Y}}^0\mathrm{)}\mathrm{,}{{\displaystyle \int }}_{{\mathit{Y}}^0}{\mathit{g}}^0\mathit{d}\mathit{y}=\bar{\mathit{G}}\}.\end{array}$$ $\bar{\mathit{Q}}$ captures the influence of the pores (and their geometry) on the effective behavior. The minimiser ${\mathit{\phi }}_{\ast }\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}$ to ${\overline{\mathit{I}}}_{\text{small}}$ takes the form ${\mathit{\phi }}_{\ast }={\mathit{g}}_{\ast }^1+{\bar{\mathit{G}}}_{\ast }$ with ${\bar{\mathit{G}}}_{\ast }:={\int }_{{\mathit{Y}}_0}{\mathit{g}}_{\ast }^0\mathrm{(}·\mathrm{,}\mathit{y}\mathrm{)}\mathit{d}\mathit{y}$ . In view of (6) the field ${\bar{\mathit{G}}}_{\ast }$ can be interpreted as the gap between the macroscopic displacement and the microscopic displacements in the pores.

In the finite strain regime, which corresponds to $\mathit{\alpha }=1$ , the displacement gradient $\mathit{\epsilon }\nabla {\mathit{\phi }}_{\mathit{\epsilon }}$ becomes infinitesimally small only in the stiff component, while large strains still may occur in the soft pores. Therefore, the Γ-limit is a non-convex (partially linearised) functional of the form $$\begin{array}{rcl}{\mathit{I}}_{\text{finite}}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}& :=& {{\displaystyle \int }}_{\mathrm{\Omega }\times \mathit{Y}}\mathcal{Q}{\mathit{W}}^0\mathrm{(}\mathit{I}+{\nabla }_{\mathit{y}}{\mathit{g}}^0\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}\mathrm{)}\mathit{d}\mathit{y}+{\mathit{Q}}_{hom}^1\mathrm{(}\nabla {\mathit{g}}^1\mathrm{(}\mathit{x}\mathrm{)}\mathrm{)}\mathit{d}\mathit{x}\\ & & -{{\displaystyle \int }}_{\mathrm{\Omega }}\mathrm{(}{{\displaystyle \int }}_{{\mathit{Y}}_0}{\mathit{g}}^0\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}\mathit{d}\mathit{y}+{\mathit{g}}^1\mathrm{(}\mathit{x}\mathrm{)}\mathrm{)}·\mathit{f}\mathrm{(}\mathit{x}\mathrm{)}\mathit{d}\mathit{x}\mathrm{,}\end{array}$$ where $\mathcal{Q}{\mathit{W}}^0$ denotes the quasiconvex envelope of ${\mathit{W}}^0$ . Similarly to the small strain regime, one can average out the fast scale y and obtain Γ-convergence of ${\mathit{I}}_{\mathit{\epsilon }}^0$ to the functional ${\overline{\mathit{I}}}_{\text{finite}}:{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}\to \mathbb{R}\cup \{+\infty \}$ given by $$\begin{array}{rcl}\text{(8)}& {\overline{\mathit{I}}}_{\mathrm{finite}}\mathrm{(}\mathit{\phi }\mathrm{)}& :=& min\left\{{{\displaystyle \int }}_{\mathrm{\Omega }}{\mathit{Q}}_{hom}^1\mathrm{\right(}\nabla {\mathit{g}}^1\mathrm{)}+\bar{\mathit{V}}\mathrm{(}\bar{\mathit{G}}\mathrm{)}:{\mathit{g}}^1\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}\mathrm{,}\bar{\mathit{G}}\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}\text{ with }{\mathit{g}}^1+\bar{\mathit{G}}=\mathit{\phi }\}& & -{{\displaystyle \int }}_{\mathrm{\Omega }}\mathit{f}·\mathit{\phi }\mathit{d}\mathit{x}\mathrm{,}\end{array}$$ with non-convex potential $\bar{\mathit{V}}:{\mathbb{R}}^{\mathit{d}}\to [0\mathrm{,}\infty \mathrm{)}$ defined by $$\begin{array}{c}\bar{\mathit{V}}\mathrm{(}\bar{\mathit{G}}\mathrm{)}:=min\left\{{{\displaystyle \int }}_{{\mathit{Y}}^0}\mathcal{Q}{\mathit{W}}^0\mathrm{\right(}\mathit{I}+\nabla {\mathit{g}}^0\mathrm{(}\mathit{y}\mathrm{)}\mathrm{)}\mathit{d}\mathit{y}:{\mathit{g}}^0\in {\mathit{H}}_0^1\mathrm{(}{\mathit{Y}}^0\mathrm{)}\mathrm{,}{{\displaystyle \int }}_{{\mathit{Y}}^0}{\mathit{g}}^0\mathit{d}\mathit{y}=\bar{\mathit{G}}\}\mathrm{,}\end{array}$$ see Remark 2. In contrast to the small strain regime, where $\bar{\mathit{Q}}$ is quadratic, the potential $\bar{\mathit{V}}$ is non-convex and expresses a nonlinear (and non-monotone) coupling between the macroscopic and microscopic displacement.

Connection to acoustic wave propagation in high-contrast materials. The sequence of functionals ${\mathit{\epsilon }}^{-2\mathit{\alpha }}{\mathcal{I}}_{\mathit{\epsilon }}$ , in either of the two regimes described above, occupies an intermediate position between a fully nonlinearly elastic composite and fully linearised models, as $\mathit{\epsilon }\to 0$ . Notably, linear models with high contrast, which are suitable for the description of small displacement fields (that often occur, say, in acoustic wave propagation) already exhibit a coupling between the macroscopic part ${\mathit{g}}^0$ and microscopic part ${\mathit{g}}^1$ of the minimiser of ${\mathit{I}}_{\mathrm{small}}$ , which in our case is obtained as a limit in the small-strain regime $\mathit{\alpha }\mathrm{>}1$ . This can be seen by considering the time-harmonic solitons to the equations of elastodynamics with the elastic part of the energy given by (4), away from the sources of the elastic motion. In this case the function f in (5) (which in our analysis we assume to be independent of the fast variable $\mathit{x}\mathrm{/}\mathit{\epsilon }$ for simplicity, an assumption that can be relaxed with no changes in the proofs needed) has to be replaced by the sum ${\mathit{g}}^0+{\mathit{g}}^1$ , with the integration in (5) carried over ${\mathit{Y}}_0$ and Q at the same time, i.e. the work of the external forces (5) is replaced by the expression for the work of “self-forces” $$\begin{array}{c}-{\mathit{\omega }}^2{{\displaystyle \int }}_{\mathrm{\Omega }}{{\displaystyle \int }}_{{\mathit{Y}}_0}\mathrm{(}{\mathit{g}}^0\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}+{\mathit{g}}^1\mathrm{(}\mathit{x}\mathrm{)}\mathrm{)}·\mathrm{(}{\mathit{g}}^0\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}+{\mathit{g}}^1\mathrm{(}\mathit{x}\mathrm{)}\mathrm{)}\mathit{d}\mathit{y}\mathit{d}\mathit{x}\mathrm{,}\end{array}$$ where ω is the frequency. The solution to the Euler–Lagrange equation for the resulting functional is a coupled system of equations for ${\mathit{g}}^0$ , ${\mathit{g}}^1$ , so that when the equation for ${\mathit{g}}^0$ is solved in terms of ${\mathit{g}}^1$ and substituted into the second equation, it takes the form (away from the sources): $$\begin{array}{}\text{(9)}& {\mathcal{A}}^{\mathrm{hom}}{\mathit{g}}^1={\mathit{\omega }}^2\mathit{\beta }\mathrm{\left(}{\mathit{\omega }}^2\mathrm{\right)}{\mathit{g}}^1\mathrm{,}\end{array}$$ for some non-negative self-adjoint differential operator ${\mathcal{A}}^{\mathrm{hom}}$ and a special nonlinear function β, which takes positive and negative values on alternating intervals of the real axis (leading to “lacunae”, or “band gaps” in the spectrum of the corresponding operator) and is obtained from the spectral decomposition of ${\mathit{g}}^0$ and the subsequent averaging over ${\mathit{Y}}_0$ , see [38]. From this point of view, the non-quadratic finite-strain functional ${\mathit{I}}_{\mathrm{finite}}$ is a “matching”, “partially quadratic”, homogenised model corresponding, e.g., to finite-amplitude, rather than small-amplitude, wave motions that can no longer be treated using a quadratic model such as ${\mathit{I}}_{\mathrm{small}}$ but can still be used in place of models of nonlinear elasticity where the elastic energy terms on both components of the composite (stiff and soft) are non-quadratic.

Methods and previous results. In this paper we appeal to analytic methods that have been developed in the last two decades in the areas of nonlinear elasticity and homogenisation. Among these are the notion of two-scale convergence introduced in [1,32] and periodic unfolding (see [18] and references therein). The convergence statements of our main results are expressed in the language of Γ-convergence (see [20] and references therein). In order to treat the geometric nonlinearity of the considered functional, we make use of the geometric rigidity estimate (see [23]). Since we consider a low energy regime, linearisation and homogenisation take place at the same time. The simultaneous treatment of both effects is inspired by recent works [24,28–31] of the third author, where various problems involving simultaneous homogenisation, linearisation and dimension reduction are studied. The homogenisation of the kind of high-contrast composites that we study is related to the homogenisation for periodically perforated domains (e.g. see [8,33]). For instance, we make use of extensions across the pores. As a side result we prove a version of the geometric rigidity estimate for perforated domains (see Lemma 4 below). We would like to remark that while the present work is one of the few papers, along with [9,16], that treat the fully nonlinear high-contrast case, during the last decade there has been a significant amount of literature devoted to the mathematical analysis of phenomena associated with, or modelled by, a high degree of contrast between the properties of the materials constituting a composite, in the linearised setting. The first contributions in this direction are due to Zhikov [38], and Bouchitté and Felbacq [6], following an earlier paper by Allaire [1] and the collection of papers by Hornung [25] (see also the references therein), where the special role of high-contrast elliptic PDE was pointed out albeit not studied in detail. These works demonstrated that the behaviour of the field variable in such models is of a two-scale type in the homogenisation limit, i.e. the limit model cannot be reduced to a one-scale formulation and fields that depend on the fast variable remain in the effective model. They also noticed that the spectrum of such materials has a band-gap structure, as in (9), and indicated how this fact could be exploited for high-resolution imaging and cloaking. It has since been an adopted approach to the theoretical construction of “negative refraction” media, or more generally “metamaterials”, which is now a hugely popular area of research in physics (see e.g. [34] and references therein). On the analytical side, a number of further works followed, in particular [3,4,7,10,11,13–15,17,26,35], where various consequences of high contrast (or, mathematically speaking, the property of non-uniform ellipticity) in the underlying equations have been explored. Among these are the “non-locality” and “micro-torsion” effects in materials with high-contrast inclusions in the shape of fibres extending in one or more directions, the “partial band-gap” wave propagation due to the high degree of anisotropy of one of the constituent media, and the localisation of energy in high-contrast media with a defect (“photonic crystal fibres”), all of which can be thought of as examples of “non-standard”, or “non-classical”, behaviour in composites, which is not available in the usual moderate-contrast materials. In the present paper we aim to develop further a rigorous high-contrast theory in the context of finite elasticity, where the underlying model is nonlinear.

With this paper we continue the multiscale theme initiated in [16], where the regime of large deformation gradients in the soft component of the composite was considered. Let us emphasise two points that contrast our contribution to some earlier work within the related field. First, we note that, apart from [9,16], a number of other articles (e.g. [5,7,12]) have treated high-contrast periodic composites in the nonlinear context. However, the related results are of limited relevance to nonlinear elasticity, due to the convexity or monotonicity assumptions made in these works. In the present paper we study a class of functionals subject to the requirement of material fame indifference (see assumption (W1) in Section 2), which makes our analysis fit the fully nonlinear elasticity framework, as opposed to the works mentioned. Second, as was discussed above, the analysis of composites with “soft” inclusions within a “stiff” matrix cannot be reduced to a “decoupled” model where the perforated medium obtained by removing the inclusions is considered first and the displacement within the inclusions is found independently, which from the physics perspective can be viewed as a kind of resonance phenomenon; cf. (9) in the linearisation regime, for which an inherent energy coupling, in the limit as $\mathit{\epsilon }\to 0$ , between the soft and stiff components of the composite is essential. On a related note, the proof of the key compactness statement (Lemma 1) involves the simultaneous analysis of the displacements on the two components. We would also like to highlight the fact that in [16] the order of the relative scaling of the displacements on the soft and stiff components of the composite are assumed from the outset, while in the present work it is the result of the above compactness argument itself.

Organisation of the paper. In Section 2 we state the assumptions on the geometry of the composite and the material law. In Section 3 we present the main results, starting with results regarding two-scale compactness, convergence results in the small strain regime and finally the convergence result in the finite strain regime. All proofs are presented in Section 4.

1.1. Notation

Here we list some notation that we use throughout the text. Additional items will be introduced whenever they are first used in the text.

$\mathit{d}⩾2$ is the (integer) dimension of the space occupied by the material.

$\mathit{p}⩾1$ is the exponent in the notation ${\mathit{L}}^{\mathit{p}}$ for a Lebesgue space.

$\mathit{Y}:={\mathrm{(}0\mathrm{,}1\mathrm{)}}^{\mathit{d}}$ the reference period cell; ${\mathit{Y}}^0$ is an open Lipschitz set whose closure is contained in Y, and ${\mathit{Y}}^1:=\mathit{Y}\setminus \overline{{\mathit{Y}}^0}$ .

Ω, ${\mathrm{\Omega }}_{\mathit{\epsilon }}^0$ and ${\mathrm{\Omega }}_{\mathit{\epsilon }}^1$ denote the reference domains of the composite, the set occupied by the pore material, and the domain occupied by the matrix material, respectively, see Section 2 for precise definition.

Unless stated otherwise, all function spaces ${\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}$ , ${\mathit{H}}^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ , ${\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ , etc. consist of functions taking values in ${\mathbb{R}}^{\mathit{d}}$ .

Function spaces whose notation contains subscript “c” consist of functions that vanish outside a compact set.

The function spaces ${\mathit{H}}_{\mathrm{\#}}^1$ , ${\mathit{H}}_0^1\mathrm{(}{\mathit{Y}}^0\mathrm{)}$ , and $\mathcal{A}\mathrm{(}{\mathit{Y}}^0\mathrm{)}$ are introduced in Section 3.1.

We write · and : for the canonical inner products in ${\mathbb{R}}^{\mathit{d}}$ and ${\mathbb{R}}^{\mathit{d}\times \mathit{d}}$ , respectively.

$\mathrm{SO}\mathrm{(}\mathit{d}\mathrm{)}$ denotes the set of rotations in ${\mathbb{R}}^{\mathit{d}\times \mathit{d}}$ .

≲ stands for ⩽ up to a multiplicative constant that only depends on d, ${\mathit{Y}}^1$ , Ω, and on p if applicable.

2. Geometric and constitutive setup

The pore geometry. The set ${\mathit{Y}}^0$ defined above describes the “pores” contained within the cell Y. Note that ${\mathit{Y}}^1$ is an open, bounded, connected set with Lipschitz boundary. Therefore, to each $\mathit{\phi }\in {\mathit{H}}^1\mathrm{(}{\mathit{Y}}^1\mathrm{)}$ we can associate (see e.g. [33]) a unique harmonic extension ${\mathit{g}}^1\in {\mathit{H}}^1\mathrm{(}\mathit{Y}\mathrm{)}$ characterised by $$\begin{array}{}\text{(10)}& {\mathit{g}}^1=\mathit{\phi }\text{in }{\mathit{Y}}^1\mathrm{,}{{\displaystyle \int }}_{{\mathit{Y}}^0}\nabla {\mathit{g}}^1:\nabla \mathit{\zeta }\mathit{d}\mathit{y}=0\forall \mathit{\zeta }\in {\mathit{H}}_0^1\mathrm{\left(}{\mathit{Y}}^0\mathrm{\right)}.\end{array}$$ For this extension the inequality $$\begin{array}{}\text{(11)}& {\Vert \nabla {\mathit{g}}^1\Vert }_{{\mathit{L}}^2\mathrm{(}{\mathit{Y}}^0\mathrm{)}}⩽\mathit{C}\Vert \nabla \mathit{\phi }{\Vert }_{{\mathit{L}}^2\mathrm{(}{\mathit{Y}}^1\mathrm{)}}\end{array}$$ holds with a constant C that only depends on ${\mathit{Y}}^1$ .

For a given domain $\mathrm{\Omega }\subset {\mathbb{R}}^{\mathit{d}}$ and $\mathit{\epsilon }\mathrm{>}0$ , we define the sets ${\mathrm{\Omega }}_{\mathit{\epsilon }}^0$ and ${\mathrm{\Omega }}_{\mathit{\epsilon }}^1$ as follows: $$\begin{array}{c}{\mathrm{\Omega }}_{\mathit{\epsilon }}^0:={\displaystyle \bigcup }\left\{\mathit{\epsilon }\mathrm{\right(}\mathit{\xi }+{\mathit{Y}}^0\mathrm{)}|\mathit{\xi }\in {\mathbb{Z}}^{\mathit{d}}\mathrm{,}\mathit{\epsilon }\mathrm{(}\mathit{\xi }+\mathit{Y}\mathrm{)}\subset \mathrm{\Omega }\}\mathrm{,}{\mathrm{\Omega }}_{\mathit{\epsilon }}^1:=\mathrm{\Omega }\setminus \overline{{\mathrm{\Omega }}_{\mathit{\epsilon }}^0}.\end{array}$$ Note that by construction ${\mathrm{\Omega }}_{\mathit{\epsilon }}^1$ is a Lipschitz domain. In particular, it is connected and $\partial \mathrm{\Omega }\subset \partial {\mathrm{\Omega }}_{\mathit{\epsilon }}^1$ . We denote by ${\mathit{\chi }}_{\mathit{\epsilon }}$ the indicator function of the set of pores: $$\begin{array}{c}{\mathit{\chi }}_{\mathit{\epsilon }}\mathrm{(}\mathit{x}\mathrm{)}:=\left\{\begin{array}{cc}1\mathrm{,}\hfill & \mathit{x}\in {\mathrm{\Omega }}_{\mathit{\epsilon }}^0\mathrm{,}\hfill \\ 0\mathrm{,}\hfill & \mathit{x}\in {\mathbb{R}}^{\mathit{d}}\setminus {\mathrm{\Omega }}_{\mathit{\epsilon }}^0.\hfill \end{array}\right.\end{array}$$

The composite. The two materials are described by energy densities ${\mathit{W}}^{\mathit{i}}:{\mathbb{R}}^{\mathit{d}\times \mathit{d}}\to [0\mathrm{,}+\infty ]$ , $\mathit{i}=0\mathrm{,}1$ . Unless stated otherwise, we assume that for $\mathit{i}=0\mathrm{,}1$ :

${\mathit{W}}^{\mathit{i}}$ is frame-indifferent, i.e. ${\mathit{W}}^{\mathit{i}}\mathrm{(}\mathit{R}\mathit{F}\mathrm{)}={\mathit{W}}^{\mathit{i}}\mathrm{(}\mathit{F}\mathrm{)}$ for all $\mathit{R}\in \mathrm{SO}\mathrm{(}\mathit{d}\mathrm{)}$ and all $\mathit{F}\in {\mathbb{R}}^{\mathit{d}\times \mathit{d}}$ ;

The identity matrix $\mathit{I}\in {\mathbb{R}}^{\mathit{d}\times \mathit{d}}$ is a “natural state”, i.e. ${\mathit{W}}^{\mathit{i}}\mathrm{(}\mathit{I}\mathrm{)}=0$ , and ${\mathit{W}}^{\mathit{i}}$ is non-degenerate: $$\begin{array}{c}{\mathit{W}}^{\mathit{i}}\mathrm{(}\mathit{F}\mathrm{)}⩾{\mathit{c}}_0{dist}^2\mathrm{(}\mathit{F}\mathrm{,}\mathrm{SO}\mathrm{(}\mathit{d}\mathrm{)}\mathrm{)}\mathrm{,}\forall \mathit{F}\in {\mathbb{R}}^{\mathit{d}\times \mathit{d}}\mathrm{,}{\mathit{c}}_0\mathrm{>}0.\end{array}$$

${\mathit{W}}^{\mathit{i}}$ has a quadratic expansion at I, i.e. there exists a non-negative quadratic form ${\mathit{Q}}^{\mathit{i}}$ on ${\mathbb{R}}^{\mathit{d}\times \mathit{d}}$ and an increasing function ${\mathit{r}}^{\mathit{i}}:[0\mathrm{,}\infty \mathrm{)}\to [0\mathrm{,}\infty ]$ with ${lim}_{\mathit{s}\downarrow 0}{\mathit{r}}^{\mathit{i}}\mathrm{(}\mathit{s}\mathrm{)}=0$ , such that $$\begin{array}{}\text{(12)}& |{\mathit{W}}^{\mathit{i}}\mathrm{(}\mathit{I}+\mathit{G}\mathrm{)}-{\mathit{Q}}^{\mathit{i}}\mathrm{(}\mathit{G}\mathrm{)}|⩽|\mathit{G}{|}^2{\mathit{r}}^{\mathit{i}}\mathrm{(}|\mathit{G}|\mathrm{)}\forall \mathit{G}\in {\mathbb{R}}^{\mathit{d}\times \mathit{d}}.\end{array}$$

As shown in [30, Lemma 2.7] the quadratic form ${\mathit{Q}}^{\mathit{i}}$ associated with ${\mathit{W}}^{\mathit{i}}$ via (W3) satisfies $$\begin{array}{}\text{(13)}& {\mathit{c}}_1|sym\mathit{G}{|}^2⩽{\mathit{Q}}^{\mathit{i}}\mathrm{(}\mathit{G}\mathrm{)}={\mathit{Q}}^{\mathit{i}}\mathrm{(}sym\mathit{G}\mathrm{)}⩽{\mathit{c}}_1^{-1}|sym\mathit{G}{|}^2\forall \mathit{G}\in {\mathbb{R}}^{\mathit{d}\times \mathit{d}}\mathrm{,}{\mathit{c}}_1\mathrm{>}0.\end{array}$$ In the finite strain regime we consider a different set of assumptions for ${\mathit{W}}^0$ , which are listed in Section 3.3.

The scaling parameter γ. Throughout the paper $\mathit{\gamma }\mathrm{>}0$ denotes a fixed scaling parameter. It is a quantitative measure of the relative contrast between the two components of the composite.

Energy functional. We define the elastic energy as a functional of the displacement, as follows: $$\begin{array}{}\text{(14)}& {\mathcal{E}}_{\mathit{\epsilon }}\mathrm{(}\mathit{u}\mathrm{)}={{\displaystyle \int }}_{\mathrm{\Omega }}\mathrm{(}{\mathit{\epsilon }}^{2\mathit{\gamma }}{\mathit{W}}^0\mathrm{(}\mathit{I}+\nabla \mathit{u}\mathrm{)}{\mathit{\chi }}_{\mathit{\epsilon }}+{\mathit{W}}^1\mathrm{(}\mathit{I}+\nabla \mathit{u}\mathrm{)}\mathrm{(}1-{\mathit{\chi }}_{\mathit{\epsilon }}\mathrm{)}\mathrm{)}\mathit{d}\mathit{x}\mathrm{,}\mathit{u}\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{,}{\mathbb{R}}^{\mathit{d}}\mathrm{)}.\end{array}$$

3. Main results

3.1. Compactness and two-scale convergence

We first present an a priori estimate and a two-scale compactness statement for sequences ${\mathit{\phi }}_{\mathit{\epsilon }}\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ whose energy is equi-bounded in the sense that $$\begin{array}{}\text{(15)}& \underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}sup}{\mathrm{\Phi }}_{\mathit{\epsilon }}^{\mathit{\gamma }}\mathrm{(}{\mathit{v}}_{\mathit{\epsilon }}\mathrm{)}\mathrm{<}\infty \mathrm{,}\end{array}$$ where $$\begin{array}{c}{\mathrm{\Phi }}_{\mathit{\epsilon }}^{\mathit{\gamma }}\mathrm{(}\mathit{v}\mathrm{)}:={{\displaystyle \int }}_{\mathrm{\Omega }}{dist}^2\mathrm{(}\mathit{I}+\nabla \mathit{v}\mathrm{(}\mathit{x}\mathrm{)}\mathrm{,}\mathrm{SO}\mathrm{(}\mathit{d}\mathrm{)}\mathrm{)}\mathrm{(}{\mathit{\epsilon }}^{2\mathit{\gamma }}{\mathit{\chi }}_{\mathit{\epsilon }}+\mathrm{(}1-{\mathit{\chi }}_{\mathit{\epsilon }}\mathrm{)}\mathrm{)}\mathit{d}\mathit{x}.\end{array}$$ Note that, by virtue of the non-degeneracy assumption (W2) the functional ${\mathrm{\Phi }}_{\mathit{\epsilon }}^{\mathit{\gamma }}\mathrm{(}·\mathrm{)}$ bounds below ${\mathcal{E}}_{\mathit{\epsilon }}\mathrm{(}\mathrm{id}+·\mathrm{)}$ , where $\mathrm{id}\mathrm{(}\mathit{x}\mathrm{)}=\mathit{x}$ , $\mathit{x}\in \mathrm{\Omega }$ . As we shall see in the upcoming Lemma 1, the inequality (15) implies that the sequence ${\mathit{\phi }}_{\mathit{\epsilon }}$ is bounded in ${\mathit{H}}^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ , and thus weakly converges (up to extracting a subsequence) to a limit displacement $\mathit{\phi }\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ . For our purpose we require a precise understanding of the oscillations that emerge along that limit. We achieve this by combining two concepts:

• We write a representation for ${\mathit{\phi }}_{\mathit{\epsilon }}$ in the spirit of an asymptotic decomposition as $\mathit{\epsilon }\downarrow 0$ .

• We study the convergence properties of the terms in this decomposition by appealing to two-scale convergence.

In the following lemma we address the first item above.

Lemma 1.

Let ${\mathit{\phi }}_{\mathit{\epsilon }}\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ and $0\mathrm{<}\mathit{\epsilon }⩽1$ .

There exists a unique pair of functions ${\mathit{g}}_{\mathit{\epsilon }}^0\in {\mathit{H}}_0^1\mathrm{(}{\mathrm{\Omega }}_{\mathit{\epsilon }}^0\mathrm{)}$ and ${\mathit{g}}_{\mathit{\epsilon }}^1\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ such that $$\begin{array}{}\text{(16)}& \begin{array}{cc}\hfill & \mathrm{(}\mathrm{i}\mathrm{)}{\mathit{\phi }}_{\mathit{\epsilon }}={\mathit{g}}_{\mathit{\epsilon }}^1+{\mathit{\epsilon }}^{1-\mathit{\gamma }}{\mathit{g}}_{\mathit{\epsilon }}^0\mathrm{,}\hfill \\ \hfill & \mathrm{(}\mathrm{ii}\mathrm{)}{{\displaystyle \int }}_{{\mathrm{\Omega }}_{\mathit{\epsilon }}^0}\nabla {\mathit{g}}_{\mathit{\epsilon }}^1:\nabla \mathit{\zeta }=0\forall \mathit{\zeta }\in {\mathit{H}}_0^1\mathrm{\left(}{\mathrm{\Omega }}_{\mathit{\epsilon }}^0\mathrm{\right)}.\hfill \end{array}\end{array}$$

There exists a positive constant C that only depends on $\mathrm{\Omega }\mathrm{,}{\mathit{Y}}^0$ such that $$\begin{array}{}\text{(17)}& {\Vert {\mathit{g}}_{\mathit{\epsilon }}^0\Vert }_{{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}}^2+{\Vert \mathit{\epsilon }\nabla {\mathit{g}}_{\mathit{\epsilon }}^0\Vert }_{{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}}^2+{\Vert {\mathit{g}}_{\mathit{\epsilon }}^1\Vert }_{{\mathit{H}}^1\mathrm{(}\mathrm{\Omega }\mathrm{)}}^2⩽\mathit{C}{\mathrm{\Phi }}_{\mathit{\epsilon }}^{\mathit{\gamma }}\mathrm{(}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}\mathrm{,}\end{array}$$ where ${\mathit{g}}_{\mathit{\epsilon }}^0$ , ${\mathit{g}}_{\mathit{\epsilon }}^1$ and ${\mathit{\phi }}_{\mathit{\epsilon }}$ are related to each other as in (a).

As already explained in the introduction, for our purpose it is convenient to appeal to two-scale convergence, see Definition 2. We use the following shorthand notation: $$\begin{array}{rcl}{\mathit{f}}_{\mathit{\epsilon }}\to {\mathit{f}}_0& :\iff & {\mathit{f}}_{\mathit{\epsilon }}\text{ strongly converges to }{\mathit{f}}_0\text{ in }{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}\mathrm{,}\\ {\mathit{f}}_{\mathit{\epsilon }}\mathrm{\rightharpoonup }{\mathit{f}}_0& :\iff & {\mathit{f}}_{\mathit{\epsilon }}\text{ weakly converges to }{\mathit{f}}_0\text{ in }{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}\mathrm{,}\\ {\mathit{f}}_{\mathit{\epsilon }}\stackrel{2}{\mathrm{\rightharpoonup }}\mathit{f}& :\iff & {\mathit{f}}_{\mathit{\epsilon }}\text{ weakly two-scale converges to }\mathit{f}\text{ in }{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\times \mathit{Y}\mathrm{)}\mathrm{,}\\ {\mathit{f}}_{\mathit{\epsilon }}\stackrel{2}{\to }\mathit{f}& :\iff & {\mathit{f}}_{\mathit{\epsilon }}\text{ strongly two-scale converges to }\mathit{f}\text{ in }{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\times \mathit{Y}\mathrm{)}.\end{array}$$ The upcoming lemma states a two-scale compactness result for the displacements ${\mathit{g}}_{\mathit{\epsilon }}^0$ and ${\mathit{g}}_{\mathit{\epsilon }}^1$ that appear in the representation (16). Due to the differential constraint satisfied by ${\mathit{g}}_{\mathit{\epsilon }}^0$ , the corresponding two-scale limits automatically satisfy certain structural properties, which can be captured with the help of the following function spaces:

${\mathit{H}}_{\mathrm{\#}}^1$ is the space of ${[0\mathrm{,}1\mathrm{)}}^{\mathit{d}}$ -periodic functions in ${\mathit{H}}_{\mathrm{loc}}^1\mathrm{(}{\mathbb{R}}^{\mathit{d}}\mathrm{)}$ .

${\mathit{H}}_0^1\mathrm{(}{\mathit{Y}}^0\mathrm{)}$ is the closed subspace of ${\mathit{H}}_{\mathrm{\#}}^1$ consisting of functions $\mathit{\psi }\in {\mathit{H}}_{\mathrm{\#}}^1$ with $\mathit{\psi }=0$ on ${\mathit{Y}}^1$ .

$\mathcal{A}\mathrm{(}{\mathit{Y}}^0\mathrm{)}$ is the closed subspace of ${\mathit{H}}_{\mathrm{\#}}^1$ consisting of functions $\mathit{\psi }\in {\mathit{H}}_{\mathrm{\#}}^1$ that satisfy the identity $$\begin{array}{c}{{\displaystyle \int }}_{{\mathit{Y}}^0}{\nabla }_{\mathit{y}}\mathit{\psi }:{\nabla }_{\mathit{y}}\mathit{\zeta }\mathit{d}\mathit{y}=0\forall \mathit{\zeta }\in {\mathit{H}}_0^1\mathrm{\left(}{\mathit{Y}}^0\mathrm{\right)}.\end{array}$$

Lemma 2.

Consider a sequence ${\mathit{\phi }}_{\mathit{\epsilon }}\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ and let $\mathrm{(}{\mathit{g}}_{\mathit{\epsilon }}^0\mathrm{,}{\mathit{g}}_{\mathit{\epsilon }}^1\mathrm{)}$ be associated with ${\mathit{\phi }}_{\mathit{\epsilon }}$ via ( 16 ). Suppose that there exists a sequence of positive numbers ${\mathit{m}}_{\mathit{\epsilon }}$ such that $$\begin{array}{c}\underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}sup}{\mathit{m}}_{\mathit{\epsilon }}^{-2}{\mathrm{\Phi }}_{\mathit{\epsilon }}^{\mathit{\gamma }}\mathrm{(}{\mathit{m}}_{\mathit{\epsilon }}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}\mathrm{<}\infty \mathit{\text{and}}{\mathit{m}}_{\mathit{\epsilon }}=\mathit{O}\mathrm{\left(}{\mathit{\epsilon }}^{\mathit{\gamma }}\mathrm{\right)}.\end{array}$$ Then there exist $$\begin{array}{}\text{(18)}& {\mathit{g}}^0\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{,}{\mathit{H}}_0^1\mathrm{(}{\mathit{Y}}^0\mathrm{\left)}\mathrm{\right)}\mathrm{,}{\mathit{g}}^1\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}\mathrm{,}\mathit{\psi }\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{,}\mathcal{A}\mathrm{(}{\mathit{Y}}^0\mathrm{\left)}\mathrm{\right)}\end{array}$$ such that, up to selecting a subsequence, one has $$\begin{array}{}\text{(19)}& \begin{array}{cc}\hfill & {\mathit{g}}_{\mathit{\epsilon }}^0\stackrel{2}{\mathrm{\rightharpoonup }}{\mathit{g}}^0\mathrm{,}\mathit{\epsilon }\nabla {\mathit{g}}_{\mathit{\epsilon }}^0\stackrel{2}{\mathrm{\rightharpoonup }}{\nabla }_{\mathit{y}}{\mathit{g}}^0\mathrm{,}\hfill \\ \hfill & {\mathit{g}}_{\mathit{\epsilon }}^1\mathrm{\rightharpoonup }{\mathit{g}}^1\mathit{\text{weakly in }}{\mathit{H}}^1\mathrm{(}\mathrm{\Omega }\mathrm{)}\mathit{\text{and}}\nabla {\mathit{g}}_{\mathit{\epsilon }}^1\stackrel{2}{\mathrm{\rightharpoonup }}\nabla {\mathit{g}}^1+{\nabla }_{\mathit{y}}\mathit{\psi }.\hfill \end{array}\end{array}$$

The identification obtained in the previous lemma is sharp, in the sense of the following statement.

Lemma 3.

Let ${\mathit{g}}^0\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{,}{\mathit{H}}_0^1\mathrm{(}{\mathit{Y}}^0\mathrm{)}\mathrm{)}$ , ${\mathit{g}}^1\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ and $\mathit{\psi }\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{,}\mathcal{A}\mathrm{(}{\mathit{Y}}^0\mathrm{)}\mathrm{)}$ . Let ${\mathit{c}}_{\mathit{\epsilon }}$ be an arbitrary sequence of positive numbers converging to zero. Then there exist function sequences ${\mathit{g}}_{\mathit{\epsilon }}^0\in {\mathit{H}}_0^1\mathrm{(}{\mathrm{\Omega }}_{\mathit{\epsilon }}^0\mathrm{)}$ , ${\mathit{g}}_{\mathit{\epsilon }}^1\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ such that $\mathrm{(}{\mathit{g}}_{\mathit{\epsilon }}^0\mathrm{,}{\mathit{g}}_{\mathit{\epsilon }}^1\mathrm{)}$ is related to ${\mathit{\phi }}_{\mathit{\epsilon }}:={\mathit{g}}_{\mathit{\epsilon }}^1+{\mathit{\epsilon }}^{1-\mathit{\gamma }}{\mathit{g}}_{\mathit{\epsilon }}^0$ as in ( 16 ), and $$\begin{array}{}\text{(20)}& \begin{array}{cc}\hfill & {\mathit{g}}_{\mathit{\epsilon }}^0\stackrel{2}{\to }{\mathit{g}}^0\mathrm{,}\mathit{\epsilon }\nabla {\mathit{g}}_{\mathit{\epsilon }}^0\stackrel{2}{\to }{\nabla }_{\mathit{y}}{\mathit{g}}^0\mathrm{,}\hfill \\ \hfill & {\mathit{g}}_{\mathit{\epsilon }}^1\mathrm{\rightharpoonup }{\mathit{g}}^1\mathit{\text{weakly in }}{\mathit{H}}^1\mathrm{(}\mathrm{\Omega }\mathrm{)}\mathit{\text{and}}\nabla {\mathit{g}}_{\mathit{\epsilon }}^1\stackrel{2}{\to }\nabla {\mathit{g}}^1+{\nabla }_{\mathit{y}}\mathit{\psi }\mathrm{,}\hfill \\ \hfill & \underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}sup}{\mathit{c}}_{\mathit{\epsilon }}\Vert \nabla {\mathit{\phi }}_{\mathit{\epsilon }}{\Vert }_{{\mathit{L}}^{\infty }\mathrm{(}\mathrm{\Omega }\mathrm{)}}=0.\hfill \end{array}\end{array}$$

Our main result is formulated in terms of the notion of convergence described in the above lemmas. For convenience we use the following notation:

Given ${\mathit{\phi }}_{\mathit{\epsilon }}\in {\mathit{H}}^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ we write ${\mathit{g}}_{\mathit{\epsilon }}^1+{\mathit{\epsilon }}^{1-\mathit{\gamma }}{\mathit{g}}_{\mathit{\epsilon }}^0\stackrel{\text{(16)}}{:=}{\mathit{\phi }}_{\mathit{\epsilon }}$ , if ${\mathit{g}}_{\mathit{\epsilon }}^1\in {\mathit{H}}^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ , ${\mathit{g}}_{\mathit{\epsilon }}^0\in {\mathit{H}}_0^1\mathrm{(}{\mathrm{\Omega }}_{\mathit{\epsilon }}^0\mathrm{)}$ , and both functions are related to ${\mathit{\phi }}_{\mathit{\epsilon }}$ as in (16).

We write ${\mathit{\phi }}_{\mathit{\epsilon }}\stackrel{2}{\mathrm{\rightharpoonup }}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}$ , if ${\mathit{g}}_{\mathit{\epsilon }}^1+{\mathit{\epsilon }}^{1-\mathit{\gamma }}{\mathit{g}}_{\mathit{\epsilon }}^0\stackrel{\text{(16)}}{:=}{\mathit{\phi }}_{\mathit{\epsilon }}$ and $$\begin{array}{}\text{(21)}& {\mathit{g}}_{\mathit{\epsilon }}^0\stackrel{2}{\mathrm{\rightharpoonup }}{\mathit{g}}^0\mathrm{,}\mathit{\epsilon }\nabla {\mathit{g}}_{\mathit{\epsilon }}^0\stackrel{2}{\mathrm{\rightharpoonup }}{\nabla }_{\mathit{y}}{\mathit{g}}^0\mathrm{,}{\mathit{g}}_{\mathit{\epsilon }}^1\mathrm{\rightharpoonup }{\mathit{g}}^1\text{weakly in }{\mathit{H}}^1\mathrm{(}\mathrm{\Omega }\mathrm{)}.\end{array}$$

We write ${\mathit{\phi }}_{\mathit{\epsilon }}\stackrel{2}{\to }\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}$ , if ${\mathit{g}}_{\mathit{\epsilon }}^1+{\mathit{\epsilon }}^{1-\mathit{\gamma }}{\mathit{g}}_{\mathit{\epsilon }}^0\stackrel{\text{(16)}}{:=}{\mathit{\phi }}_{\mathit{\epsilon }}$ and $$\begin{array}{}\text{(22)}& {\mathit{g}}_{\mathit{\epsilon }}^0\stackrel{2}{\to }{\mathit{g}}^0\mathrm{,}\mathit{\epsilon }\nabla {\mathit{g}}_{\mathit{\epsilon }}^0\stackrel{2}{\to }{\nabla }_{\mathit{y}}{\mathit{g}}^0\mathrm{,}{\mathit{g}}_{\mathit{\epsilon }}^1\mathrm{\rightharpoonup }{\mathit{g}}^1\text{weakly in }{\mathit{H}}^1\mathrm{(}\mathrm{\Omega }\mathrm{)}.\end{array}$$

3.2. Convergence in the small strain regime ${\mathit{m}}_{\mathit{\epsilon }}=\mathit{o}\mathrm{(}{\mathit{\epsilon }}^{\mathit{\gamma }}\mathrm{)}$

Throughout this section we assume that the densities ${\mathit{W}}^0$ and ${\mathit{W}}^1$ satisfy the conditions (W1)–(W3). We show that in the small strain regime the limit functional $$\begin{array}{c}{\mathcal{E}}_{\text{small}}:{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{,}{\mathit{H}}_0^1\mathrm{(}{\mathit{Y}}^0\mathrm{\left)}\mathrm{\right)}\times {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}\to [0\mathrm{,}\infty \mathrm{)}\mathrm{,}\end{array}$$ is given by $$\begin{array}{c}{\mathcal{E}}_{\text{small}}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}:={{\displaystyle \int }}_{\mathrm{\Omega }\times \mathit{Y}}\mathrm{\left(}{\mathit{Q}}^0\mathrm{\right(}{\nabla }_{\mathit{y}}{\mathit{g}}^0\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}\mathrm{)}+{\mathit{Q}}_{\mathrm{hom}}^1\mathrm{(}\nabla {\mathit{g}}^1\mathrm{(}\mathit{x}\mathrm{)}\mathrm{\left)}\mathrm{\right)}\mathit{d}\mathit{x}\mathrm{,}\end{array}$$ where $$\begin{array}{}\text{(23)}& {\mathit{Q}}_{\mathrm{hom}}^1\mathrm{(}\mathit{F}\mathrm{)}:=\underset{\mathit{\psi }\in \mathcal{A}\mathrm{(}{\mathit{Y}}^0\mathrm{)}}{min}{{\displaystyle \int }}_{{\mathit{Y}}^1}{\mathit{Q}}^1\mathrm{(}\mathit{F}+{\nabla }_{\mathit{y}}\mathit{\psi }\mathrm{(}\mathit{y}\mathrm{)}\mathrm{)}\mathit{d}\mathit{y}.\end{array}$$ More precisely, the following theorem holds.

Theorem 1.

Let ${\mathit{m}}_{\mathit{\epsilon }}$ be a sequence of positive numbers and assume that ${\mathit{m}}_{\mathit{\epsilon }}=\mathit{o}\mathrm{(}{\mathit{\epsilon }}^{\mathit{\gamma }}\mathrm{)}$ as $\mathit{\epsilon }\downarrow 0$ .

(Compactness). Suppose that ${\mathit{\phi }}_{\mathit{\epsilon }}\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ satisfy $$\begin{array}{c}\underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}sup}{\mathit{m}}_{\mathit{\epsilon }}^{-2}{\mathcal{E}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{m}}_{\mathit{\epsilon }}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}\mathrm{<}\infty .\end{array}$$ Then, up to a subsequence, one has ${\mathit{\phi }}_{\mathit{\epsilon }}\stackrel{2}{\mathrm{\rightharpoonup }}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}$ for some ${\mathit{g}}^0\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{,}{\mathit{H}}_0^1\mathrm{(}{\mathit{Y}}^0\mathrm{)}\mathrm{)}$ and ${\mathit{g}}^1\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ .

(Lower bound). Consider ${\mathit{\phi }}_{\mathit{\epsilon }}\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ and suppose that ${\mathit{\phi }}_{\mathit{\epsilon }}\stackrel{2}{\mathrm{\rightharpoonup }}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}$ for some ${\mathit{g}}^0\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{,}{\mathit{H}}_0^1\mathrm{(}{\mathit{Y}}^0\mathrm{)}\mathrm{)}$ and ${\mathit{g}}^1\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ . Then the estimate $$\begin{array}{c}\underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}inf}{\mathit{m}}_{\mathit{\epsilon }}^{-2}{\mathcal{E}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{m}}_{\mathit{\epsilon }}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}⩾{\mathcal{E}}_{\mathrm{small}}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}\end{array}$$ holds.

(Recovery sequence). For all ${\mathit{g}}^0\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{,}{\mathit{H}}_0^1\mathrm{(}{\mathit{Y}}^0\mathrm{)}\mathrm{)}$ and ${\mathit{g}}^1\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ there exists a sequence ${\mathit{\phi }}_{\mathit{\epsilon }}\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ such that ${\mathit{\phi }}_{\mathit{\epsilon }}\stackrel{2}{\to }\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}$ and $$\begin{array}{c}\underset{\mathit{\epsilon }\downarrow 0}{lim}{\mathit{m}}_{\mathit{\epsilon }}^{-2}{\mathcal{E}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{m}}_{\mathit{\epsilon }}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}={\mathcal{E}}_{\mathrm{small}}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}.\end{array}$$

In the next result we consider a minimisation problem that involves the density of the “body forces” ${\ell }_{\mathit{\epsilon }}\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}$ . We study the variational limit of the (scaled) total energy $$\begin{array}{}\text{(24)}& {\mathit{I}}_{\mathit{\epsilon }}\mathrm{(}\mathit{\phi }\mathrm{)}:={\displaystyle \frac{1}{{\mathit{m}}_{\mathit{\epsilon }}^2}}\mathrm{(}{\mathcal{E}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{m}}_{\mathit{\epsilon }}\mathit{\phi }\mathrm{)}-{{\displaystyle \int }}_{\mathrm{\Omega }}{\ell }_{\mathit{\epsilon }}·\mathrm{(}{\mathit{m}}_{\mathit{\epsilon }}\mathit{\phi }\mathrm{)}\mathit{d}\mathit{x}\mathrm{)}\mathrm{,}\end{array}$$ where the scaling factor ${\mathit{m}}_{\mathit{\epsilon }}$ is determined by the body forces via $$\begin{array}{}\text{(25a)}& {\mathit{m}}_{\mathit{\epsilon }}:={\mathit{\epsilon }}^{1-\mathit{\gamma }}\Vert {\ell }_{\mathit{\epsilon }}{\Vert }_{{\mathit{L}}^2\mathrm{(}{\mathrm{\Omega }}_{\mathit{\epsilon }}^0\mathrm{)}}+\Vert {\ell }_{\mathit{\epsilon }}{\Vert }_{{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}}.\end{array}$$ In the small strain regime we assume that the body forces are small in the sense that $$\begin{array}{}\text{(25b)}& {\mathit{m}}_{\mathit{\epsilon }}=\mathit{o}\mathrm{\left(}{\mathit{\epsilon }}^{\mathit{\gamma }}\mathrm{\right)}\mathrm{,}\mathit{\epsilon }\downarrow 0.\end{array}$$ Moreover, we assume that the (scaled) body-force densities converge, as $\mathit{\epsilon }\downarrow 0$ , in the following way: $$\begin{array}{}\text{(25c)}& {\mathit{m}}_{\mathit{\epsilon }}^{-1}{\mathit{\epsilon }}^{1-\mathit{\gamma }}{\mathit{\chi }}_{\mathit{\epsilon }}{\ell }_{\mathit{\epsilon }}\stackrel{2}{\to }{\ell }^0\mathrm{,}{\mathit{m}}_{\mathit{\epsilon }}^{-1}{\ell }_{\mathit{\epsilon }}\mathrm{\rightharpoonup }{\ell }^1.\end{array}$$ It follows from Theorem 1 that the variational limit of the total energy (24) is given by the functional $$\begin{array}{}\text{(26)}& {\mathit{I}}_{\text{small}}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}:={\mathcal{E}}_{\text{small}}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}-{{\displaystyle \int }}_{\mathrm{\Omega }}\mathrm{(}{{\displaystyle \int }}_{{\mathit{Y}}^0}{\ell }^0·{\mathit{g}}^0\mathit{d}\mathit{y}+{\ell }^1·{\mathit{g}}^1\mathrm{)}\mathit{d}\mathit{x}.\end{array}$$

Proposition 1.

Assume that ( 25a )–( 25c ) hold.

(Convergence of infima). One has $$\begin{array}{c}\underset{\mathit{\epsilon }\downarrow 0}{lim}\underset{\mathit{\phi }\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}}{inf}{\mathit{I}}_{\mathit{\epsilon }}\mathrm{(}\mathit{\phi }\mathrm{)}=min{\mathit{I}}_{\mathrm{small}}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}\mathrm{,}\end{array}$$ where the minimum on the right-hand side is taken over all ${\mathit{g}}^0\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{,}{\mathit{H}}_0^1\mathrm{(}{\mathit{Y}}^0\mathrm{)}\mathrm{)}$ and ${\mathit{g}}^1\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ . Moreover, the minimum is attained for a unique pair $\mathrm{(}{\mathit{g}}_{\ast }^0\mathrm{,}{\mathit{g}}_{\ast }^1\mathrm{)}$ .

(Convergence of minimisers). Let ${\mathit{\phi }}_{\mathit{\epsilon }}\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ be a sequence of almost minimisers, i.e. $$\begin{array}{}\text{(27)}& {\mathit{I}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}⩽\underset{\mathit{\phi }\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}}{inf}{\mathit{I}}_{\mathit{\epsilon }}\mathrm{(}\mathit{\phi }\mathrm{)}+\mathit{o}\mathrm{(}1\mathrm{)}\mathrm{,}\mathit{\epsilon }\downarrow 0.\end{array}$$ Then $$\begin{array}{c}{\mathit{\phi }}_{\mathit{\epsilon }}\stackrel{2}{\mathrm{\rightharpoonup }}\mathrm{(}{\mathit{g}}_{\ast }^0\mathrm{,}{\mathit{g}}_{\ast }^1\mathrm{)}\mathit{\text{and}}\nabla {\mathit{g}}_{\mathit{\epsilon }}^1\stackrel{2}{\mathrm{\rightharpoonup }}\nabla {\mathit{g}}_{\ast }^1+{\nabla }_{\mathit{y}}{\mathit{\psi }}_{\ast }\mathrm{,}\end{array}$$ where ${\mathit{\psi }}_{\ast }\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{,}\mathcal{A}\mathrm{(}{\mathit{Y}}^0\mathrm{)}\mathrm{)}$ denotes the unique “corrector” characterised by $$\begin{array}{}\text{(28)}& {\mathit{Q}}_{\mathrm{hom}}^1\mathrm{(}\nabla {\mathit{g}}_{\ast }^1\mathrm{(}\mathit{x}\mathrm{)}\mathrm{)}={{\displaystyle \int }}_{{\mathit{Y}}^1}{\mathit{Q}}^1\mathrm{(}\nabla {\mathit{g}}_{\ast }^1\mathrm{(}\mathit{x}\mathrm{)}+{\nabla }_{\mathit{y}}{\mathit{\psi }}_{\ast }\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}\mathrm{)}\mathit{d}\mathit{y}\mathrm{,}{{\displaystyle \int }}_{\mathit{Y}}{\mathit{\psi }}_{\ast }\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}\mathit{d}\mathit{y}=0\mathrm{,}\end{array}$$ for almost every $\mathit{x}\in \mathrm{\Omega }$ .

Next, we prove that almost minimisers ${\mathit{\phi }}_{\mathit{\epsilon }}$ satisfy the asymptotic relation $$\begin{array}{}\text{(29)}& {\mathit{\phi }}_{\mathit{\epsilon }}={\mathit{g}}_{\ast \mathrm{,}\mathit{\epsilon }}^1\mathrm{(}\mathit{x}\mathrm{)}+{\mathit{\epsilon }}^{1-\mathit{\gamma }}{\mathit{g}}_{\ast \mathrm{,}\mathit{\epsilon }}^0\mathrm{(}\mathit{x}\mathrm{)}+\mathit{o}\mathrm{(}1\mathrm{)}\mathrm{,}\mathit{\epsilon }\downarrow 0\mathrm{,}\text{in }{\mathit{W}}^{1\mathrm{,}\mathit{p}}\mathrm{(}\mathrm{\Omega }\mathrm{)}\mathrm{(}\mathit{p}\mathrm{<}2\mathrm{)}\mathrm{,}\end{array}$$ where ${\mathit{g}}_{\ast \mathrm{,}\mathit{\epsilon }}^1$ and ${\mathit{g}}_{\ast \mathrm{,}\mathit{\epsilon }}^0$ formally obey the “ansatz” $$\begin{array}{}\text{(30)}& {\mathit{g}}_{\ast \mathrm{,}\mathit{\epsilon }}^0\mathrm{(}\mathit{x}\mathrm{)}\stackrel{\text{formally}}{=}{\mathit{g}}_{\ast }^0\mathrm{(}\mathit{x}\mathrm{,}\mathit{x}\mathrm{/}\mathit{\epsilon }\mathrm{)}\mathrm{,}{\mathit{g}}_{\ast \mathrm{,}\mathit{\epsilon }}^1\mathrm{(}\mathit{x}\mathrm{)}\stackrel{\text{formally}}{=}{\mathit{g}}_{\ast }^1\mathrm{(}\mathit{x}\mathrm{)}+\mathit{\epsilon }{\mathit{\psi }}_{\ast }\mathrm{(}\mathit{x}\mathrm{,}\mathit{x}\mathrm{/}\mathit{\epsilon }\mathrm{)}.\end{array}$$ Here $\mathrm{(}{\mathit{g}}_{\ast }^0\mathrm{,}{\mathit{g}}_{\ast }^1\mathrm{)}$ and ${\mathit{\psi }}_{\ast }$ denote the minimising pair and corrector from Proposition 1. Since the functions on the right-hand sides in (30) are in general not smooth enough to define ${\mathit{g}}_{\ast \mathrm{,}\mathit{\epsilon }}^0$ and ${\mathit{g}}_{\ast \mathrm{,}\mathit{\epsilon }}^1$ by (30) directly, we use instead the approximation associated with $\mathrm{(}{\mathit{g}}_{\ast }^0\mathrm{,}{\mathit{g}}_{\ast }^1\mathrm{,}{\mathit{\psi }}_{\ast }\mathrm{)}$ via Lemma 3.

In addition to the properties of ${\mathit{Y}}^0$ assumed in Section 2, we require the following assumption on the regularity of ${\mathit{Y}}^0$ :

Assumption 1.

There exist an exponent $\mathit{p}\mathrm{<}2$ and a constant C such that for all $\mathit{\phi }\in {\mathit{H}}^1\mathrm{(}{\mathit{Y}}^1\mathrm{)}$ and ${\mathit{g}}^1\in {\mathit{H}}^1\mathrm{(}\mathit{Y}\mathrm{)}$ related via (10) we have $\Vert \nabla {\mathit{g}}^1{\Vert }_{{\mathit{L}}^{\mathit{p}}\mathrm{(}{\mathit{Y}}^0\mathrm{)}}⩽\mathit{C}\Vert \nabla \mathit{\phi }{\Vert }_{{\mathit{L}}^{\mathit{p}}\mathrm{(}{\mathit{Y}}^1\mathrm{)}}$ .

Note that Assumption 1 is satisfied if ${\mathit{Y}}^0$ can be written as the disjoint union of a finite number of Lipschitz domains ${\mathit{Y}}_1^0\mathrm{,}\dots \mathrm{,}{\mathit{Y}}_{\mathit{N}}^0$ with $\partial {\mathit{Y}}_{\mathit{i}}^0\cap \partial {\mathit{Y}}_{\mathit{j}}^0=\varnothing$ for $\mathit{i}\ne \mathit{j}$ .

Theorem 2.

Assume that ( 25a )–( 25c ) hold, and let Assumption  1 be satisfied. Let ${\mathit{\phi }}_{\mathit{\epsilon }}\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ be a sequence of almost minimisers, i.e. $$\begin{array}{c}{\mathit{I}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}⩽\underset{\mathit{\phi }\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}}{inf}{\mathit{I}}_{\mathit{\epsilon }}\mathrm{(}\mathit{\phi }\mathrm{)}+\mathit{o}\mathrm{(}1\mathrm{)}\mathrm{,}\mathit{\epsilon }\downarrow 0.\end{array}$$ Let $\mathrm{(}{\mathit{g}}_{\ast }^0\mathrm{,}{\mathit{g}}_{\ast }^1\mathrm{)}$ be the minimiser of ${\mathit{I}}_{\mathrm{small}}$ and let ${\mathit{\psi }}_{\ast }$ be defined through ( 28 ). Let $\mathrm{(}{\mathit{g}}_{\ast \mathrm{,}\mathit{\epsilon }}^0\mathrm{,}{\mathit{g}}_{\ast \mathrm{,}\mathit{\epsilon }}^1\mathrm{)}$ and ${\mathit{\phi }}_{\ast \mathrm{,}\mathit{\epsilon }}={\mathit{g}}_{\ast \mathrm{,}\mathit{\epsilon }}^1+{\mathit{\epsilon }}^{1-\mathit{\gamma }}{\mathit{g}}_{\ast \mathrm{,}\mathit{\epsilon }}^0$ be associated with $\mathrm{(}{\mathit{g}}_{\ast }^0\mathrm{,}{\mathit{g}}_{\ast }^1\mathrm{,}{\mathit{\psi }}_{\ast }\mathrm{)}$ as in Lemma  3 , i.e. ${\mathit{\phi }}_{\ast \mathrm{,}\mathit{\epsilon }}\stackrel{2}{\to }\mathrm{(}{\mathit{g}}_{\ast }^0\mathrm{,}{\mathit{g}}_{\ast }^1\mathrm{)}$ and $\nabla {\mathit{g}}_{\ast \mathrm{,}\mathit{\epsilon }}^1\stackrel{2}{\to }\nabla {\mathit{g}}_{\ast }^1+{\nabla }_{\mathit{y}}{\mathit{\psi }}_{\ast }$ . Then for ${\mathit{g}}_{\mathit{\epsilon }}^1+{\mathit{\epsilon }}^{1-\mathit{\gamma }}{\mathit{g}}_{\mathit{\epsilon }}^0\stackrel{\mathit{\text{(16)}}}{:=}{\mathit{\phi }}_{\mathit{\epsilon }}$ one has $$\begin{array}{}\text{(31)}& {\Vert {\mathit{g}}_{\mathit{\epsilon }}^0-{\mathit{g}}_{\ast \mathrm{,}\mathit{\epsilon }}^0\Vert }_{{\mathit{L}}^{\mathit{p}}\mathrm{(}{\mathrm{\Omega }}_{\mathit{\epsilon }}^0\mathrm{)}}+{\Vert \mathit{\epsilon }\nabla {\mathit{g}}_{\mathit{\epsilon }}^0-\mathit{\epsilon }\nabla {\mathit{g}}_{\ast \mathrm{,}\mathit{\epsilon }}^0\Vert }_{{\mathit{L}}^{\mathit{p}}\mathrm{(}{\mathrm{\Omega }}_{\mathit{\epsilon }}^0\mathrm{)}}+{\Vert {\mathit{g}}_{\mathit{\epsilon }}^1-{\mathit{g}}_{\ast \mathrm{,}\mathit{\epsilon }}^1\Vert }_{{\mathit{W}}^{1\mathrm{,}\mathit{p}}\mathrm{(}\mathrm{\Omega }\mathrm{)}}\to 0\mathit{\text{as }}\mathit{\epsilon }\downarrow 0.\end{array}$$

Remark 1.

To illustrate the result of Theorem 2, consider the case $\mathit{\gamma }=1$ with ${\ell }_{\mathit{\epsilon }}:={\mathit{m}}_{\mathit{\epsilon }}\ell \mathrm{(}\mathit{x}\mathrm{,}\frac{\mathit{x}}{\mathit{\epsilon }}\mathrm{)}$ , where $\ell \mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}$ is smooth both in x and y and is periodic in y. If the domain Ω and the pore set ${\mathit{Y}}^0$ are sufficiently regular, the minimisers $\mathrm{(}{\mathit{g}}_{\ast }^0\mathrm{,}{\mathit{g}}_{\ast }^1\mathrm{)}$ and ${\mathit{\psi }}_{\ast }$ are smooth, by the classical elliptic regularity theory, see e.g. [21]. In that case we may set $$\begin{array}{c}{\mathit{g}}_{\ast \mathrm{,}\mathit{\epsilon }}^0\mathrm{(}\mathit{x}\mathrm{)}:={\mathit{g}}_{\ast }^0\mathrm{(}\mathit{x}\mathrm{,}\mathit{x}\mathrm{/}\mathit{\epsilon }\mathrm{)}\text{and}{\mathit{g}}_{\ast \mathrm{,}\mathit{\epsilon }}^1:={\mathit{g}}_{\ast }^1\mathrm{(}\mathit{x}\mathrm{)}+\mathit{\epsilon }{\mathit{\psi }}_{\ast }\mathrm{(}\mathit{x}\mathrm{,}\mathit{x}\mathrm{/}\mathit{\epsilon }\mathrm{)}\mathrm{,}\end{array}$$ and the asymptotic formula for ${\mathit{\phi }}_{\mathit{\epsilon }}$ reads $$\begin{array}{c}{\mathit{\phi }}_{\mathit{\epsilon }}={\mathit{g}}_{\ast }^0\mathrm{(}\mathit{x}\mathrm{,}\mathit{x}\mathrm{/}\mathit{\epsilon }\mathrm{)}+{\mathit{g}}_{\ast }^1\mathrm{(}\mathit{x}\mathrm{)}+\mathit{\epsilon }{\mathit{\psi }}_{\ast }\mathrm{(}\mathit{x}\mathrm{,}\mathit{x}\mathrm{/}\mathit{\epsilon }\mathrm{)}+{\mathit{R}}_{\mathit{\epsilon }}\mathrm{(}\mathit{x}\mathrm{)}\mathrm{,}\end{array}$$ where $\Vert {\mathit{R}}_{\mathit{\epsilon }}{\Vert }_{{\mathit{W}}^{1\mathrm{,}\mathit{p}}\mathrm{(}\mathrm{\Omega }\mathrm{)}}\to 0$ as $\mathit{\epsilon }\downarrow 0$ .

In the remainder of this section, we restrict to the special case of the introduction (in the small strain regime), i.e. we assume that $\mathit{\alpha }\mathrm{>}1$ , $\mathit{\gamma }=1$ , ${\mathit{m}}_{\mathit{\epsilon }}={\mathit{\epsilon }}^{\mathit{\alpha }}$ and ${\ell }_{\mathit{\epsilon }}={\mathit{\epsilon }}^{\mathit{\alpha }}\mathit{f}$ for some $\mathit{f}\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}$ , so that the functionals ${\mathit{I}}_{\mathit{\epsilon }}^{\mathit{\alpha }}$ and ${\mathit{I}}_{\mathit{\epsilon }}$ defined in (2) and (24) are identical. Hence, Theorem 1 and Proposition 1 prove two-scale Γ-convergence of ${\mathit{I}}_{\mathit{\epsilon }}^{\mathit{\alpha }}$ and the convergence of the associated minimisation problems (as claimed in the Introduction), and Theorem 2 yields the two-scale expansion (6). We argue that the functionals ${\mathit{I}}_{\mathit{\epsilon }}^{\mathit{\alpha }}$ Γ-converge to the (single scale) limit ${\overline{\mathit{I}}}_{\text{small}}$ : Proposition 2.

For $\mathit{\alpha }\mathrm{>}1$ and $\mathit{f}\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}$ consider ${\mathit{I}}_{\mathit{\epsilon }}^{\mathit{\alpha }}$ and ${\overline{\mathit{I}}}_{\text{small}}$ defined in ( 2 ) and ( 7 ). We extend ${\mathit{I}}_{\mathit{\epsilon }}^{\mathit{\alpha }}$ to a functional on ${\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}$ by setting ${\mathit{I}}_{\mathit{\epsilon }}^{\mathit{\alpha }}:=+\infty$ on ${\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}\setminus {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ . Then:

(Compactness). Suppose that ${\mathit{\phi }}_{\mathit{\epsilon }}\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}$ satisfy $$\begin{array}{c}\underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}sup}{\mathit{I}}_{\mathit{\epsilon }}^{\mathit{\alpha }}\mathrm{(}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}\mathrm{<}\infty .\end{array}$$ Then, up to a subsequence, we have ${\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{\rightharpoonup }\mathit{\phi }$ weakly in ${\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}$ .

(Lower bound). For every ${\mathit{\phi }}_{\mathit{\epsilon }}\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}$ with ${\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{\rightharpoonup }\mathit{\phi }$ weakly in ${\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}$ we have $$\begin{array}{c}\underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}inf}{\mathit{I}}_{\mathit{\epsilon }}^{\mathit{\alpha }}\mathrm{(}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}⩾{\overline{\mathit{I}}}_{\text{small}}\mathrm{(}\mathit{\phi }\mathrm{)}\end{array}$$

(Upper bound). For every $\mathit{\phi }\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}$ we can find ${\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{\rightharpoonup }\mathit{\phi }$ weakly in ${\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}$ such that $$\begin{array}{c}\underset{\mathit{\epsilon }\downarrow 0}{lim}{\mathit{I}}_{\mathit{\epsilon }}^{\mathit{\alpha }}\mathrm{(}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}={\overline{\mathit{I}}}_{\text{small}}\mathrm{(}\mathit{\phi }\mathrm{)}\end{array}$$

(Convergence of the minimisation problem). Let ${\mathit{\phi }}_{\mathit{\epsilon }}\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ denote an infimizing sequence of ${\mathit{I}}_{\mathit{\epsilon }}^{\mathit{\alpha }}$ , i.e. $$\begin{array}{c}{\mathit{I}}_{\mathit{\epsilon }}^{\mathit{\alpha }}\mathrm{(}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}⩽\underset{\mathit{\phi }\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}}{inf}{\mathit{I}}_{\mathit{\epsilon }}^{\mathit{\alpha }}\mathrm{(}\mathit{\phi }\mathrm{)}+\mathit{o}\mathrm{(}1\mathrm{)}\mathrm{,}\mathit{\epsilon }\downarrow 0.\end{array}$$ Let $\mathrm{(}{\mathit{g}}_{\ast }^0\mathrm{,}{\mathit{g}}_{\ast }^1\mathrm{)}$ denote the unique minimiser of ${\mathit{I}}_{\text{small}}$ and ${\mathit{\phi }}_{\ast }$ denote the unique minimiser of ${\overline{\mathit{I}}}_{\text{small}}$ . Then ${inf}_{{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}}{\mathit{I}}_{\mathit{\epsilon }}^{\mathit{\alpha }}\to {min}_{{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}}{\overline{\mathit{I}}}_{\text{small}}$ , ${\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{\rightharpoonup }{\mathit{\phi }}_{\ast }$ weakly in ${\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}$ , and $$\begin{array}{c}{\mathit{\phi }}_{\ast }={\mathit{g}}_{\ast }^1+{{\displaystyle \int }}_{{\mathit{Y}}_0}{\mathit{g}}^0\mathit{d}\mathit{y}\mathrm{,}{\overline{\mathit{I}}}_{\text{small}}\mathrm{(}{\mathit{\phi }}_{\ast }\mathrm{)}={\mathit{I}}_{\text{small}}\mathrm{(}{\mathit{g}}_{\ast }^0\mathrm{,}{\mathit{g}}_{\ast }^1\mathrm{)}.\end{array}$$

3.3. Convergence in the finite strain regime ${\mathit{m}}_{\mathit{\epsilon }}={\mathit{\epsilon }}^{\mathit{\gamma }}$

Throughout this section we assume that

${\mathit{W}}^1$ satisfies the conditions (W1)–(W3).

${\mathit{W}}^0:{\mathbb{R}}^{\mathit{d}\times \mathit{d}}\to [0\mathrm{,}\infty \mathrm{)}$ is continuous and satisfies the growth condition $$\begin{array}{}\text{(32a)}& {\mathit{c}}_0{dist}^2\mathrm{(}\mathit{F}\mathrm{,}\mathrm{SO}\mathrm{(}\mathit{d}\mathrm{)}\mathrm{)}⩽{\mathit{W}}^0\mathrm{(}\mathit{F}\mathrm{)}⩽{\mathit{c}}_0^{-1}\mathrm{(}1+|\mathit{F}{|}^2\mathrm{)}\forall \mathit{F}\in {\mathbb{R}}^{\mathit{d}\times \mathit{d}}\mathrm{,}\end{array}$$ and the local Lipschitz condition $$\begin{array}{}\text{(32b)}& |{\mathit{W}}^0\mathrm{(}\mathit{F}+\mathit{G}\mathrm{)}-{\mathit{W}}^0\mathrm{(}\mathit{F}\mathrm{)}|⩽{\mathit{c}}_0^{-1}\mathrm{(}1+|\mathit{F}|+|\mathit{G}|\mathrm{)}|\mathit{G}|\forall \mathit{F}\mathrm{,}\mathit{G}\in {\mathbb{R}}^{\mathit{d}\times \mathit{d}}.\end{array}$$

We prove that in the finite strain regime the limit functional $$\begin{array}{c}{\mathcal{E}}_{\text{finite}}:{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{,}{\mathit{H}}_0^1\mathrm{(}{\mathit{Y}}^0\mathrm{\left)}\mathrm{\right)}\times {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}\to [0\mathrm{,}\infty \mathrm{)}\end{array}$$ is given by $$\begin{array}{c}{\mathcal{E}}_{\text{finite}}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}:={\iint }_{\mathrm{\Omega }\times {\mathit{Y}}^0}\mathcal{Q}{\mathit{W}}^0\mathrm{(}\mathit{I}+{\nabla }_{\mathit{y}}{\mathit{g}}^0\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}\mathrm{)}\mathit{d}\mathit{y}\mathit{d}\mathit{x}+{{\displaystyle \int }}_{\mathrm{\Omega }}{\mathit{Q}}_{\mathrm{hom}}^1\mathrm{(}\nabla {\mathit{g}}^1\mathrm{(}\mathit{x}\mathrm{)}\mathrm{)}\mathit{d}\mathit{x}\mathrm{,}\end{array}$$ where $\mathcal{Q}{\mathit{W}}^0$ denotes the quasiconvex envelope of ${\mathit{W}}^0$ (see e.g. [19]). The associated limit of the total energy ${\mathit{I}}^{\mathit{\epsilon }}$ , see (24), is given by (cf. (26)) $$\begin{array}{c}{\mathit{I}}_{\mathrm{finite}}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}:={\mathcal{E}}_{\mathrm{finite}}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}-{{\displaystyle \int }}_{\mathrm{\Omega }}\mathrm{(}{{\displaystyle \int }}_{{\mathit{Y}}^0}{\ell }^0·{\mathit{g}}^0\mathit{d}\mathit{y}+{\ell }^1·{\mathit{g}}^1\mathrm{)}\mathit{d}\mathit{x}\mathrm{,}\end{array}$$ where ${\ell }^0$ , ${\ell }^1$ are defined in the same way as in (25c).

Theorem 3.

(Lower bound). Consider a sequence ${\mathit{\phi }}_{\mathit{\epsilon }}\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ and the associated decomposition ${\mathit{\epsilon }}^{1-\mathit{\gamma }}{\mathit{g}}_{\mathit{\epsilon }}^0+{\mathit{g}}_{\mathit{\epsilon }}^1\stackrel{\mathit{\text{(16)}}}{:=}{\mathit{\phi }}_{\mathit{\epsilon }}$ . If ${\mathit{\phi }}_{\mathit{\epsilon }}\stackrel{2}{\mathrm{\rightharpoonup }}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}$ and ${\mathit{g}}_{\mathit{\epsilon }}^0\stackrel{2}{\to }{\mathit{g}}^0$ , then $$\begin{array}{c}\underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}inf}{\mathit{\epsilon }}^{-2\mathit{\gamma }}{\mathcal{E}}_{\mathit{\epsilon }}\mathrm{\left(}{\mathit{\epsilon }}^{\mathit{\gamma }}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{\right)}⩾{\mathcal{E}}_{\mathrm{finite}}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}.\end{array}$$

(Recovery sequence). For any ${\mathit{g}}^0\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{,}{\mathit{H}}_0^1\mathrm{(}{\mathit{Y}}^0\mathrm{)}\mathrm{)}$ and ${\mathit{g}}^1\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ there exists a sequence ${\mathit{\phi }}_{\mathit{\epsilon }}\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ such that ${\mathit{\phi }}_{\mathit{\epsilon }}\stackrel{2}{\mathrm{\rightharpoonup }}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}$ and $$\begin{array}{c}\underset{\mathit{\epsilon }\downarrow 0}{lim}{\mathit{\epsilon }}^{-2\mathit{\gamma }}{\mathcal{E}}_{\mathit{\epsilon }}\mathrm{\left(}{\mathit{\epsilon }}^{\mathit{\gamma }}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{\right)}={\mathcal{E}}_{\mathrm{finite}}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}.\end{array}$$

Suppose that the force densities ${\ell }_{\mathit{\epsilon }}\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}$ satisfy ( 25a ) and ( 25c ) with ${\mathit{m}}_{\mathit{\epsilon }}={\mathit{\epsilon }}^{\mathit{\gamma }}$ . Then the infima converge, i.e. $$\begin{array}{}\text{(33)}& \underset{\mathit{\epsilon }\downarrow 0}{lim}\underset{\mathit{\phi }\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}}{inf}{\mathit{I}}_{\mathit{\epsilon }}\mathrm{(}\mathit{\phi }\mathrm{)}=inf{\mathit{I}}_{\mathrm{finite}}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}\mathrm{,}\end{array}$$ where the infimum on the right-hand side is taken over all functions ${\mathit{g}}^0\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{,}{\mathit{H}}_0^1\mathrm{(}{\mathit{Y}}^0\mathrm{)}\mathrm{)}$ and ${\mathit{g}}^1\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ . Moreover, there exist a minimising pair $\mathrm{(}{\mathit{g}}_{\ast }^0\mathrm{,}{\mathit{g}}_{\ast }^1\mathrm{)}$ and a recovery sequence ${\mathit{\phi }}_{\mathit{\epsilon }}\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ with ${\mathit{\phi }}_{\mathit{\epsilon }}\stackrel{2}{\to }\mathrm{(}{\mathit{g}}_{\ast }^0\mathrm{,}{\mathit{g}}_{\ast }^1\mathrm{)}$ such that $$\begin{array}{}\text{(34)}& {\mathit{I}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}\to {\mathit{I}}_{\mathrm{finite}}\mathrm{(}{\mathit{g}}_{\ast }^0\mathrm{,}{\mathit{g}}_{\ast }^1\mathrm{)}=min{\mathit{I}}_{\mathrm{finite}}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}\mathit{\text{as }}\mathit{\epsilon }\downarrow 0.\end{array}$$

Remark 2 (Example in Section 1).

If we consider Theorem 3 and Proposition 1 in the case $\mathit{\gamma }=1$ , ${\mathit{m}}_{\mathit{\epsilon }}=\mathit{\epsilon }$ and ${\ell }_{\mathit{\epsilon }}=\mathit{\epsilon }\mathit{f}$ for some $\mathit{f}\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}$ , then we recover the special case (in the finite strain regime) presented in the introduction. In particular, we deduce that the functionals ${\mathit{I}}_{\mathit{\epsilon }}^{\mathit{\alpha }}$ two-scale Γ-converge (in the sense of Theorem 3) to ${\mathit{I}}_{\text{finite}}$ . Arguing as in the small-strain regime, cf. Proposition 2, we deduce that ${\mathit{I}}_{\mathit{\epsilon }}^{\mathit{\alpha }}$ Γ-converges (with respect to the weak topology in ${\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}$ ) to ${\overline{\mathit{I}}}_{\text{finite}}$ , cf. (8). We leave the details to the readers.

4. Proofs

We start by proving the auxiliary results discussed in Section 3.1. Sections 4.2 and 4.3 contain the proofs of the main statements in the small strain and finite strain cases, respectively.

4.1. Proofs of Lemma 1, Lemma 2, and Lemma 3: A priori estimate, compactness and approximation

A key ingredient in the proof of Lemma 1 is the geometric rigidity estimate by Friesecke et al. [23]:

Theorem 4 (Geometric rigidity estimate, see [23]).

Let U be an open, bounded Lipschitz domain in ${\mathbb{R}}^{\mathit{d}}$ , $\mathit{d}⩾2$ . There exists a constant $\mathit{C}\mathrm{(}\mathit{U}\mathrm{)}$ with the following property: for each $\mathit{v}\in {\mathit{H}}^1\mathrm{(}\mathit{U}\mathrm{)}$ there is a rotation $\mathit{R}\in \mathrm{SO}\mathrm{(}\mathit{d}\mathrm{)}$ such that $$\begin{array}{c}{{\displaystyle \int }}_{\mathit{U}}{|\nabla \mathit{v}\mathrm{(}\mathit{x}\mathrm{)}-\mathit{R}|}^2\mathit{d}\mathit{x}⩽\mathit{C}\mathrm{(}\mathit{U}\mathrm{)}{{\displaystyle \int }}_{\mathit{U}}{dist}^2\mathrm{(}\nabla \mathit{v}\mathrm{(}\mathit{x}\mathrm{)}\mathrm{,}\mathrm{SO}\mathrm{(}\mathit{d}\mathrm{)}\mathrm{)}\mathit{d}\mathit{x}.\end{array}$$ Moreover, the constant $\mathit{C}\mathrm{(}\mathit{U}\mathrm{)}$ is invariant under uniform scaling of U.

In fact, we need the following modified version, which is adapted to perforated domains.

Lemma 4.

There exists a constant $\mathit{C}\mathrm{>}0$ that only depends on Ω and ${\mathit{Y}}^1$ such that for all $\mathit{\epsilon }\mathrm{>}0$ and $\mathit{v}\in {\mathit{H}}^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ satisfying $$\begin{array}{}\text{(35)}& {{\displaystyle \int }}_{{\mathrm{\Omega }}_{\mathit{\epsilon }}^0}\nabla \mathit{v}:\nabla \mathit{\zeta }\mathit{d}\mathit{x}=0\forall \mathit{\zeta }\in {\mathit{H}}_0^1\mathrm{\left(}{\mathrm{\Omega }}_{\mathit{\epsilon }}^0\mathrm{\right)}\mathrm{,}\end{array}$$ the estimates $$\begin{array}{}\text{(36)}& {\Vert dist\mathrm{(}\nabla \mathit{v}\mathrm{,}\mathrm{SO}\mathrm{(}\mathit{d}\mathrm{)}\mathrm{)}\Vert }_{{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}}& ⩽& \mathit{C}{\Vert dist\mathrm{(}\nabla \mathit{v}\mathrm{,}\mathrm{SO}\mathrm{(}\mathit{d}\mathrm{)}\mathrm{)}\Vert }_{{\mathit{L}}^2\mathrm{(}{\mathrm{\Omega }}_{\mathit{\epsilon }}^1\mathrm{)}}\mathrm{,}\text{(37)}& \Vert \nabla \mathit{v}-\mathit{R}{\Vert }_{{\mathit{L}}^2\mathrm{(}{\mathrm{\Omega }}_{\mathit{\epsilon }}^1\mathrm{)}}& ⩽& \mathit{C}{\Vert dist\mathrm{(}\nabla \mathit{v}\mathrm{,}\mathrm{SO}\mathrm{(}\mathit{d}\mathrm{)}\mathrm{)}\Vert }_{{\mathit{L}}^2\mathrm{(}{\mathrm{\Omega }}_{\mathit{\epsilon }}^1\mathrm{)}}\mathrm{,}\end{array}$$ hold for some $\mathit{R}\in \mathrm{SO}\mathrm{(}\mathit{d}\mathrm{)}$ , which may depend on v. In addition, if $\mathit{v}\mathrm{(}\mathit{x}\mathrm{)}=\mathit{x}+\mathit{c}$ on $\partial \mathrm{\Omega }$ for some constant c, then we may set $\mathit{R}=\mathit{I}$ .

Proof of Lemma 4.

Step 1. The proof of the inequality (36).

Let ${\hat{\mathrm{\Omega }}}_{\mathit{\epsilon }}:=\bigcup \{\mathit{\epsilon }\mathrm{(}\mathit{\xi }+\mathit{Y}\mathrm{)}|\mathit{\xi }\in {\mathbb{Z}}^{\mathit{d}}\mathrm{,}\mathit{\epsilon }\mathrm{(}\mathit{\xi }+\mathit{Y}\mathrm{)}\subset \mathrm{\Omega }\}$ denote the union of ε-cells that are completely contained in Ω. Since $\mathrm{\Omega }\setminus {\hat{\mathrm{\Omega }}}_{\mathit{\epsilon }}\subset {\mathrm{\Omega }}_{\mathit{\epsilon }}^1$ , it suffices to prove (36) for Ω replaced by ${\hat{\mathrm{\Omega }}}_{\mathit{\epsilon }}$ , respectively. In fact we shall prove the following stronger estimate: for all $\mathit{\xi }\in {\mathbb{Z}}^{\mathit{d}}$ with $\mathit{\epsilon }\mathrm{(}\mathit{\xi }+\mathit{Y}\mathrm{)}\subset \mathrm{\Omega }$ we have $$\begin{array}{}\text{(38)}& {{\displaystyle \int }}_{\mathit{\epsilon }\mathrm{(}\mathit{\xi }+\mathit{Y}\mathrm{)}}{dist}^2\mathrm{(}\nabla \mathit{v}\mathrm{,}\mathrm{SO}\mathrm{(}\mathit{d}\mathrm{)}\mathrm{)}\mathit{d}\mathit{x}\lesssim {{\displaystyle \int }}_{\mathit{\epsilon }\mathrm{(}\mathit{\xi }+{\mathit{Y}}^1\mathrm{)}}{dist}^2\mathrm{(}\nabla \mathit{v}\mathrm{,}\mathrm{SO}\mathrm{(}\mathit{d}\mathrm{)}\mathrm{)}\mathit{d}\mathit{x}.\end{array}$$ For the argument fix an admissible $\mathit{\xi }\in {\mathbb{Z}}^{\mathit{d}}$ . Application of Theorem 4 with $\mathit{U}=\mathit{\epsilon }\mathrm{(}\mathit{\xi }+{\mathit{Y}}^1\mathrm{)}$ yields a rotation $\mathit{R}\in \mathrm{SO}\mathrm{(}\mathit{d}\mathrm{)}$ such that $$\begin{array}{}\text{(39)}& {{\displaystyle \int }}_{\mathit{\epsilon }\mathrm{(}\mathit{\xi }+{\mathit{Y}}^1\mathrm{)}}|\nabla \mathit{v}-\mathit{R}{|}^2\mathit{d}\mathit{x}\lesssim {{\displaystyle \int }}_{\mathit{\epsilon }\mathrm{(}\mathit{\xi }+{\mathit{Y}}^1\mathrm{)}}{dist}^2\mathrm{(}\nabla \mathit{v}\mathrm{,}\mathrm{SO}\mathrm{(}\mathit{d}\mathrm{)}\mathrm{)}\mathit{d}\mathit{x}.\end{array}$$ Note that the multiplicative constant in the estimate above only depends on ${\mathit{Y}}^1$ , since $\mathit{\epsilon }\mathrm{(}\mathit{\xi }+{\mathit{Y}}^1\mathrm{)}$ is a dilation and translation of ${\mathit{Y}}^1$ . On the other hand, since $\mathit{\epsilon }\mathrm{(}\mathit{\xi }+{\mathit{Y}}^0\mathrm{)}\ni \mathit{x}\mapsto \mathrm{(}\mathit{v}\mathrm{(}\mathit{x}\mathrm{)}-\mathit{R}\mathit{x}\mathrm{)}$ is harmonic, we have (cf. (11)): $$\begin{array}{c}{{\displaystyle \int }}_{\mathit{\epsilon }\mathrm{(}\mathit{\xi }+\mathit{Y}\mathrm{)}}{dist}^2\mathrm{(}\nabla \mathit{v}\mathrm{,}\mathrm{SO}\mathrm{(}\mathit{d}\mathrm{)}\mathrm{)}\mathit{d}\mathit{x}⩽{{\displaystyle \int }}_{\mathit{\epsilon }\mathrm{(}\mathit{\xi }+\mathit{Y}\mathrm{)}}|\nabla \mathit{v}-\mathit{R}{|}^2\lesssim {{\displaystyle \int }}_{\mathit{\epsilon }\mathrm{(}\mathit{\xi }+{\mathit{Y}}^1\mathrm{)}}|\nabla \mathit{v}-\mathit{R}{|}^2.\end{array}$$ Combined with (39), inequality (38) follows.

Step 2. The proof of the rigidity estimate (37).

From (36) and Theorem 4 (applied with $\mathit{U}=\mathrm{\Omega }$ ) we deduce that for some $\mathit{R}\in \mathrm{SO}\mathrm{(}\mathit{d}\mathrm{)}$ : $$\begin{array}{}\text{(40)}& \Vert \nabla \mathit{v}-\mathit{R}{\Vert }_{{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}}\lesssim {\Vert dist\mathrm{(}\nabla \mathit{v}\mathrm{,}\mathrm{SO}\mathrm{(}\mathit{d}\mathrm{)}\mathrm{)}\Vert }_{{\mathit{L}}^2\mathrm{(}{\mathrm{\Omega }}_{\mathit{\epsilon }}^1\mathrm{)}}\mathrm{,}\end{array}$$ which in particular implies (37). Finally we argue that one can set $\mathit{R}=\mathit{I}$ , if $\mathit{v}=\mathit{x}+\mathit{c}$ on $\partial \mathrm{\Omega }$ . In view of (40), it suffices to show that ${\int }_{\mathrm{\Omega }}|\nabla \mathit{v}-\mathit{I}{|}^2\mathit{d}\mathit{x}⩽{\int }_{\mathrm{\Omega }}|\nabla \mathit{v}-\mathit{R}{|}^2\mathit{d}\mathit{x}$ for all $\mathit{R}\in \mathrm{SO}\mathrm{(}\mathit{d}\mathrm{)}$ . This inequality can be seen as follows: Consider $\mathit{\phi }\mathrm{(}\mathit{x}\mathrm{)}:=\mathit{v}\mathrm{(}\mathit{x}\mathrm{)}-\mathit{x}-\mathit{c}$ and note that φ vanishes on $\partial \mathrm{\Omega }$ , so that $$\begin{array}{rcl}{{\displaystyle \int }}_{\mathrm{\Omega }}|\nabla \mathit{v}-\mathit{I}{|}^2\mathit{d}\mathit{x}& =& {{\displaystyle \int }}_{\mathrm{\Omega }}|\nabla \mathit{\phi }{|}^2\mathit{d}\mathit{x}⩽{{\displaystyle \int }}_{\mathrm{\Omega }}|\nabla \mathit{\phi }{|}^2+|\mathit{I}-\mathit{R}{|}^2\mathit{d}\mathit{x}={{\displaystyle \int }}_{\mathrm{\Omega }}{|\nabla \mathit{\phi }+\mathrm{(}\mathit{I}-\mathit{R}\mathrm{)}|}^2\mathit{d}\mathit{x}\\ & =& {{\displaystyle \int }}_{\mathrm{\Omega }}|\nabla \mathit{v}-\mathit{R}{|}^2\mathit{d}\mathit{x}\mathrm{,}\end{array}$$ which in fact holds for an arbitrary matrix R. □

We are now in position to present the proofs of Lemmas 1 and 2.

Proof of Lemma 1.

In the following, the symbol ≲ stands for ⩽ up to a multiplicative constant that only depends on ${\mathit{Y}}^1$ and Ω.

Step 1. Existence of the decomposition (16) and derivation of the estimate for ${\mathit{g}}_{\mathit{\epsilon }}^1$ .

Let ${\mathit{g}}_{\mathit{\epsilon }}^1$ denote the unique function in ${\mathit{H}}^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ characterised by ${\mathit{g}}_{\mathit{\epsilon }}^1={\mathit{\phi }}_{\mathit{\epsilon }}$ in ${\mathrm{\Omega }}_{\mathit{\epsilon }}^1$ and (16)(ii). Since $\partial {\mathrm{\Omega }}_{\mathit{\epsilon }}^0$ is Lipschitz, we deduce that ${\mathit{g}}_{\mathit{\epsilon }}^0:={\mathit{\epsilon }}^{\mathit{\gamma }-1}\mathrm{(}{\mathit{\phi }}_{\mathit{\epsilon }}-{\mathit{g}}_{\mathit{\epsilon }}^1\mathrm{)}\in {\mathit{H}}_0^1\mathrm{(}{\mathrm{\Omega }}_{\mathit{\epsilon }}^0\mathrm{)}$ . This proves the existence of the decomposition. We claim that $$\begin{array}{}\text{(41)}& {{\displaystyle \int }}_{{\mathrm{\Omega }}_{\mathit{\epsilon }}^0}{|\nabla {\mathit{g}}_{\mathit{\epsilon }}^1|}^2\mathit{d}\mathit{x}\lesssim {{\displaystyle \int }}_{{\mathrm{\Omega }}_{\mathit{\epsilon }}^1}|\nabla {\mathit{\phi }}_{\mathit{\epsilon }}{|}^2\mathit{d}\mathit{x}={{\displaystyle \int }}_{{\mathrm{\Omega }}_{\mathit{\epsilon }}^1}{|\nabla {\mathit{g}}_{\mathit{\epsilon }}^1|}^2\mathit{d}\mathit{x}.\end{array}$$ Since ${\mathrm{\Omega }}_{\mathit{\epsilon }}^0$ is defined as the union of the sets $\mathit{\epsilon }\mathrm{(}\mathit{\xi }+{\mathit{Y}}^0\mathrm{)}$ with $\mathit{\xi }\in {\mathit{Z}}_{\mathit{\epsilon }}:=\{\mathit{\xi }\in {\mathbb{Z}}^{\mathit{d}}:\mathit{\epsilon }\mathrm{(}\mathit{\xi }+\mathit{Y}\mathrm{)}\subset \mathrm{\Omega }\}$ , it suffices to prove ${\int }_{\mathit{\epsilon }\mathrm{(}\mathit{\xi }+{\mathit{Y}}^0\mathrm{)}}|\nabla {\mathit{g}}_{\mathit{\epsilon }}^1{|}^2\mathit{d}\mathit{x}\lesssim {\int }_{\mathit{\epsilon }\mathrm{(}\mathit{\xi }+{\mathit{Y}}^1\mathrm{)}}|\nabla {\mathit{\phi }}_{\mathit{\epsilon }}{|}^2\mathit{d}\mathit{x}$ . The latter follows from (11) by a scaling argument, since the rescaled functions $\mathit{y}\mapsto {\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{(}\mathit{\epsilon }\mathrm{(}\mathit{\xi }+\mathit{y}\mathrm{)}\mathrm{)}$ and $\mathit{y}\mapsto {\mathit{g}}_{\mathit{\epsilon }}^1\mathrm{(}\mathit{\epsilon }\mathrm{(}\mathit{\xi }+\mathit{y}\mathrm{)}\mathrm{)}$ satisfy (10).

Next, we prove (17). Consider ${\mathit{v}}_{\mathit{\epsilon }}\mathrm{(}\mathit{x}\mathrm{)}:=\mathit{x}+{\mathit{g}}_{\mathit{\epsilon }}^1\mathrm{(}\mathit{x}\mathrm{)}$ and note that ${\mathit{v}}_{\mathit{\epsilon }}\mathrm{(}\mathit{x}\mathrm{)}$ satisfies (35). Hence, (41) and Lemma 4 yield $$\begin{array}{rcl}{{\displaystyle \int }}_{\mathrm{\Omega }}{|\nabla {\mathit{g}}_{\mathit{\epsilon }}^1|}^2\mathit{d}\mathit{x}& \lesssim & {{\displaystyle \int }}_{{\mathrm{\Omega }}_{\mathit{\epsilon }}^1}{|\nabla {\mathit{g}}_{\mathit{\epsilon }}^1|}^2\mathit{d}\mathit{x}⩽{{\displaystyle \int }}_{\mathrm{\Omega }}|\nabla {\mathit{v}}_{\mathit{\epsilon }}-\mathit{I}{|}^2\stackrel{\text{(37)}}{⩽}\mathit{C}{{\displaystyle \int }}_{{\mathrm{\Omega }}_{\mathit{\epsilon }}^1}{dist}^2\mathrm{(}\mathit{I}+\nabla {\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{(}\mathit{x}\mathrm{)}\mathrm{,}\mathrm{SO}\mathrm{(}\mathit{d}\mathrm{)}\mathrm{)}\mathit{d}\mathit{x}\\ \text{(42)}& & ⩽& {\mathrm{\Phi }}_{\mathit{\epsilon }}^{\mathit{\gamma }}\mathrm{(}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}.\end{array}$$ Since ${\mathit{g}}_{\mathit{\epsilon }}^1$ vanishes on the boundary of Ω, the estimate upgrades (by Poincaré’s inequality) to $\Vert {\mathit{g}}_{\mathit{\epsilon }}^1{\Vert }_{{\mathit{H}}^1\mathrm{(}\mathrm{\Omega }\mathrm{)}}⩽\mathit{C}{\mathrm{\Phi }}_{\mathit{\epsilon }}^{\mathit{\gamma }}\mathrm{(}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}$ .

Step 2. Derivation of the estimate for ${\mathit{g}}_{\mathit{\epsilon }}^0$ .

Since we have an improved Poincaré inequality (see e.g. [25, Lemma 1.6]): $$\begin{array}{}\text{(43)}& \forall \mathit{g}\in {\mathit{H}}_0^1\mathrm{\left(}{\mathrm{\Omega }}_{\mathit{\epsilon }}^0\mathrm{\right)}:\Vert \mathit{g}{\Vert }_{{\mathit{L}}^2\mathrm{(}{\mathrm{\Omega }}_{\mathit{\epsilon }}^0\mathrm{)}}\lesssim \mathit{\epsilon }\Vert \nabla \mathit{g}{\Vert }_{{\mathit{L}}^2\mathrm{(}{\mathrm{\Omega }}_{\mathit{\epsilon }}^0\mathrm{)}}\mathrm{,}\end{array}$$ it suffices to prove $$\begin{array}{c}{\Vert \mathit{\epsilon }\nabla {\mathit{g}}_{\mathit{\epsilon }}^0\Vert }_{{\mathit{L}}^2\mathrm{(}{\mathrm{\Omega }}_{\mathit{\epsilon }}^0\mathrm{)}}^2\lesssim {\mathrm{\Phi }}_{\mathit{\epsilon }}^{\mathit{\gamma }}\mathrm{(}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}.\end{array}$$ To this end, notice that since ${\mathit{\phi }}_{\mathit{\epsilon }}$ vanishes on the boundary of Ω, we have $$\begin{array}{rcl}\Vert \nabla {\mathit{\phi }}_{\mathit{\epsilon }}{\Vert }_{{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}}^2& =& \underset{\mathit{R}\in \mathrm{SO}\mathrm{(}\mathit{d}\mathrm{)}}{min}\Vert \mathit{I}+\nabla {\mathit{\phi }}_{\mathit{\epsilon }}-\mathit{R}{\Vert }_{{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}}^2\stackrel{\text{Theorem 4}}{\lesssim }{{\displaystyle \int }}_{\mathrm{\Omega }}{dist}^2\mathrm{(}\mathit{I}+\nabla {\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{(}\mathit{x}\mathrm{)}\mathrm{,}\mathrm{SO}\mathrm{(}\mathit{d}\mathrm{)}\mathrm{)}\mathit{d}\mathit{x}\\ \text{(44)}& & ⩽& {\mathit{\epsilon }}^{-2\mathit{\gamma }}{\mathrm{\Phi }}_{\mathit{\epsilon }}^{\mathit{\gamma }}\mathrm{(}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}.\end{array}$$ Thanks to the first identity in (16), we get by triangle inequality: $$\begin{array}{c}{\Vert {\mathit{\epsilon }}^{1-\mathit{\gamma }}\nabla {\mathit{g}}_{\mathit{\epsilon }}^0\Vert }_{{\mathit{L}}^2\mathrm{(}{\mathrm{\Omega }}_{\mathit{\epsilon }}^0\mathrm{)}}⩽\Vert \nabla {\mathit{\phi }}_{\mathit{\epsilon }}{\Vert }_{{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}}+{\Vert \nabla {\mathit{g}}_{\mathit{\epsilon }}^1\Vert }_{{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}}.\end{array}$$ Combined with (42) and (44) we finally get $$\begin{array}{c}{\Vert \mathit{\epsilon }\nabla {\mathit{g}}_{\mathit{\epsilon }}^0\Vert }_{{\mathit{L}}^2\mathrm{(}{\mathrm{\Omega }}_{\mathit{\epsilon }}^0\mathrm{)}}^2={\mathit{\epsilon }}^{2\mathit{\gamma }}{\Vert {\mathit{\epsilon }}^{1-\mathit{\gamma }}\nabla {\mathit{g}}_{\mathit{\epsilon }}^0\Vert }_{{\mathit{L}}^2\mathrm{(}{\mathrm{\Omega }}_{\mathit{\epsilon }}^0\mathrm{)}}^2\lesssim {\mathrm{\Phi }}_{\mathit{\epsilon }}^{\mathit{\gamma }}\mathrm{(}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}.\end{array}$$ □

Proof of Lemma 2.

Step 1. A priori estimate and basic compactness.

From Lemma 1 we deduce that $$\begin{array}{}\text{(45)}& \underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}sup}\mathrm{(}{\Vert {\mathit{g}}_{\mathit{\epsilon }}^0\Vert }_{{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}}^2+{\Vert \mathit{\epsilon }\nabla {\mathit{g}}_{\mathit{\epsilon }}^0\Vert }_{{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}}^2+{\Vert {\mathit{g}}_{\mathit{\epsilon }}^1\Vert }_{{\mathit{H}}^1\mathrm{(}\mathrm{\Omega }\mathrm{)}}^2\mathrm{)}⩽\mathit{C}\underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}sup}{\mathit{m}}_{\mathit{\epsilon }}^{-2}{\mathrm{\Phi }}_{\mathit{\epsilon }}^{\mathit{\gamma }}\mathrm{(}{\mathit{m}}_{\mathit{\epsilon }}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}\mathrm{<}\infty .\end{array}$$ Hence, by standard results concerning two-scale convergence (cf. [1, Proposition 1.14] and [36, Proposition 4.2]), there exist ${\mathit{g}}^1\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ , $\mathit{\psi }\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{,}{\mathit{H}}_{\mathrm{\#}}^1\mathrm{)}$ and ${\mathit{g}}^0\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{,}{\mathit{H}}_{\mathrm{\#}}^1\mathrm{)}$ such that, up to a subsequence, one has $$\begin{array}{ll}\hfill & {\mathit{g}}_{\mathit{\epsilon }}^1\mathrm{\rightharpoonup }{\mathit{g}}^1\text{weakly in }{\mathit{H}}^1\mathrm{(}\mathrm{\Omega }\mathrm{)}\mathrm{,}\nabla {\mathit{g}}_{\mathit{\epsilon }}^1\stackrel{2}{\mathrm{\rightharpoonup }}\nabla {\mathit{g}}^1+{\nabla }_{\mathit{y}}\mathit{\psi }\mathrm{,}\\ \hfill & {\mathit{g}}_{\mathit{\epsilon }}^0\stackrel{2}{\mathrm{\rightharpoonup }}{\mathit{g}}^0\mathrm{,}\mathit{\epsilon }\nabla {\mathit{g}}_{\mathit{\epsilon }}^0\stackrel{2}{\mathrm{\rightharpoonup }}{\nabla }_{\mathit{y}}{\mathit{g}}^0.\end{array}$$

Step 2. The proof of the inclusion $\mathit{\psi }\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{,}\mathcal{A}\mathrm{(}{\mathit{Y}}^0\mathrm{)}\mathrm{)}$ .

By a density argument, it suffices to show that $$\begin{array}{}\text{(46)}& {\iint }_{\mathrm{\Omega }\times {\mathit{Y}}^0}{\nabla }_{\mathit{y}}\mathit{\psi }\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}:{\nabla }_{\mathit{y}}\mathrm{(}{\mathit{\zeta }}_1\mathrm{(}\mathit{x}\mathrm{)}{\mathit{\zeta }}_2\mathrm{(}\mathit{y}\mathrm{)}\mathrm{)}\mathit{d}\mathit{x}\mathit{d}\mathit{y}=0\end{array}$$ for all scalar functions ${\mathit{\zeta }}_1\in {\mathit{C}}_{\mathrm{c}}^{\infty }\mathrm{(}\mathrm{\Omega }\mathrm{)}$ , and all ${\mathit{\zeta }}_2\in {\mathit{C}}_{\mathrm{c}}^{\infty }\mathrm{(}{\mathit{Y}}^0\mathrm{)}$ . To this end, we identify ${\mathit{\zeta }}_2$ with its unique Y-periodic extension to ${\mathbb{R}}^{\mathit{d}}$ that vanishes on ${\mathit{Y}}^1$ , and set $$\begin{array}{c}{\mathit{\zeta }}_{\mathit{\epsilon }}\mathrm{(}\mathit{x}\mathrm{)}:=\mathit{\epsilon }{\mathit{\zeta }}_1\mathrm{(}\mathit{x}\mathrm{)}{\mathit{\zeta }}_2\mathrm{(}\mathit{x}\mathrm{/}\mathit{\epsilon }\mathrm{)}\mathrm{,}\mathit{x}\in \mathrm{\Omega }.\end{array}$$ Thanks to (16) we have $$\begin{array}{c}{{\displaystyle \int }}_{\mathrm{\Omega }}\nabla {\mathit{g}}_{\mathit{\epsilon }}^1:\nabla {\mathit{\zeta }}_{\mathit{\epsilon }}\mathit{d}\mathit{x}=0.\end{array}$$ As can be easily checked, we have $\nabla {\mathit{\zeta }}_{\mathit{\epsilon }}\stackrel{2}{\to }{\mathit{\zeta }}_1\mathrm{(}\mathit{x}\mathrm{)}{\nabla }_{\mathit{y}}{\mathit{\zeta }}_2\mathrm{(}\mathit{y}\mathrm{)}$ , so that $$\begin{array}{rcl}0& =& \underset{\mathit{\epsilon }\downarrow 0}{lim}{{\displaystyle \int }}_{\mathrm{\Omega }}\nabla {\mathit{g}}_{\mathit{\epsilon }}^1:\nabla {\mathit{\zeta }}_{\mathit{\epsilon }}\mathit{d}\mathit{x}\\ & =& {\iint }_{\mathrm{\Omega }\times {\mathit{Y}}^0}\mathrm{(}\nabla {\mathit{g}}^1\mathrm{(}\mathit{x}\mathrm{)}+{\nabla }_{\mathit{y}}\mathit{\psi }\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}\mathrm{)}:\mathrm{(}{\mathit{\zeta }}_1\mathrm{(}\mathit{x}\mathrm{)}{\nabla }_{\mathit{y}}{\mathit{\zeta }}_2\mathrm{(}\mathit{y}\mathrm{)}\mathrm{)}\mathit{d}\mathit{x}\mathit{d}\mathit{y}\\ & =& {\iint }_{\mathrm{\Omega }\times {\mathit{Y}}^0}{\nabla }_{\mathit{y}}\mathit{\psi }\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}:\mathrm{(}{\mathit{\zeta }}_1\mathrm{(}\mathit{x}\mathrm{)}{\nabla }_{\mathit{y}}{\mathit{\zeta }}_2\mathrm{(}\mathit{y}\mathrm{)}\mathrm{)}\mathit{d}\mathit{x}\mathit{d}\mathit{y}\mathrm{,}\end{array}$$ where the last identity holds thanks to the periodicity of ${\mathit{\zeta }}_2$ . This proves (46).

Step 3. The proof of the inclusion ${\mathit{g}}^0\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{,}{\mathit{H}}_0^1\mathrm{(}{\mathit{Y}}^0\mathrm{)}\mathrm{)}$ .

By a density argument, it suffices to show that $$\begin{array}{}\text{(47)}& {\iint }_{\mathrm{\Omega }\times \mathit{Y}}{\mathit{g}}^0\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}·\mathrm{(}{\mathit{\zeta }}_1\mathrm{(}\mathit{x}\mathrm{)}{\mathit{\zeta }}_2\mathrm{(}\mathit{y}\mathrm{)}\mathrm{)}\mathit{d}\mathit{x}\mathit{d}\mathit{y}=0\end{array}$$ for all scalar functions ${\mathit{\zeta }}_1\in {\mathit{C}}_{\mathrm{c}}^{\infty }\mathrm{(}\mathrm{\Omega }\mathrm{)}$ and all ${\mathit{\zeta }}_2\in {\mathit{H}}_{\mathrm{\#}}^1$ with ${\mathit{\zeta }}_2=0$ on ${\mathit{Y}}^0$ . We argue by considering the function ${\mathit{\zeta }}_{\mathit{\epsilon }}\mathrm{(}\mathit{x}\mathrm{)}:={\mathit{\zeta }}_1\mathrm{(}\mathit{x}\mathrm{)}{\mathit{\zeta }}_2\mathrm{(}\frac{\mathit{x}}{\mathit{\epsilon }}\mathrm{)}$ , $\mathit{x}\in \mathrm{\Omega }$ , the support of which is contained in ${\mathrm{\Omega }}_{\mathit{\epsilon }}^1$ for $\mathit{\epsilon }\ll 1$ . Since ${\mathit{\zeta }}_{\mathit{\epsilon }}\stackrel{2}{\to }{\mathit{\zeta }}_1\mathrm{(}\mathit{x}\mathrm{)}\mathit{\zeta }\mathrm{(}\mathit{y}\mathrm{)}$ , and since ${\mathit{g}}_{\mathit{\epsilon }}^0$ is supported in ${\mathrm{\Omega }}_{\mathit{\epsilon }}^0$ , we deduce that $$\begin{array}{c}0=\underset{\mathit{\epsilon }\downarrow 0}{lim}{{\displaystyle \int }}_{\mathrm{\Omega }}{\mathit{g}}_{\mathit{\epsilon }}^0\mathrm{(}\mathit{x}\mathrm{)}·{\mathit{\zeta }}_{\mathit{\epsilon }}\mathrm{(}\mathit{x}\mathrm{)}\mathit{d}\mathit{x}={\iint }_{\mathrm{\Omega }\times \mathit{Y}}{\mathit{g}}^0\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}·\mathrm{(}{\mathit{\zeta }}_1\mathrm{(}\mathit{x}\mathrm{)}{\mathit{\zeta }}_2\mathrm{(}\mathit{y}\mathrm{)}\mathrm{)}\mathit{d}\mathit{x}\mathit{d}\mathit{y}.\end{array}$$ This completes the argument. □

In the proof of Lemma 3 we appeal to the construction of a diagonal sequence that is due to Attouch, see [2]: Lemma 5.

For any ${\mathit{h}:[0\mathrm{,}\infty \mathrm{)}}^2\to [0\mathrm{,}+\infty ]$ , there exists a mapping $\mathrm{(}0\mathrm{,}1\mathrm{)}\ni \mathit{\epsilon }\mapsto \mathit{\delta }\mathrm{(}\mathit{\epsilon }\mathrm{)}\in \mathrm{(}0\mathrm{,}1\mathrm{)}$ such that $$\begin{array}{c}\underset{\mathit{\epsilon }\downarrow 0}{lim}\mathit{\delta }\mathrm{(}\mathit{\epsilon }\mathrm{)}=0\mathit{\text{and}}\underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}sup}\mathit{h}\mathrm{(}\mathit{\epsilon }\mathrm{,}\mathit{\delta }\mathrm{(}\mathit{\epsilon }\mathrm{)}\mathrm{)}⩽\underset{\mathit{\delta }\downarrow 0}{lim\hspace{0.17em}sup}\underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}sup}\mathit{h}\mathrm{(}\mathit{\epsilon }\mathrm{,}\mathit{\delta }\mathrm{)}.\end{array}$$

Proof of Lemma 3.

Step 1. Characterisation of strong two-scale convergence via unfolding.

For ${\mathit{f}}_{\mathit{\epsilon }}:\mathrm{\Omega }\to \mathbb{R}$ and $\mathit{f}:\mathrm{\Omega }\times \mathit{Y}\to \mathbb{R}$ define $$\begin{array}{c}{\mathit{d}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{f}}_{\mathit{\epsilon }}\mathrm{,}\mathit{f}\mathrm{)}:={{\displaystyle \int }}_{{\mathbb{R}}^{\mathit{d}}}{{\displaystyle \int }}_{\mathit{Y}}{\left|{\tilde{\mathit{f}}}_{\mathit{\epsilon }}\mathrm{\right(}\mathit{\epsilon }\lfloor \mathit{x}\mathrm{/}\mathit{\epsilon }\rfloor +\mathit{\epsilon }\mathit{y}\mathrm{)}-\tilde{\mathit{f}}\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}|}^2\mathit{d}\mathit{y}\mathit{d}\mathit{x}\mathrm{,}\end{array}$$ where ${\tilde{\mathit{f}}}_{\mathit{\epsilon }}$ denotes the extension by zero of ${\mathit{f}}_{\mathit{\epsilon }}$ to ${\mathbb{R}}^{\mathit{d}}$ , $\tilde{\mathit{f}}$ denotes the extension by zero of f to ${\mathbb{R}}^{\mathit{d}}\times \mathit{Y}$ , and $\lfloor \mathit{z}\rfloor$ denotes the unique element in ${\mathbb{Z}}^{\mathit{d}}$ with ${\mathit{z}-\lfloor \mathit{z}\rfloor \in [0\mathrm{,}1\mathrm{)}}^{\mathit{d}}$ . We recall from [36] that $$\begin{array}{}\text{(48)}& {\mathit{f}}_{\mathit{\epsilon }}\stackrel{2}{\to }\mathit{f}\iff {\mathit{d}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{f}}_{\mathit{\epsilon }}\mathrm{,}\mathit{f}\mathrm{)}\to 0.\end{array}$$ The characterisation extends in the obvious way to vector-valued functions.

Step 2. Construction of ${\mathit{g}}_{\mathit{\epsilon }}^1$ .

We claim that there exists a sequence ${\mathit{g}}_{\mathit{\epsilon }}^1$ in ${\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ whose elements satisfy (16)(ii) and $$\begin{array}{}\text{(49)}& {\mathit{g}}_{\mathit{\epsilon }}^1\mathrm{\rightharpoonup }{\mathit{g}}^1\text{weakly in }{\mathit{H}}^1\mathrm{(}\mathrm{\Omega }\mathrm{)}\mathrm{,}\nabla {\mathit{g}}_{\mathit{\epsilon }}^1\stackrel{2}{\to }{\mathit{g}}^1+{\nabla }_{\mathit{y}}\mathit{\psi }\mathrm{,}\underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}sup}{\mathit{c}}_{\mathit{\epsilon }}{\Vert \nabla {\mathit{g}}_{\mathit{\epsilon }}^1\Vert }_{{\mathit{L}}^{\infty }\mathrm{(}\mathrm{\Omega }\mathrm{)}}=0.\end{array}$$ Indeed, by a density argument there exist ${\mathit{g}}^{1\mathrm{,}\mathit{\delta }}\in {\mathit{C}}_{\mathrm{c}}^{\infty }\mathrm{(}\mathrm{\Omega }\mathrm{)}$ , $\mathit{\delta }\in \mathrm{(}0\mathrm{,}1\mathrm{)}$ , and $\mathit{\psi }\in {\mathit{C}}_{\mathrm{c}}^{\infty }\mathrm{(}\mathrm{\Omega }\mathrm{,}{\mathit{C}}_{\mathrm{\#}}^{\infty }\mathrm{(}\mathit{Y}\mathrm{)}\mathrm{)}$ such that $$\begin{array}{c}{\Vert {\mathit{g}}^{1\mathrm{,}\mathit{\delta }}-{\mathit{g}}^1\Vert }_{{\mathit{H}}^1\mathrm{(}\mathrm{\Omega }\mathrm{)}}+{\Vert {\nabla }_{\mathit{y}}{\mathit{\psi }}^{\mathit{\delta }}-{\nabla }_{\mathit{y}}\mathit{\psi }\Vert }_{{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\times \mathit{Y}\mathrm{)}}⩽\mathit{\delta }\forall \mathit{\delta }.\end{array}$$ For $\mathit{\epsilon }\mathrm{>}0$ , $\mathit{\delta }\in \mathrm{(}0\mathrm{,}1\mathrm{)}$ , define $$\begin{array}{c}{\mathit{g}}_{\mathit{\epsilon }}^{1\mathrm{,}\mathit{\delta }}\mathrm{(}\mathit{x}\mathrm{)}:={\mathit{g}}^{1\mathrm{,}\mathit{\delta }}\mathrm{(}\mathit{x}\mathrm{)}+\mathit{\epsilon }{\mathit{\psi }}^{\mathit{\delta }}\mathrm{(}\mathit{x}\mathrm{,}\mathit{x}\mathrm{/}\mathit{\epsilon }\mathrm{)}\mathrm{,}\end{array}$$ and set $$\begin{array}{c}{\mathit{d}}_{\mathit{\epsilon }}^{\mathit{\delta }}:={\mathit{d}}_{\mathit{\epsilon }}\mathrm{(}\nabla {\mathit{g}}_{\mathit{\epsilon }}^{1\mathrm{,}\mathit{\delta }}\mathrm{,}\nabla {\mathit{g}}^1+{\nabla }_{\mathit{y}}\mathit{\psi }\mathrm{)}+{\Vert {\mathit{g}}_{\mathit{\epsilon }}^{1\mathrm{,}\mathit{\delta }}-{\mathit{g}}^1\Vert }_{{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}}+{\mathit{c}}_{\mathit{\epsilon }}{\Vert \nabla {\mathit{g}}_{\mathit{\epsilon }}^{1\mathrm{,}\mathit{\delta }}\Vert }_{{\mathit{L}}^{\infty }\mathrm{(}\mathrm{\Omega }\mathrm{)}}.\end{array}$$ By construction, we have ${lim}_{\mathit{\delta }\downarrow 0}{lim\hspace{0.17em}sup}_{\mathit{\epsilon }\downarrow 0}{\mathit{d}}_{\mathit{\epsilon }}^{\mathit{\delta }}=0$ , and Lemma 5 yields a function $\mathit{\epsilon }\mapsto \mathit{\delta }\mathrm{(}\mathit{\epsilon }\mathrm{)}$ with ${lim}_{\mathit{\epsilon }\downarrow 0}{\mathit{d}}_{\mathit{\epsilon }}^{\mathit{\delta }\mathrm{(}\mathit{\epsilon }\mathrm{)}}=0$ . In view of Step 1, this implies that the diagonal sequence ${\tilde{\mathit{g}}}_{\mathit{\epsilon }}^1:={\mathit{g}}_{\mathit{\epsilon }}^{1\mathrm{,}\mathit{\delta }\mathrm{(}\mathit{\epsilon }\mathrm{)}}$ satisfies (49). Now, for each $\mathit{\epsilon }\mathrm{>}0$ , let ${\mathit{g}}_{\mathit{\epsilon }}^1$ denote the function satisfying (16)(ii) and such that ${\mathit{g}}_{\mathit{\epsilon }}^1={\tilde{\mathit{g}}}_{\mathit{\epsilon }}^1$ on ${\mathrm{\Omega }}_{\mathit{\epsilon }}^1$ . To conclude the argument, we only need to show that ${\mathit{g}}_{\mathit{\epsilon }}^1$ satisfies (49). Consider the difference ${\mathit{\eta }}_{\mathit{\epsilon }}:={\tilde{\mathit{g}}}_{\mathit{\epsilon }}^1-{\mathit{g}}_{\mathit{\epsilon }}^1$ . Since ${\mathit{\eta }}_{\mathit{\epsilon }}$ is bounded in ${\mathit{H}}^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ and ${\mathit{\eta }}_{\mathit{\epsilon }}=0$ in ${\mathrm{\Omega }}_{\mathit{\epsilon }}^1$ , we have ${\mathit{\eta }}_{\mathit{\epsilon }}\mathrm{\rightharpoonup }0$ in ${\mathit{H}}^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ , and, up to a subsequence, $\nabla {\mathit{\eta }}_{\mathit{\epsilon }}\stackrel{2}{\mathrm{\rightharpoonup }}{\nabla }_{\mathit{y}}\mathit{\phi }$ for some $\mathit{\phi }\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{,}{\mathit{H}}_0^1\mathrm{(}{\mathit{Y}}^0\mathrm{)}\mathrm{)}$ . On the other hand, since ${\mathit{g}}_{\mathit{\epsilon }}^1$ satisfies (16)(ii) and ${\tilde{\mathit{g}}}_{\mathit{\epsilon }}^1$ satisfies (49), we deduce that $$\begin{array}{r}\Vert \nabla {\mathit{\eta }}_{\mathit{\epsilon }}{\Vert }_{{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}}^2={{\displaystyle \int }}_{{\mathrm{\Omega }}_{\mathit{\epsilon }}^0}\nabla {\mathit{\eta }}_{\mathit{\epsilon }}:\nabla {\mathit{\eta }}_{\mathit{\epsilon }}\mathit{d}\mathit{x}={{\displaystyle \int }}_{{\mathrm{\Omega }}_{\mathit{\epsilon }}^0}\nabla {\tilde{\mathit{g}}}_{\mathit{\epsilon }}^1:\nabla {\mathit{\eta }}_{\mathit{\epsilon }}\mathit{d}\mathit{x}\to {\iint }_{\mathrm{\Omega }\times {\mathit{Y}}^0}\mathrm{(}\nabla {\mathit{g}}^1+{\nabla }_{\mathit{y}}\mathit{\psi }\mathrm{)}:{\nabla }_{\mathit{y}}\mathit{\phi }\mathit{d}\mathit{x}\mathit{d}\mathit{y}.\end{array}$$ Since $\nabla {\mathit{g}}^1$ is independent of y, and because $\mathit{\psi }\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{,}\mathcal{A}\mathrm{(}{\mathit{Y}}^0\mathrm{)}\mathrm{)}$ , the integral on the right-hand side vanishes. Hence, $\Vert \nabla {\mathit{\eta }}_{\mathit{\epsilon }}{\Vert }_{{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}}^2\to 0$ , and thus ${\mathit{g}}_{\mathit{\epsilon }}^1$ satisfies (49).

Step 3. Conclusion.

As can be shown by appealing to a combination of a density argument and a diagonal-sequence argument, similar to Step 1, there exists a sequence ${\mathit{g}}_{\mathit{\epsilon }}^0\in {\mathit{H}}_0^1\mathrm{(}{\mathrm{\Omega }}_{\mathit{\epsilon }}^0\mathrm{)}$ such that $$\begin{array}{c}{\mathit{g}}_{\mathit{\epsilon }}^0\stackrel{2}{\to }{\mathit{g}}^0\mathrm{,}\mathit{\epsilon }\nabla {\mathit{g}}_{\mathit{\epsilon }}^0\stackrel{2}{\to }{\nabla }_{\mathit{y}}{\mathit{g}}^0\mathrm{,}\underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}sup}{\mathit{c}}_{\mathit{\epsilon }}{\Vert {\mathit{\epsilon }}^{1-\mathit{\gamma }}\nabla {\mathit{g}}_{\mathit{\epsilon }}^0\Vert }_{{\mathit{L}}^{\infty }\mathrm{(}\mathrm{\Omega }\mathrm{)}}=0.\end{array}$$ Now define ${\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{(}\mathit{x}\mathrm{)}:={\mathit{\epsilon }}^{1-\mathit{\gamma }}{\mathit{g}}_{\mathit{\epsilon }}^0+{\mathit{g}}_{\mathit{\epsilon }}^1$ , and note that $\mathrm{(}{\mathit{g}}_{\mathit{\epsilon }}^0\mathrm{,}{\mathit{g}}_{\mathit{\epsilon }}^1\mathrm{)}$ satisfy (16). In view of the convergence of ${\mathit{g}}_{\mathit{\epsilon }}^0$ and ${\mathit{g}}_{\mathit{\epsilon }}^1$ , the sequence ${\mathit{\phi }}_{\mathit{\epsilon }}$ has the required properties. □

4.2. Proof of Theorems 1, 2 and Propositions 1, 2: Small strain regime

As a preliminary remark, we note that two different effects play a role when passing to the limit $\mathit{\epsilon }\downarrow 0$ in the small strain regime:

The non-convex energy functional is linearised at identity map (which is a stress-free state for ${\mathcal{E}}_{\mathit{\epsilon }}$ ) – this corresponds to the passage from nonlinear to linearised elasticity.

The obtained linearised, still oscillating, convex-quadratic energy is homogenised.

The following lemma is used to treat both effects simultaneously. Its proof combines convex homogenisation methods (e.g. [37, Proposition 1.3]) with a “careful Taylor expansion” in the spirit of [23, Proof of Theorem 6.2]. For notational convenience, we introduce two “linearised” functionals:

• For ${\mathit{G}}_{\mathit{\epsilon }}=\mathrm{(}{\mathit{G}}_{\mathit{\epsilon }}^0\mathrm{,}{\mathit{G}}_{\mathit{\epsilon }}^1\mathrm{)}\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{,}{\mathbb{R}}^{\mathit{d}\times \mathit{d}}\mathrm{)}\times {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{,}{\mathbb{R}}^{\mathit{d}\times \mathit{d}}\mathrm{)}$ set $$\begin{array}{c}{\mathcal{Q}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{G}}_{\mathit{\epsilon }}\mathrm{)}={\mathcal{Q}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{G}}_{\mathit{\epsilon }}^0\mathrm{,}{\mathit{G}}_{\mathit{\epsilon }}^1\mathrm{)}:={{\displaystyle \int }}_{{\mathrm{\Omega }}_{\mathit{\epsilon }}^0}{\mathit{Q}}^0\mathrm{(}{\mathit{G}}_{\mathit{\epsilon }}^0\mathrm{(}\mathit{x}\mathrm{)}\mathrm{)}\mathit{d}\mathit{x}+{{\displaystyle \int }}_{{\mathrm{\Omega }}_{\mathit{\epsilon }}^1}{\mathit{Q}}^1\mathrm{(}{\mathit{G}}_{\mathit{\epsilon }}^1\mathrm{(}\mathit{x}\mathrm{)}\mathrm{)}\mathit{d}\mathit{x}.\end{array}$$

• For $\mathit{G}=\mathrm{(}{\mathit{G}}^0\mathrm{,}{\mathit{G}}^1\mathrm{)}\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\times \mathit{Y}\mathrm{,}{\mathbb{R}}^{\mathit{d}\times \mathit{d}}\mathrm{)}\times {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\times \mathit{Y}\mathrm{,}{\mathbb{R}}^{\mathit{d}\times \mathit{d}}\mathrm{)}$ set $$\begin{array}{c}\mathcal{Q}\mathrm{(}\mathit{G}\mathrm{)}=\mathcal{Q}\mathrm{(}{\mathit{G}}^0\mathrm{,}{\mathit{G}}^1\mathrm{)}:={{\displaystyle \int }}_{\mathrm{\Omega }\times {\mathit{Y}}^0}{\mathit{Q}}^0\mathrm{(}{\mathit{G}}^0\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}\mathrm{)}\mathit{d}\mathit{x}\mathit{d}\mathit{y}+{{\displaystyle \int }}_{\mathrm{\Omega }\times {\mathit{Y}}^1}{\mathit{Q}}^1\mathrm{(}{\mathit{G}}^1\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}\mathrm{)}\mathit{d}\mathit{x}\mathit{d}\mathit{y}.\end{array}$$

Lemma 6.

Consider sequences ${\mathit{g}}_{\mathit{\epsilon }}^0\mathrm{,}{\mathit{g}}_{\mathit{\epsilon }}^1\in {\mathit{H}}^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ that satisfy $$\begin{array}{c}\underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}sup}\mathrm{(}{\Vert \mathit{\epsilon }\nabla {\mathit{g}}_{\mathit{\epsilon }}^0\Vert }_{{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}}+{\Vert \nabla {\mathit{g}}_{\mathit{\epsilon }}^1\Vert }_{{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}}\mathrm{)}\mathrm{<}\infty .\end{array}$$ Set ${\mathit{\phi }}_{\mathit{\epsilon }}:={\mathit{\epsilon }}^{1-\mathit{\gamma }}{\mathit{g}}_{\mathit{\epsilon }}^0+{\mathit{g}}_{\mathit{\epsilon }}^1$ and ${\mathit{G}}_{\mathit{\epsilon }}=\mathrm{(}{\mathit{G}}_{\mathit{\epsilon }}^0\mathrm{,}{\mathit{G}}_{\mathit{\epsilon }}^1\mathrm{)}:=\mathrm{(}\mathit{\epsilon }\nabla {\mathit{g}}_{\mathit{\epsilon }}^0\mathrm{,}\nabla {\mathit{g}}_{\mathit{\epsilon }}^1\mathrm{)}$ .

If ${\mathit{G}}_{\mathit{\epsilon }}\stackrel{2}{\mathrm{\rightharpoonup }}\mathit{G}$ and ${\mathit{m}}_{\mathit{\epsilon }}=\mathit{o}\mathrm{(}{\mathit{\epsilon }}^{\mathit{\gamma }}\mathrm{)}$ as $\mathit{\epsilon }\downarrow 0$ , then $$\begin{array}{}\text{(50)}& \begin{array}{c}\hfill \underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}inf}{\mathit{m}}_{\mathit{\epsilon }}^{-2}{\mathcal{E}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{m}}_{\mathit{\epsilon }}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}⩾\underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}inf}{\mathcal{Q}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{\theta }}_{\mathit{\epsilon }}{\mathit{G}}_{\mathit{\epsilon }}\mathrm{)}⩾\mathcal{Q}\mathrm{(}\mathit{G}\mathrm{)}\mathrm{,}\end{array}\end{array}$$ where ${\mathit{\theta }}_{\mathit{\epsilon }}:\mathrm{\Omega }\to \{0\mathrm{,}1\}$ is defined by $$\begin{array}{c}{\mathit{\theta }}_{\mathit{\epsilon }}\mathrm{(}\mathit{x}\mathrm{)}:=\left\{\begin{array}{cc}1\mathrm{,}\hfill & \mathit{\text{if }}|\nabla {\mathit{\phi }}_{\mathit{\epsilon }}|⩽{\mathrm{(}{\mathit{m}}_{\mathit{\epsilon }}{\mathit{\epsilon }}^{\mathit{\gamma }}\mathrm{)}}^{-1\mathrm{/}2}\mathrm{,}\hfill \\ 0\mathrm{,}\hfill & \mathit{\text{otherwise.}}\hfill \end{array}\right.\end{array}$$

If ${\mathit{G}}_{\mathit{\epsilon }}\stackrel{2}{\to }\mathit{G}$ , ${\mathit{m}}_{\mathit{\epsilon }}=\mathit{o}\mathrm{(}{\mathit{\epsilon }}^{\mathit{\gamma }}\mathrm{)}$ as $\mathit{\epsilon }\downarrow 0$ , and $$\begin{array}{}\text{(51)}& \underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}sup}\Vert {\mathit{m}}_{\mathit{\epsilon }}\nabla {\mathit{\phi }}_{\mathit{\epsilon }}{\Vert }_{{\mathit{L}}^{\infty }\mathrm{(}\mathrm{\Omega }\mathrm{)}}=0\mathrm{,}\end{array}$$ then $$\begin{array}{c}\begin{array}{c}\hfill \underset{\mathit{\epsilon }\downarrow 0}{lim}{\mathit{m}}_{\mathit{\epsilon }}^{-2}{\mathcal{E}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{m}}_{\mathit{\epsilon }}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}=\mathcal{Q}\mathrm{(}\mathit{G}\mathrm{)}.\end{array}\end{array}$$

Remark 3.

In Lemma 6, following [23], the function ${\mathit{\theta }}_{\mathit{\epsilon }}$ is introduced, in order to truncate the peaks of  ${\mathit{G}}_{\mathit{\epsilon }}$ . This is needed for exploiting the quadratic expansion (W3). Since both $\mathit{\epsilon }\nabla {\mathit{g}}_{\mathit{\epsilon }}^0$ and $\nabla {\mathit{g}}_{\mathit{\epsilon }}^1$ are assumed to be bounded sequences in ${\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}$ , we deduce, from the definition of ${\mathit{\theta }}_{\mathit{\epsilon }}$ , the fact that ${\mathit{m}}_{\mathit{\epsilon }}{\mathit{\epsilon }}^{-\mathit{\gamma }}=\mathit{o}\mathrm{(}1\mathrm{)}$ and the Chebyshev inequality, that $$\begin{array}{}\text{(52)}& \forall \mathit{r}\mathrm{<}\infty \Vert {\mathit{\theta }}_{\mathit{\epsilon }}-1{\Vert }_{{\mathit{L}}^{\mathit{r}}\mathrm{(}\mathrm{\Omega }\mathrm{)}}\to 0\mathrm{,}\text{and}\Vert {\mathit{\theta }}_{\mathit{\epsilon }}{\Vert }_{{\mathit{L}}^{\infty }\mathrm{(}\mathrm{\Omega }\mathrm{)}}⩽1.\end{array}$$

In the proof of Lemma 6 we need to pass to the limit in products of the form ${\mathit{f}}_{\mathit{\epsilon }}{\mathit{\theta }}_{\mathit{\epsilon }}$ , where ${\mathit{\theta }}_{\mathit{\epsilon }}$ satisfies (52), or ${\mathit{f}}_{\mathit{\epsilon }}{\mathit{\chi }}_{\mathit{\epsilon }}$ , where ${\mathit{\chi }}_{\mathit{\epsilon }}$ denotes the indicator of ${\mathrm{\Omega }}_{\mathit{\epsilon }}^0$ . This is done by appealing to the next two lemmas, the proofs of which are elementary and left to the reader.

Lemma 7.

Let ${\mathit{f}}_{\mathit{\epsilon }}\mathrm{,}{\mathit{\theta }}_{\mathit{\epsilon }}$ be sequences in ${\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}$ and assume that ${\mathit{\theta }}_{\mathit{\epsilon }}$ satisfies ( 52 ), then the following implications are valid: $$\begin{array}{ll}\hfill \underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}sup}\Vert {\mathit{f}}_{\mathit{\epsilon }}{\Vert }_{{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}}\mathrm{<}\infty & \Rightarrow \Vert {\mathit{\theta }}_{\mathit{\epsilon }}{\mathit{f}}_{\mathit{\epsilon }}-{\mathit{f}}_{\mathit{\epsilon }}{\Vert }_{{\mathit{L}}^{\mathit{p}}\mathrm{(}\mathrm{\Omega }\mathrm{)}}\to 0\forall \mathit{p}\mathrm{<}2\mathrm{,}\\ \hfill {\mathit{f}}_{\mathit{\epsilon }}\mathrm{\rightharpoonup }\mathit{f}& \Rightarrow {\mathit{\theta }}_{\mathit{\epsilon }}{\mathit{f}}_{\mathit{\epsilon }}\mathrm{\rightharpoonup }\mathit{f}\mathrm{,}\\ \hfill {\mathit{f}}_{\mathit{\epsilon }}\stackrel{2}{\mathrm{\rightharpoonup }}\mathit{f}& \Rightarrow {\mathit{\theta }}_{\mathit{\epsilon }}{\mathit{f}}_{\mathit{\epsilon }}\stackrel{2}{\mathrm{\rightharpoonup }}\mathit{f}\mathrm{,}\\ \hfill \left\{\begin{array}{cc}\hfill & {\mathit{f}}_{\mathit{\epsilon }}\mathrm{\rightharpoonup }\mathit{f}\hfill \\ \hfill & \mathrm{(}|{\mathit{f}}_{\mathit{\epsilon }}{|}^2\mathrm{)}\mathit{\text{equi-integrable}}\hfill \end{array}\right\}& \Rightarrow \Vert {\mathit{\theta }}_{\mathit{\epsilon }}{\mathit{f}}_{\mathit{\epsilon }}-{\mathit{f}}_{\mathit{\epsilon }}{\Vert }_{{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}}\to 0.\end{array}$$

Lemma 8.

Suppose that ${\mathit{f}}_{\mathit{\epsilon }}$ be a sequence in ${\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}$ and, as above, let ${\mathit{\chi }}_{\mathit{\epsilon }}$ denote the set indicator function of ${\mathrm{\Omega }}_{\mathit{\epsilon }}^0$ . Then the following implications hold: $$\begin{array}{ll}\hfill {\mathit{f}}_{\mathit{\epsilon }}\stackrel{2}{\mathrm{\rightharpoonup }}\mathit{f}& \Rightarrow {\mathit{\chi }}_{\mathit{\epsilon }}{\mathit{f}}_{\mathit{\epsilon }}\stackrel{2}{\mathrm{\rightharpoonup }}\mathit{\chi }\mathrm{(}\mathit{y}\mathrm{)}\mathit{f}\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}\mathrm{,}\\ \hfill {\mathit{f}}_{\mathit{\epsilon }}\stackrel{2}{\to }\mathit{f}& \Rightarrow {\mathit{\chi }}_{\mathit{\epsilon }}{\mathit{f}}_{\mathit{\epsilon }}\stackrel{2}{\to }\mathit{\chi }\mathrm{(}\mathit{y}\mathrm{)}\mathit{f}\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}\mathrm{,}\end{array}$$ where χ denotes the indicator function of ${\mathit{Y}}^0$ .

Proof of Lemma 6.

Step 1. Linearisation.

We claim that the following statement holds for $\mathit{i}=1\mathrm{,}2$ : Let ${\mathit{F}}_{\mathit{\epsilon }}$ denote a sequence in ${\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{,}{\mathbb{R}}^{\mathit{d}\times \mathit{d}}\mathrm{)}$ , and let ${\mathit{c}}_{\mathit{\epsilon }}$ be a sequence of positive numbers converging to zero, such that $$\begin{array}{}\text{(53)}& \underset{\mathit{\epsilon }\to 0}{lim\hspace{0.17em}sup}{\mathit{c}}_{\mathit{\epsilon }}\Vert {\mathit{F}}_{\mathit{\epsilon }}{\Vert }_{{\mathit{L}}^{\infty }\mathrm{(}\mathrm{\Omega }\mathrm{)}}=0.\end{array}$$ Then the convergence $$\begin{array}{}\text{(54)}& \underset{\mathit{\epsilon }\downarrow 0}{lim}|{\mathit{c}}_{\mathit{\epsilon }}^{-2}{{\displaystyle \int }}_{\mathrm{\Omega }}{\mathit{W}}^{\mathit{i}}\mathrm{(}\mathit{I}+{\mathit{c}}_{\mathit{\epsilon }}{\mathit{F}}_{\mathit{\epsilon }}\mathrm{)}\mathit{d}\mathit{x}-{{\displaystyle \int }}_{\mathrm{\Omega }}{\mathit{Q}}^{\mathit{i}}\mathrm{(}{\mathit{F}}_{\mathit{\epsilon }}\mathrm{(}\mathit{x}\mathrm{)}\mathrm{\left)}\mathit{d}\mathit{x}\right|=0\end{array}$$ holds. Indeed, thanks to (W3) we have $$\begin{array}{c}|{\mathit{c}}_{\mathit{\epsilon }}^{-2}{\mathit{W}}^{\mathit{i}}\mathrm{(}\mathit{I}+{\mathit{c}}_{\mathit{\epsilon }}{\mathit{F}}_{\mathit{\epsilon }}\mathrm{)}-{\mathit{Q}}^{\mathit{i}}\mathrm{(}{\mathit{F}}_{\mathit{\epsilon }}\mathrm{)}|⩽|{\mathit{F}}_{\mathit{\epsilon }}{|}^2{\mathit{r}}^{\mathit{i}}\mathrm{(}{\mathit{c}}_{\mathit{\epsilon }}|{\mathit{F}}_{\mathit{\epsilon }}|\mathrm{)}⩽|{\mathit{F}}_{\mathit{\epsilon }}{|}^2{\mathit{r}}^{\mathit{i}}\mathrm{(}{\mathit{c}}_{\mathit{\epsilon }}\Vert {\mathit{F}}_{\mathit{\epsilon }}{\Vert }_{{\mathit{L}}^{\infty }\mathrm{(}\mathrm{\Omega }\mathrm{)}}\mathrm{)}\text{a.e.}\end{array}$$ Thanks to (53), and since ${\mathit{F}}_{\mathit{\epsilon }}$ is bounded in ${\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}$ , the right-hand side converges to zero in ${\mathit{L}}^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ , and (54) follows.

Step 2. Proof of part (a).

Since the energy densities ${\mathit{W}}^0$ , ${\mathit{W}}^1$ are minimised at the identity, cf. (W2), we have $$\begin{array}{}\text{(55)}& {\mathit{m}}_{\mathit{\epsilon }}^{-2}{\mathcal{E}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}⩾{\mathrm{\left(}{\mathit{m}}_{\mathit{\epsilon }}{\mathit{\epsilon }}^{-\mathit{\gamma }}\mathrm{\right)}}^{-2}{{\displaystyle \int }}_{\mathrm{\Omega }}{\mathit{W}}^0\mathrm{(}\mathit{I}+{\mathit{m}}_{\mathit{\epsilon }}{\mathit{\epsilon }}^{-\mathit{\gamma }}{\mathit{\theta }}_{\mathit{\epsilon }}{\mathit{F}}_{\mathit{\epsilon }}^0\mathrm{(}\mathit{x}\mathrm{)}\mathrm{)}\mathit{d}\mathit{x}+{\mathit{m}}_{\mathit{\epsilon }}^{-2}{{\displaystyle \int }}_{\mathrm{\Omega }}{\mathit{W}}^1\mathrm{(}\mathit{I}+{\mathit{m}}_{\mathit{\epsilon }}{\mathit{\theta }}_{\mathit{\epsilon }}{\mathit{F}}_{\mathit{\epsilon }}^1\mathrm{(}\mathit{x}\mathrm{)}\mathrm{)}\mathit{d}\mathit{x}\mathrm{,}\end{array}$$ where $$\begin{array}{}\text{(56)}& {\mathit{F}}_{\mathit{\epsilon }}^0:={\mathit{\chi }}_{\mathit{\epsilon }}\mathrm{(}{\mathit{G}}_{\mathit{\epsilon }}^0+{\mathit{\epsilon }}^{\mathit{\gamma }}{\mathit{G}}_{\mathit{\epsilon }}^1\mathrm{)}\mathrm{,}{\mathit{F}}_{\mathit{\epsilon }}^1\mathrm{(}\mathit{x}\mathrm{)}:=\mathrm{(}1-{\mathit{\chi }}_{\mathit{\epsilon }}\mathrm{)}{\mathit{G}}_{\mathit{\epsilon }}^1.\end{array}$$ Thanks to the definition of ${\mathit{\theta }}_{\mathit{\epsilon }}$ we have $\Vert {\mathit{m}}_{\mathit{\epsilon }}{\mathit{\epsilon }}^{-\mathit{\gamma }}{\mathit{\theta }}_{\mathit{\epsilon }}{\mathit{F}}_{\mathit{\epsilon }}^0{\Vert }_{{\mathit{L}}^{\infty }\mathrm{(}\mathrm{\Omega }\mathrm{)}}+\Vert {\mathit{m}}_{\mathit{\epsilon }}{\mathit{\theta }}_{\mathit{\epsilon }}{\mathit{F}}_{\mathit{\epsilon }}^1{\Vert }_{{\mathit{L}}^{\infty }\mathrm{(}\mathrm{\Omega }\mathrm{)}}\to 0$ , so that we may apply (54) to the right-hand side in (55). We get $$\begin{array}{r}\underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}inf}{\mathit{m}}_{\mathit{\epsilon }}^{-2}{\mathcal{E}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}⩾\underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}inf}{\mathcal{Q}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{\theta }}_{\mathit{\epsilon }}{\mathit{F}}_{\mathit{\epsilon }}^0\mathrm{,}{\mathit{\theta }}_{\mathit{\epsilon }}{\mathit{F}}_{\mathit{\epsilon }}^1\mathrm{)}=\underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}inf}{\mathcal{Q}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{\theta }}_{\mathit{\epsilon }}{\mathit{G}}_{\mathit{\epsilon }}^1\mathrm{,}{\mathit{\theta }}_{\mathit{\epsilon }}{\mathit{G}}_{\mathit{\epsilon }}^0\mathrm{)}\mathrm{,}\end{array}$$ where for the last identity we used the facts that ${\mathit{F}}_{\mathit{\epsilon }}^1={\mathit{G}}_{\mathit{\epsilon }}^1$ on ${\mathrm{\Omega }}_{\mathit{\epsilon }}^1$ and $$\begin{array}{c}{\Vert {\mathit{F}}_{\mathit{\epsilon }}^0-{\mathit{G}}_{\mathit{\epsilon }}^0\Vert }_{{\mathit{L}}^2\mathrm{(}{\mathrm{\Omega }}_{\mathit{\epsilon }}^0\mathrm{)}}={\Vert {\mathit{\epsilon }}^{\mathit{\gamma }}\nabla {\mathit{g}}_{\mathit{\epsilon }}^1\Vert }_{{\mathit{L}}^2\mathrm{(}{\mathrm{\Omega }}_{\mathit{\epsilon }}^0\mathrm{)}}\to 0.\end{array}$$ It remains to argue that $$\begin{array}{c}\underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}inf}{\mathcal{Q}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{\theta }}_{\mathit{\epsilon }}{\mathit{G}}_{\mathit{\epsilon }}\mathrm{)}⩾\mathcal{Q}\mathrm{(}\mathit{G}\mathrm{)}.\end{array}$$ In order to show this, notice that $$\begin{array}{}\text{(57)}& {\mathit{Q}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{\theta }}_{\mathit{\epsilon }}{\mathit{G}}_{\mathit{\epsilon }}\mathrm{)}={{\displaystyle \int }}_{\mathrm{\Omega }}{\mathit{Q}}^0\mathrm{\left(}{\mathit{\theta }}_{\mathit{\epsilon }}{\mathit{\chi }}_{\mathit{\epsilon }}{\mathit{G}}_{\mathit{\epsilon }}^0\mathrm{\right)}\mathit{d}\mathit{x}+{{\displaystyle \int }}_{\mathrm{\Omega }}{\mathit{Q}}^1\mathrm{(}{\mathit{\theta }}_{\mathit{\epsilon }}\mathrm{(}1-{\mathit{\chi }}_{\mathit{\epsilon }}\mathrm{)}{\mathit{G}}_{\mathit{\epsilon }}^1\mathrm{)}\mathit{d}\mathit{x}.\end{array}$$ From ${\mathit{G}}_{\mathit{\epsilon }}\stackrel{2}{\mathrm{\rightharpoonup }}\mathit{G}$ we deduce, using Lemma 7, Remark 3 and Lemma 8, that $$\begin{array}{c}{\mathit{\theta }}_{\mathit{\epsilon }}{\mathit{\chi }}_{\mathit{\epsilon }}{\mathit{G}}_{\mathit{\epsilon }}^0\stackrel{2}{\mathrm{\rightharpoonup }}\mathit{\chi }\mathrm{(}\mathit{y}\mathrm{)}{\mathit{G}}^0\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}\mathrm{,}{\mathit{\theta }}_{\mathit{\epsilon }}\mathrm{(}1-{\mathit{\chi }}_{\mathit{\epsilon }}\mathrm{)}{\mathit{G}}_{\mathit{\epsilon }}^1\stackrel{2}{\mathrm{\rightharpoonup }}\mathrm{(}1-\mathit{\chi }\mathrm{(}\mathit{y}\mathrm{)}\mathrm{)}{\mathit{G}}^1\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}.\end{array}$$ By appealing to the lower semicontinuity of convex integral functionals with respect to weak two-scale convergence (cf. [37, Proposition 1.3]), we deduce that the lim inf of the right-hand side in (57) is bounded below by $\mathcal{Q}\mathrm{(}\mathit{G}\mathrm{)}$ . This completes the argument.

Step 3. Proof of part (b).

We claim that $$\begin{array}{}\text{(58)}& \underset{\mathit{\epsilon }\downarrow 0}{lim}|{\displaystyle \frac{1}{{\mathit{m}}_{\mathit{\epsilon }}^2}}{\mathcal{E}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{m}}_{\mathit{\epsilon }}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}-{\mathcal{Q}}_{\mathit{\epsilon }}\mathrm{(}\mathit{\epsilon }\nabla {\mathit{g}}_{\mathit{\epsilon }}^0\mathrm{,}\nabla {\mathit{g}}_{\mathit{\epsilon }}^1\mathrm{\left)}\right|=0.\end{array}$$ Note that $$\begin{array}{c}{\mathit{m}}_{\mathit{\epsilon }}^{-2}{\mathcal{E}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}={\mathrm{\left(}{\mathit{m}}_{\mathit{\epsilon }}{\mathit{\epsilon }}^{-\mathit{\gamma }}\mathrm{\right)}}^{-2}{{\displaystyle \int }}_{\mathrm{\Omega }}{\mathit{W}}^0\mathrm{(}\mathit{I}+{\mathit{m}}_{\mathit{\epsilon }}{\mathit{\epsilon }}^{-\mathit{\gamma }}{\mathit{F}}_{\mathit{\epsilon }}^0\mathrm{(}\mathit{x}\mathrm{)}\mathrm{)}\mathit{d}\mathit{x}+{\mathit{m}}_{\mathit{\epsilon }}^{-2}{{\displaystyle \int }}_{\mathrm{\Omega }}{\mathit{W}}^1\mathrm{(}\mathit{I}+{\mathit{m}}_{\mathit{\epsilon }}{\mathit{F}}_{\mathit{\epsilon }}^1\mathrm{(}\mathit{x}\mathrm{)}\mathrm{)}\mathit{d}\mathit{x}\mathrm{,}\end{array}$$ where ${\mathit{F}}_{\mathit{\epsilon }}^0$ and ${\mathit{F}}_{\mathit{\epsilon }}^1$ are defined in (56). By (51) we have $\Vert {\mathit{m}}_{\mathit{\epsilon }}{\mathit{\epsilon }}^{-\mathit{\gamma }}{\mathit{F}}_{\mathit{\epsilon }}^0{\Vert }_{{\mathit{L}}^{\infty }\mathrm{(}\mathrm{\Omega }\mathrm{)}}+\Vert {\mathit{m}}_{\mathit{\epsilon }}{\mathit{F}}_{\mathit{\epsilon }}^1{\Vert }_{{\mathit{L}}^{\infty }\mathrm{(}\mathrm{\Omega }\mathrm{)}}\to 0$ and (54) yields $$\begin{array}{c}|{\mathit{m}}_{\mathit{\epsilon }}^{-2}{\mathcal{E}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{m}}_{\mathit{\epsilon }}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}-{\mathcal{Q}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{F}}_{\mathit{\epsilon }}^0\mathrm{,}{\mathit{F}}_{\mathit{\epsilon }}^1\mathrm{\left)}\right|\to 0\text{as }\mathit{\epsilon }\downarrow 0.\end{array}$$ Since ${\mathit{G}}_{\mathit{\epsilon }}\stackrel{2}{\to }\mathit{G}$ we have, thanks to Lemma 8: $$\begin{array}{c}{\mathit{\chi }}_{\mathit{\epsilon }}{\mathit{F}}_{\mathit{\epsilon }}^0\stackrel{2}{\to }\mathit{\chi }\mathrm{(}\mathit{y}\mathrm{)}{\mathit{G}}^0\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}\mathrm{,}\mathrm{(}1-{\mathit{\chi }}_{\mathit{\epsilon }}\mathrm{)}{\mathit{F}}_{\mathit{\epsilon }}^1\stackrel{2}{\mathrm{\rightharpoonup }}\mathrm{(}1-\mathit{\chi }\mathrm{(}\mathit{y}\mathrm{)}\mathrm{)}{\mathit{G}}^1\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}.\end{array}$$ Hence, the continuity of convex integral functionals with respect to strong two-scale convergence (cf. [36]) yields $$\begin{array}{c}{\mathcal{Q}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{F}}_{\mathit{\epsilon }}^0\mathrm{,}{\mathit{F}}_{\mathit{\epsilon }}^1\mathrm{)}\to \mathcal{Q}\mathrm{(}\mathit{G}\mathrm{)}\mathrm{,}\end{array}$$ which completes the argument. □

We are now in a position to prove the Γ-convergence statement for the energies ${\mathcal{E}}_{\mathit{\epsilon }}$ .

Proof of Theorem 1.

Step 1. Part (a) (Compactness).

Thanks to (W2) we have $$\begin{array}{c}{\mathit{m}}_{\mathit{\epsilon }}^{-2}{\mathcal{E}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{m}}_{\mathit{\epsilon }}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}⩾{\mathit{c}}_0{\mathit{m}}_{\mathit{\epsilon }}^{-2}{\mathrm{\Phi }}_{\mathit{\epsilon }}^{\mathit{\gamma }}\mathrm{(}{\mathit{m}}_{\mathit{\epsilon }}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}.\end{array}$$ Hence, the claim of Theorem 1(a) directly follows from Lemma 2.

Step 2. Part (b) (Lower bound).

Without loss of generality we assume that $$\begin{array}{c}\underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}inf}{\mathit{m}}_{\mathit{\epsilon }}^{-2}{\mathcal{E}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{m}}_{\mathit{\epsilon }}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}=\underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}sup}{\mathit{m}}_{\mathit{\epsilon }}^{-2}{\mathcal{E}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{m}}_{\mathit{\epsilon }}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}\mathrm{<}\infty .\end{array}$$ Furthermore, thanks to Lemma 2, we can assume in addition that $\nabla {\mathit{g}}_{\mathit{\epsilon }}^1\stackrel{2}{\mathrm{\rightharpoonup }}\nabla {\mathit{g}}^1+{\nabla }_{\mathit{y}}\mathit{\psi }$ for some $\mathit{\psi }\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{,}\mathcal{A}\mathrm{(}{\mathit{Y}}^0\mathrm{)}\mathrm{)}$ , so that $$\begin{array}{c}{\mathit{G}}_{\mathit{\epsilon }}:=\mathrm{(}\mathit{\epsilon }\nabla {\mathit{g}}_{\mathit{\epsilon }}^0\mathrm{,}\nabla {\mathit{g}}_{\mathit{\epsilon }}^1\mathrm{)}\stackrel{2}{\mathrm{\rightharpoonup }}\mathrm{(}{\nabla }_{\mathit{y}}{\mathit{g}}^0\mathrm{,}\nabla {\mathit{g}}^1+{\nabla }_{\mathit{y}}\mathit{\psi }\mathrm{)}=:\mathit{G}.\end{array}$$ Applying Lemma 6(a) yields $$\begin{array}{c}\underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}inf}{\mathit{m}}_{\mathit{\epsilon }}^{-2}{\mathcal{E}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{m}}_{\mathit{\epsilon }}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}⩾{{\displaystyle \int }}_{\mathrm{\Omega }\times {\mathit{Y}}^0}{\mathit{Q}}^0\mathrm{(}{\nabla }_{\mathit{y}}{\mathit{g}}^0\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}\mathrm{)}\mathit{d}\mathit{x}\mathit{d}\mathit{y}+{{\displaystyle \int }}_{\mathrm{\Omega }\times {\mathit{Y}}^1}{\mathit{Q}}^1\mathrm{(}\nabla {\mathit{g}}^1\mathrm{(}\mathit{x}\mathrm{)}+{\nabla }_{\mathit{y}}\mathit{\psi }\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}\mathrm{)}\mathit{d}\mathit{x}\mathit{d}\mathit{y}.\end{array}$$ This completes the argument, since the right-hand side is bounded from below by ${\mathcal{E}}_{\text{small}}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}$ .

Step 3. Part (c) (Upper bound).

Choose $\mathit{\psi }\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{,}\mathcal{A}\mathrm{(}{\mathit{Y}}^0\mathrm{)}\mathrm{)}$ such that $$\begin{array}{}\text{(59)}& {{\displaystyle \int }}_{\mathrm{\Omega }\times {\mathit{Y}}^1}\mathit{Q}\mathrm{(}\nabla {\mathit{g}}^1\mathrm{(}\mathit{x}\mathrm{)}+{\nabla }_{\mathit{y}}\mathit{\psi }\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}\mathrm{)}\mathit{d}\mathit{x}\mathit{d}\mathit{y}={{\displaystyle \int }}_{\mathrm{\Omega }}{\mathit{Q}}_{hom}^1\mathrm{(}\nabla {\mathit{g}}^1\mathrm{(}\mathit{x}\mathrm{)}\mathrm{)}\mathit{d}\mathit{x}.\end{array}$$ Let ${\mathit{\phi }}_{\mathit{\epsilon }}$ denote the sequence associated with ${\mathit{g}}^0\mathrm{,}{\mathit{g}}^1$ and ψ via Lemma 3 with ${\mathit{c}}_{\mathit{\epsilon }}:={\mathit{m}}_{\mathit{\epsilon }}$ . In view of (20), applying Lemma 6(b) yields $$\begin{array}{c}\underset{\mathit{\epsilon }\downarrow 0}{lim}{\mathit{m}}_{\mathit{\epsilon }}^{-2}{\mathcal{E}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{m}}_{\mathit{\epsilon }}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}={{\displaystyle \int }}_{\mathrm{\Omega }\times \mathit{Y}}{\mathit{Q}}^0\mathrm{(}{\nabla }_{\mathit{y}}{\mathit{g}}^0\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}\mathrm{)}\mathit{d}\mathit{x}\mathit{d}\mathit{y}+{{\displaystyle \int }}_{\mathrm{\Omega }\times {\mathit{Y}}^1}{\mathit{Q}}^1\mathrm{(}\nabla {\mathit{g}}^1\mathrm{(}\mathit{x}\mathrm{)}+{\nabla }_{\mathit{y}}\mathit{\psi }\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}\mathrm{)}\mathit{d}\mathit{x}\mathit{d}\mathit{y}.\end{array}$$ It follows from (59) that the right-hand side equals ${\mathcal{E}}_{\text{small}}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}$ . □

Proof of Proposition 1.

Step 1. A priori estimate.

We claim that for every sequence ${\mathit{\phi }}_{\mathit{\epsilon }}$ in ${\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ the following implication holds: $$\begin{array}{}\text{(60)}& \underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}sup}{\mathit{I}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}\mathrm{<}\infty \Rightarrow \underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}sup}{\mathit{m}}_{\mathit{\epsilon }}^{-2}{\mathcal{E}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{m}}_{\mathit{\epsilon }}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}\mathrm{<}\infty .\end{array}$$ Indeed, we have $$\begin{array}{rcl}|{{\displaystyle \int }}_{\mathrm{\Omega }}{\ell }_{\mathit{\epsilon }}·{\mathit{m}}_{\mathit{\epsilon }}{\mathit{\phi }}_{\mathit{\epsilon }}\mathit{d}\mathit{x}|& ⩽& {\mathit{m}}_{\mathit{\epsilon }}\mathrm{(}\Vert {\ell }_{\mathit{\epsilon }}{\Vert }_{{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}}{\Vert {\mathit{g}}_{\mathit{\epsilon }}^1\Vert }_{{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}}+{\mathit{\epsilon }}^{1-\mathit{\gamma }}\Vert {\ell }_{\mathit{\epsilon }}{\Vert }_{{\mathit{L}}^2\mathrm{(}{\mathrm{\Omega }}_{\mathit{\epsilon }}^0\mathrm{)}}{\Vert {\mathit{g}}_{\mathit{\epsilon }}^0\Vert }_{{\mathit{L}}^2\mathrm{(}{\mathrm{\Omega }}_{\mathit{\epsilon }}^0\mathrm{)}}\mathrm{)}\\ & \stackrel{\text{(25a)}}{⩽}& {\mathit{m}}_{\mathit{\epsilon }}\mathrm{(}{\Vert {\mathit{m}}_{\mathit{\epsilon }}{\mathit{g}}_{\mathit{\epsilon }}^1\Vert }_{{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}}+{\Vert {\mathit{m}}_{\mathit{\epsilon }}{\mathit{g}}_{\mathit{\epsilon }}^0\Vert }_{{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}}\mathrm{)}\stackrel{\text{(17)}}{\lesssim }{\mathit{m}}_{\mathit{\epsilon }}^2\sqrt{{\mathit{m}}_{\mathit{\epsilon }}^{-2}{\mathit{\varphi }}_{\mathit{\epsilon }}^{\mathit{\gamma }}\mathrm{(}{\mathit{m}}_{\mathit{\epsilon }}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}}\\ & \stackrel{\text{(W2)}}{\lesssim }& {\mathit{m}}_{\mathit{\epsilon }}^2\sqrt{{\mathit{m}}_{\mathit{\epsilon }}^{-2}{\mathcal{E}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{m}}_{\mathit{\epsilon }}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}}.\end{array}$$ Combining this with the definition of ${\mathit{I}}_{\mathit{\epsilon }}$ we get $$\begin{array}{c}{\mathit{m}}_{\mathit{\epsilon }}^{-2}{\mathcal{E}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{m}}_{\mathit{\epsilon }}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}\lesssim {\mathit{I}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}+\sqrt{{\mathit{m}}_{\mathit{\epsilon }}^{-2}{\mathcal{E}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}}\mathrm{,}\end{array}$$ which implies (60).

Step 2. The proof of parts (a) and (b).

The existence of a minimiser to ${\mathit{I}}_{\text{small}}$ follows by the direct method. The minimiser $\mathrm{(}{\mathit{g}}_{\ast }^0\mathrm{,}{\mathit{g}}_{\ast }^1\mathrm{)}$ is unique, since the implication $$\begin{array}{c}{{\displaystyle \int }}_{\mathrm{\Omega }\times \mathit{Y}}{\mathit{Q}}^0\mathrm{(}{\nabla }_{\mathit{y}}{\tilde{\mathit{g}}}^0\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}\mathrm{)}\mathit{d}\mathit{x}\mathit{d}\mathit{y}+{{\displaystyle \int }}_{\mathrm{\Omega }}{\mathit{Q}}_{hom}^1\mathrm{(}\nabla {\tilde{\mathit{g}}}^1\mathrm{(}\mathit{x}\mathrm{)}\mathrm{)}\mathit{d}\mathit{x}=0\Rightarrow {\tilde{\mathit{g}}}^0={\tilde{\mathit{g}}}^1=0\end{array}$$ holds for all ${\tilde{\mathit{g}}}^0\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{,}{\mathit{H}}_0^1\mathrm{(}{\mathit{Y}}^0\mathrm{)}\mathrm{)}$ and ${\tilde{\mathit{g}}}^1\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ .

The remaining claims of Proposition 1 follow from the standard Γ-convergence arguments (cf. [20, Corollary 7.20]), provided the functionals ${\mathit{I}}_{\mathit{\epsilon }}$ , $\mathit{\epsilon }\mathrm{>}0$ , are equi-coercive and Γ-converge to ${\mathit{I}}_{\text{small}}$ . Indeed, thanks to (25c), it is easy to check that ${\mathit{\phi }}_{\mathit{\epsilon }}\stackrel{2}{\mathrm{\rightharpoonup }}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}$ implies $$\begin{array}{}\text{(61)}& \underset{\mathit{\epsilon }\downarrow 0}{lim}{\displaystyle \frac{1}{{\mathit{m}}_{\mathit{\epsilon }}^2}}{{\displaystyle \int }}_{\mathrm{\Omega }}{\ell }_{\mathit{\epsilon }}·{\mathit{m}}_{\mathit{\epsilon }}{\mathit{\phi }}_{\mathit{\epsilon }}\mathit{d}\mathit{x}={{\displaystyle \int }}_{\mathrm{\Omega }}\mathrm{(}{{\displaystyle \int }}_{{\mathit{Y}}^0}{\ell }^0·{\mathit{g}}^0\mathit{d}\mathit{y}+{\ell }^1·{\mathit{g}}^1\mathrm{)}\mathit{d}\mathit{x}\mathrm{,}\end{array}$$ since the integral on the left-hand side only involves products of weakly and strongly two-scale convergent factors, cf. [36, Proposition 2.8]). In combination with Theorem 1, this implies that ${\mathit{I}}_{\mathit{\epsilon }}$ Γ-converges to ${\mathit{I}}_{\text{small}}$ . In addition, the trivial inequality $$\begin{array}{c}\underset{\mathit{\phi }\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}}{inf}{\mathit{I}}_{\mathit{\epsilon }}\mathrm{(}\mathit{\phi }\mathrm{)}⩽{\mathit{I}}_{\mathit{\epsilon }}\mathrm{(}0\mathrm{)}=0\mathrm{,}\end{array}$$ combined with (60) and Lemma 2, proves that the functionals ${\mathit{I}}_{\mathit{\epsilon }}$ are equi-coercive. □

For the proof of Theorem 2 we make use of the following lemma: Lemma 9 (Decomposition Lemma, see [22,27]).

Let $\mathit{v}\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ and let ${\mathit{v}}_{\mathit{\epsilon }}\in {\mathit{H}}^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ be a sequence with ${\mathit{v}}_{\mathit{\epsilon }}\mathrm{\rightharpoonup }\mathit{v}$ in ${\mathit{H}}^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ . Then there exists a sequence ${\bar{\mathit{V}}}_{\mathit{\epsilon }}\in {\mathit{H}}^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ such that the following properties hold for a subsequence of ${\mathit{v}}_{\mathit{\epsilon }}$ (not relabelled):

${\bar{\mathit{V}}}_{\mathit{\epsilon }}\mathrm{\rightharpoonup }\mathit{v}$ in ${\mathit{H}}^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ ;

${\bar{\mathit{V}}}_{\mathit{\epsilon }}={\mathit{v}}_{\mathit{\epsilon }}$ in a neighbourhood of $\partial \mathrm{\Omega }$ ;

$\mathrm{(}|\nabla {\bar{\mathit{V}}}_{\mathit{\epsilon }}{|}^2\mathrm{)}$ is equi-integrable;

$|\{\mathit{x}\in \mathrm{\Omega }:{\mathit{v}}_{\mathit{\epsilon }}\mathrm{(}\mathit{x}\mathrm{)}\ne {\bar{\mathit{V}}}_{\mathit{\epsilon }}\mathrm{(}\mathit{x}\mathrm{)}\}|\to 0$ as $\mathit{\epsilon }\downarrow 0$ .

Proof of Theorem 2.

It suffices to prove the theorem for a subsequence. Throughout the proof we write $$\begin{array}{c}{\mathit{G}}_{\mathit{\epsilon }}:=\mathrm{(}\mathit{\epsilon }\nabla {\mathit{g}}_{\mathit{\epsilon }}^0\mathrm{,}\nabla {\mathit{g}}_{\mathit{\epsilon }}^1\mathrm{)}\mathrm{,}\mathit{G}:=\mathrm{(}{\nabla }_{\mathit{y}}{\mathit{g}}_{\ast }^0\mathrm{,}\nabla {\mathit{g}}_{\ast }^1+{\nabla }_{\mathit{y}}{\mathit{\psi }}_{\ast }\mathrm{)}.\end{array}$$ Furthermore, we make use of the functionals ${\mathcal{Q}}_{\mathit{\epsilon }}$ and $\mathcal{Q}$ introduced at the beginning of Section 4.2. Recall that $$\begin{array}{c}{\mathcal{E}}_{\text{small}}\mathrm{(}{\mathit{g}}_{\ast }^0\mathrm{,}{\mathit{g}}_{\ast }^1\mathrm{)}=\mathit{Q}\mathrm{(}\mathit{G}\mathrm{)}.\end{array}$$

Step 1. Convergence of ${\mathit{\phi }}_{\mathit{\epsilon }}$ and of the corresponding energy values.

We claim that, as $\mathit{\epsilon }\downarrow 0$ , one has $$\begin{array}{}\text{(62)}& {\mathit{G}}_{\mathit{\epsilon }}& \stackrel{2}{\mathrm{\rightharpoonup }}& \mathit{G}\mathrm{,}\text{(63)}& {\mathit{I}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}& \to & {\mathit{I}}_{\text{small}}\mathrm{(}{\mathit{g}}_{\ast }^0\mathrm{,}{\mathit{g}}_{\ast }^1\mathrm{)}\mathrm{,}\text{(64)}& {\mathit{m}}_{\mathit{\epsilon }}^{-2}{\mathcal{E}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{m}}_{\mathit{\epsilon }}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}& \to & {\mathcal{E}}_{\text{small}}\mathrm{(}{\mathit{g}}_{\ast }^0\mathrm{,}{\mathit{g}}_{\ast }^1\mathrm{)}.\end{array}$$ Indeed, from Proposition 1 we immediately deduce that ${\mathit{\phi }}_{\mathit{\epsilon }}\stackrel{2}{\mathrm{\rightharpoonup }}\mathrm{(}{\mathit{g}}_{\ast }^0\mathrm{,}{\mathit{g}}_{\ast }^1\mathrm{)}$ and (63). Furthermore, in view of the continuity of the loading term, cf. (61), this implies (64). For (62), it remains to argue that $\nabla {\mathit{g}}_{\mathit{\epsilon }}^1\stackrel{2}{\mathrm{\rightharpoonup }}\nabla {\mathit{g}}_{\ast }^1+{\nabla }_{\mathit{y}}{\mathit{\psi }}_{\ast }$ . Thanks to ${\mathit{\phi }}_{\mathit{\epsilon }}\stackrel{2}{\mathrm{\rightharpoonup }}\mathrm{(}{\mathit{g}}_{\ast }^0\mathrm{,}{\mathit{g}}_{\ast }^1\mathrm{)}$ and Lemma 2 we have, up to a subsequence, $\nabla {\mathit{g}}_{\mathit{\epsilon }}^1\stackrel{2}{\mathrm{\rightharpoonup }}\nabla {\mathit{g}}_{\ast }^1+{\nabla }_{\mathit{y}}\mathit{\psi }$ for some $\mathit{\psi }\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{,}\mathcal{A}\mathrm{(}{\mathit{Y}}^0\mathrm{)}\mathrm{)}$ . Furthermore, from (64) and Lemma 6(a) we infer that $$\begin{array}{c}{\mathit{E}}_{\text{small}}\mathrm{(}{\mathit{g}}_0^{\ast }\mathrm{,}{\mathit{g}}_1^{\ast }\mathrm{)}=\underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}inf}{\mathit{m}}_{\mathit{\epsilon }}^{-2}{\mathcal{E}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{m}}_{\mathit{\epsilon }}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}⩾{{\displaystyle \int }}_{\mathrm{\Omega }\times {\mathit{Y}}^0}{\mathit{Q}}^0\mathrm{\left(}{\nabla }_{\mathit{y}}{\mathit{g}}_{\ast }^0\mathrm{\right)}\mathit{d}\mathit{x}\mathit{d}\mathit{y}+{{\displaystyle \int }}_{\mathrm{\Omega }\times {\mathit{Y}}^1}{\mathit{Q}}^1\mathrm{(}\nabla {\mathit{g}}_{\ast }^1+{\nabla }_{\mathit{y}}\mathit{\psi }\mathrm{)}\mathit{d}\mathit{x}\mathit{d}\mathit{y}.\end{array}$$ This, in particular, implies $$\begin{array}{c}{{\displaystyle \int }}_{\mathrm{\Omega }\times {\mathit{Y}}^1}{\mathit{Q}}^1\mathrm{(}\nabla {\mathit{g}}_{\ast }^1+{\nabla }_{\mathit{y}}\mathit{\psi }\mathrm{)}\mathit{d}\mathit{x}\mathit{d}\mathit{y}={{\displaystyle \int }}_{\mathrm{\Omega }}{\mathit{Q}}_{hom}^1\mathrm{(}\nabla {\mathit{g}}_{\ast }^1\mathrm{)}\mathit{d}\mathit{x}.\end{array}$$ In view of (28) we conclude that $\mathit{\psi }={\mathit{\psi }}_{\ast }$ and (62) follows.

Step 2. Equi-integrable decomposition.

We claim that for a subsequence (not relabelled) there exist sequences ${\overline{\mathit{g}}}_{\mathit{\epsilon }}^0\mathrm{,}{\bar{\mathit{G}}}_{\mathit{\epsilon }}^1\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ such that ${\bar{\mathit{G}}}_{\mathit{\epsilon }}:=\mathrm{(}\mathit{\epsilon }\nabla {\bar{\mathit{G}}}_{\mathit{\epsilon }}^0\mathrm{,}\nabla {\bar{\mathit{G}}}_{\mathit{\epsilon }}^1\mathrm{)}$ satisfies $$\begin{array}{}\text{(65)}& {\bar{\mathit{G}}}_{\mathit{\epsilon }}-{\mathit{G}}_{\mathit{\epsilon }}\stackrel{2}{\mathrm{\rightharpoonup }}0\mathrm{,}\Vert {\bar{\mathit{G}}}_{\mathit{\epsilon }}-{\mathit{G}}_{\mathit{\epsilon }}{\Vert }_{{\mathit{L}}^{\mathit{p}}\mathrm{(}\mathrm{\Omega }\mathrm{)}}\to 0\mathrm{,}\end{array}$$ and $$\begin{array}{}\text{(66)}& {\mathcal{Q}}_{\mathit{\epsilon }}\mathrm{(}{\bar{\mathit{G}}}_{\mathit{\epsilon }}\mathrm{)}\to {\mathcal{E}}_{\text{small}}\mathrm{(}{\mathit{g}}_{\ast }^0\mathrm{,}{\mathit{g}}_{\ast }^1\mathrm{)}.\end{array}$$ To show the above, notice that thanks to Lemma 9 there exist sequences ${\bar{\mathit{G}}}_{\mathit{\epsilon }}^0\mathrm{,}{\bar{\mathit{G}}}_{\mathit{\epsilon }}^1\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ such that

$\mathrm{(}{\mathit{\epsilon }}^2|\nabla {\bar{\mathit{G}}}_{\mathit{\epsilon }}^0{|}^2\mathrm{)}$ and $\mathrm{(}|\nabla {\overline{\mathit{g}}}_{\mathit{\epsilon }}^1{|}^2\mathrm{)}$ are equi-integrable,

the indicator function ${\overline{\mathit{\theta }}}_{\mathit{\epsilon }}$ defined by $$\begin{array}{}\text{(67)}& {\overline{\mathit{\theta }}}_{\mathit{\epsilon }}\mathrm{(}\mathit{x}\mathrm{)}:=\left\{\begin{array}{cc}1\mathrm{,}\hfill & \text{if }{\bar{\mathit{G}}}_{\mathit{\epsilon }}^0\mathrm{(}\mathit{x}\mathrm{)}={\mathit{g}}_{\mathit{\epsilon }}^0\mathrm{(}\mathit{x}\mathrm{)}\text{ and }{\bar{\mathit{G}}}_{\mathit{\epsilon }}^1\mathrm{(}\mathit{x}\mathrm{)}={\mathit{g}}_{\mathit{\epsilon }}^1\mathrm{(}\mathit{x}\mathrm{)}\mathrm{,}\hfill \\ 0\mathrm{,}\hfill & \text{otherwise}\mathrm{,}\hfill \end{array}\right.\end{array}$$ satisfies (52).

Since ${\bar{\mathit{G}}}_{\mathit{\epsilon }}-{\mathit{G}}_{\mathit{\epsilon }}=\mathrm{(}1-{\overline{\mathit{\theta }}}_{\mathit{\epsilon }}\mathrm{)}\mathrm{(}{\bar{\mathit{G}}}_{\mathit{\epsilon }}-{\mathit{G}}_{\mathit{\epsilon }}\mathrm{)}$ (and $\mathit{p}\mathrm{<}2$ ), the convergence (65) follows from the boundedness of the sequence $\mathrm{(}{\bar{\mathit{G}}}_{\mathit{\epsilon }}-{\mathit{G}}_{\mathit{\epsilon }}\mathrm{)}$ in ${\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}$ , Lemma 7, and Hölder’s inequality.

We prove (66). Thanks to (62) we have ${\bar{\mathit{G}}}_{\mathit{\epsilon }}\stackrel{2}{\mathrm{\rightharpoonup }}\mathit{G}$ , so that (due to the lower semicontinuity of convex integral functionals with respect to weak two-scale convergence, cf. [37, Proposition 1.3]): $$\begin{array}{c}\underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}inf}{\mathcal{Q}}_{\mathit{\epsilon }}\mathrm{(}{\bar{\mathit{G}}}_{\mathit{\epsilon }}\mathrm{)}⩾\mathcal{Q}\mathrm{(}\mathit{G}\mathrm{)}={\mathcal{E}}_{\text{small}}\mathrm{(}{\mathit{g}}_{\ast }^0\mathrm{,}{\mathit{g}}_{\ast }^1\mathrm{)}.\end{array}$$ Hence, for (66) it suffices to prove the opposite estimate, i.e.  ${lim\hspace{0.17em}sup}_{\mathit{\epsilon }\downarrow 0}{\mathcal{Q}}_{\mathit{\epsilon }}\mathrm{(}{\bar{\mathit{G}}}_{\mathit{\epsilon }}\mathrm{)}⩽{\mathcal{E}}_{\text{small}}\mathrm{(}\mathit{G}\mathrm{)}$ , which, thanks to (50) and (64), follows from $$\begin{array}{}\text{(68)}& \underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}sup}\mathrm{(}{\mathcal{Q}}_{\mathit{\epsilon }}\mathrm{(}{\bar{\mathit{G}}}_{\mathit{\epsilon }}\mathrm{)}-{\mathcal{Q}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{\theta }}_{\mathit{\epsilon }}{\mathit{G}}_{\mathit{\epsilon }}\mathrm{)}\mathrm{)}⩽0.\end{array}$$ In order to show (68) notice that since the supports of ${\overline{\mathit{\theta }}}_{\mathit{\epsilon }}$ and $\mathrm{(}1-{\overline{\mathit{\theta }}}_{\mathit{\epsilon }}\mathrm{)}$ are disjoint, and because ${\overline{\mathit{\theta }}}_{\mathit{\epsilon }}{\mathit{\theta }}_{\mathit{\epsilon }}{\mathit{G}}_{\mathit{\epsilon }}={\overline{\mathit{\theta }}}_{\mathit{\epsilon }}{\mathit{\theta }}_{\mathit{\epsilon }}{\bar{\mathit{G}}}_{\mathit{\epsilon }}$ (cf. (67)), an expansion of the squares yields $$\begin{array}{rcl}{\mathcal{Q}}_{\mathit{\epsilon }}\mathrm{(}{\bar{\mathit{G}}}_{\mathit{\epsilon }}\mathrm{)}-{\mathcal{Q}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{\theta }}_{\mathit{\epsilon }}{\mathit{G}}_{\mathit{\epsilon }}\mathrm{)}& =& {\mathcal{Q}}_{\mathit{\epsilon }}\mathrm{(}{\overline{\mathit{\theta }}}_{\mathit{\epsilon }}{\bar{\mathit{G}}}_{\mathit{\epsilon }}\mathrm{)}+{\mathcal{Q}}_{\mathit{\epsilon }}\mathrm{(}\mathrm{(}1-{\overline{\mathit{\theta }}}_{\mathit{\epsilon }}\mathrm{)}{\bar{\mathit{G}}}_{\mathit{\epsilon }}\mathrm{)}-{\mathcal{Q}}_{\mathit{\epsilon }}\mathrm{(}{\overline{\mathit{\theta }}}_{\mathit{\epsilon }}{\mathit{\theta }}_{\mathit{\epsilon }}{\mathit{G}}_{\mathit{\epsilon }}\mathrm{)}-{\mathcal{Q}}_{\mathit{\epsilon }}\mathrm{(}\mathrm{(}1-{\overline{\mathit{\theta }}}_{\mathit{\epsilon }}\mathrm{)}{\mathit{\theta }}_{\mathit{\epsilon }}{\mathit{G}}_{\mathit{\epsilon }}\mathrm{)}\\ & =& {\mathcal{Q}}_{\mathit{\epsilon }}\mathrm{(}{\overline{\mathit{\theta }}}_{\mathit{\epsilon }}\mathrm{(}1-{\mathit{\theta }}_{\mathit{\epsilon }}\mathrm{)}{\bar{\mathit{G}}}_{\mathit{\epsilon }}\mathrm{)}+{\mathcal{Q}}_{\mathit{\epsilon }}\mathrm{(}\mathrm{(}1-{\overline{\mathit{\theta }}}_{\mathit{\epsilon }}\mathrm{)}{\bar{\mathit{G}}}_{\mathit{\epsilon }}\mathrm{)}-{\mathcal{Q}}_{\mathit{\epsilon }}\mathrm{(}\mathrm{(}1-{\overline{\mathit{\theta }}}_{\mathit{\epsilon }}\mathrm{)}{\mathit{\theta }}_{\mathit{\epsilon }}{\mathit{G}}_{\mathit{\epsilon }}\mathrm{)}\\ & ⩽& {\mathcal{Q}}_{\mathit{\epsilon }}\mathrm{(}{\overline{\mathit{\theta }}}_{\mathit{\epsilon }}\mathrm{(}1-{\mathit{\theta }}_{\mathit{\epsilon }}\mathrm{)}{\bar{\mathit{G}}}_{\mathit{\epsilon }}\mathrm{)}+{\mathcal{Q}}_{\mathit{\epsilon }}\mathrm{(}\mathrm{(}1-{\overline{\mathit{\theta }}}_{\mathit{\epsilon }}\mathrm{)}{\bar{\mathit{G}}}_{\mathit{\epsilon }}\mathrm{)}.\end{array}$$ It is easy to check that ${\overline{\mathit{\theta }}}_{\mathit{\epsilon }}\mathrm{(}1-{\mathit{\theta }}_{\mathit{\epsilon }}\mathrm{)}$ and $1-{\overline{\mathit{\theta }}}_{\mathit{\epsilon }}$ converge to zero in ${\mathit{L}}^{\mathit{r}}\mathrm{(}\mathrm{\Omega }\mathrm{)}$ for all $\mathit{r}\mathrm{<}\infty$ . Hence, since $|{\bar{\mathit{G}}}_{\mathit{\epsilon }}{|}^2$ is equi-integrable, Lemma 7 implies that the right-hand side of the previous estimate converges to zero, and (68) follows.

Step 3. Error estimate.

We claim that $$\begin{array}{}\text{(69)}& {{\displaystyle \int }}_{{\mathrm{\Omega }}_{\mathit{\epsilon }}^0}\mathit{\epsilon }{|sym\nabla \mathrm{(}{\overline{\mathit{g}}}_{\mathit{\epsilon }}^0-{\mathit{g}}_{\ast \mathrm{,}\mathit{\epsilon }}^0\mathrm{\left)}\right|}^2\mathit{d}\mathit{x}+{{\displaystyle \int }}_{{\mathrm{\Omega }}_{\mathit{\epsilon }}^1}{|sym\nabla \mathrm{(}{\overline{\mathit{g}}}_{\mathit{\epsilon }}^1-{\mathit{g}}_{\ast \mathrm{,}\mathit{\epsilon }}^1\mathrm{\left)}\right|}^2\mathit{d}\mathit{x}\to 0\text{as }\mathit{\epsilon }\downarrow 0.\end{array}$$ For the argument set ${\mathit{G}}_{\ast \mathrm{,}\mathit{\epsilon }}:=\mathrm{(}\mathit{\epsilon }\nabla {\mathit{g}}_{\ast \mathrm{,}\mathit{\epsilon }}^0\mathrm{,}\nabla {\mathit{g}}_{\ast \mathrm{,}\mathit{\epsilon }}^1\mathrm{)}$ . In view of (13) it suffices to argue that $$\begin{array}{c}{\mathcal{Q}}_{\mathit{\epsilon }}\mathrm{(}{\bar{\mathit{G}}}_{\mathit{\epsilon }}-{\mathit{G}}_{\ast \mathrm{,}\mathit{\epsilon }}\mathrm{)}\to 0\text{as }\mathit{\epsilon }\downarrow 0.\end{array}$$ The latter can be seen as follows: We have $$\begin{array}{r}{\mathcal{Q}}_{\mathit{\epsilon }}\mathrm{(}{\bar{\mathit{G}}}_{\mathit{\epsilon }}-{\mathit{G}}_{\ast \mathrm{,}\mathit{\epsilon }}\mathrm{)}={\mathcal{Q}}_{\mathit{\epsilon }}\mathrm{(}{\bar{\mathit{G}}}_{\mathit{\epsilon }}\mathrm{)}-{\mathcal{Q}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{G}}_{\ast \mathrm{,}\mathit{\epsilon }}\mathrm{)}+2{\mathcal{B}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{G}}_{\ast \mathrm{,}\mathit{\epsilon }};{\mathit{G}}_{\ast \mathrm{,}\mathit{\epsilon }}-{\bar{\mathit{G}}}_{\mathit{\epsilon }}\mathrm{)}\mathrm{,}\end{array}$$ where ${\mathcal{B}}_{\mathit{\epsilon }}$ denotes the bilinear form associated with ${\mathcal{Q}}_{\mathit{\epsilon }}$ .

The difference of the two quadratic terms on the right-hand side converges to zero, since ${\mathit{G}}_{\ast \mathrm{,}\mathit{\epsilon }}$ is associated with a recovery sequence, and thanks to (66). On the other hand, since ${\mathit{G}}_{\ast \mathrm{,}\mathit{\epsilon }}$ strongly two-scale converges, and ${\mathit{G}}_{\ast \mathrm{,}\mathit{\epsilon }}-{\overline{\mathit{G}}}_{\mathit{\epsilon }}\stackrel{2}{\mathrm{\rightharpoonup }}0$ by (65), we deduce that $$\begin{array}{c}\underset{\mathit{\epsilon }\downarrow 0}{lim}{\mathcal{B}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{G}}_{\ast \mathrm{,}\mathit{\epsilon }};{\mathit{G}}_{\ast \mathrm{,}\mathit{\epsilon }}-{\bar{\mathit{G}}}_{\mathit{\epsilon }}\mathrm{)}=0\mathrm{,}\end{array}$$ as ${\mathcal{B}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{G}}_{\ast \mathrm{,}\mathit{\epsilon }};{\mathit{G}}_{\ast \mathrm{,}\mathit{\epsilon }}-{\bar{\mathit{G}}}_{\mathit{\epsilon }}\mathrm{)}$ only involves products between a weakly and a strongly two-scale convergent factor (cf. [36, Proposition 2.8]).

Step 4. Conclusion (Proof of (31)).

We split the estimate into $$\begin{array}{}\text{(70)}& {{\displaystyle \int }}_{\mathrm{\Omega }}\mathrm{(}{|{\mathit{g}}_{\mathit{\epsilon }}^0-{\mathit{g}}_{\ast \mathrm{,}\mathit{\epsilon }}^0|}^{\mathit{p}}+{|\mathit{\epsilon }\nabla {\mathit{g}}_{\mathit{\epsilon }}^0-\mathit{\epsilon }\nabla {\mathit{g}}_{\ast \mathrm{,}\mathit{\epsilon }}^0|}^{\mathit{p}}\mathrm{)}\mathit{d}\mathit{x}& \to & 0\text{as }\mathit{\epsilon }\downarrow 0\mathrm{,}\text{(71)}& {{\displaystyle \int }}_{\mathrm{\Omega }}\mathrm{(}{|{\mathit{g}}_{\mathit{\epsilon }}^1-{\mathit{g}}_{\ast \mathrm{,}\mathit{\epsilon }}^1|}^{\mathit{p}}+{|\nabla {\mathit{g}}_{\mathit{\epsilon }}^1-\nabla {\mathit{g}}_{\ast \mathrm{,}\mathit{\epsilon }}^1|}^{\mathit{p}}\mathrm{)}\mathit{d}\mathit{x}& \to & 0\text{as }\mathit{\epsilon }\downarrow 0.\end{array}$$ Thanks to (65) and Step 3 we have $$\begin{array}{}\text{(72)}& {{\displaystyle \int }}_{{\mathrm{\Omega }}_{\mathit{\epsilon }}^0}\mathit{\epsilon }{|sym\nabla \mathrm{(}{\mathit{g}}_{\mathit{\epsilon }}^0-{\mathit{g}}_{\ast \mathrm{,}\mathit{\epsilon }}^0\mathrm{\left)}\right|}^{\mathit{p}}\mathit{d}\mathit{x}+{{\displaystyle \int }}_{{\mathrm{\Omega }}_{\mathit{\epsilon }}^1}{|sym\nabla \mathrm{(}{\mathit{g}}_{\mathit{\epsilon }}^1-{\mathit{g}}_{\ast \mathrm{,}\mathit{\epsilon }}^1\mathrm{\left)}\right|}^{\mathit{p}}\mathit{d}\mathit{x}\to 0\text{as }\mathit{\epsilon }\downarrow 0.\end{array}$$ Argument for (70): Set ${\mathit{\eta }}_{\mathit{\epsilon }}^0:={\mathit{g}}_{\mathit{\epsilon }}^0-{\mathit{g}}_{\ast \mathrm{,}\mathit{\epsilon }}^0$ . Since ${\mathit{\eta }}_{\mathit{\epsilon }}^0\in {\mathit{H}}_0^1\mathrm{(}{\mathrm{\Omega }}_{\mathit{\epsilon }}^0\mathrm{)}\subset {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ , Korn’s inequality yields $$\begin{array}{c}{{\displaystyle \int }}_{\mathrm{\Omega }}{|\nabla {\mathit{\eta }}_{\mathit{\epsilon }}^0|}^{\mathit{p}}\mathit{d}\mathit{x}⩽\mathit{C}{{\displaystyle \int }}_{\mathrm{\Omega }}{|sym\nabla {\mathit{\eta }}_{\mathit{\epsilon }}^0|}^{\mathit{p}}\mathit{d}\mathit{x}=\mathit{C}{{\displaystyle \int }}_{{\mathrm{\Omega }}_{\mathit{\epsilon }}^0}{|sym\nabla {\mathit{\eta }}_{\mathit{\epsilon }}^0|}^{\mathit{p}}\mathit{d}\mathit{x}\mathrm{,}\end{array}$$ where $\mathit{C}\mathrm{>}0$ only depends on Ω, p, and d. Combined with the improved Poincaré inequality (43) and (72), (70) follows.

Argument for (71): We claim that (71) follows from $$\begin{array}{}\text{(73)}& {\Vert sym\nabla {\mathit{\eta }}_{\mathit{\epsilon }}^1\Vert }_{{\mathit{L}}^{\mathit{p}}\mathrm{(}{\mathrm{\Omega }}_{\mathit{\epsilon }}^0\mathrm{)}}\to 0\mathrm{,}\end{array}$$ where ${\mathit{\eta }}_{\mathit{\epsilon }}^1:={\mathit{g}}_{\mathit{\epsilon }}^1-{\mathit{g}}_{\ast \mathrm{,}\mathit{\epsilon }}^1$ . Indeed, since ${\mathit{\eta }}_{\mathit{\epsilon }}^1$ vanishes on $\partial \mathrm{\Omega }$ , (73), (72) and Korn’s first inequality yield (71).

Thanks to the definition of ${\mathrm{\Omega }}_{\mathit{\epsilon }}^0$ , the argument for (73) can be reduced to the following statement: For all $\mathit{\xi }\in {\mathit{Z}}_{\mathit{\epsilon }}:=\{\mathit{\xi }\in {\mathbb{Z}}^{\mathit{d}}:\mathit{\epsilon }\mathrm{(}\mathit{\xi }+\mathit{Y}\mathrm{)}\subset \mathrm{\Omega }\}$ we have $$\begin{array}{}\text{(74)}& {{\displaystyle \int }}_{\mathit{\epsilon }\mathrm{(}\mathit{\xi }+{\mathit{Y}}^0\mathrm{)}}|sym\nabla {\mathit{\eta }}_{\mathit{\epsilon }}{|}^{\mathit{p}}\mathit{d}\mathit{x}\lesssim {{\displaystyle \int }}_{\mathit{\epsilon }\mathrm{(}\mathit{\xi }+{\mathit{Y}}^1\mathrm{)}}|sym\nabla {\mathit{\eta }}_{\mathit{\epsilon }}{|}^{\mathit{p}}\mathit{d}\mathit{x}.\end{array}$$ For the argument consider the rescaled function $$\begin{array}{c}{\hat{\mathit{\eta }}}_{\mathit{\epsilon }}:\mathit{Y}\mapsto {\mathbb{R}}^{\mathit{d}}\mathrm{,}{\hat{\mathit{\eta }}}_{\mathit{\epsilon }}\mathrm{(}\mathit{y}\mathrm{)}:={\mathit{\eta }}_{\mathit{\epsilon }}\mathrm{(}\mathit{\epsilon }\mathrm{(}\mathit{\xi }+\mathit{y}\mathrm{)}\mathrm{)}-\mathit{S}\mathit{y}+\mathit{c}\mathrm{,}\end{array}$$ where $\mathit{S}\in {\mathbb{R}}_{\mathrm{skew}}^{\mathit{d}\times \mathit{d}}$ and $\mathit{c}\in {\mathbb{R}}^{\mathit{d}}$ are chosen such that the Poincaré and Korn inequalities yield $$\begin{array}{}\text{(75)}& {{\displaystyle \int }}_{{\mathit{Y}}^1}\mathrm{(}|{\hat{\mathit{\eta }}}_{\mathit{\epsilon }}{|}^{\mathit{p}}+|\nabla {\hat{\mathit{\eta }}}_{\mathit{\epsilon }}{|}^{\mathit{p}}\mathrm{)}\mathit{d}\mathit{y}\lesssim {{\displaystyle \int }}_{{\mathit{Y}}^1}|sym\nabla {\hat{\mathit{\eta }}}_{\mathit{\epsilon }}{|}^{\mathit{p}}\mathit{d}\mathit{y}.\end{array}$$ Since both ${\mathit{g}}_{\mathit{\epsilon }}^1$ and ${\mathit{g}}_{\ast \mathrm{,}\mathit{\epsilon }}^1$ satisfy (16)(ii), we have $-△{\hat{\mathit{\eta }}}_{\mathit{\epsilon }}=0$ in ${\mathit{Y}}^0$ in the distributional sense. Hence, thanks to Assumption 1 and (75), we have $$\begin{array}{c}{{\displaystyle \int }}_{{\mathit{Y}}^0}|\nabla {\hat{\mathit{\eta }}}_{\mathit{\epsilon }}{|}^{\mathit{p}}\mathit{d}\mathit{y}\lesssim {{\displaystyle \int }}_{{\mathit{Y}}^1}|sym\nabla {\hat{\mathit{\eta }}}_{\mathit{\epsilon }}{|}^{\mathit{p}}\mathit{d}\mathit{y}\mathrm{,}\end{array}$$ and thus $$\begin{array}{rcl}{{\displaystyle \int }}_{\mathit{\epsilon }\mathrm{(}\mathit{\xi }+{\mathit{Y}}^0\mathrm{)}}|sym\nabla {\mathit{\eta }}_{\mathit{\epsilon }}{|}^{\mathit{p}}\mathit{d}\mathit{x}& =& {\mathit{\epsilon }}^{\mathit{d}-\mathit{p}}{{\displaystyle \int }}_{{\mathit{Y}}^0}|sym\nabla {\hat{\mathit{\eta }}}_{\mathit{\epsilon }}{|}^{\mathit{p}}\mathit{d}\mathit{y}\lesssim {\mathit{\epsilon }}^{\mathit{d}-\mathit{p}}{{\displaystyle \int }}_{{\mathit{Y}}^1}|sym\nabla {\hat{\mathit{\eta }}}_{\mathit{\epsilon }}{|}^{\mathit{p}}\mathit{d}\mathit{y}\\ & =& {{\displaystyle \int }}_{\mathit{\epsilon }\mathrm{(}\mathit{\xi }+{\mathit{Y}}^1\mathrm{)}}|sym\nabla {\mathit{\eta }}_{\mathit{\epsilon }}{|}^{\mathit{p}}\mathit{d}\mathit{x}.\end{array}$$ □

Proof of Proposition 2.

Step 1. Proof of (a).

Arguing as in Step 1 in the proof of Proposition 1 we find that ${lim\hspace{0.17em}sup}_{\mathit{\epsilon }\downarrow 0}{\mathit{\epsilon }}^{-2\mathit{\alpha }}{\mathcal{E}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{\epsilon }}^{\mathit{\alpha }}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}\mathrm{<}\infty$ , and thus the compactness part of Theorem 1 (and the fact that two-scale convergence implies weak convergence) yields (for a subsequence, which we  do not relabel): $$\begin{array}{}\text{(76)}& {\mathit{\phi }}_{\mathit{\epsilon }}\stackrel{2}{\mathrm{\rightharpoonup }}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}\text{and}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{\rightharpoonup }\mathit{\phi }:={\mathit{g}}^1+{{\displaystyle \int }}_{{\mathit{Y}}^0}{\mathit{g}}^0\mathit{d}\mathit{y}.\end{array}$$

Step 2. Proof of (b).

We may restrict to the case $$\begin{array}{c}\underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}inf}{\mathit{I}}_{\mathit{\epsilon }}^{\mathit{\alpha }}\mathrm{(}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}=\underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}sup}{\mathit{I}}_{\mathit{\epsilon }}^{\mathit{\alpha }}\mathrm{(}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}\mathrm{<}\infty .\end{array}$$ Thanks to Step 1 we may assume without loss of generality that (76) holds. From the lower bound part of Theorem 1 (and the fact that ${\int }_{\mathrm{\Omega }}\mathit{f}·{\mathit{\phi }}_{\mathit{\epsilon }}\mathit{d}\mathit{x}\to {\int }_{\mathrm{\Omega }}\mathit{f}·\mathit{\phi }\mathit{d}\mathit{x}$ ) we thus deduce that $$\begin{array}{c}\underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}inf}{\mathit{I}}_{\mathit{\epsilon }}^{\mathit{\alpha }}\mathrm{(}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}⩾{\mathit{I}}_{\text{small}}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}.\end{array}$$ With the definition of $\bar{\mathit{Q}}$ from the introduction and $\bar{\mathit{G}}:={\int }_{{\mathit{Y}}_0}{\mathit{g}}^0\mathit{d}\mathit{y}$ we get the inequality: $$\begin{array}{rcl}{\mathit{I}}_{\text{small}}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}& ⩾& {{\displaystyle \int }}_{\mathrm{\Omega }}{\mathit{Q}}_{hom}^1\mathrm{(}\nabla {\mathit{g}}^1\mathrm{(}\mathit{x}\mathrm{)}\mathrm{)}-\mathrm{(}{\mathit{g}}^1\mathrm{(}\mathit{x}\mathrm{)}+\bar{\mathit{G}}\mathrm{(}\mathit{x}\mathrm{)}\mathrm{)}·\mathit{f}\mathrm{(}\mathit{x}\mathrm{)}\mathit{d}\mathit{x}\\ & & +{{\displaystyle \int }}_{\mathrm{\Omega }}min\left\{{{\displaystyle \int }}_{{\mathit{Y}}^0}{\mathit{Q}}^0\mathrm{\right(}{\nabla }_{\mathit{y}}{\tilde{\mathit{g}}}^0\mathrm{(}\mathit{y}\mathrm{)}\mathrm{)}:{\tilde{\mathit{g}}}^0\in {\mathit{H}}_0^1\mathrm{(}{\mathit{Y}}^0\mathrm{)}\mathrm{,}{{\displaystyle \int }}_{{\mathit{Y}}^0}{\tilde{\mathit{g}}}^0\mathrm{(}\mathit{y}\mathrm{)}=\bar{\mathit{G}}\mathrm{(}\mathit{x}\mathrm{)}\}\\ & =& {{\displaystyle \int }}_{\mathrm{\Omega }}{\mathit{Q}}_{hom}^1\mathrm{(}\nabla {\mathit{g}}^1\mathrm{(}\mathit{x}\mathrm{)}\mathrm{)}+\bar{\mathit{Q}}\mathrm{(}\bar{\mathit{G}}\mathrm{(}\mathit{x}\mathrm{)}\mathrm{)}-\mathrm{(}{\mathit{g}}^1\mathrm{(}\mathit{x}\mathrm{)}+\bar{\mathit{G}}\mathrm{(}\mathit{x}\mathrm{)}\mathrm{)}·\mathit{f}\mathrm{(}\mathit{x}\mathrm{)}\mathit{d}\mathit{x}.\end{array}$$ Since $\mathit{\phi }={\mathit{g}}^1+\overline{\mathit{G}}$ by (76), the right-hand side is bounded from below by ${\overline{\mathit{I}}}_{\text{small}}\mathrm{(}\mathit{\phi }\mathrm{)}$ , which completes the argument for (b).

Step 3. Proof of (c).

Let $\mathit{\phi }\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}$ . It suffices to argue that there exist ${\mathit{g}}^0\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{,}{\mathit{H}}_0^1\mathrm{(}{\mathit{Y}}^0\mathrm{)}\mathrm{)}$ and ${\mathit{g}}^1\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ with $\mathit{\phi }={\mathit{g}}^1+{\int }_{{\mathit{Y}}^0}{\mathit{g}}^0\mathit{d}\mathit{y}$ and $$\begin{array}{}\text{(77)}& {\overline{\mathit{I}}}_{\text{small}}\mathrm{(}\mathit{\phi }\mathrm{)}={\mathit{I}}_{\text{small}}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}.\end{array}$$ Indeed, in that case, we can find by part (c) of Theorem 1 a sequence ${\mathit{\phi }}_{\mathit{\epsilon }}\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ such that ${\mathit{\phi }}_{\mathit{\epsilon }}\stackrel{2}{\mathrm{\rightharpoonup }}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}$ and ${lim}_{\mathit{\epsilon }\downarrow 0}{\mathit{I}}_{\mathit{\epsilon }}^{\mathit{\alpha }}\mathrm{(}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}={\mathit{I}}_{\text{small}}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}$ . Since ${\mathit{\phi }}_{\mathit{\epsilon }}\stackrel{2}{\mathrm{\rightharpoonup }}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}$ implies ${\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{\rightharpoonup }\mathit{\phi }={\mathit{g}}^1+{\int }_{{\mathit{Y}}^0}{\mathit{g}}^0\mathit{d}\mathit{y}$ weakly in ${\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}$ , we deduce from (77) that ${\mathit{\phi }}_{\mathit{\epsilon }}$ is the sought for recovery sequence. In order to prove (77), let ${\mathit{\rho }}_{\mathit{i}}$ ( $\mathit{i}=1\mathrm{,}2\mathrm{,}3$ ) denote the unique minimiser in ${\mathit{H}}_0^1\mathrm{(}{\mathit{Y}}^0\mathrm{)}$ to the functional $$\begin{array}{c}{\mathit{H}}_0^1\mathrm{\left(}{\mathit{Y}}^0\mathrm{\right)}\ni \mathit{\rho }\mapsto {{\displaystyle \int }}_{{\mathit{Y}}^0}{\mathit{Q}}^0\mathrm{(}\nabla \mathit{\rho }\mathrm{(}\mathit{y}\mathrm{)}\mathrm{)}\mathit{d}\mathit{y}\text{subject to }{{\displaystyle \int }}_{{\mathit{Y}}^0}\mathit{\rho }\mathit{d}\mathit{y}={\mathit{e}}_{\mathit{i}}\mathrm{,}\end{array}$$ and set $$\begin{array}{c}{\overline{\mathit{q}}}_{\mathit{i}}:={{\displaystyle \int }}_{{\mathit{Y}}^0}{\mathit{Q}}^0\mathrm{(}\nabla {\mathit{\rho }}_{\mathit{i}}\mathrm{(}\mathit{y}\mathrm{)}\mathrm{)}\mathit{d}\mathit{y}.\end{array}$$ Then it is easy to check that $$\begin{array}{c}\bar{\mathit{Q}}\mathrm{(}\mathit{G}\mathrm{)}={\displaystyle \sum }_{\mathit{i}=1}^3{\mathit{G}}_{\mathit{i}}^2{\overline{\mathit{q}}}_{\mathit{i}}\text{for all }\mathit{G}\in {\mathbb{R}}^3\mathrm{,}\end{array}$$ and since ${\overline{\mathit{q}}}_1\mathrm{,}\dots \mathrm{,}{\overline{\mathit{q}}}_3\mathrm{>}0$ , we deduce that $\bar{\mathit{Q}}$ is a positive definite quadratic form. Hence, since ${\mathit{Q}}_{hom}^1\mathrm{(}\mathit{F}\mathrm{)}⩾\mathit{c}|sym\mathit{F}{|}^2$ for some $\mathit{c}\mathrm{>}0$ , we deduce that we can find a unique function ${\mathit{g}}^1\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ that minimizes the functional $$\begin{array}{}\text{(78)}& {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}\ni \mathit{g}\mapsto {{\displaystyle \int }}_{\mathrm{\Omega }}{\mathit{Q}}_{hom}^1\mathrm{(}\nabla \mathit{g}\mathrm{(}\mathit{x}\mathrm{)}\mathrm{)}+\bar{\mathit{Q}}\mathrm{(}\mathit{\phi }\mathrm{(}\mathit{x}\mathrm{)}-\mathit{g}\mathrm{(}\mathit{x}\mathrm{)}\mathrm{)}\mathit{d}\mathit{x}.\end{array}$$ Setting ${\mathit{g}}^0\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}:={\sum }_{\mathit{i}=1}^3\mathrm{(}\mathrm{(}\mathit{\phi }\mathrm{(}\mathit{x}\mathrm{)}-\mathit{g}\mathrm{(}\mathit{x}\mathrm{)}\mathrm{)}·{\mathit{e}}_{\mathit{i}}\mathrm{)}{\mathit{\rho }}_{\mathit{i}}\mathrm{(}\mathit{y}\mathrm{)}$ we deduce that ${\int }_{\mathrm{\Omega }}\bar{\mathit{Q}}\mathrm{(}\mathit{\phi }\mathrm{(}\mathit{x}\mathrm{)}-\mathit{g}\mathrm{(}\mathit{x}\mathrm{)}\mathrm{)}\mathit{d}\mathit{x}={\int }_{\mathrm{\Omega }\times {\mathit{Y}}^0}{\mathit{Q}}^0\mathrm{(}{\nabla }_{\mathit{y}}{\mathit{g}}^0\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}\mathrm{)}\mathit{d}\mathit{y}$ , and thus (thanks to the definition of ${\overline{\mathit{I}}}_{\text{small}}$ ) (77) follows.

Step 4. Proof of (d).

By Proposition 1 we have ${\mathit{\phi }}_{\mathit{\epsilon }}\stackrel{2}{\mathrm{\rightharpoonup }}\mathrm{(}{\mathit{g}}_{\ast }^0\mathrm{,}{\mathit{g}}_{\ast }^1\mathrm{)}$ and thus ${\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{\rightharpoonup }\tilde{\mathit{\phi }}:={\mathit{g}}_{\ast }^1+{\int }_{{\mathit{Y}}^0}{\mathit{g}}_{\ast }^0\mathit{d}\mathit{y}$ . By definition of ${\overline{\mathit{I}}}_{\text{small}}$ we have $$\begin{array}{}\text{(79)}& {\mathit{I}}_{\text{small}}\mathrm{(}{\mathit{g}}_{\ast }^0\mathrm{,}{\mathit{g}}_{\ast }^1\mathrm{)}⩾{\overline{\mathit{I}}}_{\text{small}}\mathrm{(}\tilde{\mathit{\phi }}\mathrm{)}⩾\underset{{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}}{inf}{\overline{\mathit{I}}}_{\text{small}}.\end{array}$$ On the other hand, the map ${\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}\ni \mathit{\phi }\mapsto \mathit{g}\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ with g minimizing the functional in (78) is linear and bounded; hence, we deduce that ${\overline{\mathit{I}}}_{\text{small}}$ is quadratic and strictly convex. It thus admits a unique minimiser ${\mathit{\phi }}_{\ast }\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}$ . By Step 3 (cf. (77)) we associated with ${\mathit{\phi }}_{\ast }$ a pair $\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}$ such that ${\overline{\mathit{I}}}_{\text{small}}\mathrm{(}{\mathit{\phi }}_{\ast }\mathrm{)}={\mathit{I}}_{\text{small}}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}$ . Combined with (79) we get $$\begin{array}{c}\underset{{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}}{min}{\overline{\mathit{I}}}_{\text{small}}={\overline{\mathit{I}}}_{\text{small}}\mathrm{(}{\mathit{\phi }}_{\ast }\mathrm{)}={\mathit{I}}_{\text{small}}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}⩾{\mathit{I}}_{\text{small}}\mathrm{(}{\mathit{g}}_{\ast }^0\mathrm{,}{\mathit{g}}_{\ast }^1\mathrm{)}⩾{\overline{\mathit{I}}}_{\text{small}}\mathrm{(}\tilde{\mathit{\phi }}\mathrm{)}⩾\underset{{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}}{min}{\overline{\mathit{I}}}_{\text{small}}.\end{array}$$ Hence, equality holds everywhere and the claimed identities follow from the strict convexity of ${\overline{\mathit{I}}}_{\text{small}}$ . The convergence of ${inf}_{{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}}{\mathit{I}}_{\mathit{\epsilon }}^{\mathit{\alpha }}\to {min}_{{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}}{\overline{\mathit{I}}}_{\text{small}}$ follows from part (b) and (c) by standard arguments from theory of Γ-convergence. □

4.3. Proof of Theorem 3: Finite strain regime

We define, for ${\mathit{g}}_{\mathit{\epsilon }}^0\in {\mathit{H}}_0^1\mathrm{(}{\mathrm{\Omega }}_{\mathit{\epsilon }}^0\mathrm{)}$ , ${\mathit{g}}^0\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{,}{\mathit{H}}_0^1\mathrm{(}{\mathit{Y}}^0\mathrm{)}\mathrm{)}$ , and ${\mathit{g}}_{\mathit{\epsilon }}^1\mathrm{,}{\mathit{g}}^1\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ , the following functionals: $$\begin{array}{rcl}{\mathit{I}}_{\mathit{\epsilon }}^0\mathrm{\left(}{\mathit{g}}_{\mathit{\epsilon }}^0\mathrm{\right)}& :=& {{\displaystyle \int }}_{{\mathrm{\Omega }}_{\mathit{\epsilon }}^0}\mathrm{\left(}{\mathit{W}}^0\mathrm{\right(}\mathit{I}+\mathit{\epsilon }\nabla {\mathit{g}}_{\mathit{\epsilon }}^0\mathrm{(}\mathit{x}\mathrm{)}\mathrm{)}-{\mathit{\epsilon }}^{1-2\mathit{\gamma }}{\ell }_{\mathit{\epsilon }}·{\mathit{g}}_{\mathit{\epsilon }}^0\mathrm{)}\mathit{d}\mathit{x}\mathrm{,}\\ {\mathit{I}}_0^0\mathrm{\left(}{\mathit{g}}^0\mathrm{\right)}& :=& {{\displaystyle \int }}_{\mathrm{\Omega }\times {\mathit{Y}}^0}\mathrm{\left(}\mathcal{Q}{\mathit{W}}^0\mathrm{\right(}\mathit{I}+{\nabla }_{\mathit{y}}{\mathit{g}}^0\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}\mathrm{)}-{\ell }^0·{\mathit{g}}^0\mathrm{)}\mathit{d}\mathit{x}\mathit{d}\mathit{y}\mathrm{,}\\ {\mathit{I}}_{\mathit{\epsilon }}^1\mathrm{\left(}{\mathit{g}}_{\mathit{\epsilon }}^1\mathrm{\right)}& :=& {\mathit{\epsilon }}^{-2\mathit{\gamma }}{{\displaystyle \int }}_{{\mathrm{\Omega }}_{\mathit{\epsilon }}^1}\mathrm{\left(}{\mathit{W}}^1\mathrm{\right(}\mathit{I}+{\mathit{\epsilon }}^{\mathit{\gamma }}\nabla {\mathit{g}}_{\mathit{\epsilon }}^1\mathrm{(}\mathit{x}\mathrm{)}\mathrm{)}-{{\displaystyle \int }}_{\mathrm{\Omega }}{\mathit{\epsilon }}^{-\mathit{\gamma }}{\ell }_{\mathit{\epsilon }}·{\mathit{g}}_{\mathit{\epsilon }}^1\mathrm{)}\mathit{d}\mathit{x}\mathrm{,}\\ {\mathit{I}}_0^1\mathrm{\left(}{\mathit{g}}^1\mathrm{\right)}& :=& {{\displaystyle \int }}_{\mathrm{\Omega }}\mathrm{\left(}{\mathit{Q}}_{hom}^1\mathrm{\right(}\nabla {\mathit{g}}^1\mathrm{(}\mathit{x}\mathrm{)}\mathrm{)}-{\ell }^1·{\mathit{g}}^1\mathrm{)}\mathit{d}\mathit{x}.\end{array}$$ Thanks to the Lipschitz condition (32b), we can decompose ${\mathit{I}}_{\mathit{\epsilon }}$ into the sum ${\mathit{I}}_{\mathit{\epsilon }}^0+{\mathit{I}}_{\mathit{\epsilon }}^1$ at the expense of a small error. More precisely, the following lemma holds.

Lemma 10.

Suppose that ${\mathit{m}}_{\mathit{\epsilon }}={\mathit{\epsilon }}^{\mathit{\gamma }}$ . There exists a constant $\mathit{C}\mathrm{>}0$ such that for all $\mathit{\epsilon }\mathrm{>}0$ and ${\mathit{\phi }}_{\mathit{\epsilon }}\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ , ${\mathit{\epsilon }}^{1-\mathit{\gamma }}{\mathit{g}}_{\mathit{\epsilon }}^0+{\mathit{g}}_{\mathit{\epsilon }}^1\stackrel{\mathit{\text{(16)}}}{:=}{\mathit{\phi }}_{\mathit{\epsilon }}$ we have $$\begin{array}{c}|{\mathit{I}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}-\mathrm{(}{\mathit{I}}_{\mathit{\epsilon }}^0\mathrm{\left(}{\mathit{g}}_{\mathit{\epsilon }}^0\mathrm{\right)}+{\mathit{I}}_{\mathit{\epsilon }}^1\mathrm{\left(}{\mathit{g}}_{\mathit{\epsilon }}^1\mathrm{\right)}\mathrm{\left)}\right|⩽\mathit{C}{\mathit{\epsilon }}^{\mathit{\gamma }}\mathrm{(}1+{\mathrm{\Phi }}_{\mathit{\epsilon }}^{\mathit{\gamma }}\mathrm{(}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}\mathrm{)}.\end{array}$$

Proof.

Note that $$\begin{array}{ll}\hfill & |{\mathit{I}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}-\mathrm{(}{\mathit{I}}_{\mathit{\epsilon }}^0\mathrm{\left(}{\mathit{g}}_{\mathit{\epsilon }}^0\mathrm{\right)}+{\mathit{I}}_{\mathit{\epsilon }}^1\mathrm{\left(}{\mathit{g}}_{\mathit{\epsilon }}^1\mathrm{\right)}\mathrm{\left)}\right|\\ \hfill & =\left|{{\displaystyle \int }}_{{\mathrm{\Omega }}_{\mathit{\epsilon }}^0}{\mathit{W}}^0\mathrm{\right(}\mathit{I}+{\mathit{\epsilon }}^{\mathit{\gamma }}\mathrm{(}{\mathit{\epsilon }}^{1-\mathit{\gamma }}\nabla {\mathit{g}}_{\mathit{\epsilon }}^0+\nabla {\mathit{g}}_{\mathit{\epsilon }}^1\mathrm{)}\mathrm{)}-{\mathit{W}}^0\mathrm{(}\mathit{I}+\mathit{\epsilon }\nabla {\mathit{g}}_{\mathit{\epsilon }}^0\mathrm{\left)}\mathit{d}\mathit{x}\right|.\end{array}$$ In view of (32b) and (17), the statement follows. □

The following lemma is a simple consequence of [16, Lemma 21, Lemma 22] and (25c):

Lemma 11.

Assume ( 25a ), ( 25c ) and ${\mathit{m}}_{\mathit{\epsilon }}={\mathit{e}}^{\mathit{\gamma }}$ .

Consider a sequence ${\mathit{g}}_{\mathit{\epsilon }}^0\in {\mathit{H}}_0^1\mathrm{(}{\mathrm{\Omega }}_{\mathit{\epsilon }}^0\mathrm{)}$ . If ${\mathit{g}}_{\mathit{\epsilon }}^0\stackrel{2}{\to }{\mathit{g}}^0$ and $\mathit{\epsilon }\nabla {\mathit{g}}_0^1\stackrel{2}{\mathrm{\rightharpoonup }}{\nabla }_{\mathit{y}}{\mathit{g}}^0$ , then $$\begin{array}{c}\underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}inf}{\mathit{I}}_{\mathit{\epsilon }}^0\mathrm{\left(}{\mathit{g}}_{\mathit{\epsilon }}^0\mathrm{\right)}⩾{\mathit{I}}_0^0\mathrm{\left(}{\mathit{g}}^0\mathrm{\right)}.\end{array}$$

For all ${\mathit{g}}^0\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{,}{\mathit{H}}_0^1\mathrm{(}{\mathit{Y}}^0\mathrm{)}\mathrm{)}$ there exists a sequence ${\mathit{g}}_{\mathit{\epsilon }}^0\in {\mathit{H}}_0^1\mathrm{(}{\mathrm{\Omega }}_{\mathit{\epsilon }}^0\mathrm{)}$ such that ${\mathit{g}}_{\mathit{\epsilon }}^0\stackrel{2}{\to }{\mathit{g}}^0$ , $\mathit{\epsilon }\nabla {\mathit{g}}_{\mathit{\epsilon }}^0\stackrel{2}{\mathrm{\rightharpoonup }}{\nabla }_{\mathit{y}}{\mathit{g}}^0$ , and $$\begin{array}{c}\underset{\mathit{\epsilon }\downarrow 0}{lim}{\mathit{I}}_{\mathit{\epsilon }}^0\mathrm{\left(}{\mathit{g}}_{\mathit{\epsilon }}^0\mathrm{\right)}={\mathit{I}}_0^0\mathrm{\left(}{\mathit{g}}^0\mathrm{\right)}.\end{array}$$

For the stiff part one can prove (similar to Lemma 6) the following lemma:

Lemma 12.

Assume that ( 25a )–( 25c ) hold.

Consider ${\mathit{g}}_{\mathit{\epsilon }}^1\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ . If ${\mathit{g}}_{\mathit{\epsilon }}^1\mathrm{\rightharpoonup }{\mathit{g}}^1$ weakly in ${\mathit{H}}^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ , then $$\begin{array}{c}\underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}inf}{\mathit{I}}_{\mathit{\epsilon }}^1\mathrm{\left(}{\mathit{g}}_{\mathit{\epsilon }}^1\mathrm{\right)}⩾{\mathit{I}}_0^1\mathrm{\left(}{\mathit{g}}^1\mathrm{\right)}.\end{array}$$

For all ${\mathit{g}}^1\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ there exists a sequence ${\mathit{g}}_{\mathit{\epsilon }}^1\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ such that $$\begin{array}{c}{\mathit{g}}_{\mathit{\epsilon }}^1\mathrm{\rightharpoonup }{\mathit{g}}^1\mathit{\text{weakly in }}{\mathit{H}}^1\mathrm{(}\mathrm{\Omega }\mathrm{)}\mathrm{,}\mathit{\text{and}}\underset{\mathit{\epsilon }\downarrow 0}{lim}{\mathit{I}}_{\mathit{\epsilon }}^1\mathrm{\left(}{\mathit{g}}_{\mathit{\epsilon }}^1\mathrm{\right)}={\mathit{I}}_0^1\mathrm{\left(}{\mathit{g}}^1\mathrm{\right)}.\end{array}$$

We proceed to the proof of Theorem 3.

Proof of Theorem 3.

Step 1. The proof of parts (a) and (b).

Statement (a) and (b) directly follow from Lemma 10, Lemma 11 and Lemma 12.

Step 2. Proof of (33): convergence of the minima.

For brevity set $$\begin{array}{l}\hfill {\mathit{e}}_{\mathit{\epsilon }}:=\underset{\mathit{\phi }\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}}{inf}{\mathit{\epsilon }}^{-2\mathit{\gamma }}{\mathit{I}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}\mathrm{,}{\mathit{e}}_0:=\underset{{\mathit{g}}^1\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}}{\underset{{\mathit{g}}^0\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{,}{\mathit{H}}_0^1\mathrm{(}{\mathit{Y}}^0\mathrm{)}\mathrm{)}}{inf}}{\mathit{I}}_{\text{finite}}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}.\end{array}$$ We prove (33) in the form of the two inequalities $$\begin{array}{}\text{(80)}& \underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}sup}{\mathit{e}}_{\mathit{\epsilon }}⩽{\mathit{e}}_0\mathrm{,}\underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}inf}{\mathit{e}}_{\mathit{\epsilon }}⩾{\mathit{e}}_0.\end{array}$$ The argument for the first inequality in (80) is standard: for $\mathit{\delta }\mathrm{>}0$ choose $\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}$ with ${\mathit{I}}_{\text{finite}}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}⩽{\mathit{e}}_0+\mathit{\delta }$ . By part (b) there exists a recovery sequence ${\mathit{\phi }}_{\mathit{\epsilon }}$ , so that ${\mathit{I}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}\to {\mathit{I}}_{\text{finite}}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}$ . Hence $$\begin{array}{c}\underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}sup}{\mathit{e}}_{\mathit{\epsilon }}⩽\underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}sup}{\mathit{I}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}={\mathit{I}}_{\text{finite}}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}⩽{\mathit{e}}_0+\mathit{\delta }.\end{array}$$ Since this is valid for all $\mathit{\delta }\mathrm{>}0$ , the first inequality in (80) follows.

Next, we prove the second inequality in (80). Let ${\mathit{\phi }}_{\mathit{\epsilon }}$ denote a sequence with the property ${lim\hspace{0.17em}inf}_{\mathit{\epsilon }\downarrow 0}{\mathit{e}}_{\mathit{\epsilon }}={lim\hspace{0.17em}inf}_{\mathit{\epsilon }\downarrow 0}{\mathit{I}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}$ ; e.g. choose ${\mathit{\phi }}_{\mathit{\epsilon }}\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}$ such that ${\mathit{I}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}⩽{\mathit{e}}_{\mathit{\epsilon }}+\mathit{\epsilon }$ . Combining this with Lemma 10 we deduce that $$\begin{array}{c}\underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}inf}{\mathit{e}}_{\mathit{\epsilon }}=\underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}inf}{\mathit{I}}_{\mathit{\epsilon }}\mathrm{(}{\mathit{\phi }}_{\mathit{\epsilon }}\mathrm{)}=\underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}inf}\mathrm{\left(}{\mathit{I}}_{\mathit{\epsilon }}^0\mathrm{\right(}{\mathit{g}}_{\mathit{\epsilon }}^0\mathrm{)}+{\mathit{I}}_{\mathit{\epsilon }}^1\mathrm{(}{\mathit{g}}_{\mathit{\epsilon }}^1\mathrm{\left)}\mathrm{\right)}⩾\underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}inf}{\mathit{I}}_{\mathit{\epsilon }}^0\mathrm{\left(}{\mathit{g}}_{\mathit{\epsilon }}^0\mathrm{\right)}+\underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}inf}{\mathit{I}}_{\mathit{\epsilon }}^1\mathrm{\left(}{\mathit{g}}_{\mathit{\epsilon }}^1\mathrm{\right)}\mathrm{,}\end{array}$$ where ${\mathit{\epsilon }}^{1-\mathit{\gamma }}{\mathit{g}}_{\mathit{\epsilon }}^0+{\mathit{g}}_{\mathit{\epsilon }}^1\stackrel{\text{(16)}}{:=}{\mathit{\phi }}_{\mathit{\epsilon }}$ . By passing to a subsequence, we assume without loss of generality that ${\mathit{\phi }}_{\mathit{\epsilon }}\stackrel{2}{\mathrm{\rightharpoonup }}\mathrm{(}{\mathit{g}}^0\mathrm{,}{\mathit{g}}^1\mathrm{)}$ . Since, thanks to Lemma 12, we have $$\begin{array}{c}\underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}inf}{\mathit{I}}_{\mathit{\epsilon }}^1\mathrm{\left(}{\mathit{g}}_{\mathit{\epsilon }}^1\mathrm{\right)}⩾{\mathit{I}}_0^1\mathrm{\left(}{\mathit{g}}^1\mathrm{\right)}⩾\underset{{\tilde{\mathit{g}}}^1\in {\mathit{H}}_0^1\mathrm{(}\mathrm{\Omega }\mathrm{)}}{inf}{\mathit{I}}_0^1\mathrm{\left(}{\tilde{\mathit{g}}}^1\mathrm{\right)}\mathrm{,}\end{array}$$ it remains to argue that $$\begin{array}{}\text{(81)}& \underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}inf}{\mathit{I}}_{\mathit{\epsilon }}^0\mathrm{\left(}{\mathit{g}}_{\mathit{\epsilon }}^0\mathrm{\right)}⩾\underset{{\tilde{\mathit{g}}}^0\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{,}{\mathit{H}}_0^1\mathrm{(}{\mathit{Y}}^0\mathrm{)}\mathrm{)}}{inf}{\mathit{I}}_0^0\mathrm{\left(}{\tilde{\mathit{g}}}^0\mathrm{\right)}.\end{array}$$ We identify ${\mathit{g}}_{\mathit{\epsilon }}^0$ with its extension by zero to ${\mathbb{R}}^{\mathit{d}}$ , and consider the periodic unfolding of ${\mathit{g}}_{\mathit{\epsilon }}^0$ defined as $$\begin{array}{c}{\mathit{G}}_{\mathit{\epsilon }}^0:\mathrm{\Omega }\times \mathit{Y}\to {\mathbb{R}}^{\mathit{d}}\mathrm{,}{\mathit{G}}_{\mathit{\epsilon }}^0\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}:={\mathit{g}}_{\mathit{\epsilon }}^0\mathrm{(}\mathit{\epsilon }\lfloor \mathit{x}\mathrm{/}\mathit{\epsilon }\rfloor +\mathit{\epsilon }\mathit{y}\mathrm{)}\mathrm{,}\end{array}$$ where $\lfloor \mathit{z}\rfloor$ stands the unique vector in ${\mathbb{Z}}^{\mathit{d}}$ such that ${\mathit{z}-\lfloor \mathit{z}\rfloor \in [0\mathrm{,}1\mathrm{)}}^{\mathit{d}}$ . Further, note that ${\mathit{G}}_{\mathit{\epsilon }}^0\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{,}{\mathit{H}}_0^1\mathrm{(}{\mathit{Y}}^0\mathrm{)}\mathrm{)}$ , and for $\mathit{\xi }\in {\mathit{Z}}_{\mathit{\epsilon }}:=\{\mathit{\xi }\in {\mathbb{Z}}^{\mathit{d}}:\mathit{\epsilon }\mathrm{(}\mathit{\xi }+\mathit{Y}\mathrm{)}\subset \mathrm{\Omega }\}$ and $\mathit{y}\in \mathit{Y}$ one has $$\begin{array}{c}{\mathit{g}}_{\mathit{\epsilon }}^0\mathrm{(}\mathit{\epsilon }\mathrm{(}\mathit{\xi }+\mathit{y}\mathrm{)}\mathrm{)}={\mathit{G}}_{\mathit{\epsilon }}^0\mathrm{(}\mathit{\epsilon }\mathrm{(}\mathit{\xi }+\mathit{y}\mathrm{)}\mathrm{,}\mathit{y}\mathrm{)}\mathrm{,}\mathit{\epsilon }\nabla {\mathit{g}}_{\mathit{\epsilon }}^0\mathrm{(}\mathit{\epsilon }\mathrm{(}\mathit{\xi }+\mathit{y}\mathrm{)}\mathrm{)}={\nabla }_{\mathit{y}}{\mathit{G}}_{\mathit{\epsilon }}^0\mathrm{(}\mathit{\epsilon }\mathrm{(}\mathit{\xi }+\mathit{y}\mathrm{)}\mathrm{,}\mathit{y}\mathrm{)}.\end{array}$$ Now we consider ${\mathit{I}}_{\mathit{\epsilon }}^0\mathrm{(}{\mathit{g}}_{\mathit{\epsilon }}^0\mathrm{)}$ , which involves an integral over the set ${\mathrm{\Omega }}_{\mathit{\epsilon }}^0$ . Since the latter can be written as a union of sets of the form $\mathit{\epsilon }\mathrm{(}\mathit{\xi }+{\mathit{Y}}^0\mathrm{)}$ with $\mathit{\xi }\in {\mathit{Z}}_{\mathit{\epsilon }}$ , an elementary calculation shows that $$\begin{array}{rcl}{\mathit{I}}_{\mathit{\epsilon }}\mathrm{\left(}{\mathit{g}}_{\mathit{\epsilon }}^0\mathrm{\right)}& =& {\displaystyle \sum }_{\mathit{\xi }\in {\mathit{Z}}_{\mathit{\epsilon }}}{{\displaystyle \int }}_{\mathit{\epsilon }\mathrm{(}\mathit{\xi }+{\mathit{Y}}^0\mathrm{)}}\mathrm{\left(}{\mathit{W}}^0\mathrm{\right(}\mathit{I}+\mathit{\epsilon }\nabla {\mathit{g}}_{\mathit{\epsilon }}^0\mathrm{(}\mathit{x}\mathrm{)}\mathrm{)}-{\mathit{\epsilon }}^{2\mathit{\gamma }-1}{\ell }_{\mathit{\epsilon }}\mathrm{(}\mathit{x}\mathrm{)}·{\mathit{g}}_{\mathit{\epsilon }}^0\mathrm{(}\mathit{x}\mathrm{)}\mathrm{)}\mathit{d}\mathit{x}\\ & =& {\displaystyle \sum }_{\mathit{\xi }\in {\mathit{Z}}_{\mathit{\epsilon }}}{\mathit{\epsilon }}^{\mathit{d}}{{\displaystyle \int }}_{{\mathit{Y}}^0}\mathrm{\left(}{\mathit{W}}^0\mathrm{\right(}\mathit{I}+{\nabla }_{\mathit{y}}{\mathit{G}}_{\mathit{\epsilon }}^0\mathrm{(}\mathit{\epsilon }\mathrm{(}\mathit{\xi }+\mathit{y}\mathrm{)}\mathrm{,}\mathit{y}\mathrm{)}\mathrm{)}-{\mathit{\epsilon }}^{2\mathit{\gamma }-1}{\ell }_{\mathit{\epsilon }}\mathrm{(}\mathit{\epsilon }\mathrm{(}\mathit{\xi }+\mathit{y}\mathrm{)}\mathrm{)}·{\mathit{G}}_{\mathit{\epsilon }}^0\mathrm{(}\mathit{\epsilon }\mathrm{(}\mathit{\xi }+\mathit{y}\mathrm{)}\mathrm{,}\mathit{y}\mathrm{\left)}\mathrm{\right)}\mathit{d}\mathit{y}\\ & =& {{\displaystyle \int }}_{\mathrm{\Omega }}{{\displaystyle \int }}_{{\mathit{Y}}^0}\mathrm{\left(}{\mathit{W}}^0\mathrm{\right(}\mathit{I}+{\nabla }_{\mathit{y}}{\mathit{G}}_{\mathit{\epsilon }}^0\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}\mathrm{)}-{\mathit{\epsilon }}^{2\mathit{\gamma }-1}{\ell }_{\mathit{\epsilon }}\mathrm{(}\mathit{\epsilon }\mathrm{(}\lfloor \mathit{x}\mathrm{/}\mathit{\epsilon }\rfloor +\mathit{y}\mathrm{)}\mathrm{)}·{\mathit{G}}_{\mathit{\epsilon }}^0\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}\mathrm{)}\mathit{d}\mathit{y}\mathit{d}\mathit{x}.\end{array}$$ Thanks to (25c), the characterisation of strong two-scale convergence introduced in Step 1 of the proof of Lemma 3, and the fact that ${lim\hspace{0.17em}sup}_{\mathit{\epsilon }\downarrow 0}\Vert {\mathit{g}}_{\mathit{\epsilon }}^0{\Vert }_{{\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{)}}\mathrm{<}\infty$ , we have $$\begin{array}{c}\underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}sup}\left|{{\displaystyle \int }}_{\mathrm{\Omega }\times {\mathit{Y}}^0}\mathrm{\right(}{\mathit{\epsilon }}^{2\mathit{\gamma }-1}{\ell }_{\mathit{\epsilon }}\mathrm{\left(}\mathit{\epsilon }\mathrm{\right(}\lfloor \mathit{x}\mathrm{/}\mathit{\epsilon }\rfloor +\mathit{y}\mathrm{\left)}\mathrm{\right)}-{\ell }^0\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}\mathrm{)}·{\mathit{G}}_{\mathit{\epsilon }}^0\mathrm{(}\mathit{\epsilon }\mathrm{(}\lfloor \mathit{x}\mathrm{/}\mathit{\epsilon }\rfloor +\mathit{y}\mathrm{)}\mathrm{,}\mathit{y}\mathrm{\left)}\mathit{d}\mathit{x}\mathit{d}\mathit{y}\right|=0.\end{array}$$ Hence, one has $$\begin{array}{rcl}\underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}inf}{\mathit{I}}_{\mathit{\epsilon }}\mathrm{\left(}{\mathit{g}}_{\mathit{\epsilon }}^0\mathrm{\right)}& =& \underset{\mathit{\epsilon }\downarrow 0}{lim\hspace{0.17em}inf}{{\displaystyle \int }}_{\mathrm{\Omega }}{{\displaystyle \int }}_{{\mathit{Y}}^0}\mathrm{\left(}{\mathit{W}}^0\mathrm{\right(}\mathit{I}+{\nabla }_{\mathit{y}}{\mathit{G}}_{\mathit{\epsilon }}^0\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}\mathrm{)}-{\ell }^0\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}·{\mathit{G}}_{\mathit{\epsilon }}^0\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}\mathrm{)}\mathit{d}\mathit{y}\mathit{d}\mathit{x}\\ & ⩾& \underset{{\mathit{g}}^0\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{,}{\mathit{H}}_0^1\mathrm{(}{\mathit{Y}}^0\mathrm{)}\mathrm{)}}{inf}{{\displaystyle \int }}_{\mathrm{\Omega }}{{\displaystyle \int }}_{{\mathit{Y}}^0}\mathrm{\left(}{\mathit{W}}^0\mathrm{\right(}\mathit{I}+{\nabla }_{\mathit{y}}{\mathit{g}}^0\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}\mathrm{)}-{\ell }^0\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}·{\mathit{g}}^0\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}\mathrm{)}\mathit{d}\mathit{y}\mathit{d}\mathit{x}\\ & ⩾& \underset{{\mathit{g}}^0\in {\mathit{L}}^2\mathrm{(}\mathrm{\Omega }\mathrm{,}{\mathit{H}}_0^1\mathrm{(}{\mathit{Y}}^0\mathrm{)}\mathrm{)}}{inf}{{\displaystyle \int }}_{\mathrm{\Omega }}{{\displaystyle \int }}_{{\mathit{Y}}^0}\mathrm{\left(}\mathcal{Q}{\mathit{W}}^0\mathrm{\right(}\mathit{I}+{\nabla }_{\mathit{y}}{\mathit{g}}^0\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}\mathrm{)}-{\ell }^0\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}·{\mathit{g}}^0\mathrm{(}\mathit{x}\mathrm{,}\mathit{y}\mathrm{)}\mathrm{)}\mathit{d}\mathit{y}\mathit{d}\mathit{x}\mathrm{,}\end{array}$$ which proves (81).

Step 3. The proof of the convergence (34).

Since $\mathcal{Q}{\mathit{W}}^0$ is quasiconvex and ${\mathit{Q}}_{hom}^1$ quadratic, there exists a pair $\mathrm{(}{\mathit{g}}_{\ast }^0\mathrm{,}{\mathit{g}}_{\ast }^1\mathrm{)}$ that minimises ${\mathit{I}}_{\text{finite}}$ . Now the sequence associated with $\mathrm{(}{\mathit{g}}_{\ast }^0\mathrm{,}{\mathit{g}}_{\ast }^1\mathrm{)}$ via Theorem 3(b) satisfies (34). □