$$ \newcommand{\C}{{\mathbb{{C}}}} \newcommand{\R}{{\mathbb{{R}}}} \newcommand{\Q}{{\mathbb{{Q}}}} \newcommand{\Z}{{\mathbb{{Z}}}} \newcommand{\N}{{\mathbb{{N}}}} \newcommand{\uu}[1]{{\boldsymbol{{#1}}}} \newcommand{\uuuu}[1]{{\symbb{{#1}}}} \newcommand{\uv}[1]{{\underline{{#1}}}} \newcommand{\ve}[1]{{\uv{{e}}_{{#1}}}} \newcommand{\x}{{\uv{{x}}}} \newcommand{\n}{{\uv{{n}}}} \newcommand{\eps}{{\uu{{\varepsilon}}}} \newcommand{\E}{{\uu{{E}}}} \newcommand{\sig}{{\uu{{\sigma}}}} \newcommand{\Sig}{{\uu{{\Sigma}}}} \newcommand{\cod}{{\uv{{\symscr{b}}}}} \newcommand{\trans}[1]{{{}^{t}{#1}}} \newcommand{\sotimes}{{\stackrel{s}{\otimes}}} \newcommand{\sboxtimes}{\stackrel{s}{\boxtimes}} \newcommand{\norm}[1]{{\lVert{{#1}}\rVert}} \newcommand{\ud}{{\,\mathrm{d}}} \DeclareMathOperator{\arcosh}{arcosh} \DeclareMathOperator{\divz}{div} \DeclareMathOperator{\divu}{\uv{div}} \DeclareMathOperator{\hess}{hess} \DeclareMathOperator{\gradu}{\uv{grad}} \DeclareMathOperator{\graduu}{\uu{grad}} \DeclareMathOperator{\Mat}{Mat} \DeclareMathOperator{\tr}{tr} \DeclareMathOperator{\ISO}{ISO} \newcommand{\jump}[1]{[\hspace*{-.15em}[\hspace*{.1em}{#1}% \hspace*{.1em}]\hspace*{-.15em}]} $$

1  Basic problem of elasticity homogenization

1.1 System of equations

Consider a representative volume element (RVE) \(\Omega\) composed of a heterogeneous material. Neglecting body forces in a problem posed at the scale of a RVE is consistent with the fact that the order of magnitude of mechanical effects induced by body forces is in general much lower than that of the macroscopic strain \(\E\) or stress \(\Sig\) effects accounting for interactions with particles surrounding the RVE (see (Dormieux et al., 2006)). The hypothesis of quasi-static equilibrium is also invoked here to write the balance law involving the Cauchy stress field \(\sig\) \[ \divu{\sig}=\uv{0} \quad(\Omega) \tag{1.1}\]

In the sequel, the small perturbation hypothesis is adopted so that the strain field \(\eps\) derives from the displacement one \(\uv{u}\) as the symmetrical part of its gradient \[ \eps=\frac{\graduu{\uv{u}}+\trans{\graduu{\uv{u}}}}{2} \quad(\Omega) \tag{1.2}\]

In the framework of random media homogenization, two types of conditions applied at the boundary \(\partial\Omega\) of a RVE \(\Omega\) are usually considered:

  • homogeneous strain boundary conditions corresponding to prescribed displacements \(\uv{u}^g\) at \(\partial\Omega\) \[ \uv{u}^g=\E\cdot\x \quad(\partial\Omega) \tag{1.3}\] It is noticeable that in this case the divergence theorem implies the following relationship between the microscopic and macroscopic strain tensors \[ <\eps>_{\Omega}=\frac{1}{\lvert \Omega \rvert}\int_\Omega \eps\ud \Omega =\frac{1}{\lvert \Omega \rvert}\int_{\partial\Omega} \uv{u}\sotimes\n \ud S =\E \tag{1.4}\] where the spatial average over a domain \(\omega\) is denoted by \(<\bullet>_\omega\) and \(\n\) is the unit outward normal at the boundary. The macroscopic stress tensor is then simply defined as the average \[ \Sig= <\sig>_{\Omega}=\frac{1}{\lvert \Omega \rvert}\int_\Omega \sig\ud \Omega \tag{1.5}\]

  • homogeneous stress boundary conditions corresponding to prescribed surface tractions \(\uv{T}^g\) at \(\partial\Omega\) \[ \uv{T}^g=\Sig\cdot\n \quad(\partial\Omega) \tag{1.6}\] Now owing to the remarkable identity \((x_i\sigma_{jk})_{,k}=\sigma_{ij}\) resulting from (1.1) and the symmetry of \(\sig\), the relationship between the microscopic and macroscopic stress is ensured by the divergence theorem \[ <\sig>_{\Omega}=\frac{1}{\lvert \Omega \rvert}\int_\Omega \sig\ud \Omega =\frac{1}{\lvert \Omega \rvert}\int_{\partial\Omega} \x\sotimes(\sig\cdot\n) \ud S =\Sig \tag{1.7}\] The macroscopic stress tensor is then simply defined as the average \[ \E= <\eps>_{\Omega}=\frac{1}{\lvert \Omega \rvert}\int_\Omega \eps\ud \Omega \tag{1.8}\]

Hill lemma

Note that whatever the choice of boundary conditions between (1.3) and (1.6), the consistency between the microscopic and macroscopic works is ensured by \[ <\sig:\eps>_{\Omega} =\frac{1}{\lvert \Omega \rvert}\int_\Omega \sig:\eps\ud \Omega =\frac{1}{\lvert \Omega \rvert}\int_{\partial\Omega} \uv{u}\cdot\sig\cdot\n \ud S =\Sig:\E \tag{1.9}\] which results from the application of the divergence theorem to \((u_i\sigma_{ij})_{,k}=u_{i,j}\sigma_{ij}=\varepsilon_{ij}\sigma_{ij}\).

The set of equations defining the problem posed on the RVE is finally completed by the local constitutive law relating the strain and stress fields. The hypothesis of linear elasticity is adopted in this part so that \[ \sig=\uuuu{c}:\eps \quad(\Omega) \tag{1.10}\] where \(\uuuu{c}(\x)\) denotes the heterogeneous (positive definite fourth-order) stiffness tensor field satisfying the conditions of minor (\(c_{jikl}=c_{ijlk}=c_{ijkl}\)) and major (\(c_{klij}=c_{ijkl}\)) symmetries. The compliance tensor field is introduced as the inverse \(\uuuu{s}=\uuuu{c}^{-1}\) in the sense of fourth-order tensors operating over symmetrical second-order tensors.

In short, the system of equations posed on the RVE is given by (1.1), (1.2), (1.3) or (1.6) and (1.10).

1.2 Macroscopic stiffness or compliance tensors

Whatever the boundary condition of homogeneous strain or stress type (1.3) or (1.6), the linearity of the problem allows to invoke the existence of concentration tensors relating the microscopic strain \(\eps\) and stress \(\sig\) fields to the macroscopic strain \(\E\) or stress \(\Sig\) tensors \[\begin{aligned} \eps & = \uuuu{A}_E:\E \\ \sig & = \uuuu{B}_E:\E \quad\textrm{with}\quad \uuuu{B}_E=\uuuu{c}:\uuuu{A}_E \\ \sig & = \uuuu{B}_\Sigma:\Sig \\ \eps & = \uuuu{A}_\Sigma:\Sig \quad\textrm{with}\quad \uuuu{A}_\Sigma=\uuuu{s}:\uuuu{B}_\Sigma \end{aligned} \tag{1.11}\]