An asynchronous variational integrator for the phase field approach to dynamic fracture

Zongwu Niu,Vahid Ziaei-Rad,Zongyuan Wu,Yongxing Shen

DOI: https://doi.org/10.48550/arXiv.2203.08797

2022-03-17

Abstract:The phase field approach is widely used to model fracture behaviors due to the absence of the need to track the crack topology and the ability to predict crack nucleation and branching. In this work, the asynchronous variational integrators (AVI) is adapted for the phase field approach of dynamic brittle fractures. The AVI is derived from Hamilton's principle and allows each element in the mesh to have its own local time step that may be different from others'. While the displacement field is explicitly updated, the phase field is implicitly solved, with upper and lower bounds strictly and conveniently enforced. In particular, the AT1 and AT2 variants are equally easily implemented. Three benchmark problems are used to study the performances of both AT1 and AT2 models and the results show that AVI for phase field approach significantly speeds up the computational efficiency and successfully captures the complicated dynamic fracture behavior.

Numerical Analysis

What problem does this paper attempt to address?

The problems that this paper attempts to solve are the computational efficiency and the capture of complex crack behaviors in dynamic fracture simulation. Specifically:

1. **Low computational efficiency**: When the traditional phase - field method is used to simulate dynamic brittle fracture, a very fine spatial discretization is required to resolve the regularization length scale, resulting in a very small critical time step and thus a very high computational cost.

2. **Capture of complex crack behaviors**: Under dynamic loading, the crack development process in materials is accompanied by rapidly changing external loads and rapid crack propagation, and inertial forces play an important role in this process. Therefore, how to accurately capture complex behaviors such as crack nucleation, propagation, and branching is a challenge.

To solve these problems, the author introduced the Asynchronous Variational Integrator (AVI) and applied it to the phase - field method to simulate dynamic brittle fracture. The main advantages of AVI include:

- **Allowing each finite element to have its own local time step**: This enables elements with smaller time steps to be updated more frequently, thus improving computational efficiency.

- **Explicitly updating the displacement field and implicitly solving the phase field**: The displacement field is updated by an explicit scheme, while the phase field is solved implicitly for the inequality - constrained optimization problem, ensuring that the upper and lower bounds of the phase - field values are strictly and conveniently imposed.

- **Adapting to AT1 and AT2 models**: These two variants can be implemented equally easily, and the performance on three benchmark problems shows that AVI significantly improves computational efficiency and successfully captures complex dynamic fracture behaviors.

Through these improvements, the paper demonstrates the potential of AVI in improving computational efficiency and capturing complex dynamic fracture behaviors.