Fracture and healing of elastomers: A phase-transition theory and numerical implementation

From an applications point of view, the proposed theory amounts to solving a system of two coupled and nonlinear PDEs for the deformation field and an order parameter, or phase field. A numerical scheme is presented to generate solutions for these PDEs in N=2 and 3 space dimensions. This is based on an efficient non-conforming finite-element discretization, which remains stable for large deformations and elastomers of any compressibility, together with an implicit gradient flow solver, which is able to deal with the large changes in the deformation field that can ensue locally in space and time from the nucleation of fracture. The last part of this paper is devoted to presenting sample simulations of the so-called Gent-Park experiment. Those are confronted with recent experimental results for various types of silicone elastomers.