DeepAI AI Chat
Log In Sign Up

A "parallel universe" scheme for crack nucleation in the phase field approach to fracture

by   Yihao Chen, et al.

Crack nucleation is crucial in many industrial applications. The phase field method for fracture transforms the crack nucleation problem into a minimization problem of the sum of the elastic potential energy and the crack surface energy. Due to the polyconvexity of the formulation, starting from a crackless solid, a standard Newton iteration may lead to a solution with no crack, even though a cracked solution has a lower total energy. As such, the critical load for cracking is highly overestimated. Here, we propose an algorithm termed "parallel universe" algorithm to capture the global minimum. This algorithm has two key ingredients: (a) a necessary condition for cracking solely based on the current crackless solution, and (b) beginning from when this condition is met, Newton iteration with two initial guesses, a crackles one and a cracked one, will both be performed and the converged candidate solution with lower energy is accepted as the solution at that load step. Once the cracked candidate solution is accepted, the crackless one is discarded, i.e., only one universe is retained. This cracked initial guess is obtained only once for all load steps by solving a series of similar minimization problems with a progressively reduced critical crack energy release rate. Numerical examples with isotropic and anisotropic critical crack energy release rates indicate that the proposed algorithm is more reliable (as there is no need to retrace) and more efficient than the standard Newton iteration and a well-known backtracking algorithm.


page 12

page 14

page 19


Truncated Nonsmooth Newton Multigrid for phase-field brittle-fracture problems

We propose the Truncated Nonsmooth Newton Multigrid Method (TNNMG) as a ...

Parallel energy stable phase field simulations of Ni-based alloys system

In this paper, we investigate numerical methods for solving Nickel-based...

Munchausen Iteration

We present a method for solving polynomial equations over idempotent ome...

Efficient numerical methods for computing the stationary states of phase field crystal models

Finding the stationary states of a free energy functional is an importan...

A Global Newton-Type Scheme Based on a Simplified Newton-Type Approach

Globalization concepts for Newton-type iteration schemes are widely used...