Premium
Configurational‐Force‐Based Adaptive FE Solver for a Phase Field Model of Fracture
Author(s) -
Welschinger Fabian,
Hofacker Martina,
Miehe Christian
Publication year - 2010
Publication title -
pamm
Language(s) - English
Resource type - Journals
ISSN - 1617-7061
DOI - 10.1002/pamm.201010330
Subject(s) - classification of discontinuities , solver , finite element method , adaptive mesh refinement , complex fracture , extended finite element method , computer science , fracture (geology) , fracture mechanics , field (mathematics) , phase (matter) , mathematical optimization , mathematics , computational science , materials science , structural engineering , physics , mathematical analysis , engineering , quantum mechanics , pure mathematics , composite material
The computational modeling of failure mechanisms in solids due to fracture based on sharp crack discontinuities suffers in situations with complex crack topologies. This can be overcome by diffusive crack modeling, based on the introduction of a crack phase field as outlined in [1, 2]. Following these formulations, we outline a thermodynamically consistent framework for phase field models of crack propagation in elastic solids, develop incremental variational principles and, as an extension to [1, 2], consider their numerical implementations by an efficient h ‐adaptive finite element method. A key problem of the phase field formulation is the mesh density, which is required for the resolution of the diffusive crack patterns. To this end, we embed the computational framework into an adaptive mesh refinement strategy that resolves the fracture process zones. We construct a configurational‐force‐based framework for h ‐adaptive finite element discretizations of the gradient‐type diffusive fracture model. We develop a staggered computational scheme for the solution of the coupled balances in physical and material space. The balance in the material space is then used to set up indicators for the quality of the finite element mesh and accounts for a subsequent h ‐type mesh refinement. The capability of the proposed method is demonstrated by means of a numerical example. (© 2010 Wiley‐VCH Verlag GmbH & Co. KGaA, Weinheim)