• Geotechnical Engineering
• /
• v.14 no.4
• /
• pp.129-140
• /
• 1998

Elasto-plastic finite element analysis of geotechnical boundary value problems necessitate the stress integration for the known strain increments. For the elasto-plastic constitutive model, the stress integration is generally achieved by numerical schemes, because analytical integration is impossible for general strain path. In this case, the accuracy of numerical stress integration has an important role on the overall accuracy of nonlinear finite element solution. In this study, the Sloan's substepping method which is one of explicit integration methods has been adopted and iris applicability has been checked. The unstability and inaccuracy of ifs results initiated from initial stress level were revealed. So. a new modified numerical integration method which employs the basic concept of modified Euler scheme for error control is proposed and accuracy and stability of the solutions are confirmed by triaxial test simulation.