A modified combined active-set Newton method for solving phase-field fracture into the monolithic limit
arXiv:2306.03501 · doi:10.1016/j.cma.2023.116170
Abstract
In this work, we examine a numerical phase-field fracture framework in which the crack irreversibility constraint is treated with a primal-dual active set method and a linearization is used in the degradation function to enhance the numerical stability. The first goal is to carefully derive from a complementarity system our primal-dual active set formulation, which has been used in the literature in numerous studies, but for phase-field fracture without its detailed mathematical derivation yet. Based on the latter, we formulate a modified combined active-set Newton approach that significantly reduces the computational cost in comparison to comparable prior algorithms for quasi-monolithic settings. For many practical problems, Newton converges fast, but active set needs many iterations, for which three different efficiency improvements are suggested in this paper. Afterwards, we design an iteration on the linearization in order to iterate the problem to the monolithic limit. Our new algorithms are implemented in the programming framework pfm-cracks [T. Heister, T. Wick; pfm-cracks: A parallel-adaptive framework for phase-field fracture propagation, Software Impacts, Vol. 6 (2020), 100045]. In the numerical examples, we conduct performance studies and investigate efficiency enhancements. The main emphasis is on the cost complexity by keeping the accuracy of numerical solutions and goal functionals. Our algorithmic suggestions are substantiated with the help of several benchmarks in two and three spatial dimensions. Therein, predictor-corrector adaptivity and parallel performance studies are explored as well.
arXiv comment49 pages, 45 figures, 9 tables
References in corpus (15)
- A physics-informed variational DeepONet for predicting the crack path in brittle materials
- The deal.II finite element library: design, features, and insights
- A unified Abaqus implementation of the phase field fracture method using only a user material subroutine
- An Adaptive Global-Local Approach for Phase-Field Modeling of Anisotropic Brittle Fracture
- An accelerated staggered scheme for phase-field modeling of brittle fracture
- Interior-point methods for the phase-field approach to brittle and ductile fracture
- An iterative staggered scheme for phase field brittle fracture propagation with stabilizing parameters
- An efficient and robust monolithic approach to phase-field quasi-static brittle fracture using a modified Newton method
- Recursive multilevel trust region method with application to fully monolithic phase-field models of brittle fracture
- Multilevel Global-Local techniques for adaptive ductile phase-field fracture
- Matrix-free multigrid solvers for phase-field fracture problems
- A least-squares Galerkin gradient recovery method for fully nonlinear elliptic equations
- Fast staggered schemes for the phase-field model of brittle fracture based on the fixed-stress concept
- Nonlinear Field-split Preconditioners for Solving Monolithic Phase-field Models of Brittle Fracture
- Truncated Nonsmooth Newton Multigrid for phase-field brittle-fracture problems, with analysis