jievoo/cracks

Finite element code for crack propagation

Crack propagation

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

How to run

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".

Notes

The code is published under GPL v2 or newer.

Authors: Timo Heister, Thomas Wick.