Details
Original language | English |
---|---|
Pages (from-to) | 577-611 |
Number of pages | 35 |
Journal | Computer Methods in Applied Mechanics and Engineering |
Volume | 325 |
Publication status | Published - 1 Oct 2017 |
Externally published | Yes |
Abstract
Our goal in this work is to develop and compare modified Newton methods for fully monolithic quasi-static brittle phase-field fracture propagation. In variational phase-field fracture, a smoothed phase-field indicator variable denoting the crack path is coupled to elasticity. Moreover, a crack irreversibility condition is incorporated. To develop a fully monolithic scheme is an extremely challenging task since the underlying problem is non-convex and the Jacobian of Newton's method is indefinite. To split the problem using alternating minimization, thus a partitioned approach, is a possible resort. However, there are good reasons to consider monolithic approaches such as for example robustness and efficiency. Although an error-oriented Newton method can cope with a larger variety of configurations, it appears that this method is not always robust and also not always efficient. Inspired by nonlinear flow problems, as alternative, we develop a modified Newton scheme in which globalization is based on a dynamic modification of the Jacobian matrix rather than utilizing line-search or trust-region strategies. This variation switches smoothly between full Newton and Newton-like steps. In several 2D and 3D numerical examples, all of them with different characteristic features, our modified Newton solver is compared to a backtracking line-search Newton method, another line-search method monitoring the global energy and allowing for negative curvatures, and to already published results of an error-oriented version. These computations also include further modifications of Newton's method and detailed discussions why certain schemes either work or fail. Revisiting all findings, the main outcome of this paper is that the modified Newton scheme with Jacobian modification is currently the only method that works in a robust and efficient way for all provided examples, whereas line-search schemes or the error-oriented scheme show deficiencies for certain configurations.
Keywords
- Benchmark tests, Inexact augmented Lagrangian, Jacobian modification, Line-search, Modified Newton's method, Phase-field fracture propagation
ASJC Scopus subject areas
- Engineering(all)
- Computational Mechanics
- Engineering(all)
- Mechanics of Materials
- Engineering(all)
- Mechanical Engineering
- Physics and Astronomy(all)
- Computer Science(all)
- Computer Science Applications
Cite this
- Standard
- Harvard
- Apa
- Vancouver
- BibTeX
- RIS
In: Computer Methods in Applied Mechanics and Engineering, Vol. 325, 01.10.2017, p. 577-611.
Research output: Contribution to journal › Article › Research › peer review
}
TY - JOUR
T1 - Modified Newton methods for solving fully monolithic phase-field quasi-static brittle fracture propagation
AU - Wick, Thomas
N1 - Publisher Copyright: © 2017 Elsevier B.V. Copyright: Copyright 2018 Elsevier B.V., All rights reserved.
PY - 2017/10/1
Y1 - 2017/10/1
N2 - Our goal in this work is to develop and compare modified Newton methods for fully monolithic quasi-static brittle phase-field fracture propagation. In variational phase-field fracture, a smoothed phase-field indicator variable denoting the crack path is coupled to elasticity. Moreover, a crack irreversibility condition is incorporated. To develop a fully monolithic scheme is an extremely challenging task since the underlying problem is non-convex and the Jacobian of Newton's method is indefinite. To split the problem using alternating minimization, thus a partitioned approach, is a possible resort. However, there are good reasons to consider monolithic approaches such as for example robustness and efficiency. Although an error-oriented Newton method can cope with a larger variety of configurations, it appears that this method is not always robust and also not always efficient. Inspired by nonlinear flow problems, as alternative, we develop a modified Newton scheme in which globalization is based on a dynamic modification of the Jacobian matrix rather than utilizing line-search or trust-region strategies. This variation switches smoothly between full Newton and Newton-like steps. In several 2D and 3D numerical examples, all of them with different characteristic features, our modified Newton solver is compared to a backtracking line-search Newton method, another line-search method monitoring the global energy and allowing for negative curvatures, and to already published results of an error-oriented version. These computations also include further modifications of Newton's method and detailed discussions why certain schemes either work or fail. Revisiting all findings, the main outcome of this paper is that the modified Newton scheme with Jacobian modification is currently the only method that works in a robust and efficient way for all provided examples, whereas line-search schemes or the error-oriented scheme show deficiencies for certain configurations.
AB - Our goal in this work is to develop and compare modified Newton methods for fully monolithic quasi-static brittle phase-field fracture propagation. In variational phase-field fracture, a smoothed phase-field indicator variable denoting the crack path is coupled to elasticity. Moreover, a crack irreversibility condition is incorporated. To develop a fully monolithic scheme is an extremely challenging task since the underlying problem is non-convex and the Jacobian of Newton's method is indefinite. To split the problem using alternating minimization, thus a partitioned approach, is a possible resort. However, there are good reasons to consider monolithic approaches such as for example robustness and efficiency. Although an error-oriented Newton method can cope with a larger variety of configurations, it appears that this method is not always robust and also not always efficient. Inspired by nonlinear flow problems, as alternative, we develop a modified Newton scheme in which globalization is based on a dynamic modification of the Jacobian matrix rather than utilizing line-search or trust-region strategies. This variation switches smoothly between full Newton and Newton-like steps. In several 2D and 3D numerical examples, all of them with different characteristic features, our modified Newton solver is compared to a backtracking line-search Newton method, another line-search method monitoring the global energy and allowing for negative curvatures, and to already published results of an error-oriented version. These computations also include further modifications of Newton's method and detailed discussions why certain schemes either work or fail. Revisiting all findings, the main outcome of this paper is that the modified Newton scheme with Jacobian modification is currently the only method that works in a robust and efficient way for all provided examples, whereas line-search schemes or the error-oriented scheme show deficiencies for certain configurations.
KW - Benchmark tests
KW - Inexact augmented Lagrangian
KW - Jacobian modification
KW - Line-search
KW - Modified Newton's method
KW - Phase-field fracture propagation
UR - http://www.scopus.com/inward/record.url?scp=85027698782&partnerID=8YFLogxK
U2 - 10.1016/j.cma.2017.07.026
DO - 10.1016/j.cma.2017.07.026
M3 - Article
AN - SCOPUS:85027698782
VL - 325
SP - 577
EP - 611
JO - Computer Methods in Applied Mechanics and Engineering
JF - Computer Methods in Applied Mechanics and Engineering
SN - 0045-7825
ER -