Processing math: 100%
Research article Special Issues

Homogenisation of high-contrast brittle materials

  • Received: 13 October 2019 Accepted: 26 November 2019 Published: 31 December 2019
  • This paper is an overview on some recent results concerning the variational analysis of static fracture in the so-called high-contrast brittle composite materials. The paper is divided into two main parts. The first part is devoted to establish a compactness result for a general class of free-discontinuity functionals with degenerate (or high-contrast) integrands. The second part is focussed on some specific examples which show that the degeneracy of the integrands may lead to non-standard limit effects, which are specific to this high-contrast setting.

    Citation: Caterina Ida Zeppieri. Homogenisation of high-contrast brittle materials[J]. Mathematics in Engineering, 2020, 2(1): 174-202. doi: 10.3934/mine.2020009

    Related Papers:

    [1] Manuel Friedrich . Griffith energies as small strain limit of nonlinear models for nonsimple brittle materials. Mathematics in Engineering, 2020, 2(1): 75-100. doi: 10.3934/mine.2020005
    [2] Stefano Almi, Giuliano Lazzaroni, Ilaria Lucardesi . Crack growth by vanishing viscosity in planar elasticity. Mathematics in Engineering, 2020, 2(1): 141-173. doi: 10.3934/mine.2020008
    [3] Giacomo Canevari, Arghir Zarnescu . Polydispersity and surface energy strength in nematic colloids. Mathematics in Engineering, 2020, 2(2): 290-312. doi: 10.3934/mine.2020015
    [4] Zhiwen Zhao . Exact solutions for the insulated and perfect conductivity problems with concentric balls. Mathematics in Engineering, 2023, 5(3): 1-11. doi: 10.3934/mine.2023060
    [5] Paola F. Antonietti, Chiara Facciolà, Marco Verani . Unified analysis of discontinuous Galerkin approximations of flows in fractured porous media on polygonal and polyhedral grids. Mathematics in Engineering, 2020, 2(2): 340-385. doi: 10.3934/mine.2020017
    [6] Raimondo Penta, Ariel Ramírez-Torres, José Merodio, Reinaldo Rodríguez-Ramos . Effective governing equations for heterogenous porous media subject to inhomogeneous body forces. Mathematics in Engineering, 2021, 3(4): 1-17. doi: 10.3934/mine.2021033
    [7] Patrick Dondl, Martin Jesenko, Martin Kružík, Jan Valdman . Linearization and computation for large-strain visco-elasticity. Mathematics in Engineering, 2023, 5(2): 1-15. doi: 10.3934/mine.2023030
    [8] Federico Cluni, Vittorio Gusella, Dimitri Mugnai, Edoardo Proietti Lippi, Patrizia Pucci . A mixed operator approach to peridynamics. Mathematics in Engineering, 2023, 5(5): 1-22. doi: 10.3934/mine.2023082
    [9] Chiara Gavioli, Pavel Krejčí . Deformable porous media with degenerate hysteresis in gravity field. Mathematics in Engineering, 2025, 7(1): 35-60. doi: 10.3934/mine.2025003
    [10] Luca Formaggia, Alessio Fumagalli, Anna Scotti . A multi-layer reactive transport model for fractured porous media. Mathematics in Engineering, 2022, 4(1): 1-32. doi: 10.3934/mine.2022008
  • This paper is an overview on some recent results concerning the variational analysis of static fracture in the so-called high-contrast brittle composite materials. The paper is divided into two main parts. The first part is devoted to establish a compactness result for a general class of free-discontinuity functionals with degenerate (or high-contrast) integrands. The second part is focussed on some specific examples which show that the degeneracy of the integrands may lead to non-standard limit effects, which are specific to this high-contrast setting.


    In this note we analyse the large-scale behaviour of high-contrast composite materials which can undergo fracture. In a variational setting, the microscopic behaviour of high-contrast composites is typically described by means of scale-dependent energy functionals with "degenerate" integrands. For brittle materials the scale-dependent energies are of the general form

    Fε(u)=Ωfε(x,u)dx+Sugε(x,νu)dHn1, (1.1)

    where ε>0 describes both the composite-microstructure and the degeneracy of the mechanical properties of the material (cf. (1.3)). In (1.1) the variable u:ΩRnR belongs to SBV(Ω), the space of special functions of bounded variation in Ω. In this simplified scalar setting, u represents an anti-plane displacement and Ω is the cross-section of an infinite cylindrical body. Being u an SBV-function, discontinuities are allowed and the discontinuity set of u, denoted by Su, models the cracks in the material. The deformation gradient Du can be decomposed into the sum of a bulk part udx and a surface part (u+u)νuHn1Su, where u denotes the approximate gradient of u, u+ and u the traces of u on both sides of Su, and νu denotes the (generalised) normal to Su. The volume term in Fε represents the elastic energy stored in the unfractured part of the material, whereas the surface term in Fε accounts for the presence of cracks. According to the Griffith criterion, in brittle materials, already for the smallest crack-amplitude, there is no interaction between the two lips of the crack, so that the corresponding fracture energy does not depend on [u]=u+u.

    For finite-contrast brittle materials, the limit behaviour of energies of type (1.1) is by-now well-understood and the corresponding theory provides a rigorous micro-to-macro upscaling for brittle fracture. In fact, if fε and gε satisfy (mild regularity assumptions and) standard growth and coercivity conditions of type

    c1|ξ|pfε(x,ξ)c2(1+|ξ|p)andc3gε(x,ν)c4, (1.2)

    for every ε>0, x,ξRn, νSn1, for some p>1, and 0<c1c2<+, 0<c3c4<+, then in [18] Giacomini and Ponsiglione showed, among other, that the limit behaviour of Fε is captured by a scale-independent free-discontinuity functional of the same type as Fε; i.e.,

    F0(u)=Ωf0(x,u)dx+Sug0(x,νu)dHn1,

    with f0 and g0 also satisfying (1.2). Under these assumptions, Giacomini and Ponsiglione also showed that volume and surface energy decouple in the limit, so that the energy density f0 is not affected by the presence of the surface term in Fε, whereas the surface energy density g0 is not affected by the volume term in Fε. In a recent work, Cagnetti, Dal Maso, Scardia and Zeppieri [12] generalised the asymptotic analysis carried out in [18] and devised (nearly optimal) sufficient conditions which ensure a macroscopic bulk-surface energy decoupling for a wide class of finite-contrast vectorial free-discontinuity functionals which may also depend on [u]. The class of periodic free-discontinuity functionals originally analysed by Braides, Defranceschi and Vitali [9] satisfy the sufficient conditions provided in [12]. Moreover, random free-discontinuity functionals with stationary finite-contrast integrands can be also seen as a special instance of those treated in [12], as shown by Cagnetti, Dal Maso, Scardia and Zeppieri in [13]. Therefore, a volume-surface interaction can be ruled out for a large class of finite-contrast free-discontinuity functionals. In this setting, in particular, microscopic brittle energies always converge to macroscopic brittle energies. However, the general theory established in [9,12,13,18] is not well-suited for studying the large-scale behaviour of those brittle composites whose different constituents have very different mechanical properties from one another. Indeed, in this case the integrands fε and gε in (1.1) may exhibit a so-called high-contrast behaviour and satisfy (1.2) only in a subset Ωε of Ω.

    In the last decade there has been an ever increasing interest in the study of high-contrast free-discontinuity functionals and in the derivation of their effective properties. In particular, the case where (at least) one of the conditions in (1.2) is violated in "many small" periodically distributed regions inside Ω has been considered (see, e.g., [3,4,5,6,11,14,16,17,19,20,21]). Depending on the type of degeneracy of fε and gε, nonstandard limit effects have been also observed. These nonstandard effects are typical of the high-contrast setting and arise from a nontrivial volume-surface limit interaction, which cannot be excluded in this degenerate setting. In fact, in the two companion papers [4,16], Barchiesi, Dal Maso and Zeppieri show that when only gε is degenerate, already for very simple free-discontinuity functionals of Mumford-Shah type, a bulk-surface interaction cannot be ruled out. Namely, a volume-surface coupling can be observed when homogenising a material made of "many" purely brittle inclusions periodically distributed in a connected unbreakable structure, whose fracture-resistance is assumed to be infinite. This coupling produces a homogeneous material whose overall behaviour is of ductile (or cohesive) type; in other words, the homogenised surface energy explicitly depends on [u]. A similar phenomenon is also observed by Barchiesi, Lazzaroni and Zeppieri [6] who show that a ductile behaviour can be seen as the macroscopic effect of a nontrivial volume-surface interaction in the homogenisation of two purely brittle materials with a high-contrast bulk energy. Moreover, in the recent work [19] Pellet, Scardia and Zeppieri prove, instead, that nonstandard constitutive laws may arise when homogenising two purely brittle materials with a high-contrast surface energy. The functionals analysed in [6] and [19] are both of type

    Fαε,βεε(u)=Ωaε(xε)|u|2dx+Subε(xε)dHn1 (1.3)

    where the elastic modulus aε and the fracture resistance (or fracture thoughness) bε are Q-periodic functions and in the unit periodicity cell Q:=(1/2,1/2)n are defined as

    aε(y)={αεify¯Qr1ifyQ¯Qrbε(y)={βεify¯Qr1ifyQ¯Qr

    with αε,βε[0,1], r(0,1), and Qr:=(r/2,r/2)n. Since αε,βε are not bounded away from zero, the functions aε and bε can be degenerate. In their turn, the integrands fε(y,ξ)=aε(y)|ξ|2 and gε(y,ν)=bε(y) in (1.3) will not satisfy, in general, the coercivity conditions in (1.2).

    The limit case αε=βε=0 corresponds to the case of periodically perforated brittle materials studied by Cagnetti and Scardia [14] and by Focardi, Gelli, and Ponsiglione [17] (see also Barchiesi and Focardi [5] for more general free-discontinuity functionals). In spite of the strong degeneracy of the coefficients aε and bε, which in this case are equal to zero in a "large" portion of Ω, in this case it can be proven that the functionals F0,0ε exhibit a limit behaviour which is qualitatively similar to that of free-discontinuity functionals with coercive integrands. Namely, in this case bulk and surface terms do not interact in the limit.

    The aim of this note is to show that, contrary to the coercive case, where general homogenisation results can be proven to describe the limit behaviour of a large class of free-discontinuity functionals, in the non-coercive setting, already for special functionals of type (1.3), a unified homogenisation theory cannot be established. In fact, the limit behaviour of Fαε,βεε is highly sensitive both to the choice of the parameters αε and βε and to their vanishing rate compared to the period of the microstructure ε.

    This note is divided into two main parts and organised as follows. In first part we will deal with sequences of general free-discontinuity functionals of type (1.1) whose coefficients fε and gε are "weakly coercive" or "degenerate"; i.e., they satisfy the lower bounds in (1.2) only in a set Ωε which is obtained removing from Ω many small periodically distributed connected regions. We will use the localisation method of Γ-convergence [8,15] to prove that these kind of functionals are (pre)compact. That is, up to subsequences, they always Γ-converge to a free-discontinuity functional of type

    Ωf(x,u)dx+Sug(x,[u],νu)dHn1.

    Moreover, the limit integrands f and g are non-degenerate and satisfy coercivity conditions of type (1.2) for some positive constants ˆc1,ˆc3 which are strictly smaller than c1,c3, respectively. In this part of the analysis a pivotal role is played by an extension result for SBV-functions defined in periodically perforated domains, proved by Cagnetti and Scardia [14] (see also the later variant in [5]).

    In the second part of this note we will specialise the general theory to some prototypical and yet relevant model cases. Namely, we will briefly review the case of perforated (or porous) brittle materials studied by Cagnetti and Scardia [14] and by Focardi, Gelli and Ponsiglione [17] (see also [5]), the case of high-contrast brittle materials with soft inclusions treated by Barchiesi, Lazzaroni and Zeppieri in [6], and eventually the case of high-contrast brittle materials with weak inclusions analysed by Pellet, Scardia and Zeppieri in [19]. In particular we will show that the choice of the integrands fε and gε in (1.1) strongly affects the form of the Γ-limit which can give rise to macroscopic models accounting for damage as well as to models accounting for cohesive fracture.

    In this part we will use the localisation method of Γ-convergence [8,15] to prove a convergence result for a general class of free-discontinuity functionals of brittle type, with degenerate coefficients.

    In the choice of the convergence to compute the Γ-limit, a crucial role will be played by an extension result for SBV-functions defined in periodically perforated domains due to Cagnetti and Scardia [14,Theorem 1.3] and by a later variant due to Barchiesi and Focardi [5,Theorem 1].

    We list below a few notation which will be used throughout the paper.

    ΩRn denotes an open and bounded set with Lipschitz boundary. The set A(Ω) denotes the collection of all open subsets of Ω;

    Q denotes the open unit cube of Rn centred at the the origin, whereas for xRn and r>0 we set Qr(x):=rQ+x;

    ● for νSn1 we denote with Qν the open unit cube of Rn centred at the the origin, with one face orthogonal to ν and for xRn and r>0 we set Qνr(x):=rQν+x;

    ● for xRn and νSn1 we denote by Πν(x) the hyperplane through x and perpendicular to ν; i.e., Πν(x):={yRn:(yx)ν=0}. If x=0 we simply write Πν;

    ● For uL1(Ω) and m>0 the function um denotes the truncated function of u at level m; i.e., um:=(um)(m);

    ● For ξRn we denote by uξ the linear function with gradient equal to ξ; i.e., uξ(x):=ξx, for every xRn;

    ● For xRn, tR, and νSn1 we denote with uν,tx the piecewise constant function taking values 0,t and jumping across the hyperplane Πν(x); i.e.,

    uν,tx(y):={tif(yx)ν0,0if(yx)ν<0.

    The functional setting we are going to consider in this note is that of SBV, the space of special functions of bounded variation. We recall here only the definition of the spaces which are relevant for our analysis and we refer the reader to [2] for a comprehensive treatment on the subject. We set

    SBV(Ω):={uBV(Ω): Du=uLn+(u+u)νudHn1Su}.

    Here Su denotes the approximate discontinuity set of u, νu is the generalised normal to Su, u+ and u are the traces of u on both sides of Su. In this paper we work with the following vector subspace of SBV(Ω)

    SBVp(Ω):={uSBV(Ω): uLp(Ω) and Hn1(Su)<+},

    where p>1. We consider also the larger space of generalised special functions of bounded variation in Ω,

    GSBV(Ω):={uL1(Ω): umSBV(Ω) for all mN},

    as well as

    GSBVp(Ω):={uGSBV(Ω): uLp(Ω) and Hn1(Su)<+}.

    We consider also

    SBVpc(Ω):={uSBV(Ω):u=0Lna.e.Hn1(Su)<+};

    it is known (see [2,Theorem 4.23]) that every u in SBVpc(Ω)L(Ω) is piecewise constant in the sense of [2,Definition 4.21], namely there exists a Caccioppoli partition (Ei) of Ω such that u is constant Ln-a.e. in each set Ei. Moreover, we set

    P(Ω):={uSBVpc(Ω):u(x){0,1} Lna.e.inΩ}.

    For u,wL1(U), in what follows, by "u=w near U" we mean that there exists a neighbourhood V of U in Rn such that u=w Ln-a.e. in VU. Let fk:Rn×Rn[0,+) be Carathéodory functions such that

    (H1) there exist p>1 and 0<c1c2<+ such that for every (x,ξ)Rn×Rn and for every kN

    c1|ξ|pfk(x,ξ)c2(1+|ξ|p); (2.1)

    (H2) fk(x,0)=0 for every xRn and for every kN.

    Let moreover gk:Rn×Sn1(0,+) be Borel functions such that

    (H3) there exist 0<c3c4<+ such that for every (x,ν)Rn×Sn1 and every kN

    c3gk(x,ν)c4; (2.2)

    (H4) gk(x,ν)=gk(x,ν), for every (x,ν)Rn×Sn1 and every kN.

    Let ΩRn be open bounded and with Lipschitz boundary and let KQ be compact and such that QK has a Lipschitz boundary. We define

    E:=RniZn(K+i);

    the set E is open, connected, Q-periodic, and has a Lipschitz boundary. Let εk be a sequence of positive numbers such that εk0 as k+ and denote by Ωk the εkQ-periodic set defined as Ωk:=ΩεkE (see Figure 1).

    Figure 1.  Schematic of a high-contrast composite material.

    Let moreover αk,βk[0,1] and consider the sequence of functionals Fk:L1(Ω)[0,+] defined as

    Fk(u):={Ωkfk(x,u)dx+αkΩΩkfk(x,u)dx+SuΩkgk(x,νu)dHn1+βkSu(ΩΩk)gk(x,νu)dHn1ifuSBVp(Ω),+otherwise in L1(Ω). (2.3)

    We observe that thanks to assumption (H2) the functionals Fk decrease by truncation, wherease they do not satisfy the standard coercivity conditions required, e.g., in [9,12,18] since the coefficients αk,βk are not bounded away from zero.

    Due to the possible degeneracy of the coefficients αk and βk, the functionals Fk are not, in general, equi-coercive with respect to the strong L1(Ω)-convergence. Similarly as in [5,11,17,19], in what follows we give a notion of convergence on L1(Ω) which is weaker that the L1(Ω)-convergence and ensures the equi-coercivity of the functionals Fk. This will be done by appealing to [5,Theorem 1]. For the readers' convenience we recall here a slightly simplified version of this result which is useful for our purposes.

    Theorem 2.1 (cf. Theorem 1 in [5]). Let (uk)SBVp(Ωk) be such that

    supk(Ωk|uk|pdx+Ωk|uk|pdx+Hn1(SukΩk))<+. (2.4)

    Then, there exist (˜uk)SBVp(Ω), with ˜uk=uk a.e. in Ωk, and a function uGSBVp(Ω)Lp(Ω) such that (up to subsequences) ˜uku in L1(Ω).

    If moreover supkukL(Ωk)<+ then uSBVp(Ω)L(Ω) and ˜uku in Lp(Ω).

    Let (uk)L1(Ω) be a sequence satisfying

    supkukLp(Ωk)<+andsupkFk(uk)<+.

    Then, clearly (uk)SBVp(Ωk); moreover in view of (H1) and (H3) the sequence (uk) satisfies the uniform bound (2.4). Therefore invoking Theorem 2.1 immediately yields the existence of a function uGSBVp(Ω)Lp(Ω) and a sequence (˜uk)SBVp(Ω) with ˜uk=uk a.e. in Ωk, such that (up to subsequences not relabelled) ˜uku in L1(Ω).

    This observation motivates the choice of the following notion of convergence on L1(Ω).

    Definition 2.2 (Convergence). Let (uk) be a sequence in L1(Ω). We say that (uk) converges to a function uL1(Ω), and we write uku, if there exists a sequence (˜uk)L1(Ω) such that ˜uk=uk a.e. in Ωk, and ˜uk converges to u in L1(Ω).

    Remark 2.3 (Uniqueness of the limit). We observe that since C(K):=Ln(QK)>0, then the limit in the sense of Definition 2.2 is well-defined. Indeed, assume that uku1 and uku2. Then by definition there exist (˜u1,k),(˜u2,k)L1(Ω) such that ˜u1,k=˜u2,k=uk in Ωk and ˜u1,ku1 and ˜u2,ku2 in L1(Ω). Therefore

    0=limk+Ωk|˜u1,k˜u2,k|dx=limk+Ω|˜u1,k˜u2,k|χΩkdx=C(K)Ω|u1u2|dx,

    where the last inequality follows by the Riemann-Lebesgue Theorem applied to the εkQ-periodic function χΩk. Then, since C(K)>0 we necessarily have u1=u2 a.e. in Ω.

    We notice moreover that the convergence uku readily implies

    limk+ukuL1(Ωk)=0.

    Remark 2.4 (Convergence of truncated functions). Let (uk)L1(Ω) be such that uku for some uL1(Ω). Let mN and denote by (umk) the sequence of truncated functions of uk at level m, then umkum where um denotes the truncated function of u at level m. Indeed, set vk:=(˜uk)m, then vk=umk a.e. in Ωk, moreover since ˜uku in L1(Ω) then (˜uk)mum in L1(Ω), and actually in any Lp(Ω).

    In what follows we study the Γ-convergence of the functionals Fk with respect to the convergence as in Definition 2.2. To this end we give the following sequential notion of Γ-convergence.

    Definition 2.5 (Sequential Γ-convergence). Let Fk,F:L1(Ω)[0,+]; we say that the functionals Fε -converge to F with respect to the convergence as in Definition 2.2 if for every uL1(Ω) the two following conditions are satisfied:

    (i) (Ansatz-free lower bound) For every (uk)L1(Ω) with uku we have

    F(u)lim infk+Fk(uk);

    (ii) (Existence of a recovery sequence) There exists (ˉuk)L1(Ω) with ˉuku such that

    F(u)lim supk+Fk(ˉuk).

    Remark 2.6. It is standard to show that F is lower semicontinuous with respect to the convergence as in Definition 2.2 and hence with respect to the strong L1(Ω)-convergence.

    For every uL1(Ω) we consider the functionals

    Γ-lim infk+Fk(u):=inf{lim infk+Fk(uk):uku} (2.5)

    and

    Γ-lim supk+Fk(u):=inf{lim supk+Fk(uk):uku}. (2.6)

    It is easy to show that the infima in (2.5) and (2.6) are actually attained.

    In what follows we also use the compact notation

    F(u):=Γ-lim infk+Fk(u) and F"(u):=Γ-lim supk+Fk(u). (2.7)

    It is immediate to see that Definition 2.5 is equivalent to F=F"=F in L1(Ω).

    Remark 2.7 (The case αk,βk=0). In the case of porous brittle materials [5,14,17], which corresponds to the parameter choice αk,βk=0, the Γ-convergence of the functionals Fk can be equivalently studied with respect to the strong L1(Ω)-convergence. Indeed, in this case a sequence (uk) with equibounded energy can be replaced by the L1(Ω)-converging sequence (˜uk) given by Theorem 2.1, without changing the energy.

    The following proposition shows that the domain of the Γ-limit of Fk (if it exists) is GSBVp(Ω).

    Proposition 2.8 (Domain of the Γ-limit). Let F and F" be as in (2.7); then

    domF=domF"=GSBVp(Ω).

    Proof. We first show that GSBVp(Ω)domF". By the growth conditions (2.1) and (2.2) we have Fk(u)G(u) where

    G(u):={c2Ω(1+|u|p)dx+c3Hn1(SuΩ)in GSBVp(Ω)+otherwise in L1(Ω) (2.8)

    The functional G is lower semicontinuous with respect to the strong L1(Ω)-convergence, hence we have

    inf{lim supk+Fk(uk):ukuinL1(Ω)}G(u).

    Since the convergence in Definition 2.2 is weaker than the L1(Ω)-convergence we then have F"≤G, and thus the desired inclusion.

    We now prove that domFGSBVp(Ω). To this end, let udomF then there exists (uk)L1(Ω) with uku such that lim infkFk(uk)=F(u)<+. Then, up to subsequences (not relabelled) we have supkF(uk)<+, thus in particular (uk)SBVp(Ω).

    Let mN and let umk be the truncated function of uk at level m; then (umk)SBVp(Ω)L(Ω). Since the functionals Fk decrease by truncation, for every fixed mN it also holds supkFk(umk)<+. Therefore, for mN fixed we can appeal to Theorem 2.1 to deduce the existence of a sequence (vk)L1(Ω) such that vk=umk a.e. in Ωk and of a function vSBVp(Ω) such that up to subsequence (not relabelled) vkv in L1(Ω). Since uku we have

    0=limk+Ωk|vkumk|dx=limk+Ω|vk(˜uk)m|χΩkdx=C(K)Ω|vum|dx,

    therefore v=um a.e. in Ω. Eventually, the arbitrariness of mN yields uGSBVp(Ω).

    In this section we show that, up to subsequences, the functionals Fk Γ-converge to a free-discontinuity functional of the form

    F(u)=Ωf(x,u)dx+Sug(x,u+u,νu)dHn1.

    for some f and g. Moreover, we show that, despite the degeneracy of the coefficients αk,βk, the limit integrands f and g satisfy standard coercivity conditions similar to (2.1) and (2.2), respectively.

    If not otherwise specified, in what follows the Γ-convergence of the functionals Fk is always understood in the sense of Definition 2.5.

    To prove the existence of a Γ-convergent subsequence of Fk we make use of the so-called localisation method [8,15] which we adapt to the sequential notion of Γ-convergence as in Definition 2.5.

    We start by localising the functionals Fk; that is we consider Fk:L1(Ω)×A(Ω)[0,+] defined as

    Fk(u,U):={Ukfk(x,u)dx+αkUUkfk(x,u)dx+SuUkgk(x,νu)dHn1+βkSu(UUk)gk(x,νu)dHn1ifuSBVp(U),+otherwise in L1(Ω), (2.9)

    where Uk:=UεkE.

    We also define the localised versions of (2.5) and (2.6); i.e., for every UA(Ω) we consider the functionals defined as

    F(,U):=Γ-lim infk+Fk(,U),F"(,U):=Γ-lim supk+Fk(,U). (2.10)

    Remark 2.9 (Properties of F,F"). It is easy to show that F and F" are lower semicontinuous with respect to the convergence in Definition 2.2, local, and that they decrease by truncation. Moreover, as set functions they are both increasing, whereas F is also superadditive.

    Remark 2.10 (On assumption (H2)). If we drop assumption (H2) the functionals Fk will not decrease by truncation, but rather satisfy

    Fk(um,U)Fk(u,U)+c2Ln(Uk{|u|m})+αkc2Ln(UUk{|u|m}). (2.11)

    If αk is infinitesimal, the inequality in (2.11) implies

    F(um,U)F(u,U)+c2muL1(Ω), (2.12)

    (and analogously for F"). In fact, by definition of Γ-liminf there exists a sequence (uk)L1(Ω) such that uku and F(u,U)=lim infk+Fk(uk,U). Then if umk is the truncated function of uk at level m, by (2.11) we get

    Fk(umk,U)Fk(uk,U)+c2Ln(Uk{|uk|m})+αkc2Ln(UUk{|uk|m})Fk(uk,U)+c2Ln(Uk{|˜uk|m})+αkc2Ln(Ω)Fk(uk,U)+c2Ln(U{|˜uk|m})+αkc2Ln(Ω)

    where ˜uk is as in Definition 2.2 and thus ˜uku in L1(Ω). Therefore, taking the liminf as k+ gives

    lim infk+Fk(umk,U)F(u,U)+c2Ln(U{|u|m}),

    hence (2.12) follows by the definition of Γ-liminf, taking into account that umkum, and by the Chebyshev inequality. Therefore, (2.12) ensures that F "almost" decreases by truncation up to an error which becomes small for m large. Inequality (2.12) is then enough to to carry out the Γ-convergence analysis below (cf. [12]). Hence, if αk is infinitesimal assumption (H2) can be dropped.

    However, if the sequence αk is uniformly bounded from below, we have no control on the term αkc2Ln(UUk{|uk|m}), therefore from (2.11) we cannot infer (2.12). Since with we want to study the Γ convergence of Fk for any choices of αk[0,1], assumption (H2) is actually necessary.

    In general the set functions F(u,) and F"(u,) are not inner regular. Then we consider their inner regular envelopes defined as:

    F(u,U):=sup{F(u,V):V⊂⊂U,VA(Ω)}.

    and

    F"(u,U):=sup{F"(u,V):V⊂⊂U,VA(Ω)}.

    Remark 2.11 (Properties of F,F"). The functionals F and F" are lower semicontinuous with respect to the convergence in Definition 2.2 [15,Remark 15.10], local [15,Remark 15.25], and it is immediate to check that they decrease by truncation. Furthermore, as set functions, they are both increasing and F is superadditive [15,Remark 15.10].

    The following compactness result is the analogue of [15,Theorem 16.9], when the sequential notion of Γ-convergence in Definition 2.5 is considered. We omit its proof since it is standard.

    Proposition 2.12 (Compactness by Γ-convergence). Let Fk be the localised functionals as in (2.9). Then there exists a subsequence (Fkj)(Fk) such that the corresponding functionals F and F" defined in (2.10) satisfy F=F".

    We now set

    F:=F=F". (2.13)

    In what follows we show that actually F coincides with the Γ-limit of the subsequence (Fkj). To this end we start noticing that by monotonicity we always have F"=FFF". Therefore, if we show that F"=F"; i.e., that F" is inner regular, we immediately get F=F"=F and therefore that Fkj(,U) Γ convergesto {\mathcal F}(\cdot, U) forevery U\in A(\Omega) $, as desired.

    A crucial preliminary result needed to prove the inner-regularity of F" is the so-called fundamental estimate, which has to hold uniformly in k. Since the Γ-limit is computed with respect to the convergence in Definition 2.2, the fundamental estimate we need is non-standard. Namely, we have to prove that the error in the fundamental estimate tends to zero when uku. This is achieved by first showing that the error goes like ukuLp(Ωk) and then by resorting to a truncation argument.

    We notice that an analogous estimate for degenerate functionals defined in Sobolev spaces can be found in [10,Proposition 3.3]. Whereas in the SBV-setting, for functionals of Mumford-Shah type with degenerate surface energy it can be found in the recent [19,Lemma 4.4].

    Following [10] we start showing how to construct suitable cut-off functions which are constant in iZnεk(K+i). To this end let δ>0 be small enough so that the set Kδ:={xRn:dist(x,K)<δ} satisfies Kδ⊂⊂Q. Let ψC0(Q) be a cut-off function between K and Kδ (that is 0ψ1, ψ1 on K, and sptψKδ) such that |ψ|2η.

    For kN and iZn, we define the operator Rki:W1,loc(Rn)W1,loc(Rn) as

    Rki(ϕ)(x):=(1ψ(xεki))ϕ(x)+ψ(xεki)εkKδ+εkiϕ(y)dy.

    By definition we have that

    Rki(ϕ)(x)=ϕ(x)if xεkKδ+εki,

    while Rki is constant in εkK+εki, namely we have

    Rki(ϕ)(x)=εkKδ+εkiϕ(y)dyif xεkK+εki.

    Finally, we consider the operator Rk:W1,loc(Rn)W1,loc(Rn) defined as

    Rk(ϕ)(x):={Rki(ϕ)(x)ifxεkKδ+εki,iZn,ϕ(x)otherwise.

    Let URn be open and bounded and let ϕW1,(U) then Rk(ϕ) is uniformly bounded in k. More precisely, we have

    Rk(ϕ)L(U;Rn)(2δd+1)ϕL(U;Rn), (2.14)

    where d denotes the diameter of Kδ. In fact,

    Rk(ϕ)L(U;Rn)2εkδsupiϕεKδ+εiϕ(y)dyL(εkKδ+εki;Rn)+ϕL(U;Rn)

    and

    ϕεkKδ+εkiϕ(y)dyL(εkKδ+εki;Rn)εkdϕL(U;Rn).

    In the next proposition we make use of the operator Rk to construct cut-off functions whose gradient vanishes in RnεkE; these cut-off functions are then used to prove the desired fundamental estimate.

    Proposition 2.13 (Fundamental estimate). For every η>0, and for every U, U", VA(Ω), with U⊂⊂U", there exist two constants M(η)>0 and kηN satisfying the following property: For every k>kη, for every uL1(Ω) with uSBVp(U"), and for every vL1(Ω) with vSBVp(V), there exists a function φC0(Ω) with φ=1 in a neighbourhood of U, sptφU" and 0φ1 such that

    Fk(φu+(1φ)v,UV)(1+η)(Fk(u,U")+Fk(v,V))+M(η)||uv||Lp(SεkE) (2.15)

    where S:=(U"U)V.

    Proof. Let U,U",VA(Ω) be as in the statement. Let η>0 be fixed and choose NN in a way such that

    1Nmax{c2Ln((U"U)V),3p1max{c2c1,c4c3}}<η. (2.16)

    Let moreover UA(Ω) be such that U⊂⊂U⊂⊂U" and consider the open sets

    U⊂⊂U1⊂⊂⊂⊂U3N⊂⊂U"

    where

    Ul:={x:dist(x,U)<dist(U,U)3Nl},for everyl=1,,3N.

    We notice that by definition of Ul we have that

    dist(Ul,Ul+1)=13N,for everyl=1,,3N1. (2.17)

    For every j=0,,N1 let ϕj be a cut-off function between U3j+1 and U3j+2 with |ϕj|<4N.

    Let kηN be such that

    εkd<13Nfor everyk>kη, (2.18)

    where d:=diam(Kδ)<2.

    If iZn is such that (εkK+εki)Ul for every k>kη, then thanks to (2.17)–(2.18) we can deduce that (εkK+εki)(RnUl+1)=. Therefore the functions φj:=Rk(ϕj) are cut-off functions between the sets U3j and U3(j+1), for every j=0,,N1 (where we have set U0:=U).

    Now let uSBVp(U") and vSBVp(V); for every j=0,,N1 fixed we have

    Fk(φju+(1φj)v,UV)=Fk(u,(UV)¯U3j)+Fk(v,VU3(j+1))+Fk(φju+(1φj)v,V(U3(j+1)¯U3j))Fk(u,U")+Fk(v,V)+Fk(φju+(1φj)v,V(U3(j+1)¯U3j)).

    We set

    wj:=φju+(1φj)v,Sj:=V(U3(j+1)¯U3j)

    and estimate the term Fk(wj,Sj). We clearly have

    Fk(wj,Sj)=Fk(wj,SjεkE)+Fk(wj,SjεkE). (2.19)

    By construction φj=0 in RnεkE, therefore appealing to (2.1) and (2.2) we deduce

    Fk(wj,SjεkE)=αkSjεkEfk(x,φju+(1φj)v)dx+βk(SjεkE)Swjgk(x,νwj)dHn1c2αk(Ln(SjεkE)+SjεkE|u|pdx+SjεkE|v|pdx)+c4βk(Hn1((SjεkE)Su)+Hn1((SjεkE)Sv))c2αkLn(SjεkE)+c2c1(αkSjεkEfk(x,u)dx+αkSjεkEfk(x,v)dx)+c4c3(βk(SjεkE)Sugk(x,νu)dHn1+βk(SjεkE)Svgk(x,νv)dHn1)c2αkLn(SjεkE)+max{c2c1,c4c3}(Fk(u,SjεkE)+Fk(v,SjεkE)). (2.20)

    Moreover, again invoking (2.1) and (2.2), in εkE we have

    Fk(wj,SjεkE)=SjεkEfk(x,wj)dx+SjεkESwjgk(x,νwj)dHn1c2(Ln(SjεkE)+SjεkE|φj(uv)+φju+(1φj)v|pdx)+c4(Hn1((SjεkE)Su)+Hn1((SjεkE)Sv))c2Ln(SjεkE)+c23p1φjpL(U;Rn)SjεkE|uv|pdx+3p1max{c2c1,c4c3}(Fk(u,SjεkE)+Fk(v,SjεkE)). (2.21)

    Since |ϕj|4N, combining the definition of φj with (2.14) gives

    φjL(U;Rn)(2δd+1)4N. (2.22)

    In view of (2.19), by gathering (2.20)–(2.22) we then obtain for every j=0,,N1

    Fk(wj,Sj)c2Ln(Sj)+3p1max{c2c1,c4c3}(Fk(u,Sj)+Fk(v,Sj))+M(η)SjεkE|uv|pdx,

    where

    M(η):=c23p1(2δd+1)p(4N)p.

    Therefore there exists j{0,,N1} such that

    Fk(wj,Sj)1NN1i=0Fk(wj,Sj)c2NLn((U"U)V)+3p1Nmax{c2c1,c4c3}(Fk(u,U")+Fk(v,V))+M(η)(U"U)VεkE|uv|pdx.

    Finally the thesis follows from (2.16) by choosing φj as a cut-off function and setting S:=(U"U)V.

    Thanks to the fundamental estimate Proposition 2.13 we are now able to prove the following abstract Γ-convergence result for the sequence of localised functionals Fkj(,U).

    Theorem 2.14 (Abstract Γ-convergence and properties of the Γ-limit). Let F be as in (2.13), then:

    a. (locality and lower semicontinuity) for every UA(Ω), the functional F(,U) is local and lower semicontinuous with respect to the L1(Ω)-convergence;

    b. (measure property) for every uGSBVp(Ω), the set function F(u,) is the restriction to A(Ω) of a Radon measure on Ω;

    c. (Γ-convergence) for every UA(Ω) it holds F(,U)=F(,U)=F"(,U)on GSBVp(Ω)

    d. (translational invariance in u) for every uL1(Ω) and UA(Ω) there holds F(u+s,U)=F(u,U) for every sR.

    Proof. Since the L1(Ω)-convergence implies the convergence in the sense of Definition 2.2, property 1 immediately follows from Remark 2.11. In view of Remark 2.11, property 2 follows by the De Giorgi and Letta criterion (see, e.g., [15,Theorem 14.23]) once we show that for every uGSBVp(Ω) the set function F(u,) is subadditive. In its turn, the subadditivity of F(u,) follows from Proposition 2.13. Since in our setting this proof is not entirely standard, we discuss it in detail for the readers' convenience.

    We start observing that on GSBVp(Ω) the functional F satisfies the following limsup-type inequality: For every uGSBVp(Ω) and for every U,UA(Ω) with U⊂⊂U, there exists a sequence (uj)GSBVp(U)L1(Ω) with uju such that

    lim supj+Fkj(uj,U)F(u,U)

    (see, e.g., [15,Proposition 16.4 and Remark 16.5] also recalling that the infimum in the definition of F" is actually attained).

    Now let U,VA(Ω) and let uGSBVp(Ω)L(Ω). Fix any U⊂⊂U, V⊂⊂V, U,VA(Ω). Choose an open set U" such that U⊂⊂U"⊂⊂U and two sequences (uj)GSBVp(U")L1(Ω) and (vj)GSBVp(V)L1(Ω), with uju and vju such that

    lim supj+Fkj(uj,U")F(u,U),lim supj+Fkj(vj,V)F(u,V). (2.23)

    Since the functionals Fk decrease by truncation, we can additionally assume that ujL(Ω),vjL(Ω)uL(Ω); clearly,

    limj+ujuLp(Ωkj)=limj+vjuLp(Ωkj)=0. (2.24)

    Let η>0 be fixed and arbitrary. The fundamental estimate Proposition 2.13 provides us with constants M(η)>0 and jηN and with a sequence (φj) of cut-off functions between U and U" such that

    Fkj(φjuj+(1φj)vj,UV)(1+η)(Fkj(uj,U")+Fkj(vj,V))+M(η)ujvjpLp(Ωkj)

    for every jjη. Hence appealing to (2.23), to the convergence φjuj+(1φj)vju, and to the obvious inequality FF, by taking the limit as j+, we get

    F(u,UV)(1+η)(F(u,U)+F(u,V)).

    Now letting η0, and then UU, VV in view of the inner-regularity of F we get

    F(u,UV)F(u,U)+F(u,V), (2.25)

    hence the subadditivity of F(u,) for uGSBVp(Ω)L(Ω).

    Now let uGSBVp(Ω) and, for every mN, set um:=(um)(m). Then, since F decreases by truncation (2.25) immediately gives

    F(um,UV)F(u,U)+F(u,V).

    Then, taking the limit as m+, in view of the convergence umu in L1(Ω) and the lower semicontinuity of F we obtain

    F(u,UV)lim infm+F(um,UV)F(u,U)+F(u,V),

    and thus the subadditivity of F(u,) for every uGSBVp(Ω).

    The proof of property 3 is achieved by showing that F" is inner-regular. Indeed, this is equivalent to F"=F", which by definition of F implies F"≤FF. Since clearly FF", we actually deduce that F is the Γ-limit of Fkj.

    The inner regularity of F" follows from the fundamental estimate Proposition 2.13. To see this, for every UA(Ω) let G(,U) be the localised version of the functional G defined in (2.8); i.e.,

    G(u,U):={c2U(1+|u|p)dx+c4Hn1(SuU)ifuGSBVp(U),+otherwise in L1(Ω). (2.26)

    Now fix WA(Ω) and uGSBVp(Ω); since G(u,) is the restriction to A(Ω) of a Radon measure, for every η>0 there exists a compact set ˜WW such that and MS(u,W˜W)<η.

    Now choose U,UA(Ω) satisfying ˜WU⊂⊂U⊂⊂W and set V:=W˜W. Recalling that F"(u,) is increasing, appealing to Proposition 2.13 easily gives

    F"(u,W)F"(u,UV)F"(u,U)+F"(u,V)=F"(u,U)+F"(u,W˜W).

    Recalling that F"≤G, by taking the sup on U⊂⊂W we get

    F"(u,W)F"(u,W)+G(u,W˜W)F"(u,W)+η.

    Hence, by the arbitrariness of η>0 we get F"(u,W)F"(u,W) for every WA(Ω) and every uGSBVp(Ω). Since the opposite inequality is always satisfied, we readily deduce the inner regularity of F"(u,), as desired.

    Eventually, the proof of property 4 is standard and follows as in, e.g., [9,Lemma 3.7].

    In the following theorem we show that the Γ-limit F can be represented in an integral form as a free-discontinuity functional. Moreover, thanks to [5,Theorem 4] the functional F turns out to be non-degenerate, unlike the functionals Fk.

    Theorem 2.15 (Integral representation of the Γ-limit). Let F be as in Theorem 2.14. Then, there exist a Carathéodory function f:Rn×Rn[0,+) and a Borel function g:Rn×Sn1(0,+) such that

    F(u,U)=Uf(x,u)dx+SuUg(x,[u],νu)dHn1 (2.27)

    for every uGSBVp(Ω) and every UA(Ω).

    Furthermore, the function f:Rn×Rn[0,+) satisfies the following properties:

    i) (convexity in ξ) for a.e. xRn, f(x,) is convex;

    ii) (p growth and coercivity) there exists ˜c1>0 such that for a.e. xRn and for every ξRn it holds

    ˜c1|ξ|pf(x,ξ)c2(1+|ξ|p), (2.28)

    where c2 is as in (2.1).

    The function g:Rn×R×Sn1[0,+) satisfies the following properties:

    iii) (monotonicity in t and symmetry) for a.e. xRn and for every νSn1, g(x,,ν) is nondecreasing on (0,+) and satisfies the symmetry condition g(x,t,ν)=g(x,t,ν) for every tR;

    iv) (subadditivity in t) for a.e. xRn and for every νSn1

    g(x,t1+t2,ν)g(x,t1,ν)+g(x,t2,ν),

    for every t1,t2R;

    v) (convexity in ν) for a.e. xRn and for every tR, the 1-homogeneous extension of g(x,t,) to Rn is convex. Equivalently, for a.e. xRn and for every tR the function g satisfies

    g(x,t,ν)λ1g(x,t,ν1)+λ2g(x,t,ν2),

    for every ν,ν1,ν2Sn1, λ1,λ20 such that λ1ν1+λ2ν2=ν;

    vi) (bounds) there exists ˜c3>0 such that for a.e. xRn, for every tR, and every νSn1 it holds

    ˜c3g(x,t,ν)c4, (2.29)

    where c4 is as in (2.2).

    Proof. Let Ek:L1(Ω)×A(Ω)[0,+] be the functionals defined as

    Ek(u,U):={c1Uk|u|pdx+c3Hn1(SuUk)ifuSBVp(U)+otherwise in L1(Ω), (2.30)

    with c1 and c3 as in (2.1) and (2.2), respectively. Appealing to [5,Theorem 4] and also noticing that the Lp-convergence in the statement can be equivalently replaced by the convergence in Definition 2.2, we deduce that Ek(,U) Γ-converges to E(,U) for every UA(Ω), where

    E(u,U)=Uˆf(u)dx+SuUˆg(ν)dHn1

    with ˆf and ˆg as in [5,Theorem 4] formulas (40) and (41), respectively. Moreover ˆf and ˆg satisfy

    ˜c1|ξ|pˆf(ξ)for every ξRnand˜c3ˆg(ν)for every νSn1,

    for some ˜c1,˜c3>0. Then, since EkFk, we may deduce that for every uSBVp(Ω) and every UA(Ω) we have

    E(u,U)F(u,U). (2.31)

    We recall that for every uSBVp(Ω) and every UA(Ω) we also have

    F(u,U)G(u,U), (2.32)

    where G is as in (2.26).

    Now let σ>0 and for every uSBVp(Ω) and UA(Ω) set

    Fσ(u,U):=F(u,U)+σSuU|[u]|dHn1.

    For every fixed σ>0 the functional Fσ satisfies properties 1, 2, and 4 in Theorem 2.14. Moreover, in view of (2.31)–(2.32) it holds

    ˜c1U|u|pdx+SuU(˜c3+σ|[u]|)dHn1Fσ(u,U)c2U(1+|u|p)dx+SuU(c4+σ|[u]|)dHn1

    Therefore, we can invoke the integral representation result [7,Theorem 1] to deduce that for every uSBVp(Ω) and every UA(Ω) we have

    Fσ(u,U)=Ufσ(x,u)dx+SuUgσ(x,[u],νu)dHn1,

    where fσ and gσ are given by the following derivation formulas

    fσ(x,ξ):=lim supρ0+1ρninf{Fσ(u,Qρ(x)):uSBVp(Qρ(x)),u=uξnearQρ(x)} (2.33)

    and

    gσ(x,t,ν):=lim supρ0+1ρn1inf{Fσ(u,Qνρ(x)):uSBVp(Qνρ(x)),u=ut,νxnearQνρ(x)}. (2.34)

    By (2.33) and (2.34) the sequences (fσ)σ>0 and (gσ)σ>0 are decreasing as σ decreases, therefore by setting f:=limσ0+fσ and g:=limσ0+gσ, by the pointwise convergence of (Fσ)σ>0 to F and the Monotone Convergence Theorem, we get

    F(u,U)=Uf(x,u)dx+SuUg(x,[u],νu)dHn1,

    for every uSBVp(Ω) and UA(Ω). Eventually, a standard truncation and continuity argument allows to extend this integral representation to the whole space GSBVp(Ω) and thus to get exactly (2.27).

    The measurability properties of f and g follow from the derivation formulas (2.33) and (2.34), arguing as in the appendix of [12]. The convexity of f in ξ, the subadditivity of g in t, and the convexity in ν of its 1-homogeneous extension are immediate consequences of the L1(Ω)-lower semicontinuity of F.

    To show that f and g satisfy, respectively, the lower bounds as in ii) and vi) we argue as follows. Set

    Φσ(u,U):={˜c1U|u|pdx+SuU(˜c3+σ|[u]|)dHn1ifuSBVp(U)+otherwise in L1(Ω),

    and for every xRn and ξRn define

    ϕσ(x,ξ):=lim supρ0+1ρninf{Φσ(u,Qρ(x)):uSBVp(Qρ(x)),u=uξnearQρ(x)},

    while for every xRn, tR, and νSn1 set

    ψσ(x,t,ν):=lim supρ0+1ρn1inf{Φσ(u,Qνρ(x)):uSBVp(Qνρ(x)),u=uν,txnearQνρ(x)}.

    Since ΦσFσ on SBVp(Ω) we clearly have both ϕσfσ and ψσgσ. We now show that ϕσ(x,ξ)=˜c1|ξ|p for every xRn and every ξRn and ψσ(x,t,ν)=˜c3+σt. To do so we notice that by the homogeneity in x of Φσ, we have both ϕσ(x,ξ)=ϕσ(0,ξ) for every xRn and every ξRn and ψσ(x,t,ν)=ψσ(0,t,ν). We can now apply the integral representation result [7,Theorem 1] to Φσ so that choosing u=uξ and U=Q we obtain

    ˜c1|ξ|p=Φσ(uξ,Q)=Qϕσ(y,ξ)dy=ϕσ(0,ξ)=ϕσ(x,ξ),

    while choosing u=ut,ν0 and U=Qν we obtain

    ˜c3+σt=Φσ(ut,ν0,Qν)=ΠνQνψσ(y,t,ν)dHn1=ψσ(0,t,ν)=ψσ(x,t,ν),

    and hence the desired equalities. Therefore we deduce

    ˜c1|ξ|p=ϕσ(x,ξ)fσ(x,ξ)for every x,ξRn

    which immediately gives the lower bound of f; moreover there holds

    ˜c3˜c3+σt=ψσ(x,t,ν)gσ(x,t,ν)for every xRn,tR,νSn1

    hence, taking the inf on σ>0 yields the the lower bound on g.

    The upper bound in ii) immediately follows from (2.33) and the obvious inquality Fσ(uξ,Qρ(x))ρnc2(1+|ξ|p), while the upper bound in vi) follows from (2.34) and

    F(ut,νx,Qνρ(x))Fσ(ut,νx,Qνρ(x))ρn1(c4+σt),

    which holds true for every σ>0 and hence also in the limit as σ0+.

    Finally, the monotonicity in t and the symmetry of g easily follow from (2.34).

    Theorem 2.16 (Γ-convergence). Let Fk be the functionals defined in (2.3). Then, there exists a subsequence kj+ such that (Fkj) Γ-converges to the functional F given by (2.27), for some Carathéodory function f:Rn×Rn[0,+) and some Borel function g:Rn×Sn1(0,+) satisfying properties (i)(vi) as in Theorem 2.15.

    Proof. The proof is an immediate consequence of Theorem 2.14 and Theorem 2.15.

    Corollary 2.17 (Γ-convergence of porous brittle materials). Let αk=βk=0 and let Fk be the corresponding functionals given by (2.3). Then, there exists a subsequence kj+ such that (Fkj) Γ-converges with respect to the L1(Ω)-convergence to the functional F given by (2.27), for some Carathéodory function f:Rn×Rn[0,+) and some Borel function g:Rn×Sn1(0,+) satisfying properties (i)(vi) as in Theorem 2.15.

    Proof. Since the L1(Ω)-convergence implies the convergence in Definition 2.2, the proof of the liminf inequality is immediate from Theorem 2.16. Now let uGSBVp(Ω), then by Theorem 2.16 there exists (uj)L1(Ω) such that uju and limjFkj(uj)=F(u). In view of Definition 2.2 this means that there exists a sequence (˜uj)L1(Ω) such that ˜uj=uj a.e. in Ωkj and ˜uju in L1(Ω). Then, since the choice αkj=βkj=0 implies the equality Fkj(˜uj)=Fkj(uj), the sequence (˜uj) is the desired recovery sequence.

    On account of the Γ-convergence result Theorem 2.14 in this section we establish a convergence result for minimisation problems associated to a suitable perturbation of the functionals Fk. To this end, let hL(Ω) and for every k set

    Mk:=inf{Fk(u)+uhpLp(Ωk):uL1(Ω)}.

    By a standard truncation argument it is immediate to show that

    Mk=inf{Fk(u)+uhpLp(Ωk):uSBVp(Ω),uL(Ω)hL(Ω)}. (2.35)

    Proposition 2.18. Let F=Γ-limjFkj and let (uj)SBVp(Ω) be such that

    limj+(Fkj(uj)+ujhpLp(Ωkj)Mj)=0. (2.36)

    Then, up to subsequences (not relabelled), uj converges in the sense of Definition 2.2 to a function ˉuSBVp(Ω)L(Ω) which solves

    M:=min{F(u)+C(K)uhpLp(Ω):uSBVp(Ω),uL(Ω)hL(Ω)},

    where C(K):=Ln(QK). Moreover it holds MjM, as j+.

    Proof. Let (uj)SBVp(Ω) be as in (2.36). Then, in view of (2.35), (H1), and (H3) we have

    supj(ujL(Ω)+Ωkj|uj|pdx+Hn1(SujΩkj))<+.

    Therefore Theorem 2.1 yields the existence of a function ˉuSBVp(Ω)L(Ω) and of a sequence (˜uj)SBVp(Ω) with ˜uj=uj a.e. in Ωkj such that (up to subsequences) ˜ujˉu in Lp(Ω), moreover ˉuL(Ω)hL(Ω). We have

    C(K)ˉuhpLp(Ω)=limj+(˜ujh)χΩkjpLp(Ω)=limj+ujhpLp(Ωkj),

    thus by Theorem 2.14 we get

    F(ˉu)+C(K)ˉuhpLp(Ω)lim infj+(Fkj(uj)+ujhpLp(Ωkj)).

    Therefore, by definition of uj we obtain

    F(ˉu)+C(K)ˉuhpLp(Ω)lim infj+Mj. (2.37)

    Now let wSBVp(Ω)L(Ω) be an arbitrary function such that wL(Ω)hL(Ω). Again appealing to Theorem 2.14 we can find (wj)L1(Ω) such that wjw and limjFj(wj)=F(w). Now let ˜wj be as in Definition 2.2, let m:=hL(Ω) and denote with (˜wmj) the sequence of truncated functions of (˜wj) at level m. We clearly have ˜wmj=wmj a.e. in Ωkj and ˜wmjw in Lp(Ω). Hence

    limj+wmjhpLp(Ωkj)=limj+(˜wmjh)χΩkjpLp(Ω)=C(K)whpLp(Ω).

    Moreover, since lim supjFj(wmj)F(w), we immediately deduce

    lim supj+MjF(w)+C(K)whpLp(Ω). (2.38)

    Finally, by gathering (2.37) and (2.38) we obtain

    F(ˉu)+C(K)ˉuhpLp(Ω)lim infj+Mjlim supj+MjF(w)+C(K)whpLp(Ω),

    hence by the arbitrariness of w we deduce that ˉu is a minimiser for F+C(K)hpLp(Ω). Finally, taking w=ˉu also implies MjM. Since moreover this limit does not depend on the subsequence, the convergence holds true for the whole (Mj).

    In this section we restrict the analysis to the case of εk-periodic integrands fk and gk. That is, we consider the functionals Fαk,βkk:L1(Ω)[0,+] defined as

    Fαk,βkk(u):={Ωkf(xεk,u)dx+αkΩΩkf(xεk,u)dx+SuΩkg(xεk,νu)dHn1+βkSu(ΩΩk)g(xεk,νu)dHn1ifuSBVp(Ω)+otherwise, (3.1)

    where f and g are Q-periodic in the first variable and satisfy (H1)–(H2) and (H3)–(H4), respectively.

    With the help of some specific examples, which correspond to some specific choices of f, g, and Ωk, we show that the Γ-limit of Fαk,βkk is highly sensitive both to the choice of the coefficients αk, βk and to the asymptotic behaviour of αk, βk compared to the period of the microstructure εk. The examples we are going to discuss are taken from Barchiesi and Focardi [5] (see also Cagnetti and Scardia [14] and Focardi, Gelli and Ponsiglione [17]), from Barchiesi, Lazzaroni and Zeppieri [6], and from Pellet, Scardia and Zeppieri [19]. For the corresponding proofs we refer the reader to the aforementioned papers.

    In this subsection we consider the limit case αk=βk=0; i.e., we consider the functionals

    F0,0k(u):={Ωkf(xεk,u)dx+SuΩkg(xεk,νu)dHn1ifuSBVp(Ω),+otherwise in L1(Ω). (3.2)

    Loosely speaking, in this case the soft or weak inclusions in the material are replaced by perforations [5,14,17].

    Theorem 3.1 (Homogenisation of periodic porous brittle materials). Let F0,0k be the functionals as in (3.2). Then (F0,0k) Γ-converges both with respect to the convergence in Definition 2.2 and with respect to the L1(Ω)-convergence to the functional F0 which is finite on GSBVp(Ω) and given by

    F0(u)=Ωf0(u)dx+Sug0(νu)dHn1, (3.3)

    where f0 and g0 are, respectively, given by the following homogenisation formulas

    f0(ξ)=inf{QEf(y,u)dx:uW1,p(QE),u=uξnearQ}, (3.4)

    for every ξRn, whereas

    g0(ν)=limT+1Tn1inf{SuTQνEg(y,νu)dHn1:uP(TQνE),u=uν,10nearTQν}, (3.5)

    for every νSn1.

    Proof. Theorem 2.16 and Corollary 2.17 yield the existence of a subsequence kj+ such that the corresponding functionals F0,0kj Γ-converge to F as in (2.27), both with respect to the convergence in Definition 2.2 and to the L1(Ω)-convergence. Then, the homogenisation formulas (3.4) and (3.5) together with the identity F=F0 follow from [5,Theorem 4]. Finally, since (3.4) and (3.5) are subsequence-independent, invoking the Urysohn property [15,Proposition 8.3] readily implies the Γ-convergence of the whole sequence (F0,0k) to F0.

    The following result is an immediate consequence of Theorem 3.1 and of an adaptation of the Cagnetti and Scardia extension result [14,Theorem 1.3] to the case of a general exponent p>1.

    Corollary 3.2. Let αk,βk0 and let Fαk,βkk be the corresponding functionals as in (2.3). Then, the sequence (Fαk,βkk) Γ-converges to the functional F0 given by (3.3).

    Proof. By Theorem 2.16 (up to subsequences not relabelled) the functionals Fαk,βkk Γ-converge to F. Since F0,0kFαk,βkk by Theorem 3.1 we immediately get F0F.

    We now prove the converse inequality. To this end let uSBVp(Ω)L(Ω) and let (uk)SBVp(Ωk) be a recovery sequence for F0,0k. That is uku and limkF0,0k(uk)=F(u). Since the functionals F0,0k decrease by truncation it is not restrictive to assume that ukL(Ω)uL(Ω). Starting from (uk) we now want to construct a sequence (vk) which both satisfies vku in L1(Ω) and limkFαk,βkk(vk)=F0(u). To this end, we start noticing that the bounds (2.1) and (2.2) readily imply

    supk(Ωk|uk|pdx+Hn1(SukΩk))<+. (3.6)

    For every fixed k let vk:=TkukSBVp(Ω) be the extended function of uk to Ω whose existence is given by [14,Theorem 1.3]; i.e., vk is such that vk=uk a.e. in Ωk, ukL(Ω)uL(Ω), and

    Ω|vk|pdx+Hn1(SvkΩ)C(Ωk|uk|pdx+Hn1(SukΩk)) (3.7)

    for some C>0 independent of k. By definition of vk, also invoking the Ambrosio compactness Theorem, it is immediate to check that vku in L1(Ω).

    By (3.7) we get that

    αkΩΩk|vk|pdx+βkHn1(Svk(ΩΩk))αkΩ|vk|pdx+αkHn1(SvkΩ)Cmax{αk,βk}(Ωk|uk|pdx+Hn1(SukΩk)), (3.8)

    where (3.8) is infinitesimal thanks to (3.6), since max{αk,βk}0 as k+. Thus eventually

    limk+Fαk,βkk(vk)=limk+F0,0k(uk)=F0(u),

    hence (vk) is the desired sequence. Therefore, by the Γ-convergence of Fαk,βkk to F we can deduce that for every uSBVp(Ω)L(Ω) it holds F(u)F0(u).

    Now let uGSBVp(Ω) and denote with um its truncated function at level m>0. We clearly have F(um)F0(um)F0(u), hence the desired inequality follows by the L1(Ω) convergence um to u and by the lower semicontinuity of F.

    The following remarks are in order.

    Remark 3.3. In view of Remark 2.10, both in Theorem 3.1 and in Corollary 3.2 assumption (H2) on f can be dropped.

    Remark 3.4 (On f0). The homogenised volume energy density f0 given by (3.4) is the same as that obtained by Acerbi, Chiadó-Piat, Dal Maso and Percivale [1] in the case of elastic perforated materials. Moreover, it is easy to check that if f is p-homogeneous then the corresponding f0 given by (3.4) is also p-homogeneous.

    Remark 3.5 (Energy decoupling). In spite of the strong degeneracy of the integrands in (3.2) (resp. in (3.1)), which in this case are identically equal to zero (resp. both infinitesimal) in the εk-periodic set ΩΩk, Theorem 3.1 (resp. Corollary 3.2) shows that the functionals F0,0k (resp. Fαk,βkk) exhibit a limit behaviour which is qualitatively similar to that of free-discontinuity functionals with coercive integrands [9,13,18]. Namely, in the homogenised limit there is no interaction between bulk and surface term. As a consequence the homogenised surface energy density g0 does not depend on t, and therefore the Γ-limit is of brittle type.

    In this section we show that if only one of the coefficients αk and βk is infinitesimal (while the other stays uniformly bounded from below), then the asymptotic behaviour of the functionals Fαk,βkk can be very different from that of F0,0k (or of Fαk0,βk0k). In particular, we show that in this case a volume-surface energy coupling cannot be excluded in general. To do so we exhibit coefficients αk,βk, integrands f,g and a geometry for the periodic set E which give rise to the desired limit coupling. This is done by resorting to the analysis of Barchiesi, Lazzaroni and Zeppieri [6] and Pellet, Scardia and Zeppieri [19], which is briefly reviewed in Subsection 3.2.1 and Subsection 3.2.2, respectively.

    The functionals analysed in [6] and [19] are both of Mumford-Shah type and can be written in the form

    MSαk,βkk(u)=Ωak(xεk)|u|2dx+Subk(xεk)dHn1,uSBV2(Ω) (3.9)

    where ak,bk:Rn[0,1] are Q-periodic functions and in the periodicity cell Q are defined as

    ak(y)={αkify¯Qr1ifyQ¯Qrbk(y)={βkify¯Qr1ifyQ¯Qr (3.10)

    with r(0,1). From (3.9)–(3.10) we infer that in this case f=f(ξ)=|ξ|2, g1, and Ωk=ΩεkE with E=RniZn(¯Qr+i).

    Remark 3.6 (Mumford-Shah functional in perforated domains). The choice αk=βk=0 corresponds to the Mumford-Shah functionals in perforated domain. The functional MS0,0k is a special instance of (3.2) and its homogenised limit is treated in [14,17] for general sets E. In this case the homogenised integrands (3.4) and (3.5) reduce, respectively, to

    f0(ξ)=inf{Q¯Qr|u|2dx:uW1,2(Q¯Qr),u=uξ near Q}, (3.11)

    for every ξRn, and to

    g0(ν)=limT+1Tn1inf{Hn1(SuTQνE):uP(TQνE),u=uν,10 near TQν}, (3.12)

    for every νSn1. From (3.11) and (3.12) it is easy to check that f0(ξ)=A0ξξ, for some A0Rn×n which satisfies ˜c1IA0I, in the sense of quadratic forms (cf. (2.28)). Hence, f0 is a positive quadratic form. Moreover, it holds g0(ei)=1rn1, for every i=1,,n.

    We consider the case αk0 and βk=1 which models the situation where the periodic set ΩΩk is occupied by a brittle material with a very small elastic modulus. For this reason, we refer to the set ΩΩk as the set of soft inclusions. With this choice the functionals in (3.9) become

    MSαk,1k(u)=Ωk|u|2dx+αkΩΩk|u|2dx+Hn1(Su),uSBV2(Ω). (3.13)

    In [6] Barchiesi, Lazzaroni and Zeppieri showed that the asymptotic behaviour of MSαk,1k heavily depends on the mutual vanishing rate of αk and εk; that is, it depends on the parameter

    :=limk+αkεk[0,+]. (3.14)

    For the proof of Theorem 3.7 below we refer the reader to [6,Theorems 1,4,and Remark 6].

    Theorem 3.7 (Homogenisation of periodic brittle materials with soft inclusions). Let MSαk,1k be the functionals defined in (3.13) and let [0,+] be as in (3.14). Then, up to subsequences not relabelled, (MSαk,1k) Γ-converges to the functional F which is finite on GSBV2(Ω) and given by

    F(u)=Ωf0(u)dx+Sug([u],νu)dHn1, (3.15)

    where f0 is as in (3.11) and for every tR, νSn1

    g(t,ν)={g0(ν)if=01if=+.

    Moreover for every (0,+) it holds

    min{g0(ei)+ct2,1}g(t,ei)min{g0(ei)+ˆct,1} (3.16)

    for every t>0, i=1,,n, and for some c,ˆc>0, with lim0+c=lim0+ˆc=0.

    Remark 3.8. The following remarks are in order.

    (i) As far as the homogenised volume energy is concerned, the soft inclusions are (energetically) equivalent to the perforations in the material.

    (ii) For =0, which corresponds to αkεk, the functionals MSαk,1k are equivalent to the functionals MS0,0k, in the sense of Γ-convergence.

    (iii) For (0,+) the bounds in (3.16) imply that, along the coordinate directions, g depends on t. Moreover it becomes constant (and equal to 1) above a certain threshold t0>0; i.e., g is of cohesive type. Being the microscopic energies MSαk,1k of brittle type, the cohesive behaviour of g can only be explained as the result of a non trivial bulk-surface coupling by homogenisation. This interaction is particularly apparent from the upper-bound construction in [6] which we briefly illustrate here in the case n=2.

    For i=1,2 we have g(t,ei)=F(uei,t0,Q), moreover it is immediate to check that g(t,e1)=g(t,e2). Clearly g(t,e2)1 for every t>0. Then, to get the upper bound in (3.16) it suffices to show that g(t,e2)g0(e2)+ˆct for some ˆc>0. Let RQR2 be the open rectangle defined as

    R:=(r2,r2)×(τ2,τ2),

    with τ(0,r) to be determined. Set

    Rk:=QiZ(εkR+(iεk,0))

    and let (uk)SBV2(Q) be the sequence of functions defined as

    uk(x):={tifxQRk and x20,t2+tτεkxnifxRk,0ifxQRk and x2<0,

    (see Figure 2). We clearly have ukue2,t0 in L1(Q); moreover

    Rk|uk|2dx(1εk+1)t2τandH1(Suk)εk(1εk+1)(1r+2τ),
    Figure 2.  Construction of the recovery sequence in an εk-cell, across the interface x2=0.

    therefore

    g(t,e2)=F(ue2,t0,Q)lim supk+MSαk,1k(uk,Q).1r+2τ+t2τ

    Hence, by optimising on τ we get

    g(t,e2)1r+22t (3.17)

    thus the desired estimate follows with ˆc=22, by recalling that g0(e2)=1r.

    Loosely speaking, the construction as above shows that, the cost of an elastic deformation of the soft inclusions is of the same order of the energy spent to create a microscopic crack. Since the former depends linearly on t (while the latter is constant in t) for small values of t, to approximate a macroscopic crack it can be convenient to combine microscopic deformations of the soft inclusions (with high gradients) and microscopic jumps.

    (iv) Even if not immediately apparent from the homogenisation formulas, a volume-surface interaction takes place for =0, as well. Indeed, in this case g=g0 whereas in MSαk,1k the surface energy density is identically equal to one. In this case in fact, the cost of an elastic deformation of the soft inclusions is negligible (cf. (3.17) for =0) so that to approximate a macroscopic crack it is never convenient to introduce microscopic cracks inside the soft material. On the contrary, in the regime =+, which corresponds to αkεk, there is a complete volume-surface decoupling, as in the coercive case.

    We consider the case αk=1 and αk0 which models the situation where the periodic set ΩΩk is occupied by a brittle material with a very small fracture resistance. For this reason, we refer to the set ΩΩk as the set of weak inclusions. With this choice the functionals in (3.9) become

    MS1,βkk(u)=Ω|u|2dx+Hn1(SuΩk)+βkHn1(Su(ΩΩk)),uSBV2(Ω). (3.18)

    In [19] Pellet, Scardia and Zeppieri showed that the asymptotic behaviour of MS1,βkk heavily depends on the mutual vanishing rate of βk and εk, that is on the parameter

    :=limkβkεk[0,+]. (3.19)

    For the proof of Theorem 3.9 below we refer the reader to [19].

    Theorem 3.9 (Homogenisation of periodic brittle materials with weak inclusions). Let MS1,βkk be the functionals defined in (3.18) and let [0,+] be as in (3.19). Then, up to subsequences not relabelled, (MS1,βkk) Γ-converges to the functional F which is finite on GSBV2(Ω) and given by

    F(u)=Ωf(u)dx+Sug0(νu)dHn1, (3.20)

    where g0 is as in (3.12) and for every ξRn

    f(ξ)={f0(ξ)if=0|ξ|2if=+.

    Moreover for every (0,+) it holds

    f0(ξ)f(ξ)min{|ξ|2,f0(ξ)+C} (3.21)

    for every ξRn and for some C>0.

    Remark 3.10 The following remarks are in order.

    (i) As far as the homogenised surface energy is concerned, the weak inclusions are (energetically) equivalent to the perforations in the material.

    (ii) For =0, which corresponds to βkεk, the functionals MS1,βkk are equivalent to the functionals MS0,0k, in the sense of Γ-convergence. Indeed, "removing the weak inclusions from the material" has an infinitesimal cost of order βk/εk given by the perimeter of the weak inclusions (proportional to βkεn1k) multiplied by εnk (the number of εk-cells contained in Ω). In this case a volume-surface energy coupling takes place since the elastic energy can be lowered by introducing cracks in the materials.

    (iii) For (0,+) the bounds in (3.21) hold true (see [19,Lemma 6.1]). The bound from below is immediate and it is a consequence of the trivial bound MS0,0kMS1,βkk. The bound from above shows that for large deformations; i.e., for large |ξ|, to approximate a macroscopic elastic deformation is energetically favourable to mix elastic deformations and jumps in the weak inclusions. Moreover, (3.21) implies that for |ξ| large it holds f(ξ)<|ξ|2. The latter shows that a stiffness degradation occurs in the homogenised limit, and that the macroscopic energy F describes a damaged material (the same being true for =0).

    (iv) The bounds in (3.21) combined with an easy scaling argument show that in the regime (0,+) the homogenised volume energy density f is not 2-homogeneous. Indeed, assume by contradiction that this is not the case and let λ0. Taking into account that f0 is 2-homogeneous (see Remark 3.4), we can replace in (3.21) ξ with λξ and divide by λ2 to get

    f0(ξ)f(ξ)min{|ξ|2,f0(ξ)+Cλ2}.

    Therefore by letting |λ|+ we get ff0 which leads to a contradiction in view of [19,Proposition 6.10].

    (v) In the regime =+, which corresponds to βkεk, there is a complete volume-surface decoupling, as in the coercive case. Loosely speaking, in this case the fracture resistance of the weak inclusions is not small enough to make cracks energetically more convenient than (or at least comparable to) elastic deformations.

    This work was supported by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under the Germany's Excellence Strategy EXC 2044-390685587, Mathematics Münster: Dynamics–Geometry–Structure.

    The author declares no conflict of interest.



    [1] Acerbi E, Chiadó Piat V, Dal Maso G, et al. (1992) An extension theorem from connected sets, and homogenization in general periodic domains. Nonlinear Anal Theor 18: 481-496. doi: 10.1016/0362-546X(92)90015-7
    [2] Ambrosio L, Fusco N, Pallara D (2000) Functions of Bounded Variations and Free Discontinuity Problems. Oxford: Clarendon Press.
    [3] Barchiesi M (2018) Toughening by crack deflection in the homogenization of brittle composites with soft inclusions. Arch Ration Mech Anal 227: 749-766. doi: 10.1007/s00205-017-1173-5
    [4] Barchiesi M, Dal Maso G (2009) Homogenization of fiber reinforced brittle materials: The extremal cases. SIAM J Math Anal 41: 1874-1889. doi: 10.1137/080744372
    [5] Barchiesi M, Focardi M (2011) Homogenization of the Neumann problem in perforated domains: An alternative approach. Calc Var Partial Dif 42: 257-288. doi: 10.1007/s00526-010-0387-2
    [6] Barchiesi M, Lazzaroni G, Zeppieri CI (2016) A bridging mechanism in the homogenization of brittle composites with soft inclusions. SIAM J Math Anal 48: 1178-1209. doi: 10.1137/15M1007343
    [7] Bouchitté G, Fonseca I, Leoni G, et al. (2002) A global method for relaxation in W1,p and in SBVp. Arch Ration Mech Anal 165: 187-242. doi: 10.1007/s00205-002-0220-y
    [8] Braides A, Defranceschi A (1998) Homogenization of Multiple Integrals. New York: Oxford University Press.
    [9] Braides A, Defranceschi A, Vitali E (1996) Homogenization of free discontinuity problems. Arch Ration Mech Anal 135: 297-356. doi: 10.1007/BF02198476
    [10] Braides A, Garroni A (1995) Homogenization of periodic nonlinear media with stiff and soft inclusions. Math Mod Meth Appl Sci 5: 543-564. doi: 10.1142/S0218202595000322
    [11] Braides A, Solci M (2013) Multi-scale free-discontinuity problems with soft inclusions. Boll Unione Mat Ital 1: 29-51.
    [12] Cagnetti F, Dal Maso G, Scardia L, et al. (2019) Γ-convergence of free-discontinuity problems. Ann Inst H Poincaré Anal Non Linéaire 36: 1035-1079. doi: 10.1016/j.anihpc.2018.11.003
    [13] Cagnetti F, Dal Maso G, Scardia L, et al. (2019) Homogenisation of stochastic free-discontinuity problems. Arch Ration Mech Anal 233: 935-974. doi: 10.1007/s00205-019-01372-x
    [14] Cagnetti F, Scardia L (2011) An extension theorem in SBV and an application to the homogenization of the Mumford-Shah functional in perforated domains. J Math Pur Appl 95: 349-381. doi: 10.1016/j.matpur.2010.03.002
    [15] Dal Maso G (1993) An Introduction to Γ-Convergence. Boston: Birkhäuser.
    [16] Dal Maso G, Zeppieri CI (2010) Homogenization of fiber reinforced brittle materials: The intermediate case. Adv Calc Var 3: 345-370.
    [17] Focardi M, Gelli MS, Ponsiglione M (2009) Fracture mechanics in perforated domains: A variational model for brittle porous media. Math Mod Meth Appl Sci 19: 2065-2100. doi: 10.1142/S0218202509004042
    [18] Giacomini A, Ponsiglione M (2006) A Γ-convergence approach to stability of unilateral minimality properties. Arch Ration Mech Anal 180: 399-447. doi: 10.1007/s00205-005-0392-3
    [19] Pellet X, Scardia L, Zeppieri CI (2019) Homogenization of high-contrast Mumford-Shah energies. SIAM J Math Anal 51: 1696-1729. doi: 10.1137/18M1189804
    [20] Scardia L (2008) Damage as Γ-limit of microfractures in anti-plane linearized elasticity. Math Mod Meth Appl Sci 18: 1703-1740. doi: 10.1142/S0218202508003170
    [21] Scardia L (2010) Damage as the Γ-limit of microfractures in linearized elasticity under the non-interpenetration constraint. Adv Calc Var 3: 423-458.
  • This article has been cited by:

    1. L. De Luca, M. Ponsiglione, Variational models in elasticity, 2021, 3, 2640-3501, 1, 10.3934/mine.2021015
    2. A.G. Kolpakov, S.I. Rakin, Local stresses in the reinforced plate with orthogonal systems of fibers. The homogenized strength criterion, 2023, 275, 00457949, 106922, 10.1016/j.compstruc.2022.106922
    3. Xavier Pellet, Lucia Scardia, Caterina Ida Zeppieri, Stochastic homogenisation of free-discontinuity functionals in randomly perforated domains, 2024, 17, 1864-8258, 643, 10.1515/acv-2022-0052
  • Reader Comments
  • © 2020 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(3492) PDF downloads(375) Cited by(3)

Figures and Tables

Figures(2)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog