This repository contains an implementation of our paper on Computer Methods in Applied Mechanics and Engineering: "Mapped phase field method for brittle fracture". We proposed to use a local reparametrization of the physical domain near crack, so that without modifying the finite element mesh we obtain higher accuracy for crack simulation based on the phase field method.

To run the cases, do

python -m src.cases.driver

under the root directory (/crack).

You may need to create folders to hold data, e.g.,

crack/data/pvd/brittle

For educational purposes of building phase field models with FEniCS, you may find the file

crack/src/exp/unit_test.py

useful.

The following vedios correspond to the four MPFM numerical examples in "Section 5. Numerical examples" in the manuscript.

The left image shows x10 displacement, and the right image shows the reconstruced discretization of the physical domain.

The left image shows x10 displacement, and the right image shows the reconstruced discretization of the physical domain.

The left image shows x5 displacement, and the right image shows the reconstruced discretization of the physical domain.

The left image shows x100 displacement, and the right image shows the reconstruced discretization of the physical domain.

If you find our work useful in your research, please consider citing us.