Numerical methods for implementing the bounding surface plasticity model for clays

Previous work on the development and numerical implementation of the Bounding Surface Plasticity model for clays is discussed. Modifications are made to the hardening relationship to improve the numerical performance in the tensile range. A rate equation for the loading surface is developed. Modifications are made to the invariant description of the bounding surface to avoid numerical difficulties in evaluating the derivatives.; The Closest Point Projection method is described for simple and general internal variable plasticity models. The method is developed within the classical plasticity framework and uses the Newton-Raphson method to satisfy the implicit integration of the rate equations and the consistency condition. An explicit treatment of the internal variables is discussed. Application of this method for the Bounding Surface Plasticity model for clays is developed by adding an internal variable and using the rate equation for the loading surface.; A new algorithm coined "the Reduced Newton method" is developed for the Bounding Surface Plasticity model for clays. It involves mapping the stress rate equations, internal variable rate equations and the consistency condition into two nonlinear equations and integrating them with a backwards Euler formula using Newton-Raphson iteration.; Comparison of predictions for a number of sample problems is made using the trapezoidal integration, Closest Point and Reduced Newton methods. "Exact" solutions for stress points that start on the bounding surface are developed by assigning an arbitrary stress path and calculating the corresponding strains using numerical integration with a tight tolerance. These "exact" solutions are used to evaluate the effectiveness of the proposed general numerical implementation of the Bounding Surface model.; A standard effective stress interface is proposed for finite element programs that use the Newton-Raphson method. The Reduced Newton model is implemented within the DYSAC2 finite element program and is used to analyze an earth embankment subjected to earthquake and shock loads.