General constitutive updating for finite strain formulations based on assumed strains and the Jacobian

P. Areias, R. Natal Jorge, J. César de Sá, Timon Rabczuk, J. Reinoso

Research output: Contribution to journalArticle

2 Citations (Scopus)

Abstract

Compatibility between element technology featuring assumed (finite)-strains based on least-squares and current constitutive formulations employed in elastic and inelastic contexts is a demanding task. Local frames are required for anisotropic and cohesive laws, some assumed-strain element technologies do not explicitly provide the deformation gradient, and total Lagrangian approaches are often inadequate for advanced plasticity models. Kirchhoff stress-based Fe Fp decompositions are also not convenient for ductile damage models. In addition, if rotational degrees-of-freedom are used, as is the case in beams and shells, the adoption of a fixed undeformed configuration causes implementation brittleness. An additional aspect to consider is remeshing by element partitioning, which precludes the storage of constitutive tensors in local frames, invalidating the stored quantities. Based on seven algorithmic requirements and the corresponding design solutions, we introduce a general constitutive updating algorithm based on the strain and the Jacobian provided by the element. This allows the use of virtually any constitutive law with any finite-strain element formulation while satisfying the seven requirements. In addition, Newton-Raphson convergence properties are extraordinary, at the cost of precision in the strain rate estimation. As a prototype element implementation, we present a stable hexahedron based on least-squares strains. A BFGS secant estimation is employed for the weight in the least-squares so that softening constitutive laws can be adopted without stability issues at the element level.

Original languageEnglish
Pages (from-to)32-45
Number of pages14
JournalFinite Elements in Analysis and Design
Volume143
DOIs
Publication statusPublished - 2018 Apr 1
Externally publishedYes

Fingerprint

Finite Strain
Updating
Least Squares
Constitutive Law
Formulation
Hexahedron
Remeshing
Newton-Raphson
Requirements
Softening
Strain Rate
Chord or secant line
Plasticity
Convergence Properties
Compatibility
Partitioning
Shell
Tensor
Damage
Degree of freedom

Keywords

  • Anisotropy
  • Constitutive integration
  • Finite strains
  • Local frame
  • Löwdin frame
  • Softening

ASJC Scopus subject areas

  • Analysis
  • Engineering(all)
  • Computer Graphics and Computer-Aided Design
  • Applied Mathematics

Cite this

General constitutive updating for finite strain formulations based on assumed strains and the Jacobian. / Areias, P.; Natal Jorge, R.; César de Sá, J.; Rabczuk, Timon; Reinoso, J.

In: Finite Elements in Analysis and Design, Vol. 143, 01.04.2018, p. 32-45.

Research output: Contribution to journalArticle

Areias, P. ; Natal Jorge, R. ; César de Sá, J. ; Rabczuk, Timon ; Reinoso, J. / General constitutive updating for finite strain formulations based on assumed strains and the Jacobian. In: Finite Elements in Analysis and Design. 2018 ; Vol. 143. pp. 32-45.
@article{20504f030d394dfeb965b9b24fbd6951,
title = "General constitutive updating for finite strain formulations based on assumed strains and the Jacobian",
abstract = "Compatibility between element technology featuring assumed (finite)-strains based on least-squares and current constitutive formulations employed in elastic and inelastic contexts is a demanding task. Local frames are required for anisotropic and cohesive laws, some assumed-strain element technologies do not explicitly provide the deformation gradient, and total Lagrangian approaches are often inadequate for advanced plasticity models. Kirchhoff stress-based Fe Fp decompositions are also not convenient for ductile damage models. In addition, if rotational degrees-of-freedom are used, as is the case in beams and shells, the adoption of a fixed undeformed configuration causes implementation brittleness. An additional aspect to consider is remeshing by element partitioning, which precludes the storage of constitutive tensors in local frames, invalidating the stored quantities. Based on seven algorithmic requirements and the corresponding design solutions, we introduce a general constitutive updating algorithm based on the strain and the Jacobian provided by the element. This allows the use of virtually any constitutive law with any finite-strain element formulation while satisfying the seven requirements. In addition, Newton-Raphson convergence properties are extraordinary, at the cost of precision in the strain rate estimation. As a prototype element implementation, we present a stable hexahedron based on least-squares strains. A BFGS secant estimation is employed for the weight in the least-squares so that softening constitutive laws can be adopted without stability issues at the element level.",
keywords = "Anisotropy, Constitutive integration, Finite strains, Local frame, L{\"o}wdin frame, Softening",
author = "P. Areias and {Natal Jorge}, R. and {C{\'e}sar de S{\'a}}, J. and Timon Rabczuk and J. Reinoso",
year = "2018",
month = "4",
day = "1",
doi = "10.1016/j.finel.2018.01.006",
language = "English",
volume = "143",
pages = "32--45",
journal = "Finite Elements in Analysis and Design",
issn = "0168-874X",
publisher = "Elsevier",

}

TY - JOUR

T1 - General constitutive updating for finite strain formulations based on assumed strains and the Jacobian

AU - Areias, P.

AU - Natal Jorge, R.

AU - César de Sá, J.

AU - Rabczuk, Timon

AU - Reinoso, J.

PY - 2018/4/1

Y1 - 2018/4/1

N2 - Compatibility between element technology featuring assumed (finite)-strains based on least-squares and current constitutive formulations employed in elastic and inelastic contexts is a demanding task. Local frames are required for anisotropic and cohesive laws, some assumed-strain element technologies do not explicitly provide the deformation gradient, and total Lagrangian approaches are often inadequate for advanced plasticity models. Kirchhoff stress-based Fe Fp decompositions are also not convenient for ductile damage models. In addition, if rotational degrees-of-freedom are used, as is the case in beams and shells, the adoption of a fixed undeformed configuration causes implementation brittleness. An additional aspect to consider is remeshing by element partitioning, which precludes the storage of constitutive tensors in local frames, invalidating the stored quantities. Based on seven algorithmic requirements and the corresponding design solutions, we introduce a general constitutive updating algorithm based on the strain and the Jacobian provided by the element. This allows the use of virtually any constitutive law with any finite-strain element formulation while satisfying the seven requirements. In addition, Newton-Raphson convergence properties are extraordinary, at the cost of precision in the strain rate estimation. As a prototype element implementation, we present a stable hexahedron based on least-squares strains. A BFGS secant estimation is employed for the weight in the least-squares so that softening constitutive laws can be adopted without stability issues at the element level.

AB - Compatibility between element technology featuring assumed (finite)-strains based on least-squares and current constitutive formulations employed in elastic and inelastic contexts is a demanding task. Local frames are required for anisotropic and cohesive laws, some assumed-strain element technologies do not explicitly provide the deformation gradient, and total Lagrangian approaches are often inadequate for advanced plasticity models. Kirchhoff stress-based Fe Fp decompositions are also not convenient for ductile damage models. In addition, if rotational degrees-of-freedom are used, as is the case in beams and shells, the adoption of a fixed undeformed configuration causes implementation brittleness. An additional aspect to consider is remeshing by element partitioning, which precludes the storage of constitutive tensors in local frames, invalidating the stored quantities. Based on seven algorithmic requirements and the corresponding design solutions, we introduce a general constitutive updating algorithm based on the strain and the Jacobian provided by the element. This allows the use of virtually any constitutive law with any finite-strain element formulation while satisfying the seven requirements. In addition, Newton-Raphson convergence properties are extraordinary, at the cost of precision in the strain rate estimation. As a prototype element implementation, we present a stable hexahedron based on least-squares strains. A BFGS secant estimation is employed for the weight in the least-squares so that softening constitutive laws can be adopted without stability issues at the element level.

KW - Anisotropy

KW - Constitutive integration

KW - Finite strains

KW - Local frame

KW - Löwdin frame

KW - Softening

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

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

U2 - 10.1016/j.finel.2018.01.006

DO - 10.1016/j.finel.2018.01.006

M3 - Article

AN - SCOPUS:85042180892

VL - 143

SP - 32

EP - 45

JO - Finite Elements in Analysis and Design

JF - Finite Elements in Analysis and Design

SN - 0168-874X

ER -