Energy consistent framework for continuously evolving 3D crack propagation

Łukasz Kaczmarczyk, Zahur Ullah, Chris J. Pearce

Research output: Contribution to journalArticle

7 Citations (Scopus)

Abstract

This paper presents an enhanced theoretical formulation and associated computational framework for brittle fracture in elastic solids within the context of configurational mechanics, building on the authors’ previous paper, Kaczmarczyk et al. (2014). The local form of the first law of thermodynamics provides an equilibrium condition for the crack front, expressed in terms of the configurational forces. Applying the principle of maximal energy dissipation, it is shown that the direction of the crack propagation is given by the direction of the configurational forces. In combination with a fracture criterion, these are utilised to determine the position of the continuously evolving crack front. This exploitation of the crack front equilibrium condition leads to a completely new, implicit, crack propagation formulation. A monolithic solution strategy is adopted, solving simultaneously for both the material displacements (i.e. crack extension) and the spatial displacements. The resulting crack path is resolved as a discrete displacement discontinuity, where the material displacements of the nodes on the crack front change continuously, without the need for mesh splitting or the use of enrichment techniques. In order to trace the dissipative loading path, an arc-length procedure is adopted that controls the incremental crack area growth. In order to maintain mesh quality, smoothing of the mesh is undertaken as a continuous process, together with face flipping, node merging and edge splitting where necessary. Hierarchical basis functions of arbitrary polynomial order are adopted to increase the order of approximation without the need to change the finite element mesh. Performance of the formulation is demonstrated by means of three representative numerical simulations, demonstrating both accuracy and robustness.
LanguageEnglish
Pages54-73
JournalComputer Methods in Applied Mechanics and Engineering
Volume324
DOIs
Publication statusPublished - 15 Jun 2017

Fingerprint

Crack propagation
Cracks
Brittle fracture
Merging
Energy dissipation
Mechanics
Polynomials
Thermodynamics
Computer simulation

Keywords

  • Brittle fracture
  • Configurational forces
  • Crack tracking
  • Mesh smoothing

Cite this

Kaczmarczyk, Łukasz ; Ullah, Zahur ; Pearce, Chris J. / Energy consistent framework for continuously evolving 3D crack propagation. In: Computer Methods in Applied Mechanics and Engineering. 2017 ; Vol. 324. pp. 54-73.
@article{83dafcb5005e41c980f9d396cd64f3d8,
title = "Energy consistent framework for continuously evolving 3D crack propagation",
abstract = "This paper presents an enhanced theoretical formulation and associated computational framework for brittle fracture in elastic solids within the context of configurational mechanics, building on the authors’ previous paper, Kaczmarczyk et al. (2014). The local form of the first law of thermodynamics provides an equilibrium condition for the crack front, expressed in terms of the configurational forces. Applying the principle of maximal energy dissipation, it is shown that the direction of the crack propagation is given by the direction of the configurational forces. In combination with a fracture criterion, these are utilised to determine the position of the continuously evolving crack front. This exploitation of the crack front equilibrium condition leads to a completely new, implicit, crack propagation formulation. A monolithic solution strategy is adopted, solving simultaneously for both the material displacements (i.e. crack extension) and the spatial displacements. The resulting crack path is resolved as a discrete displacement discontinuity, where the material displacements of the nodes on the crack front change continuously, without the need for mesh splitting or the use of enrichment techniques. In order to trace the dissipative loading path, an arc-length procedure is adopted that controls the incremental crack area growth. In order to maintain mesh quality, smoothing of the mesh is undertaken as a continuous process, together with face flipping, node merging and edge splitting where necessary. Hierarchical basis functions of arbitrary polynomial order are adopted to increase the order of approximation without the need to change the finite element mesh. Performance of the formulation is demonstrated by means of three representative numerical simulations, demonstrating both accuracy and robustness.",
keywords = "Brittle fracture, Configurational forces, Crack tracking, Mesh smoothing",
author = "Łukasz Kaczmarczyk and Zahur Ullah and Pearce, {Chris J.}",
year = "2017",
month = "6",
day = "15",
doi = "10.1016/j.cma.2017.06.001",
language = "English",
volume = "324",
pages = "54--73",
journal = "Computer Methods in Applied Mechanics and Engineering",
issn = "0045-7825",
publisher = "Elsevier",

}

Energy consistent framework for continuously evolving 3D crack propagation. / Kaczmarczyk, Łukasz; Ullah, Zahur; Pearce, Chris J.

In: Computer Methods in Applied Mechanics and Engineering, Vol. 324, 15.06.2017, p. 54-73.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Energy consistent framework for continuously evolving 3D crack propagation

AU - Kaczmarczyk, Łukasz

AU - Ullah, Zahur

AU - Pearce, Chris J.

PY - 2017/6/15

Y1 - 2017/6/15

N2 - This paper presents an enhanced theoretical formulation and associated computational framework for brittle fracture in elastic solids within the context of configurational mechanics, building on the authors’ previous paper, Kaczmarczyk et al. (2014). The local form of the first law of thermodynamics provides an equilibrium condition for the crack front, expressed in terms of the configurational forces. Applying the principle of maximal energy dissipation, it is shown that the direction of the crack propagation is given by the direction of the configurational forces. In combination with a fracture criterion, these are utilised to determine the position of the continuously evolving crack front. This exploitation of the crack front equilibrium condition leads to a completely new, implicit, crack propagation formulation. A monolithic solution strategy is adopted, solving simultaneously for both the material displacements (i.e. crack extension) and the spatial displacements. The resulting crack path is resolved as a discrete displacement discontinuity, where the material displacements of the nodes on the crack front change continuously, without the need for mesh splitting or the use of enrichment techniques. In order to trace the dissipative loading path, an arc-length procedure is adopted that controls the incremental crack area growth. In order to maintain mesh quality, smoothing of the mesh is undertaken as a continuous process, together with face flipping, node merging and edge splitting where necessary. Hierarchical basis functions of arbitrary polynomial order are adopted to increase the order of approximation without the need to change the finite element mesh. Performance of the formulation is demonstrated by means of three representative numerical simulations, demonstrating both accuracy and robustness.

AB - This paper presents an enhanced theoretical formulation and associated computational framework for brittle fracture in elastic solids within the context of configurational mechanics, building on the authors’ previous paper, Kaczmarczyk et al. (2014). The local form of the first law of thermodynamics provides an equilibrium condition for the crack front, expressed in terms of the configurational forces. Applying the principle of maximal energy dissipation, it is shown that the direction of the crack propagation is given by the direction of the configurational forces. In combination with a fracture criterion, these are utilised to determine the position of the continuously evolving crack front. This exploitation of the crack front equilibrium condition leads to a completely new, implicit, crack propagation formulation. A monolithic solution strategy is adopted, solving simultaneously for both the material displacements (i.e. crack extension) and the spatial displacements. The resulting crack path is resolved as a discrete displacement discontinuity, where the material displacements of the nodes on the crack front change continuously, without the need for mesh splitting or the use of enrichment techniques. In order to trace the dissipative loading path, an arc-length procedure is adopted that controls the incremental crack area growth. In order to maintain mesh quality, smoothing of the mesh is undertaken as a continuous process, together with face flipping, node merging and edge splitting where necessary. Hierarchical basis functions of arbitrary polynomial order are adopted to increase the order of approximation without the need to change the finite element mesh. Performance of the formulation is demonstrated by means of three representative numerical simulations, demonstrating both accuracy and robustness.

KW - Brittle fracture

KW - Configurational forces

KW - Crack tracking

KW - Mesh smoothing

U2 - 10.1016/j.cma.2017.06.001

DO - 10.1016/j.cma.2017.06.001

M3 - Article

VL - 324

SP - 54

EP - 73

JO - Computer Methods in Applied Mechanics and Engineering

T2 - Computer Methods in Applied Mechanics and Engineering

JF - Computer Methods in Applied Mechanics and Engineering

SN - 0045-7825

ER -