In this paper, the semi-implicit return mapping algorithm has been developed to cope the integration of the Lade-Duncan elasto-plastic constitutive model.