Skip to content
BY 4.0 license Open Access Published by De Gruyter November 7, 2020

Integral representation for energies in linear elasticity with surface discontinuities

  • Vito Crismale ORCID logo , Manuel Friedrich and Francesco Solombrino ORCID logo EMAIL logo

Abstract

In this paper we prove an integral representation formula for a general class of energies defined on the space of generalized special functions of bounded deformation ( GSBD p ) in arbitrary space dimensions. Functionals of this type naturally arise in the modeling of linear elastic solids with surface discontinuities including phenomena as fracture, damage, surface tension between different elastic phases, or material voids. Our approach is based on the global method for relaxation devised in [G. Bouchitté, I. Fonseca and L. Mascarenhas, A global method for relaxation, Arch. Ration. Mech. Anal. 145 1998, 1, 51–98] and a recent Korn-type inequality in GSBD p , cf. [F. Cagnetti, A. Chambolle and L. Scardia, Korn and Poincaré–Korn inequalities for functions with a small jump set, preprint 2020]. Our general strategy also allows to generalize integral representation results in SBD p , obtained in dimension two [S. Conti, M. Focardi and F. Iurlano, Integral representation for functionals defined on SBD p in dimension two, Arch. Ration. Mech. Anal. 223 2017, 3, 1337–1374], to higher dimensions, and to revisit results in the framework of generalized special functions of bounded variation ( GSBV p ).

1 Introduction

Integral representation results are a fundamental tool in the abstract theory of variational limits by Γ-convergence or in relaxation problems (see [32]). The topic has attracted widespread attention in the mathematical community over the last decades, with applications in various contexts, such as homogenization, dimension reduction, or atomistic-to-continuum approximations. In this paper we contribute to this topic by proving an integral representation result for a general class of energies arising in the modeling of linear elastic solids with surface discontinuities.

Integral representation theorems have been provided with increasing generality, ranging from functionals defined on Sobolev spaces [1, 17, 18, 19, 33, 49] to those defined on spaces of functions of bounded variation [12, 22, 30, 14], in particular on the subspace SBV of special functions of bounded variation [13, 15, 16] and on piecewise constant functions [2]. In recent years, this analysis has been further improved to deal with functionals and variational limits on GSBV p (generalized special functions of bounded variation with p-integrable bulk density), which is the natural energy space for the variational description of many problems with free discontinuities, see among others [6, 7, 8, 9, 21, 37, 41]. A very general method for dealing with all the abovementioned classes of functionals, the so-called global method for relaxation, has been developed by Bouchitté, Fonseca, Leoni, and Mascarenhas in [13, 14]. It essentially consists in comparing asymptotic Dirichlet problems on small balls with different boundary data depending on the local properties of the functions and allows to characterize energy densities in terms of cell formulas.

When coming to the variational description of rupture phenomena in general linearly elastic materials, however, the functional setting to be considered becomes weaker. Indeed, problems need to be formulated in suitable subspaces of functions of bounded deformation ( BD functions) for which the distributional symmetrized gradient is a bounded Radon measure.

In the mathematical description of linear elasticity, the elastic properties are determined by the elastic strain. For a solid in a (bounded) reference configuration Ω d , whose displacement field with respect to the equilibrium is u : Ω d , the elastic strain is given by the symmetrized gradient e ( u ) = 1 2 ( u + ( u ) T ) . In standard models, the corresponding linear elastic energy is a suitable quadratic form of e ( u ) , possibly depending on the material point, see, e.g., [38, Section 2.1]. However, this is often generalized to the case of p-growth for a power p > 1 , see [46, Sections 10–11]. The presence of surface discontinuities is related to several dissipative phenomena, such as cracks, surface tension between different elastic phases, or internal cavities. In the energetic description, this is represented by a term concentrated on the jump set J u . This set is characterized by the property that for x J u , when blowing up around x, the jump set approximates a hyperplane with normal ν u ( x ) 𝕊 d - 1 and the displacement field is close to two suitable values u + ( x ) , u - ( x ) d on the two sides of the material with respect to this hyperplane.

Prototypical examples of functionals described above are energies which are controlled from above and below by suitable multiples of

(1) Ω | e ( u ) | p d x + J u Ω ( 1 + | [ u ] | ) d d - 1 ,

where [ u ] ( x ) = u + ( x ) - u - ( x ) denotes the jump opening, or which are controlled by multiples of Griffith’s energy [45]

(2) Ω | e ( u ) | p d x + d - 1 ( J u Ω ) .

Whereas in case (1) the energy space is be given by SBD p , a subspace of BD , problems with control of type (2) are naturally formulated on generalized special functions of bounded deformation GSBD p , introduced by Dal Maso [31]. (We refer to Section 3.1 for more details.) The only available integral representation result in this context is due to Conti, Focardi, and Iurlano [28] who considered variational functionals controlled locally in terms of (1) in dimension d = 2 . Let us mention that the behavior is quite different if linear growth on the symmetrized gradient is assumed (corresponding to p = 1 ), as suited for the description of plasticity. In that case, representation results in the framework of BD have been obtained, for instance, in [10, 36] and [23] (see also [35, 48], containing essential tools for the proof).

The goal of the present article is twofold: we generalize the results of [28] for energies with control of type (1) to arbitrary space dimensions and, more importantly, we extend the theory to encompass also problems of the form (2), which are most relevant from an applicative viewpoint. Indeed, already in dimension two, the extension of [28] to the case where only a control of type (2) is available is no straightforward task. This is a fundamental difference with respect to the BV -theory where problems for generalized functions of bounded variation can be reconducted to SBV by a perturbation trick (see for instance [21]): one considers a small perturbation of the functional, depending on the jump opening, to represent functionals on SBV p . Then, by letting the perturbation parameter vanish and by truncating functions suitably, the representation can be extended to GSBV p . Unfortunately, the trick of reducing problem (2) to (1) is not expedient in the linearly elastic context and does not allow to deduce an integral representation result in GSBD p from the one in SBD p . This is mainly due to the fact that, given a control only on the symmetrized gradient, it is in principle not possible to use smooth truncations to decrease the energy up to a small error.

Let us also remark that, while in the majority of integral representation results in BV and BD the L 1 -topology was considered, this is not the right choice when only a lower bound of the form (2) is at hand. Indeed, in this case, the available compactness results [25, 31] have been established with respect to the topology of the convergence in measure. This latter is also the topology where recently an integral representation result for the subspace P R ( Ω ) of piecewise rigid functions has been proved in [44].

In our main result (Theorem 2.1), we prove an integral representation for variational functionals

: GSBD p ( Ω ) × ( Ω ) [ 0 , + )

( ( Ω ) denoting the Borel subsets of Ω) that satisfy the standard abstract conditions to be Borel measures in the second argument, lower semicontinuous with respect to convergence in measure, and local in the first argument. Moreover, we require control of type (2), localized to any B ( Ω ) .

Let us comment on the proof strategy. We follow the general approach of the global method for relaxation provided in [13, 14] for variational functionals in BV . The proof strategy recovers the integral bulk and surface densities as blow-up limits of cell minimization formulas. The steps to be performed are the following:

  1. one first shows that, for fixed u GSBD p ( Ω ) , the set function ( u , ) is asymptotically equivalent to its minimum 𝐦 ( u , ) over competitors attaining the same boundary conditions as u on the boundaries of small balls centered in x 0 Ω with vanishing radii. With this we mean that the two quantities have the same Radon–Nikodym derivative with respect to μ := d Ω + d - 1 J u Ω (Lemma 4.1),

  2. one then proves that the Radon–Nikodym derivative d 𝐦 ( u , ) d μ only depends on x 0 , the value u ( x 0 ) , and the (approximate) gradient u ( x 0 ) at a Lebesgue point x 0 , while at a jump point x 0 it is uniquely determined by the one-sided traces u + ( x 0 ) , u - ( x 0 ) and the normal vector ν u ( x 0 ) to J u in x 0 (Lemmas 4.2 and 4.3).

When dealing with all of the abovementioned issues, a key ingredient is given by a Korn-type inequality for special functions of bounded deformation, established recently by Cagnetti, Chambolle, and Scardia [20], which generalizes a two-dimensional result in [28] (see also [39]) to arbitrary dimension. It provides a control of the full gradient in terms of the symmetrized gradient, up to an exceptional set whose perimeter has a surface measure comparable to that of the discontinuity set. In particular, this estimate is used to approximate the function u with functions u ε , which have Sobolev regularity in a ball (around a Lebesgue point), or in half-balls oriented by the jump normal (around a jump point), and which converge to the purely elastic competitor u ( x 0 ) + u ( x 0 ) ( - x 0 ) , or the two-valued function with values u - ( x 0 ) and u + ( x 0 ) , respectively. This is done in Lemmas 5.1 and 6.1, respectively, and is used for proving Lemmas 4.2 and 4.3. Let us mention that this application of the Korn-type inequality is similar to the one in dimension two [28] (with the topology of convergence in measure in place of L 1 ), and constitutes the counterpart of the SBV -Poincaré inequality [34] used in the SBV -case [13]. We also point out that our construction for approximating two-valued functions in Lemma 4.3 slightly differs from the ones in [13, 28] in order to fix a possible flaw contained in these proofs, see Remark 6.2 for details.

In contrast to [28], the Korn inequality is also used in the proof of Lemma 4.1: at this point, one needs to show that functions of the form

v δ := v i δ χ B i δ

approximate u in the topology of the convergence in measure, where B i δ is a fine cover of a given set with disjoint balls of radius smaller than δ and v i δ denote minimizers for 𝐦 ( u , B i δ ) . In [28], the lower bound in (1) allows to control the distributional symmetrized gradient E u which along with a scaling argument and the classical Korn–Poincaré inequality in BD (see [50, Theorem 2.2]) shows that v i δ is close to u on each B i δ . (In [13], the SBV -Poincaré inequality is used.) Our weaker lower bound of the form (2), however, calls for novel arguments and we use the Korn-type inequality to show that v i δ are close to u in L p up to exceptional sets ω i δ whose volumes scale like δ ( ( u , B i δ ) + μ ( B i δ ) ) .

We also point out that, if instead a control of the type (1) is assumed, the arguments leading to Theorem 2.1 can be successfully adapted to extend the result for functionals on SBD p (see [28]) to arbitrary space dimensions, see Theorem 7.1. This is done by exploiting the stronger blow-up properties of SBD functions. We note that, in principle, this result could be also obtained by adapting the arguments in [28] to higher dimension by employing the Korn inequality [20]. We however preferred to give a self-contained proof of Theorem 7.1, which requires only slight modifications of the arguments used for Theorem 2.1 and nicely illustrates the differences between SBD p and its generalized space.

For a related purpose, in Section 8 we discuss how our arguments can also provide a direct proof for integral representation results on GSBV p , if a local control on the full deformation gradient of the form

Ω | u | p d x + d - 1 ( J u Ω )

is given, see Theorem 8.1. In particular, no perturbation or truncation arguments are needed in the proof. Therefore, we believe that this provides a new perspective and a slightly simpler approach to integral representation results in GSBV p without necessity of the perturbation trick discussed before, relying on the SBV result. Let us, however, mention that in [21] a more general growth condition from above is considered: dealing with such a condition would instead require a truncation method in the proof.

We close the introduction by mentioning that in a subsequent work [42] we use the present result to obtain integral representation of Γ-limits for sequences of energies in linear elasticity with surface discontinuities. There, we additionally characterize the bulk and surface densities as blow-up limits of cell minimization formulas where the minimization is not performed on GSBD p but more specifically on Sobolev functions (bulk density) and piecewise rigid functions [44] (surface density). The latter characterization particularly allows to identify integrands of relaxed functionals and to treat homogenization problems.

The paper is organized as follows. In Section 2 we present our main integral representation result in GSBD p . Section 3 is devoted to some preliminaries about the function space. In particular, we present the Korn-type inequality established in [20] and prove a fundamental estimate. Section 4 contains the general strategy and the proof of Lemma 4.1. The identifications of the bulk and surface density (Lemmas 4.2 and 4.3) are postponed to Sections 5 and 6, respectively. In Section 7 we describe the modifications necessary to obtain the SBD p -case. Finally, in Section 8 we explain how our method can be used to establish an integral representation result in GSBV p .

2 The integral representation result

In this section we present our main result. We start with some basic notation. Let Ω d be open, bounded with Lipschitz boundary. Let 𝒜 ( Ω ) be the family of open subsets of Ω, and denote by ( Ω ) the family of Borel sets contained in Ω. For every x d and ε > 0 we indicate by B ε ( x ) d the open ball with center x and radius ε. For x, y d , we use the notation x y for the scalar product and | x | for the Euclidean norm. Moreover, we let 𝕊 d - 1 := { x d : | x | = 1 } and we denote by 𝕄 d × d the set of d × d matrices. The m-dimensional Lebesgue measure of the unit ball in m is indicated by γ m for every m . We denote by d and k the d-dimensional Lebesgue measure and the k-dimensional Hausdorff measure, respectively.

For definition and properties of the space GSBD p ( Ω ) , 1 < p < , we refer the reader to [31]. Some relevant properties are collected in Section 3 below. In particular, the approximate gradient is denoted by u (it is well defined, see Lemma 3.5) and the (approximate) jump set is denoted by J u with corresponding normal ν u and one-sided limits u + and u - . We also define e ( u ) = 1 2 ( u + ( u ) T ) .

We consider functionals : GSBD p ( Ω ) × ( Ω ) [ 0 , + ) with the following general assumptions:

  1. ( u , ) is a Borel measure for any u GSBD p ( Ω ) ,

  2. ( , A ) is lower semicontinuous with respect to convergence in measure on Ω for any A 𝒜 ( Ω ) ,

  3. ( , A ) is local for any A 𝒜 ( Ω ) , in the sense that if u , v GSBD p ( Ω ) satisfy u = v a.e. in A, then ( u , A ) = ( v , A ) ,

  4. there exist 0 < α < β such that for any u GSBD p ( Ω ) and B ( Ω ) we have

    α ( B | e ( u ) | p d x + d - 1 ( J u B ) ) ( u , B ) β ( B ( 1 + | e ( u ) | p ) d x + d - 1 ( J u B ) ) .

We now formulate the main result of this article addressing integral representation of functionals satisfying (H1)(H4). To this end, we introduce some further notation: for every u GSBD p ( Ω ) and A 𝒜 ( Ω ) we define

(2.1) 𝐦 ( u , A ) = inf v GSBD p ( Ω ) { ( v , A ) : v = u  in a neighborhood of  A } .

For x 0 Ω , u 0 d , and ξ 𝕄 d × d we introduce the functions x 0 , u 0 , ξ : d d by

(2.2) x 0 , u 0 , ξ ( x ) = u 0 + ξ ( x - x 0 ) .

Moreover, for x 0 Ω , a , b d , and ν 𝕊 d - 1 we introduce u x 0 , a , b , ν : d d by

(2.3) u x 0 , a , b , ν ( x ) = { a if  ( x - x 0 ) ν > 0 , b if  ( x - x 0 ) ν < 0 .

In this paper, we will prove the following result.

Theorem 2.1 (Integral representation in GSBD p ).

Let Ω R d be open, bounded with Lipschitz boundary and suppose that F : GSBD p ( Ω ) × B ( Ω ) [ 0 , + ) satisfies (H1)(H4). Then

( u , B ) = B f ( x , u ( x ) , u ( x ) ) d x + J u B g ( x , u + ( x ) , u - ( x ) , ν u ( x ) ) d d - 1 ( x )

for all u GSBD p ( Ω ) and B B ( Ω ) , where f is given by

(2.4) f ( x 0 , u 0 , ξ ) = lim sup ε 0 𝐦 ( x 0 , u 0 , ξ , B ε ( x 0 ) ) γ d ε d

for all x 0 Ω , u 0 R d , ξ M d × d , and g is given by

(2.5) g ( x 0 , a , b , ν ) = lim sup ε 0 𝐦 ( u x 0 , a , b , ν , B ε ( x 0 ) ) γ d - 1 ε d - 1

for all x 0 Ω , a , b R d , and ν S d - 1 .

Remark 2.2.

We proceed with some remarks on the result.

  1. In general, if f is not convex in ξ, in spite of the growth conditions (H4), the functional may fully depend on u and not just on the symmetric part e ( u ) . We refer to [28, Remark 4.14] for an example in this direction.

  2. As is lower semicontinuous on W 1 , p with respect to weak convergence, the integrand f is quasiconvex [47]. Since is lower semicontinuous on piecewise rigid functions, the integrand g is BD -elliptic [43] (at least if one can ensure, for instance, that g has a continuous dependence in x). A fortiori, g is BV -elliptic [3].

  3. If the functional additionally satisfies ( u + a , A ) = ( u , A ) for all affine functions a : d d with e ( a ) = 0 , then there are two functions f : Ω × 𝕄 d × d [ 0 , + ) and g : Ω × d × 𝕊 d - 1 [ 0 , + ) such that

    ( u , B ) = B f ( x , e ( u ) ( x ) ) d x + J u B g ( x , [ u ] ( x ) , ν u ( x ) ) d d - 1 ( x ) ,

    where [ u ] ( x ) := u + ( x ) - u - ( x ) .

  4. A variant of the proof shows that, in the minimization problems (2.4)–(2.5), one may replace balls B ε ( x 0 ) by cubes Q ε ν ( x 0 ) with sidelength ε, centered at x 0 , and two faces orthogonal to ν = ν u ( x 0 ) .

  5. An analogous result holds on the space GSBV p ( Ω ; m ) for m . We refer to Section 8 for details.

We will additionally discuss the minor modifications needed in order to deal with functionals

: SBD p ( Ω ) × ( Ω ) [ 0 , + )

satisfying (H1)(H3) and

  1. there exist 0 < α < β such that for any u SBD p ( Ω ) and B ( Ω ) we have

    α ( B | e ( u ) | p d x + J u B ( 1 + | [ u ] | ) d d - 1 ) ( u , B ) β ( B ( 1 + | e ( u ) | p ) d x + J u B ( 1 + | [ u ] | ) d d - 1 ) .

In this case, SBD p ( Ω ) (see Section 3.1) is the natural energy space for . Furthermore, sequences of competitors with bounded energy, which are converging in measure, are additionally L 1 -convergent if we assume (H4’), due to the classical Korn–Poincaré inequality in BD (see [50, Theorem II.2.2]). Hence, in this latter case, (H2) is equivalent to requiring lower semicontinuity with respect to the L 1 -convergence. The statement of the result in this setting, as well as of the changes needed in the proofs, will be given in Section 7.

3 Preliminaries

We start this preliminary section by introducing some further notation. For E d , ε > 0 , and x 0 d we set

(3.1) E ε , x 0 := x 0 + ε ( E - x 0 ) .

The diameter of E is indicated by diam ( E ) . Given two sets E 1 , E 2 d , we denote their symmetric difference by E 1 E 2 . We write χ E for the characteristic function of any E d , which is 1 on E and 0 otherwise. If E is a set of finite perimeter, we denote its essential boundary by * E , see [5, Definition 3.60]. We denote the set of symmetric and skew-symmetric matrices by 𝕄 sym d × d and 𝕄 skew d × d , respectively.

3.1 BD and GBD functions

Let U d be open. A function v L 1 ( U ; d ) belongs to the space of functions of bounded deformation, denoted by BD ( U ) , if the distribution E v := 1 2 ( D v + ( D v ) T ) is a bounded 𝕄 sym d × d -valued Radon measure on U, where D v = ( D 1 v , , D d v ) is the distributional differential. It is well known (see [4, 50]) that for v BD ( U ) the jump set J v is countably d - 1 -rectifiable (in the sense of [5, Definition 2.57]), and that

E v = E a v + E c v + E j v ,

where E a v is absolutely continuous with respect to d , E c v is singular with respect to d and such that | E c v | ( B ) = 0 if d - 1 ( B ) < , while E j v is concentrated on J v . The density of E a v with respect to d is denoted by e ( v ) .

The space SBD ( U ) is the subspace of all functions v BD ( U ) such that E c v = 0 . For p ( 1 , ) , we define SBD p ( U ) := { v SBD ( U ) : e ( v ) L p ( U ; 𝕄 sym d × d ) , d - 1 ( J v ) < } . For a complete treatment of BD and SBD functions, we refer to to [4, 11, 50].

The spaces GBD ( U ) of generalized functions of bounded deformation and GSBD ( U ) GBD ( U ) of generalized special functions of bounded deformation have been introduced in [31] (cf. [31, Definitions 4.1 and 4.2]), and are defined as follows.

Definition 3.1.

Let U d be a bounded open set, and let v : U d be measurable. We introduce the notation

Π ξ := { y d : y ξ = 0 } , B y ξ := { t : y + t ξ B }  for any  y d  and  B d

for fixed ξ 𝕊 d - 1 , and for every t B y ξ we let

v y ξ ( t ) := v ( y + t ξ ) , v ^ y ξ ( t ) := v y ξ ( t ) ξ .

Then v GBD ( U ) if there exists λ v b + ( U ) such that v ^ y ξ BV loc ( U y ξ ) for d - 1 -a.e. y Π ξ , and for every Borel set B U ,

Π ξ ( | D v ^ y ξ | ( B y ξ J v ^ y ξ 1 ) + 0 ( B y ξ J v ^ y ξ 1 ) ) d d - 1 ( y ) λ v ( B ) , J v ^ y ξ 1 := { t J v ^ y ξ : | [ v ^ y ξ ] | ( t ) 1 } .

Moreover, the function v belongs to GSBD ( U ) if v GBD ( U ) and v ^ y ξ SBV loc ( U y ξ ) for every ξ 𝕊 d - 1 and for d - 1 -a.e. y Π ξ .

We recall that every v GBD ( U ) has an approximate symmetric gradient e ( v ) L 1 ( U ; 𝕄 sym d × d ) and an approximate jump set J v which is still countably d - 1 -rectifiable (cf. [31, Theorem 9.1, Theorem 6.2]).

The notation for e ( v ) and J v , which is the same as that one in the SBD case, is consistent: in fact, if v lies in SBD ( U ) , the objects coincide (up to negligible sets of points with respect to d and d - 1 , respectively). For x J v there exist v + ( x ) , v - ( x ) d and ν v ( x ) 𝕊 d - 1 such that

(3.2) lim ε 0 ε - d d ( { y B ε ( x ) : ± ( y - x ) ν v ( x ) > 0 } { | v - v ± ( x ) | > ϱ } ) = 0

for every ϱ > 0 , and the function [ v ] := v + - v - : J v d is measurable. For 1 < p < , the space GSBD p ( U ) is given by

GSBD p ( U ) := { v GSBD ( U ) : e ( v ) L p ( U ; 𝕄 sym d × d ) , d - 1 ( J v ) < } .

Any function v GSBD ( U ) with [ v ] integrable belongs to SBD ( U ) , as follows from [27, Theorem 2.9] for 𝔸 v = E v (see [27, Remark 2.5]). This corresponds to the following proposition.

Proposition 3.2.

If v GSBD p ( U ) is such that [ v ] L 1 ( J v ; R d ) , then v SBD p ( U ) .

If U has Lipschitz boundary, for each v GBD ( U ) the traces on U are well defined (see [31, Theorem 5.5]), in the sense that for d - 1 -a.e. x U there exists tr ( v ) ( x ) d such that

(3.3) lim ε 0 ε - d d ( U B ε ( x ) { | v - tr ( v ) ( x ) | > ϱ } ) = 0    for all  ϱ > 0 .

3.2 Korn’s inequality and fundamental estimate

In this subsection we discuss two important tools which will be instrumental for the proof of Theorem 2.1. We start by the following Korn and Korn–Poincaré inequalities in GSBD for functions with small jump sets, see [20, Theorem 1.1, Theorem 1.2]. In the following, we say that a : d d is an infinitesimal rigid motion if a is affine with e ( a ) = 1 2 ( a + ( a ) T ) = 0 .

Theorem 3.3 (Korn inequality for functions with small jump set).

Let Ω R d be a bounded Lipschitz domain and let 1 < p < + . Then there exists a constant c = c ( Ω , p ) > 0 such that for all u GSBD p ( Ω ) there is a set of finite perimeter ω Ω with

(3.4) d - 1 ( * ω ) c d - 1 ( J u ) , d ( ω ) c ( d - 1 ( J u ) ) d d - 1

and an infinitesimal rigid motion a such that

(3.5) u - a L p ( Ω ω ) + u - a L p ( Ω ω ) c e ( u ) L p ( Ω ) .

Moreover, there exists v W 1 , p ( Ω ; R d ) such that v = u on Ω ω and

e ( v ) L p ( Ω ) c e ( u ) L p ( Ω ) .

Note that the result is indeed only relevant for functions with sufficiently small jump set, as otherwise one can choose ω = Ω , and (3.5) trivially holds. Note that, in [20], d ( ω ) c ( d - 1 ( J u ) ) d d - 1 has not been shown, but it readily follows from d - 1 ( * ω ) c d - 1 ( J u ) by the isoperimetric inequality.

Remark 3.4 (Almost Sobolev regularity, constants, and scaling invariance).

(i) More precisely, in [20] it is proved that there exists v W 1 , p ( Ω ; d ) such that v = u on Ω ω and e ( v ) L p ( Ω ) c e ( u ) L p ( Ω ) , whence by Korn’s and Poincaré’s inequality in W 1 , p ( Ω ; d ) we get

v - a L p ( Ω ) + v - a L p ( Ω ) c e ( u ) L p ( Ω )

for an infinitesimal rigid motion a. This directly implies (3.5), see [20, Theorem 4.1, Theorem 4.4].

(ii) Given a collection of bounded Lipschitz domains ( Ω k ) k which are related through bi-Lipschitzian homeomorphisms with Lipschitz constants of both the homeomorphism itself and its inverse bounded uniformly in k, in Theorem 3.3 we can choose a constant c uniformly for all Ω k , see [20, Remark 4.2].

(iii) Recall definition (3.1). Consider a bounded Lipschitz domain Ω, ε > 0 , and x 0 d . Then for each u GSBD p ( Ω ε , x 0 ) we find ω Ω ε , x 0 and a rigid motion a such that

d - 1 ( * ω ) C d - 1 ( J u ) , d ( ω ) C ( d - 1 ( J u ) ) d d - 1

and

ε - 1 u - a L p ( Ω ε , x 0 ω ) + u - a L p ( Ω ε , x 0 ω ) C e ( u ) L p ( Ω ε , x 0 ) ,

where C = C ( Ω , p ) > 0 is independent of ε. This follows by a standard rescaling argument.

From Theorem 3.3, one can also deduce that for u GSBD p ( Ω ) the approximate gradient u exists d -a.e. in Ω, see [20, Corollary 5.2].

Lemma 3.5 (Approximate gradient).

Let Ω R d be open, bounded with Lipschitz boundary, let 1 < p < + , and u GSBD p ( Ω ) . Then for L d -a.e. x 0 Ω there exists a matrix in M d × d , denoted by u ( x 0 ) , such that

lim ε 0 ε - d d ( { x B ε ( x 0 ) : | u ( x ) - u ( x 0 ) - u ( x 0 ) ( x - x 0 ) | | x - x 0 | > ϱ } ) = 0 for all  ϱ > 0 .

We point out that the result in Lemma 3.5 has already been obtained in [40] for p = 2 , as a consequence of the embedding GSBD 2 ( Ω ) ( GBV ( Ω ) ) d , see [40, Theorem 2.9].

To control the affine mappings appearing in Theorem 3.3, we will make use of the following elementary lemma on affine mappings, see, e.g., [44, Lemma 3.4] or [29, Lemmas 4.3] for similar statements. (It is obtained by the equivalence of norms in finite dimensions and by standard rescaling arguments.)

Lemma 3.6.

Let 1 p < + , let x 0 R d , and let R , θ > 0 . Let a : R d R d be affine, defined by a ( x ) = A x + b for x R d , and let E B R ( x 0 ) R d with L d ( E ) θ L d ( B R ( x 0 ) ) . Then there exists a constant c 0 > 0 only depending on p and θ such that

a L p ( B R ( x 0 ) ) γ d 1 p R d p a L ( B R ( x 0 ) ) c 0 a L p ( E ) , | A | c 0 R - 1 - d p a L p ( E ) .

We now proceed with another consequence of Theorem 3.3.

Corollary 3.7.

Let Ω R d be a bounded Lipschitz domain and let 1 < p < + . Then there exists a constant C = C ( Ω , p ) > 0 such that for all u GSBD p ( Ω ) with trace tr ( u ) = 0 on Ω (see (3.3)) there is a set of finite perimeter ω R d with

(3.6) d - 1 ( * ω ) C d - 1 ( J u ) , d ( ω ) C ( d - 1 ( J u ) ) d d - 1

such that

(3.7) u L p ( Ω ω ) + u L p ( Ω ω ) C e ( u ) L p ( Ω ) .

Proof.

We start by choosing a bounded Lipschitz domain Ω d with Ω Ω . Each u GSBD p ( Ω ) with tr ( u ) = 0 on Ω can be extended to a function u ~ GSBD p ( Ω ) by u ~ = 0 on Ω Ω such that J u ~ = J u . We first note that it is not restrictive to assume that

(3.8) d - 1 ( J u ) ( d ( Ω Ω ) 2 c ) d - 1 d ,

where c = c ( Ω , p ) > 0 is the constant of Theorem 3.3. In fact, otherwise we could take ω = Ω and the statement would be trivially satisfied since (3.7) is clearly trivial and for (3.6) we use that

d - 1 ( * Ω ) C ( d ( Ω Ω ) 2 c ) d - 1 d , d ( Ω ) C d ( Ω Ω ) 2 c

for a sufficiently large constant C > 0 depending only on Ω and Ω .

Now, consider a function u satisfying (3.8). We apply Theorem 3.3 on u ~ GSBD p ( Ω ) and obtain a set ω Ω d satisfying (3.4) as well as an infinitesimal rigid motion a such that

(3.9) u ~ - a L p ( Ω ω ) + u ~ - a L p ( Ω ω ) c e ( u ~ ) L p ( Ω ) = c e ( u ) L p ( Ω ) .

In particular, u ~ = 0 on Ω Ω implies

(3.10) a L p ( Ω ( Ω ω ) ) c e ( u ) L p ( Ω ) .

By (3.4) and (3.8) we get d ( ω ) 1 2 d ( Ω Ω ) . In view of (3.10), we apply Lemma 3.6 on E = Ω ( Ω ω ) with R = diam ( Ω ) and θ = 1 2 d ( Ω Ω ) / γ d R d to get

a L p ( Ω ) c a L p ( Ω ( Ω ω ) ) c e ( u ) L p ( Ω ) ,

and, in a similar fashion, | a | c e ( u ) L p ( Ω ) , where c > 0 depends on Ω, Ω , and p. Then (3.7) follows from (3.9), the triangle inequality, and the fact that u = u ~ on Ω. ∎

We conclude this subsection with another important tool in the proof of the integral representation, namely a fundamental estimate in GSBD p .

Lemma 3.8 (Fundamental estimate in GSBD p ).

Let Ω R d be open, bounded with Lipschitz boundary, and let 1 < p < + . Let η > 0 and let A , A , A ′′ A ( Ω ) with A A . For every functional F satisfying (H1), (H3), and (H4) and for every u GSBD p ( A ) , v GSBD p ( A ′′ ) there exists a function φ C ( R d ; [ 0 , 1 ] ) such that w := φ u + ( 1 - φ ) v GSBD p ( A A ′′ ) satisfies

(3.11) { (i) ( w , A A ′′ ) ( 1 + η ) ( ( u , A ) + ( v , A ′′ ) ) + M u - v L p ( ( A A ) A ′′ ) p + η d ( A A ′′ ) , (ii) w = u on  A    𝑎𝑛𝑑    w = v on  A ′′ A ,

where M = M ( A , A , A ′′ , p , η ) > 0 depends only on A , A , A ′′ , p , η , but is independent of u and v. Moreover, if for ε > 0 and x 0 R d we have A ε , x 0 , A ε , x 0 , A ε , x 0 ′′ Ω , then

(3.12) M ( A ε , x 0 , A ε , x 0 , A ε , x 0 ′′ , p , η ) = ε - p M ( A , A , A ′′ , p , η ) ,

where we used the notation introduced in (3.1).

The same statement holds if F satisfies (H4’), u SBD p ( A ) , and v SBD p ( A ′′ ) .

In the statement above, we intend that u - v L p ( ( A A ) A ′′ ) p = + if u - v L p ( ( A A ) A ′′ ) .

Proof.

The proof follows the lines of [16, Proposition 3.1]. Choose k such that

(3.13) k max { 3 p - 1 β η α , β η } .

Let A 1 , , A k + 1 be open subsets of d with

A A 1 A k + 1 A .

For i = 1 , , k let φ i C 0 ( A i + 1 ; [ 0 , 1 ] ) with φ i = 1 in a neighborhood V i of A i ¯ .

Consider u GSBD p ( A ) and v GSBD p ( A ′′ ) . We can clearly assume that u - v L p ( ( A A ) A ′′ ) as otherwise the result is trivial. We define the function w i = φ i u + ( 1 - φ i ) v GSBD p ( A A ′′ ) , where u and v are extended arbitrarily outside A and A ′′ , respectively. Letting T i = A ′′ ( A i + 1 A i ¯ ) we get by (H1) and (H3)

( w i , A A ′′ ) ( u , ( A A ′′ ) V i ) + ( v , A ′′ supp ( φ i ) ) + ( w i , T i )
(3.14) ( u , A ) + ( v , A ′′ ) + ( w i , T i ) .

For the last term, we compute using (H4) ( denotes the symmetrized vector product)

( w i , T i ) β T i ( 1 + | e ( w i ) | p ) d x + β d - 1 ( J w i T i )
β T i ( 1 + | φ i e ( u ) + ( 1 - φ i ) e ( v ) + φ i ( u - v ) | p ) + β d - 1 ( ( J u J v ) T i )
β d ( T i ) + 3 p - 1 β T i ( | e ( u ) | p + | e ( v ) | p + | φ i | p | u - v | p ) + β d - 1 ( J u T i ) + β d - 1 ( J v T i )
3 p - 1 β α - 1 ( ( u , T i ) + ( v , T i ) ) + 3 p - 1 β ϕ i p u - v L p ( T i ) p + β d ( T i ) .

Notice that we can obtain the same estimate also if satisfies (H4’), u SBD p ( A ) , and v SBD p ( A ′′ ) . (We refer to [16, proof of Proposition 3.1] for details.) Consequently, recalling (3.13) and using (H1), we find i 0 { 1 , , k } such that

( w i 0 , T i 0 ) 1 k i = 1 k ( w i , T i ) η ( ( u , A ) + ( v , A ′′ ) ) + M u - v L p ( ( A A ) A ′′ ) p + η d ( ( A A ) A ′′ ) ,

where M := 3 p - 1 β k - 1 max i = 1 , , k φ i p . This along with (3.14) concludes the proof of (3.11) by setting w = w i 0 . To see the scaling property (3.12), it suffices to use the cut-off functions φ i ε C 0 ( ( A i + 1 ) ε , x 0 ; [ 0 , 1 ] ) , i = 1 , , k , defined by φ i ε ( x ) = φ i ( x 0 + 1 ε ( x - x 0 ) ) for x ( A i + 1 ) ε , x 0 . This concludes the proof. ∎

4 The global method

This section is devoted to the proof of Theorem 2.1 which is based on three ingredients. First, we show that is equivalent to 𝐦 (see (2.1)) in the sense that the two quantities have the same Radon–Nikodym derivative with respect to μ := d Ω + d - 1 J u Ω .

Lemma 4.1.

Suppose that F satisfies (H1)(H4). Let u GSBD p ( Ω ) and μ = L d Ω + H d - 1 J u Ω . Then for μ-a.e. x 0 Ω we have

lim ε 0 ( u , B ε ( x 0 ) ) μ ( B ε ( x 0 ) ) = lim ε 0 𝐦 ( u , B ε ( x 0 ) ) μ ( B ε ( x 0 ) ) .

We prove this lemma in the final part of this section. The second ingredient is that, asymptotically as ε 0 , the minimization problems 𝐦 ( u , B ε ( x 0 ) ) and 𝐦 ( u ¯ x 0 bulk , B ε ( x 0 ) ) coincide for d -a.e. x 0 Ω , where we write u ¯ x 0 bulk := x 0 , u ( x 0 ) , u ( x 0 ) for brevity, see (2.2).

Lemma 4.2.

Suppose that F satisfies (H1) and (H3)(H4) and let u GSBD p ( Ω ) . Then for L d -a.e. x 0 Ω we have

(4.1) lim ε 0 𝐦 ( u , B ε ( x 0 ) ) γ d ε d = lim sup ε 0 𝐦 ( u ¯ x 0 bulk , B ε ( x 0 ) ) γ d ε d .

We defer the proof of Lemma 4.2 to Section 5. The third ingredient is that, asymptotically as ε 0 , the minimization problems 𝐦 ( u , B ε ( x 0 ) ) and 𝐦 ( u ¯ x 0 surf , B ε ( x 0 ) ) coincide for d - 1 -a.e. x 0 J u , where we write u ¯ x 0 surf := u x 0 , u + ( x 0 ) , u - ( x 0 ) , ν u ( x 0 ) for brevity, see (2.3).

Lemma 4.3.

Suppose that F satisfies (H1) and (H3)(H4) and let u GSBD p ( Ω ) . Then for H d - 1 -a.e. x 0 J u we have

(4.2) lim ε 0 𝐦 ( u , B ε ( x 0 ) ) γ d - 1 ε d - 1 = lim sup ε 0 𝐦 ( u ¯ x 0 surf , B ε ( x 0 ) ) γ d - 1 ε d - 1 .

We defer the proof of Lemma 4.3 to Section 6, and now proceed to prove Theorem 2.1.

Proof of Theorem 2.1.

In view of (H4) on and of the Besicovitch derivation theorem (cf. [5, Theorem 2.22]), we need to show that for d -a.e. x 0 Ω one has

(4.3) d ( u , ) d d ( x 0 ) = f ( x 0 , u ( x 0 ) , u ( x 0 ) ) ,

where f was defined in (2.4), and that for d - 1 -a.e. x 0 J u one has

(4.4) d ( u , ) d d - 1 J u ( x 0 ) = g ( x 0 , u + ( x 0 ) , u - ( x 0 ) , ν u ( x 0 ) ) ,

where g was defined in (2.5).

By Lemma 4.1 and the fact that lim ε 0 ( γ d ε d ) - 1 μ ( B ε ( x 0 ) ) = 1 for d -a.e. x 0 Ω we deduce

d ( u , ) d d ( x 0 ) = lim ε 0 ( u , B ε ( x 0 ) ) μ ( B ε ( x 0 ) ) = lim ε 0 𝐦 ( u , B ε ( x 0 ) ) μ ( B ε ( x 0 ) ) = lim ε 0 𝐦 ( u , B ε ( x 0 ) ) γ d ε d <

for d -a.e. x 0 Ω . Then (4.3) follows from (2.4) and Lemma 4.2. By Lemma 4.1 and the fact that

lim ε 0 ( γ d - 1 ε d - 1 ) - 1 μ ( B ε ( x 0 ) ) = 1

for d - 1 -a.e. x 0 J u we deduce

d ( u , ) d d - 1 J u ( x 0 ) = lim ε 0 ( u , B ε ( x 0 ) ) μ ( B ε ( x 0 ) ) = lim ε 0 𝐦 ( u , B ε ( x 0 ) ) μ ( B ε ( x 0 ) ) = lim ε 0 𝐦 ( u , B ε ( x 0 ) ) γ d - 1 ε d - 1 <

for d - 1 -a.e. x 0 J u . Now, (4.4) follows from (2.5) and Lemma 4.3. ∎

In the remaining part of the section we prove Lemma 4.1. We basically follow the lines of [13, 14, 28], with the difference that the required compactness results are more delicate due to the weaker growth condition from below (see H4) compared to [13, 14, 28]. We start with some notation. For δ > 0 and A 𝒜 ( Ω ) , we define

𝐦 δ ( u , A ) = inf { i = 1 𝐦 ( u , B i ) : B i A  pairwise disjoint balls , diam ( B i ) δ , μ ( A i = 1 B i ) = 0 } ,

where, as before, μ = d Ω + d - 1 J u Ω . As 𝐦 δ ( u , A ) is decreasing in δ, we can also introduce

(4.5) 𝐦 * ( u , A ) = lim δ 0 𝐦 δ ( u , A ) .

In the following lemma, we prove that and 𝐦 * coincide under our assumptions.

Lemma 4.4.

Suppose that F satisfies (H1)(H4) and let u GSBD p ( Ω ) . Then for all A A ( Ω ) there holds F ( u , A ) = m F * ( u , A ) .

Proof.

We follow the lines of the proof of [28, Lemma 4.1] focusing on the necessary adaptions due to the weaker growth condition from below (see H4) compared to [28]. For each ball B A , 𝐦 ( u , B ) ( u , B ) by definition. By (H1) we get 𝐦 δ ( u , A ) ( u , A ) for all δ > 0 . This shows 𝐦 * ( u , A ) ( u , A ) , cf. (4.5).

We now address the reverse inequality. We fix A 𝒜 ( Ω ) and δ > 0 . Let ( B i δ ) i be balls as in the definition of 𝐦 δ ( u , A ) such that

(4.6) i = 1 𝐦 ( u , B i δ ) 𝐦 δ ( u , A ) + δ .

By the definition of 𝐦 , we find v i δ GSBD p ( B i δ ) such that v i δ = u in a neighborhood of B i δ and

(4.7) ( v i δ , B i δ ) 𝐦 ( u , B i δ ) + δ d ( B i δ ) .

We define

(4.8) v δ , n := i = 1 n v i δ χ B i δ + u χ N 0 δ , n for  n , v δ := i = 1 v i δ χ B i δ + u χ N 0 δ ,

where N 0 δ , n := Ω i = 1 n B i δ and N 0 δ := Ω i = 1 B i δ . By construction, we have that each v δ , n lies in GSBD p ( Ω ) and that sup n ( e ( v δ , n ) L p ( Ω ) + d - 1 ( J v δ , n ) ) < + by (4.6)–(4.7) and (H4). Moreover, v δ , n v δ pointwise a.e. in Ω. Then [25, Theorem 1.1] yields v δ GSBD p ( Ω ) .

( v δ , A ) = i = 1 ( v i δ , B i δ ) + ( u , N 0 δ A )
i = 1 ( 𝐦 ( u , B i δ ) + δ d ( B i δ ) )
(4.9) 𝐦 δ ( u , A ) + δ ( 1 + d ( A ) ) ,

where we also used the fact that μ ( N 0 δ A ) = ( u , N 0 δ A ) = 0 by the definition of ( B i δ ) i and (H4). For later purpose, we also note by (H4) that this implies

(4.10) e ( v δ ) L p ( A ) p + d - 1 ( J v δ A ) α - 1 ( 𝐦 δ ( u , A ) + δ ( 1 + d ( A ) ) ) .

We now claim that v δ u in measure on A. To this end, we apply Remark 3.4 (iii) and Corollary 3.7 on each B i δ for the function u - v i δ and we get sets of finite perimeter ω i δ B i δ such that

(4.11) { (i) ( d ( ω i δ ) ) d - 1 d C d - 1 ( ( J u v δ ) B i δ ) , (ii) u - v i δ L p ( B i δ ω i δ ) p C δ p ( e ( u ) L p ( B i δ ) p + e ( v δ ) L p ( B i δ ) p )

for a constant C > 0 only depending on p. Here, we used that diam ( B i δ ) δ and ( u - v δ ) B i δ GSBD p ( B i δ ) with trace zero on B i δ . We define ψ : [ 0 , + ) [ 0 , + ) by ψ ( t ) = min { t p , 1 } and observe that v δ u in measure on A is equivalent to A ψ ( | u - v δ | ) d x 0 as δ 0 . In view of (4.8), we compute

(4.12) A ψ ( | u - v δ | ) d x = i = 1 B i δ ψ ( | u - v i δ | ) d x i = 1 ( u - v i δ L p ( B i δ ω i δ ) p + d ( ω i δ ) ) .

By (4.11) (ii) and the fact that the balls ( B i δ ) i are pairwise disjoint we get

(4.13) i = 1 u - v i δ L p ( B i δ ω i δ ) p C δ p ( e ( u ) L p ( A ) p + e ( v δ ) L p ( A ) p ) .

As ω i δ B δ i and diam ( B i δ ) δ , we further get by (4.11) (i)

(4.14) i = 1 d ( ω i δ ) γ d 1 d δ i = 1 ( d ( ω i δ ) ) d - 1 d γ d 1 / d C δ d - 1 ( ( J u J v δ ) A ) .

Now, combining (4.12)–(4.14) and using (4.10), we find A ψ ( | u - v δ | ) d x 0 as δ 0 . With this, using (H2), (4.5), and (4.9) we get the required inequality 𝐦 * ( u , A ) ( u , A ) in the limit as δ 0 . This concludes the proof. ∎

Proof of Lemma 4.1.

The statement follows by repeating exactly the arguments in [28, proofs of Lemma 4.2 and Lemma 4.3]. We report a sketch for the reader’s convenience.

The definition of 𝐦 gives readily that for every x 0 Ω ,

lim sup ε 0 𝐦 ( u , B ε ( x 0 ) ) μ ( B ε ( x 0 ) ) lim sup ε 0 ( u , B ε ( x 0 ) ) μ ( B ε ( x 0 ) ) .

The converse inequality follows by proving that, for every t > 0 , the set

E t := { x Ω : lim inf ε 0 ( u , B ε ( x ) ) - 𝐦 ( u , B ε ( x ) ) μ ( B ε ( x ) ) > t }

satisfies μ ( E t ) = 0 . To this end, we fix t > 0 . We introduce the family of balls (depending on t)

X δ := { B ε ( x ) : ε < δ , B ε ( x ) ¯ Ω , μ ( B ε ( x ) ) = 0 , ( u , B ε ( x ) ) > 𝐦 ( u , B ε ( x ) ) + t μ ( B ε ( x ) ) } ,

and we show that

U * := δ > 0 { x Ω : B ε ( x ) X δ  for some  ε > 0 }

satisfies

(4.15) E t U * , μ ( U * ) = 0 .

Then μ ( E t ) = 0 indeed holds true and the proof is concluded.

We now confirm (4.15). The inclusion E t U * follows from the definition of E t which permits to find, for any x E t , ε < δ such that ( u , B ε ( x ) ) > 𝐦 ( u , B ε ( x ) ) + t μ ( B ε ( x ) ) . Note that, due to the left continuity of ε 𝐦 ( u , B ε ( x ) ) (cf. [28, Lemma 4.2]) one can also ensure the additional property μ ( B ε ( x ) ) = 0 by slightly varying ε. In order to prove that μ ( U * ) = 0 , one fixes a compact set K U * , two positive numbers δ < η , and defines

U η := { B ε ( x ) : B ε ( x ) X η } , Y δ := { B ε ( x ) : ε < δ , B ε ( x ) ¯ U η K , μ ( B ε ( x ) ) = 0 } .

By recalling also the definition of U * , we see that X δ and Y δ are fine covers of K and U η K , respectively. Thus there exists countable many pairwise disjoint B i X δ , B ^ j Y δ , and a set N with μ ( N ) = 0 such that U η = i B i j B ^ j N . In view of assumptions (H1), the definitions of X δ , Y δ (in particular the balls have radii smaller than δ) give that

( u , U η ) i 𝐦 ( u , B i ) + j 𝐦 ( u , B ^ j ) + t μ ( i B i ) 𝐦 δ ( u , U η ) + t μ ( K ) .

Passing to the limit in δ, (4.5) and Lemma 4.4 imply

( u , U η ) 𝐦 * ( u , U η ) + t μ ( K ) = ( u , U η ) + t μ ( K ) ,

so that μ ( K ) = 0 . Then μ ( U * ) = 0 by the regularity of μ. ∎

To conclude the proof of Theorem 2.1, it remains to prove Lemmas 4.2 and 4.3. This is the subject of the following two sections.

5 The bulk density

This section is devoted to the proof of Lemma 4.2. We start by analyzing the blow-up at points with approximate gradient. The latter exists for d -a.e. point in Ω by Lemma 3.5.

Lemma 5.1 (Blow-up at points with approximate gradient).

Let u GSBD p ( Ω ) . Let θ ( 0 , 1 ) . Then for L d -a.e. x 0 Ω there exists a family u ε GSBD p ( B ε ( x 0 ) ) such that

(5.1) { ( i ) u ε = u in a neighborhood of  B ε ( x 0 ) , lim ε 0 ε - ( d + 1 ) d ( { u ε u } ) = 0 , ( ii ) lim ε 0 ε - ( d + p ) B ( 1 - θ ) ε ( x 0 ) | u ε ( x ) - u ( x 0 ) - u ( x 0 ) ( x - x 0 ) | p d x = 0 , ( iii ) lim ε 0 ε - d B ε ( x 0 ) | e ( u ε ) ( x ) - e ( u ) ( x 0 ) | p d x = 0 , ( iv ) lim ε 0 ε - d d - 1 ( J u ε ) = 0 .

Proof.

Let x 0 Ω be such that

(5.2) { ( i ) lim ε 0 ε - d B ε ( x 0 ) | e ( u ) ( x ) - e ( u ) ( x 0 ) | p d x = 0 , ( ii ) lim ε 0 ε - d d - 1 ( J u B ε ( x 0 ) ) = 0 , ( iii ) lim ε 0 ε - d d ( { x B ε ( x 0 ) : | u ( x ) - u ( x 0 ) - u ( x 0 ) ( x - x 0 ) | | x - x 0 | > ϱ } ) = 0  for all  ϱ > 0 .

These properties hold for d -a.e. x 0 Ω by Lemma 3.5 and the facts that | e ( u ) | p L 1 ( Ω ) and J u is countably d - 1 -rectifiable. We use again the notation u ¯ x 0 bulk = x 0 , u 0 , u ( x 0 ) = u ( x 0 ) + u ( x 0 ) ( - x 0 ) for brevity, see (2.2).

Fix θ > 0 . We apply Theorem 3.3 and Remark 3.4 (i) for the function u - u ¯ x 0 bulk on the set B ( 1 - θ ) ε ( x 0 ) to obtain a set of finite perimeter ω ε B ( 1 - θ ) ε ( x 0 ) , a function v ε W 1 , p ( B ( 1 - θ ) ε ( x 0 ) ; d ) with v ε = u - u ¯ x 0 bulk in B ( 1 - θ ) ε ( x 0 ) ω ε , and an infinitesimal rigid motion a ε such that

(5.3) { ( i ) d - 1 ( * ω ε ) c d - 1 ( J u B ε ( x 0 ) ) , d ( ω ε ) c ( d - 1 ( J u B ε ( x 0 ) ) ) d d - 1 , ( ii ) v ε - a ε L p ( B ( 1 - θ ) ε ( x 0 ) ) c ε e ( u - u ¯ x 0 bulk ) L p ( B ε ( x 0 ) ) , ( iii ) e ( v ε ) L p ( B ( 1 - θ ) ε ( x 0 ) ) c e ( u - u ¯ x 0 bulk ) L p ( B ε ( x 0 ) ) ,

where c > 0 depends only on p, cf. also Remark 3.4 (iii). We directly note by (5.2) (ii) and (5.3) (i) that

(5.4) lim ε 0 ε - d 2 d - 1 d ( ω ε ) = 0 .

We define u ε GSBD p ( B ε ( x 0 ) ) as

(5.5) u ε := u χ B ε ( x 0 ) B ( 1 - θ ) ε ( x 0 ) + ( v ε + u ¯ x 0 bulk ) χ B ( 1 - θ ) ε ( x 0 ) ,

and proceed by confirming the properties stated in (5.1). Notice that, by construction, u ε = u in B ε ( x 0 ) ω ε . First, (5.1) (i) follows directly from the fact that ω ε B ( 1 - θ ) ε ( x 0 ) , as well as (5.4)–(5.5). Moreover, (5.3) (i) and (5.2) (ii) imply (5.1) (iv). As for (5.1) (iii), we notice that by (5.3) (iii) and (5.2) (i) we have

lim ε 0 ε - d B ( 1 - θ ) ε ( x 0 ) | e ( v ε ) ( x ) | p d x = 0 .

Since, by a direct computation, e ( u ε ) ( x ) - e ( u ) ( x 0 ) = e ( v ε ) ( x ) for x B ( 1 - θ ) ε ( x 0 ) , see (5.5), in combination with (5.2) (i) we obtain (5.1) (iii). It therefore remains to prove (5.1) (ii).

To this end, fix ϱ > 0 and define ω ^ ε := { x B ε ( x 0 ) : | u ( x ) - u ¯ x 0 bulk ( x ) | > ϱ ε } . In view of (5.2) (iii) and (5.4), we can choose ε 0 > 0 sufficiently small such that for all 0 < ε ε 0 we have

(5.6) d ( ω ε ω ^ ε ) 1 2 d ( B ( 1 - θ ) ε ( x 0 ) ) .

By the definition of ω ^ ε and the fact that v ε = u - u ¯ x 0 bulk in B ( 1 - θ ) ε ( x 0 ) ω ε , we have | v ε ( x ) | ϱ ε for all x B ( 1 - θ ) ε ( x 0 ) ( ω ε ω ^ ε ) . Hence, (5.3) (ii) and the triangle inequality give

a ε L p ( B ( 1 - θ ) ε ( x 0 ) ( ω ε ω ^ ε ) ) p C ε p e ( u - u ¯ x 0 bulk ) L p ( B ε ( x 0 ) ) p + C d ( B ε ( x 0 ) ) ϱ p ε p ,

where C > 0 depends only on p. By (5.6) and Lemma 3.6 we get

a ε L p ( B ( 1 - θ ) ε ( x 0 ) ) p C ε p e ( u - u ¯ x 0 bulk ) L p ( B ε ( x 0 ) ) p + C d ( B ε ( x 0 ) ) ϱ p ε p .

Therefore, by using also (5.2) (i), we derive lim sup ε 0 ε - ( d + p ) a ε L p ( B ( 1 - θ ) ε ( x 0 ) ) p C γ d ϱ p . As ϱ > 0 was arbitrary, we get

(5.7) lim ε 0 ε - ( d + p ) B ( 1 - θ ) ε ( x 0 ) | a ε | p d x = 0 .

Now, (5.2) (i) and (5.3) (ii) give that

lim ε 0 ε - ( d + p ) v ε - a ε L p ( B ( 1 - θ ) ε ( x 0 ) ) p c ε - d e ( u - u ¯ x 0 bulk ) L p ( B ε ( x 0 ) ) p = 0 .

As u ε - u ¯ x 0 bulk = v ε in B ( 1 - θ ) ε ( x 0 ) , this shows (5.1) (ii) by (5.7). ∎

We are now in a position to prove Lemma 4.2.

Proof of Lemma 4.2.

It suffices to prove (4.1) for points x 0 Ω where the statement of Lemma 5.1 holds and we have lim ε 0 ε - d μ ( B ε ( x 0 ) ) = γ d . This holds true for d -a.e. x 0 Ω . Then also lim ε 0 ε - d 𝐦 ( u , B ε ( x 0 ) ) exists, see Lemma 4.1. As before, we write u ¯ x 0 bulk = u ( x 0 ) + u ( x 0 ) ( - x 0 ) for shorthand.

Step 1 (Inequality “ ” in (4.1)). We fix η > 0 and θ > 0 . Choose z ε GSBD p ( B ( 1 - 3 θ ) ε ( x 0 ) ) with z ε = u ¯ x 0 bulk in a neighborhood of B ( 1 - 3 θ ) ε ( x 0 ) and

(5.8) ( z ε , B ( 1 - 3 θ ) ε ( x 0 ) ) 𝐦 ( u ¯ x 0 bulk , B ( 1 - 3 θ ) ε ( x 0 ) ) + ε d + 1 .

We extend z ε to a function in GSBD p ( B ε ( x 0 ) ) by setting z ε = u ¯ x 0 bulk outside B ( 1 - 3 θ ) ε ( x 0 ) . Let ( u ε ) ε be the family given by Lemma 5.1. We apply Lemma 3.8 on z ε (in place of u) and u ε (in place of v) for η as above and the sets

(5.9) A = B 1 - 2 θ ( x 0 ) , A = B 1 - θ ( x 0 ) , A ′′ = B 1 ( x 0 ) B 1 - 4 θ ( x 0 ) ¯ .

By (3.11)–(3.12) there exist functions w ε GSBD p ( B ε ( x 0 ) ) such that w ε = u ε on B ε ( x 0 ) B ( 1 - θ ) ε ( x 0 ) and

(5.10) ( w ε , B ε ( x 0 ) ) ( 1 + η ) ( ( z ε , A ε , x 0 ) + ( u ε , A ε , x 0 ′′ ) ) + M ε p z ε - u ε L p ( ( A A ) ε , x 0 ) p + d ( B ε ( x 0 ) ) η ,

where M > 0 depends on θ and η, but is independent of ε. Here and in the following, we use notation (3.1). In particular, we have w ε = u ε = u in a neighborhood of B ε ( x 0 ) by (5.1) (i). By (5.1) (ii), (5.9), and the fact that z ε = u ¯ x 0 bulk outside B ( 1 - 3 θ ) ε ( x 0 ) we find

(5.11) lim ε 0 ε - ( d + p ) z ε - u ε L p ( ( A A ) ε , x 0 ) p = lim ε 0 ε - ( d + p ) u ε - u ¯ x 0 bulk L p ( B ( 1 - θ ) ε ( x 0 ) ) p = 0 .

This along with (5.10) shows that there exists a sequence ( ρ ε ) ε ( 0 , + ) with ρ ε 0 such that

(5.12) ( w ε , B ε ( x 0 ) ) ( 1 + η ) ( ( z ε , A ε , x 0 ) + ( u ε , A ε , x 0 ′′ ) ) + ε d ρ ε + γ d ε d η .

On the one hand, by using that z ε = u ¯ x 0 bulk on B ε ( x 0 ) B ( 1 - 3 θ ) ε ( x 0 ) A ε , x 0 ′′ , (H1), (H4), and (5.8) we compute

lim sup ε 0 ( z ε , A ε , x 0 ) ε d lim sup ε 0 ( z ε , B ( 1 - 3 θ ) ε ( x 0 ) ) ε d + lim sup ε 0 ( u ¯ x 0 bulk , A ε , x 0 ′′ ) ε d
lim sup ε 0 𝐦 ( u ¯ x 0 bulk , B ( 1 - 3 θ ) ε ( x 0 ) ) ε d + β γ d [ 1 - ( 1 - 4 θ ) d ] ( 1 + | e ( u ) ( x 0 ) | p )
(5.13) ( 1 - 3 θ ) d lim sup ε 0 𝐦 ( u ¯ x 0 bulk , B ε ( x 0 ) ) ( ε ) d + β γ d [ 1 - ( 1 - 4 θ ) d ] ( 1 + | e ( u ) ( x 0 ) | p ) ,

where in the last step we substituted ( 1 - 3 θ ) ε by ε . On the other hand, by (H4) and (5.9) we also find

( u ε , A ε , x 0 ′′ ) β A ε , x 0 ′′ ( 1 + | e ( u ε ) | p ) + β d - 1 ( J u ε A ε , x 0 ′′ )
γ d ε d β [ 1 - ( 1 - 4 θ ) d ] ( 1 + 2 p - 1 | e ( u ) ( x 0 ) | p ) + 2 p - 1 β e ( u ε ) - e ( u ) ( x 0 ) L p ( B ε ( x 0 ) ) p + β d - 1 ( J u ε ) .

By (5.1) (iii)–(iv) this implies

(5.14) lim sup ε 0 ( u ε , A ε , x 0 ′′ ) ε d β γ d [ 1 - ( 1 - 4 θ ) d ] ( 1 + 2 p - 1 | e ( u ) ( x 0 ) | p ) .

Recall that w ε = u in a neighborhood of B ε ( x 0 ) . This along with (5.12)–(5.14) and ρ ε 0 yields

lim ε 0 𝐦 ( u , B ε ( x 0 ) ) γ d ε d lim sup ε 0 ( w ε , B ε ( x 0 ) ) γ d ε d
( 1 + η ) ( 1 - 3 θ ) d lim sup ε 0 𝐦 ( u ¯ x 0 bulk , B ε ( x 0 ) ) γ d ε d
+ 2 ( 1 + η ) β [ 1 - ( 1 - 4 θ ) d ] ( 1 + 2 p - 1 | e ( u ) ( x 0 ) | p ) + η .

Passing to η , θ 0 , we obtain inequality “ ” in (4.1).

Step 2 (Inequality “ ” in (4.1)). We fix η , θ > 0 and let ( u ε ) ε be again the family from Lemma 5.1. By (5.1) (i) and Fubini’s Theorem, for each ε > 0 we can find s ε ( 1 - 4 θ , 1 - 3 θ ) ε such that

(5.15) { ( i ) lim ε 0 ε - d d - 1 ( { u u ε } B s ε ( x 0 ) ) = 0 , ( ii ) d - 1 ( ( J u J u ε ) B s ε ( x 0 ) ) = 0 for all  ε > 0 .

We consider z ε GSBD p ( B s ε ( x 0 ) ) such that z ε = u in a neighborhood of B s ε ( x 0 ) , and

(5.16) ( z ε , B s ε ( x 0 ) ) 𝐦 ( u , B s ε ( x 0 ) ) + ε d + 1 .

We extend z ε to a function in GSBD p ( B ε ( x 0 ) ) by setting

(5.17) z ε = u ε in  B ε ( x 0 ) B s ε ( x 0 ) .

We apply Lemma 3.8 on z ε (in place of u) and u ¯ x 0 bulk (in place of v) for the sets indicated in equation (5.9). By (3.11)–(3.12) there exist functions w ε GSBD p ( B ε ( x 0 ) ) such that w ε = u ¯ x 0 bulk on B ε ( x 0 ) B ( 1 - θ ) ε ( x 0 ) and

( w ε , B ε ( x 0 ) ) ( 1 + η ) ( ( z ε , A ε , x 0 ) + ( u ¯ x 0 bulk , A ε , x 0 ′′ ) ) + M ε p z ε - u ¯ x 0 bulk L p ( ( A A ) ε , x 0 ) p + d ( B ε ( x 0 ) ) η .

By (5.17) and the choice of s ε we get that z ε = u ε outside B ( 1 - 3 θ ) ε ( x 0 ) . Thus, similar to Step 1, cf. (5.11) and (5.12), we find a sequence ( ρ ε ) ε ( 0 , + ) with ρ ε 0 such that

(5.18) ( w ε , B ε ( x 0 ) ) ( 1 + η ) ( ( z ε , A ε , x 0 ) + ( u ¯ x 0 bulk , A ε , x 0 ′′ ) ) + ε d ρ ε + γ d ε d η .

Let us now estimate the terms in (5.18). We get by (H1), (H4), (5.16)–(5.17), and the choice of s ε that

(5.19) ( z ε , A ε , x 0 ) 𝐦 ( u , B s ε ( x 0 ) ) + ε d + 1 + β d - 1 ( ( { u u ε } J u J u ε ) B s ε ( x 0 ) ) + ( u ε , A ε , x 0 ′′ ) .

Therefore, by (5.14), (5.15), and the fact that s ε ( 1 - 3 θ ) ε we derive

lim sup ε 0 ( z ε , A ε , x 0 ) ε d ( s ε ε ) d lim sup ε 0 𝐦 ( u , B s ε ( x 0 ) ) s ε d + β γ d [ 1 - ( 1 - 4 θ ) d ] ( 1 + 2 p - 1 | e ( u ) ( x 0 ) | p )
(5.20) ( 1 - 3 θ ) d lim sup ε 0 𝐦 ( u , B ε ( x 0 ) ) ε d + β γ d [ 1 - ( 1 - 4 θ ) d ] ( 1 + 2 p - 1 | e ( u ) ( x 0 ) | p ) .

Estimating ( u ¯ x 0 bulk , A ε , x 0 ′′ ) as in (5.13), with (5.18)–(5.20) and ρ ε 0 , we then obtain

lim sup ε 0 ( w ε , B ε ( x 0 ) ) ε d ( 1 + η ) ( 1 - 3 θ ) d lim sup ε 0 𝐦 ( u , B ε ( x 0 ) ) ε d
+ 2 ( 1 + η ) β γ d [ 1 - ( 1 - 4 θ ) d ] ( 1 + 2 p - 1 | e ( u ) ( x 0 ) | p ) + γ d η .

Passing to η , θ 0 and recalling that w ε = u ¯ x 0 bulk in a neighborhood of B ε ( x 0 ) , we derive

lim sup ε 0 𝐦 ( u ¯ x 0 bulk , B ε ( x 0 ) ) γ d ε d lim sup ε 0 ( w ε , B ε ( x 0 ) ) γ d ε d
lim sup ε 0 𝐦 ( u , B ε ( x 0 ) ) γ d ε d = lim ε 0 𝐦 ( u , B ε ( x 0 ) ) γ d ε d .

This shows inequality “ ” in (4.1) and concludes the proof. ∎

6 The surface density

This section is devoted to the proof of Lemma 4.3. We start by analyzing the blow-up at jump points. In the following, for any x 0 J u we adopt the notation u ¯ x 0 surf for the function u x 0 , u + ( x 0 ) , u - ( x 0 ) , ν u ( x 0 ) , see (2.3), with ν u ( x 0 ) 𝕊 d - 1 and u ± ( x 0 ) d being the approximate normal to J u and the traces on both sides of J u at x 0 , respectively. Recall also notation (3.1).

Lemma 6.1 (Blow-up at jump points).

Let u GSBD p ( Ω ) and let θ ( 0 , 1 ) . For H d - 1 -a.e. x 0 J u there exists a family u ε GSBD p ( B ε ( x 0 ) ) such that

(6.1) { ( i ) u ε = u in a neighborhood of  B ε ( x 0 ) lim ε 0 ε - d d ( { u ε u } ) = 0 , ( ii ) lim ε 0 ε - ( d - 1 + p ) B ( 1 - θ ) ε ( x 0 ) | u ε ( x ) - u ¯ x 0 surf | p d x = 0 , ( iii ) lim ε 0 ε - ( d - 1 ) d - 1 ( J u ε E ε , x 0 ) = d - 1 ( Π 0 E ) for all Borel sets  E B 1 ( x 0 ) , ( iv ) lim ε 0 ε - ( d - 1 ) B ε ( x 0 ) | e ( u ε ) ( x ) | p d x = 0 ,

where Π 0 denotes the hyperplane passing through x 0 with normal ν u ( x 0 ) .

Proof.

We start by using the fact that J u is countably d - 1 -rectifiable and the blow-up properties of GSBD p functions. Arguing as in, e.g., [24, proof of Theorem 2], [26, proof of Theorem 1.1], [28, Lemma 3.4], we infer that for d - 1 -a.e. x 0 J u there exist ε ¯ > 0 , ν u ( x 0 ) 𝕊 d - 1 , u ± ( x 0 ) d , and a hypersurface Γ which is a graph of a function h defined on Π 0 , being C 1 and Lipschitz, such that x 0 Γ , Π 0 is tangent to Γ in x 0 , Γ B ε ( x 0 ) separates B ε ( x 0 ) in two open connected components B ε Γ , ± ( x 0 ) for each ε < ε ¯ , and

(6.2) { ( i ) lim ε 0 ε - ( d - 1 ) d - 1 ( ( J u Γ ) B ε ( x 0 ) ) = 0 , lim ε 0 ε - ( d - 1 ) d - 1 ( Γ E ε , x 0 ) = d - 1 ( Π 0 E ) for all Borel sets  E B 1 ( x 0 ) , ( ii ) lim ε 0 ε - ( d - 1 ) B ε ( x 0 ) | e ( u ) | p d x = 0 , ( iii ) lim ε 0 ε - d d ( { x B ε ( x 0 ) : | u - u ¯ x 0 surf | > ϱ } ) = 0 for all  ϱ > 0 .

In particular, (ii) follows from the fact that | e ( u ) | p L 1 ( Ω ) and (iii) from (3.2). Then, since Π 0 is tangent to Γ in x 0 , Γ B ε ( x 0 ) is the graph of a Lipschitz function h ε defined on a subset of Π 0 , with Lipschitz constant L ε such that lim ε 0 L ε = 0 . Therefore, it holds that

(6.3) lim ε 0 ε - d d ( B ε Γ , ± ( x 0 ) B ε ± ( x 0 ) ) = 0 ,

where B ε ± ( x 0 ) := { y B ε ( x 0 ) : ± ( y - x 0 ) ν u ( x 0 ) > 0 } . By this and Fubini’s Theorem, for each ε > 0 we can find s ε ( 1 - θ , 1 - θ 2 ) ε such that

(6.4) lim ε 0 ε - ( d - 1 ) d ( ( B ε Γ , ± ( x 0 ) B ε ± ( x 0 ) ) B s ε ( x 0 ) ) = 0 .

For any ε > 0 , we apply Theorem 3.3 and Remark 3.4 (i) on u in the two connected components B s ε Γ , ± ( x 0 ) for ε < ε ¯ . This gives two functions v ε ± W 1 , p ( B s ε Γ , ± ( x 0 ) ; d ) , two sets of finite perimeter ω ε ± B s ε Γ , ± ( x 0 ) , and two infinitesimal rigid motions a ε ± such that

(6.5) { ( i ) v ε ± = u in  B s ε Γ , ± ( x 0 ) ω ε ± , ( ii ) d - 1 ( * ω ε ± ) c d - 1 ( J u B ε Γ , ± ( x 0 ) ) , d ( ω ε ± ) c ( d - 1 ( J u B ε Γ , ± ( x 0 ) ) ) d d - 1 , ( iii ) v ε ± - a ε ± L p ( B s ε Γ , ± ( x 0 ) ) c ε e ( u ) L p ( B ε Γ , ± ( x 0 ) ) , ( iv ) v ε ± - a ε ± L p ( B s ε Γ , ± ( x 0 ) ) c e ( u ) L p ( B ε Γ , ± ( x 0 ) ) ,

where c > 0 is independent of ε. (See Remark 3.4 (ii)–(iii) and recall that the Lipschitz constant of h ε vanishes as ε 0 .) By the Sobolev extension theorem we extend v ε ± to v ^ ε ± W 1 , p ( B s ε ( x 0 ) ; d ) , and (6.5) (iii)–(iv) along with the linearity of the extension operator yield

(6.6) ε - 1 v ^ ε ± - a ε ± L p ( B s ε ( x 0 ) ) + v ^ ε ± - a ε ± L p ( B s ε ( x 0 ) ) c e ( u ) L p ( B ε Γ , ± ( x 0 ) ) ,

where, as before, the constant is independent of ε. (Here, we used again the properties of the functions h ε recalled below (6.2).) We define u ε GSBD p ( B ε ( x 0 ) ) as

(6.7) u ε := { v ^ ε + in  B s ε + ( x 0 ) , v ^ ε - in  B s ε - ( x 0 ) , u in  B ε ( x 0 ) B s ε ( x 0 ) ,

where B s ε ± ( x 0 ) is defined below (6.3). We now prove the properties in (6.1). First, by definition we have that u ε = u in a neighborhood of B ε ( x 0 ) . By B ( 1 - θ ) ε Γ , ± ( x 0 ) Γ = , (6.2) (i), (6.3), and (6.5) (i)–(ii) we obtain lim ε 0 ε - d d ( { u ε u } ) = 0 . This concludes (6.1) (i). Moreover, (6.2) (ii) and (6.6) imply (6.1) (iv). By the definition of u ε and (6.5) (i) it holds that

J u ε ( Π 0 B s ε ( x 0 ) ¯ ) ( J u ( B ε ( x 0 ) B s ε ( x 0 ) ) ) * ω ε + * ω ε - ( ( B ε Γ , ± ( x 0 ) B ε ± ( x 0 ) ) B s ε ( x 0 ) ) .

We now show (6.1) (iii). Concerning the ” ” inequality, for a fixed Borel set E B 1 ( x 0 ) we have to estimate the measure of the intersection with E ε , x 0 and any of the five sets in the right-hand side above: it holds that

d - 1 ( Π 0 B s ε ( x 0 ) ¯ E ε , x 0 ) = ε d - 1 d - 1 ( Π 0 B s ε / ε ( x 0 ) ¯ E )

for any ε > 0 by rescaling, that

lim ε 0 | ε - ( d - 1 ) d - 1 ( J u ( E ε , x 0 B s ε ( x 0 ) ) ) - d - 1 ( Π 0 ( E B s ε / ε ( x 0 ) ) ) | = 0

by (6.2) (i), while the last three terms are estimated by (6.5) (ii) and (6.4). To see the converse, we first apply [25, Theorem 1.1] to the functions u ε ( x 0 + ε ) , which converge in measure to u ¯ x 0 surf in B 1 ( 0 ) by (6.2), (6.3), and (6.5) (ii). Then we scale back to B ε ( x 0 ) . Hence, (6.1) (iii) holds.

It remains to prove (6.1) (ii). We notice that this easily follows from

(6.8) lim ε 0 ε - ( d - 1 + p ) B s ε ± ( x 0 ) | a ε ± - u ¯ x 0 surf | p d x = 0 .

In fact, (6.2) (ii) and (6.6) give that

(6.9) lim ε 0 ε - ( d - 1 + p ) B s ε ± ( x 0 ) | v ^ ε ± - a ε ± | p d x = 0 .

Then (6.8), (6.9), the triangle inequality, s ε ( 1 - θ ) ε , and (6.7) imply (6.1) (ii).

Therefore, let us now confirm (6.8). We only address the “ + ” case, for the “ - ” case is analogous. We first observe that by (6.2) (iii), (6.3), and a diagonal argument, we may find a sequence ( ϱ ε ) ε ( 0 , + ) with lim ε 0 ϱ ε = 0 such that the sets

ω ^ ε + := { x B ε ( x 0 ) : | u ( x ) - u + ( x 0 ) | > ϱ ε } B s ε Γ , + ( x 0 )

satisfy

(6.10) lim ε 0 ε - d d ( ω ^ ε + ) = 0 .

In view of (6.5) (i), (iii), we have that

(6.11) u - a ε + L p ( B s ε Γ , + ( x 0 ) ω ε + ) c ε e ( u ) L p ( B ε Γ , + ( x 0 ) ) .

Then, by (6.11), the definition of ω ^ ε + , and the triangle inequality we get that

(6.12) u + ( x 0 ) - a ε + L p ( B s ε Γ , + ( x 0 ) ( ω ε + ω ^ ε + ) ) c ε e ( u ) L p ( B ε Γ , + ( x 0 ) ) + γ d 1 p ε d p ϱ ε .

By (6.2) (i), (6.3), (6.5) (ii), and (6.10) we obtain d ( ω ε + ω ^ ε + ) 1 2 d ( B s ε Γ , + ( x 0 ) ) for ε sufficiently small. Then, by Lemma 3.6, we have that γ d ε d p u + ( x 0 ) - a ε + L ( B ε ( x 0 ) ) is less or equal than the right-hand side of (6.12), up to multiplication with a constant. This along with (6.2) (ii), p 1 , and the fact that ϱ ε 0 implies

(6.13) lim ε 0 u + ( x 0 ) - a ε + L ( B ε ( x 0 ) ) = 0 .

Let us consider A ε + 𝕄 skew d × d and b ε + d such that a ε + ( x ) = A ε + ( x - x 0 ) + b ε + . Then (6.8) follows by

(6.14)

(6.14a) lim ε 0 ε | A ε + | p = 0 ,
(6.14b) lim ε 0 ε 1 - p p | b ε + - u + ( x 0 ) | = 0 .

So we are left to prove (6.14) which corresponds to [28, equations (3.18)–(3.19)]. The proof goes in the same way with slight modifications that we indicate below. For fixed δ > 0 small, by (6.2) (ii) there exists ε ^ > 0 , depending on δ, such that

(6.15) ε - ( d - 1 ) B ε ( x 0 ) | e ( u ) | p d x δ p for all  ε ε ^ .

For ε ~ < ε < ε ^ , we set ε k := min { 2 k ε ~ , ε } and adopt the notation k in place of ε k in the subscripts. We then obtain

(6.16) a k + - a k + 1 + L ( B ε k + ( x 0 ) ) c γ d - 1 p ε k - d p a k + - a k + 1 + L p ( B s ε k Γ , + ( x 0 ) ( ω k + ω k + 1 + ) ) c δ ε k p - 1 p .

In fact, the first inequality follows from (6.5) (ii) and Lemma 3.6, and the second one from (6.11), (6.15), and the triangle inequality. Similarly, employing (6.5) (i), (iv) in place of (6.11), and recalling a ε k + = A k + , we obtain

(6.17) | A k + - A k + 1 + | c δ ε k - 1 p .

At this stage, (6.14) follow exactly as in [28]: for k ^ being the first index such that ε k ^ = ε , recalling ε ~ = ε 0 and summing (6.17) gives

ε ~ | A ε ~ + | p ε ~ ( | A k ^ + | + k = 0 k ^ - 1 | A k + - A k + 1 + | ) p c δ p + c ε ~ | A k ^ + | p .

The right-hand side vanishes as ε ~ 0 and δ 0 , and this proves (6.14a). Moreover, summing (6.16) (and since | b k + - b k + 1 + | a k + - a k + 1 + L ( B ε k + ( x 0 ) ) ), we obtain

| b ε ~ + - b ε + | c δ ε p - 1 p

for all 0 < ε ~ ε . By passing to the limit as ε ~ 0 together with (6.13), we get

ε 1 - p p | u + ( x 0 ) - b ε + | c δ

for ε < ε ^ = ε ^ ( δ ) . Thus, (6.14b) follows by the arbitrariness of δ > 0 , concluding the proof. ∎

Remark 6.2 (Construction of u ε ).

We point out that our definition of u ε in (6.7) differs from the corresponding constructions in [13, Lemma 3] and [28, Lemma 3.4] in order to fix a possible flaw contained in these proofs. Roughly speaking, in our notation, in [13, 28], u ε on B s ε ( x 0 ) is defined as

(6.18) u ε = { v ε + in  B s ε Γ , + ( x 0 ) , v ε - in  B s ε Γ , - ( x 0 ) .

Then, instead of (6.8), one needs to check

lim ε 0 ε - ( d - 1 + p ) B s ε Γ ± ( x 0 ) | a ε ± - u ¯ x 0 surf | p d x = 0 .

This, however, is in general false if lim inf ε 0 ε - ( d - 1 + p ) d ( B ε Γ , ± ( x 0 ) B ε ± ( x 0 ) ) > 0 (which is clearly possible). Let us also remark that, in contrast to our construction, (6.18) allows to prove an estimate of the form

(6.19) lim ε 0 ε - ( d - 1 ) d - 1 ( J u ε J u ) = 0 ,

see [13, equation (24)] and [28, Lemma 3.4 (i)]. It is not clear to us if it is possible that for u ε satisfying the fundamental blow-up property (6.1) (ii) one may still have an estimate of the form (6.19). The latter, however, is not needed for our proofs.

We now proceed with the proof of Lemma 4.3.

Proof of Lemma 4.3.

The proof follows the same strategy of the proof of Lemma 4.2. We fix x 0 J u such that the statement of Lemma 6.1 holds at x 0 and lim ε 0 ε - ( d - 1 ) μ ( B ε ( x 0 ) ) = γ d - 1 . This is possible for d - 1 -a.e. x 0 J u . Then also lim ε 0 ε - ( d - 1 ) 𝐦 ( u , B ε ( x 0 ) ) exists, see Lemma 4.1. We prove (4.2) for x 0 of this type.

Step 1 (Inequality “ ” in (4.2)). We fix η , θ > 0 , and consider z ε GSBD p ( B ( 1 - 3 θ ) ε ( x 0 ) ) with z ε = u ¯ x 0 surf in a neighborhood of B ( 1 - 3 θ ) ε ( x 0 ) and

(6.20) ( z ε , B ( 1 - 3 θ ) ε ( x 0 ) ) 𝐦 ( u ¯ x 0 surf , B ( 1 - 3 θ ) ε ( x 0 ) ) + ε d .

We extend z ε to a function in GSBD p ( B ε ( x 0 ) ) by setting z ε = u ¯ x 0 surf outside B ( 1 - 3 θ ) ε ( x 0 ) . Let ( u ε ) ε be the family given by Lemma 6.1. As in the proof of Lemma 4.2, we apply Lemma 3.8 on z ε (in place of u) and u ε (in place of v) for η fixed above and the sets

A = B 1 - 2 θ ( x 0 ) , A = B 1 - θ ( x 0 ) , A ′′ = B 1 ( x 0 ) B 1 - 4 θ ( x 0 ) ¯ .

Recalling notation (3.1) and (3.11)–(3.12), we find a function w ε GSBD p ( B ε ( x 0 ) ) such that w ε = u ε on B ε ( x 0 ) B ( 1 - θ ) ε ( x 0 ) and

(6.21) ( w ε , B ε ( x 0 ) ) ( 1 + η ) ( ( z ε , A ε , x 0 ) + ( u ε , A ε , x 0 ′′ ) ) + M ε p z ε - u ε L p ( ( A A ) ε , x 0 ) p + d ( B ε ( x 0 ) ) η ,

where M > 0 depends on θ and η, but is independent of ε. In particular, we have w ε = u ε = u in a neighborhood of B ε ( x 0 ) by (6.1) (i). By (6.1) (ii) and the fact that z ε = u ¯ x 0 surf outside B ( 1 - 3 θ ) ε ( x 0 ) we find

lim ε 0 ε - ( d - 1 + p ) z ε - u ε L p ( ( A A ) ε , x 0 ) p = lim ε 0 ε - ( d - 1 + p ) u ε - u ¯ x 0 surf L p ( B ( 1 - θ ) ε ( x 0 ) ) p = 0 .

Inserting this in (6.21), we find that, for a suitable sequence ( ρ ε ) ε ( 0 , + ) with ρ ε 0 ,

(6.22) ( w ε , B ε ( x 0 ) ) ( 1 + η ) ( ( z ε , A ε , x 0 ) + ( u ε , A ε , x 0 ′′ ) ) + ε d - 1 ρ ε + γ d ε d η .

We evaluate the first terms on the right-hand side of (6.22): since z ε = u ¯ x 0 bulk on B ε ( x 0 ) B ( 1 - 3 θ ) ε ( x 0 ) A ε , x 0 ′′ , by (H1), (H4), and (6.20) we have

lim sup ε 0 ( z ε , A ε , x 0 ) ε d - 1 lim sup ε 0 ( z ε , B ( 1 - 3 θ ) ε ( x 0 ) ) ε d - 1 + lim sup ε 0 ( u ¯ x 0 surf , A ε , x 0 ′′ ) ε d - 1
lim sup ε 0 𝐦 ( u ¯ x 0 surf , B ( 1 - 3 θ ) ε ( x 0 ) ) ε d - 1 + β d - 1 ( A ′′ Π 0 )
(6.23) ( 1 - 3 θ ) d - 1 lim sup ε 0 𝐦 ( u ¯ x 0 surf , B ε ( x 0 ) ) ε d - 1 + β γ d - 1 ( 1 - ( 1 - 4 θ ) d - 1 ) ,

where, as in Lemma 6.1, we denote by Π 0 the hyperplane passing through x 0 with normal ν u ( x 0 ) . By (H4) and (6.1) (iii), (iv) we get

lim sup ε 0 ( u ε , A ε , x 0 ′′ ) β ε d - 1 lim sup ε 0 A ε , x 0 ′′ ( 1 + | e ( u ε ) | p ) d x ε d - 1 + lim sup ε 0 d - 1 ( J u ε A ε , x 0 ′′ ) ε d - 1
(6.24) = γ d - 1 ( 1 - ( 1 - 4 θ ) d - 1 ) .

Collecting (6.22), (6.23), (6.24), and recalling ρ ε 0 , as well as the fact that w ε = u in a neighborhood of B ε ( x 0 ) , we obtain

lim ε 0 𝐦 ( u , B ε ( x 0 ) ) γ d - 1 ε d - 1 lim sup ε 0 ( w ε , B ε ( x 0 ) ) γ d - 1 ε d - 1 ( 1 + η ) ( 1 - 3 θ ) d - 1 lim sup ε 0 𝐦 ( u ¯ x 0 surf , B ε ( x 0 ) ) γ d - 1 ε d - 1 + 2 ( 1 + η ) β ( 1 - ( 1 - 4 θ ) d - 1 ) .

Passing to the limit as η , θ 0 , we conclude inequality “ ” in (4.2).

Step 2 (Inequality “ ” in (4.2)). We fix η , θ > 0 and let, as in Step 1, ( u ε ) ε be the family given by Lemma 6.1. By (6.1) (i) and Fubini’s Theorem, for each ε > 0 we can find s ε ( 1 - 4 θ , 1 - 3 θ ) ε such that

(6.25) { ( i ) lim ε 0 ε - ( d - 1 ) d - 1 ( { u u ε } B s ε ( x 0 ) ) = 0 , ( ii ) d - 1 ( ( J u J u ε ) B s ε ( x 0 ) ) = 0 for all  ε > 0 .

We consider z ε GSBD p ( B s ε ( x 0 ) ) such that z ε = u in a neighborhood of B s ε ( x 0 ) , and

(6.26) ( z ε , B s ε ( x 0 ) ) 𝐦 ( u , B s ε ( x 0 ) ) + ε d .

We extend z ε to a function in GSBD p ( B ε ( x 0 ) ) by setting

(6.27) z ε = u ε in  B ε ( x 0 ) B s ε ( x 0 ) .

We apply Lemma 3.8 for z ε (in place of u), u ¯ x 0 surf (in place of v), and for the sets A, A , B as in Step 1, in correspondence to ε. By (3.11)–(3.12), there exists w ε GSBD p ( B ε ( x 0 ) ) such that w ε = u ¯ x 0 surf on B ε ( x 0 ) B ( 1 - θ ) ε ( x 0 ) , and

( w ε , B ε ( x 0 ) ) ( 1 + η ) ( ( z ε , A ε , x 0 ) + ( u ¯ x 0 surf , A ε , x 0 ′′ ) ) + M ε p z ε - u ¯ x 0 surf L p ( ( A A ) ε , x 0 ) p + d ( B ε ( x 0 ) ) η .

We observe that z ε = u ε outside B ( 1 - 3 θ ) ε ( x 0 ) , by (6.27) and the choice of s ε . Then, as done in Step 1, we may employ (6.1) (ii). This gives us a sequence ( ρ ε ) ε ( 0 , + ) with ρ ε 0 as ε 0 such that

(6.28) ( w ε , B ε ( x 0 ) ) ( 1 + η ) ( ( z ε , A ε , x 0 ) + ( u ¯ x 0 surf , A ε , x 0 ′′ ) ) + ε d - 1 ρ ε + γ d ε d η .

We estimate the first terms in (6.28). We get by (H1), (H4), (6.26)–(6.27), and the choice of s ε that

(6.29) ( z ε , A ε , x 0 ) 𝐦 ( u , B s ε ( x 0 ) ) + ε d + β d - 1 ( ( { u u ε } J u J u ε ) B s ε ( x 0 ) ) + ( u ε , A ε , x 0 ′′ ) .

By (6.24), (6.25), and the fact that s ε ( 1 - 3 θ ) ε we thus deduce that

lim sup ε 0 ( z ε , A ε , x 0 ) ε d - 1 lim sup ε 0 𝐦 ( u , B s ε ( x 0 ) ) ε d - 1 + β γ d - 1 ( 1 - ( 1 - 4 θ ) d - 1 )
(6.30) ( 1 - 3 θ ) d - 1 lim sup ε 0 𝐦 ( u , B ε ( x 0 ) ) ε d - 1 + β γ d - 1 ( 1 - ( 1 - 4 θ ) d - 1 ) ,

and, similarly to (6.23),

(6.31) lim sup ε 0 ( u ¯ x 0 surf , A ε , x 0 ′′ ) ε d - 1 β γ d - 1 ( 1 - ( 1 - 4 θ ) d - 1 ) .

Collecting (6.28), (6.30), (6.31) and using ρ ε 0 , we derive

lim sup ε 0 ( w ε , B ε ( x 0 ) ) ε d - 1 ( 1 + η ) ( ( 1 - 3 θ ) d - 1 lim sup ε 0 𝐦 ( u , B ε ( x 0 ) ) ε d - 1 + 2 β γ d - 1 ( 1 - ( 1 - 4 θ ) d - 1 ) ) .

Finally, recalling that w ε = u ¯ x 0 surf in a neighborhood of B ε ( x 0 ) , and using the arbitrariness of η, θ > 0 , we obtain

lim sup ε 0 𝐦 ( u ¯ x 0 surf , B ε ( x 0 ) ) γ d - 1 ε d - 1 lim sup ε 0 ( w ε , B ε ( x 0 ) ) γ d - 1 ε d - 1 lim sup ε 0 𝐦 ( u , B ε ( x 0 ) ) γ d - 1 ε d - 1 = lim ε 0 𝐦 ( u , B ε ( x 0 ) ) γ d - 1 ε d - 1 .

This shows “ ” in (4.2) and concludes the proof. ∎

7 The SBD p case

This section is devoted to the analysis of the integral representation result for : SBD p ( Ω ) × ( Ω ) [ 0 , + ) satisfying (H1)(H3) and (H4’). This case has been addressed, for d = 2 , in [28]. On the one hand, the arguments there could be now generalized to general dimension by virtue of Theorem 3.3. On the other hand, as we are going to show, the result can also be obtained with minor changes of our more general strategy.

We start by pointing out that, under (H4’), only competitors in SBD p may have finite energy. In fact, in view of Proposition 3.2, in the present setting definition (2.1) reads as

(7.1) 𝐦 ( u , A ) = inf v SBD p ( Ω ) { ( v , A ) : v = u  in a neighborhood of  A } .

Then the following integral representation result holds.

Theorem 7.1 (Integral representation in SBD p ).

Let Ω R d be open, bounded with Lipschitz boundary and suppose that F : SBD p ( Ω ) × B ( Ω ) [ 0 , + ) satisfies (H1)(H3) and (H4’). Then

( u , B ) = B f ( x , u ( x ) , u ( x ) ) d x + J u B g ( x , u + ( x ) , u - ( x ) , ν u ( x ) ) d d - 1 ( x )

for all u SBD p ( Ω ) and B B ( Ω ) , where f is given by

f ( x 0 , u 0 , ξ ) = lim sup ε 0 𝐦 ( x 0 , u 0 , ξ , B ε ( x 0 ) ) γ d ε d

for all x 0 Ω , u 0 R d , ξ M d × d , and x 0 , u 0 , ξ as in (2.2), and g is given by

g ( x 0 , a , b , ν ) = lim sup ε 0 𝐦 ( u x 0 , a , b , ν , B ε ( x 0 ) ) γ d - 1 ε d - 1

for all x 0 Ω , a , b R d , ν S d - 1 , and u x 0 , a , b , ν as in (2.3).

The remainder of this section is devoted to the proof of Theorem 7.1 which follows along the lines of the proof of Theorem 2.1 devised in Section 4. First, the analogue of Lemma 4.1 holds essentially with the same proof.

Lemma 7.2.

Suppose that F satisfies (H1)(H3) and (H4’). Let u SBD p ( Ω ) and let μ = L d Ω + H d - 1 J u Ω . Then for μ-a.e. x 0 Ω we have

lim ε 0 ( u , B ε ( x 0 ) ) μ ( B ε ( x 0 ) ) = lim ε 0 𝐦 ( u , B ε ( x 0 ) ) μ ( B ε ( x 0 ) ) .

Proof.

One can follow the same argument used to prove Lemma 4.1 through Lemma 4.4.

First, we remark that [28, Lemmas 4.2 and 4.3] is proved under the assumptions (H4’) and u SBD p ( Ω ) , hence it can be used to derive (4.15).

Concerning Lemma 4.4, as the lower bound of (H4’) is stronger than the one of (H4), the GSBD compactness result [25, Theorem 1.1] is still applicable. First, this shows v δ GSBD p ( Ω ) . Additionally, (4.9), (H4’), and Proposition 3.2 imply that the function v δ belongs indeed to SBD p ( Ω ) . The rest of the proof remains unchanged, upon noticing that, under assumption (H4’), (4.10) still holds. ∎

We now address the adaptions necessary for the bulk density. When u SBD p ( Ω ) , we show that the approximating sequence constructed in Lemma 5.1 satisfies some additional properties.

Lemma 7.3.

Let u SBD p ( Ω ) . Let θ ( 0 , 1 ) . For L d -a.e. x 0 Ω there exists a family u ε SBD p ( B ε ( x 0 ) ) such that (5.1) holds, and additionally

(7.2) { ( i ) lim ε 0 ε - ( d + 1 ) B ε ( x 0 ) | u ε - u | d x = 0 , ( ii ) lim ε 0 ε - d J u ε | [ u ε ] | d d - 1 = 0 .

Proof.

Since u SBD ( Ω ) , for d -a.e. x 0 it holds that

(7.3) lim ε 0 ε - d J u B ε ( x 0 ) | [ u ] | d d - 1 = 0 ,

and (see [4, Theorem 7.4]) that

(7.4) lim ε 0 ε - ( d + 1 ) B ε ( x 0 ) | u - u ¯ x 0 bulk | d x = 0 ,

where for brevity we again let u ¯ x 0 bulk = x 0 , u ( x 0 ) , u ( x 0 ) . Hence, with Fubini’s Theorem we can fix s ε ( 1 - θ , 1 - θ 2 ) ε so that d - 1 ( J u B s ε ( x 0 ) ) = 0 and

(7.5) lim ε 0 ε - d B s ε ( x 0 ) | u - u ¯ x 0 bulk | d d - 1 = 0 .

We can now perform the same construction as in (5.5) with s ε in place of ( 1 - θ ) ε . Notice that in this case u ε SBD p ( B ε ( x 0 ) ) . By arguing exactly as in the proof of Lemma 5.1, we derive (5.1) (i), (iii), (iv) while (ii) holds in B s ε ( x 0 ) and a fortiori in B ( 1 - θ ) ε ( x 0 ) . In particular, this in combination with Hölder’s inequality, (7.4), and u = u ε in B ε ( x 0 ) B s ε ( x 0 ) yields (7.2) (i).

To see (7.2) (ii), observe that, since s ε ε is bounded from above and from below, (5.1) (ii), (iii), the fact that u ε B s ε ( x 0 ) W 1 , p ( B s ε ( x 0 ) ; d ) (see (5.5)), and Korn’s inequality imply that

u ε ( x 0 + s ε ) - u ¯ x 0 bulk ( x 0 + s ε ) s ε 0 in  W 1 , p ( B 1 ( 0 ) ; d ) .

Hence, by the trace inequality and by scaling back to B s ε ( x 0 ) , we obtain

lim ε 0 s ε - ( d - 1 ) B s ε ( x 0 ) 1 s ε | u ε - u ¯ x 0 bulk | d d - 1 = 0 .

With this, (7.5), and the fact that s ε ε is bounded from below we then have

lim ε 0 ε - d B s ε ( x 0 ) | u ε - u | d d - 1 = 0 .

Hence, we get by the construction of u ε and (7.3) that

lim ε 0 ε - d J u ε | [ u ε ] | d d - 1 = lim ε 0 ε - d J u ( B ε ( x 0 ) B s ε ( x 0 ) ) | [ u ] | d d - 1 = 0 ,

which concludes the proof. ∎

With the above lemma at our disposal, we can deduce the asymptotic equivalence of the minimization problems (7.1) for u and u ¯ x 0 bulk = x 0 , u ( x 0 ) , u ( x 0 ) .

Lemma 7.4.

Suppose that F satisfies (H1), (H3), and (H4’). Let u SBD p ( Ω ) . Then for L d -a.e. x 0 Ω we have

(7.6) lim ε 0 𝐦 ( u , B ε ( x 0 ) ) γ d ε d = lim sup ε 0 𝐦 ( u ¯ x 0 bulk , B ε ( x 0 ) ) γ d ε d .

Proof.

We argue as in the proof of Lemma 4.2.

For the ” ” inequality, take u ε satisfying (5.1) and (7.2), and perform the same construction as in Lemma 4.2. (Observe that the fundamental estimate also holds in this case, see Lemma 3.8.) Notice that, in this case, we have by (H4’) and (5.9) that

( u ε , A ε , x 0 ′′ ) β A ε , x 0 ′′ ( 1 + | e ( u ε ) | p ) d x + β J u ε A ε , x 0 ′′ ( 1 + | [ u ε ] | ) d d - 1 .

Thus, using (5.1) (iii), (iv) and (7.2) (ii), we still get (5.14), and may deduce inequality “ ” in (7.6).

For the “ ” inequality, we start by observing that, if u ε also satisfies (7.2), in addition to (5.15) we may require that

(7.7) lim ε 0 ε - d B s ε ( x 0 ) | u + - u ε - | d d - 1 = 0 ,

where u ε - and u + indicate the inner and outer traces at B s ε ( x 0 ) , respectively. We also have that (5.14) holds, as seen in the previous step. We then perform the same construction as in Lemma 4.2. In this case, inequality (5.19) is replaced by

( z ε , A ε , x 0 ) 𝐦 ( u , B s ε ( x 0 ) ) + ε d + 1 + β d - 1 ( ( { u u ε } J u J u ε ) B s ε ( x 0 ) ) + β B s ε ( x 0 ) | u + - u ε - | d d - 1 + ( u ε , A ε , x 0 ′′ ) ,

so that, using (5.14), (5.15), (7.7), and the fact that s ε ( 1 - 3 θ ) ε , we are still in a position to deduce (5.20). The rest of the argument remains unchanged and we obtain inequality “ ” in (7.6). ∎

Similar changes have to be performed also for the surface density. We first deduce the analogue of Lemma 6.1. We again set u ¯ x 0 surf = u x 0 , u + ( x 0 ) , u - ( x 0 ) , ν u ( x 0 ) for brevity, see (2.3). We also recall the notation in (3.1).

Lemma 7.5.

Let u SBD p ( Ω ) and θ ( 0 , 1 ) . For H d - 1 -a.e. x 0 J u there exists a family u ε SBD p ( B ε ( x 0 ) ) such that (6.1) holds, and additionally

(7.8) { ( i ) lim ε 0 ε - d B ε ( x 0 ) | u ε - u | d x = 0 , ( ii ) lim ε 0 ε - ( d - 1 ) J u ε E ε , x 0 | [ u ε ] | d d - 1 = | [ u ¯ x 0 surf ] | d - 1 ( Π 0 E ) for all Borel sets  E B 1 ( x 0 ) ,

where Π 0 denotes the hyperplane passing through x 0 with normal ν u ( x 0 ) .

Proof.

Since u SBD ( Ω ) , for d - 1 -a.e. x 0 J u it holds that

(7.9) lim ε 0 ε - d B ε ( x 0 ) | u - u ¯ x 0 surf | d x = 0 .

Hence, by Fubini’s Theorem we find s ε ( 1 - θ , 1 - θ 2 ) ε such that (6.4) holds, we have

d - 1 ( J u B s ε ( x 0 ) ) = 0

and

(7.10) lim ε 0 ε - ( d - 1 ) B s ε ( x 0 ) | u - u ¯ x 0 surf | d d - 1 = 0 .

We can now perform the same construction as in (6.7) and derive (6.1). In particular, this in combination with Hölder’s inequality, p > 1 , (7.9), and u = u ε in B ε ( x 0 ) B s ε ( x 0 ) yields (7.8) (i).

To see (7.8) (ii), observe that, since s ε ε is bounded from above and from below, (6.1) (ii), (iv), p > 1 , the fact that u ε B s ε ± ( x 0 ) W 1 , p ( B s ε ± ( x 0 ) ; d ) (see (6.7)), and Korn’s inequality imply that

s ε - 1 ( B s ε ± ( x 0 ) - x 0 ) | u ε ( x 0 + s ε y ) - u ¯ x 0 surf | p d y + s ε - 1 ( B s ε ± ( x 0 ) - x 0 ) | y u ε ( x 0 + s ε y ) | p d y 0 .

Hence, by the trace inequality and by scaling back to B s ε ± ( x 0 ) , we obtain

lim ε 0 s ε - ( d - 1 ) B s ε ± ( x 0 ) | u ε - u ± ( x 0 ) | d d - 1 = 0 .

Since s ε ε is bounded from below, we then get that

(7.11) lim ε 0 ε - ( d - 1 ) B s ε ± ( x 0 ) | u ε - u ± ( x 0 ) | d d - 1 = 0 .

Given a Borel set E B 1 ( x 0 ) , we define E ε = E B s ε ε ( x 0 ) for every ε > 0 . Then by (6.7) and (7.11) we get

(7.12) lim ε 0 ( ε - ( d - 1 ) J u ε E ε , x 0 ε | [ u ε ] | d d - 1 - | [ u ¯ x 0 surf ] | d - 1 ( Π 0 E ε ) ) = 0 .

By (7.10) and (7.11) we also have

lim ε 0 ε - ( d - 1 ) B s ε ( x 0 ) | u ε - - u + | d d - 1 = 0 ,

where u ε - and u + indicate the inner and outer traces at B s ε ( x 0 ) , respectively. Hence, by construction of u ε in (6.7) and since u SBD ( Ω ) , we obtain

lim ε 0 ( ε - ( d - 1 ) J u ε ( E E ε ) ε , x 0 | [ u ε ] | d d - 1 - | [ u ¯ x 0 surf ] | d - 1 ( Π 0 ( E E ε ) ) ) = lim ε 0 ( ε - ( d - 1 ) J u ( E E ε ) ε , x 0 | [ u ] | d d - 1 - | [ u ¯ x 0 surf ] | d - 1 ( Π 0 ( E E ε ) ) ) = 0 .

By combining with (7.12), this concludes the proof of (7.8) (ii). ∎

With this lemma at hand, we can address the equivalence of minimization problems for the surface scaling.

Lemma 7.6.

Suppose that F satisfies (H1), (H3), and (H4’). Let u SBD p ( Ω ) . Then for H d - 1 -a.e. x 0 J u we have

(7.13) lim ε 0 𝐦 ( u , B ε ( x 0 ) ) γ d - 1 ε d - 1 = lim sup ε 0 𝐦 ( u ¯ x 0 surf , B ε ( x 0 ) ) γ d - 1 ε d - 1 .

Proof.

We argue as in the proof of Lemma 4.3.

For the “ ” inequality, take u ε satisfying (6.1) and (7.8), and perform the same construction as in Lemma 4.3. Estimates (6.22) and (6.23) continue to hold, provided one replaces β with the larger constant β ( 1 + | [ u ¯ x 0 surf ] | ) . Then, with (H4’), (6.1) (iii), (iv), and (7.8) (ii) we get

(7.14) lim sup ε 0 ( u ε , A ε , x 0 ′′ ) ε d - 1 β ( 1 + | [ u ¯ x 0 surf ] | ) γ d - 1 ( 1 - ( 1 - 4 θ ) d - 1 ) ,

which is the analogue of (6.24). This is enough to derive inequality “ ” in (7.13).

For the reverse one, take again u ε satisfying (6.1) and (7.8). Then, by (7.8) (i), in addition to (6.25) we may also require that

(7.15) lim ε 0 ε - ( d - 1 ) B s ε ( x 0 ) | u + - u ε - | d d - 1 = 0 ,

where u ε - and u + indicate the inner and outer traces at B s ε ( x 0 ) , respectively. We perform the same construction as in Lemma 4.3. In this case, inequality (6.29) is replaced by

( z ε , A ε , x 0 ) 𝐦 ( u , B s ε ( x 0 ) ) + ε d + β d - 1 ( ( { u u ε } J u J u ε ) B s ε ( x 0 ) ) + β B s ε ( x 0 ) | u + - u ε - | d d - 1 + ( u ε , A ε , x 0 ′′ ) ,

so that, using (6.25), (7.14), (7.15), and the fact that s ε ( 1 - 3 θ ) ε we deduce that

lim sup ε 0 ( z ε , A ε , x 0 ) ε d - 1 ( 1 - 3 θ ) d - 1 lim sup ε 0 𝐦 ( u , B ε ( x 0 ) ) ε d - 1 + β ( 1 + | [ u ¯ x 0 surf ] | ) γ d - 1 ( 1 - ( 1 - 4 θ ) d - 1 ) .

Then one concludes exactly in the same way, upon replacing β in equation (6.31) with the larger constant β ( 1 + | [ u ¯ x 0 surf ] | ) . ∎

Proof of Theorem 7.1.

The result follows from Lemmas 7.2, 7.4, and 7.6, arguing exactly as in the proof of Theorem 2.1. ∎

8 The GSBV p case

In this section we briefly remark that the strategy devised in this paper allows also to establish an integral representation result in the space GSBV p ( Ω ; m ) for m . We consider functionals

: GSBV p ( Ω ; m ) × ( Ω ) [ 0 , + )

with the following general assumptions:

  1. ( u , ) is a Borel measure for any u GSBV p ( Ω ; m ) ,

  2. ( , A ) is lower semicontinuous with respect to convergence in measure on Ω for any A 𝒜 ( Ω ) ,

  3. ( , A ) is local for any A 𝒜 ( Ω ) , in the sense that if u , v GSBV p ( Ω ; m ) satisfy u = v a.e. in A, then ( u , A ) = ( v , A ) ,

  4. there exist 0 < α < β such that for any u GSBV p ( Ω ; m ) and B ( Ω ) we have

    α ( B | u | p d x + d - 1 ( J u B ) ) ( u , B ) β ( B ( 1 + | u | p ) d x + d - 1 ( J u B ) ) .

In this setting, we replace definition (2.1) by

𝐦 ( u , A ) = inf v GSBV p ( Ω ; m ) { ( v , A ) : v = u  in a neighborhood of  A } .

Moreover, as in (2.2)–(2.3), we define the functions x 0 , u 0 , ξ ( x ) = u 0 + ξ ( x - x 0 ) and u x 0 , a , b , ν ( x ) = a on the set { ( x - x 0 ) ν > 0 } and u x 0 , a , b , ν ( x ) = b on the set { ( x - x 0 ) ν < 0 } for x 0 Ω , u 0 m , ξ 𝕄 m × d , a , b m , and ν 𝕊 d - 1 .

Theorem 8.1 (Integral representation in GSBV p ).

Let Ω R d be open, bounded with Lipschitz boundary, let m N , and suppose that F : GSBV p ( Ω ; R m ) × B ( Ω ) [ 0 , + ) satisfies ($\mathbb{H}$1)($\mathbb{H}$4). Then

( u , B ) = B f ( x , u ( x ) , u ( x ) ) d x + J u B g ( x , u + ( x ) , u - ( x ) , ν u ( x ) ) d d - 1 ( x )

for all u GSBV p ( Ω ; R m ) and B B ( Ω ) , where f is given by

f ( x 0 , u 0 , ξ ) = lim sup ε 0 𝐦 ( x 0 , u 0 , ξ , B ε ( x 0 ) ) γ d ε d

for all x 0 Ω , u 0 R m , ξ M m × d , and g is given by

g ( x 0 , a , b , ν ) = lim sup ε 0 𝐦 ( u x 0 , a , b , ν , B ε ( x 0 ) ) γ d - 1 ε d - 1

for all x 0 Ω , a , b R m , and ν S d - 1 .

We point out that integral representation results in GSBV p have been used in several contributions, see, e.g., [6, 7, 8, 9, 21, 37]. They all rely on [13] along with a perturbation and truncation argument as follows: first, one considers the regularization

σ ( u ) := ( u ) + σ J u | [ u ] | d d - 1 ,

restricted to u SBV p ( Ω ; m ) . Then the assumptions of the integral representation result in SBV p [13] are satisfied and one obtains a representation of σ . In a second step, this representation is extended to GSBV p by a truncation argument which allows to approximate GSBV p functions by SBV p functions. Eventually, by sending σ 0 , an integral representation result for the original functional can be obtained. We refer to [21, Theorem 4.3, Theorem 5.1] for details on this procedure. (We notice that in [21] a more general growth condition from above is allowed in the surface energy density, cf. [21, assumption (1.4)], analogous to the one in (H4’).) With our result at hand, this method can be considerably simplified since no perturbation and truncation arguments are needed.

For the proof of the theorem we need the following Poincaré-type inequality, which can be directly deduced from Theorem 3.3.

Theorem 8.2 (Poincaré inequality for functions with small jump set).

Assume Ω R d is a bounded Lipschitz domain and let 1 < p < + . Then there exists a constant c = c ( Ω , p , m ) > 0 such that for all u GSBV p ( Ω ; R m ) there is a set of finite perimeter ω Ω with

d - 1 ( * ω ) c d - 1 ( J u ) , d ( ω ) c ( d - 1 ( J u ) ) d d - 1

and v W 1 , p ( Ω ; R m ) such that v = u on Ω ω and

(8.1) v L p ( Ω ) c u L p ( Ω ) .

In particular, for all u GSBV p ( Ω ; R m ) there is a constant b R m such that

u - b L p ( Ω ω ) c u L p ( Ω ) .

Proof.

It suffices to consider the case m = 1 and to prove (8.1). This can be obtained for instance by applying Theorem 3.3 to the function u ¯ : Ω d defined as u ¯ := ( u , 0 , , 0 ) and using the Sobolev–Korn inequality to get v on the left-hand side. ∎

Proof of Theorem 8.1.

We follow the proof of Theorem 2.1 and only indicate briefly the necessary adaptions. First, we observe that a version of the fundamental estimate in Lemma 3.8 holds true in GSBV p ( Ω ; m ) by repeating the proof with ($\mathbb{H}$4) in place of (H4). (We refer also to [16, Proposition 3.1].) Recall that the result follows by combining Lemmas 4.1, 4.2, and 4.3.

Lemma 4.1: The result is proved via (4.15) and Lemma 4.4. Notice indeed that the derivation of (4.15), as well as the argument ensuing therefrom are the same, up to using the growth condition ($\mathbb{H}$4) instead of (H4). (We also refer to [13, Lemma 6] for the corresponding argument in SBV p .) In the proof of Lemma 4.4, due to the (stronger) lower bound in ($\mathbb{H}$4) and Ambrosio’s compactness theorem in GSBV p (see [5, Theorem 4.36]), one can ensure that the function v δ defined in (4.8) now belongs to GSBV p ( Ω ; m ) . Then the result follows with the same argument, up to using Theorem 8.2 in place of Theorem 3.3.

Lemma 4.2: With the fundamental estimate in GSBV p at hand, we can follow the proof of Lemma 4.2 for each u GSBV p ( Ω ; m ) with ($\mathbb{H}$4) instead of (H4). The family ( u ε ) ε is defined as in Lemma 5.1 (see (5.5)) using Theorem 8.2 in place of Theorem 3.3. First, u ε GSBV p ( B ε ( x 0 ) ; m ) since u GSBV p ( Ω ; m ) , and u ε B ( 1 - θ ) ε ( x 0 ) W 1 , p ( B ( 1 - θ ) ε ( x 0 ) ; m ) . Observing that

(8.2) lim ε 0 ε - d B ε ( x 0 ) | u ( x ) - u ( x 0 ) | p d x = 0

for d -a.e. x 0 Ω (as u L p ( Ω ; 𝕄 m × d ) ) and using (8.1), we further get

lim ε 0 ε - d B ( 1 - θ ) ε ( x 0 ) | u ε ( x ) - u ( x 0 ) | p d x = 0 .

With (5.5), and using again (8.2), we deduce that (5.1) (iii) can be improved to

lim ε 0 ε - d B ε ( x 0 ) | u ε ( x ) - u ( x 0 ) | p d x = 0 .

This adaption is enough to redo the proof of Lemma 4.2 in the present situation.

Lemma 4.3: Here, we can follow the proof of Lemma 4.3 for each u GSBV p ( Ω ; m ) with ($\mathbb{H}$4) instead of (H4), and Theorem 8.2 in place of Theorem 3.3. The family ( u ε ) ε defined in Lemma 6.1 needs to satisfy u ε GSBV p ( B ε ( x 0 ) ; m ) and (6.1) (iv) needs to be improved to

(8.3) lim ε 0 ε - ( d - 1 ) B ε ( x 0 ) | u ε | p d x = 0 .

First, we use (6.7) to see that u ε GSBV p ( B ε ( x 0 ) ; m ) . Arguing as in the proof of (6.1) (iv), with Theorem 8.2 at hand, we obtain

lim ε 0 ε - ( d - 1 ) B s ε ± ( x 0 ) | u ε | p d x = 0 .

This along with lim ε 0 ε - ( d - 1 ) B ε ( x 0 ) | u | p d x = 0 for d - 1 -a.e. x 0 J u and (6.7) concludes the proof of equation (8.3). ∎


Communicated by Jan Kristensen


Award Identifier / Grant number: 793018

Award Identifier / Grant number: FR 4083/1-1

Award Identifier / Grant number: EXC 2044-390685587

Award Identifier / Grant number: PRIN 2017

Funding statement: For this project Vito Crismale has received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No. 793018. Manuel Friedrich is supported by the DFG project FR 4083/1-1 and by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy EXC 2044-390685587, Mathematics Münster: Dynamics–Geometry–Structure. The work of Francesco Solombrino is part of the project “Variational methods for stationary and evolution problems with singularities and interfaces” PRIN 2017 financed by the Italian Ministry of Education, University, and Research.

References

[1] G. Alberti, Integral representation of local functionals, Ann. Mat. Pura Appl. (4) 165 (1993), 49–86. 10.1007/BF01765841Search in Google Scholar

[2] L. Ambrosio and A. Braides, Functionals defined on partitions in sets of finite perimeter. I. Integral representation and Γ-convergence, J. Math. Pures Appl. (9) 69 (1990), no. 3, 285–305. Search in Google Scholar

[3] L. Ambrosio and A. Braides, Functionals defined on partitions in sets of finite perimeter. II. Semicontinuity, relaxation and homogenization, J. Math. Pures Appl. (9) 69 (1990), no. 3, 307–333. Search in Google Scholar

[4] L. Ambrosio, A. Coscia and G. Dal Maso, Fine properties of functions with bounded deformation, Arch. Ration. Mech. Anal. 139 (1997), no. 3, 201–238. 10.1007/s002050050051Search in Google Scholar

[5] L. Ambrosio, N. Fusco and D. Pallara, Functions of Bounded Variation and Free Discontinuity Problems, Oxford Math. Monogr., Clarendon Press, New York, 2000. Search in Google Scholar

[6] A. Bach, A. Braides and C. I. Zeppieri, Quantitative analysis of finite-difference approximations of free-discontinuity problems, Interfaces Free Bound. 22 (2020), no. 3, 317–381. 10.4171/IFB/443Search in Google Scholar

[7] A. Bach, M. Cicalese and M. Ruf, Random finite-difference discretizations of the Ambrosio–Tortorelli functional with optimal mesh size, preprint (2019), https://arxiv.org/abs/1902.08437. 10.1137/20M1312927Search in Google Scholar

[8] M. Barchiesi and M. Focardi, Homogenization of the Neumann problem in perforated domains: An alternative approach, Calc. Var. Partial Differential Equations 42 (2011), no. 1–2, 257–288. 10.1007/s00526-010-0387-2Search in Google Scholar

[9] M. Barchiesi, G. Lazzaroni and C. I. Zeppieri, A bridging mechanism in the homogenization of brittle composites with soft inclusions, SIAM J. Math. Anal. 48 (2016), no. 2, 1178–1209. 10.1137/15M1007343Search in Google Scholar

[10] A. C. Barroso, I. Fonseca and R. Toader, A relaxation theorem in the space of functions of bounded deformation, Ann. Sc. Norm. Super. Pisa Cl. Sci. (4) 29 (2000), no. 1, 19–49. Search in Google Scholar

[11] G. Bellettini, A. Coscia and G. Dal Maso, Compactness and lower semicontinuity properties in SBD ( Ω ) , Math. Z. 228 (1998), no. 2, 337–351. 10.1007/PL00004617Search in Google Scholar

[12] G. Bouchitté and G. Dal Maso, Integral representation and relaxation of convex local functionals on BV ( Ω ) , Ann. Sc. Norm. Super. Pisa Cl. Sci. (4) 20 (1993), no. 4, 483–533. Search in Google Scholar

[13] G. Bouchitté, I. Fonseca, G. Leoni and L. Mascarenhas, A global method for relaxation in W 1 , p and in SBV p , Arch. Ration. Mech. Anal. 165 (2002), no. 3, 187–242. 10.1007/s00205-002-0220-ySearch in Google Scholar

[14] G. Bouchitté, I. Fonseca and L. Mascarenhas, A global method for relaxation, Arch. Ration. Mech. Anal. 145 (1998), no. 1, 51–98. 10.1007/s002050050124Search in Google Scholar

[15] A. Braides and V. Chiadò Piat, Integral representation results for functionals defined on SBV ( Ω ; m ) , J. Math. Pures Appl. (9) 75 (1996), no. 6, 595–626. Search in Google Scholar

[16] A. Braides, A. Defranceschi and E. Vitali, Homogenization of free discontinuity problems, Arch. Ration. Mech. Anal. 135 (1996), no. 4, 297–356. 10.1007/BF02198476Search in Google Scholar

[17] G. Buttazzo and G. Dal Maso, Γ-limits of integral functionals, J. Anal. Math. 37 (1980), 145–185. 10.1007/BF02797684Search in Google Scholar

[18] G. Buttazzo and G. Dal Maso, A characterization of nonlinear functionals on Sobolev spaces which admit an integral representation with a Carathéodory integrand, J. Math. Pures Appl. (9) 64 (1985), no. 4, 337–361. Search in Google Scholar

[19] G. Buttazzo and G. Dal Maso, Integral representation and relaxation of local functionals, Nonlinear Anal. 9 (1985), no. 6, 515–532. 10.1016/0362-546X(85)90038-0Search in Google Scholar

[20] F. Cagnetti, A. Chambolle and L. Scardia, Korn and Poincaré–Korn inequalities for functions with a small jump set, preprint (2020). 10.1007/s00208-021-02210-wSearch in Google Scholar

[21] F. Cagnetti, G. Dal Maso, L. Scardia and C. I. Zeppieri, Γ-convergence of free-discontinuity problems, Ann. Inst. H. Poincaré Anal. Non Linéaire 36 (2019), no. 4, 1035–1079. 10.1016/j.anihpc.2018.11.003Search in Google Scholar

[22] L. Carbone and C. Sbordone, Some properties of Γ-limits of integral functionals, Ann. Mat. Pura Appl. (4) 122 (1979), 1–60. 10.1007/BF02411687Search in Google Scholar

[23] M. Caroccia, M. Focardi and N. Van Goethem, On the integral representation of variational functionals on BD , SIAM J. Math. Anal. 52 (2020), no. 4, 4022–4067. 10.1137/19M1277564Search in Google Scholar

[24] A. Chambolle, An approximation result for special functions with bounded deformation, J. Math. Pures Appl. (9) 83 (2004), no. 7, 929–954. 10.1016/j.matpur.2004.02.004Search in Google Scholar

[25] A. Chambolle and V. Crismale, Compactness and lower semicontinuity in GSBD , preprint (2018), https://arxiv.org/abs/1802.03302; to appear in J. Eur. Math. Soc. (JEMS). Search in Google Scholar

[26] A. Chambolle and V. Crismale, A density result in GSBD p with applications to the approximation of brittle fracture energies, Arch. Ration. Mech. Anal. 232 (2019), no. 3, 1329–1378. 10.1007/s00205-018-01344-7Search in Google Scholar

[27] A. Chambolle and V. Crismale, Phase-field approximation for a class of cohesive fracture energies with an activation threshold, Adv. Calc. Var. (2020), 10.1515/acv-2019-0018. 10.1515/acv-2019-0018Search in Google Scholar

[28] S. Conti, M. Focardi and F. Iurlano, Integral representation for functionals defined on SBD p in dimension two, Arch. Ration. Mech. Anal. 223 (2017), no. 3, 1337–1374. 10.1007/s00205-016-1059-ySearch in Google Scholar

[29] S. Conti, M. Focardi and F. Iurlano, Which special functions of bounded deformation have bounded variation?, Proc. Roy. Soc. Edinburgh Sect. A 148 (2018), no. 1, 33–50. 10.1017/S030821051700004XSearch in Google Scholar

[30] G. Dal Maso, Integral representation on BV ( Ω ) of Γ-limits of variational integrals, Manuscripta Math. 30 (1979/80), no. 4, 387–416. 10.1007/BF01301259Search in Google Scholar

[31] G. Dal Maso, Generalised functions of bounded deformation, J. Eur. Math. Soc. (JEMS) 15 (2013), no. 5, 1943–1997. 10.4171/JEMS/410Search in Google Scholar

[32] G. Dal Maso and L. Modica, A general theory of variational functionals, Topics in Functional Analysis 1980–81, Quaderni, Scuola Normale Superiore, Pisa (1981), 149–221. Search in Google Scholar

[33] E. De Giorgi, Sulla convergenza di alcune successioni d’integrali del tipo dell’area, Rend. Mat. (6) 8 (1975), 277–294. Search in Google Scholar

[34] E. De Giorgi, M. Carriero and A. Leaci, Existence theorem for a minimum problem with free discontinuity set, Arch. Ration. Mech. Anal. 108 (1989), no. 3, 195–218. 10.1007/BF01052971Search in Google Scholar

[35] G. De Philippis and F. Rindler, On the structure of 𝒜 -free measures and applications, Ann. of Math. (2) 184 (2016), no. 3, 1017–1039. 10.4007/annals.2016.184.3.10Search in Google Scholar

[36] F. Ebobisse and R. Toader, A note on the integral representation of functionals in the space SBD Ω, Rend. Mat. Appl. (7) 23 (2003), no. 2, 189–201. Search in Google Scholar

[37] M. Focardi, M. S. Gelli and M. Ponsiglione, Fracture mechanics in perforated domains: A variational model for brittle porous media, Math. Models Methods Appl. Sci. 19 (2009), no. 11, 2065–2100. 10.1142/S0218202509004042Search in Google Scholar

[38] G. A. Francfort and J.-J. Marigo, Revisiting brittle fracture as an energy minimization problem, J. Mech. Phys. Solids 46 (1998), no. 8, 1319–1342. 10.1016/S0022-5096(98)00034-9Search in Google Scholar

[39] M. Friedrich, A Korn-type inequality in SBD for functions with small jump sets, Math. Models Methods Appl. Sci. 27 (2017), no. 13, 2461–2484. 10.1142/S021820251750049XSearch in Google Scholar

[40] M. Friedrich, A piecewise Korn inequality in SBD and applications to embedding and density results, SIAM J. Math. Anal. 50 (2018), no. 4, 3842–3918. 10.1137/17M1129982Search in Google Scholar

[41] M. Friedrich, A compactness result in GSBV p and applications to Γ-convergence for free discontinuity problems, Calc. Var. Partial Differential Equations 58 (2019), no. 3, Paper No. 86. Search in Google Scholar

[42] M. Friedrich, M. Perugini and F. Solombrino, Γ-convergence for free-discontinuity problems in linear elasticity: Homogenization and relaxation, preprint (2020), https://arxiv.org/abs/2010.05461. Search in Google Scholar

[43] M. Friedrich, M. Perugini and F. Solombrino, Lower semicontinuity for functionals defined on piecewise rigid functions and on GSBD , preprint (2020), https://arxiv.org/abs/2002.08133. Search in Google Scholar

[44] M. Friedrich and F. Solombrino, Functionals defined on piecewise rigid functions: Integral representation and Γ-convergence, Arch. Ration. Mech. Anal. 236 (2020), no. 3, 1325–1387. 10.1007/s00205-020-01493-8Search in Google Scholar

[45] A. A. Griffith, The phenomena of rupture and flow in solids, Philos. Trans. Roy. Soc. London Ser. A 221 (1920), 163–198. Search in Google Scholar

[46] J. W. Hutchinson, A course on nonlinear fracture mechanics, Department of Solid Mechanics, Technical University of Denmark, Lyngby, 1989. Search in Google Scholar

[47] C. B. Morrey, Jr., Quasi-convexity and the lower semicontinuity of multiple integrals, Pacific J. Math. 2 (1952), 25–53. 10.2140/pjm.1952.2.25Search in Google Scholar

[48] F. Rindler, Lower semicontinuity for integral functionals in the space of functions of bounded deformation via rigidity and Young measures, Arch. Ration. Mech. Anal. 202 (2011), no. 1, 63–113. 10.1007/s00205-011-0408-0Search in Google Scholar

[49] C. Sbordone, Su alcune applicazioni di un tipo di convergenza variazionale, Ann. Sc. Norm. Super. Pisa Cl. Sci. (4) 2 (1975), no. 4, 617–638. Search in Google Scholar

[50] R. Temam, Mathematical Problems in Plasticity, Gauthier-Villars, Paris, 1983. Search in Google Scholar

Received: 2020-05-14
Revised: 2020-10-13
Accepted: 2020-10-20
Published Online: 2020-11-07
Published in Print: 2022-10-01

© 2020 Walter de Gruyter GmbH, Berlin/Boston

This work is licensed under the Creative Commons Attribution 4.0 International License.

Downloaded on 25.4.2024 from https://www.degruyter.com/document/doi/10.1515/acv-2020-0047/html
Scroll to top button