This is a finite element code based on the finite element library deal.II to simulate crack propagation in elastic and porous media.
Features:
- phase field approach for the crack location
- support for 2d and 3d computations
- primal-dual active set method for the irreversibility constraint
- novel adaptive mesh refinement technique
This project was originally developed for the papers mentioned below, but has been extended considerably since then. Please cite these papers if you make use of our work. Thank you!
The original code versions are available as separate branches in this repository:
T. Heister, M. F. Wheeler, T. Wick: A primal-dual active set method and predictor-corrector mesh adaptivity for computing fracture propagation using a phase-field approach. Comp. Meth. Appl. Mech. Engrg., Vol. 290 (2015), pp. 466-495 http://dx.doi.org/10.1016/j.cma.2015.03.009
A preprint is available here: http://www.math.clemson.edu/~heister/preprints/HeWheWi15_CMAME_accepted.pdf
T. Heister, T. Wick: Parallel solution, adaptivity, computational convergence, and open-source code of 2d and 3d pressurized phase-field fracture problems Proc. Appl. Math. Mech., 2018, e201800353 https://doi.org/10.1002/pamm.201800353
A preprint is available here: https://arxiv.org/abs/1806.09924
You need to install deal.II (see http://www.dealii.org) with external dependencies p4est and Trilinos. Then configure with:
cmake -D DEAL_II_DIR=/your/dealii-installation/ .
Compile with:
make
and finally run with:
mpirun -n 2 ./cracks parameters_sneddon_2d.prm
The code is published under GPL v2 or newer.
Authors: Timo Heister, Thomas Wick.