An anisotropic cohesive fracture model: advantages and limitations of length-scale insensitive phase-field damage models

TL;DR

This paper develops anisotropic phase-field damage models for fracture that are insensitive to length scale parameters by incorporating direction-dependent fracture energy and strength, enabling efficient simulations of anisotropic materials.

Contribution

It introduces a novel anisotropic cohesive phase-field damage model that accounts for direction-dependent fracture properties and remains insensitive to length scale variations.

Findings

Models are almost insensitive to length scale parameter.

Mesh size can be increased, reducing computational time.

The approach captures anisotropic fracture behavior effectively.

Abstract

The goal of the current work is to explore direction-dependent damage initiation and propagation within an arbitrary anisotropic solid. In particular, we aim at developing anisotropic cohesive phase-field (PF) damage models by extending the idea introduced in \cite{REZAEI2021a} for direction-dependent fracture energy and also anisotropic PF damage models based on structural tensors. The cohesive PF damage formulation used in the current contribution is motivated by the works of \cite{LORENTZ201120, wu2018, GEELEN2019}. The results of the latter models are shown to be insensitive with respect to the length scale parameter for the isotropic case. This is because they manage to formulate the fracture energy as a function of diffuse displacement jumps in the localized damaged zone. In the present paper, we discuss numerical examples and details on finite element implementations where the fracture energy, as well as the material strength, are introduced as an arbitrary function of the crack direction. Using the current formulation for anisotropic cohesive fracture, the obtained results are almost insensitive with respect to the length scale parameter. The latter is achieved by including the direction-dependent strength of the material in addition to its fracture energy. Utilizing the current formulation, one can increase the mesh size which reduces the computational time significantly without any severe change in the predicted crack path and overall obtained load-displacement curves. We also argue that these models still lack to capture mode-dependent fracture properties. Open issues and possible remedies for future developments are finally discussed as well.