A strain-gradient isotropic elastoplastic damage model with J3 dependence

INTERNATIONAL JOURNAL OF SOLIDS AND STRUCTURES

卷:174

Article

Papadioti, I¹  Aravas, N.^1,2  Lian, J.^3,4  Muenstermann, S.⁵ 

[1] Univ Thessaly, Dept Mech Engn, Volos 38334, Greece

[2] Kyushu Univ, Int Inst Carbon Neutral Energy Res WPI I2CNER, Nishi Ku, 744 Moto Oka, Fukuoka, Fukuoka 8190395, Japan

[3] Aalto Univ, Dept Mech Engn, Puumiehenkuja 3, Espoo 02150, Finland

[4] MIT, Dept Mech Engn, Impact & Crashworthiness Lab, 77 Massachusetts Ave, Cambridge, MA 02139 USA

[5] Rhein Westfal TH Aachen, Dept Ferrous Met, Intzestr 1, D-52072 Aachen, Germany

关键词: Strain gradient plasticity;    J(3) dependence;    Damage mechanics;    Finite elements;    Numerical integration of elastoplastic equations;    Hourglass control;   

DOI  :  10.1016/j.ijsolstr.2019.05.027

来源: Elsevier

【 摘 要 】

A plastic strain gradient version of an isotropic elastoplastic damage model that depends on the third invariant J(3) of the stress deviator is developed. The model is based on the non-local equivalent plastic strain e(p) defined by Peerlings et al. (2001) and Engelen et al. (2003) and introduces a material length l to the constitutive equations. It is shown that the non-local equivalent plastic strain e(p) at a material point P can be identified with the average value of the local von Mises equivalent plastic strain (epsilon) over bar (p) over a sphere centered at P and of radius approximately equal to 3 l. A methodology for the numerical integration of the constitutive equations is presented. The algorithm is appropriate for rate-independent as well as rate-dependent (viscoplastic) models. The model is implemented in the ABAQUS general-purpose finite element program and both quasi-static and dynamic problems are solved. Two possible ABAQUS implementations are discussed. First,user elements are developed, which can be used for the solution of both quasi-static and dynamic problems. Reduced 1-point Gauss integration is discussed in 8-node hexahedral elements and the physical stabilization method of Puso (2000) is used to remove the resulting numerical singularities (hourglass control). Second, the implementation of the model via user material subroutines is discussed. Quasi-static problems can be solved with ABAQUS/Standard using a *COUPLED TEMPERATURE-DISPLACEMENT, STEADY STATE analysis together with user subroutine UMAT, in which temperature is identified with the non-local equivalent plastic strain e(p); the solution of dynamic problems requires use of ABAQUS/Explicit together with a *DYNAMIC TEMPERATURE-DISPLACEMENT analysis option and user subroutines VUMAT and DFLUX. Several example problems are solved. (C) 2019 Elsevier Ltd. All rights reserved.

【 授权许可】

Free