Element-wise fracture algorithm based on rotation of edges

P. Areias, Timon Rabczuk, D. Dias-da-Costa

Research output: Contribution to journalArticle

152 Citations (Scopus)

Abstract

We propose an alternative, simpler algorithm for FEM-based computational fracture in brittle, quasi-brittle and ductile materials based on edge rotations. Rotation axes are the crack front edges (respectively nodes in surface discretizations) and each rotated edge affects the position of only one or two nodes. Modified positions of the entities minimize the difference between the predicted crack path (which depends on the specific propagation theory in use) and the edge or face orientation. The construction of all many-to-many relations between geometrical entities in a finite element code motivates operations on existing entities retaining most of the relations, in contrast with remeshing (even tip remeshing) and enrichment which alter the structure of the relations and introduce additional entities to the relation graph (in the case of XFEM, enriched elements which can be significantly different than classical FEM elements and still pose challenges for ductile fracture or large amplitude sliding). In this sense, the proposed solution has algorithmic and generality advantages. The propagation algorithm is simpler than the aforementioned alternatives and the approach is independent of the underlying element used for discretization. For history-dependent materials, there are still some transfer of relevant quantities between meshes. However, diffusion of results is more limited than with tip or full remeshing. To illustrate the advantages of our approach, two prototype models are used: tip energy dissipation (LEFM) and cohesive-zone approaches. The Sutton crack path criterion is employed. Traditional fracture benchmarks and newly proposed verification tests are solved. These were found to be very good in terms of crack path and load/deflection accuracy.

Original languageEnglish
Pages (from-to)113-137
Number of pages25
JournalEngineering Fracture Mechanics
Volume110
DOIs
Publication statusPublished - 2013 Sep 1

Fingerprint

Cracks
Finite element method
Ductile fracture
Energy dissipation

Keywords

  • Ductile materials
  • Edge rotations
  • Finite strains
  • Fracture
  • Quasi-brittle materials

ASJC Scopus subject areas

  • Mechanical Engineering
  • Mechanics of Materials
  • Materials Science(all)

Cite this

Element-wise fracture algorithm based on rotation of edges. / Areias, P.; Rabczuk, Timon; Dias-da-Costa, D.

In: Engineering Fracture Mechanics, Vol. 110, 01.09.2013, p. 113-137.

Research output: Contribution to journalArticle

Areias, P. ; Rabczuk, Timon ; Dias-da-Costa, D. / Element-wise fracture algorithm based on rotation of edges. In: Engineering Fracture Mechanics. 2013 ; Vol. 110. pp. 113-137.
@article{fc875bf9fcd047e3b2986893a00153d1,
title = "Element-wise fracture algorithm based on rotation of edges",
abstract = "We propose an alternative, simpler algorithm for FEM-based computational fracture in brittle, quasi-brittle and ductile materials based on edge rotations. Rotation axes are the crack front edges (respectively nodes in surface discretizations) and each rotated edge affects the position of only one or two nodes. Modified positions of the entities minimize the difference between the predicted crack path (which depends on the specific propagation theory in use) and the edge or face orientation. The construction of all many-to-many relations between geometrical entities in a finite element code motivates operations on existing entities retaining most of the relations, in contrast with remeshing (even tip remeshing) and enrichment which alter the structure of the relations and introduce additional entities to the relation graph (in the case of XFEM, enriched elements which can be significantly different than classical FEM elements and still pose challenges for ductile fracture or large amplitude sliding). In this sense, the proposed solution has algorithmic and generality advantages. The propagation algorithm is simpler than the aforementioned alternatives and the approach is independent of the underlying element used for discretization. For history-dependent materials, there are still some transfer of relevant quantities between meshes. However, diffusion of results is more limited than with tip or full remeshing. To illustrate the advantages of our approach, two prototype models are used: tip energy dissipation (LEFM) and cohesive-zone approaches. The Sutton crack path criterion is employed. Traditional fracture benchmarks and newly proposed verification tests are solved. These were found to be very good in terms of crack path and load/deflection accuracy.",
keywords = "Ductile materials, Edge rotations, Finite strains, Fracture, Quasi-brittle materials",
author = "P. Areias and Timon Rabczuk and D. Dias-da-Costa",
year = "2013",
month = "9",
day = "1",
doi = "10.1016/j.engfracmech.2013.06.006",
language = "English",
volume = "110",
pages = "113--137",
journal = "Engineering Fracture Mechanics",
issn = "0013-7944",
publisher = "Elsevier BV",

}

TY - JOUR

T1 - Element-wise fracture algorithm based on rotation of edges

AU - Areias, P.

AU - Rabczuk, Timon

AU - Dias-da-Costa, D.

PY - 2013/9/1

Y1 - 2013/9/1

N2 - We propose an alternative, simpler algorithm for FEM-based computational fracture in brittle, quasi-brittle and ductile materials based on edge rotations. Rotation axes are the crack front edges (respectively nodes in surface discretizations) and each rotated edge affects the position of only one or two nodes. Modified positions of the entities minimize the difference between the predicted crack path (which depends on the specific propagation theory in use) and the edge or face orientation. The construction of all many-to-many relations between geometrical entities in a finite element code motivates operations on existing entities retaining most of the relations, in contrast with remeshing (even tip remeshing) and enrichment which alter the structure of the relations and introduce additional entities to the relation graph (in the case of XFEM, enriched elements which can be significantly different than classical FEM elements and still pose challenges for ductile fracture or large amplitude sliding). In this sense, the proposed solution has algorithmic and generality advantages. The propagation algorithm is simpler than the aforementioned alternatives and the approach is independent of the underlying element used for discretization. For history-dependent materials, there are still some transfer of relevant quantities between meshes. However, diffusion of results is more limited than with tip or full remeshing. To illustrate the advantages of our approach, two prototype models are used: tip energy dissipation (LEFM) and cohesive-zone approaches. The Sutton crack path criterion is employed. Traditional fracture benchmarks and newly proposed verification tests are solved. These were found to be very good in terms of crack path and load/deflection accuracy.

AB - We propose an alternative, simpler algorithm for FEM-based computational fracture in brittle, quasi-brittle and ductile materials based on edge rotations. Rotation axes are the crack front edges (respectively nodes in surface discretizations) and each rotated edge affects the position of only one or two nodes. Modified positions of the entities minimize the difference between the predicted crack path (which depends on the specific propagation theory in use) and the edge or face orientation. The construction of all many-to-many relations between geometrical entities in a finite element code motivates operations on existing entities retaining most of the relations, in contrast with remeshing (even tip remeshing) and enrichment which alter the structure of the relations and introduce additional entities to the relation graph (in the case of XFEM, enriched elements which can be significantly different than classical FEM elements and still pose challenges for ductile fracture or large amplitude sliding). In this sense, the proposed solution has algorithmic and generality advantages. The propagation algorithm is simpler than the aforementioned alternatives and the approach is independent of the underlying element used for discretization. For history-dependent materials, there are still some transfer of relevant quantities between meshes. However, diffusion of results is more limited than with tip or full remeshing. To illustrate the advantages of our approach, two prototype models are used: tip energy dissipation (LEFM) and cohesive-zone approaches. The Sutton crack path criterion is employed. Traditional fracture benchmarks and newly proposed verification tests are solved. These were found to be very good in terms of crack path and load/deflection accuracy.

KW - Ductile materials

KW - Edge rotations

KW - Finite strains

KW - Fracture

KW - Quasi-brittle materials

UR - http://www.scopus.com/inward/record.url?scp=84885862548&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84885862548&partnerID=8YFLogxK

U2 - 10.1016/j.engfracmech.2013.06.006

DO - 10.1016/j.engfracmech.2013.06.006

M3 - Article

AN - SCOPUS:84885862548

VL - 110

SP - 113

EP - 137

JO - Engineering Fracture Mechanics

JF - Engineering Fracture Mechanics

SN - 0013-7944

ER -