• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2007.05a
• /
• pp.272-275
• /
• 2007

This paper is concerned with modified integration algorithm on the strain-space for rate and temperature dependent elasto-plastic constitutive relations in order to obtain more accurate results in numerical implementation. The proposed algorithm is integrated analytically using integration by part and chain rule and then is applied to the 2-stage Lobatto IIIA with second-order accuracy. It has advantage that is able to consider the convective stress rates on the yield surface of the strain-space. Also this paper is carried out the iteration procedure using the Newton-Raphson method to enforce consistency at the end of the step. And the performance of the proposed algorithm for rate and temperature dependent constitutive relation is illustrated by means of analysis of adiabatic shear bands.

• Journal of the Korean Geotechnical Society
• /
• v.22 no.9
• /
• pp.45-59
• /
• 2006

The successful performance of any numerical geotechnical simulation depends on the accuracy and efficiency of the numerical implementation of constitutive model used to simulate the stress-strain (constitutive) response of the soil. The corner stone of the numerical implementation of constitutive models is the numerical integration of the incremental form of soil-plasticity constitutive equations over a discrete sequence of time steps. In this paper a well known two-surface soil plasticity model is implemented using a generalized implicit return mapping algorithm to arbitrary convex yield surfaces referred to as the Closest-Point-Projection method (CPPM). The two-surface model describes the nonlinear behavior of coarse-grained materials by incorporating a bounding surface concept together with isotropic and kinematic hardening as well as fabric formulation to account for the effect of fabric formation on the unloading response. In the course of investigating the performance of the CPPM integration method, it is proven that the algorithm is an accurate, robust, and efficient integration technique useful in finite element contexts. It is also shown that the algorithm produces a consistent tangent operator $\frac{d\sigma}{d\varepsilon}$ during the iterative process with quadratic convergence rate of the global iteration process.