Iterative convergence in phase-field brittle fracture computations: exact line search is all you need

Abstract

Variational phase-field models of brittle fracture pose a local constrained minimization problem of a non-convex energy functional. In the discrete setting, the problem is most often solved by alternate minimization, exploiting the separate convexity of the energy with respect to the two unknowns. This approach is theoretically guaranteed to converge, provided each of the individual subproblems is solved successfully. However, strong non-linearities of the energy functional may lead to failure of iterative convergence within one or both subproblems. We analyze and visualize the energy along Newton directions to illustrate why Newton's method without line search fails. Motivated by this, we propose to employ an exact line search algorithm based on bisection, which (under certain conditions) can guarantee global convergence of Newton's method for each subproblem and consequently the successful determination of critical points of the energy through the alternate minimization scheme. Through several benchmark tests computed with various strain energy decompositions and two strategies for the enforcement of the irreversibility constraint in two and three dimensions, we demonstrate the robustness of the approach and assess its efficiency in comparison with other commonly used line search algorithms. With the outlined approach, we are able to compute the especially demanding Brazilian test featuring contact in 3D with the star-convex model.

References in corpus (13)

• Construction of arbitrary order finite element degree-of-freedom maps on polygonal and polyhedral cell meshes
• On penalization in variational phase-field models of brittle fracture
• Composing Scalable Nonlinear Algebraic Solvers
• 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
• Recursive multilevel trust region method with application to fully monolithic phase-field models of brittle fracture
• An efficient and robust monolithic approach to phase-field quasi-static brittle fracture using a modified Newton method
• 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
• Variational phase-field modeling of cohesive fracture with flexibly tunable strength surface
• Truncated Nonsmooth Newton Multigrid for phase-field brittle-fracture problems, with analysis
• Greedy Newton: Newton's Method with Exact Line Search