Mapped phase field method for simulating quasi-static brittle fracture based on FEniCS

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.

Dependency and Usage

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.

Demos

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

mode-I.mp4

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

mode-II.mp4

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

TPB.mp4

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

L_shape.mp4

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

Citation

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