A stress-strain integration algorithm for unsaturated soil elastoplasticity with automatic error control

Stress-strain integration algorithms are a very important component in the development of finite element codes. The use of accurate, robust and fast stress-strain integration algorithms accounts for a significant part of the performance of a finite element code, especially when complex elasto-plastic constitutive models are used. This paper presents the formulation of an algorithm for the stress-strain integration of the Barcelona Basic Model, an elasto-plastic volumetric hardening constitutive model for unsaturated soils. The proposed algorithm is based on some earlier ideas of substepping explicit integration with automatic error control. Stress-strain integration in the plastic domain is performed by using an explicit algorithm which accuracy is improved by dividing the initial strain increment in a number of substeps. The number of substeps depends on an estimation of the integration error that is obtained by using a modified Euler procedure. The performance of the stress-strain algorithm is presented for different types of stress paths to assess its dependency on factors such as, for example, the initial size of the strain increment, the error tolerance and the initial stress state. A drift correction algorithm is also proposed and its influence on the results is evaluated. 2 INTEGRATION OF ELASTO-PLASTIC STRAIN

[1]  Richard A. Regueiro,et al.  Implicit numerical integration of a three-invariant, isotropic/kinematic hardening cap plasticity model for geomaterials , 2005 .

[2]  Ronaldo I. Borja,et al.  Cam-Clay plasticity. Part V: A mathematical framework for three-phase deformation and strain localization analyses of partially saturated porous media , 2004 .

[3]  Lidija Zdravković,et al.  Finite element analysis in geotechnical engineering , 1999 .

[4]  Scott W. Sloan,et al.  Aspects of finite element implementation of critical state models , 2000 .

[5]  S. Sloan,et al.  Refined explicit integration of elastoplastic models with automatic error control , 2001 .

[6]  Antonio Gens,et al.  Finite element formulation and algorithms for unsaturated soils. Part I: Theory , 2003 .

[7]  Antonio Gens,et al.  Finite element formulation and algorithms for unsaturated soils. Part II: Verification and application , 2003 .

[8]  Ronaldo I. Borja,et al.  Cam-Clay plasticity, Part II: implicit integration of constitutive equation based a nonlinear elastic stress predictor , 1991 .

[9]  Scott W. Sloan,et al.  Substepping schemes for the numerical integration of elastoplastic stress–strain relations , 1987 .

[10]  A. Gens,et al.  An elasto-plastic model for unsaturated soil incorporating the effects of suction and degree of saturation on mechanical behaviour , 2003 .

[11]  Ronaldo I. Borja,et al.  Cam-Clay plasticity, part I: implicit integration of elastoplastic constitutive relations , 1990 .

[12]  Antonio Gens,et al.  A constitutive model for partially saturated soils , 1990 .

[13]  Antonio Gens,et al.  A stress point algorithm for an elastoplastic model in unsaturated soils , 2000 .

[14]  L. Wang,et al.  Formulation of the return mapping algorithm for elastoplastic soil models , 2004 .

[15]  J. M. Watt Numerical Initial Value Problems in Ordinary Differential Equations , 1972 .

[16]  R. de Borst,et al.  Implicit integration of a generalized plasticity constitutive model for partially saturated soil , 2001 .