On the numerical integration of three-invariant elastoplastic constitutive models

We investigate the performance of a numerical algorithm for the integration of isotropically hardening three-invariant elastoplastic constitutive models with convex yield surfaces. The algorithm is based on a spectral representation of stresses and strains for infinitesimal and finite deformation pl...

Full description

Saved in:
Bibliographic Details
Published in:Computer methods in applied mechanics and engineering Vol. 192; no. 9; pp. 1227 - 1258
Main Authors: Borja, Ronaldo I., Sama, Kossi M., Sanz, Pablo F.
Format: Journal Article
Language:English
Published: Amsterdam Elsevier B.V 28.02.2003
Elsevier
Subjects:
Computational techniques
Exact sciences and technology
Finite-element and galerkin methods
Fundamental areas of phenomenology (including applications)
Inelasticity (thermoplasticity, viscoplasticity...)
Mathematical methods in physics
Physics
Solid mechanics
Structural and continuum mechanics
Viscoelasticity, plasticity, viscoplasticity
33
34
46
60
65
Constitutive equation
Yield criterion
Modeling
Soil mechanics
Finite element method
Return mapping
Inelasticity
Boundary value problem
Elastoplasticity
Isotropic hardening
Variational calculus
Principal stress
Strain hardening
ISSN:0045-7825, 1879-2138
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:We investigate the performance of a numerical algorithm for the integration of isotropically hardening three-invariant elastoplastic constitutive models with convex yield surfaces. The algorithm is based on a spectral representation of stresses and strains for infinitesimal and finite deformation plasticity, and a return mapping in principal stress directions. Smooth three-invariant representations of the Mohr–Coulomb model, such as the Lade–Duncan and Matsuoka–Nakai models, are implemented within the framework of the proposed algorithm. Among the specific features incorporated into the formulation are the hardening/softening responses and the tapering of the yield surfaces toward the hydrostatic axis with increasing confining pressure. Isoerror maps are generated to study the local accuracy of the numerical integration algorithm. Finally, a boundary-value problem involving loading of a strip foundation on a soil is analyzed with and without finite deformation effects to investigate the performance of the integration algorithm in a full-scale non-linear finite element simulation.
Bibliography:ObjectType-Article-2
SourceType-Scholarly Journals-1
ObjectType-Feature-1
content type line 23
ISSN:0045-7825
1879-2138
DOI:10.1016/S0045-7825(02)00620-5