Explicit Runge–Kutta methods for the integration of rate-type constitutive equations

Modern constitutive models have the potential to improve the quality of engineering calculations involving non-linear anisotropic materials. The adoption of complex models in practice, however, depends on the availability of reliable and accurate solution methods for the stress point integration problem. This paper presents a modular implementation of explicit Runge–Kutta methods with error control, that is suitable for use, without change, with any rate-type constitutive model. The paper also shows how the complications caused by the algebraic constraint of conventional plasticity are resolved through a simple subloading modification. With this modification any rate-independent model can be implemented without difficulty, using the integration module as an accurate and robust standard procedure. The effectiveness and efficiency of the method are demonstrated through a comparative evaluation of second and fifth-order formulas, applied to a complex constitutive model for natural clay, full details of which are given.

[1]  J. Dormand,et al.  A family of embedded Runge-Kutta formulae , 1980 .

[2]  D. Muir Wood,et al.  A new approach to anisotropic, bounding surface plasticity: general formulation and simulations of natural and reconstituted clay behaviour , 2001 .

[3]  Lidija Zdravković,et al.  General Formulation of Two Kinematic Hardening Constitutive Models with a Smooth Elastoplastic Transition , 2006 .

[4]  Yannis F. Dafalias,et al.  BOUNDING SURFACE PLASTICITY, I: MATHEMATICAL FOUNDATION AND HYPOPLASTICITY , 1986 .

[5]  O. Scherf,et al.  A BDF2 integration method with step size control for elasto-plasticity , 2004 .

[6]  Koichi Hashiguchi,et al.  Constitutive Equations of Elastoplastic Materials With Anisotropic Hardening and Elastic-Plastic Transition , 1981 .

[7]  E. Hairer,et al.  Solving Ordinary Differential Equations II , 2010 .

[8]  Michael Ortiz,et al.  An analysis of a new class of integration algorithms for elastoplastic constitutive relations , 1986 .

[9]  Takashi Okayasu,et al.  Evaluation of typical conventional and unconventional plasticity models for prediction of softening behaviour of soils , 2002 .

[10]  Bernd Simeon,et al.  Runge---Kutta methods in elastoplasticity , 2002 .

[11]  Mohamed Rouainia,et al.  Implicit numerical integration for a kinematic hardening soil plasticity model , 2001 .

[12]  Mohamed Rouainia,et al.  An implicit constitutive algorithm for finite strain Cam-clay elasto-plastic model , 2000 .

[13]  L. Shampine,et al.  Some practical Runge-Kutta formulas , 1986 .

[14]  Otto T. Bruhns,et al.  On the simultaneous estimation of model parameters used in constitutive laws for inelastic material behaviour , 1998 .

[15]  Mohamed Rouainia,et al.  A kinematic hardening constitutive model for natural clays with loss of structure , 2000 .

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

[17]  Jerrold L. Wagener,et al.  Fortran 95 Handbook , 1997 .

[18]  Koichi Hashiguchi,et al.  DESCRIPTION OF INHERENT/INDUCED ANISOTROPY OF SOILS : ROTATIONAL HARDENING RULE WITH OBJECTIVITY , 2001 .

[19]  William H. Press,et al.  Numerical Recipes: FORTRAN , 1988 .

[20]  D. L. Parnas,et al.  On the criteria to be used in decomposing systems into modules , 1972, Software Pioneers.

[21]  Koichi Hashiguchi.  Constitutive Equations of Elastoplastic Materials With Elastic-Plastic Transition , 1980 .

[22]  L. Shampine,et al.  Fundamentals of Numerical Computing , 1997 .

[23]  E. Hairer,et al.  Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems , 2010 .

[24]  Robert L. Taylor,et al.  Finite element implementation of non-linear elastoplastic constitutive laws using local and global explicit algorithms with automatic error control† , 2001 .