Overall behavior of elastic materials with embedded cohesive zone models:rigorous bounds and mesh size-to-cohesive properties relationships

Nawfal Blal (1, 2) , Loïc Daridon (3) , Yann Monerie (2, 4) , Stéphane Pagano (1, 2)

1 LMGC - Laboratoire de Mécanique et Génie Civil
2 MIST - Laboratoire de micromécanique et intégrité des structures
3 M3 - Modélisation Mathématique en Mécanique
4 IRSN - Institut de Radioprotection et de Sûreté Nucléaire

Incorporating softening behaviors in scales transition problems is a complex task, especially due to strain and stress localization. Micromechanical approaches (Gurson-type models) have to invoke some ad hoc coalescence models in order to reproduce relevant overall fracture [1]. Numerical simulations based on continuum damage mechanics lead to a vanishing energy (energy dissipated due to strain-softening damage) when the mesh tends to zero and nonlocal damage theories have to be formulated [2]. Recently, a renewed interest is noticed in the literature for these nonlocal (explicit, integral or variational) damage models in particular concerning the diffuse microcracking observed in quasi-brittle materials before the onset of a localized macrocrack (see [3], [4], [5] and references therein). The new micromechanical approach, presented here, allows multiscale transitions in softening problems (up to fracture) without any nonlocal or first gradient theory. This approach, first proposed by [6], is based on a cohesive-volumetric finite element method together with rigorous theoretical upper bounds on the overall strain potential. The surface damage is defined through a modified secant modulus of any intrinsic cohesive zone model (traction-separation law with initial stiffness). Considering each face of a three dimensional mesh as a penny shaped damageable inclusion, a linear comparison composite (LCC) is defined using the variational approach of [7]. A convenient estimate of the LCC properties and an asymptotic analysis lead to a micromechanical damage model based on the cohesive zone approach. This model includes a dependence to the mesh size allowing finite values of the overall properties (stiffness, dissipated energy, maximal stress, etc.) at mesh convergence. This model used in complex multiscale problems as fracture of functionally graded materials looks very promising.

Domaines: Mécanique des solides [physics.class-ph]; Mécanique des solides [physics.class-ph]

Nawfal Blal, Loïc Daridon, Yann Monerie, Stéphane Pagano. Overall behavior of elastic materials with embedded cohesive zone models:rigorous bounds and mesh size-to-cohesive properties relationships. ESIS-TC Workshop on Numerical Methods for Materials Failure, 2011, France. ⟨hal-00628017⟩