Damage and fracture evolution in brittle materials by shape optimization methods

Authors:
Grégoire Allaire;François Jouve;Nicolas Van Goethem
Affiliations:
Centre de Mathématiques Appliquées, ícole Polytechnique, 91128 Palaiseau, France;Laboratoire J.-L. Lions, Université Paris Diderot (Paris 7), 75205 Paris, France;Centre de Mathématiques Appliquées, ícole Polytechnique, 91128 Palaiseau, France and Universidade de Lisboa, Faculdade de Ciências, Departamento de Matemática, CMAF, Av. P ...
Venue:
Journal of Computational Physics
Year:
2011

Citing 7
Cited 1

Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations

Journal of Computational Physics
Image segmentation by variational methods: Mumford and Shah functional and the discrete approximations

SIAM Journal on Applied Mathematics
On the Topological Derivative in Shape Optimization

SIAM Journal on Control and Optimization
A remark on computing distance functions

Journal of Computational Physics
The Topological Asymptotic for PDE Systems: The Elasticity Case

SIAM Journal on Control and Optimization
Etude de Problème d'Optimal Design

Proceedings of the 7th IFIP Conference on Optimization Techniques: Modeling and Optimization in the Service of Man, Part 2
Direct Elastic Imaging of a Small Inclusion

SIAM Journal on Imaging Sciences

Phase field method to optimize dielectric devices for electromagnetic wave propagation

Journal of Computational Physics

Quantified Score

Hi-index	31.45

Visualization

Abstract

This paper is devoted to a numerical implementation of the Francfort-Marigo model of damage evolution in brittle materials. This quasi-static model is based, at each time step, on the minimization of a total energy which is the sum of an elastic energy and a Griffith-type dissipated energy. Such a minimization is carried over all geometric mixtures of the two, healthy and damaged, elastic phases, respecting an irreversibility constraint. Numerically, we consider a situation where two well-separated phases coexist, and model their interface by a level set function that is transported according to the shape derivative of the minimized total energy. In the context of interface variations (Hadamard method) and using a steepest descent algorithm, we compute local minimizers of this quasi-static damage model. Initially, the damaged zone is nucleated by using the so-called topological derivative. We show that, when the damaged phase is very weak, our numerical method is able to predict crack propagation, including kinking and branching. Several numerical examples in 2d and 3d are discussed.