Crack propagation in anisotropic brittle materials: from a phase-field model to a shape optimization approach

Tim Suchan,Chaitanya Kandekar,Wolfgang E. Weber,Kathrin Welker

2023-09-06

Abstract:The phase-field method is based on the energy minimization principle which is a geometric method for modeling diffusive cracks that are popularly implemented with irreversibility based on Griffith's criterion. This method requires a length-scale parameter that smooths the sharp discontinuity, which influences the diffuse band and results in mesh-sensitive fracture propagation results. Recently, a novel approach based on the optimization on Riemannian shape spaces has been proposed, where the crack path is realized by techniques from shape optimization. This approach requires the shape derivative, which is derived in a continuous sense and used for a gradient-based algorithm to minimize the energy of the system. Due to the continuous derivation of the shape derivative, this approach yields mesh-independent results. In this paper, the novel approach based on shape optimization is presented, followed by an assessment of the predicted crack path in anisotropic brittle material using numerical calculations from a phase-field model.

Optimization and Control

What problem does this paper attempt to address?

### Problems the Paper Aims to Solve

This paper aims to address the issue of simulating crack propagation in anisotropic brittle materials and to compare the performance of a new method based on shape optimization with the phase-field method in numerical calculations.

#### Main Research Content:
1. **Phase-Field Method**: This method is based on the principle of energy minimization and uses a continuous scalar field variable to smooth out sharp discontinuities, thereby simulating diffusive cracks. This method requires a length scale parameter to smooth the discontinuous band, but it results in mesh-sensitive outcomes.
2. **Shape Optimization Method**: This is a method based on shape space optimization on Riemannian manifolds, which minimizes the system energy through a gradient descent algorithm. This method does not require a length scale parameter and can achieve mesh-independent results.

#### Research Objectives:
1. **Extend Gradient-Based Methods**: Extend gradient-based methods to anisotropic material behavior and compare them with the phase-field method.
2. **Numerical Tests**: Validate the effectiveness of the shape optimization method through numerical calculations on two boundary value problems (isotropic and anisotropic materials) and analyze the impact of anisotropy on crack paths.

### Summary
The main purpose of the paper is to propose a new shape optimization-based method to simulate crack propagation in anisotropic brittle materials and to validate its effectiveness through numerical experiments. Additionally, by comparing the new method with the traditional phase-field method, the paper demonstrates the performance of both methods under different material behaviors.