Element-wise algorithm for modeling ductile fracture with the Rousselier yield function

P. Areias, D. Dias-Da-Costa, J. M. Sargado, Timon Rabczuk

Research output: Contribution to journalArticle

34 Citations (Scopus)

Abstract

Within the theme of ductile fracture in metals, we propose an algorithm for FEM-based computational fracture based on edge rotations and smoothing of complementarity conditions. Rotation axes are the crack front nodes in surface discretizations and each rotated edge affects the position of only one or two nodes. Modified edge positions correspond to the predicted crack path. To represent softening, porous plasticity in the form of the Rousselier yield function is used. The finite strain integration algorithm makes use of a consistent updated Lagrangian formulation which makes use of polar decomposition between each increment. Constitutive updating is based on the implicit integration of a regularized non-smooth problem. The proposed alternative is advantageous when compared with enriched elements that can be significantly different than classical FEM elements and still pose challenges for ductile fracture or large amplitude sliding. 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, fracture examples making use of the Rousselier yield function are presented. The Ma-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/displacement accuracy.

Original languageEnglish
Pages (from-to)1429-1443
Number of pages15
JournalComputational Mechanics
Volume52
Issue number6
DOIs
Publication statusPublished - 2013 Dec 1
Externally publishedYes

Fingerprint

Ductile Fracture
Ductile fracture
Crack
Cracks
Modeling
Path
Finite element method
Polar decomposition
Finite Strain
Remeshing
Softening
Complementarity
Vertex of a graph
Plasticity
Increment
Updating
Smoothing
Metals
Discretization
Mesh

Keywords

  • Ductile
  • Fracture
  • Plasticity
  • Porous

ASJC Scopus subject areas

  • Computational Theory and Mathematics
  • Mechanical Engineering
  • Ocean Engineering
  • Applied Mathematics
  • Computational Mathematics

Cite this

Element-wise algorithm for modeling ductile fracture with the Rousselier yield function. / Areias, P.; Dias-Da-Costa, D.; Sargado, J. M.; Rabczuk, Timon.

In: Computational Mechanics, Vol. 52, No. 6, 01.12.2013, p. 1429-1443.

Research output: Contribution to journalArticle

Areias, P. ; Dias-Da-Costa, D. ; Sargado, J. M. ; Rabczuk, Timon. / Element-wise algorithm for modeling ductile fracture with the Rousselier yield function. In: Computational Mechanics. 2013 ; Vol. 52, No. 6. pp. 1429-1443.
@article{0b6a61de4e6a48deabfd606d5b3416a8,
title = "Element-wise algorithm for modeling ductile fracture with the Rousselier yield function",
abstract = "Within the theme of ductile fracture in metals, we propose an algorithm for FEM-based computational fracture based on edge rotations and smoothing of complementarity conditions. Rotation axes are the crack front nodes in surface discretizations and each rotated edge affects the position of only one or two nodes. Modified edge positions correspond to the predicted crack path. To represent softening, porous plasticity in the form of the Rousselier yield function is used. The finite strain integration algorithm makes use of a consistent updated Lagrangian formulation which makes use of polar decomposition between each increment. Constitutive updating is based on the implicit integration of a regularized non-smooth problem. The proposed alternative is advantageous when compared with enriched elements that can be significantly different than classical FEM elements and still pose challenges for ductile fracture or large amplitude sliding. 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, fracture examples making use of the Rousselier yield function are presented. The Ma-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/displacement accuracy.",
keywords = "Ductile, Fracture, Plasticity, Porous",
author = "P. Areias and D. Dias-Da-Costa and Sargado, {J. M.} and Timon Rabczuk",
year = "2013",
month = "12",
day = "1",
doi = "10.1007/s00466-013-0885-0",
language = "English",
volume = "52",
pages = "1429--1443",
journal = "Computational Mechanics",
issn = "0178-7675",
publisher = "Springer Verlag",
number = "6",

}

TY - JOUR

T1 - Element-wise algorithm for modeling ductile fracture with the Rousselier yield function

AU - Areias, P.

AU - Dias-Da-Costa, D.

AU - Sargado, J. M.

AU - Rabczuk, Timon

PY - 2013/12/1

Y1 - 2013/12/1

N2 - Within the theme of ductile fracture in metals, we propose an algorithm for FEM-based computational fracture based on edge rotations and smoothing of complementarity conditions. Rotation axes are the crack front nodes in surface discretizations and each rotated edge affects the position of only one or two nodes. Modified edge positions correspond to the predicted crack path. To represent softening, porous plasticity in the form of the Rousselier yield function is used. The finite strain integration algorithm makes use of a consistent updated Lagrangian formulation which makes use of polar decomposition between each increment. Constitutive updating is based on the implicit integration of a regularized non-smooth problem. The proposed alternative is advantageous when compared with enriched elements that can be significantly different than classical FEM elements and still pose challenges for ductile fracture or large amplitude sliding. 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, fracture examples making use of the Rousselier yield function are presented. The Ma-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/displacement accuracy.

AB - Within the theme of ductile fracture in metals, we propose an algorithm for FEM-based computational fracture based on edge rotations and smoothing of complementarity conditions. Rotation axes are the crack front nodes in surface discretizations and each rotated edge affects the position of only one or two nodes. Modified edge positions correspond to the predicted crack path. To represent softening, porous plasticity in the form of the Rousselier yield function is used. The finite strain integration algorithm makes use of a consistent updated Lagrangian formulation which makes use of polar decomposition between each increment. Constitutive updating is based on the implicit integration of a regularized non-smooth problem. The proposed alternative is advantageous when compared with enriched elements that can be significantly different than classical FEM elements and still pose challenges for ductile fracture or large amplitude sliding. 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, fracture examples making use of the Rousselier yield function are presented. The Ma-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/displacement accuracy.

KW - Ductile

KW - Fracture

KW - Plasticity

KW - Porous

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

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

U2 - 10.1007/s00466-013-0885-0

DO - 10.1007/s00466-013-0885-0

M3 - Article

AN - SCOPUS:84889881976

VL - 52

SP - 1429

EP - 1443

JO - Computational Mechanics

JF - Computational Mechanics

SN - 0178-7675

IS - 6

ER -