Crack growth by vanishing viscosity in planar elasticity

TL;DR

This paper establishes the existence of quasistatic crack evolutions in brittle materials using a vanishing viscosity method within planar linearized elasticity, allowing the crack path to emerge naturally.

Contribution

It introduces a novel vanishing viscosity approach to model crack growth without prescribing the crack path in planar elasticity.

Findings

Existence of quasistatic evolutions proven

Crack paths are determined as part of the solution

Analysis of energy release rate properties

Abstract

We show the existence of quasistatic evolutions in a fracture model for brittle materials by a vanishing viscosity approach, in the setting of planar linearized elasticity. The crack is not prescribed a priori and is selected in a class of (unions of) regular curves. To prove the result, it is crucial to analyze the properties of the energy release rate.

Equations

$$\Omega \times \mathbb{R}\ni (x_1,x_2,x_3)\mapsto (x_1+u_1(x_1,x_2),x_2+u_2(x_1,x_2),x_3),\text{where}u=(u_1,u_2):\Omega \to {\mathbb{R}}^2.$$

$${\overline{\mathrm{B}}}_{\eta }^1\cap {\overline{\mathrm{B}}}_{\eta }^2=\{x\}\text{and}({\mathrm{B}}_{\eta }^1\cup {\mathrm{B}}_{\eta }^2)\cap \Gamma =\varnothing .$$

$$\mathrm{B}_{2\eta}(x)\cap\Big{(}\partial\Omega\cup\bigcup_{l\neq m}\Gamma^{l}\Big{)}=\emptyset\,.$$

$$d_H({\Gamma }^{\prime };\Gamma ):=\max \left\{\underset{x\in {\Gamma }^{\prime }}{\sup }\mathrm{dist}(x,\Gamma ),\underset{x\in \Gamma }{\sup }\mathrm{dist}(x,{\Gamma }^{\prime })\right\},$$

$$\min \left\{\frac{1}{2}{\int }_{\Omega \setminus \Gamma }\mathbb{C}\mathrm{E}v:\mathrm{E}v\mathrm{d}x-{\int }_{\Omega \setminus \Gamma }f(t)\cdot v\mathrm{d}x-{\int }_{{\partial }_S\Omega }g(t)\cdot v\mathrm{d}{\mathcal{H}}^1:v\in H^1(\Omega \setminus \Gamma ;{\mathbb{R}}^2),v=w(t)\text{on}{\partial }_D\Omega \right\}$$

$$\mathcal{E}(t;\Gamma ):=\frac{1}{2}{\int }_{\Omega \setminus {\Gamma }_s}\mathbb{C}\mathrm{E}u(t;\Gamma ):\mathrm{E}u(t;\Gamma )\mathrm{d}x-{\int }_{\Omega \setminus \Gamma }f(t)\cdot u(t;\Gamma )\mathrm{d}x-{\int }_{{\partial }_S\Omega }g(t)\cdot u(t;\Gamma )\mathrm{d}{\mathcal{H}}^1.$$

$$\mathcal{F}(t;\Gamma ):=\mathcal{E}(t;\Gamma )+\mathcal{K}(\Gamma ),\mathcal{K}(\Gamma ):={\int }_{\Gamma }\kappa \mathrm{d}{\mathcal{H}}^1,$$

$$0=t_{k,0}<t_{k,1}<\cdots <t_{k,k}=T\text{and}\underset{k\to \infty }{\lim }\underset{1\le i\le k}{\max }(t_{k,i}-t_{k,i-1})=0.$$

$$\min \left\{\mathcal{E}(t_{k,i};\Gamma )+{\mathcal{H}}^1(\Gamma )+\frac{\epsilon }{2}\sum _{m=1}^M\frac{{\mathcal{H}}^1({\Gamma }^m\setminus {\Gamma }_{\epsilon ,k,i-1}^m{)}^2}{t_{k,i}-t_{k,i-1}}:\Gamma \in {\mathcal{R}}_{\eta },\Gamma \supseteq {\Gamma }_{\epsilon ,k,i-1}\right\},$$

$${\Gamma }_{\epsilon ,k}(0):={\Gamma }_0,{\Gamma }_{\epsilon ,k}(t):={\Gamma }_{\epsilon ,k,i}\text{for}t\in (t_{k,i-1},t_{k,i}].$$

$$\frac{\mathrm{d}\mathcal{E}(t;\Gamma_{\sigma})}{\mathrm{d}\sigma}\bigg{|}_{\sigma=s}=\frac{\mathrm{d}\mathcal{E}(t;\hat{\Gamma}_{\sigma})}{\mathrm{d}\sigma}\bigg{|}_{\sigma=s}.$$

$$\begin{array}{rl}\mathcal{F}(\tilde{t}(s);\tilde{\Gamma }(s))=& \mathcal{F}(0;{\Gamma }_0)+{\int }_0^s{\int }_{\Omega }\mathbb{C}\mathrm{E}\tilde{u}(\sigma ):\mathrm{E}\dot{w}(\tilde{t}(\sigma )){\tilde{t}}^{\prime }(\sigma )\mathrm{d}x\mathrm{d}\sigma \\ & -\sum _{m=1}^M{\int }_0^s({\tilde{\mathcal{G}}}^m(\sigma )-\kappa ({\tilde{P}}^m(\sigma )))({\tilde{l}}^m{)}^{\prime }(\sigma )\mathrm{d}\sigma \\ & -{\int }_0^s{\int }_{\Omega }\dot{f}(\tilde{t}(\sigma ))\cdot \tilde{u}(\sigma ){\tilde{t}}^{\prime }(\sigma )\mathrm{d}x\mathrm{d}\sigma -{\int }_0^s{\int }_{\Omega }f(\tilde{t}(\sigma ))\cdot \dot{w}(\tilde{t}(\sigma )){\tilde{t}}^{\prime }(\sigma )\mathrm{d}x\mathrm{d}\sigma \\ & -{\int }_0^s{\int }_{{\partial }_S\Omega }\dot{g}(\tilde{t}(\sigma ))\cdot \tilde{u}(\sigma ){\tilde{t}}^{\prime }(\sigma )\mathrm{d}{\mathcal{H}}^1\mathrm{d}\sigma -{\int }_0^s{\int }_{{\partial }_S\Omega }g(\tilde{t}(\sigma ))\cdot \dot{w}(\tilde{t}(\sigma )){\tilde{t}}^{\prime }(\sigma )\mathrm{d}{\mathcal{H}}^1\mathrm{d}\sigma .\end{array}$$

$$\text{if}{\tilde{t}}^{\prime }(\sigma )>0,\text{then}\tilde{\Gamma }(\sigma )=\Gamma (\tilde{t}(\sigma )),$$

$$\Vert \nabla u{\Vert }_2\le C(\Vert u{\Vert }_2+\Vert \mathrm{E}u{\Vert }_2).$$

$$\Vert \nabla u{\Vert }_2\le C(\Vert u{\Vert }_2+\Vert \mathrm{E}u{\Vert }_2)\text{for every u∈H1(Ω∖Γn;R2)}.$$

$$\Vert \nabla u{\Vert }_2\le C\Vert \mathrm{E}u{\Vert }_2\text{and}\Vert u{\Vert }_2\le C\Vert \mathrm{E}u{\Vert }_2.$$

$$\Vert \nabla u{\Vert }_{2,{\Omega }^{\prime }}\le C^{\prime }(\Vert u{\Vert }_{2,{\Omega }^{\prime }}+\Vert \mathrm{E}u{\Vert }_{2,{\Omega }^{\prime }})\text{for every u∈H1(Ω∖Γn;R2)}.$$

$$\Vert \nabla u{\Vert }_{2,{\Omega }_1^n\setminus \overline{{\Omega }^{\prime }}}\le \tilde{C}(\Vert u{\Vert }_{2,{\Omega }_1^n\setminus \overline{{\Omega }^{\prime }}}+\Vert \mathrm{E}u{\Vert }_{2,{\Omega }_1^n\setminus \overline{{\Omega }^{\prime }}})\text{for every u∈H1(Ω∖Γn;R2)}.$$

$$\Vert u{\Vert }_2\le C\Vert \mathrm{E}u{\Vert }_2\text{for every u∈H1(Ω∖Γn;R2) with u=0 H1-a.e. on ∂DΩ, n large enough}.$$

$${\mathcal{E}}_n=E_n(u_n)=\frac{1}{2}{\int }_{\Omega }\mathbb{C}\mathrm{E}{u}_n:\mathrm{E}{u}_n\mathrm{d}x-{\int }_{\Omega }{f}_n\cdot u_n\mathrm{d}x-{\int }_{{\partial }_S\Omega }{g}_n\cdot u_n\mathrm{d}{\mathcal{H}}^1\text{for n∈N∪{∞}},$$

$${\Lambda }_n(x):=x+{\lambda }_n(|x-{\Pi }_{{\Gamma }_{\infty }}(x)|)({\hat{\gamma }}_n(s(x))-{\gamma }_{\infty }(s(x)))\text{for x∈R2}.$$

$$E_n(u_n)\le E_n(v_n)<C.$$

$$\underset{n\to \infty }{\lim }{\int }_{\Omega }{f}_n\cdot u_n\mathrm{d}x+{\int }_{{\partial }_S\Omega }{g}_n\cdot u_n\mathrm{d}{\mathcal{H}}^1={\int }_{\Omega }{f}_{\infty }\cdot \phi \mathrm{d}x+{\int }_{{\partial }_S\Omega }{g}_{\infty }\cdot \phi \mathrm{d}{\mathcal{H}}^1.$$

$$E_{\infty }(\phi )\le E_{\infty }(v)\text{for every v∈H1(Ω∖Γ∞;R2) with v=w∞ H1-a.e. on ∂DΩ}.$$

$$E_{\infty }(u_{\infty })\le \underset{n}{\mathrm{liminf}}{E}_n(u_n)\le \underset{n}{\mathrm{limsup}}{E}_n(u_n)\le \underset{n}{\lim }{E}_n(u_{\infty }\circ {\Lambda }_n^{-1}+w_n-w_{\infty })=E_{\infty }(u_{\infty }),$$

$$\frac{\mathrm{d}\mathcal{E}(\Gamma_{\sigma})}{\mathrm{d}\sigma}\bigg{|}_{\sigma=s}\coloneqq\lim_{\sigma\to s}\frac{\mathcal{E}(\Gamma_{\sigma})-\mathcal{E}(\Gamma_{s})}{\sigma-s}\,.$$

$$\frac{\mathrm{d}\mathcal{E}(\Gamma_{\sigma})}{\mathrm{d}\sigma}\bigg{|}_{\sigma=s}=\frac{\mathrm{d}\mathcal{E}(\hat{\Gamma}_{\sigma})}{\mathrm{d}\sigma}\bigg{|}_{\sigma=s}\,.$$

$${\Gamma }_s:=\{\gamma (\sigma ):0\le \sigma \le s\},$$

$$F_{s,\delta}(x)\coloneqq x+\left(\begin{array}[]{cc}(\gamma_{1}(s+\delta)-\gamma_{1}(s))\varphi(x)\\[2.84526pt] \zeta(x_{1}+(\gamma_{1}(s+\delta)-\gamma_{1}(s))\varphi(x))-\zeta(x_{1})\end{array}\right)\,,$$

$$\displaystyle\rho_{s}(x)\coloneqq\partial_{\delta}(F_{s,\delta}(x))|_{\delta=0}=\gamma^{\prime}_{1}(s)\varphi(x)\left(\begin{array}[]{cc}1\\ \zeta^{\prime}(x_{1})\end{array}\right)$$

Full text

Introduction

In many applications of engineering, it is crucial to predict the propagation of fracture in structures and to understand whether cracks are stable. When the external loading is very slow if compared with the time scale of internal oscillations (such as, e.g., in a building in standard conditions), it is possible to ignore inertia and to assume that the system is always at equilibrium: the resulting model si called quasistatic. Quasistatic (or rate-independent) processes have been extensively analyzed in the mathematical literature both in the context of fracture and of other models (see [32] and references therein).

The first difficulties in modeling fracture are related to identifying equilibrium configurations. In fact, in order to state that a configuration is stable, one would have to use a derivative of the mechanical energy with respect to the crack set, which is not well defined. Thus one may prefer a derivative-free formulation where equilibria are restricted to global minimizers (of the sum of the mechanical energy and of the dissipated energy due to crack growth), in the context of energetic solutions to rate-independent systems, see e.g. [18, 14, 4, 17, 11, 12, 19].

A second approach allows one to take into account of more equilibria by restricting the set of the admissible cracks. In fact, the problem is to select a class of regular curves and to prove the existence of a derivative of the mechanical energy with respect to the elongation of a crack in that class. The opposite of this derivative is called energy release rate and represents the gain in stored elastic energy due to an infinitesimal crack growth. Griffith’s criterion [20] allows crack growth only when the energy release rate reaches the toughness of the material (i.e., the energy spent to produce an infinitesimal crack).

In this context, some existence results for crack evolution were given in the case of antiplane linear elasticity, where the deformation is represented by a scalar function (that is the vertical displacement, depending on the two horizontal components, while the horizontal displacement is zero). The results of [22] and [34] deal with the case of a prescribed crack, i.e., before the evolution starts one already knows the set which is going to crack. An algorithm for predicting a stable crack path (chosen from a class of regular curves) was proposed in [28] and extended in [8] to a class of curves with branches and kinks. A nonlinear model with vector displacements (still in dimension two) was studied in [25] for a prescribed crack path.

In this paper we prove an existence result for crack evolution based on Griffith’s criterion, in the context of planar linear elasticity, in dimension two. In this case the displacement is a vector (with two components). In our model, the path followed by the crack is not a priori known. In fact, the crack is assumed to be the union of a fixed number of $C^{1,1}$ curves and is selected among a class of (unions of) curves with bounded curvature, with no self-intersections, and with at most one point meeting the boundary of the domain (in the reference configuration). Some geometric constraints guarantee that this class is compact with respect to the Hausdorff convergence of sets. The same class of admissible cracks was considered in [28].

In order to write the flow rule for crack propagation, we need the expression of the energy release rate. The first step is to prove that, when the crack is a prescribed curve, then the mechanical energy (i.e., the sum of the stored elastic energy and of the work of external volume and surface forces) is differentiable with respect to the arc length of the curve, and its derivative can be written as a surface integral depending on the deformation gradient (Proposition 3.1). Since we want a model predicting the crack path (not prescribed a priori), we need to prove that the energy release rate is independent of the extension of the crack (in the class of $C^{1,1}$ curves). This is done in Theorem 3.6. Moreover, the energy release rate is continuous with respect to the Hausdorff convergence of cracks (see Remark 3.11). When there are more curves, there is an energy release rate for each crack tip.

Proving such properties of the energy release rate(s) is fundamental to study quasistatic crack evolution and is the major technical difficulty of this work. In fact, the strategy of the proof differs from the method used in the corresponding results in the antiplane case, cf. [27, 28]. In planar elasticity, assuming that the crack is $C^{\infty}$, that there are no external forces, and that the elasticity tensor is constant, it was proven in [2] that the energy release rate can be expressed in terms of two stress intensity factors, which characterize the singularities of the elastic equilibrium; since the stress intensity factors only depend on the current crack, it turns out that the energy release rate is independent of the crack’s extension. In this paper we need a corresponding property for $C^{1,1}$ cracks (in the class where we have compactness with respect to Hausdorff convergence) and for energies with external forces and nonconstant elasticity tensor. The same strategy does not apply to the nonsmooth case, in particular we do not prove the existence of the stress intensity factors; nonetheless, we prove that the energy release rate is stable under Hausdorff convergence in the class of $C^{1,1}$ cracks, so we can employ the results of [2] via some approximation arguments (see Section 3).

We point out that an energy release rate associated with a crack tip does exist also under much weaker regularity conditions on the crack set. For instance, the results of [3] apply to cracks that are merely closed and connected. However, in this setting energy release rates can be characterized just up to subsequences through a blow-up limit, thus uniqueness is not guaranteed and, ultimately, the independence on extensions may not hold. On the other hand, the results of [5] do not have this limitation, but the initial crack needs to be straight, which makes it impossible to use such characterization in the context of an evolution problem. (We also refer to [7] for related results in antiplane elasticity.) For these reasons in this paper we resort to the class of (unions of) $C^{1,1}$ cracks where, as mentioned, better properties can be proven.

This allows us to employ the well known vanishing viscosity method for finding balanced viscosity solutions to rate-independent systems, see [22, 28, 8] for fracture in antiplane elasticity and [32] for further references. We fix a time discretization and solve some incremental problems where we minimize the sum of the mechanical energy and of the dissipated energy. Notice that in the present work the dissipated energy density is nonconstant and depends on the position of the crack tip in the reference configuration. In the minimum problems, the total energy is perturbed with a term penalizing brutal propagations between energy wells, multiplied by a parameter $\varepsilon$. Passing to the continuous time, we obtain a viscous version of Griffith’s criterion, with a regularizing term multiplied by $\varepsilon$; a second passage to the limit as $\varepsilon\to 0$ leads to rate-independent solutions. It is also possible to characterize the time discontinuities of the resulting evolution using the reparametrization technique first proposed in [15] and then refined in [29, 30, 31, 33].

The main result of this paper, extending the results of [28] to planar elasticity, is the existence of a quasistatic evolution (more precisely, a balanced viscosity evolution) fulfilling Griffith’s criterion: the length of each component of the crack is a nondecreasing function of time; at all continuity points of these functions, the energy release rate at each tip is less than or equal to the material’s toughness at that tip (which is a stability condition); the length is increasing only if the energy release rate reaches the toughness. Moreover, time discontinuities (corresponding to brutal propagation) can be interpolated by a transition, characterized by a viscous flow rule, where the energy release rates are larger than or equal to the toughness (see e.g. [9, 10, 23, 24] for corresponding results in damage and plasticity).

Notation

Given two vectors $a,b\in\mathbb{R}^{d}$, their scalar product is denoted by $a\cdot b$. We set $\mathbb{M}^{d}$ the space of $d\times d$ square matrices, and we denote by $\mathbb{M}^{d}_{sym}$ and $\mathbb{M}^{d}_{skw}$ the subsets of symmetric and skew-symmetric ones, respectively. We set $\mathbf{I}$ the identity matrix in $\mathbb{M}^{d}$. Given $A$ and $B$ in $\mathbb{M}^{d}$, we write $A:B$ to denote their Euclidean scalar product, namely $A:B=A_{ij}B_{ij}$. Here and in the rest of the paper we adopt the convention of summation over repeated indices. For every $p\geq 1$ we define the $p$-norm in $\mathbb{R}^{d}$ as |x|_{p}:=\big{(}\sum_{i=1}^{d}|x_{i}|^{p}\big{)}^{1/p}. The 2-norm will be simply denoted by $|\cdot|$. The latter induces the distance $\mathop{\rm dist}\nolimits(C,D):=\inf\{|x-y|\,:\,\,x\in C,\,y\in D\}$ between two sets $C$ and $D$. The maximal distance between two points of a set $E$, namely its diameter, is denoted by $\mathop{\mathrm{diam}}(C)$.

The symbol $\mathrm{B}_{\rho}(x)$ denotes the open ball of radius $\rho$ in $\mathbb{R}^{2}$, centred at $x$. The support of a function $f$, namely the closure of $\{f\neq 0\}$, is denoted by $\mathop{\rm spt}\nolimits(f)$. For a tensor field $V\in C^{1}(\mathbb{R}^{d};\mathbb{M}^{d})$, by $\mathop{\rm div}\nolimits V$ we mean its divergence with respect to lines, namely $(\mathop{\rm div}\nolimits V)_{i}:=\partial_{j}V_{ij}$. The symmetric gradient of a vector field $u\in C^{1}(\mathbb{R}^{d};\mathbb{R}^{d})$ is denoted by $\mathrm{E}u$, namely $(\mathrm{E}u)_{ij}:=(\partial_{i}u_{j}+\partial_{j}u_{i})/2$.

We adopt standard notations for Lebesgue and Sobolev spaces on a bounded open set of $\mathbb{R}^{d}$. The boundary values of a Sobolev function are always intended in the sense of traces. Boundary integrals on Lipschitz curves are done with respect to the 1-dimensional Hausdorff measure $\mathcal{H}^{1}$. Given an interval $I\subset\mathbb{R}$ and a Banach space $X$, $L^{p}(I;X)$ is the space of $L^{p}$ functions from $I$ to $X$. Similarly, the sets of continuous and absolutely continuous functions from $I$ to $X$ are denoted by $C^{0}(I;X)$ and $AC(I;X)$, respectively. Derivatives of functions depending on one variable are denoted by a prime or, when the variable is time, by a dot.

The identity map in a vector space is denoted by $id$. Given a normed vector space $X$ the norm in $X$ is denoted by $\|\cdot\|_{X}$. We adopt the same notation also for vector valued functions in $X$. For brevity, the norm in $L^{p}$ over an open set $\Omega$ of $\mathbb{R}^{d}$ is denoted by $\|\cdot\|_{p,\Omega}$ or, when no ambiguity may arise, simply by $\|\cdot\|_{p}$.

1. Description of the model and existence results

We describe a crack model in planar elasticity for a brittle body. The body is represented in its reference configuration by an infinite cylinder $\Omega\times\mathbb{R}$, where $\Omega\subset\mathbb{R}^{2}$ is a bounded connected open set, with Lipschitz boundary. By assumption, the displacement $u$ produced by the external loading is horizontal and depends only on the two horizontal components: the deformation is then given by

[TABLE]

1.1. Admissible cracks

The set of possible discontinuity points of $u$ (the crack) is assumed to lie in a class of admissible regular cracks. We now define such class following [27, 28]. It depends on a parameter $\eta>0$ that is thought as small, but is fixed throughout the paper.

Definition 1.1**.**

Fixed $\eta>0$, the set $\mathcal{R}_{\eta}^{0}$ contains all closed subsets $\Gamma\subseteq\overline{\Omega}$ such that

• (a)

$\Gamma$ is a union of a finite number of arcs of $C^{1,1}$ curves, each of them intersecting $\partial\Omega$ in at most one endpoint,

• (b)

$\mathcal{H}^{1}(\Gamma\cap\Omega)>0$ and $\Omega\setminus\Gamma$ is a connected open set, union of a finite number of Lipschitz domains,

• (c)

for every $x\in\Gamma$ there exist two open balls $\mathrm{B}_{\eta}^{1},\mathrm{B}_{\eta}^{2}\subseteq\mathbb{R}^{2}$ of radius $\eta$ such that

[TABLE]

Furthermore, we denote with $\mathcal{R}^{0,1}_{\eta}$ the class of curves $\Gamma\in\mathcal{R}_{\eta}^{0}$ such that $\Gamma$ is one arc of curve of class $C^{1,1}$ intersecting $\partial\Omega$ in exactly one endpoint.

Notice that $\mathcal{R}_{\eta}^{0}\subseteq\mathcal{R}_{\eta^{\prime}}^{0}$ if $\eta>\eta^{\prime}$. The role of (1.1) is twofold: on the one hand it gives a uniform bound (depending on $1/\eta$) on the curvature of each connected component of any set $\Gamma\in\mathcal{R}_{\eta}^{0}$, on the other hand it ensures that each of these components is an arc of a simple curve, i.e., a curve with no self-intersections. Because of $(a)$, each of the arcs has one or two endpoints contained in $\Omega$; we say that these points are the crack tips.

Since quasistatic models are in general unable to predict crack initiation [6], i.e., nucleation of a new crack from sound material, we assume that there is an initial crack $\Gamma_{0}\in\mathcal{R}_{\eta}^{0}$. Each connected component of an admissible crack $\Gamma$ will be the extension of a connected component of $\Gamma_{0}$, starting from its crack tips. Let $M$ be the number of crack tips of $\Gamma_{0}$; notice that $M$ may be larger than the number of connected components of $\Gamma_{0}$. We parametrize $\Gamma_{0}$ by introducing $M$ injective functions $\gamma^{m}$ of class $C^{1,1}$, for $m=1,\dots,M$, in the following way:

• •

If a connected component $\Gamma_{0}^{m}$ of $\Gamma_{0}$ intersects $\partial\Omega$ in a single endpoint $x_{0}$, we consider its arc-length parametrization $\gamma^{m}\colon[0,\mathcal{H}^{1}(\Gamma_{0}^{m})]\to\Gamma_{0}^{m}$ such that $\gamma^{m}(0)\in\partial\Omega$ and $\gamma^{m}(\mathcal{H}^{1}(\Gamma_{0}^{m}))$ is the crack tip of $\Gamma_{0}^{m}$. In particular, a crack in $\mathcal{R}^{0,1}_{\eta}$ has exactly one tip.

• •

If a connected component of $\Gamma_{0}$ is contained in $\Omega$, we see it as the union of two curves $\Gamma_{0}^{m}$, $\Gamma_{0}^{m+1}$, intersecting at a single point $\bar{x}$ (that is not a tip of $\Gamma_{0}$); then we consider two arc-length parametrizations $\gamma^{m}\colon[0,\mathcal{H}^{1}(\Gamma_{0}^{m})]\to\Gamma_{0}^{m}$, $\gamma^{m+1}\colon[0,\mathcal{H}^{1}(\Gamma_{0}^{m+1})]\to\Gamma_{0}^{m+1}$, such that $\gamma^{m}(0)=\gamma^{m+1}(0)=\bar{x}$ and $\gamma^{m}(\mathcal{H}^{1}(\Gamma_{0}^{m}))$ and $\gamma^{m+1}(\mathcal{H}^{1}(\Gamma_{0}^{m+1}))$ are the two crack tips.

We then have $\Gamma_{0}=\bigcup_{m=1}^{M}\Gamma_{0}^{m}$. Analogous parametrizations will be used for the extensions of $\Gamma_{0}$. In the next definition, $M$ is the number fixed above.

Definition 1.2**.**

The set $\mathcal{R}_{\eta}$ contains all subsets $\Gamma\in\mathcal{R}_{\eta}^{0}$ such that

• (d)

$\Gamma$ is the union of $M$ connected subsets $\Gamma^{1},\dots,\Gamma^{M}$, such that any two of them intersect in up to a point,

• (e)

$\Gamma^{m}\supseteq\Gamma_{0}^{m}$ for every $m=1,\dots,M$,

• (f)

for every $m=1,\dots,M$ and for every $x\in\Gamma^{m}{\setminus}\Gamma_{0}^{m}$

[TABLE]

Given a set $\Gamma=\bigcup_{m=1}^{M}\Gamma^{m}\in\mathcal{R}_{\eta}$, we extend the functions $\gamma^{m}$, $m=1,\dots,M$, defined above, to arc-length parametrizations $\gamma^{m}\colon[0,\mathcal{H}^{1}(\Gamma^{m})]\to\Gamma^{m}$; it turns out that they are injective and of class $C^{1,1}$. Properties (a)–(f) ensure that the class $\mathcal{R}_{\eta}$ is sequentially compact with respect to Hausdorff convergence (see the next section for details), induced by the following distance.

Definition 1.3**.**

Given two compact subsets $\Gamma,\Gamma^{\prime}\subset\overline{\Omega}$, their Hausdorff distance is given by

[TABLE]

with the conventions $d_{H}(x;\emptyset)=\mathop{\mathrm{diam}}{\Omega}$ and $\sup\emptyset=0$. A sequence $(\Gamma_{n})_{n\in\mathbb{N}}$ of compact subsets of $\overline{\Omega}$ converges to $\Gamma$ in the Hausdorff metric if $d_{H}(\Gamma_{n};\Gamma)\to 0$ as $n\to\infty$.

Remark 1.4**.**

There are choices of $\Gamma_{0}$ such that $\mathcal{R}_{\eta}$ contains no elements different from $\Gamma_{0}$: we mention a few examples with $\Omega=[-1,1]^{2}$. Let $\Gamma_{0}=[-1,0]\times\{0,\frac{1}{4}\}$: then $\Gamma_{0}\in\mathcal{R}_{\eta}^{0}$ only if $\eta\leq 1/8$, thus $\mathcal{R}_{\eta}=\emptyset$ if and only if $\eta>1/8$. If instead $\Gamma_{0}=[-1,0]\times\{0\}$, we have $\mathcal{R}_{\eta}=\{\Gamma_{0}\}$ if and only if $\eta\geq 1/2$. However, given $\Gamma_{0}$ such that $\mathcal{R}_{\eta}$ is trivial, one can find $\eta^{\prime}<\eta$ such that $\mathcal{R}_{\eta^{\prime}}$ contains nontrivial extensions of $\Gamma_{0}$. Starting from an initial crack with nontrivial extensions, the model described in this paper is reliable as long as our algorithm finds a current configuration $\Gamma(t)$ such that there are nontrivial extensions. If, during the evolution, some tip becomes $(2\eta)$-close to $\partial\Omega$ or to other connected components of the crack, the results should not be regarded as meaningful.

1.2. The mechanical energy and the incremental scheme

Since the body is brittle, the uncracked part $\Omega\setminus\Gamma$ is elastic; we assume that the displacements are small (so we adopt the setting of linear elasticity) and the crack is traction-free. We look for evolutions in the time interval $[0,T]$, produced by the time-dependent external loading:

• (H1)

a boundary condition $w\in AC([0,T];H^{1}(\Omega{\setminus}\Gamma_{0};\mathbb{R}^{2}))$, to be satisfied on a relatively open subset $\partial_{D}\Omega$ of $\partial\Omega$, with a finite number of connected components,

• (H2)

a volume force $f\in AC([0,T];L^{2}(\Omega{\setminus}\Gamma;\mathbb{R}^{2}))$ and a surface force $g\in AC([0,T];L^{2}(\partial_{S}\Omega;\mathbb{R}^{2}))$, where $\partial_{S}\Omega$ is a relatively open subset of $\partial\Omega$ such that $\partial_{S}\Omega\Subset\partial\Omega{\setminus}\overline{\partial_{D}\Omega}$.

Without loss of generality, we assume that $\mathop{\rm spt}\nolimits(w)\subseteq\{x\in\overline{\Omega}\ :\ \mathop{\rm dist}\nolimits(x,\partial\Omega)\leq\eta\}$, so $w\equiv 0$ around any crack tip.

At each point $x\in\Omega$, the stress tensor is $\mathbb{C}(x)\colon\mathbb{M}^{2}_{sym}\to\mathbb{M}^{2}_{sym}$, where

• (H3)

$\mathbb{C}(x)A=\lambda(x)\mathrm{tr}(A)\mathbf{I}+2\mu(x)A$ for every $A\in\mathbb{M}^{2}_{sym}$, with $\lambda,\mu\in C^{0,1}(\overline{\Omega})$ such that $\mu(x)>0$ and $\lambda(x)+\mu(x)>0$ for every $x\in\overline{\Omega}$.

Notice that the standard conditions $\mu(x)>0$ and $\lambda(x)+\mu(x)>0$ ensure the positive definiteness of $\mathbb{C}(x)$, uniformly in $x$.

Given $t\in[0,T]$ and $\Gamma\in\mathcal{R}_{\eta}$, the minimum problem

[TABLE]

has a unique solution, denoted by $u(t;\Gamma)\colon\Omega{\setminus}\Gamma\to\mathbb{R}^{2}$, with elastic energy

[TABLE]

According to the assumption of brittle behavior, in order to produce a crack the system employs an energy depending (only) on the geometry of the crack itself, in the context of Griffith’s theory [20]. The total energy of the configuration corresponding to a crack $\Gamma$ at time $t$ is

[TABLE]

where the surface energy density satisfies

• (H4)

$\kappa\in C^{0}(\overline{\Omega};[\kappa_{1},\kappa_{2}])$,

where $0<\kappa_{1}<\kappa_{2}$.

Starting from the initial condition $\Gamma_{0}$ fixed above, we define a discrete-time evolution of stable states by solving some incremental minimum problems. For every $k\in\mathbb{N}$ we consider a subdivision of the time interval $[0,T]$ in nodes $\{t_{k,i}\}_{0\leq i\leq k}$ such that

[TABLE]

Fixed $\varepsilon>0$, we define by recursion the sets $\Gamma_{\varepsilon,k,i}$, $i=0,\dots,k$, as follows. We set $\Gamma_{\varepsilon,k,0}:=\Gamma_{0}$; for $i\geq 1$, $\Gamma_{\varepsilon,k,i}$ is a solution to the minimum problem

[TABLE]

where the role of the term multiplied by $\varepsilon$ is to penalize transitions between energy wells. The existence of solutions to (1.4) is proven in Corollary 2.5 exploiting the compactness properties of $\mathcal{R}_{\eta}$ with respect to the Hausdorff convergence, see Section 2 for details.

We define a piecewise constant interpolation on $[0,T]$ by

[TABLE]

The unilateral constraint $\Gamma\supseteq\Gamma_{\varepsilon,k,i-1}$ in (1.4) enforces irreversibility of the crack growth, indeed the set function $t\mapsto\Gamma_{\varepsilon,k}(t)$ is nondecreasing with respect to the inclusion.

1.3. Existence results

Passing to the limit as $k\to\infty$ and exploiting again the compactness of $\mathcal{R}_{\eta}$, we obtain a time-continuous evolution $t\mapsto\Gamma_{\varepsilon}(t)$. In order to understand its properties, we need to define the energy release rate associated to a crack.

For simplicity, let us first consider the case of a prescribed curve with only one tip. Given an increasing family of cracks $\Gamma_{\sigma}\in\mathcal{R}_{\eta}^{0,1}$ parametrized by their arc length $\sigma\in[0,S]$, we will prove that the map $\sigma\mapsto\mathcal{E}(t;\Gamma_{\sigma})$ is differentiable for every fixed $t$. Moreover, we will show that the derivative only depends on the current configuration $\Gamma_{s}$, and not on its possible extensions, i.e., if $\Gamma_{\sigma}=\hat{\Gamma}_{\sigma}$ for $\sigma\leq s$, then

[TABLE]

In particular, we are allowed to write -\frac{\mathrm{d}\mathcal{E}(t;\Gamma_{\sigma})}{\mathrm{d}\sigma}\bigg{|}_{\sigma=s}=:\mathcal{G}(t;\Gamma_{s}) with no ambiguity. The quantity $\mathcal{G}(t;\Gamma_{s})$ is the energy release rate corresponding to the crack $\Gamma_{s}$ and represents the (partial) derivative of the energy $\mathcal{E}$ with respect to variations of crack in the set of all admissible curves $\mathcal{R}^{0,1}_{\eta}$ larger than $\Gamma_{s}$. For the details of these results, we refer to Section 3 below.

In the case of a curve with several connected components $\Gamma\in\mathcal{R}_{\eta}$, for every tip indexed by $m$ we define the $m$-th energy release rate $\mathcal{G}^{m}(t;\Gamma)$ as above, with respect to variations of the sole component $\Gamma^{m}$ of $\Gamma$. The energy release rate will be in this case a vector $\mathcal{G}(t;\Gamma)\coloneqq(\mathcal{G}^{1}(t;\Gamma),\ldots,\mathcal{G}^{M}(t;\Gamma))$.

The properties of the evolution $t\mapsto\Gamma_{\varepsilon}(t)$ are summarized in the next proposition, whose proof is postponed to Section 4.

Proposition 1.5**.**

Fix $\eta>0$, $\Gamma_{0}\in\mathcal{R}_{\eta}^{0}$, and $\varepsilon>0$. Assume (H1)–(H4). Let $\Gamma_{\varepsilon,k}$ be as in (1.5). Then there are a subsequence (not relabeled) of $\Gamma_{\varepsilon,k}$ and a set function $t\mapsto\Gamma_{\varepsilon}(t)\in\mathcal{R}_{\eta}$ such that $\Gamma_{\varepsilon,k}(t)$ converges to $\Gamma_{\varepsilon}(t)$ in the Hausdorff metric for every $t\in[0,T]$.

Set $\Gamma_{\varepsilon}(t)=\bigcup_{m=1}^{M}\Gamma_{\varepsilon}^{m}(t)$, with the conventions of Definition 1.2, and $l_{\varepsilon}^{m}(t):=\mathcal{H}^{1}(\Gamma_{\varepsilon}^{m}(t))$. Then for every $m=1,\dots,M$, and for a.e. $t\in[0,T]$

1. (G1)ε

$\dot{l}_{\varepsilon}^{m}(t)\geq 0$; 2. (G2)ε

$\kappa(P_{\varepsilon}^{m}(t))-\mathcal{G}_{\varepsilon}^{m}(t)+\varepsilon\,\dot{l}_{\varepsilon}^{m}(t)\geq 0$; 3. (G3)ε

$\dot{l}_{\varepsilon}^{m}(t)\,[\kappa(P_{\varepsilon}^{m}(t))-\mathcal{G}_{\varepsilon}^{m}(t)+\varepsilon\,\dot{l}_{\varepsilon}^{m}(t)]=0$,

where $\mathcal{G}_{\varepsilon}^{m}(t)$ is the energy release rate corresponding to $\Gamma_{\varepsilon}^{m}(t)$.

Moreover, along a suitable $\varepsilon$-subsequence, $\varepsilon\|\dot{l}^{m}_{\varepsilon}\|_{2}^{2}$ is bounded uniformly w.r.t. $\varepsilon$.

Properties (G1)ε–(G3)ε show that the term multiplied by $\varepsilon$ in (1.4) has a regularizing effect, indeed the flow rule for the evolution of $l_{\varepsilon}:=(l_{\varepsilon}^{1},\dots,l_{\varepsilon}^{M})$ features a time derivative of the unknown. For this reason the corresponding solutions are called viscous.

In the passage to the limit as $\varepsilon\to 0$, such viscous regularizing term vanishes, so the system follows an evolution of stable states. We thus obtain a balanced viscosity evolution. The next result is proven in Section 4.

Theorem 1.6**.**

Fix $\eta>0$ and $\Gamma_{0}\in\mathcal{R}_{\eta}^{0}$. Assume (H1)–(H4). For every $\varepsilon>0$, let $\Gamma_{\varepsilon}$ be the evolution found in Proposition 1.5. Then there are a subsequence (not relabeled) of $\Gamma_{\varepsilon}$ and a set function $t\mapsto\Gamma(t)\in\mathcal{R}_{\eta}$ such that $\Gamma_{\varepsilon}(t)$ converges to $\Gamma(t)$ in the Hausdorff metric for every $t\in[0,T]$.

Set $\Gamma(t)=\bigcup_{m=1}^{M}\Gamma^{m}(t)$, with the conventions of Definition 1.2, and $l^{m}(t):=\mathcal{H}^{1}(\Gamma^{m}(t))$. Then for every $m=1,\dots,M$

1. (G1)

for a.e. $t\in[0,T]$, $\dot{l}^{m}(t)\geq 0$; 2. (G2)

for every $t\in[0,T]$ of continuity for $l^{m}$, $\kappa(P^{m}(t))-\mathcal{G}^{m}(t)\geq 0$; 3. (G3)

for a.e. $t\in[0,T]$, $\dot{l}^{m}(t)\,[\kappa(P^{m}(t))-\mathcal{G}^{m}(t)]=0$,

where $\mathcal{G}^{m}(t)$ is the energy release rate corresponding to $\Gamma^{m}(t)$.

Properties (G1)–(G3) are a formulation of Griffith’s criterion for crack growth and show the stability of the evolution $t\mapsto l(t):=(l^{1}(t),\dots,l^{M}(t))$ in its continuity points. However, the function $t\mapsto l(t)$ may have discontinuities and Theorem 1.6 does not provide a characterization of jumps in time. The existence result is refined in the following theorem, where we show that there are a reparametrization of the time interval and a parametrized evolution, continuous in time, that interpolates $l$ and follows a viscous flow rule in the intervals corresponding to the discontinuities of $l$. The next theorem is proven in Section 5.

Theorem 1.7** (Griffith’s criterion).**

Fix $\eta>0$ and $\Gamma_{0}\in\mathcal{R}_{\eta}^{0}$. Assume (H1)–(H4). There are absolutely continuous functions $\tilde{t}\colon[0,S]\to[0,T]$ and $\widetilde{\Gamma}^{m}\colon[0,S]\to\mathcal{R}_{\eta}$, $m\in\{1,\ldots,M\}$, such that for a.e. $\sigma\in[0,S]$, setting $\widetilde{\Gamma}(\sigma)=\bigcup_{m=1}^{M}\widetilde{\Gamma}^{m}(\sigma)$, with the conventions of Definition 1.2, and $\tilde{l}^{m}(\sigma):=\mathcal{H}^{1}(\widetilde{\Gamma}^{m}(\sigma))$,

• (pG1)

$\tilde{t}^{\prime}(\sigma)\geq 0$* and $(\tilde{l}^{m})^{\prime}(\sigma)\geq 0$ for every $m=1,\dots,M$;*

• (pG2)

if $\tilde{t}^{\prime}(\sigma)>0$, then $\widetilde{\mathcal{G}}^{m}(\sigma)\leq\kappa(\widetilde{P}^{m}(\sigma))$ for every $m=1,\dots,M$;

• (pG3)

if $\tilde{t}^{\prime}(\sigma)>0$ and $(\tilde{l}^{m})^{\prime}(\sigma)>0$ for some $m\in\{1,\dots,M\}$, then $\widetilde{\mathcal{G}}^{m}(\sigma)=\kappa(\widetilde{P}^{m}(\sigma))$;

• (pG4)

if $\tilde{t}^{\prime}(\sigma)=0$, then there is $m\in\{1,\dots,M\}$ such that $(\tilde{l}^{m})^{\prime}(\sigma)>0$; moreover, for every $m$ with this property, we have $\widetilde{\mathcal{G}}^{m}(\sigma)\geq\kappa(\widetilde{P}^{m}(\sigma))$,

where $\widetilde{\mathcal{G}}^{m}(\sigma)$ is the energy release rate corresponding to $\widetilde{\Gamma}^{m}(\sigma)$. Moreover, denoting with $\tilde{u}(\sigma)$ the solution of (1.2) at time $\tilde{t}(\sigma)$ with a crack $\widetilde{\Gamma}(\sigma)$, for every $s\in[0,S]$ it holds

[TABLE]

Finally,

[TABLE]

where $\Gamma$ is the balanced viscosity evolution found in Theorem 1.6.

2. Preliminary results

In this section we collect some properties of the class of admissible cracks $\mathcal{R}_{\eta}$ and of the associated displacements. We recall that, given a crack, the associated displacement is the unique solution to the corresponding minimum problem (1.2).

As already mentioned in the previous section, the elements of $\mathcal{R}_{\eta}$ have no self-intersections, and, during the evolution, their crack tips stay uniformly far from the boundary and from the other connected components of the crack set. Moreover, it is easy to show that the curvature and the $\mathcal{H}^{1}$ measure of the elements of $\mathcal{R}_{\eta}$ are controlled from above by $\eta^{-1}$ and by some constant $C(\Omega,\Gamma_{0},\eta)$, respectively. Finally, as proven in [27, Proposition 2.9 and Remark 2.10], the class of admissible cracks $\mathcal{R}_{\eta}$ is sequentially compact with respect to the Hausdorff convergence introduced in Definition 1.3.

Theorem 2.1**.**

Every sequence $(\Gamma_{n})_{n\in\mathbb{N}}\subset\mathcal{R}_{\eta}$ admits (up to a subsequence) a limit $\Gamma_{\infty}\in\mathcal{R}_{\eta}$ in the Hausdorff metric. Moreover, along the subsequence (not relabeled), we have $\mathcal{H}^{1}(\Gamma_{n})\to\mathcal{H}^{1}(\Gamma)$ as $n\to\infty$.

In what follows we show the continuity of the elastic energy $\mathcal{E}$ w.r.t. Hausdorff convergence of the crack set $\Gamma\in\mathcal{R}_{\eta}$. This will in particular imply the existence of solutions for the incremental minimum problems (1.4).

We start with recalling in Proposition 2.2 a Korn inequality for $\Omega\setminus\Gamma$. In Proposition 2.3, instead, we show that, along sequences of cracks $\Gamma_{n}\in\mathcal{R}_{\eta}$ converging in the Hausdorff metric, such an inequality is independent of $n$. The study is carried out disregarding the time variable, which for brevity is omitted. Accordingly, the elastic energy associated to a fracture $\Gamma$ writes $\mathcal{E}(\Gamma)$. Furthermore, when explicitly needed, we highlight the dependence on the data by writing $\mathcal{E}(f,g,w,\mathbb{C};\Gamma)$ for $\mathcal{E}(\Gamma)$.

Proposition 2.2**.**

Let $\Gamma\in\mathcal{R}_{\eta}$. Then, there exists a positive constant $C=C(\Omega,\Gamma)$ such that for every $u\in H^{1}(\Omega\setminus\Gamma;\mathbb{R}^{2})$

[TABLE]

Proof.

Being $\Omega\setminus\Gamma$ connected by arcs (see Definition 1.2), it is possible to fix $\widehat{\Gamma}\supset\Gamma$ such that $\Omega\setminus\widehat{\Gamma}$ is the union of $N$ disjoint open sets $\Omega_{i}$ with Lipschitz boundaries $\partial\Omega_{i}$ such that $\mathcal{H}^{1}(\partial_{D}\Omega\cap\partial\Omega_{i})>0$ for $i\in\{1,\ldots,N\}$, and apply the usual Korn inequality to $u$ restricted to $\Omega_{i}$. ∎

Proposition 2.3**.**

Let $\Gamma_{n},\Gamma_{\infty}\in\mathcal{R}_{\eta}$ be such that $\Gamma_{n}$ converges to $\Gamma_{\infty}$ in the Hausdorff metric as $n\to\infty$. Then, there exists a positive constant $C=C(\Omega)$ (independent of $n$) such that for $n$ sufficiently large

[TABLE]

Moreover, for $u\in H^{1}(\Omega\setminus\Gamma_{n};\mathbb{R}^{2})$ with $u=0$ $\mathcal{H}^{1}$-a.e. on $\partial_{D}\Omega$ we have

[TABLE]

Proof.

At least for $n$ sufficiently large, we may assume that there exists an extension $\widehat{\Gamma}_{n}$ of $\Gamma_{n}$ such that $\Omega\setminus\widehat{\Gamma}_{n}=\bigcup_{i=1}^{N}\Omega_{i}^{n}$, where $\Omega_{i}^{n}$ ($i=1,\dots,N$) are open bounded disjoint sets with Lipschitz boundaries and Lipschitz constant $L$ independent of $n$. Moreover, we can assume that $\mathcal{H}^{1}(\partial_{D}\Omega\cap\partial\Omega^{n}_{i})>0$ for $i\in\{1,\ldots,N\}$ and every $n$. The same construction can be repeated for $n=\infty$ in such a way that $\Omega_{i}^{n}$ converges to $\Omega_{i}^{\infty}$ in the Hausdorff metric as $n\to\infty$.

Let us now fix $\Omega^{\prime}\Subset\Omega_{1}^{\infty}$. For $n$ large enough (including the case $n=\infty$), we have that $\Omega^{\prime}\Subset\Omega_{1}^{n}$. Hence, applying Proposition 2.2 in $\Omega^{\prime}$ we deduce that there exists a positive constant $C^{\prime}$ independent of $n$ such that

[TABLE]

Since $\Omega_{1}^{n}$ and $\Omega_{1}^{\infty}$ share the same Lipschitz constant $L$, applying locally, close to the boundary of $\Omega_{1}^{n}$ (resp. $\Omega_{1}^{\infty}$), the results of [16, Theorem 4.2], we also obtain that there exists a positive constant $\widetilde{C}$ such that

[TABLE]

The same inequality can be proven for $\Omega^{n}_{i}$, $i\geq 2$. Therefore, combining (2.3) and (2.4) we get (2.1) for some positive constant $C$ independent of $n\in\mathbb{N}\cup\{\infty\}$, $n$ large enough.

To prove (2.2) it is enough to show that

[TABLE]

We proceed with the usual contradiction argument. Assume that (2.5) is false. Then, for every $k\in\mathbb{N}$ there exist $n_{k}>n_{k-1}$ and $u_{k}\in H^{1}(\Omega\setminus\Gamma_{n_{k}};\mathbb{R}^{2})$ such that $\|u_{k}\|_{2}>k\|\mathrm{E}u_{k}\|_{2}$. Without loss of generality, we may assume that $\|u_{k}\|_{2}=1$. By (2.1) we deduce that $\|\nabla{u}_{k}\|_{2}$ is bounded. Hence, up to a subsequence, $\nabla{u}_{k}\rightharpoonup\varphi$ weakly in $L^{2}(\Omega;\mathbb{M}^{2})$ and $u_{k}\to u$ in $L^{2}_{loc}(\Omega\setminus\Gamma_{\infty};\mathbb{R}^{2})$, which implies that $u\in H^{1}_{loc}(\Omega\setminus\Gamma_{\infty};\mathbb{R}^{2})$ with $\nabla{u}=\varphi$. Since $\varphi\in L^{2}(\Omega;\mathbb{M}^{2})$, applying [13, Proposition 7.1] we deduce that $u\in H^{1}(\Omega\setminus\Gamma_{\infty};\mathbb{R}^{2})$. Since $\mathrm{E}u_{k}$ converges to [math] in $L^{2}(\Omega;\mathbb{M}^{2})$, we get that $\mathrm{E}u=0$ in $\Omega$. Thus, $u$ is a rigid movement in $\Omega$, i.e., there exist $\mathrm{A}\in\mathbb{M}^{2}_{skw}$ and $b\in\mathbb{R}^{2}$ such that $u=\mathrm{A}x+b$ for $x\in\Omega$. Moreover, setting $\Omega_{\eta}:=\{x\in\Omega:\,\mathrm{dist}(x,\partial\Omega)<\eta\}$, by Definition 1.2 we have $(\Gamma_{n}\cap\Omega_{\eta})\setminus\Gamma_{0}=\emptyset$ and $u_{k}\rightharpoonup u$ in $H^{1}(\Omega_{\eta}\setminus\Gamma_{0};\mathbb{R}^{2})$. Therefore, $u=0$ $\mathcal{H}^{1}$-a.e. on $\partial_{D}\Omega$, which implies that $u=0$. We claim that $\|u_{k}\|_{2}\to\|u\|_{2}$. Indeed, $\|u_{k}\|_{2,\Omega^{\prime}}\to\|u\|_{2,\Omega^{\prime}}$ for every $\Omega^{\prime}\Subset\Omega\setminus\Gamma_{\infty}$. By a simple reflection argument applied on both sides of the crack set $\Gamma_{n}$, we instead obtain that $\|u_{k}\|_{2,\Omega\setminus\overline{\Omega}^{\prime}}\to\|u\|_{2,\Omega\setminus\overline{\Omega}^{\prime}}$. Thus, $1=\|u_{k}\|_{2}\to\|u\|_{2}=0$, which is a contradiction. This concludes the proof of (2.2). ∎

We are now ready to prove the continuity of the energy $\mathcal{E}$ with respect to the crack set. The following lemma is actually stated in a more general setting. Indeed, we show the continuity of the displacement $u$ solution of (1.2) not only w.r.t. the Hausdorff convergence of sets in $\mathcal{R}_{\eta}$, but also w.r.t. the data of the problem, i.e., the applied forces, the boundary datum, and the elasticity tensor. Such a continuity result will be useful in the next section, where we prove the differentiability of $\mathcal{E}$ w.r.t. crack elongations by using some approximations.

Lemma 2.4**.**

Let $f_{n},f_{\infty}\in L^{2}(\Omega;\mathbb{R}^{2})$, $g_{n},g_{\infty}\in L^{2}(\partial_{S}\Omega;\mathbb{R}^{2})$, $w_{n},w_{\infty}\in H^{1}(\partial_{S}\Omega;\mathbb{R}^{2})$, $\mathbb{C}_{n},\mathbb{C}_{\infty}\in C^{0,1}(\overline{\Omega})$, $\Gamma_{n},\Gamma_{\infty}\in\mathcal{R}_{\eta}$, and $n\in\mathbb{N}$ be such that $f_{n}\to f_{\infty}$ strongly in $L^{2}(\Omega;\mathbb{R}^{2})$, $w_{n}\to w_{\infty}$ in $H^{1}(\Omega\setminus\Gamma_{0};\mathbb{R}^{2})$, $g_{n}\rightharpoonup g_{\infty}$ weakly in $L^{2}(\partial_{S}\Omega;\mathbb{R}^{2})$, $\mathbb{C}_{n}\to\mathbb{C}_{\infty}$ uniformly in $\overline{\Omega}$, and $\Gamma_{n}\to\Gamma_{\infty}$ in the Hausdorff metric, as $n\to\infty$.

Then, the energies $\mathcal{E}(f_{n},g_{n},w_{n},\mathbb{C}_{n};\Gamma_{n})$ converge to $\mathcal{E}(f_{\infty},g_{\infty},w_{\infty},\mathbb{C}_{\infty};\Gamma_{\infty})$ in the limit as $n\to\infty$. Moreover, the corresponding displacements $u_{n}$ and $u_{\infty}$, solutions to the associated minimum problems (1.2), satisfy $\nabla u_{n}\to\nabla u_{\infty}$ strongly in $L^{2}(\Omega;\mathbb{M}^{2})$.

Proof.

The proof is carried out following the steps of [35, Lemma 3.7]. The letter $C$ will denote a positive constant, which can possibly change from line to line.

For the sake of clarity, we consider cracks in $\mathcal{R}_{\eta}^{0,1}$. The proof can be easily generalized to the whole class $\mathcal{R}_{\eta}$. For brevity, we set $\mathcal{E}_{n}\coloneqq\mathcal{E}(f_{n},g_{n},w_{n},\mathbb{C}_{n};\Gamma_{n})$ and $\mathcal{E}_{\infty}\coloneqq\mathcal{E}(f_{\infty},g_{\infty},w_{\infty},\mathbb{C}_{\infty};\Gamma_{\infty})$; furthermore, along the proof we denote with $E_{n}$ and $E_{\infty}$ the functionals appearing in the minimization (1.2) with data $\{f_{n},g_{n},w_{n},\mathbb{C}_{n},\Gamma_{n}\}$ and $\{f_{\infty},g_{\infty},w_{\infty},\mathbb{C}_{\infty},\Gamma_{\infty}\}$, respectively. Clearly, we have

[TABLE]

where $\mathrm{E}u_{n}$ are interpreted as functions defined a.e. in $\Omega$.

Let $\gamma_{n}\in C^{1,1}([0,\ell_{n}];\mathbb{R}^{2})$ and $\gamma_{\infty}\in C^{1,1}([0,\ell_{\infty}];\mathbb{R}^{2})$ be the arc-length parametrizations of $\Gamma_{n}$ and $\Gamma_{\infty}$, respectively, where $\ell_{n}$ and $\ell_{\infty}$ denote the $\mathcal{H}^{1}$ measures of the crack sets. By a simple rescaling of $\gamma_{n}$, we construct a $C^{1,1}$ parametrization $\hat{\gamma}_{n}$ of $\Gamma_{n}$, defined in $[0,\ell_{\infty}]$. The new parametrization, by definition of $\mathcal{R}_{\eta}^{0,1}$, belongs to $W^{2,\infty}([0,\ell_{\infty}];\mathbb{R}^{2})$ and its norm is bounded above by a constant independent of $n$. From the Hausdorff convergence of $\Gamma_{n}$ to $\Gamma_{\infty}$, we deduce that $\hat{\gamma}_{n}$ converges to $\gamma_{\infty}$ weakly* in $W^{2,\infty}([0,\ell_{\infty}];\mathbb{R}^{2})$ and strongly in $W^{1,\infty}([0,\ell_{\infty}];\mathbb{R}^{2})$.

Let us fix $\rho>0$ sufficiently small, so that the projection $\Pi_{\Gamma_{\infty}}$ over $\Gamma_{\infty}$ is well defined in $\mathcal{I}_{\rho}({\Gamma}_{\infty}):=\{x\in\Omega:\,d(x,\Gamma_{\infty})<\rho)\}$. For $n$ large enough we have ${\Gamma}_{n}\subseteq\mathcal{I}_{\rho}({\Gamma}_{\infty})$. We want to construct a Lipschitz function $\Lambda_{n}$ such that $\Lambda_{n}(\Gamma_{\infty})=\Gamma_{n}$ and $\Lambda_{n}(x)=x$ for $x\in\mathbb{R}^{2}\setminus\mathcal{I}_{\rho}(\Gamma_{\infty})$. For every $x\in\mathcal{I}_{\rho}(\Gamma_{\infty})$ we define $s(x)\in[0,\ell_{\infty}]$ in such a way that $\gamma_{\infty}(s(x))=\Pi_{\Gamma_{\infty}}(x)$. We notice that the map $x\mapsto s(x)$ is locally Lipschitz, while $\Pi_{\Gamma_{\infty}}$ is Lipschitz on $\mathcal{I}_{\rho}(\Gamma_{\infty})$. Moreover, we set $d_{n}:=\|\hat{\gamma}_{n}-\gamma_{\infty}\|^{1/2}_{W^{1,\infty}}$ and \lambda_{n}(t):=\big{(}1-\tfrac{|t|}{d_{n}}\big{)}_{+}, where $(\cdot)_{+}$ stands for the positive part. With this notation at hand, we define

[TABLE]

In particular, $\Lambda_{n}$ is Lipschitz, $\|\Lambda_{n}-id\|_{W^{1,\infty}}\leq Cd_{n}\to 0$ as $n\to\infty$, and, for $n$ large enough, $\Lambda_{n}(\Gamma_{\infty})=\Gamma_{n}$ and $\Lambda_{n}=id$ out of $\mathcal{I}_{d_{n}}(\Gamma_{\infty})$. Applying the Hadamard Theorem [26, Theorem 6.2.3], we deduce that $\Lambda_{n}$ is globally invertible with $\|\Lambda_{n}^{-1}-id\|_{W^{1,\infty}}\to 0$ as $n\to\infty$.

Given $v\in H^{1}(\Omega\setminus\Gamma_{\infty};\mathbb{R}^{2})$ with $v=w_{\infty}$ on $\partial_{D}\Omega\setminus\Gamma_{0}$, we have that the function $v_{n}:=v\circ\Lambda_{n}^{-1}+w_{n}-w_{\infty}$ belongs to $H^{1}(\Omega\setminus\Gamma_{n};\mathbb{R}^{2})$ and satisfies $v_{n}=g_{n}$ on $\partial_{D}\Omega\setminus\Gamma_{0}$. Moreover, $\nabla{v_{n}}\to\nabla{v}$ in $L^{2}(\Omega;\mathbb{M}^{2})$, $v_{n}\to v$ in $L^{2}(\Omega;\mathbb{R}^{2})$, and, by the continuity of the trace operator, $v_{n}\to v$ strongly in $L^{2}(\partial_{S}\Omega;\mathbb{R}^{2})$. This asymptotic analysis implies that the sequence \big{(}E_{n}(v_{n})\big{)}_{n\in\mathbb{N}} is bounded and converges to $E_{\infty}(v)$ as $n\to\infty$.

By the minimality of $u_{n}$ for $E_{n}$, we have

[TABLE]

It is easy to see that the functionals $E_{n}$ are equi-coercive in $H^{1}(\Omega\setminus\Gamma_{n};\mathbb{R}^{2})$, so that inequality (2.6), together with Proposition 2.3, provides a uniform bound on the $L^{2}$ norm of $u_{n}$, of its gradient, and of its trace. Therefore, up to a subsequence (not relabeled), we have $u_{n}\rightharpoonup\varphi$ weakly in $L^{2}(\Omega;\mathbb{R}^{2})$ for some $\varphi\in L^{2}(\Omega;\mathbb{R}^{2})$. Moreover, in a suitably small neighborhood $\mathcal{U}$ of the boundary, this convergence is stronger, since $(\Omega\setminus\Gamma_{n})\cap\mathcal{U}=(\Omega\setminus\Gamma_{0})\cap\mathcal{U}$ for every $n$. More precisely, we have $u_{n}\rightharpoonup\varphi$ weakly in $H^{1}((\Omega\setminus\Gamma_{0})\cap\mathcal{U};\mathbb{R}^{2})$ and, therefore, $u_{n}\to\varphi$ strongly in $L^{2}(\partial\Omega;\mathbb{R}^{2})$ and $\varphi=w_{\infty}$ on $\partial_{D}\Omega$. The above convergences imply that

[TABLE]

Hence, passing to the liminf in (2.6) we get

[TABLE]

Thus, $\varphi$ is a minimizer of $E_{\infty}$ in $H^{1}(\Omega\setminus\Gamma_{\infty};\mathbb{R}^{2})$ with boundary condition $w_{\infty}$. Therefore, by uniqueness of the minimizer, $\varphi=u_{\infty}$. The strong convergence of the gradients follows by considering (2.6) for $v=u_{\infty}$. Indeed, we have

[TABLE]

which implies, together with (2.7), that $\mathcal{E}_{n}\to\mathcal{E}_{\infty}$ and $\mathrm{E}u_{n}\to\mathrm{E}u_{\infty}$ in $L^{2}(\Omega;\mathbb{M}^{2}_{sym})$. Applying Proposition 2.3 and recalling that $w_{n}\to w_{\infty}$ in $H^{1}(\Omega\setminus\Gamma_{0};\mathbb{R}^{2})$, we also obtain the strong convergence of $\nabla{u}_{n}$ to $\nabla{u}_{\infty}$ in $L^{2}(\Omega;\mathbb{M}^{2})$. This concludes the proof of the lemma. ∎

As a corollary of Lemma 2.4 we deduce the existence of solutions of the incremental minimum problems (1.4).

Corollary 2.5**.**

Fix $\varepsilon>0$, $k\in\mathbb{N}$, and $i\in\{1,\dots,k\}$. Then the minimum problem (1.4) admits a solution.

Proof.

It is sufficient to apply the direct method. Let $(\Gamma_{n})_{n\in\mathbb{N}}\subseteq\mathcal{R}_{\eta}$ be a minimizing sequence for (1.4). By Theorem 2.1, $\Gamma_{n}$ converges in the Hausdorff metric, up to a subsequence (not relabeled), to some $\Gamma_{\infty}\in\mathcal{R}_{\eta}$ such that the constraint $\Gamma_{\infty}\supseteq\Gamma_{\varepsilon,k,i-1}$ is preserved; moreover we have $\mathcal{H}^{1}(\Gamma_{n})\to\mathcal{H}^{1}(\Gamma_{\infty})$. Applying Lemma 2.4 with $\mathbb{C}_{n}=\mathbb{C}_{\infty}=\mathbb{C}$, $f_{n}=f_{\infty}=f(t_{k,i})$, $g_{n}=g_{\infty}=g(t_{k,i})$, and $w_{n}=w_{\infty}=w(t_{k,i})$, we obtain the convergence of the corresponding energies $\mathcal{E}(t_{k,i};\Gamma_{n})\to\mathcal{E}(t_{k,i};\Gamma_{\infty})$. Hence, $\Gamma_{\infty}$ is a solution to the minimum problem. ∎

3. The energy release rate

This section is devoted to the definition of the energy release rate, i.e., the opposite of the derivative of the energy $\mathcal{E}(t;\cdot)$ with respect to the crack elongation. The problem is clearly time-independent, therefore we omit the variable $t$, which is kept fixed. As in the previous section, the energy in (1.2) simply writes $\mathcal{E}(\Gamma)$.

Our aim is to generalize the results obtained in [2], where the energy release rate has been computed only in presence of smooth cracks $\Gamma$, in the absence of forces, and with a spatially constant elasticity tensor. Here we extend its definition to the case $\Gamma\in\mathcal{R}_{\eta}$, non-zero volume and boundary forces $f\in L^{2}(\Omega;\mathbb{R}^{2})$ and $g\in L^{2}(\partial_{S}\Omega;\mathbb{R}^{2})$, boundary datum $w\in H^{1}(\Omega\setminus\Gamma_{0};\mathbb{R}^{2})$, and non-constant tensor $\mathbb{C}\in C^{0,1}(\overline{\Omega})$.

As in [2], the fundamental steps are the following:

• $(i)$

Given an increasing family of cracks $\Gamma_{\sigma}\in\mathcal{R}_{\eta}^{0,1}$ parametrized by their arc length $\sigma\in[0,S]$, we prove that the map $\sigma\mapsto\mathcal{E}(\Gamma_{\sigma})$ is differentiable, thus

[TABLE]

• $(ii)$

We show that the above derivative only depends on the current configuration $\Gamma_{s}$, and not on its possible extensions, i.e., if $\Gamma_{\sigma}=\hat{\Gamma}_{\sigma}$ for $\sigma\leq s$, then

[TABLE]

In particular, we are allowed to write -\frac{\mathrm{d}\mathcal{E}(\Gamma_{\sigma})}{\mathrm{d}\sigma}\bigg{|}_{\sigma=s}=:\mathcal{G}(\Gamma_{s}) with no ambiguity.

We point out a difference of our strategy with respect to the proof of [27] for the antiplane case. In that case, the energy release rate is first characterized via the stress intensity factor assuming that the volume force is null in a neighborhood of the crack tip; then, one treats general forces by approximation, using the property that the stress intensity factor is continuous with respect to the force. In this paper, in the planar case we do not prove the existence of stress intensity factors for nonsmooth curves. Hence, when expressing the energy release rate via integral forms, we have to deal carefully with the terms containing the external force. Once the existence of the energy release rate is guaranteed, we will reduce to the case of forces that are null close to the tip via some approximation arguments, see Lemma 3.8 below.

In order to rigorously proceed with $(i)$, we first restrict our attention to cracks $\Gamma_{s}\in\mathcal{R}_{\eta}^{0,1}$. We write $\Gamma_{s}$ as

[TABLE]

where $\gamma\in C^{1,1}$ is the arc-length parametrization of $\Gamma_{s}$. We will discuss in Remark 3.10 how to tackle the general case $\Gamma\in\mathcal{R}_{\eta}$. For brevity, we denote with $u_{s}\in H^{1}(\Omega\setminus\Gamma_{s};\mathbb{R}^{2})$ the minimizer of (1.2). As in the previous section, when explicitly needed, we will highlight the dependence on the data by writing $\mathcal{E}(f,g,w,\mathbb{C};\Gamma_{s})$ for $\mathcal{E}(\Gamma_{s})$.

In order to make explicit computations, for every $s\in(0,S)$ and $\delta\in\mathbb{R}$ with $|\delta|$ small, we need to introduce a $C^{1,1}$ diffeomorphism $F_{s,\delta}$ that transforms $\Gamma_{s+\delta}$ in $\Gamma_{s}$ and maps $\Omega$ into itself. Precisely, for $r>0$ small enough we may assume that the curve $\Gamma_{\sigma}\cap\mathrm{B}_{r}(\gamma(s))$, for $|s-\sigma|$ small , is the graph of a $C^{1,1}$-function $\zeta$, with $\zeta^{\prime}(\gamma_{1}(s))=0$, where we have denoted with $\gamma_{1}$ the first component of the arc-length parametrization $\gamma$. For $\delta\in\mathbb{R}$ small in modulus , we define the function $F_{s,\delta}\colon\mathrm{B}_{r/2}(\gamma(s))\to\mathbb{R}^{2}$ by

[TABLE]

where $\varphi\in C^{\infty}_{c}(\mathrm{B}_{r/2}(\gamma(s)))$ is a suitable cut-off function equal to $1$ close to $\gamma(s)$. Notice that, for $r$ small enough, $\mathop{\rm spt}\nolimits(\varphi)\cap\mathop{\rm spt}\nolimits(w)=\emptyset$. We extend $F_{s,\delta}$ to the identity in $\mathbb{R}^{2}\setminus\mathrm{B}_{r/2}(\gamma(s))$.

By the regularity of $\zeta$, $F_{s,\delta}$ is a $C^{1,1}$ diffeomorphism of $\mathbb{R}^{2}$ such that $F_{s,\delta}(\Gamma_{s})=\Gamma_{s+\delta}$ and $F_{s,0}=id$. Moreover, the following equalities hold:

[TABLE]

With this notation at hand, we are in a position to prove the differentiability of $s\mapsto\mathcal{E}(\Gamma_{s})$.

Proposition 3.1**.**

Let $\{\Gamma_{\sigma}\}_{\sigma\geq 0}\subseteq\mathcal{R}^{0,1}_{\eta}$ be parametrized as in (3.1). Let $f\in L^{2}(\Omega;\mathbb{R}^{2})$, $g\in L^{2}(\partial_{S}\Omega;\mathbb{R}^{2})$, $w\in H^{1}(\Omega\setminus\Gamma_{0};\mathbb{R}^{2})$, and $\mathbb{C}\in C^{0,1}(\overline{\Omega})$. Then, the map $\sigma\mapsto\mathcal{E}(\Gamma_{\sigma})$ is differentiable and

[TABLE]

where $\mathrm{D}\mathbb{C}$ denotes the fourth order tensor

[TABLE]

Proof.

To prove (3.5), we compute explicitly the limits

[TABLE]

and show that the two limits coincide.

Let us start with (3.6). For every $\delta>0$, the function $u_{s}\circ F_{s,\delta}^{-1}$ belongs to $H^{1}(\Omega\setminus\Gamma_{s+\delta};\mathbb{R}^{2})$ and $u_{s}\circ F_{s,\delta}^{-1}=w$ on $\partial_{D}\Omega$. Hence,

[TABLE]

By a change of coordinate in the first integral in (3.8) we deduce that

[TABLE]

By a simple computation, we can rewrite (3.9) as

[TABLE]

Since

[TABLE]

where the limits are uniform in $\delta$, we obtain

[TABLE]

Applying e.g. [1, Lemma 3.8] (see also [21, Lemma 4.1]), we have that $\delta^{-1}(u_{s}{\,\circ\,}F_{s,\delta}^{-1}-u_{s})\to-\nabla{u_{s}}\,\rho_{s}$ in $L^{2}(\Omega)$ as $\delta\to 0$. Thus,

[TABLE]

Combining (3.10)-(3.14) we get

[TABLE]

In order to obtain the opposite inequality, we consider the function $u_{s+\delta}\circ F_{s,\delta}\in H^{1}(\Omega\setminus\Gamma_{s};\mathbb{R}^{2})$. By the minimality of $u_{s}$ we have that

[TABLE]

For simplicity of notation, we denote with $U_{s,\delta}\coloneqq u_{s+\delta}\circ F_{s,\delta}$. By a change of variable in the first integral in (3.17) and we deduce that

[TABLE]

Repeating the computations of (3.10)-(3.17) and taking into account that $\delta^{-1}(u_{s+\delta}-U_{s,\delta})\rightharpoonup-\nabla{u}\,\rho_{s}$ weakly in $L^{2}(\Omega;\mathbb{R}^{2})$ (see, e.g., [1, Lemma 3.8]), we infer that

[TABLE]

which, together with (3.17) implies that

[TABLE]

Adapting the above argument to the case $\delta<0$, cf. (3.7), it is also possible to prove that

[TABLE]

This concludes the proof of (3.5). ∎

The following corollary states the continuity of the derivative (3.5) w.r.t. the data $f$, $g$, $w$, $\mathbb{C}$, and $\Gamma_{s}$.

Corollary 3.2**.**

Let $f_{n},f\in L^{2}(\Omega;\mathbb{R}^{2})$, $g_{n},g\in L^{2}(\partial_{S}\Omega;\mathbb{R}^{2})$, $w_{n},w\in H^{1}(\Omega\setminus\Gamma_{0};\mathbb{R}^{2})$, and $\mathbb{C}_{n},\mathbb{C}\in C^{0,1}(\overline{\Omega})$ be such that $f_{n}\to f$ strongly in $L^{2}(\Omega;\mathbb{R}^{2})$, $g_{n}\rightharpoonup g$ weakly in $L^{2}(\partial_{S}\Omega;\mathbb{R}^{2})$, $w_{n}\rightharpoonup w$ weakly in $H^{1}(\Omega\setminus\Gamma_{0};\mathbb{R}^{2})$, and $\mathbb{C}_{n}\rightharpoonup\mathbb{C}$ weakly in $W^{1,\infty}(\Omega)$. Moreover, let $S>0$, let $\{\Gamma_{s}\}_{s\in[0,S]}\subseteq\mathcal{R}^{0,1}_{\eta}$ be as in (3.1), and assume that there exists a sequence $\{\Gamma^{n}_{s}\}_{s\in[0,S]}\subseteq\mathcal{R}^{0,1}_{\eta}$ such that $\Gamma^{n}_{s}$ converges to $\Gamma_{s}$ in the Hausdorff metric of sets for every $s\in[0,S]$. Then, for every $s\in(0,S)$ we have*

[TABLE]

Proof.

Let us denote with $u_{n}$ and $u$ the displacements associated to $\mathcal{E}(f_{n},g_{n},w_{n},\mathbb{C}_{n};\Gamma^{n}_{s})$ and to $\mathcal{E}(f,g,w,\mathbb{C};\Gamma_{s})$, respectively. By Lemma 2.4 and by the hypotheses, it follows that $\nabla u_{n}$ converges to $\nabla u$ strongly in $L^{2}(\Omega;\mathbb{M}^{2})$. Let us denote by $\rho^{n}_{s}$ the quantity defined as in (3.4) and corresponding to $\Gamma^{n}_{s}$. Since $\Gamma_{s}^{n}$ converges to $\Gamma_{s}$ in the Hausdorff metric of sets for every $s\in[0,S]$, we have that $\rho^{n}_{s}\to\rho_{s}$ uniformly in $\Omega$ and weakly* in $W^{1,\infty}(\Omega;\mathbb{R}^{2})$ for every $s\in(0,S)$. Thus (3.18) follows by (3.5). ∎

We notice that the dependence of \frac{\mathrm{d}\mathcal{E}(\Gamma_{\sigma})}{\mathrm{d}\sigma}\big{|}_{\sigma=s} on $\{\Gamma_{\sigma}\}_{\sigma\in[0,S]}$ is encoded by the quantity $\rho_{s}$ introduced in (3.4). The rest of this section is devoted to step (ii), namely at proving that the above derivative only depends on the current fracture $\Gamma_{s}$, and not on its possible extensions, i.e., on the choice of the family $\{\Gamma_{\sigma}\}_{\sigma\in[0,S]}$. We start by recalling a result of [2] stating that this very same property holds for $C^{\infty}$ cracks in absence of external forces and with constant elasticity tensor.

Theorem 3.3** ([2, Theorem 4.1]).**

Let $f=g=0$ and let $\mathbb{C}$ be constant in $\Omega$. Let $\{\Gamma_{\sigma}\}_{\sigma\in[0,S]}\subseteq\mathcal{R}^{0,1}_{\eta}$ be as in (3.1) and assume that there exists $s\in(0,S)$ such that $\Gamma_{\sigma}$ is of class $C^{\infty}$ for every $\sigma\in(0,\bar{s}]$. Then, for every $s\in(0,\bar{s}]$ there exist two constants $Q_{1}(\Gamma_{s})$ and $Q_{2}(\Gamma_{s})$ (independent of $\Gamma_{\sigma}$ for $\sigma>s$) such that

[TABLE]

where $C(\lambda,\mu)$ is a constant which depends only on the Lamé coefficients $\lambda$ and $\mu$.

Remark 3.4**.**

The constants $Q_{1}(\Gamma_{s})$ and $Q_{2}(\Gamma_{s})$ are the so called stress intensity factors. Indeed, it has been proven in [2, Theorem 2.5] that, in the condition of Theorem 3.3, the displacement $u_{s}$ can be written as

[TABLE]

for suitable functions $u_{R}\in H^{2}(\Omega\setminus\Gamma_{s};\mathbb{R}^{2})$ and $\Phi_{1},\Phi_{2}\in H^{1}(\Omega\setminus\Gamma_{s};\mathbb{R}^{2})\setminus H^{2}(\Omega\setminus\Gamma_{s};\mathbb{R}^{2})$. Moreover, the proof of formula (3.19) follows from the above decomposition.

The next proposition is a simple localization of Theorem 3.3.

Proposition 3.5**.**

Let $\{\Gamma_{s}\}_{s\in[0,S]}\subseteq\mathcal{R}^{0,1}_{\eta}$ be as in (3.1). Let $s\in(0,S)$, $f\in L^{2}(\Omega;\mathbb{R}^{2})$, $g\in L^{2}(\partial_{S}\Omega;\mathbb{R}^{2})$, $w\in H^{1}(\Omega\setminus\Gamma_{0};\mathbb{R}^{2})$, and $\mathbb{C}\in C^{0,1}(\overline{\Omega})$ be such that $\Gamma_{s}$ is $C^{\infty}$, $f=0$, and $\mathbb{C}$ is constant in a neighborhood of the tip $\gamma(s)$ of $\Gamma_{s}$. Then, there exist two constants $Q_{1}(\Gamma_{s})$ and $Q_{2}(\Gamma_{s})$ (independent of $\Gamma_{\sigma}$ for $\sigma>s$) such that

[TABLE]

where $C(\lambda_{s},\mu_{s})$ coincides with the constant appearing in (3.19) and $\lambda_{s},\mu_{s}$ denote the Lamé coefficients of $\mathbb{C}$ in $\gamma(s)$.

Proof.

As mentioned in Remark 3.4, the proof of formula (3.21) follows directly from a splitting of the form (3.20) for the displacement $u_{s}$ solution of

[TABLE]

close to the tip $\gamma(s)$ of $\Gamma_{s}$. Indeed, given (3.20), we can simply repeat step by step the proof of [2, Theorem 4.1] and get (3.21). In order to obtain such a decomposition in a neighborhood of $\gamma(s)$, we note that $u_{s}$ is also solution of

[TABLE]

with $\ell$ chosen in such a way that $\Gamma_{s}$ is smooth, $f=0$, and $\mathbb{C}$ is constant in $\mathrm{B}_{\ell}(\gamma(s))$. This enables us to apply [2, Theorem 2.5] on the domain $\mathrm{B}_{\ell}(\gamma(s))$ and to deduce the decomposition (3.21) on $\mathrm{B}_{\ell}(\gamma(s))$. ∎

We are now in a position to state and prove the main result of this section.

Theorem 3.6**.**

Let $\{\Gamma_{\sigma}\}_{\sigma\in[0,S]},\{\hat{\Gamma}_{\sigma}\}_{\sigma\in[0,S]}\subseteq\mathcal{R}^{0,1}_{\eta}$ be as in (3.1). Let $f\in L^{2}(\Omega;\mathbb{R}^{2})$, $g\in L^{2}(\partial_{S}\Omega;\mathbb{R}^{2})$, $w\in H^{1}(\Omega\setminus\Gamma_{0};\mathbb{R}^{2})$, and $\mathbb{C}\in C^{0,1}(\Omega)$. Let $s\in(0,S)$ and assume that $\Gamma_{\sigma}=\hat{\Gamma}_{\sigma}$ for $\sigma\leq s$. Then,

[TABLE]

Remark 3.7**.**

The previous theorem states that the derivative \frac{\mathrm{d}\mathcal{E}(\Gamma_{\sigma})}{\mathrm{d}\sigma}\big{|}_{\sigma=s} computed in Proposition 3.1 does not depend on the possible extension of $\Gamma_{s}$ in the class $\mathcal{R}^{0,1}_{\eta}$. Hence, it represents the slope of the energy $\mathcal{E}$ with respect to variations of crack in the set of admissible curves $\mathcal{R}^{0,1}_{\eta}$.

The proof of Theorem 3.6 is a corollary of the following lemma, where we use an approximation argument to reduce ourselves to the case of smooth cracks, constant elasticity tensor, and forces that are null close to the crack tip. In the latter case, the relation between the energy release rate and the stress intensity factors shows (3.22), cf. [2].

Lemma 3.8**.**

Let $\{\Gamma_{\sigma}\}_{\sigma\in[0,S]}$, $f$, and $\mathbb{C}$ be as in the statement of Theorem 3.6, and let $s\in(0,S)$. Then, there exist $\delta>0$, $f_{n}\in L^{2}(\Omega;\mathbb{R}^{2})$, $\mathbb{C}_{n}\in C^{0,1}(\overline{\Omega})$, and $\{\Gamma_{\sigma}^{n}\}_{\sigma\in[0,s+\delta]}\subseteq\mathcal{R}^{0,1}_{\eta}$ such that $f_{n}\to f$ strongly in $L^{2}(\Omega;\mathbb{R}^{2})$, $\mathbb{C}_{n}\rightharpoonup\mathbb{C}$ weakly in $W^{1,\infty}(\Omega)$, $\Gamma^{n}_{\sigma}\to\Gamma_{\sigma}$ in the Hausdorff metric of sets for every $\sigma\in[0,s+\delta]$, and, close to the tip of $\Gamma_{s}^{n}$, $\Gamma_{s}^{n}$ is smooth, $f_{n}=0$, and $\mathbb{C}_{n}$ is constant.*

Moreover, if $\{\hat{\Gamma}_{\sigma}\}_{\sigma\in[0,S]}$ is another family of curves in $\mathcal{R}_{\eta}^{0,1}$ with $\hat{\Gamma}_{\sigma}=\Gamma_{\sigma}$ for $\sigma\leq s$, then the sequences $\{\Gamma_{\sigma}^{n}\}_{\sigma\in[0,s+\delta]}$, $\{\hat{\Gamma}_{\sigma}^{n}\}_{\sigma\in[0,s+\delta]}$, $\mathbb{C}_{n}$, $\hat{\mathbb{C}}_{n}$, $f_{n}$, and $\hat{f}_{n}$ can be chosen in such a way that $\hat{\Gamma}_{\sigma}^{n}=\Gamma_{\sigma}^{n}$ for $\sigma\leq s$, $\mathbb{C}_{n}=\hat{\mathbb{C}}_{n}$, and $f_{n}=\hat{f}_{n}$.

Proof.

We start with the construction of an approximating sequence for $\Gamma_{\sigma}$. Let $\delta>0$ be such that $s+\delta<S$. Let us fix a sequence $s_{n}\nearrow s$. By definition of the class $\mathcal{R}_{\eta}^{0,1}$, for every $n$ there exist two open balls $\mathrm{B}_{\eta,n}^{1}$ and $\mathrm{B}_{\eta,n}^{2}$ of radius $\eta$ such that $\overline{\mathrm{B}}_{\eta,n}^{1}\cap\overline{\mathrm{B}}_{\eta,n}^{2}=\{\gamma(s_{n})\}$ and $(\mathrm{B}_{\eta,n}^{1}\cup\mathrm{B}_{\eta,n}^{2})\cap\Gamma_{s}=\emptyset$. Up to a redefinition of $\delta$, for $n$ large enough we may assume that the portion of curve $\{\gamma(\sigma):\,s_{n}-\delta\leq\sigma\leq s+\delta\}\subseteq\Gamma_{s+\delta}$ can be represented, in a suitable coordinate system $(x_{1},x_{2})$ possibly dependent on $n$, as graph of a function $\psi_{n}$ of class $C^{1,1}$ with $\psi^{\prime}_{n}(x_{1}^{s_{n}})=0$, where the point $(x_{1}^{s_{n}},\psi_{n}(x_{1}^{s_{n}}))$ coincides with $\gamma(s_{n})$. A similar notation is used for $\gamma(s)=(x_{1}^{s},\psi_{n}(x_{1}^{s}))$. Without loss of generality, we assume that $\psi_{n}^{\prime}(x_{1}^{s})\geq 0$.

The idea of our construction is to extend the curve $\Gamma_{s_{n}}$ with the arc of circumference of equation

[TABLE]

where $\bar{x}_{1}$ is the smallest $x_{1}\geq x_{1}^{s_{n}}$ such that $\psi_{n}^{\prime}(\bar{x}_{1})=\psi_{n}^{\prime}(x_{1}^{s})$. We notice that (3.23) is the equation of the boundary of one of the two open balls $\mathrm{B}_{\eta,n}^{i}$ and that $\bar{x}_{1}=x_{1}^{s_{n}}$ whenever $\psi_{n}^{\prime}(x_{1}^{s})=0$. We denote with $\Lambda^{n}$ the extension of $\Gamma_{s_{n}}$ with the arc (3.23) and its tip with $P_{n}$. We also use the symbol $\Lambda_{\sigma}^{n}$, $\sigma\in[0,\mathcal{H}^{1}(\Lambda^{n})]$, to indicate the piece of curve contained in $\Lambda^{n}$ of length $\sigma$.

A direct computation gives $\mathcal{H}^{1}(\Lambda^{n})=s_{n}+\eta\arctan(\psi_{n}^{\prime}(x_{1}^{s}))$, which can also be written as follows:

[TABLE]

On the other hand, exploiting the upper bound $\eta^{-1}$ on the curvature of the crack set, which reads $|\psi_{n}^{\prime\prime}|\,[1+(\psi_{n}^{\prime})^{2}]^{-3/2}\leq\eta^{-1}$ in terms of the graph parametrization, we get

[TABLE]

Comparing (3.24) and (3.25), we conclude that $\mathcal{H}^{1}(\Lambda^{n})\leq\mathcal{H}^{1}(\Gamma_{s})=s$.

If $\mathcal{H}^{1}(\Lambda^{n})<s$, we denote with $\alpha^{n}_{\ell}$ the segment of length $\ell$, initial point $P_{n}$ and parallel to $\gamma^{\prime}(s)$, and we define $\{\Gamma^{n}_{\sigma}\}_{\sigma\in[0,s+\delta]}$ as follows:

[TABLE]

where we have used the notation $v+E\coloneqq\{v+e:\>e\in E\}$ for $v\in\mathbb{R}^{2}$ and $E\subseteq\mathbb{R}^{2}$.

If $\mathcal{H}^{1}(\Lambda^{n})=s$, we simply set

[TABLE]

In both cases, we have that $\{\Gamma^{n}_{\sigma}\}_{\sigma\in[0,s+\delta]}\subseteq\mathcal{R}^{0,1}_{\eta}$, $\Gamma^{n}_{\sigma}\to\Gamma_{\sigma}$ in the Hausdorff metric of sets for every $\sigma\in[0,s+\delta]$, and $\Gamma^{n}_{s}$ is of class $C^{\infty}$ close to its tip.

The construction of $f_{n}$ is trivial, since the set of functions in $L^{2}(\Omega;\mathbb{R}^{2})$ that vanish close to $\gamma(s)$ is dense in $L^{2}(\Omega;\mathbb{R}^{2})$ w.r.t. the $L^{2}$-norm. We only have to ensure that $f_{n}$ is also null close to the tip of $\Gamma^{n}_{s}$, which is still possible because of the Hausdorff convergence.

As for the elasticity tensor $\mathbb{C}$, for every $r>0$ we consider a cut off function $\varphi_{r}$ in $\mathrm{B}_{r}(\gamma(s))$ with $\varphi_{r}(x)=\varphi_{r}(|x-\gamma(s)|)$, $\varphi_{r}=1$ in $\overline{\mathrm{B}}_{r/2}(\gamma(s))$, and $|\nabla\varphi_{r}|\leq C/r$ for some positive constant $C$ independent of $r$. Let us set $\mathbb{C}_{s}\coloneqq\mathbb{C}(\gamma(s))$ and $\mathbb{C}_{r}\coloneqq\varphi_{r}\mathbb{C}+(1-\varphi_{r})\mathbb{C}_{s}$. It is easy to see that $\mathbb{C}_{r}\in C^{0,1}(\overline{\Omega})$ with Lipschitz constant bounded by $C(\mathrm{Lip}(\mathbb{C})+1)$. Hence, $\mathbb{C}_{r}\rightharpoonup\mathbb{C}$ weakly* in $W^{1,\infty}(\Omega)$ as $r\searrow 0$. To conclude, it is enough to choose a suitable sequence $r_{n}\searrow 0$ in such a way that $\mathbb{C}_{n}\coloneqq\mathbb{C}_{r_{n}}$ is constant close to the tip of $\Gamma^{n}_{s}$. This is possible thanks to the Hausdorff convergence of $\Gamma^{n}_{s}$ to $\Gamma_{s}$.

The last part of the lemma is a trivial consequence of the above construction. ∎

Proof of Theorem 3.6.

To prove (3.22), we apply Lemma 3.8 to both $\Gamma_{\sigma}$ and $\hat{\Gamma}_{\sigma}$. Fixed $s\in(0,S)$ and $\delta>0$ small, let $\{\Gamma^{n}_{\sigma}\}_{\sigma\in[0,s+\delta]}$, $\{\hat{\Gamma}_{\sigma}\}_{\{\sigma\in[0,s+\delta]}$, $\mathbb{C}_{n}$, and $f_{n}$ be as in Lemma 3.8. By Corollary 3.2 we have that

[TABLE]

Taking into account Proposition 3.5, we have that

[TABLE]

and we deduce (3.22). ∎

We are now in a position to give the precise definition of energy release rate for a crack of the form (3.1). We stress that this is now possible thanks to Theorem 3.6.

Definition 3.9**.**

Let $\Gamma\in\mathcal{R}^{0,1}_{\eta}$, $s\coloneqq\mathcal{H}^{1}(\Gamma)$, $f\in L^{2}(\Omega;\mathbb{R}^{2})$, $g\in L^{2}(\partial_{S}\Omega;\mathbb{R}^{2})$, $w\in H^{1}(\Omega\setminus\Gamma_{0};\mathbb{R}^{2})$, and $\mathbb{C}\in C^{0,1}(\overline{\Omega})$. Let $S>s$ and let $\{\Gamma_{\sigma}\}_{\sigma\in[0,S]}\subseteq\mathcal{R}^{0,1}_{\eta}$ be such $\Gamma_{s}=\Gamma$. We define the energy release $\mathcal{G}(\Gamma)$ as

[TABLE]

Remark 3.10**.**

Definition 3.9, stated for a curve $\Gamma\in\mathcal{R}^{0,1}_{\eta}$, can be further generalized in order to consider general cracks in the class $\mathcal{R}_{\eta}$. Indeed, given $\Gamma\in\mathcal{R}_{\eta}$, it is enough to represent it as union of arcs of $C^{1,1}$ curves $\Gamma^{m}$, $m=1,\ldots,M$. In particular, each component belongs to $\mathcal{R}^{0,1}_{\eta}$ and can be written as in (3.1). Hence, for every $m$ we define the $m$-th energy release rate $\mathcal{G}^{m}(\Gamma)$ as in Definition 3.9 w.r.t. variations of the sole component $\Gamma^{m}$ of $\Gamma$. The energy release rate will be in this case the vector

[TABLE]

Remark 3.11**.**

We collect here the main properties of the energy release rate $\mathcal{G}$.

• $(a)$

$\mathcal{G}$ is continuous w.r.t. the Hausdorff convergence of cracks $\Gamma\in\mathcal{R}_{\eta}$, strong convergence of volume forces $f\in L^{2}(\Omega;\mathbb{R}^{2})$, weak convergence of surface forces $g\in L^{2}(\partial_{S}\Omega;\mathbb{R}^{2})$, and convergence of Dirichlet boundary data $w\in H^{1}(\Omega\setminus\Gamma_{0};\mathbb{R}^{2})$;

• $(b)$

there exists a positive constant $C=C(\mathbb{C},\eta)$ such that for every $\Gamma\in\mathcal{R}_{\eta}$, every $f\in L^{2}(\Omega;\mathbb{R}^{2})$, every $g\in L^{2}(\partial_{S}\Omega;\mathbb{R}^{2})$, and every $w\in H^{1}(\Omega\setminus\Gamma_{0};\mathbb{R}^{2})$,

[TABLE]

where $u\in H^{1}(\Omega\setminus\Gamma;\mathbb{R}^{2})$ is the solution of (1.2) with data $\Gamma$, $f$, $g$, $w$, and $\mathbb{C}$.

We will make use of these two properties in the proofs of Proposition 1.5 and Theorems 1.6 and 1.7.

4. Vanishing viscosity evolutions

In this section we focus on the proofs of existence of a viscous evolution $\Gamma_{\varepsilon}$ (see Proposition 1.5) and of a balanced viscosity evolution $\Gamma$ (Theorem 1.6), the latter obtained as limit of $\Gamma_{\varepsilon}$ as the viscosity parameter $\varepsilon$ tends to 0. To this end, we follow the method employed in a wide literature on rate-independent processes [32]. However, we point out that the abstract results of [29, 30, 31] do not directly apply to our setting.

Since the problem we analyze depends explicitly on time $t$ through the applied loads $f$, $g$, and $w$, from now on we denote with $\mathcal{G}(t;\Gamma)$ the energy release rate defined in Definition 3.9 and Remark 3.10 for a crack $\Gamma\in\mathcal{R}_{\eta}$ at time $t\in[0,T]$.

As anticipated in Section 1, the proofs of Proposition 1.5 and of Theorem 1.6 are based on a time-discretization procedure. Let the initial crack $\Gamma_{0}\in\mathcal{R}_{\eta}^{0}$ and the viscosity parameter $\varepsilon>0$ be fixed, and let us set $l^{m}_{0}\coloneqq\mathcal{H}^{1}(\Gamma^{m}_{0})$, where $\Gamma_{0}=\bigcup_{m=1}^{M}\Gamma^{m}_{0}$, according to Definition 1.1. For every $k\in\mathbb{N}$ we fix a partition $\{t_{k,i}\}_{i=0}^{k}$ of the time interval $[0,T]$ as in (1.3). For $i=0$ we set $\Gamma_{\varepsilon,k,0}\coloneqq\Gamma_{0}$. For $i\in\{1,\ldots,k\}$ we denote with $\Gamma_{\varepsilon,k,i}$ a minimizer of the incremental minimum problem (1.4), whose existence is provided by Corollary 2.5. Recalling the conventions of Definition 1.2, we write $\Gamma_{\varepsilon,k,i}=\bigcup_{m=1}^{M}\Gamma_{\varepsilon,k,i}^{m}$, we set $l^{m}_{\varepsilon,k,i}\coloneqq\mathcal{H}^{1}(\Gamma^{m}_{\varepsilon,k,i})$, and we denote with $P^{m}_{\varepsilon,k,i}$ the tip of $\Gamma^{m}_{\varepsilon,k,i}$. Furthermore, we define the interpolation functions

[TABLE]

where we denoted with $u^{\varepsilon}_{k,i}$ the function $u(t_{k,i};\Gamma_{\varepsilon,k,i})\in H^{1}(\Omega\setminus\Gamma_{\varepsilon,k,i})$.

In the following proposition we state a time discrete version of the Griffith’s criterion (G1)ε–(G3)ε.

Proposition 4.1**.**

For every $\varepsilon>0$, every $k\in\mathbb{N}$, every $m\in\{1,\ldots,M\}$, and a.e. $t\in(0,T]$ it holds:

• (G1)k

$\dot{l}^{m}_{\varepsilon,k}(t)\geq 0$;

• (G2)k

$\kappa(P^{m}_{\varepsilon,k}(t))-\mathcal{G}^{m}_{\varepsilon,k}(t)+\varepsilon\dot{l}^{m}_{\varepsilon,k}(t)\geq 0$;

• (G3)k

$\dot{l}^{m}_{\varepsilon,k}(t)(\kappa(P^{m}_{\varepsilon,k}(t))-\mathcal{G}^{m}_{\varepsilon,k}(t)+\varepsilon\dot{l}_{\varepsilon,k}(t))=0$.

Proof.

By construction, $l^{m}_{\varepsilon,k}$ is a non-decreasing function, so that (G1)k is clearly satisfied. In order to show (G2)k–(G3)k we take into account the minimality of $\Gamma_{\varepsilon,k,i}$. Let us fix $i\in\{1,\ldots\,k\}$. For every $\overline{m}\in\{1,\ldots,M\}$, let $\Gamma_{\overline{m}}\in\mathcal{R}_{\eta}$ be such that $\Gamma_{\overline{m}}=\bigcup_{m\neq\overline{m}}\Gamma_{\varepsilon,k,i}^{m}\cup\Lambda^{\overline{m}}$ with $\Lambda^{\overline{m}}\supseteq\Gamma^{\overline{m}}_{\varepsilon,k,i}$, and let us set $\lambda\coloneqq\mathcal{H}^{1}(\Lambda^{\overline{m}})\geq l^{\overline{m}}_{\varepsilon,k,i}$. Then,

[TABLE]

which implies

[TABLE]

We divide (4.1) by $\lambda-l^{\overline{m}}_{\varepsilon,k,i}$ and pass to the limit as $\lambda\to l^{\overline{m}}_{\varepsilon,k,i}$, obtaining (G2)k as a consequence of (H4) and of Definition 3.9. If, moreover, $\Gamma^{\overline{m}}_{\varepsilon,k,i-1}\subsetneq\Gamma^{\overline{m}}_{\varepsilon,k,i}$, we can consider as a competitor a set $\Gamma_{\overline{m}}\in\mathcal{R}_{\eta}$ as above, with $\Gamma^{\overline{m}}_{\varepsilon,k,i-1}\subseteq\Lambda^{\overline{m}}\subseteq\Gamma^{\overline{m}}_{\varepsilon,k,i}$, so that $l^{\overline{m}}_{\varepsilon,k,i-1}\leq\lambda\leq l^{\overline{m}}_{\varepsilon,k,i}$. Repeating the above computation we obtain

[TABLE]

This concludes the proof of (G3)k. ∎

We now show an a priori bound on $l_{\varepsilon,k}$ and on $u_{\varepsilon,k}$.

Proposition 4.2**.**

The following facts hold:

• $(a)$

there exist two positive constants $c$ and $C$ independent of $\varepsilon$, $k$, and $i$ such that for every $\varepsilon>0$, every $k\in\mathbb{N}$, and every $t\in[0,T]$,

[TABLE]

• $(b)$

for every $\varepsilon>0$, along a suitable (not relabeled) subsequence, $\|u_{\varepsilon,k}(t)\|_{2}$ and $\|\nabla{u}_{\varepsilon,k}(t)\|_{2}$ are bounded uniformly w.r.t. $t\in[0,T]$ and $k\in\mathbb{N}$;

• $(c)$

for every $\varepsilon>0$, along a suitable (not relabeled) subsequence, $\varepsilon\|\dot{l}^{m}_{\varepsilon,k}\|^{2}_{2}$ is bounded uniformly w.r.t. $k\in\mathbb{N}$ and $m\in\{1,\ldots,M\}$.

Proof.

For the sake of simplicity, let us denote with $w_{k,i}$, $f_{k,i}$, and $g_{k,i}$ the functions $w(t_{k,i})$, $f(t_{k,i})$, and $g(t_{k,i})$, respectively.

By definition of $u_{k,i}^{\varepsilon}$, by hypothesis (H3), and by the regularity of the data of the problem $f$, $g$, and $w$, we have that

[TABLE]

for some positive constants $C_{1}$ and $C_{2}$ depending only on $f$, $g$, $w$, and $\mathbb{C}$.

Since, for every $\varepsilon$ and $k$, the set function $\Gamma_{\varepsilon,k}\colon[0,T]\to\mathcal{R}_{\eta}$ is nondecreasing, we have that $u_{k,i}^{\varepsilon}\in H^{1}(\Omega\setminus\Gamma_{\varepsilon,k}(T);\mathbb{R}^{2})$. By definition of the class $\mathcal{R}_{\eta}$, the curves $\Gamma_{\varepsilon,k}(T)$ have bounded length uniformly w.r.t. $\varepsilon$ and $k$. Hence, we may assume that, up to a not relabeled subsequence, $\Gamma_{\varepsilon,k}(T)\to\widehat{\Gamma}_{\varepsilon}\in\mathcal{R}_{\eta}$ in the Hausdorff metric of sets as $k\to\infty$. We are therefore in a position to apply Proposition 2.3 to $\Gamma_{\varepsilon,k}(T)$, $\widehat{\Gamma}_{\varepsilon}$, and $u^{\varepsilon}_{k,i}$, which, together with (4.3), implies (b).

By definition of $\Gamma_{\varepsilon,k,i}$ and of the energy $\mathcal{E}(t_{k,i};\Gamma)$ we have that

[TABLE]

Iterating the above chain of inequalities for $i\in\{1,\ldots,k\}$ and using (H2) we deduce (4.2), which, together with (b), implies (c). ∎

In the following proposition we discuss the properties of the limit of the sequence $\Gamma_{\varepsilon,k}$ as $k\to\infty$.

Proposition 4.3**.**

For every $\varepsilon>0$ there exists a subsequence (not relabeled) of $\Gamma_{\varepsilon,k}$ and a set function $t\mapsto\Gamma_{\varepsilon}(t)\in\mathcal{R}_{\eta}$ such that $\Gamma^{m}_{\varepsilon,k}(t)$ converges to $\Gamma^{m}_{\varepsilon}(t)$ in the Hausdorff metric for every $t\in[0,T]$ and every $m\in\{1,\ldots,M\}$, and

• $(a)$

$\Gamma_{\varepsilon}$* is nondecreasing in time;*

• $(b)$

$l^{m}_{\varepsilon,k}\rightharpoonup l^{m}_{\varepsilon}$* weakly in $H^{1}(0,T)$ and $l^{m}_{\varepsilon,k}(t)\to l^{m}_{\varepsilon}(t)$ for every $t\in[0,T]$ and every $m$, where $l^{m}_{\varepsilon}(t)\coloneqq\mathcal{H}^{1}(\Gamma^{m}_{\varepsilon}(t))$;*

• $(c)$

$\nabla u_{\varepsilon,k}(t)\to\nabla u_{\varepsilon}(t)$* strongly in $L^{2}(\Omega;\mathbb{M}^{2})$ for every $t\in[0,T]$, where $u_{\varepsilon}(t)\coloneqq u(t;\Gamma_{\varepsilon}(t))$;*

• $(d)$

$\mathcal{G}_{\varepsilon,k}(t)\to\mathcal{G}_{\varepsilon}(t)$* for every $t\in[0,T]$, where $\mathcal{G}_{\varepsilon}(t)\coloneqq\mathcal{G}(t;\Gamma_{\varepsilon}(t))$;*

• $(e)$

$\mathcal{G}_{\varepsilon,k}\to\mathcal{G}_{\varepsilon}$* in $L^{2}(0,T)$.*

Moreover, along a suitable (not relabeled) subsequence, we have

• $(f)$

$\varepsilon\|{\dot{l}^{m}_{\varepsilon}}\|_{2}^{2}$* is uniformly bounded in $\varepsilon$ for every $m\in\{1,\ldots,M\}$;*

• $(g)$

$\|{\nabla u_{\varepsilon}(t)}\|_{2}$* is uniformly bounded in $\varepsilon$ and $t$. *

Proof.

For brevity, in the following we will not relabel subsequences. For $\varepsilon>0$ let us consider the subsequence $\Gamma_{\varepsilon,k}$ detected in (b) and (c) of Proposition 4.2. Since $\Gamma_{\varepsilon,k}$ is a sequence of increasing set functions with uniformly bounded length, there exists a nondecreasing set function $\Gamma_{\varepsilon}\colon[0,T]\to\mathcal{R}_{\eta}$ such that, up to a further subsequence, $\Gamma_{\varepsilon,k}(t)$ converges to $\Gamma_{\varepsilon}(t)$ in the Hausdorff metric of sets for every $t\in[0,T]$. Hence, for every $m\in\{1,\ldots,M\}$ and every $t\in[0,T]$ it holds $\Gamma^{m}_{\varepsilon,k}(t)\to\Gamma^{m}_{\varepsilon}(t)$.

For $\varepsilon>0$ fixed, by Proposition 4.2 we have that $l^{m}_{\varepsilon,k}\in H^{1}(0,T)$ is bounded w.r.t. $k$ and $m$. Therefore, for every $m\in\{1,\ldots,M\}$ the sequence $l^{m}_{\varepsilon,k}$ converges weakly in $H^{1}(0,T)$ to a nondecreasing function $l^{m}_{\varepsilon}$. Up to a further subsequence, we may assume that $l^{m}_{\varepsilon,k}(t)\to l^{m}_{\varepsilon}(t)$ for every $t\in[0,T]$ and every $m$. In particular, $l^{m}_{\varepsilon}(t)=\mathcal{H}^{1}(\Gamma^{m}_{\varepsilon}(t))$, so that (b) is proven. We also notice that, because of the continuity of $l^{m}_{\varepsilon}$, we have that $\underline{\Gamma}_{\varepsilon,k}(t)\to\Gamma_{\varepsilon}(t)$ in the Hausdorff metric as $k\to\infty$.

The $L^{2}$-convergence of $\nabla u_{\varepsilon,k}(t)$ to $\nabla u_{\varepsilon}(t)$ is a consequence of the convergence of $\Gamma_{\varepsilon,k}(t)$ to $\Gamma_{\varepsilon}(t)$ and of Lemma 2.4. In a similar way, since $\underline{\Gamma}_{\varepsilon,k}(t)$ converges to $\Gamma_{\varepsilon}(t)$, $\nabla\underline{u}_{\varepsilon,k}(t)\to\nabla u_{\varepsilon}(t)$ in $L^{2}(\Omega;\mathbb{M}^{2})$. Moreover, by Remark 3.11 we have that $\mathcal{G}_{\varepsilon,k}(t)\to\mathcal{G}(t;\Gamma_{\varepsilon}(t))=:\mathcal{G}_{\varepsilon}(t)$ for every $t\in[0,T]$, so that (d) holds. Being $\|\nabla{u}_{\varepsilon,k}(t)\|_{2}$ and $\|u_{\varepsilon,k}(t)\|_{2}$ bounded uniformly w.r.t. $t$ and $k$, again by Remark 3.11 we infer that $\mathcal{G}_{\varepsilon,k}(t)$ is bounded, so that $\mathcal{G}_{\varepsilon,k}\to\mathcal{G}_{\varepsilon}$ in $L^{2}(0,T)$ and (e) is concluded.

In order to prove (f) and (g), we employ Proposition 4.2, obtaining

[TABLE]

where $C>0$ is independent of $t$ and $\varepsilon$. Arguing as in the proof of Proposition 4.2, we have that $u_{\varepsilon}(t)\in H^{1}(\Omega\setminus\Gamma_{\varepsilon}(T);\mathbb{R}^{2})$ for every $t\in[0,T]$. Since $\Gamma_{\varepsilon}(T)\in\mathcal{R}_{\eta}$ has a uniformly bounded length, we may assume that, up to a subsequence, $\Gamma_{\varepsilon}(T)\to\widehat{\Gamma}\in\mathcal{R}_{\eta}$ in the Hausdorff metric of sets. Thus, we can apply Proposition 2.3 to $\Gamma_{\varepsilon}(T)$, $\widehat{\Gamma}$, and $u_{\varepsilon}(t)$, to deduce from (4.4) that $\|\nabla{u}_{\varepsilon}(t)\|_{2}$ is bounded uniformly w.r.t. $\varepsilon$ and $t$, so that (g) holds.

Finally, we pass to the liminf in (4.2) for $t=T$, obtaining

[TABLE]

By the boundedness of $\|\nabla{u}_{\varepsilon}(t)\|_{2}$ and of $\|u_{\varepsilon}(t)\|_{2}$ we immediately get (f), and the proof is thus concluded. ∎

We are now in a position to prove Proposition 1.5.

Proof of Proposition 1.5.

Let $\Gamma_{\varepsilon}$, $l^{m}_{\varepsilon}$, and $\mathcal{G}^{m}_{\varepsilon}$ be the functions determined in Proposition 4.3. Since $l^{m}_{\varepsilon}$ is nondecreasing, (G1)ε is satisfied. In order to prove (G2)ε let us consider $\psi\in L^{2}(0,T)$ with $\psi\geq 0$. By (G2)k we have

[TABLE]

From the Hausdorff convergence of $\Gamma^{m}_{\varepsilon,k}(t)$ to $\Gamma^{m}_{\varepsilon}(t)$ it follows that $P^{m}_{\varepsilon,k}(t)\to P^{m}_{\varepsilon}(t)$ for every $t\in[0,T]$ and every $m$, where $P^{m}_{\varepsilon}(t)$ stands for the tip of $\Gamma^{m}_{\varepsilon}(t)$. By hypothesis (H4) we have that $\kappa(P^{m}_{\varepsilon,k})\to\kappa(P^{m}_{\varepsilon})$ in $L^{2}(0,T)$. Hence, passing to the limit in (4.5) as $k\to\infty$ and taking into account (e) of Proposition 4.3 we get

[TABLE]

By the arbitrariness of $\psi\in L^{2}(0,T)$, $\psi\geq 0$, we infer (G2)ε.

As for (G3)ε, we integrate (G3)k over $[0,T]$ and pass to the liminf as $k\to\infty$. By (b) and (e) of Proposition 4.3 and by the convergence of $\kappa(P^{m}_{\varepsilon,k})$ to $\kappa(P^{m}_{\varepsilon})$ we obtain

[TABLE]

Combining the previous inequality with (G1)ε and (G2)ε we deduce (G3)ε. Finally, the uniform boundedness of $\varepsilon\|\dot{l}^{m}_{\varepsilon}\|_{2}^{2}$ has been stated in (f) of Proposition 4.3. ∎

Remark 4.4**.**

Let $\Gamma_{\varepsilon}$ be as in Proposition 1.5. Then, for every $t\in[0,T]$ it holds

[TABLE]

Indeed, being $l_{\varepsilon}\in H^{1}([0,T];\mathbb{R}^{M})$, the function $t\mapsto\mathcal{F}(t,\Gamma_{\varepsilon}(t))=\mathcal{E}(t;\Gamma_{\varepsilon}(t))+\mathcal{K}(\Gamma_{\varepsilon}(t))$ belongs to $H^{1}(0,T)$ with

[TABLE]

We conclude with the proof of Theorem 1.6.

Proof of Theorem 1.6.

For $\varepsilon>0$ and $m\in\{1,\ldots,M\}$ let $\Gamma_{\varepsilon}$, $\Gamma^{m}_{\varepsilon}$, and $l^{m}_{\varepsilon}$ be the viscous evolutions determined in Proposition 1.5. Let us consider, without relabeling, the $\varepsilon$-subsequence satisfying (f) and (g) of Proposition 4.3. Since $\Gamma_{\varepsilon}$ is a sequence of nondecreasing set functions and $\mathcal{H}^{1}(\Gamma_{\varepsilon}(t))$ is uniformly bounded w.r.t. $t\in[0,T]$ and $\varepsilon>0$, there exists a nondecreasing set function $\Gamma\colon[0,T]\to\mathcal{R}_{\eta}$ such that $\Gamma_{\varepsilon}(t)\to\Gamma(t)$ in the Hausdorff metric of sets for every $t\in[0,T]$. In particular, $\Gamma^{m}_{\varepsilon}(t)\to\Gamma^{m}(t)$ for every $t$ and every $m\in\{1,\ldots,M\}$, where $\Gamma(t)=\bigcup_{m=1}^{M}\Gamma^{m}(t)$. Moreover, being $l^{m}_{\varepsilon}$ a sequence of bounded nondecreasing functions, we may assume that, up to a further subsequence, $l^{m}_{\varepsilon}(t)\to l^{m}(t)$ for every $t\in[0,T]$ and $l^{m}_{\varepsilon}\to l^{m}$ in $L^{2}(0,T)$. In particular, $l^{m}(t)=\mathcal{H}^{1}(\Gamma^{m}(t))$ and (G1) is proven.

In order to show (G2), let us consider $\psi\in L^{2}(0,T)$ with $\psi\geq 0$. In view of (G2)ε we have

[TABLE]

Since $\Gamma^{m}_{\varepsilon}(t)\to\Gamma^{m}(t)$, we have that $P^{m}_{\varepsilon}(t)\to P^{m}(t)$ for every $t$ and every $m$, where $P^{m}(t)$ is the tip of $\Gamma^{m}(t)$. Thus, by hypothesis (H4) we get that $\kappa(P^{m}_{\varepsilon})\to\kappa(P^{m})$ in $L^{2}(0,T)$ for every $m$. From (e) and (f) of Proposition 4.3 we deduce that $\varepsilon\dot{l}^{m}_{\varepsilon}\to 0$ and $\mathcal{G}^{m}_{\varepsilon}\to\mathcal{G}^{m}$ in $L^{2}(0,T)$. Hence, passing to the limit in (4.7) we get

[TABLE]

As a consequence, $\kappa(P^{m}(t))-\mathcal{G}^{m}(t)\geq 0$ for a.e. $t\in[0,T]$. By continuity, this inequality holds in all the continuity points of $\Gamma^{m}(t)$. Hence, (G2) is proven.

As for (G3), we integrate (G3)ε over the interval $[0,T]$ and notice that the term $\varepsilon(\dot{l}^{m}_{\varepsilon})^{2}$ is positive, so that

[TABLE]

Passing to the limit in the previous inequality we get

[TABLE]

Combining (4.8) with (G1) and (G2) we deduce (G3). ∎

5. Parametrized evolutions

This section is devoted to the proof of Theorem 1.7. The strategy is to perform a change of variables which transforms the lengths $l^{m}_{\varepsilon}$ obtained in Proposition 1.5 in absolutely continuous functions. Roughly speaking, this is done by a parametrization of time on the jump points of the viscous solution $l^{m}_{\varepsilon}$.

Let us fix the sequence $\varepsilon\to 0$ determined in Proposition 4.3 and Theorem 1.6. For $t\in[0,T]$ we set

[TABLE]

Thanks to the properties of $l^{m}_{\varepsilon}$ (see Proposition 1.5), $\sigma_{\varepsilon}$ is strictly increasing, continuous, and $\dot{\sigma}_{\varepsilon}(t)\geq 1$ for every $\varepsilon>0$ and a.e. $t\in[0,T]$. Therefore, $\sigma_{\varepsilon}$ is invertible and we denote its inverse with $\tilde{t}_{\varepsilon}\colon[0,\sigma_{\varepsilon}(T)]\to[0,T]$. We deduce that $\tilde{t}_{\varepsilon}$ is strictly increasing, continuous, and $0<\tilde{t}^{\prime}_{\varepsilon}(\sigma)\leq 1$ for every $\varepsilon>0$ and a.e. $\sigma\in[0,\sigma_{\varepsilon}(T)]$, where the symbol ′ stands for the derivative with respect to $\sigma$.

For $m=1,\ldots,M$ and $\sigma\in[0,\sigma_{\varepsilon}(T)]$, we set

[TABLE]

By (5.1) we have $\sigma=\tilde{t}_{\varepsilon}(\sigma)+|\tilde{l}_{\varepsilon}(\sigma)|_{1}-|l_{0}|_{1}$. Differentiating this relation we get

[TABLE]

for every $\varepsilon>0$ and a.e. $\sigma\in[0,\sigma_{\varepsilon}(T)]$. By (5.2) and the monotonicity of $\tilde{l}_{\varepsilon}^{m}$ we have $0\leq(\tilde{l}_{\varepsilon}^{m})^{\prime}(\sigma)\leq 1$ for every $\varepsilon>0$, every $m=1,\ldots,M$, and a.e. $\sigma\in[0,\sigma_{\varepsilon}(T)]$. Moreover, $\tilde{t}_{\varepsilon}$ and $\tilde{l}_{\varepsilon}$ are Lipschitz continuous.

We define $\widetilde{\mathcal{G}}_{\varepsilon}^{m}(\sigma)\coloneqq\mathcal{G}^{m}(\tilde{t}_{\varepsilon}(\sigma);\widetilde{\Gamma}_{\varepsilon}(\sigma))$ for $\sigma\in[0,\sigma_{\varepsilon}(T)]$ and $\bar{S}\coloneqq\sup_{\varepsilon>0}\sigma_{\varepsilon}(T)$, which is bounded by a constant depending on $T$ and on the class $\mathcal{R}_{\eta}$. In order to deal with functions defined on the same interval, we extend $\tilde{t}_{\varepsilon}$, $\tilde{l}_{\varepsilon}$, $\widetilde{\Gamma}_{\varepsilon}$, $\widetilde{\Gamma}^{m}_{\varepsilon}$, $\tilde{t}_{\varepsilon}^{\prime}$, and $\tilde{l}_{\varepsilon}^{\prime}$ on $(\sigma_{\varepsilon}(T),\bar{S}]$ by $\tilde{t}_{\varepsilon}(\sigma)\coloneqq\tilde{t}_{\varepsilon}(\sigma_{\varepsilon}(T))$, $\tilde{l}_{\varepsilon}(\sigma)\coloneqq\tilde{l}_{\varepsilon}(\sigma_{\varepsilon}(T))$, $\widetilde{\Gamma}_{\varepsilon}(\sigma)\coloneqq\widetilde{\Gamma}_{\varepsilon}(\sigma_{\varepsilon}(T))$, $\widetilde{\Gamma}^{m}_{\varepsilon}(\sigma)\coloneqq\widetilde{\Gamma}^{m}_{\varepsilon}(\sigma_{\varepsilon}(T))$, $\tilde{t}_{\varepsilon}^{\prime}(\sigma)\coloneqq 0$, and $\tilde{s}_{\varepsilon}^{\prime}(\sigma)\coloneqq 0$.

Recalling that $\tilde{t}_{\varepsilon}^{\prime}(\sigma)>0$ on $[0,\sigma_{\varepsilon}(T)]$, the Griffith’s criterion stated in Proposition 1.5 reads in the new variables as

[TABLE]

for every $m$, every $\varepsilon$, and a.e. $\sigma\in[0,\bar{S}]$.

Finally, we observe that by (f) of Proposition 4.3

[TABLE]

uniformly in $\varepsilon$ and $m\in\{1,\ldots,M\}$. Therefore, $\varepsilon(\tilde{l}^{m}_{\varepsilon})^{\prime}\to 0$ in $L^{2}(0,\bar{S})$.

Passing to the limit as $\varepsilon\to 0$, we are now able to prove Theorem 1.7, showing that the parametrized solution $\widetilde{\Gamma}$ satisfies a generalized Griffith’s criterion.

Proof of Theorem 1.7.

Since $\widetilde{\Gamma}_{\varepsilon}\colon[0,\bar{S}]\to\mathcal{R}_{\eta}$ is a nondecreasing set function with uniformly bounded length, there exists $\widetilde{\Gamma}\colon[0,\bar{S}]\to\mathcal{R}_{\eta}$ such that, up to a not relabeled subsequence, $\widetilde{\Gamma}_{\varepsilon}(\sigma)\to\widetilde{\Gamma}(\sigma)$ and $\widetilde{\Gamma}^{m}_{\varepsilon}(\sigma)\to\widetilde{\Gamma}^{m}(\sigma)$ in the Hausdorff metric of sets for every $\sigma\in[0,\bar{S}]$ and every $m\in\{1,\ldots,M\}$. We denote with $\widetilde{P}^{m}(\sigma)$ the tip of $\widetilde{\Gamma}^{m}(\sigma)$ and we notice that $\widetilde{P}^{m}_{\varepsilon}(\sigma)\to\widetilde{P}^{m}(\sigma)$ for every $\sigma\in[0,\bar{S}]$ and every $m\in\{1,\ldots,M\}$.

Being $\tilde{t}_{\varepsilon},\tilde{l}^{m}_{\varepsilon}$ bounded in $W^{1,\infty}(0,\bar{S})$, up to a further subsequence we have that $\tilde{t}_{\varepsilon}$ and $\tilde{l}^{m}_{\varepsilon}$ converge weakly* in $W^{1,\infty}(0,\bar{S})$ to some functions $\tilde{t}$ and $\tilde{l}^{m}$, respectively, As a consequence, we have that $\tilde{l}^{m}(\sigma)=\mathcal{H}^{1}(\widetilde{\Gamma}^{m}(\sigma))$, so that $\widetilde{\Gamma}^{m}\colon[0,S]\to\mathcal{R}_{\eta}$ is continuous in the Hausdorff metric of sets. We can also assume that $\sigma_{\varepsilon}(T)\to S$ and $\tilde{t},\tilde{l}^{m}\in W^{1,\infty}(0,S)$. Moreover, writing (5.2) in an integral form and passing to the limit, we deduce that for a.e. $\sigma\in[0,S]$

[TABLE]

where we have set $\tilde{l}(\sigma)\coloneqq(l^{1}(\sigma),\ldots,\tilde{l}^{M}(\sigma))$. For $m=1\ldots,M$ and $\sigma\in[0,S]$ we define,

[TABLE]

We notice that, by Remark 3.11, $\widetilde{\mathcal{G}}_{\varepsilon}(\sigma)$ converges to $\widetilde{\mathcal{G}}(\sigma)$ for every $\sigma\in[0,S]$ and $\widetilde{\mathcal{G}}_{\varepsilon}\to\widetilde{\mathcal{G}}$ in $L^{2}(0,S)$, as $\varepsilon\to 0$.

By the monotonicity of $\tilde{t}$ and $\tilde{l}^{m}$, we have $\tilde{t}^{\prime}(\sigma)\geq 0$ and $(\tilde{l}^{m})^{\prime}(\sigma)\geq 0$ for every $m$ and a.e. $\sigma\in[0,S]$. Moreover, by (5.4) they can not be simultaneously zero.

Let us fix $m\in\{1,\ldots,M\}$ and $\psi\in L^{2}(0,S)$ with $\psi\geq 0$. Thanks to (5.3), for every $\varepsilon$ we have

[TABLE]

Since $\tilde{t}^{\prime}_{\varepsilon}$ converges to $\tilde{t}^{\prime}$ weakly* in $L^{\infty}(0,S)$, $\varepsilon(\tilde{l}^{m}_{\varepsilon})^{\prime}\to 0$ in $L^{2}(0,S)$, $\widetilde{P}^{m}_{\varepsilon}(\sigma)\to\widetilde{P}^{m}(\sigma)$ for $\sigma\in[0,S]$, and $\widetilde{\mathcal{G}}^{m}_{\varepsilon}\to\widetilde{\mathcal{G}}^{m}$ in $L^{2}(0,S)$, passing to the limit in (5.5) as $\varepsilon\to 0$ we get

[TABLE]

which implies (pG2).

We notice that if (pG1), (pG2) and (5.4) hold, then (pG3) and (pG4) are equivalent to the following property:

[TABLE]

Let us therefore assume that $\widetilde{\mathcal{G}}^{m}(\bar{\sigma})<\kappa(\widetilde{P}^{m}(\bar{\sigma}))$. We first claim that there exist $\bar{\varepsilon}>0$ and $\delta>0$ such that $\widetilde{\mathcal{G}}^{m}_{\varepsilon}(\sigma)<\kappa(\widetilde{P}^{m}_{\varepsilon}(\sigma))$ for every $\sigma\in(\bar{\sigma}-\delta,\bar{\sigma}+\delta)$ and every $\varepsilon\leq\bar{\varepsilon}$. By contradiction, suppose that this is not the case. Then, there exist $\sigma_{k}\to\bar{\sigma}$ and $\varepsilon_{k}\to 0$ such that $\widetilde{\mathcal{G}}^{m}_{\varepsilon_{k}}(\sigma_{k})\geq\kappa(\widetilde{P}^{m}_{\varepsilon_{k}}(\sigma_{k}))$. By continuity and monotonicity of $\widetilde{\Gamma}^{m}$, we have that $\widetilde{\Gamma}^{m}_{\varepsilon_{k}}(\sigma_{k})\to\widetilde{\Gamma}^{m}(\sigma)$ in the Hausdorff metric of sets and $\widetilde{P}^{m}_{\varepsilon_{k}}(\sigma_{k})\to\widetilde{P}^{m}(\sigma)$ as $k\to\infty$. Hence, the continuity of the energy release rate and the hypothesis (H4) lead us to the contradiction $\widetilde{\mathcal{G}}^{m}(\bar{\sigma})\geq\kappa(\widetilde{P}^{m}(\bar{\sigma}))$.

Let $\delta$ and $\bar{\varepsilon}$ be as above. We deduce from the Griffith’s criterion (5.3) that $\tilde{l}^{m}_{\varepsilon}$ is constant in $(\bar{\sigma}-\delta,\bar{\sigma}+\delta)$ for every $\varepsilon\leq\bar{\varepsilon}$. Since $\tilde{l}^{m}_{\varepsilon}$ converges to $\tilde{l}^{m}$ weakly* in $W^{1,\infty}(0,S)$, we get that $\tilde{l}^{m}$ is locally constant around $\bar{\sigma}$, and this concludes the proof of (pG3) and (pG4).

In order to show that $\Gamma(\tilde{t}(\sigma))=\widetilde{\Gamma}(\sigma)$ for every $\sigma\in[0,S]$ such that $\tilde{t}^{\prime}(\sigma)>0$, we define

[TABLE]

If $\tilde{t}^{\prime}(\sigma)>0$, then we have $s(\tilde{t}(\sigma))=\sigma$ and $s(\tilde{t}(\bar{\sigma}))\neq s(\tilde{t}(\sigma))$ for $\bar{\sigma}\neq\sigma$. Let us prove that $\tilde{t}(\sigma)$ is a continuity point for $\Gamma$, where the map $t\mapsto\Gamma(t)$ has been determined in Theorem 1.6. By contradiction, assume that $\tilde{t}(\sigma)$ is a discontinuity point of $\Gamma$. Then, there exist $t_{\varepsilon}^{1}<t_{\varepsilon}^{2}$ such that $t_{\varepsilon}^{1},t_{\varepsilon}^{2}\to\tilde{t}(\sigma)$ and $\Gamma_{\varepsilon}(t^{1}_{\varepsilon})\to\Gamma^{-}(\tilde{t}(\sigma))$ and $\Gamma_{\varepsilon}(t^{2}_{\varepsilon})\to\Gamma^{+}(\tilde{t}(\sigma))$ in the Hausdorff metric of sets, where we have denoted with $\Gamma^{\pm}(\tilde{t}(\sigma))$ the left and right limits of $\Gamma(t)$ in $\tilde{t}(\sigma)$. As a consequence, $s(\tilde{t}(\bar{\sigma}))=s(\tilde{t}(\sigma))$ for $\bar{\sigma}\in(\sigma-\mathcal{H}^{1}(\Gamma^{+}(\tilde{t}(\sigma))\setminus\Gamma^{-}(\tilde{t}(\sigma))),\sigma+\mathcal{H}^{1}(\Gamma^{+}(\tilde{t}(\sigma))\setminus\Gamma^{-}(\tilde{t}(\sigma))))$, which is a contradiction. Hence, $\tilde{t}(\sigma)$ is a continuity point of $t\mapsto\Gamma(t)$. Therefore, $\widetilde{\Gamma}_{\varepsilon}(\sigma)=\Gamma_{\varepsilon}(\tilde{t}_{\varepsilon}(\sigma))$ converges to $\Gamma(\tilde{t}(\sigma))$ in the Hausdorff metric of sets. This implies that $\Gamma(\tilde{t}(\sigma))=\widetilde{\Gamma}(\sigma)$.

We conclude with the energy balance (1.6). Let $s\in[0,S]$. By the change of variable $t=t_{\varepsilon}(\sigma)$ in (4.6), for $\varepsilon>0$ we have

[TABLE]

where we have set $\tilde{u}_{\varepsilon}(\sigma)\coloneqq u_{\varepsilon}(\tilde{t}_{\varepsilon}(\sigma))$. Since $\tilde{t}_{\varepsilon}$ and $\tilde{l}^{m}_{\varepsilon}$ converge weakly* in $W^{1,\infty}(0,S)$ to $\tilde{t}$ and $\tilde{l}^{m}$ and $\widetilde{\mathcal{G}}^{m}_{\varepsilon}\to\widetilde{\mathcal{G}}^{m}$ in $L^{2}(0,S)$, passing to the limit as $\varepsilon\to 0$ in (5.6) we get (1.6). This concludes the proof of the theorem. ∎

Acknowledgments

The authors would like to acknowledge the kind hospitality of the Erwin Schrödinger International Institute for Mathematics and Physics (ESI), where part of this research was developed during the workshop New trends in the variational modeling of failure phenomena. All authors would like to acknowledge the kind hospitality of the University of Naples Federico II, to which GL was affiliated when this research was initiated. SA wishes to thank the Technical University of Munich, where he worked during the preparation of this paper, with partial support from the SFB project TRR109 Shearlet approximation of brittle fracture evolutions. GL and IL are members of the Gruppo Nazionale per l’Analisi Matematica, la Probabilità e le loro Applicazioni (GNAMPA) of the Istituto Nazionale di Alta Matematica (INdAM). GL received support from the INdAM-GNAMPA 2018 Project Analisi variazionale per difetti e interfacce nei materiali and through the 2019 Project Modellazione e studio di proprietà asintotiche per problemi variazionali in fenomeni anelastici.