Implicit and explicit integration schemes in the anisotropic bounding surface plasticity model for cyclic behaviours of saturated clay

Two integration algorithms, namely the implicit return mapping and explicit sub-stepping schemes, are adopted in the anisotropic bounding surface plasticity model for cyclic behaviours of saturated clay and are implemented into finite element code. The model is a representative of a series of bounding surface models that have typical characteristics, including isotropic and kinematic hardening rules and a rotational bounding surface to capture complex but important cyclic behaviours of soils, such as cyclic shakedown and degradation. However, there is no explicit current yield surface in the model to which the conventional implicit algorithm returns the stress state back or the sub-stepping integration corrects the drift of the stress state. Hence, necessary modifications have been made for both of the integration schemes. First, the image stress point is mapped or corrected to the bounding surface instead of mapping back or correcting the stress state to the yield surface. Second, the unloading–loading criterion is checked to determine the image stress point rather than checking the yield criterion after giving the trial stress state in a conventional way. Comparative studies on the accuracy, stability and efficiency of the two integration schemes are conducted not only at the element level but also in solving boundary value problems of monotonic and cyclic bearing behaviours of rigid footings on saturated clay. For smaller strain increments, there is no significant difference in the accuracy between the two integration schemes, but the explicit integration shows a higher efficiency and accuracy. For relatively larger increments, the implicit return mapping algorithm presents good accuracy and more robustness, while the sub-stepping algorithm shows deteriorating accuracy and suffers the convergence problem. With the tolerance used in the present model, the bearing capacity of the rigid footing predicted by the return mapping algorithm is closer to the available analytical and numerical solutions, while the bearing capacity predicted by the sub-stepping algorithm shows a marginal increase.