This is a finite element code based in deal.II to simulate crack propagation in elastic and porous media.
Features:
- phase field approach for the crack location
- primal-dual active set method for the irreversibility constraint
- novel adaptive mesh refinement technique
This is the example program to accompany the paper
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
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
If the code crashes with an exception "ExcIO", create an empty directory called "output".
The code is published under GPL v2 or newer.
Authors: Timo Heister, Thomas Wick.