arxiv: v3 [physics.class-ph] 18 Dec 2017

Transcription

1 Micromechanics based framework with second-order damage tensors R. Desmorat a, B. Desmorat b,c, M. Olive a, B. Kolev d a LMT-Cachan (ENS Cachan, CNRS, Université Paris Saclay), F-9425 Cachan Cedex, France b Sorbonne Université, UMPC Univ Paris 06, CNRS, UMR 7190, Institut d Alembert,F Paris Cedex 05, France c Université Paris SUD 11, Orsay, France d Aix Marseille Université, CNRS, Centrale Marseille, I2M, UMR 77, 145 Marseille, France arxiv: v [physics.class-ph] 18 Dec 2017 Abstract The harmonic product of tensors leading to the concept of harmonic factorization has been defined in a previous work (Olive et al, 2017). In the practical case of D crack density measurements on thin or thick walled structures, this mathematical tool allows us to factorize the harmonic (irreducible) part of the fourth-order damage tensor as an harmonic square: an exact harmonic square in 2D, an harmonic square over the set of so-called mechanically accessible directions for measurements in the D case. The corresponding micro-mechanics framework based on second instead of fourth order damage tensors is derived. An illustrating example is provided showing how the proposed framework allows for the modeling of the so-called hydrostatic sensitivity up to high damage levels. Keywords: anisotropic damage, crack density, harmonic decomposition PACS: a, x, Ba 1. Introduction The damage anisotropy encountered in quasi-brittle materials is induced by the loading direction and multiaxiality. From a micro-mechanics point of view, it is the consequence of an oriented microcracking pattern. From the Continuum Damage Mechanics (CDM) point of view, the anisotropic damage state is represented by a tensorial thermodynamics variable, either an eightorder tensor (Chaboche, 1978, 1979), a fourth-order damage tensor D (Chaboche, 1978, 1979; Leckie and Onat, 1980; Chaboche, 1984; Lemaitre and Chaboche, 1985; Andrieux et al., 1986; Ju, 1989; Kachanov, 199; Zheng and Collins, 1998; Cormery and Welemane, 2010; Dormieux and Kondo, 2016) or a symmetric second-order damage tensor d (Vakulenko and Kachanov, 1971; Murakami and Ohno, 1978; Cordebois and Sidoroff, 1982; Ladevèze, 198; Murakami, 1988). There exist many second-order anisotropic damage frameworks (Murakami, 1988; Kattan and Voyiadjis, 1990; Ramtani et al., 1992; Papa and Taliercio, 1996; Halm and Dragon, 1998; Steinmann and Carol, 1998; Lemaitre et al., 2000; Carol et al., 2001; Menzel and Steinmann, 2001; Menzel et al., 2002; Brunig, 200; Lemaitre and Desmorat, 2005; Desmorat et al., 2007; Badel et al., 2007; Desmorat and Otin, 2008; Desmorat, 2016), as their unification into a single model is partial (Ladevèze, 198, 1995). A link with the theory of second order fabric tensors has been made in (Zysset and Curnier, 1995; Voyiadjis and Kattan, 2006). From a theoretical point of view (Leckie and Onat, 1980; Onat, 1984), second order damage frameworks are usually seen to be restrictive compared to the fourth-order tensorial one. Nevertheless, the interpretation of a damage variable being simpler when a second-order tensor is considered (the three principal values d i of d naturally correspond to orthogonal families of microcracks), less damage parameters are introduced and the second-order frameworks have been widely used for either ductile or quasi-brittle materials. The recent analysis of 2D cracked media with both open and closed microcraks has shown that the so-called irreducible (harmonic) part H 2D of the damage tensor can be decomposed by means of a second-order damage tensor (Desmorat and Desmorat, 2016). More precisely, the standard second-order crack density tensor of Vakulenko and Kachanov (1971) still represents the open cracks contribution when a novel (deviatoric) second-order damage tensor represents the closed sliding cracks (previously represented by a fourth-order tensor, Andrieux et al. (1986); Kachanov (199)). This can be achieved using Verchery s polar decomposition of 2D fourth-order tensors (Verchery, 1979; Vannucci, 2005), which includes both (Desmorat and Desmorat, 2015): addresses: desmorat@lmt.ens-cachan.fr (R. Desmorat), boris.desmorat@upmc.fr (B. Desmorat), marc.olive@math.cnrs.fr (M. Olive), boris.kolev@math.cnrs.fr (B. Kolev) the harmonic decomposition of considered tensor; the harmonic factorization of its fourth-order irre- Preprint submitted to European Journal of Mechanics A/Solids (accepted, 24, 2018

2 ducible (harmonic) part H 2D, by means of a deviatoric second-order tensor h 2D : H 2D = h 2D * h 2D. The harmonic product between harmonic tensors, written as h 1 *h 2, is defined as the projection of the (totally) symmetric tensor product h 1 h 2 onto the space of harmonic tensors (see Sections and 4.4). The question arises then as how to extend these results in D? We know from (Olive et al., 2017) that any D harmonic fourth-order tensor can be factorized into H = h 1 * h 2, i.e., represented by two (deviatoric) second-order tensors h 1, h 2. However, the factorization is far from being unique. To overpass these difficulties, we point out that triaxial mechanical testing is of high complexity, both from the experimental set-up needed (a triaxial machine) and from the difficulty to measure mechanical properties in different space directions (Calloch, 1997; Calloch and Marquis, 1999). We propose, here, to restrict ourselves to a simpler, but still sufficiently general, situation: the case of measurements of a D crack density function on structures. Well-known cases are the thin walled structures, such as plates, tubes and shells for which the thinner direction is the normal ν. But the present work also applies to D thick structures (as the cube of Section 8), as long as an out-of-plane normal can be locally defined. Instead of considering the representation of crack density in any direction n, we shall then consider, in Section 5, its representation to a restricted set of directions R(ν) := {τ; τ = 1 and τ ν = 0} {ν}, i.e., the in-plane directions τ, orthogonal to ν, and the out-of-plane direction ν, normal to the structure itself. In the present work, we consider these directions as the mechanically accessible directions for measurements. After recalling the required mathematical tools (harmonic decomposition, harmonic product and Sylvester s theorem), we revisit the link between crack density function and the tensorial nature of the damage variables. This will allow us to derive a general micro-mechanics based D framework with second instead of fourth order damage tensors. Definitions We denote by T s the totally symmetric part of a possibly non symmetric tensor T. More precisely T s (x 1,..., x n ) := 1 T(x σ(1),..., x σ(n) ), n! σ S n where S n is the permutation group on the indices {1,..., n}. The symmetric tensor product of two tensors 2 T 1 and T 2, of respective orders n 1 and n 2, is the symmetrization of T 1 T 2, defining a totally symmetric tensor of order n = n 1 + n 2 : T 1 T 2 := (T 1 T 2 ) s. Contracting two indices i, j of a tensor T of order n defines a new tensor of order n 2 denoted as tr ij T. For a totally symmetric tensor T, this operation does not depend on a particular choice of the pair i, j. Thus, we can refer to this contraction just as the trace of T and we will denote it as tr T. It is a totally symmetric tensor of order n 2. Iterating the process, we define tr k T = tr(tr( (tr T))), which is a totally symmetric tensor of order n 2k. In D, a totally symmetric fourth-order tensor T has no more than 15 independent components, instead of 21 for a triclinic elasticity tensor (i.e. a tensor E having minor symmetry E ijkl = E jikl = E ijlk and major symmetry E ijkl = E klij ). Totally symmetric elasticity tensors were called rari-constant in the nineteenth century (Navier, 1827; Cauchy, 1828b,a; Poisson, 1829; Love, 1905; Vannucci and Desmorat, 2016). 2. Harmonic decomposition The harmonic decomposition of tensors (Schouten, 1954; Spencer, 1970), introduced in geophysics by Backus (1970), has been popularized by Leckie and Onat (1980) and Onat (1984) when deriving fourth-order damage tensor and later by Forte and Vianello (1996) when classifying elasticity symmetries Harmonic tensors and corresponding polynomials An harmonic tensor is a traceless, totally symmetric tensor, i.e. H = H s, and tr H = 0. To every totally symmetric tensor H of order n, with components H i1i 2 i n, corresponds a unique homogenous polynomial (and conversely). More precisely, h(x) := H(x, x,..., x) = H i1i 2 i n x i1 x i2 x in is a homogeneous polynomial h(x) = h(x 1, x 2, x ), of degree n in the spacial coordinates x 1, x 2, x. It is harmonic since 2 h = 0, due to the traceless property tr H = 0.

3 2.2. Harmonic decomposition of a symmetric tensor Any totally symmetric tensor T of order n can be decomposed uniquely as T = H H r 1 H r r H r (2.1) where r = [n/2] is the integer part of n/2, H k is an harmonic tensor of degree n 2k and 1 k = 1 1 means the symmetrized tensorial product of k copies of the (second-order) identity tensor. For n even (n = 2r), one has: H r = H n = 1 2 n + 1 tr n 2 T, (2.2) where H r = H r is a scalar in that case. Moreover, H r 1,..., H 0 are obtained inductively (Olive et al., 2017) as follows: H k = μ(k, n) tr k [ T where μ(k, n) = r j=k+1 1 j H j ] (2n 4k + 1)!(n k)!n! (2n 2k + 1)!k!(n 2k)!(n 2k)!. (2.) Remark 2.1. It is worth emphasizing the fact that this harmonic decomposition is just a generalization to higher order symmetric tensors of the well-known decomposition of a symmetric second-order tensor into its deviatoric/spheric parts: d = d + 1 (tr d)1. Decomposition (2.1) is an orthogonal decomposition (relative to the natural Euclidean product on the space of symmetric tensors). The projection H 0 onto the space of highest order harmonic tensors (same order as T) will be called the harmonic part of T and denoted by (T) 0 : (T) 0 := H 0 = T 1 H 1 1 r H r. (2.4) 2.. Harmonic decomposition of the elasticity tensor The harmonic decomposition of an elasticity tensor E, a fourth-order tensor having both minor and major symmetries (E ijkl = E jikl = E ijlk and E ijkl = E klij ), was first obtained by Backus (1970), as: E = α 1 (4) 1+β 1 (2,2) 1+1 (4) a +1 (2,2) b +H (2.5) where ( ) = ( ) 1 tr( ) 1 denotes the deviatoric part of a second-order tensor. In formula (2.5), the Young-symmetrized tensor products (4) and (2,2), between two symmetric second-order tensors y, z, are defined as follows: { ) y (4) z = 6( 1 y z + z y + 2 y z + 2 z y, ) y (2,2) z = ( 1 y z + z y y z z y, (2.6) where (y z) ijkl = 1 2 (y ikz jl + y il z jk ) so that (4) is the same as the totally symmetric tensor product : y (4) z = y z. In the harmonic decomposition (2.5), H is a fourth-order harmonic tensor, α, β are scalars, and a, b are secondorder harmonic tensors (deviators) related to the dilatation tensor di = tr 12 E and the Voigt tensor vo = tr 1 E by the formulas: and α = 1 15 (tr di + 2 tr vo), β = 1 (tr di tr vo), (2.7) 6 a = 2 7 ( di + 2vo ), b = 2 ( di vo ). (2.8) The harmonic part of E is defined as: or similarly as: (E) 0 := H = E 1 (4) a 1 (2,2) b (2.9) (E) 0 := E 1 a 1 (1 b + b 1 1 b b 1), where a = a + α1 and b = b + β1. The scalars α, β and the second-order deviators a, b are given by (2.7) and (2.8).. The harmonic product and Sylvester s theorem The harmonic product of two harmonic tensors of order n 1 and n 2, defining an harmonic tensor of order n = n 1 + n 2, has been introduced in (Olive et al., 2017) as the harmonic part of the symmetric tensor product: H 1 * H 2 := (H 1 H 2 ) 0. Note that this product is associative: and commutative: H 1 * (H 2 * H ) = (H 1 * H 2 ) * H, H 1 * H 2 = H 2 * H 1. For two vectors w 1, w 2, we have w 1 * w 2 =(w 1 w 2 ) = 1 2 (w 1 w 2 + w 2 w 1 ) 1 (w 1 w 2 ) 1, where w 1 w 2 = w T 1 w 2 is the scalar product. For two second-order harmonic tensors (deviators) h 1, h 2, we have h 1 * h 2 = h 1 h (h 1h 2 + h 2 h 1 ) tr(h 1h 2 ) 1 1. (.1)

4 Sylvester s theorem (Sylvester, 1909; Backus, 1970; Baerheim, 1998) states that any harmonic tensor H of order n can be factorized as H = w 1 * w 2 * * w n, i.e. as the harmonic products of n (real) vectors w k, the so-called Sylvester-Maxwell multipoles. Note however, that this factorization is far from being unique, as discussed in (Olive et al., 2017). Setting h 1 = w 1 * w 2 and h 2 = w * w 4 which are harmonic second-order tensors (deviators), we obtain the non unique harmonic factorization of H by means of two second-order tensors: H = h 1 * h 2, (.2) as detailed in (Desmorat and Desmorat, 2016; Olive et al., 2017). 4. Link between fourth-order crack density and damage tensors Before formulating our main result, Theorem 5.1, we summarize, in this section, the present state of the art in Continuum Mechanics leading to the representation of damage of cracked media by a fourth-order tensor (Chaboche, 1979). We make an explicit link with the harmonic decomposition and we present, by comparison to the 2D case, the problem of representation of damage by second-order tensors in D Crack density function and tensors The damage state of a microcracked material is classically defined by spatial arrangement, orientation and geometry of the cracks present at the microscale (Kachanov, 1972; Leckie and Onat, 1980; Ladevèze, 198; Onat, 1984; Lemaitre and Chaboche, 1985; Murakami, 1988; Kachanov, 199). The crack density, related to any possible D direction defined by a unit vector n, refers to a dimensionless scalar property defined in a continuous manner at the Representative Volume Element scale as a spatial crack density function Ω = Ω(n). Owing to the property Ω(n) = Ω( n), it is expressed by means of a totally symmetric tensor F (the so-called fabric tensor) of even order n = 2r (Kanatani, 1984) as: Ω(n) = F (n n n) (4.1) where means the contraction over the n subscripts. Note that Ω(n) corresponds to a homogeneous polynomial (see section 2.2) h(n 1, n 2, n ) = F(n, n,..., n). The fabric tensor F, which is totally symmetric, can be determined as the least square error approximation of an 4 experimental (measured) density distribution Ω(n), F being thus solution of 2 4π min F 2n + 1 F 1 n Ω(n) n n ds, with solid angle S and where x =1 n k := n n n. Moreover, the following equality has been used: 1 n 2n 1 ds = 4π 2n n. x =1 Note that n k = n k is a totally symmetric tensor. Comparative studies of the tensorial order, needed to represent given microcracking patterns, can be found in (Lubarda and Krajcinovic, 199; Krajcinovic, 1996; Tikhomirov et al., 2001). Expression (4.1) is often rewritten into the finite expansion (Kanatani, 1984; Onat, 1984; Krajcinovic, 1996): Ω(n) = F 4 (n n n n) + Ω 6 (n 6 ) + with fourth-order part + Ω 2k (n 2k ) + + Ω n (n n ) (4.2) F 4 (n n n n) =Ω 0 + Ω 2 (n n) + Ω 4 (n n n n) (4.) with n = 2r even and where Ω 2k are totally symmetric traceless (harmonic) tensors of order 2k. The scalar term Ω 0 is the crack density within considered Continuum Mechanics representative volume element Ω 0 = 1 4π x =1 Ω(n) ds. Crack density tensors Ω 2, Ω 4,..., Ω n are harmonic tensors of even order 2, 4,..., n. They constitute independent crack density variables representative of the microcraking pattern (and anisotropy), determined uniquely up to order n from the knowledge of the D spatial crack density distribution Ω(n) Derivation of the crack density tensors from the harmonic decomposition Let us point out that the harmonic tensors Ω 2k correspond to the tensors H r k issued from the harmonic decomposition (2.1) of the fabric tensor F: F = H H r 1 H r r H r, with r = n/2, where H r = H n and the harmonic tensors 2 H k of degree n 2k are given by (2.2) and (2.). Observe, moreover, that: (1 k H k ) n n = H k n n 2k,

5 and we get thus: Ω(n) = H r + H r 1 (n n) + H r 2 (n n n n) + + H 0 (n n ), which is the finite expansion (4.2), where Ω 0 = H n 2, Ω 2k = H r k. 4.. Fourth-order damage tensor Using the decomposition (4.2) and assuming open microcracks in an initially D isotropic medium, Leckie and Onat (1980) and Onat (1984) have shown that the damage variable defined by the coupling microcraking/elasticity is at most a fourth-order tensor, built from F 4 only, see (4.2). This result holds for non interacting closed sliding without friction pennyshaped microcracks (Kachanov, 199) and, as pointed out by Cormery and Welemane (2010), for many stress based homogenization schemes, as long as all the microcracks are in the same state, either open or closed. Setting: J = I 1 1 1, the following general definition of a fourth-order damage tensor has then been derived for initially isotropic materials (Kachanov, 199; Zheng and Collins, 1998; Cormery and Welemane, 2010): D = p 0 Ω p 1 Ω 0 J + p 2 (1 Ω 2 + Ω 2 1) + p (1 Ω 2 + Ω 2 1) + p 4 Ω 4, (4.4) where Ω 0 should be interpreted as a scalar damage variable, the symmetric deviator Ω 2 as a second-order damage variable, and the harmonic tensor Ω 4 as a fourth-order damage variable. The expression of the scalars p i depends on the initial elasticity parameters, on the homogenization scheme and on the microcracks state (simultaneously open or simultaneously closed for all cracks). Remark 4.1. The scalars p i do not depend on Ω 0, Ω 2, Ω 4. Remark 4.2. (4.4) is the harmonic decomposition (2.5) of the fourth-order damage tensor D, which has the major and the minor indicial symmetries (D ijkl = D klij = D jikl ) as an elasticity tensor. The deviatoric parts of the dilatation and Voigt tensors are both proportional to the secondorder harmonic tensor Ω 2, with the scalar factors κ di and κ vo depending only on the initial elastic parameters of the undamaged isotropic material: di (D) = (tr 12 D) = κ di Ω 2, vo (D) = (tr 1 D) = κ vo Ω 2. The traces of the dilatation and the Voigt tensors are both proportional to the scalar crack density Ω 0, with scalar factors k di and k vo depending only on the elastic parameters of virgin (undamaged) isotropic material: tr di(d) = tr(tr 12 D) = k di Ω 0, tr vo(d) = tr(tr 1 D) = k vo Ω 0. 5 Remark 4.. An alternative framework is due to Voyiadjis and Kattan (2006). These authors extend to anisotropic damage the framework of Zysset and Curnier (1995) for the representation of microstructure morphology of granular materials (the considered framework neglects the fourth order contribution Ω 4 ). They propose, then, a nonlinear link between the second order fabric tensor Ω 2 and a fourth order tensorial damage variable D, setting for the effective (damaged) elasticity tensor Ẽ = Ẽ = λ φ φ + 2μ φ φ, φ = (Ω Ω 2 ) k (4.5) This means that they replace the identity tensor 1 in the usual isotropic elasticity law by a negative power k = 0.2 of the second order crack density tensor Ω Ω 2 (powers being taken in terms of the principal values, λ, μ considered as Lamé-like constants). The crack density Ω 0 must have a non zero initial value. Due to the presence of quadratic terms in crack densities, this framework does not satisfy previous proportionality properties. Fourth order damage tensor D is such that Ẽ = (1 D) : E D case In 2D, cracks are represented by 2D straight lines. Expression (4.2) for crack density holds, recovering a Fourier finite expansion (Kanatani, 1984; Burr et al., 1995): Ω(n) = ω 2D + ω 2D (n n) + H 2D (n n n n) +, (4.6) where the unit vector n is related to the possible planar direction ω 2D = 1 2π Ω(n) dθ 2π is the 2D crack density, and where ω 2D, H 2D are respectively the 2D harmonic second and the fourth-order crack density tensors. Verchery s decomposition (Verchery, 1979; Vannucci, 2005), and its rewriting into a tensorial form (Desmorat and Desmorat, 2015), shows that any 2D harmonic fourthorder tensor is an harmonic square. Applied to H 2D, this gives: H 2D = h 2D * h 2D, (4.7) where h 2D is an harmonic second-order tensor (deviator), and the notation * denotes the harmonic product defined as the orthogonal projection of the symmetrized product h 2D h 2D onto 2D fourth-order harmonic tensors space: h 2D * h 2D := (h 2D h 2D ) 0. For second-order harmonic tensors, this reads: h 2D * h 2D = h 2D h 2D 1 4 (tr h2 2D) 1 1. This means that in 2D, any anisotropic microcracking pattern can be expressed, up to order 4, exactly by means 0

6 of the scalar ω 2D and the two independent second-order deviatoric damage variables ω = Ω 2D 2 and h 2D = h 2D. This result is consistent with the fact that the micro-mechanics of 2D media with open and closed (sliding without friction) microcracks can be represented by two second-order damage tensors only (Desmorat and Desmorat, 2015). The question arises then whether the expansion (4.6) (4.7) holds in D, i.e. with ω and h (now D) deviatoric second-order tensors. The answer is negative in the general triclinic case for which we have only Ω 4 = h 1 * h 2 (.2) with usually different second-order tensors h 1, h 2 (Olive et al., 2017). Furthermore, the factorization is not unique, forbidding to interpret h 1 and h 2 as damage variables. 5. D second-order damage tensors from walled structures It was noticed by Lubarda and Krajcinovic (199) and Krajcinovic (1996) that the fourth-order crack density tensor Ω 4 = (F 4 ) 0 (the harmonic part of the totally symmetric tensor F 4 ) induced by particular loadings is related, sometimes, as a square of the second-order harmonic contribution Ω 2. More precisely, Ω 4 is proportional to the harmonic square Ω 2 * Ω 2 in the particular situation that occurs for a family of parallel penny shaped microcracks having identical therefore coplanar normal m (induced for example by uniaxial tension on quasi-brittle materials). An harmonic square is also present in the work of Voyiadjis and Kattan (2006). Let us show that it corresponds to the fourth order harmonic tensor φ * φ with φ = (Ω Ω 2 ) 0.2 : The harmonic part (Ẽ) 0 of the effective elasticity tensor (4.5) is (Ẽ) 0 = (λ (φ φ) s + 2μ (φ φ) s ) 0, (5.1) where (.) s means the totally symmetric part and (.) 0 is the harmonic projection defined in Eq. (2.9). By definition (φ φ) s = (φ φ) s = φ φ. The harmonic fourth order part of Ẽ is thus (Ẽ) 0 = (λ + 2μ)(φ φ) 0 = (λ + 2μ) φ * φ, where * is the harmonic product (.1). See also Appendix B. such as plates. This allows us to introduce the unit normal ν ( ν = 1) to the walled structure and the set R(ν) := {τ, τ = 1 and τ ν = 0} {ν}, of directions n restricted to so-called mechanically accessible directions for measurements. As an extension of both the previous remark on fourthorder harmonic squares and the 2D result (4.6) (4.7), we have in D the following theorem (the proof of which is given at the end of the present section): Theorem 5.1. For a given unit vector ν, any density function Ω(n) is represented, up to fourth-order, for all directions n R(ν) by means of a scalar ω m and two harmonic (symmetric deviatoric) second-order tensors ω and h as: Ω(n) = ω m + ω (n n) + (h * h) (n n n n) + (5.2) for all n R(ν). This representation is unique, up to ±h, if (ω ν) ν = h ν = 0. Remark 5.2. If we set e = ν, the conditions (ω ν) ν = 0 (which is equivalent to (ω ν) = λ ν) and h ν = 0 mean that ω ω 11 ω 12 0 = ω 12 ω 22 0, (5.) 0 0 (ω 11 + ω 22) and h = h 11 h 12 0 h 12 h (5.4) Applied to thinned walled structured, for which n R(ν) are the accessible directions for mechanical measurements, Theorem 5.1 states that any microcracking pattern, possibly triclinic, can be represented, up to order 4, by means of only two symmetric second-order crack density tensors, ω = ω + ω m 1 and h, the second one being a deviator. In that case, we can recast (5.2) as: Ω(n) = ω (n n) + (h * h) (n n n n) + for all n R(ν), where ω is the crack density tensor introduced by Vakulenko and Kachanov (1971) (Kachanov, 1972, 199). The second and fourth-order crack density variables Ω 2 and Ω 4 are independent in general case. We prefer to keep the fourth order contribution Ω 4 in our further analyses and we will use next the standard expression (4.4) (for instance instead of Eq. (4.5) which neglects Ω 4 ). In order to built a micromechanics based framework with second-order damage tensors, one considers as the general case except for soils that the measurements of D crack density Ω(n) is performed on thin or thick walled structures, i.e. on thin or thick tubes or on 2D structures 6 Practical formulas Proof of Theorem 5.1 Recall that, up to order four, the crack density function Ω(n) is represented by the fabric tensor F 4 (4.2) (4.). Consider now an orthonormal frame {e 1, e 2, ν} and let (ω m, ω, h) be a triplet as in (5.) (5.4). Set h 11 + ih 12 = 1 2[ (F4 ) (F 4 ) (F 4 ) i ( (F 4 ) 1112 (F 4 ) 1222 ) ] 1/2 (5.5)

7 and ω m = 1 4 (F 4) (F 4) (F 4) (F 4) 1 15 (h h 2 12), ω 11 = 5 8 (F 4) (F 4) (F 4) (F 4) (h h 2 12), ω 22 = 8 (F 4) (F 4) (F 4) (F 4) (h h 2 12), ω 12 =(F 4 ) (F 4 ) 1112, where i = 1 is the pure imaginary number. It can be checked by a direct computation that (ω m, ω, h) is a solution of (5.2) in Theorem 5.1. Remark 5.. Because of the square root in (5.5), both h and h are solutions. We will now show the uniqueness of the solution, up to a sign, and under the assumption that: (ω ν) ν = h ν = 0. Suppose thus, that a second solution (ω m, * ω *, h * ) to (5.2) exists, with ω * and h * as in (5.) (5.4). Equaling, for different directions n R(ν), the density function Ω(n) defined by (5.2), calculated first with (ω m, ω, h) and then with (ω m, * ω *, h * ), we obtain for n = ν: 5(ω * 11 + ω * 22 ω 11 ω 22 + ω m ω * m) + 4h h h *2 11 4h *2 12 = 0. (5.6) Then, for n = (cos θ, sin θ, 0), we get: a 0 + a 2 cos 2θ + b 2 sin 2θ + a 4 cos 4θ + b 4 sin 4θ = 0, (5.7) where a 0 = 2(ω m ω m) * + ω 11 ω 11 * + ω 22 ω 22 * + 5 (h2 11 h * h 2 12 h *2 12), a 2 = ω 11 ω 11 * ω 22 + ω 22, * b 2 = 2(ω 12 ω 12), * a 4 = h 2 11 h 2 12 h * h *2 12, b 4 = 2(h 11 h 12 h * 11h * 12). Since (5.7) holds for all θ, we have a 0 = a 2 = b 2 = a 4 = b 4 = 0. Since a 4 = b 4 = 0, we get { h *2 11 h *2 12 =h 2 11 h 2 12, h * 11h * 12 =h 11 h 12, 7 from which we deduce that (h * 11 + ih * 12) 2 = (h 11 + ih 12 ) 2 and therefore that h * = ±h (in accordance with Remark 5.). From (5.6) and a 0 = a 2 = b 2 = 0, we get ω m ω m * + ω 11 * ω 11 + ω 22 * ω 22 = 0, 2(ω m ω m) * + ω 11 ω 11 * + ω 22 ω 22 * = 0, ω 11 ω 11 * ω 22 + ω 22 * = 0, ω 12 ω 12 * = 0, i.e. ω * m = ω m and ω * = ω, which achieves the proof. 6. General micro-mechanics based framework with two second-order damage variables Using the results from Section 4., we deduce from (5.2) that the representation by means of two symmetric secondorder tensors holds for the damage tensor itself, at least when the microcracks are all in the same state, all open or all closed. This means that, disposing from sufficiently many in-plane measurements (along directions n orthogonal to ν) and an out-of-plane measurement (along n = ν), the general fourth-order damage tensor of Chaboche Leckie Onat can be expressed by means of two symmetric second-order damage variables only, for example ω and h (the second one being a deviator). A general damage framework using this feature is derived next, clarifying the link between Cordebois and Sidoroff (1982) and Ladevèze (198, 1995) phenomenological second-order damage models and micro-mechanics based framework. We shall assume that the homogenization result (4.4) holds, where the constants p i are given (refer to the works of Kachanov (199) and Dormieux and Kondo (2016) for comparison of different homogenization schemes). Gibbs free enthalpy density writes ρψ = 1 18K (tr σ) G σ : σ + 1 σ : D : σ, 2E E E (1 2ν) where ρ is the density and E, G = 2(1+ν), K = are, respectively, the Young, shear and bulk moduli. The elasticity law, coupled with the anisotropic damage, writes then as ε e = ρ ψ σ = 1 2G σ + 1 9K (tr σ) E D : σ, or, in a more compact form, as ε e = S : σ, where ε e is the elastic strain tensor and S, the effective fourth-order compliance tensor S = 1 9K G J + 1 E D, J = I (6.1) Having many in-plane and possibly one out-of-plane measurements, allows us to use remark 4.1 and (5.2)

8 instead of (4.4), within the considered homogenization scheme. We can thus recast the fourth-order damage tensor D by substituting the scalar Ω 0 by ω m, the secondorder tensor Ω 2 by the deviatoric tensor ω and the fourthorder tensor Ω 4 by the harmonic (i.e totally symmetric and traceless) tensor h * h. More precisely, we get D = p 0 ω m p 1 ω m J + p 2 (1 ω + ω 1) + p (1 ω + ω 1) + p 4 h * h. (6.2) Using (.1), the term h * h expands as h * h = 1 h h + 2 h h 2 21 ( 1 h 2 + h (1 h 2 + h 2 1) ) (tr h2 ) ( ), so that the enthalpic contribution, due to the microcracks, writes σ : D : σ = p 0 ω m (tr σ) 2 + p 1 ω m tr(σ 2 ) + 2p 2 tr(ω σ) tr σ + p tr(ω σ 2 ) + p 4 [ 1 (tr(hσ )) tr(σ hσ h) 8 21 tr(h2 σ 2 ) tr h2 tr σ 2]. Again, as in Remark 4.2, (6.2) is nothing else but the harmonic decomposition of the fourth-order damage tensor D, but with the following particularities. Let di(d) := (tr 12 D) be the dilatation tensor of D and vo(d) := (tr 1 D) be the Voigt tensor of D. Then: the harmonic part H = h * h of D is factorized as an harmonic square; the following proportionality relations hold between the deviatoric parts of di(d) and vo(d): 7. A second-order anisotropic damage model in micro-mechanics based framework Following Cordebois and Sidoroff (1982) and Ladevèze (198), a symmetric second-order, unbounded damage variable Φ is introduced in the Gibbs free enthalpy (with initial value Φ = 1 for a virgin material, and with damage growth d dtφ positive definite). The usual second-order damage tensor writes as: d = 1 Φ 2 (with initial value d = 0). A general but phenomenological coupling of elasticity with second-order anisotropic damage is described in (Desmorat, 2006). It reads ρψ = g(φ) 18K (tr σ) G tr(φ σ Φ σ ), (7.1) where σ = σ 1 (tr σ)1 is the stress deviator. The function g was chosen as g(φ) := 1 1 tr d = 1 tr Φ 2 2 for metals in (Lemaitre et al., 2000), and as g(φ) := 1 tr Φ2 (7.2) for concrete in (Desmorat, 2016). In both models, the convexity with respect to σ and the positivity of the intrinsic dissipation are satisfied (see also Chambart et al. (2014)). The elasticity law writes ε e = ρ ψ σ = g(φ) 1 (tr σ)1 + 9K 2G (Φ σ Φ) (7.) with effective compliance tensor S = 1 9K G J + 1 D, (7.4) E where 1 E D has for harmonic decomposition (proof given in Appendix B) di (D) vo (D) ω, (6.) which is equivalent for the effective compliance tensor S to satisfy the same conditions: di ( S) vo ( S) ω ; (6.4) the following proportionality relations hold between the traces of di(d) and vo(d): tr di(d) tr vo(d) ω m. (6.5) Following Cormery and Welemane (2010), who consider the scalar constants p i as material parameters, conditions 1 to above, are the conditions for a damage model for instance built in a phenomenological manner which should be considered as micro-mechanics based. 8 1 E D = g(φ) 1 9K β 1 2G J 2 G (1 b + b 1) with + 1 G (1 b + b 1) + 1 2G Φ * Φ (7.5) β = 1 [ tr(φ 2 ) + (tr Φ) 2], b = 1 [ (tr Φ)Φ (Φ 2 ) ] The harmonic part of the fourth order damage tensor is the harmonic square, (D) 0 = (1 + ν) Φ * Φ, (7.6)

9 where ν is Poisson ratio of undamaged material. It satisfies the first condition on the effective compliance S to be of micro-mechanics based form (6.2), with Moreover, we have (tr 12 D) = 0, p 4 = 1 + ν, h = Φ. (tr 1 D) = 7(1 + ν) b. Both deviatoric parts (tr 12 D) and (tr 1 D) are obviously proportional, so that the phenomenological anisotropic damage model satisfies the second condition on the effective compliance S to be of micro-mechanics based form (6.2), with p 2 ω 2(1 + ν) ( = (Φ 2 ) (tr Φ) Φ ), 21 p ω = 2 p 2ω. By identification of the isotropic part of the harmonic decomposition of D, we get p 0 ω m = 1 2ν p 1 ω m =(1 + ν) (g(φ) 1), ( 1 10 (tr Φ) tr Φ2 1 ). Recall that the material constants p 0, p 1 are independent of Φ, they are considered as material parameters, so that the proportionality requirement p 0 ω m p 1 ω m can be satisfied if following Burr et al. (1995) and Lemaitre et al. (2000) we define the hydrostatic sensitivity parameter η as the material constant and set g(φ) = (1 η) + η η = p 0(1 + ν) p 1 (1 2ν), ( 1 10 (tr Φ)2 + 1 ) tr Φ2. (7.7) 0 Condition for the model (7.1) to be considered as micromechanics based is then fulfilled. For η 0, the Gibbs free enthalpy density (7.1) is furthermore convex with respect to both the stress tensor σ and the damage tensor Φ. Note that a full anisotropic damage model including damage evolution for quasi brittle materials can be naturally derived following (Desmorat, 2016). The hydrostatic sensitivity just obtained in Eq. (7.7) is quantified in next Section. 8. Hydrostatic sensitivity Let us consider quasi-brittle materials such as concrete and the micro-mechanically based damage model just derived (Section 7, Gibbs free enthalpy density (7.1), elasticity law (7.) and function g(φ), given by Eq. (7.7)). The 9 second order damage variable, which vanishes for virgin material, is d = 1 Φ 2. Even when the damage state is anisotropic, the constitutive equations considered in Section 7 is such that the dilatation tensor of the effective compliance tensor is null. This allows for the definition of an effective bulk modulus K for damaged materials (i.e. a modulus function of damage tensor and of material parameters of the undamaged material) such as tr σ = K tr ε e Calculating the trace of elasticity law (7.), we get K = K g(φ) = K (1 η) + η ( 1 10 (tr Φ)2 + 1 (8.1) 0 tr Φ2), E (1 2ν) where K = is the bulk modulus for the virgin (undamaged) material. Recall that g(1) = 1. The measurement and therefore the identification of the material parameter η are not easy tasks for quasi-brittle materials. This is why Discrete Elements (Cundall and Strack, 1979; Herrmann and Roux, 1990; Schlangen and Garboczi, 1997; Van Mier et al, 2002; Olivier-Leblond et al, 2015) are often used as numerical experimentation for the tensile states of stresses of quasi-brittle materials. In such numerical tests the material is described as a particles assembly representative of the material heterogeneity, the particles being here linked by elastic-brittle beams. The size mm of Representative Volume Element of a micro-concrete is considered (Fig. 1); it is representative of a micro-concrete of Young s modulus E = 0000 MPa and Poisson s ration ν = 0.2. The number of particles is,072 and the number of degrees of freedom 24,576. The crack pattern obtained at the end of the equi-triaxial loading (quasi-rupture) is the one given in Fig. 1, with a number of beams to break before failure of 8,000. It has been shown from such computations (Fig. 2) that at low damage the dependence K(damage) is close to be the same linear function K K(1 1.2 d H ), d H = 1 tr d, (8.2) of the hydrostatic damage d H whatever the stress triaxiality. For each mark of the Figures the components of the damage tensor d = 1 Φ 2 have been measured by means of repeated numerical elastic loading-unloading sequences performed in uniaxial tension on the mm cube (even for the triaxial loading), using then the coupling of elasticity with anisotropic damage given by elasticity law (7.) with one non zero principal stress σ i = σ, the two others σ j i = 0. Note that considering Eq. (8.2) at high damage level means that in uniaxial tension the bulk modulus K fully vanishes at tr d = /1.2 = 2.5, i.e. at maximum principal damage max d i larger than 1 ( K cannot vanish then

10 D n fact versus D H for comparison) and versus the norm D ig. 1 Crack patterns for the two samples (left 8 8 8, right 16 16) Figure 1: Discrete Element sample considered as Representative Volume Element ( mm cube, from Delaplace and Desmorat (2008)) and crack pattern at the end of equi-triaxial loading (high lue η = damage 0correspondstounphysicalresponsewith level). o damage developed in tritension. Conclusion Triaxial loading Uniaxial loading he popularity and the use of a constitutive model epend on its robustness, 0.6 its simplicity, and its easiness implement in a numerical K code. Concerning simplicy, the number, thek physical 0.4 meaning and the identificaon easiness of material parameters is an important feare to be considered. 0.2Both the discrete and anisotropic amage models have been developed with respect to ese considerations with quite a reduced number of arameters introduced. 12 ~ d H = 1 tr d Figure 2: Effective bulk modulus K from Discrete Element computations as a function of hydrostatic damage d H (from Delaplace and Desmorat (2008)). K K Uniaxial Triaxial loading Uniaxial loading End of uniaxial loading at d1=1, d2=d=0, dh=1/ ~ End of equi-triaxial loading at d1=d2=d=1, dh=1/1.2 as principal damages d i therefore here d 1 are always bounded by 1, see also Fig. (5) for η = 1.2). This corresponds to a quite high (spurious) elastic stiffness K which is kept at rupture (see Fig. ). Enforcing then gradually K 0 but allowing for damage tensor d to evolve up to second order unit tensor 1 in an adequate procedure for the numerical control of rupture is a solution which leads to numerical difficulties in Finite Element computations (Badel et al., 2007; Ragueneau et al, 2008; Leroux, 2012). The relation (8.1) relates the effective bulk modulus to damage tensor Φ and hydrostatic sensitivity parmeter η (obtained from proposed micro-mechanically based second order damage framework by function g defined in Eq. (7.7)). This relation implies that the effective (damaged) bulk modulus K vanishes exactly when the maximum eigenvalue of damage tensor d = 1 Φ 2 is equal to 1, whatever the stress multiaxiality and without the need of a procedure bounding the damage eigenvalues to 1. To illustrate this property, we describe below the three particular cases of uniaxial, equi-biaxial and equi-triaxial tension loadings. It is shown that, at low damage in those three loading cases, one recovers the expression K = K(1 ηd H ) due to Lemaitre et al. (2000). In uniaxial tension the damage tensor of quasi-brittle materials is classically d = diag[d 1, 0, 0] (Lubarda and Krajcinovic, 199; Krajcinovic, 1996) and Φ 1 = (1 d 1 ) 1/2 1, Φ 2 = Φ = 1 so that the effective bulk modulus (8.1) has for expression K = 5K(1 d H ) 5 (15 8η)d H 2η ( 1 1 d H ). (8.) In this uniaxial loading, the maximum principal damage is d 1 = d H. The equi-triaxial tension case corresponds to spherical damage tensors d = d H 1, Φ = Φ H 1 = (1 d H ) 1/2 1, with thus tr Φ 2 = 1 (tr Φ)2 = /(1 d H ), and so Eq. (8.1) rewrites as K = K(1 d H) 1 (1 η)d H. (8.4) In this equi-triaxial loading, the maximum principal damages are d i = d H. Note that, then, the value η = 1 leads to the linear law K = K(1 d H ) over the whole range of damage. 0.2 Equi-triaxial d H = 1 tr d Figure : Effective bulk modulus K from initial modeling K = K(1 1.2 d H ) (dashed lines) as a function of hydrostatic damage. 10 We can also consider equi-biaxial tension for which d 1 = d 2 0, d = 0, Φ 1 = Φ 2 = (1 d 1 ) 1/2 0, Φ = 1. The effective bulk modulus (8.1) becomes then 5K(2 d H ) K = 10 (15 1η)d H 2η ( 2 ). 4 6d H (8.5) In this equi-biaxial loading, the maximum principal damages are d 1 = d 2 = 2 d H.

11 K K Triaxial loading Uniaxial loading ~ Figure 4 shows that, for the whole range of hydrostatic damage, the loss of bulk modulus (8.1) with η = 1.2 behave in a similar manner as the one of Fig. 2 obtained from Discrete Element computations. As expected, the effective (damaged) bulk modulus K vanishes exactly when maximum principal damage(s) are equal to 1 (solid lines in Fig. 5) in all these loading cases. The first order expansions in d H (at small damage) of the three Equations (8.), (8.4) and (8.5) gives for all three uniaxial, equi-biaxial and equi-triaxial loading cases K K(1 ηd H ) (8.6) 0.2 Uniaxial Equi-triaxial (dashed lines in Fig. 5) i.e. the exact expression introduced by Lemaitre et al. (2000) in a fully phenomenological manner (recovering Eq. (8.2) when η = 1.2 is set) d H = 1 tr d Figure 4: Effective bulk modulus K from (8.1) as a function of hydrostatic damage d H with d = 1 Φ 2 (for value η = 1.2 representative of a micro-concrete). 9. Conclusion Some mathematical tools such as the harmonic product and the harmonic factorization into lower order tensors have been presented. Together with the notion of mechanically accessible directions for measurements, this has allowed us to derive, at harmonic order 4, both a crack density expansion Ω(n) and a micro-mechanics based damage framework that makes use of second-order tensorial variables only, instead of fourth-order in standard micromechanics based approaches. The hydrostatic sensitivity obtained from such a second-order damage framework is shown to have the sought property of leading to the vanishing of the effective (damaged) bulk modulus at maximum principal damage max d i exactly equal to 1. K K Uniaxial Equi-biaxial Appendix A. Spherical/Deviatoric harmonic decomposition Different practical expressions of the harmonic decomposition of a fourth order tensor of elasticity type exist (Backus, 1970; Onat, 1984; Auffray, 2017). We use here the spherical/deviatoric harmonic decomposition introduced by Auffray (2017) of a fourth order tensor T of the elasticity type 0.2 Equi-triaxial max d i Figure 5: Effective bulk modulus K from (8.1) (solid lines) and from initial modeling K = K(1 η d H ) (dashed lines) as a function of hydrostatic damage d H with d = 1 Φ 2 (for value η = 1.2 representative of a micro-concrete). T = α β J + 1 c + c 1 ( + 2 (1 b + b 1) 2 ) (1 b + b 1) + H, (A.1) where, di and vo being the dilatation and voigt tensors of T and di and vo their deviatoric parts, di = tr 12 T, vo = tr 1 T, α = tr di, β = ( tr di + tr vo), 0 c = 1 di, b = 1 7 ( 2di + vo ). 11

12 Introducing two second symmetric second order tensors y and x such that y = T : x, the following equalities are obtained : and x : T : x = α(tr x) 2 + 2βx : x + 2(tr x)(c : x ) + 4b : x 2 + x : H : x, y = α(tr x)1 + 2βx + (c : x )1 + (tr x)c tr y = (α tr x + c : x ), + 2 [(b.x ) + (x.b ) ] + H : x, y = (tr x)c + 2βx + 2 [(b.x ) + (x.b ) ] + H : x. with 1 E D =g(φ) 1 9K β 1 2G J + 1 ( 1 b + b 1 2 ) G (1 b + b 1) + 1 2G Φ * Φ. A straightforward calculation leads to the property (tr 12 D) = 0, (tr 1 D) = 7(1 + ν) b. Appendix B. Harmonic decomposition of D References References From equation (7.1), the elasticity law is obtained as with fourth order tensor S = g(φ) 9K G G, (B.1) G = Φ Φ (tr Φ2 ) (1 Φ2 + Φ 2 1). (B.2) The harmonic decomposition (A.1) of Φ Φ reads Φ Φ = α β J + 1 c + c 1 ( + 2 (1 b + b 1) 2 ) (1 b + b 1) + H with di = Φ 2, vo = 1 [ (tr Φ)Φ + Φ 2], 2 α = 1 9 tr(φ2 ), β = 1 [ tr(φ 2 ) + (tr Φ) 2], 60 c = 1 (Φ2 ), b = 1 [ (tr Φ)Φ (Φ 2 ) ]. 14 As (Φ Φ) S = Φ Φ, we get Noting that H = ((Φ Φ) S ) 0 = (Φ Φ) 0 = Φ * Φ. α c + c 1 tr Φ2 = ( 1 (Φ 2 ) + (Φ 2 ) 1 ) 9 tr Φ2 = ( 1 (Φ 2 ) + (Φ 2 ) 1 ) 9 equations (B.1) and (B.2) lead to S = 1 9K G J + 1 E D, (B.) 12 Auffray N., Géométrie des espaces de tenseurs, application à l élasticité anisotrope classique et généralisée, Habilitation à diriger des recherches, Université Paris-Est. Andrieux, S., Bamberger, Y., Marigo, J.-J., A model for microcracked material for concrete and rocks. J. Theo. Appl. Mech. 5, Backus, G., A geometrical picture of anisotropic elastic tensors. Reviews of Geophysics 8 (), Badel, P., Godard, V., Leblond, J., Application of some anisotropic damage model to the prediction of the failure of some complex industrial concrete structure. International Journal of Solids and Structures 44, Baerheim, R., Classification of symmetry by means of maxwell multipoles. Q. J. Mech. Appl. Math. 51, Brunig, M., 200. An anisotropic ductile damage model based on irreversible thermodynamics. Int. J. Plasticity 19, Burr, A., Hild, F., Leckie, F. A., Micro-mechanics and continuum damage mechanics. Archive of Applied Mechanics 65, Calloch, S., Essais triaxiaux non proportionnels et ingénierie des modéles de plasticité cyclique. Ph.D. thesis. Calloch, S., Marquis, D., Triaxial tension-compression tests for multiaxial cyclic plasticity. International Journal of Plasticity 15, Carol, I., Rizzi, E., Willam, K., On the formulation of anisotropic elastic degradation. part i: Theory based on a pseudologarithmic damage tensor rate. Int. J. of Solids and Structures 8 (4), Cauchy, A. L., 1828a. De la pression ou tension dans un système de points matériels. Exercices de mathématiques, Cauchy, A. L., 1828b. Sur l équilibre et le mouvement d un système de points matériels sollicités par des forces d attraction ou de repulsion mutuelle. Exercices de mathématiques, Chaboche, J.-L., Description thermodynamique et phénoménologique de la viscoplasticité cyclique avec endommagement., Doctorat d Etat Université Paris VI and Technical report Onera Chaboche, J.-L., Le concept de contrainte effective appliqué à l élasticité et à la viscoplasticité en présence d un endommagement anisotrope. In: Boehler, J.-P. (Ed.), Colloque Int. CNRS 295, Villard de Lans. Martinus Nijhoff Publishers and Editions du CNRS, 1982, pp Chaboche, J., Anisotropic creep damage in the framework of continuum damage mechanics. Nuclear Engineering and Design 79 (), Chambart, M., Desmorat, R., Gatuingt, F., Intrinsic dissipation of a modular anisotropic damage model: Application to concrete under impact. Eng. Fract. Mech. 127,

13 Cordebois, J., Sidoroff, F., Endommagement anisotrope en élasticité et plasticité. J. Meca. Th. Appl., Special Volume, Cormery, F., Welemane, H., A stress-based macroscopic approach for microcracks unilateral effect. Comp. Mat. Sci. 47, Cundall P.A., Strack O.D.L. (1979) A discrete numerical model for granular assemblies. Geotechnique, 29: Delaplace, A. and Desmorat, R. (2008). Discrete D model as complimentary numerical testing for anisotropic damage. International Journal of Fracture, 148, Desmorat, B., Desmorat, R., Tensorial polar decomposition of 2D fourth order tensors. C. R. Mécanique 4, Desmorat, B., Desmorat, R., Second order tensorial framework for 2D medium with open and closed cracks. European Journal of Mechanics A/Solids 58, Desmorat, R., Positivity of intrinsic dissipation of a class of nonstandard anisotropic damage models. C. R. Mecanique 4 (10), Desmorat, R., Anisotropic damage modeling of concrete materials. Int. J. Dam. Mech. 6, Desmorat, R., Gatuingt, F., Ragueneau, F., Nonlocal anisotropic damage model and related computational aspects for quasi-brittle materials. Engineering Fracture Mechanics 74 (10), Desmorat, R., Otin, S., Cross-identification isotropic/anisotropic damage and application to anisothermal structural failure. Engineering Fracture Mechanics 75 (11), Dormieux, L., Kondo, D., Micromechanics of fracture and damage. Wiley-ISTE. Forte, S., Vianello, M., Symmetry classes for elasticity tensors. J. Elasticity 4 (2), Halm, D., Dragon, A., An anisotropic model of damage and frictional sliding for brittle materials. Eur. J. Mech. A/Solids 17, Herrmann H.J., Roux S. (1990), Statistical models for the fracture of disordered media. Elsevier Science Publishers, Amsterdam. Ju, J., On energy-based coupled elastoplastic damage theories: Constitutive modeling and computational aspects. Int. J. Sol. Struct. 25, Kachanov, M., On continuum theory of medium with cracks. Mech. Solids 7, Kachanov, M., 199. Elastic solids with many cracks and related problems. Vol. 1. J. Hutchinson and T. Wu Ed., Academic Press Pub., pp Kanatani, K., Distribution of directional data and fabric tensors. Int. J. Eng. Sci. 22, Kattan, P., Voyiadjis, G., A coupled theory of damage mechanics and finite strain elasto-plasticity i. damage and elastic deformations. Int. J. Eng. Sci. 28, Krajcinovic, D., Damage Mechanics. Applied Mathematics and Mechanics. North Holland. Ladevèze, P., 198. Sur une théorie de l endommagement anisotrope, internal report 4 of lmt-cachan. Ladevèze, P., Modeling and simulation of the mechanical behavior of cmcs. High-Temperature Ceramic-Matrix Composites 47, 5 6. Leckie, F. A., Onat, E. T., Tensorial nature of damage measuring internal variables. J. Hult and J. Lemaitre eds, Springer Berlin, Ch. Physical Non-Linearities in Structural Analysis, pp Lemaitre, J., Chaboche, J.-L., Mécanique des matériaux solides. Dunod, english translation 1990 Mechanics of Solid Materials Cambridge University Press. Lemaitre, J., Desmorat, R., Engineering Damage Mechanics : Ductile, Creep, Fatigue and Brittle Failures. Springer. Lemaitre, J., Desmorat, R., Sauzay, M., Anisotropic damage law of evolution. Eur. J. Mech., A/ Solids 19, Leroux A., Modèle multiaxial dendommagement anisotrope : gestion numérique de la rupture et application la ruine de structures en béton armé sous impacts. PhD ENS-Cachan, France. Love, A., A Treatise on the Mathematical Theory of Elasticity. 2 nd edition, Cambridge University Press. Lubarda, V., Krajcinovic, D., 199. Damage tensors and the crack density distribution. Int. J. Solids Structures 0, Menzel, A., Ekh, M., Steinmann, P., Runesson, K., Anisotropic damage coupled to plasticity: Modelling based on the effective configuration concept. Int. J. Numer. Meth. Eng. 54, Menzel, A., Steinmann, P., A theoretical and computational setting for anisotropic continuum damage mechanics at large strains. Int. J. Solids Struct. 8, Murakami S. and Ohno M., A constitutive equation of creep damage in polycrystalline metals, IUTAM Coll. Euromech 111 Constitutive Modelling in Inelasticity, Marienbad, Czechoslovakia, Murakami, S., Mechanical modeling of material damage. ASME J. Appl. Mech. 55, Navier, L., Mémoire sur les lois de l équilibre et du mouvement des solides élastiques. Mémoires de l Académie Royale des Sciences de l Institut National 7, Olive, M., Kolev, B., Desmorat, B., Desmorat, R., Harmonic factorization and reconstruction of the elasticity tensor, Journal of Elasticity, in Press. DOI: /s y. Also URL (2016). Olivier-Leblond C., Delaplace A., Ragueneau F., Modelling of three-dimensional crack patterns in deep reinforced concrete structures, Engineering Structures, 8, Onat, E. T., Effective properties of elastic materials that contain penny shaped voids. Int. J. Eng. Sci. 22, Papa, E., Taliercio, A., Anisotropic damage model for the multi-axial static and fatigue behaviour of plain concrete. Engineering Fracture Mechanics 55, Poisson, S. D., Mémoire sur l équilibre et le mouvement des corps élastiques. Mémoires de l Académie Royale des Sciences de l Institut National 8, Ramtani, S., Berthaud, Y., Mazars, J., Orthotropic behaviour of concrete with directional aspects: modelling and experiments. Nucl. Eng. Design 1, Ragueneau F., Desmorat R., Gatuingt F., Anisotropic damage modelling of biaxial behaviour and ruprure of concrete structures, Computers and Concrete, 5(4), 2008, pp Schlangen E., Garboczi E.J. (1997) Fracture simulations of con- crete using lattice models: computational aspects. Eng Fract Mech 57(2/):19 2. Schouten, J., Tensor Analysis for Physicists. Clarendon Press. Spencer, A., A note on the decomposition of tensors into traceless symmetric tensors. Int. J. Engng Sci. 8, Steinmann, P., Carol, I., A framework for geometrically nonlinear continuum damage mechanics. Int. J. Eng. Sci. 6, Sylvester, J. J., Note on spherical harmonics. In: Collected Mathematical Papers. Vol.. Cambridge University Press, pp Tikhomirov, D., Niekamp, R., Stein, E., On three-dimensional microcrack density distribution. Z. Angew. Math. Mech. 81, 16. Vakulenko, P., Kachanov, M., Continuum theory of medium with cracks. Mech. of Solids (Eng. Transl. of Izvestia AN SSSR, Mekhanika Tverdogo Tela) 6, Van Mier J.G.M., Van Vliet M.R.A., Wang T.K. (2002) Fracture mechanisms in particle composites: statistical aspects in lattice type analysis. Mech Mater, 4: Vannucci, P., Plane anisotropy by the polar method. Meccanica 40, Vannucci, P., Desmorat, B., Plane anisotropic rari-constant materials. Math. Meth. Appl. Sci. 9, Verchery, G., Les invariants des tenseurs d ordre 4 du type de l élasticité. In: Boehler, J.-P. (Ed.), Colloque Int. CNRS 295, Villard de Lans. Martinus Nijhoff Publishers and Editions du CNRS, 1982, pp Voyiadjis, G.Z., Kattan, P.I., Damage mechanics with fabric tensors. Mechanics of Advanced Materials and Structures 1 (4),