Stochastic homogenization of rate-dependent models of monotone type in plasticity

TL;DR

This paper studies the stochastic homogenization of rate-dependent monotone models in plasticity, using Fitzpatrick functions and two-scale convergence to analyze the asymptotic behavior of inelastic materials with microstructure.

Contribution

It introduces a novel approach combining Fitzpatrick functions and stochastic two-scale convergence to analyze rate-dependent monotone models in plasticity.

Findings

Established existence of weak solutions for the models.

Derived homogenization results using stochastic two-scale convergence.

Connected homogenization with classical $\Gamma$-convergence theory.

Abstract

In this work we deal with the stochastic homogenization of the initial boundary value problems of monotone type. The models of monotone type under consideration describe the deformation behaviour of inelastic materials with a microstructure which can be characterised by random measures. Based on the Fitzpatrick function concept we reduce the study of the asymptotic behaviour of monotone operators associated with our models to the problem of the stochastic homogenization of convex functionals within an ergodic and stationary setting. The concept of Fitzpatrick's function helps us to introduce and show the existence of the weak solutions for rate-dependent systems. The derivations of the homogenization results presented in this work are based on the stochastic two-scale convergence in Sobolev spaces. For completeness, we also present some two-scale homogenization results for convex functionals, which are related to the classical $\Gamma$-convergence theory.

Equations

$${\mathbf{Q}}_t=\mathbf{Q}\times (0,t).$$

$$-{\mathrm{div}}_x{\sigma }_{\eta }(x,t)$$

$${\sigma }_{\eta }(x,t)$$

$${\partial }_t{z}_{\eta }(x,t)$$

$$u_{\eta }(x,t)=0,(x,t)\in \partial \mathbf{Q}\times (0,\infty ),$$

$$z_{\eta }(x,0)=z_{\eta }^{(0)}(x),x\in \mathbf{Q}.$$

$$\epsilon ({\nabla }_x{u}_{\eta }(x,t))=\frac{1}{2}({\nabla }_x{u}_{\eta }(x,t)+({\nabla }_x{u}_{\eta }(x,t){)}^T)\in {\mathcal{S}}^3$$

$$B^{T}v=(\hat{z},0{)}^T$$

$$\alpha |\xi {|}^2\le {\mathbb{C}}_{\eta }[x]\xi \cdot \xi \le \beta |\xi {|}^2\text{for}\text{any}\xi \in {\mathcal{S}}^3.$$

$$0\in g_{\eta }(x,0),x\in \mathbf{Q}.$$

$$g_{\eta }(x,\Sigma )=[|\Sigma |-{\sigma }_{\mathrm{y}}(x){]}_{+}^{r_{\eta }(x)}\frac{\Sigma }{|\Sigma |},\Sigma \in {\mathcal{S}}^3,x\in \mathbf{Q},$$

$$\begin{array}{rl}\left[f(\omega )=f({\tau }_x\omega )\forall x\in {\mathbb{R}}^n,a.e.\omega \in \Omega \right]& \Rightarrow \left[f(\omega )=const\mathcal{P}\text{-a.e.}\omega \in \Omega \right].\end{array}$$

$$\mathbb{C}[x,\omega ]=\tilde{\mathbb{C}}[{\tau }_x\omega ],L[x,\omega ]=\tilde{L}[{\tau }_x\omega ],$$

$$g(x,\omega ,v)=\tilde{g}({\tau }_x\omega ,v),$$

$${\mathbb{C}}_{\eta }[x]=\mathbb{C}\left[\frac{x}{\eta },\omega \right],L_{\eta }[x]=L\left[\frac{x}{\eta },\omega \right],$$

$$g_{\eta }(x,v)=g\left(\frac{x}{\eta },\omega ,v\right).$$

$$z_{\eta }^{(0)}(x)={\tilde{z}}^{(0)}\left(x,{\tau }_{\frac{x}{\eta }}\omega \right),x\in \mathbf{Q}.$$

$$(\xi ,\zeta {)}_S={\int }_S(\xi (s),\zeta (s))ds.$$

$$\mathrm{dom}(\varphi )=\{v\in V\mid \varphi (v)<\infty \},\mathrm{epi}(\varphi )=\{(v,t)\in V\times \mathcal{R}\mid \varphi (v)\le t\}$$

$${\varphi }^{\ast }(v^{\ast })=\underset{v\in V}{\sup }\{⟨v^{\ast },v⟩-\varphi (v)\}.$$

$$\forall v\in V,\forall v^{\ast }\in V^{\ast }:{\varphi }^{\ast }(v^{\ast })+\varphi (v)\ge ⟨v^{\ast },v⟩,$$

$$\mathrm{int}\mathrm{dom}(\varphi )=\mathrm{int}\mathrm{dom}(\partial \varphi ),$$

$$I_{\varphi }(v)=\{\begin{array}{ll}{\int }_G\varphi (v(x))dx,& \varphi (v)\in L^1(G,{\mathbb{R}}^k)\\ +\infty ,& \mathrm{otherwise}\end{array},$$

$$v^{\ast }\in L^q(G,{\mathbb{R}}^k),v\in L^p(G,{\mathbb{R}}^k)\mathrm{and}{v}^{\ast }(x)\in \partial \varphi (v(x)),a.e.$$

$$\big{(}I_{\phi}\big{)}^{*}=I_{\phi^{*}},$$

$$D(A)=\{v\in V\mid Av\ne \varnothing \},GrA=\{[v,v^{\ast }]\in V\times V^{\ast }\mid v\in D(A),v^{\ast }\in Av\}$$

$$⟨v^{\ast }-u^{\ast },v-u⟩\ge 0\forall [v,v^{\ast }],[u,u^{\ast }]\in GrA.$$

$$⟨v^{\ast }-u^{\ast },v-u⟩\ge 0\forall [u,u^{\ast }]\in GrA$$

$$\{x\in S\mid A(x)\cap U\ne \varnothing \}$$

$$A(x)=\overline{{\cup }_{n\in \mathbb{N}}{v}_n(x)}.$$

$$Gr{\mathcal{A}}_p=\{[v,v^{\ast }]\in L^p(S,V)\times L^q(S,V^{\ast })\mid [v(x),v^{\ast }(x)]\in GrA(x)fora.e.x\in S\}.$$

Full text

Key words: stochastic homogenization, random measures, plasticity, stochastic two-scale convergence, $\Gamma$-convergence, monotone operator method, Fitzpatrick’s function, Palm measures, random microstructure.

AMS 2000 subject classification: 74Q15, 74C05, 74C10, 74D10, 35J25, 34G20, 34G25, 47H04, 47H05

1 Introduction

In this work we are concerned with the homogenization of the initial boundary value problem describing the deformation behavior of inelastic materials with a microstructure which can be characterised by random measures.

While the periodic homogenization theory for elasto/visco-plastic models is sufficiently well established (see [2, 11, 17, 18, 19, 26, 27, 30, 31] and references therein), some improvement in the development of the techniques for the stochastic homogenization of the quasi-static initial boundary value problems of monotone type has to be achieved yet. To the best knowledge of the authors, there are only two works ([13, 14]) available on the market which are concerned with the homogenization problem of rate-independent systems in plasticity within an ergodic and stationary setting. In this work we extend the results obtained in [14] for perfectly elasto-plastic models to rate-dependent plasticity. Our main ingredient in the construction of the stochastic homogenization theory for rate-dependent models of monotone type is the combination of the Fitzpatrick function concept and the two-scale convergence technique in spaces equipped with random measures due to V.V. Zhikov and A.L. Pyatnitskii (see [34]). The Fitzpatrick function is used here to reduce the study of the asymptotic behaviour of monotone operators associated with the models under consideration to the problem of the stochastic homogenization of convex functionals defined on Sobolev spaces with random measures.

Setting of the problem.

Let $\boldsymbol{Q}\subset\mathbb{R}^{3}$ be an open bounded set, the set of material points of the solid body, with a Lipschitz boundary $\partial\boldsymbol{Q}$, the number $\eta>0$ denote the scaling parameter of the microstructure and $T_{e}$ be some positive number (time of existence). For $0<t\leq T_{e}$

[TABLE]

Let ${\cal S}^{3}$ denote the set of symmetric $3\times 3$-matrices, and let $u_{\eta}(x,t)\in{\mathbb{R}}^{3}$ be the unknown displacement of the material point $x$ at time $t$, $\sigma_{\eta}(x,t)\in{\cal S}^{3}$ be the unknown Cauchy stress tensor and $z_{\eta}(x,t)\in{\mathbb{R}}^{N}$ denote the unknown vector of internal variables. The model equations of the problem (the microscopic problem) are

[TABLE]

together with the homogeneous Dirichlet boundary condition

[TABLE]

and the initial condition

[TABLE]

In model equations (1) - (5)

[TABLE]

denotes the strain tensor (the measure of deformation), $B:{\mathbb{R}}^{N}\rightarrow{\cal S}^{3}$ is a linear mapping, which assigns to each vector of internal variables $z_{\eta}(x,t)$ the plastic strain tensor ${\varepsilon}_{p,\eta}(x,t)\in{\cal S}^{3}$, i.e. the following relation ${\varepsilon}_{p,\eta}(x,t)=Bz_{\eta}(x,t)$ holds. We recall that the space ${\cal S}^{3}$ can be isomorphically identified with the space ${\mathbb{R}}^{6}$ (see [1, p. 31]). Therefore, the linear mapping $B:{\mathbb{R}}^{N}\to{\cal S}^{3}$ is defined as a composition of a projector from ${\mathbb{R}}^{N}$ onto ${\mathbb{R}}^{6}$ and the isomorphism between ${\mathbb{R}}^{6}$ and ${\cal S}^{3}$. The transpose $B^{T}:{\cal S}^{3}\rightarrow{\mathbb{R}}^{N}$ is given by

[TABLE]

for $v\in{\cal S}^{3}$ and $z=(\hat{z},\tilde{z})^{T}\in{\mathbb{R}}^{N},$ $\hat{z}\in{\mathbb{R}}^{6},$ $\tilde{z}\in{\mathbb{R}}^{N-6}$.

For every $x\in\boldsymbol{Q}$ we denote by ${\mathbb{C}}_{\eta}[x]:{\cal S}^{3}\rightarrow{\cal S}^{3}$ a linear symmetric mapping, the elasticity tensor. It is assumed that the mapping $x\rightarrow{\mathbb{C}}_{\eta}[x]$ is measurable. Further, we suppose that there exist two positive constants $0<\alpha<\beta$ such that the two-sided inequality

[TABLE]

is satisfied uniformly with respect to $x\in\boldsymbol{Q}$ and $\eta>0$. The given function $b:\boldsymbol{Q}\times[0,\infty)\rightarrow{\mathbb{R}}^{3}$ is the volume force. The ($N\times N$)-matrix $L_{\eta}[x]$ represents hardening effects. It is assumed to be positive semi-definite, only. For all $x\in\boldsymbol{Q}$ the function $z\rightarrow g_{\eta}(x,z):{\mathbb{R}}^{N}\rightarrow 2^{{\mathbb{R}}^{N}}$ is maximal monotone and satisfies the following condition

[TABLE]

The mapping $x\rightarrow\left(L_{\eta}[x],g_{\eta}(x,\cdot)\right)$ is measurable.

Remark 1.1**.**

Visco-plasticity is typically included in the former conditions by choosing the function $g_{\eta}$ to be in Norton-Hoff form, i.e.

[TABLE]

where $\sigma_{\mathrm{y}}:\boldsymbol{Q}\to(0,\infty)$ is the flow stress function and $r_{\eta}:\boldsymbol{Q}\to(0,\infty)$ is some material function together with $[x]_{+}:=\max(x,0)$.

In order to specify the dependence of the model coefficients in (1) - (5) on the microstructure scaling parameter $\eta>0$, we introduce the concept of a spatial dynamical system. Throughout this paper, we follow the setting of Papanicolaou and Varadhan [22] and make the following assumptions.

Assumption 1.1**.**

*Let $(\Omega,{\cal F}_{\Omega},{\cal P})$ be a probability space with countably generated $\sigma$-algebra ${\cal F}_{\Omega}$. Further, we assume we are given a family $(\tau_{x})_{x\in\mathbb{R}^{n}}$ of measurable bijective mappings $\tau_{x}:\Omega\mapsto\Omega$, having the properties of a *dynamical system on $(\Omega,{\cal F}_{\Omega},{\cal P})$, i.e. they satisfy (i)-(iii):

1. (i)

$\tau_{x}\circ\tau_{y}=\tau_{x+y}$* , $\tau_{0}=id$ (Group property)* 2. (ii)

${\cal P}(\tau_{-x}B)={\cal P}(B)\quad\forall x\in\mathbb{R}^{n},\,\,B\in{\cal F}_{\Omega}$* (Measure preserving)* 3. (iii)

$A:\mathbb{R}^{n}\times\Omega\rightarrow\Omega,\,(x,\omega)\mapsto\tau_{x}\omega$* is measurable (Measurablility of evaluation)*

We finally assume that the system $(\tau_{x})_{x\in\mathbb{R}^{n}}$ is ergodic. This means that for every measurable function $f:\Omega\rightarrow\mathbb{R}$ there holds

[TABLE]

For reader’s convenience, we recall the following well-known result (see [9, Section VI.15]).

Lemma 1.1**.**

Let $(A,{\cal F},\mu)$ be a finite measure space with countably generated $\sigma$-algebra ${\cal F}$. Then, for all $1\leq p<\infty$, $L^{p}(A;\mu)$ contains a countable dense set of simple functions.

The coefficients in (1) - (5) are defined as follows. First, we define the stationary random fields through the relations

[TABLE]

and for every fixed $v\in{\mathbb{R}}^{N}$

[TABLE]

where $\tilde{\mathbb{C}}$, $\tilde{L}$ are measurable functions over $\Omega$ and $\omega\mapsto\tilde{g}(\omega,\cdot)$ is measurable in the sense of Definition 2.2. Then, given the specified assumptions on the random fields, the coefficients ${\mathbb{C}}_{\eta}[x]$, $L_{\eta}[x]$ and the mapping $x\mapsto g_{\eta}(x,\cdot)$ are defined as

[TABLE]

and for each fixed $v\in{\mathbb{R}}^{N}$

[TABLE]

Furthermore, we assume that

[TABLE]

for some ergodic function $\tilde{z}^{(0)}\in L^{2}(\boldsymbol{Q}\times\Omega;{\cal L}\otimes\mu)$.

From a modelling perspective, this construction is equivalent to the assumption that the coefficients and the given functions in (1) - (5) are statistically homogeneous (see [7], for example).

Notation.

The symbols $|\cdot|$ and $(\cdot,\cdot)$ will denote a norm and a scalar product in ${\mathbb{R}}^{k}$, respectively. Let $S$ be a measurable set in ${\mathbb{R}}^{s}$. For $m\in\mathbb{N}$, $q\in[1,\infty]$, we denote by $W^{m,q}(S,{\mathbb{R}}^{k})$ the Banach space of Lebesgue integrable functions having $q$-integrable weak derivatives up to order $m$. This space is equipped with the norm $\|\cdot\|_{m,q,S}$. If $m=0$, we write $\|\cdot\|_{q,S}$; and if (additionally) $q=2$, we also write $\|\cdot\|_{S}$. We set $H^{m}(S,{\mathbb{R}}^{k})=W^{m,2}(S,{\mathbb{R}}^{k})$. We choose the numbers $p,q$ satisfying $1<p,q<\infty$ and $1/p+1/q=1$. For such $p$ and $q$ one can define the bilinear form on the product space $L^{p}(S,{\mathbb{R}}^{k})\times L^{q}(S,{\mathbb{R}}^{k})$ by

[TABLE]

For functions $v$ defined on $\Omega\times[0,\infty)$ we denote by $v(t)$ the mapping $x\mapsto v(x,t)$, which is defined on $\Omega$. The space $L^{q}(0,T_{e};X)$ denotes the Banach space of all Bochner-measurable functions $u:[0,T_{e})\to X$ such that $t\mapsto\|u(t)\|^{q}_{X}$ is integrable on $[0,T_{e})$. Finally, we frequently use the spaces $W^{m,q}(0,T_{e};X)$, which consist of Bochner measurable functions having $q$-integrable weak derivatives up to order $m$.

2 Preliminaries.

In this section we briefly recall some basic facts from convex analysis and nonlinear functional analysis which are needed for further discussions. For more details see [5, 15, 23, 33], for example.

Let $V$ be a reflexive Banach space with the norm $\|\cdot\|$, $V^{*}$ be its dual space with the norm $\|\cdot\|_{*}$. The brackets $\left<\cdot,\cdot\right>$ denote the duality pairing between $V$ and $V^{*}$. By $V$ we shall always mean a reflexive Banach space throughout this section.

For a function $\phi:V\to\overline{\mathbb{R}}$ the sets

[TABLE]

are called the effective domain and the epigraph of $\phi$, respectively. One says that the function $\phi$ is proper if $\operatorname{dom}(\phi)\not=\emptyset$ and $\phi(v)>-\infty$ for every $v\in V$. The epigraph is a non-empty closed convex set iff $\phi$ is a proper lower semi-continuous convex function or, equivalently, iff $\phi$ is a proper weakly lower semi-continuous convex function (see [33, Theorem 2.2.1]).

The Legendre-Fenchel conjugate of a proper convex lower semi-continuous function $\phi:V\to\overline{\mathbb{R}}$ is the function $\phi^{*}$ defined for each $v^{*}\in V^{*}$ by

[TABLE]

The Legendre-Fenchel conjugate $\phi^{*}$ is convex, lower semi-continuous and proper on the dual space $V^{*}$. Moreover, the Young-Fenchel inequality holds

[TABLE]

and the inequality $\phi\leq\psi$ implies $\psi^{*}\leq\phi^{*}$ for any two proper convex lower semi-continuous functions $\psi,\phi:V\to\overline{\mathbb{R}}$ (see [33, Theorem 2.3.1]).

Due to Proposition II.2.5 in [5] a proper convex lower semi-continuous function $\phi$ satisfies the following identity

[TABLE]

where $\partial\phi:V\to 2^{V^{*}}$ denotes the subdifferential of the function $\phi$. We note that the equality in (7) holds iff $v^{*}\in\partial\phi(v)$.

Remark 2.1**.**

We recall that the subdifferential of a lower semi-continuous proper and convex function is maximal monotone (see [5, Theorem II.2.1]) in the sense of Definition 2.1 below.

Convex integrands.

Let the numbers $p,q$ satisfy $1<q\leq 2\leq p<\infty,1/p+1/q=1.$ For a proper convex lower semi-continuous function $\phi:{\mathbb{R}}^{k}\to\overline{\mathbb{R}}$ we define a functional $I_{\phi}$ on $L^{p}(G,{\mathbb{R}}^{k})$ by

[TABLE]

where $G$ is a bounded domain in ${\mathbb{R}}^{N}$ with some $N\in{\mathbb{N}}$. Due to Proposition II.8.1 in [28], the functional $I_{\phi}$ is proper, convex, lower semi-continuous, and $v^{*}\in\partial I_{\phi}(v)$ iff

[TABLE]

Due to the result of Rockafellar in [24, Theorem 2], the Legendre-Fenchel conjugate of $I_{\phi}$ is equal to $I_{\phi^{*}}$, i.e.

[TABLE]

where $\phi^{*}$ is the Legendre-Fenchel conjugate of $\phi$.

Maximal monotone operators.

For a multivalued mapping $A:V\to 2^{V^{*}}$ the sets

[TABLE]

are called the effective domain and the graph of $A$, respectively.

Definition 2.1**.**

A mapping $A:V\to 2^{V^{*}}$ is called monotone if and only if the following inequality holds

[TABLE]

A monotone mapping $A:V\to 2^{V^{*}}$ is called maximal monotone iff the inequality

[TABLE]

implies $[v,v^{*}]\in GrA$.

It is well known ([23, p. 105]) that if $A$ is a maximal monotone operator, then for any $v\in D(A)$ the image $Av$ is a closed convex subset of $V^{*}$ and the graph $GrA$ is demi-closed111A set $A\in V\times V^{*}$ is demi-closed if $v_{n}$ converges strongly to $v_{0}$ in $V$ and $v^{*}_{n}$ converges weakly to $v^{*}_{0}$ in $V^{*}$ (or $v_{n}$ converges weakly to $v_{0}$ in $V$ and $v^{*}_{n}$ converges strongly to $v^{*}_{0}$ in $V^{*}$) and $[v_{n},v_{n}^{*}]\in GrA$, then $[v,v^{*}]\in GrA$.

Canonical extensions of maximal monotone operators.

In this subsection we briefly present some facts about measurable multi-valued mappings (see [4, 6, 15, 21], for example). We assume that $V$, and hence $V^{*}$, is separable and denote the set of maximal monotone operators from $V$ to $V^{*}$ by ${\mathfrak{M}}(V\times V^{*})$. Further, let $(S,\Sigma(S),\mu)$ be a $\sigma-$finite $\mu-$complete measurable space.

Definition 2.2**.**

A mapping $A:S\to{\mathfrak{M}}(V\times V^{*})$ is measurable iff for every open set $U\in V\times V^{*}$ (respectively closed set, Borel set, open ball, closed ball),

[TABLE]

is measurable in $S$.

The fact that the closed or Borel sets can be equivalently used in Definition 2.2 follows from the closedness of the values of the mapping $A:S\to{\mathfrak{M}}(V\times V^{*})$ (see [4, Theorem 8.1.4]).

Remark 2.2**.**

Theorem 8.1.4 in [4] also implies that under the above conditions the measurability of a mapping $A:S\to{\mathfrak{M}}(V\times V^{*})$ is equivalent to the existence of a countable dense subset consisting of measurable selectors, i.e. there exists a sequence of measurable functions $\{v_{n}\}_{n\in\mathbb{N}}:\,S\to V\times V^{*}$ such that for any $x\in S$ the image $A(x)$ can be represented as follows

[TABLE]

The following lemma will be used in the sequel (see [32, Lemma 3.1]).

Lemma 2.1**.**

Let a mapping $A:S\to{\mathfrak{M}}(V\times V^{*})$ be measurable. For any ${\cal L}(S)$-measurable function $v:S\to V$, the multivalued mapping $\hat{A}:x\mapsto A(x,v(x))$ is then closed-valued and measurable.

Given a mapping $A:S\to{\mathfrak{M}}(V\times V^{*})$, one can define a monotone graph from $L^{p}(S,V)$ to $L^{q}(S,V^{*})$, where $1/p+1/q=1$, as follows:

Definition 2.3**.**

Let $A:S\to{\mathfrak{M}}(V\times V^{*})$. The canonical extension of $A$ from $L^{p}(S,V)$ to $L^{q}(S,V^{*})$, where $1/p+1/q=1$, is defined by:

[TABLE]

In the following, we will drop the index $p$ for readability. Since we always work fix $p$ at the beginning of a statement, there cannot occur confusion with this notation. Monotonicity of ${\cal A}$ defined in Definition 2.3 is obvious, while its maximality follows from the next proposition (see [8, Proposition 2.13]).

Proposition 2.1**.**

Let $A:S\to{\mathfrak{M}}(V\times V^{*})$ be measurable. If $Gr{\cal A}\not=\emptyset$, then ${\cal A}$ is maximal monotone.

Remark 2.3**.**

We point out that the maximality of $A(x)$ for almost every $x\in S$ does not imply the maximality of ${\cal A}$ as the latter can be empty (see [8]).

Fitzpatrick’s function.

For a proper operator $\beta:V\to 2^{V^{*}}$ the Fitzpatrick function is defined as the convex and lower semicontinuous function given by

[TABLE]

It is known ([10]) that, whenever $\beta$ is maximal monotone,

[TABLE]

Any measurable maximal monotone operator $A:S\to{\mathfrak{M}}(V\times V^{*})$ can be represented by its Fitzpatrick function $f_{A}:S\times V\times V^{*}\to\overline{\mathbb{R}}$, which is $\Sigma(S)\otimes{\cal B}(V\times V^{*})$-measurable. Namely, the graph of a mapping $A:S\to{\mathfrak{M}}(V\times V^{*})$ can be written in the form (see [32, Proposition 3.2])

[TABLE]

We note that the measurability of the Fitzpatrick function $f_{A}:S\times V\times V^{*}\to\overline{\mathbb{R}}$ follows directly from its definition and Remark 2.2.

The graph of the canonical extension of a measurable operator $A:S\to{\mathfrak{M}}(V\times V^{*})$ can be equivalently represented in terms of its Fitzpatrick function $F_{{\cal A}_{p}}:L^{p}(S,V)\times L^{q}(S,V^{*})\to\overline{\mathbb{R}}$, i.e.

[TABLE]

Again, we omit $p$ if no confusion occurs. Moreover, the following result holds (see [32, Proposition 3.3])

• •

the functional $F_{\cal A}$ is convex and lower semi-continuous;

• •

for any $[v,v^{*}]\in L^{p}(S,V)\times L^{q}(S,V^{*})$, the integral

[TABLE]

exists either finite or equal to $+\infty$;

• •

if there exists a pair $[v,v^{*}]\in L^{p}(S,V)\times L^{q}(S,V^{*})$ such that $F_{\cal A}(v,v^{*})<+\infty$, then

[TABLE]

holds for all $[v^{*},v]\in L^{q}(S,V^{*})\times L^{p}(S,V)$.

3 Existence of solutions

In this section we introduce and show the existence of weak solutions for the initial boundary value (1) - (5). To simplify the notations, throughout the whole section we ignore the fact the coefficients and the given functions in (1) - (5) depend on $\omega\in\Omega$. The results proved below hold for a.e. $\omega\in\Omega$.

Solvability concept.

We start this section with the presentation of the intuitive ideas which lead to the definition of weak solutions for the initial boundary value problem (1) - (5). To give a meaning for the solvability of problem (1) - (5) we are going to use the concept of Fitzpatrick functions defined in (9).

We assume first that a triple of functions $(u_{\eta},\sigma_{\eta},z_{\eta})$ is given with the following properties: for every $t\in(0,T_{e})$ the function $(u_{\eta}(t),\sigma_{\eta}(t))$ is a weak solution of the boundary value problem

[TABLE]

This particularly holds for $z_{\eta}(0)=z_{\eta}^{(0)}$ and the corresponding initial values $(u_{\eta}(0),\sigma_{\eta}(0))=(u_{\eta}^{(0)},\sigma_{\eta}^{(0)})$. The equations (3) - (5) are satisfied pointwise for almost every $(x,t)$, and $b$ as well as $(u_{\eta},\sigma_{\eta},z_{\eta})$ are smooth enough. Then, based on equivalence (11), we can rewrite equation (3) as follows

[TABLE]

[TABLE]

which holds for almost every $(x,t)\in{\boldsymbol{Q}}\times(0,T_{e})$. Integrating the last equality over ${\boldsymbol{Q}}$ gives

[TABLE]

Using (1), (2) and (4) the right hand side in (15) becomes (${\mathbb{A}}_{\eta}:={\mathbb{C}}_{\eta}^{-1}$)

[TABLE]

Integrating relations (15) and (3) with respect to $t$ leads to

[TABLE]

Taking into account the inequality (10), we conclude that the triple of functions $(u_{\eta},\sigma_{\eta},z_{\eta})$ satisfies equality (3) if and only if the inequality

[TABLE]

holds for all $t\in(0,T_{e})$ and some function $\sigma^{(0)}_{\eta}\in L^{2}({\boldsymbol{Q}},{\cal S}^{3})$ solving the elliptic boundary value problem (12) - (14).

The above computations suggest the following notion of weak solutions for the initial boundary value problem (1) - (5).

Definition 3.1**.**

Let the numbers $p,q$ satisfy $1<q\leq 2\leq p<\infty,1/p+1/q=1.$ A function $(u_{\eta},\sigma_{\eta},z_{\eta})$ such that

[TABLE]

[TABLE]

with

[TABLE]

is called a weak solution of the initial boundary value problem (1) - (5), if for every $t\in(0,T_{e})$ the function $(u_{\eta}(t),\sigma_{\eta}(t))$ is a weak solution of the boundary value problem (1) - (2), (4) for every given $Bz_{\eta}(t)\in L^{q}({\boldsymbol{Q}},{\cal S}^{3})$, the initial condition (5) is satisfied pointwise for almost every $(x,t)$ and the inequality (3) holds for all $t\in(0,T_{e})$ and the function $\sigma_{\eta}^{(0)}\in L^{2}({\boldsymbol{Q}},{\cal S}^{3})$ determined by equations (12) - (14).

Now, we show that the above definition of weak solutions for (1) - (5) is consistent. Namely, we are going to prove that if a triple of functions $(u_{\eta},\sigma_{\eta},z_{\eta})$ is a weak solution of (1) - (5) in the sense of Definition 3.1 and possesses additional regularity, then this triple of functions is a solution of the initial boundary value problem (1) - (5), i.e. the constitutive inclusion (3) is satisfied pointwise for a.e. $(x,t)\in{\boldsymbol{Q}}_{T_{e}}$. To this end, we assume that the weak solution $(u_{\eta},\sigma_{\eta},z_{\eta})$ has the following regularity

[TABLE]

[TABLE]

Then, it is immediately seen that the function $\sigma_{\eta}^{(0)}\in L^{2}({\boldsymbol{Q}},{\cal S}^{3})$ as a unique solution of the problem (12) - (14) satisfies the relation $\sigma^{(0)}_{\eta}(x)=\sigma_{\eta}(x,0)$ for a.e. $x\in{\boldsymbol{Q}}$ and the following identity

[TABLE]

Moreover, we have that

[TABLE]

Then, the inequality (3) can be rewritten as follows

[TABLE]

Handling the equations (1) - (2) as above we obtain that the last inequality takes the following form

[TABLE]

or, equivalently,

[TABLE]

Therefore, by (10) and the standard localization argument we get that

[TABLE]

[TABLE]

which holds for a.e. $(x,t)\in{\boldsymbol{Q}}_{T_{e}}$. Now, based on the equivalence result (11) we conclude that the inclusion (3) is satisfied pointwise from the assumed temporal regularity of $(u_{\eta},\sigma_{\eta},z_{\eta})$. The pointwise meaning of (5) follows.

Existence result.

First, we define a class of maximal monotone functions we deal with in this work.

Definition 3.2**.**

Let $S$ be a measurable set in ${\mathbb{R}}^{s}$ and $m\in L^{1}(S,\mathbb{R})$. For $\alpha_{1},\alpha_{2}\in{\mathbb{R}}_{+}$, ${\cal M}(S,{\mathbb{R}}^{k},\alpha_{1},\alpha_{2},m)$ is the set of measurable multi-valued functions $h:S\to{\mathfrak{M}}({\mathbb{R}}^{k}\times{\mathbb{R}}^{k})$ (in the sense of Definition 2.2) such that with the following inequality

[TABLE]

holds for a.e. $x\in S$ and every $v^{*}\in h(x,v)$, where $p$ and $q$ satisfy the relations $2\leq p<\infty$ and $q=p/(p-1)$.

The main properties of the class ${\cal M}(S,{\mathbb{R}}^{k},\alpha_{1},\alpha_{2},m)$ are collected in the following proposition (see [8, Corollary 2.15]).

Proposition 3.1**.**

Let ${\cal H}$ be a canonical extension of a function $h:S\to{\mathfrak{M}}({\mathbb{R}}^{k}\times{\mathbb{R}}^{k})$ in the sense of Definition 2.3, which belongs to ${\cal M}(S,{\mathbb{R}}^{k},\alpha_{1},\alpha_{2},m)$. Then ${\cal H}$ is maximal monotone, surjective and $D({\cal H})=L^{p}(S,{\mathbb{R}}^{k})$.

Now, we can state the main result of this section.

Theorem 3.1**.**

Assume that $L_{\eta}$ is positive semi-definite, ${\mathbb{C}}_{\eta}$ is uniformly positive definite and ${\mathbb{C}}_{\eta}\in C(\bar{\boldsymbol{Q}},{\cal L}({\cal S}^{3},{\cal S}^{3}))$, the mappings $g_{\eta}\in{\cal M}({\boldsymbol{Q}},{\mathbb{R}}^{N},\alpha_{1},\alpha_{2},m)$ with a function $m$ from $L^{1}({\boldsymbol{Q}},{\mathbb{R}})$. Suppose that $b\in W^{1,p}(0,T_{e};W^{-1,p}({\boldsymbol{Q}},{\mathbb{R}}^{3}))$ and $z_{\eta}^{(0)}\in L^{2}({\boldsymbol{Q}},{\mathbb{R}}^{N})$.

*Then the initial boundary value problem (1) - (5) has at least one weak solution $(u_{\eta},T_{\eta},z_{\eta})$ in the sense of Definition 3.1. *

Remark 3.1**.**

We point out that the requirement of the continuity of ${\mathbb{C}}_{\eta}$ is superfluous and is only made to simplify the proof of Theorem 3.1. The proof itself works for the measurable function ${\mathbb{C}}_{\eta}$ as well. The continuity assumption allows us to apply the $L^{p}$-regularity theory for linear elliptic systems in [12] directly to our problem. In case of ${\mathbb{C}}_{\eta}\in L^{\infty}({\boldsymbol{Q}},{\cal L}({\cal S}^{3},{\cal S}^{3}))$, some extra technical work has to be done before one can use the $L^{p}$-regularity theory for linear elliptic systems (this strategy is realized in [20]). To avoid the technicalities we assume the continuity of ${\mathbb{C}}_{\eta}$ here.

Proof.

To simplify the notations we drop $\eta$. The proof of the theorem is presented in [19]. Therefore, we only sketch it here. We show this by the Rothe method (a time-discretization method, see [25] for details). In order to introduce a time-discretized problem, let us fix any $m\in{\mathbb{N}}$ and set

[TABLE]

We are looking for functions $u^{n}_{m}\in H^{1}_{0}({\boldsymbol{Q}},{\mathbb{R}}^{3})$, $\sigma^{n}_{m}\in L^{2}({\boldsymbol{Q}},{\cal S}^{3})$ and $z^{n}_{m}\in L^{2}({\boldsymbol{Q}},{\mathbb{R}}^{N})$ with

[TABLE]

solving the following problem

[TABLE]

together with the boundary conditions

[TABLE]

The proof of the existence of the triple

[TABLE]

satisfying (20) - (23) can be found in [19].

A-priori estimates. Multiplying (20) by $(u^{n}_{m}-u^{n-1}_{m})/h$ and then integrating over ${\boldsymbol{Q}}$ we get

[TABLE]

Applying $g^{-1}(x)$ to both sides of (22), multiplying by $w^{n}_{m}:=(z^{n}_{m}-z^{n-1}_{m})/h$ and then integrate over ${\boldsymbol{Q}}$ to obtain

[TABLE]

With (24) we get that

[TABLE]

[TABLE]

Multiplying by $h$ and summing the obtained relation for $n=1,...,l$ for any fixed $l\in[1,2^{m}]$ we derive the following inequality (${\mathbb{A}}={\mathbb{C}}^{-1}$)

[TABLE]

where

[TABLE]

We estimate now the right hand side of the last inequality. Since $u^{n}_{m}$ is a solution of the linear elliptic problem formed by the equations (20), (21) and (23), it satisfies (see [12]) the inequality

[TABLE]

where $C$ is a positive constant independent of $n$ and $m$. Therefore, using the linearity of the problem formed by (20), (21) and (23), the inequality (26) and Young’s inequality with $\epsilon>0$ we get that

[TABLE]

where $C_{\epsilon}$ is a positive constant appearing in the Young inequality. Combining the inequalities (25) and (3), applying (10) and (19) and choosing an appropriate value for $\epsilon>0$ we obtain the following estimate

[TABLE]

where $\tilde{C},\tilde{C}_{\epsilon}$ and $\hat{C}_{\epsilon}$ are some positive constants. Now, using the definition of Rothe’s approximation functions (see (68)) we rewrite (28) as follows

[TABLE]

From the estimate (29) we get then that

[TABLE]

In particular, the uniform boundness of the sequences in (30) - (33) yields

[TABLE]

Employing (69), the estimates (31) - (34) further imply that the sequnces $\{\sigma_{m}\}_{m}$, $\{L^{1/2}z_{m}\}_{m}$, $\left\{z_{m}/\sqrt{m}\right\}_{m}$ and $\left\{\Sigma_{m}\right\}_{m}$ are also uniformly bounded in the spaces $L^{\infty}(0,{T_{e}};L^{2}({\boldsymbol{Q}},{\cal S}^{3}))$, $L^{\infty}(0,{T_{e}};L^{2}({\boldsymbol{Q}},{\mathbb{R}}^{N}))$, $L^{\infty}(0,{T_{e}};L^{2}({\boldsymbol{Q}},{\mathbb{R}}^{N}))$ and $L^{p}(0,{T_{e}};L^{p}({\boldsymbol{Q}},{\mathbb{R}}^{N}))$, respectively. Moreover, due to (30) and the following obvious relation

[TABLE]

we may conclude that $\{\bar{z}_{m}\}_{m}$ is uniformly bounded in $L^{q}(0,{T_{e}};L^{q}({\boldsymbol{Q}},{\mathbb{R}}^{N}))$.

In [19] it is shown that the limit functions denoted by $u,T,z$ and $\Sigma$ of the corresponding weakly convergent sequences have the following properties

[TABLE]

and

[TABLE]

To prove that the weak limit of $(u_{m},T_{m},z_{m})$ is a weak solution of the problem (1) - (5), we are going to employ the concept of the Fitzpatrick function again. To this end, we rewrite (25) as follows

[TABLE]

Next, using the lower semi-continuity of convex functionals we get (3) after passing to the weak limit in (36). This completes the proof of Theorem 3.1. ∎

4 Stochastic homogenization

Throughout this section, we follow the setting for stochastic homogenization proposed in [14] for rate-independent systems.

Remark 4.1**.**

In the following, we introduce the concept of Palm measures. Note that we will need this concept only in the context of the results in Section 5. For the main results proved in Section 6 we will restrict to the case $\mu_{\omega}^{\eta}={\cal L}$ which implies $\mu_{{\cal P}}={\cal P}$ (this follows from the translation invariance and Fubini’s theorem). In this case, we will omit $\mu_{{\cal P}}$ and every integral over $\Omega$ is meant with respect to ${\cal P}$. In particular, we will write $\int_{\Omega}f:=\int_{\Omega}f(\omega)d{\cal P}(\omega)$.

4.1 Concept of Palm measures

Let $(\Omega,{\cal F}_{\Omega},{\cal P},\tau)$ be a probability space with dynamical system satisfying Assumption 1.1 and let ${\cal M}(\mathbb{R}^{n})$ be the set of Radon measures on $\mathbb{R}^{n}$ equipped with the Vague topology.

Definition 4.1**.**

Let $(\Omega,{\cal F}_{\Omega},{\cal P},\tau)$ satisfy Assumtion 1.1. A random measure is a mapping $\mu_{\bullet}:\,\Omega\to{\cal M}(\mathbb{R}^{n})$, $\omega\mapsto\mu_{\omega}$ such that $\omega\mapsto\mu_{\omega}(A)$ is measurable for all Borel sets $A\subset\mathbb{R}^{n}$. A random measure is called stationary, if $\mu_{\tau_{x}\omega}(A)=\mu_{\omega}(A+x)$ for all Borel sets $A\subset\mathbb{R}^{n}$. The intensity $\lambda(\mu_{\omega})$ is defined by:

[TABLE]

Theorem 4.1** (Mecke [16, 7]: Existence of Palm measure).**

Let $\omega\mapsto\mu_{\omega}$ be a stationary random measure. Then there exists a unique measure $\mu_{{\cal P}}$ on $\Omega$ such that

[TABLE]

for all ${\cal L}\times\mu_{{\cal P}}$-measurable non negative functions and all ${\cal L}\times\mu_{{\cal P}}$- integrable functions $f$. Furthermore for all $A\subset\Omega$, $u\in L^{1}(\Omega,\mu_{{\cal P}})$ there holds

[TABLE]

for an arbitrary $g\in L^{1}(\mathbb{R}^{n},{\cal L})$ with $\int_{\mathbb{R}^{n}}g(x)dx=1$ and $\mu_{{\cal P}}$ is $\sigma$-finite.

Remark 4.2**.**

a) Setting $g(s):=\chi_{[0,1]^{n}}(s)$, the Palm measure can equally be defined through (38).

b) For the constant measure $\omega\mapsto{\cal L}$, we simply find $\mu_{{\cal P}}={\cal P}$, the original probability measure. This is a direct consequence of (38), Fubini’s theorem and Assumption 1.1 (ii).

For a random measure $\mu_{\omega}$, we define

[TABLE]

Theorem 4.2** (Ergodic Theorem [7]).**

Let Assumption 1.1 hold for $(\Omega,{\cal F}_{\Omega},{\cal P},\tau)$. Let $\mu_{\omega}$ be a stationary random measure with finite intensity and Palm measure $\mu_{{\cal P}}$. Then, for all $g\in L^{1}(\Omega,\mu_{{\cal P}})$ there holds ${\cal P}$ almost surely

[TABLE]

for all bounded Borel sets $A$.

The ergodic theorem only holds for functions on $\Omega$. Nevertheless, it motivates the following generalization of the concept of ergodicity:

Definition 4.2**.**

Let $f\in L^{p}(\boldsymbol{Q}\times\Omega;{\cal L}\otimes\mu_{{\cal P}})$ for some $1\leq p<\infty$. We say that $f$ is an ergodic function if it has a ${\cal B}(\boldsymbol{Q})\otimes{\cal F}_{\Omega}$-measurable representative $\tilde{f}$ such that for $a=\tilde{f}$ and $a=|\tilde{f}|^{p}$ it holds

[TABLE]

Finally, we say $f\in L^{p}(\boldsymbol{Q}\times\Omega;{\cal L}\otimes\mu_{{\cal P}})$ is an ergodic function if (42) holds for almost all $\omega$.

By Theorem 4.2, we find that every $g\in L^{1}(\Omega,\mu_{{\cal P}})$ is ergodic. In [14, Section 2.5] a larger set of ergodic functions was identified:

Lemma 4.1**.**

Let Assumption 1.1 hold for $(\Omega,{\cal F}_{\Omega},{\cal P},\tau)$. Let $\boldsymbol{Q}\subset\mathbb{R}^{n}$ be a bounded domain and let $f\in L^{\infty}(\boldsymbol{Q}\times\Omega;{\cal L}\otimes\mu_{{\cal P}})$. Then, $f$ is an ergodic function.

4.2 Potentials and Solenoidals

Let ${\cal BD}(\mathbb{R}^{n})$ denote the set of bounded domains in $\mathbb{R}^{n}$. For every $p$ with $1<p<\infty$, we introduce the following spaces:

[TABLE]

On $\Omega$, we introduce the corresponding spaces

[TABLE]

Then, there holds the following orthogonal decomposition.

Lemma 4.2** ([21, Theorem 3.1.2]).**

Let $1<p<\infty$ and $p^{-1}+q^{-1}=1$ and let Assumption 1.1 hold for $(\Omega,{\cal F}_{\Omega},{\cal P},\tau)$. Then the following relations hold in the sense of duality between the spaces $L^{p}(\Omega,{\cal P})$ and $L^{q}(\Omega,{\cal P})$:

[TABLE]

Every $L_{\rm pot}^{p}(\Omega)$ function can be optained as the ergodic limit of a sequence of gradients with vanishing potentials. The following result can be proved like in [13, Section 2.3].

Lemma 4.3**.**

Let $\upsilon\in L_{\rm pot}^{p}(\Omega)$, $1<p<\infty$. For almost every $\omega$ there exists $C>0$ such that the following holds: For every $\eta>0$ there exists a unique $V_{\eta}^{\omega}$ with $\nabla V_{\eta}^{\omega}(x)=\upsilon(\tau_{\frac{x}{\eta}}\omega)$ and $\|V_{\eta}\|_{H^{1}(\boldsymbol{Q})}\leq C\|\upsilon\|_{L_{\rm pot}^{p}(\Omega)}$ for all $\eta>0$. Furthermore,

[TABLE]

Furthermore, we find the following important Korn inequality, which can be proved like in [13, Section 2.3].

Lemma 4.4**.**

Let Assumption 1.1 hold for $(\Omega,{\cal F}_{\Omega},{\cal P},\tau)$. For every $1<p<\infty$ there exists $C_{p}>0$ such that for all $\upsilon\in L_{\rm pot}^{p}(\Omega;\mathbb{R}^{n})$

[TABLE]

4.3 Two-scale convergence: time independent case

Let Assumption 1.1 hold for $(\Omega,{\cal F}_{\Omega},{\cal P},\tau)$ and let $\omega\mapsto\mu_{\omega}$ be a stationary random measure with $\mu_{\omega}^{\eta}$ and $\mu_{{\cal P}}$ defined through (40) and (38). For the case $\mu_{\omega}={\cal L}$ we recall Remark 4.1.

Remark 4.3**.**

The product $\sigma$-algebra ${\cal B}_{Q}\otimes{\cal F}_{\Omega}$ is countably generated and therefore, the space $L^{p}(\boldsymbol{Q}\times\Omega)$ is separable (see [9, Section VI.15, p. 92]) for every $1\leq p<\infty$. In particular, for every $1\leq p<\infty$, there exists a countable dense subset $\Phi_{p}$ of simple functions in $L^{p}(\boldsymbol{Q}\times\Omega;{\cal L}\otimes\mu_{{\cal P}})$. By Lemma 4.1, every $\phi\in\Phi_{p}$ is an ergodic function and there exists a set $\Omega_{\Phi_{p}}\subset\Omega$ with ${\cal P}(\Omega_{\Phi_{p}})=1$ such that all $\phi\in\Phi_{p}$ satisfy (42) (i.e. they admit ergodic realizations) for all $\omega\in\Omega_{\Phi_{p}}$. This corresponds to the setting of [14].

Definition 4.3**.**

Let $1<p,q<\infty$ with $1/p+1/q=1$. Let $\Phi_{p}$ be the set of Remark 4.3 and let $\omega\in\Omega_{\Phi_{p}}$. Let $u^{\eta}\in L^{q}(\boldsymbol{Q};\mu_{\omega}^{\eta})$ for all $\eta>0$. We say that $u^{\eta}$ converges (weakly) in two scales to $u\in L^{q}(\boldsymbol{Q};L^{q}(\Omega,\mu_{{\cal P}}))$ and write $u^{\eta}\stackrel{{\scriptstyle 2s}}{{\rightharpoonup}}u$ if $\sup_{\eta}\|u^{\eta}\|_{L^{q}(\boldsymbol{Q};\mu_{\omega}^{\eta})}<\infty$ and if for all $\phi\in\Phi_{p}$ there holds with $\phi_{\omega,\eta}(x):=\phi(x,\tau_{\frac{x}{\eta}}\omega)$

[TABLE]

Furthermore, we say that $u^{\eta}$ converges strongly in two scales to $u$, written $u^{\eta}\stackrel{{\scriptstyle 2s}}{{\to}}u$, if for all weakly two-scale converging sequences $v^{\eta}\in L^{p}(\boldsymbol{Q};\mu_{\omega}^{\eta})$ with $v^{\eta}\stackrel{{\scriptstyle 2s}}{{\rightharpoonup}}v\in L^{p}(\boldsymbol{Q};L^{p}(\Omega,\mu_{{\cal P}}))$ as $\eta\to 0$ there holds

[TABLE]

Lemma 4.5** (Existence of two-scale limits [14]).**

Let $\omega\in\Omega_{\Phi_{p}}$, $1<p<\infty$ and $u^{\eta}\in L^{p}(\boldsymbol{Q};\mu_{\omega}^{\eta})$ be a sequence of functions such that $\left\|u^{\eta}\right\|_{L^{p}(\boldsymbol{Q};\mu_{\omega}^{\eta})}\leq C$ for some $C>0$ independent of $\eta$. Then there exists a subsequence of $u^{\eta^{\prime}}$ and $u\in L^{p}(\boldsymbol{Q};L^{p}(\Omega,\mu_{{\cal P}}))$ such that $u^{\eta^{\prime}}\stackrel{{\scriptstyle 2s}}{{\rightharpoonup}}u$ and

[TABLE]

Closely connected with the definition of two-scale convergence and Lemma 4.5 is the following result.

Lemma 4.6**.**

Let $f\in L^{1}(\Omega,\mu_{{\cal P}})$. For almost all $\omega\in\Omega$ it holds: for all $\psi\in C_{0}(\boldsymbol{Q})$ and $u^{\eta}\in L^{q}(\boldsymbol{Q};\mu_{\omega}^{\eta})$, $\eta>0$, and $u\in L^{q}(\boldsymbol{Q};L^{q}(\Omega,\mu_{{\cal P}}))$ with $u^{\eta}\stackrel{{\scriptstyle 2s}}{{\rightharpoonup}}u$ it holds

[TABLE]

The following Lemma is well known in the periodic case ([3]) but also in the stochastic setting ([14, 34] for $p=2$). The following version can be proofed along the same lines as Lemma 6.2 in [14].

Lemma 4.7**.**

Let $1<p<\infty$. If $u^{\eta}\in W_{0}^{1,p}(\boldsymbol{Q};\mathbb{R}^{n})$ for all $\eta$ with $\|\nabla u^{\eta}\|_{L^{p}(\boldsymbol{Q})}<C$ for $C$ independent from $\eta>0$ then there exists a subsequence $u^{\eta^{\prime}}$ and functions $u\in W_{0}^{1,p}(\boldsymbol{Q};\mathbb{R}^{n})$ and $v\in L^{p}(\boldsymbol{Q};L_{pot}^{p}(\Omega;\mathbb{R}^{n}))$ such that

[TABLE]

We finally collect some usefull results.

Lemma 4.8**.**

Let $u\in L^{p}(\boldsymbol{Q};L^{p}(\Omega;\mu_{{\cal P}}))$. Then, for almost every $\omega\in\Omega$, there exists a sequence $u^{\eta}\in L^{p}(\boldsymbol{Q};\mu_{\omega}^{\eta})$ such that $u^{\eta}\stackrel{{\scriptstyle 2s}}{{\rightharpoonup}}u$ as $\eta\to 0$.

Lemma 4.9**.**

Let $N\in\mathbb{N}$ and let $A\in L^{\infty}(\boldsymbol{Q};L^{\infty}(\Omega;\mathscr{L}(\mathbb{R}^{N},\mathbb{R}^{N})))$ be symmetric and assume $A$ is ${\cal B}_{Q}\otimes{\cal F}_{\Omega}$ -measurable. We furthermore assume the existence of a constant $\alpha>0$ such that

[TABLE]

Then, for almost all $\omega\in\Omega_{\Phi_{p}}$ there holds: For all sequences $u^{\eta}\in L^{2}(\boldsymbol{Q};\mu_{\omega}^{\eta};\mathbb{R}^{N})$ with weak two-scale limit $u\in L^{2}(\boldsymbol{Q};L^{2}(\Omega;\mu_{{\cal P}};\mathbb{R}^{N}))$ there holds with $A_{\eta,\omega}(x):=A(x,\tau_{\frac{x}{\eta}}\omega)$

[TABLE]

4.4 Two-scale convergence: time dependent case

We are also interested in the convergence behavior of functions $u^{\eta}:\,[0,T]\to L^{p}(\boldsymbol{Q},\mu_{\omega}^{\eta})$.

Definition 4.4**.**

Let $1<r,r^{\prime},p,q<\infty$ with $\frac{1}{p}+\frac{1}{q}=1$ and $\frac{1}{r^{\prime}}+\frac{1}{r}=1$. Let $\Phi_{q}$ be the set of Remark 4.3 and let $\omega\in\Omega_{\Phi_{q}}$. Let $u^{\eta}\in L^{r}(0,T;L^{p}(\boldsymbol{Q};\mu_{\omega}^{\eta}))$ for all $\eta>0$. We say that $u^{\eta}$ converges (weakly) in two scales to $u\in L^{r}(0,T;L^{p}(\boldsymbol{Q};L^{p}(\Omega,\mu_{{\cal P}})))$ and write $u^{\eta}\stackrel{{\scriptstyle 2s}}{{\rightharpoonup}}u$ if for all continuous and piecewise affine functions $\phi:\,[0,T]\to\Phi_{q}$ there holds with $\phi_{\omega,\eta}(t,x):=\phi(t,x,\tau_{\frac{x}{\eta}}\omega)$

[TABLE]

The following two lemmas where proved in [14].

Lemma 4.10**.**

Asssume that $1<p<\infty$ and $1<r\leq\infty$. Then, every sequence of functions $\left(u^{\eta}\in L^{r}(0,T;L^{p}(\boldsymbol{Q};\mu_{\omega}^{\eta}))\right)_{\varepsilon>0}$ satisfying

[TABLE]

for some $C>0$ independent from $\eta$ has a weakly two-scale convergent subsequence with limit function $u\in L^{r}(0,T;L^{p}(\boldsymbol{Q};L^{p}(\Omega,\mu_{{\cal P}})))$. Furthermore, if

[TABLE]

uniformly for $1<p\leq\infty$, then also $\left\|\partial_{t}u\right\|_{L^{r}(0,T;L^{p}(\boldsymbol{Q};L^{p}(\Omega,\mu_{{\cal P}})))}\leq C$ and $\partial_{t}u^{\eta}\stackrel{{\scriptstyle 2s}}{{\rightharpoonup}}\partial_{t}u$ in the sense of Definition 4.4 as well as $u^{\eta}(t)\stackrel{{\scriptstyle 2s}}{{\rightharpoonup}}u(t)$ for all $t\in[0,T]$.

5 Homogenization of convex functionals

Lemma 5.1**.**

Let Assumption 1.1 hold for $(\Omega,{\cal F}_{\Omega},{\cal P},\tau)$ and the random measure $\mu_{\omega}$. Let $f:\,\boldsymbol{Q}\times\Omega\times\mathbb{R}^{N}\to\mathbb{R}$ be a convex functional in $\mathbb{R}^{N}$. For almost all $\omega\in\Omega_{\Phi_{p}}$ the following holds: Let $u^{\eta}\in L^{q}(\boldsymbol{Q};\mathbb{R}^{N})$ be a sequence such that $\left\|u^{\eta}\right\|_{L^{q}(\boldsymbol{Q})}\leq C$ for some $0<C<\infty$ and such that $u^{\eta}\stackrel{{\scriptstyle{\scriptstyle 2s}}}{{\rightharpoonup}}u\in L^{q}(\boldsymbol{Q}\times\Omega;{\cal L}\otimes\mu_{{\cal P}};\mathbb{R}^{N})$. Then, it holds

[TABLE]

The proof of Lemma 5.1 is literally the same as for Theorem 7.1 in [35]. However, we provide it here for completeness.

Proof.

Let $\omega\in\Omega_{\Phi_{p}}$ and let $f^{\ast}$ denote the Fenchel conjugate of $f$ in the third variable. Without loss of generality, we may assume that

[TABLE]

for all $\psi\in\Phi_{p}^{N}$ and all $\omega\in\Omega_{\Phi_{p}}$. We first consider the case

[TABLE]

for almost every $(x,\tilde{\omega})\in\boldsymbol{Q}\times\Omega$ and all $\xi\in\mathbb{R}^{N}$. We then find for every $\psi\in\Phi_{p}^{N}$

[TABLE]

Due to $u^{\eta}\stackrel{{\scriptstyle{\scriptstyle 2s}}}{{\rightharpoonup}}u$ and (46) we find

[TABLE]

for all $\psi\in\Phi_{p}^{N}$. Since (47) holds, $f^{\ast}$ is continuous in $\xi$ and the last inequality implies

[TABLE]

In the general case, let

[TABLE]

Then, $0<F_{\eta}^{\delta}-F_{\eta}\leq C\delta$ and (48) implies that

[TABLE]

Hence the claim follows.∎

Lemma 5.2**.**

Let Assumption 1.1 hold for $(\Omega,{\cal F}_{\Omega},{\cal P},\tau)$ and let $\mu_{\omega}$ be a random measure. Let $f:\,\boldsymbol{Q}\times\Omega\times\mathbb{R}^{N}\to\mathbb{R}$ be such that for a.e. $(x,\tilde{\omega})$ the function $f(x,\tilde{\omega},\cdot)$ is convex in $\mathbb{R}^{N}$. Then, for almost every $\omega\in\Omega_{\Phi_{p}}$ the following holds: If $u^{\eta}_{\omega}\in L^{q}(\boldsymbol{Q};\mathbb{R}^{N})$ is a sequence of minimizers of the functionals

[TABLE]

and if $\sup_{\eta>0}\left\|u_{\omega}^{\eta}\right\|_{L^{q}(\boldsymbol{Q};\mu_{\omega}^{\eta})}<\infty$, then there exists $u_{0,\omega}\in L^{q}(\boldsymbol{Q}\times\Omega;{\cal L}\otimes\mu_{{\cal P}};\mathbb{R}^{N})$ such that $u^{\eta}_{\omega}\stackrel{{\scriptstyle{\scriptstyle 2s}}}{{\rightharpoonup}}u_{0,\omega}$ along a subsequence and $u_{0,\omega}$ is a minimizer of

[TABLE]

Proof.

Let $u_{0}\in L^{q}(\boldsymbol{Q}\times\Omega;{\cal L}\otimes\mu_{{\cal P}})$ be a minimizer of $F_{0}$. By [29, Theorem III-39] we can assume that $u_{0}(x,\tilde{\omega})$ minimizes $f(x,\tilde{\omega})$ for almost every $(x,\tilde{\omega})$. Then, for almost all $\omega\in\Omega$ it holds $u^{\eta}_{0}(x):=u_{0}(x,\tau_{\frac{x}{\eta}}\omega)\in L^{q}(\boldsymbol{Q};\mathbb{R}^{N})$ and

[TABLE]

We chose a subsequence $u_{\omega}^{\eta^{\prime}}$ and $u_{0,\omega}\in L^{q}(\boldsymbol{Q}\times\Omega;{\cal L}\otimes\mu_{{\cal P}};\mathbb{R}^{N})$ such that $u_{\omega}^{\eta^{\prime}}\stackrel{{\scriptstyle{\scriptstyle 2s}}}{{\rightharpoonup}}u_{0,\omega}$. Since $F_{\eta,\omega}(u^{\eta}_{0})\to F_{0}(u_{0})$, we find

[TABLE]

Hence, $u_{0,\omega}$ is a minimizer of $F_{0}$. ∎

Theorem 5.1**.**

Let Assumption 1.1 hold for $(\Omega,{\cal F}_{\Omega},{\cal P},\tau)$ and let $\mu_{\omega}$ be a random measure and let $1<p,q<\infty$ with $\frac{1}{p}+\frac{1}{q}=1$. Let $\boldsymbol{Q}\subset\mathbb{R}^{n}$ be a bounded domain and $f:\,\boldsymbol{Q}\times\Omega\times\mathbb{R}^{d}\to\mathbb{R}$ be measurable and for all $(x,\omega)$. Let $f(x,\omega,\cdot)$ be continuous and convex in $\mathbb{R}^{d}$ with $f(x,\omega,\xi)\geq|\xi|^{q}$. Then, for almost every $\omega\in\Omega_{\Phi_{p}}$ it holds: If $u^{\eta}\in W^{1,q}(\boldsymbol{Q})$ is a sequence of minimizers of the functional

[TABLE]

such that $\sup_{\eta>0}\left\|u_{\omega}^{\eta}\right\|_{L^{q}(\boldsymbol{Q};\mu_{\omega}^{\eta})}<\infty$, then there exist $u_{\omega}\in W_{0}^{1,q}(\boldsymbol{Q})$ and $\upsilon_{\omega}\in L^{q}(\boldsymbol{Q};L_{\rm pot}^{q}(\Omega))$ and a subsequence $u^{\eta^{\prime}}$ such that $u^{\eta^{\prime}}\to u_{\omega}$ strongly in $L^{q}(\boldsymbol{Q})$ and $\nabla u^{\eta^{\prime}}\stackrel{{\scriptstyle{\scriptstyle 2s}}}{{\rightharpoonup}}\nabla u_{\omega}+\upsilon_{\omega}$ as $\eta\to 0$ and $(u_{\omega},\upsilon_{\omega})$ is a minimizer of the functional

[TABLE]

Proof.

Let $\Phi_{\rm pot}$ be a countable dense subset of $L^{q}_{\rm pot}(\Omega)$ and let $\tilde{\Omega}\subset\Omega_{\Phi_{q}}$ be a set of full measure such that Lemma lem:valid-ts-test-function holds for all $\omega\in\tilde{\Omega}$. By ${\rm span}\Phi_{\rm pot}$ we denote finite linear combinations of elements of $\Phi_{\rm pot}$. In what follows we restrict to the case $\omega\in\tilde{\Omega}$.

Due to Lemma 4.7 there exist $u_{\omega}\in W^{1,q}(\boldsymbol{Q})$ and $\upsilon_{\omega}\in L^{q}(\boldsymbol{Q};L_{\rm pot}^{q}(\Omega))$ such that $\nabla u^{\eta}\stackrel{{\scriptstyle{\scriptstyle 2s}}}{{\rightharpoonup}}\nabla u_{\omega}+\upsilon_{\omega}$ and $u^{\eta}\to u_{\omega}$ along a subsequence, which we denote $u^{\eta}$ for simplicity. Let $u_{0}\in W_{0}^{1,q}(\boldsymbol{Q})$ and $\upsilon_{0}\in L^{q}(\boldsymbol{Q};L_{\rm pot}^{q}(\Omega))$ be a minimizer of the functional $F_{0}$.

Now, let $\delta>0$. There exists $\upsilon_{\delta}\in L^{q}(\boldsymbol{Q};\mathbb{R}\Phi_{\rm pot})$ which is simple and has compact support in $\boldsymbol{Q}$ such that $\|\upsilon_{0}-\upsilon_{\delta}\|_{L^{q}(\boldsymbol{Q};L_{\rm pot}^{q}(\Omega))}<\delta$. In particular, we find sets $A_{i}\subset{\boldsymbol{Q}}$, $1\leq i\leq K_{\delta}$ and functions $\hat{\upsilon}_{i}\in{\rm span}\Phi_{\rm pot}$ such that

[TABLE]

Let $\left(\varphi_{\varepsilon}\right)_{\varepsilon>0}\subset C^{\infty}_{0}({\rm B}_{\varepsilon})$ be a family of mollifiers. For $\varepsilon>0$ we denote $\upsilon_{\varepsilon,\delta}(\cdot,\omega):=\varphi_{\varepsilon}\ast\upsilon(\cdot,\omega)$, where $\ast$ is the convolution with respect to the $\boldsymbol{Q}$-coordinate. Then $\upsilon_{\varepsilon,\delta}\in C^{1}_{0}(\overline{\boldsymbol{Q}};{\rm span}\Phi_{\rm pot})$ for $\varepsilon>0$ small enough and $\|\upsilon_{\varepsilon,\delta}-\upsilon_{\delta}\|_{L^{q}(\boldsymbol{Q};L_{\rm pot}^{q}(\Omega))}\to 0$ as $\varepsilon\to 0$.

Given $x\in\boldsymbol{Q}$ we apply Lemma 4.3 and denote $V_{\eta,\varepsilon,\delta}^{\omega}(x,\cdot)\in H^{1}(\boldsymbol{Q})$ the $\eta$-potential to $\upsilon_{\varepsilon,\delta}(x)$ and $\eta$ and $V_{\eta,\delta}^{\omega}(x,\cdot)\in H^{1}(\boldsymbol{Q})$ the potential to $\upsilon_{\delta}(x)$ and $\eta$. Further, if $\hat{V}_{i,\eta}^{\omega}$ is the corresponding $\eta$-potential to $\hat{\upsilon}_{i}$, we find

[TABLE]

Since the mapping $\hat{\upsilon}_{i}\mapsto\hat{V}_{i,\eta}^{\omega}$ is linear, we find $V_{\eta,\varepsilon,\delta}^{\omega}\in C^{1}_{0}(\boldsymbol{Q};H^{1}(\boldsymbol{Q}))$ with

[TABLE]

For the first term on the right hand side we obtain

[TABLE]

Since the last expression on the right hand side converges to [math] as $\eta\to 0$ by Lemma 4.3, we find that $\nabla V_{\eta,\varepsilon,\delta}(x,x)\stackrel{{\scriptstyle{\scriptstyle 2s}}}{{\rightharpoonup}}\upsilon_{\varepsilon,\delta}(x,\omega)$.

Hence, we find for $\varepsilon$ small enough that $V_{\eta,\varepsilon,\delta}^{\omega}(x,x)$ is a valid point of evalutation for $F_{\eta}$ and

[TABLE]

On the other hand, due to Lemma 5.1, we have

[TABLE]

Since $f$ is continuous in $\xi$, we obtain from successively passing to the limits $\varepsilon\to 0$ and $\delta\to 0$ that

[TABLE]

∎

6 Homogenized system of equations

In this section, we are in the setting $\mu_{\omega}={\cal L}$ for all $\omega$. Hence, we frequently use the notations introduced in Remark 4.1.

The model equations of the problem (the microscopic problem) are

[TABLE]

together with the homogeneous Dirichlet boundary condition

[TABLE]

and the initial condition

[TABLE]

Now, we state the main result on the stochastic homogenization of the weak solution $(u_{\eta},\sigma_{\eta},z_{\eta})$ of problem (49) - (53).

Theorem 6.1**.**

Suppose that all assumptions of Theorem 3.1 are satisfied. Let $(u_{\eta},\sigma_{\eta},z_{\eta})$ be a weak solution of the initial-boundary value problem (49) - (53). Then, there exist

[TABLE]

such that (up to a subsequence)

[TABLE]

The weak two-scale limit function $(u_{0},u_{1},\sigma_{0},z_{0})$ solves the following homogenized problem:

[TABLE]

which hold for $(x,\omega,t)\in\boldsymbol{Q}\times\Omega\times[0,T_{e}]$, and with the boundary condition

[TABLE]

Moreover, the following variational inequality holds ($\tilde{\mathbb{A}}:=\tilde{\mathbb{C}}^{-1}$)

[TABLE]

where $(v^{(0)},v_{1})\in H^{1}(\boldsymbol{Q},{\mathbb{R}}^{3})\times L^{2}(\boldsymbol{Q},L^{2}_{pot}(\Omega))$ and $\sigma^{(0)}\in L^{2}(\boldsymbol{Q},L^{2}_{sol}(\Omega))$ solve the linear elasticity problem

[TABLE]

The careful reading of the proof of Theorem 3.1 suggests the following result.

Proposition 6.1**.**

Suppose that all assumptions of Theorem 3.1 are satisfied. Then, the weak solution $(u_{\eta},\sigma_{\eta},z_{\eta})$ of problem (49) - (53) (in the sense of Definition 3.1) fullfills the uniform estimates

[TABLE]

The result of Proposition 6.1 plays an important role in the proof of Theorem 6.1 below.

Proof of Theorem 6.1

Proof.

Proposition 6.1 provides the required uniform estimates for the solution of the microscopic problem (49) - (53). Therefore, due to Lemma 4.10 there exist functions $u_{0}$, $u_{1}$, $\sigma_{0}$ and $z_{0}$ with the prescribed regularities in Theorem 6.1 such that the convergence results in (54) hold. We note that (50) gives equation (56), namely

[TABLE]

Next, we test equation (49) with a function $\phi\in C_{0}^{\infty}(\boldsymbol{Q},{\mathbb{R}}^{3})$. Passing to the stochastic two-scale limit in the integral identity corresponding to (49) yields

[TABLE]

Now, we consider $\phi_{\tilde{\omega},\eta}(x,t)=\eta\phi(x,t)V_{\eta}^{\tilde{\omega}}(x)$, where $\phi\in C_{0}^{\infty}(\boldsymbol{Q}_{T_{e}},{\mathbb{R}})$ and $\upsilon\in L^{q}_{\rm pot}(\Omega)$ with potential $V_{\eta}^{\tilde{\omega}}$ given by Lemma 4.3, as another test function in (49) and obtain

[TABLE]

The stochastic two-scale limit in equation (6) yields

[TABLE]

Equation (65) implies that the integral identity

[TABLE]

holds for ever $\upsilon\in L^{q}_{\rm pot}(\Omega)$ for a.e. $(x,t)\in{\boldsymbol{Q}}\times(0,T_{e}).$ Integral equality (66) yields that $\sigma_{0}(x,\cdot,t)\in L^{2}_{\rm sol}(\Omega;{\cal S}^{3})$ for a.e. $(x,t)\in{\boldsymbol{Q}}\times(0,T_{e})$.

To pass to the stochastic two-scale limit in the inequality

[TABLE]

we use the results of Lemma 4.9 and Lemma 5.1 and obtain

[TABLE]

where $\sigma^{(0)}\in L^{2}(\boldsymbol{Q},L^{2}(\Omega))$ solves the following linear elasticity problem (59)–(61), which is obtained by the passage to the stochastic two-scale limit in equations (12) - (14). Here, $v^{(0)}\in H^{1}(\boldsymbol{Q},{\mathbb{R}}^{3})$ and $v_{1}\in L^{2}(\boldsymbol{Q},L^{2}_{\rm pot}(\Omega))$.

Therefore, we conclude that the limit function $(u_{0},u_{1},\sigma_{0},z_{0})$ satisfies the homogenized problem (55) - (57) and the variational inequality (6.1). ∎

Appendix: Rothe’s approximation functions

Here we recall the definition of Rothe’s approximation functions. For any family $\{\xi^{n}_{m}\}_{n=0,...,2m}$ of functions in a reflexive Banach space $X$ and for $h=T_{e}/m$, we define the piecewise affine interpolant $\xi_{m}\in C([0,T_{e}],X)$ by

[TABLE]

and the piecewise constant interpolant $\bar{\xi}_{m}\in L^{\infty}(0,T_{e};X)$ by

[TABLE]

For the further analysis we recall the following property of $\bar{\xi}_{m}$ and $\xi_{m}$:

[TABLE]

where $\bar{\xi}_{m}$ is formally extended to $t\leq 0$ by $\xi^{0}_{m}$ and $1\leq p\leq\infty$ (see [25]).

Acknowledgement

M.H. is financed by DFG within research project SFB 1114 – C05. N. S. acknowledges the financial support by the DFG within the research project NE 1498/4-1.