Abstract:On the basis of the proposed model which can describe the whole damage progress from hardening to softening for sandy soils and the UMAT materials program interface provided by Abaqus, a program is implemented to embed the proposed model to the Abaqus in which the forward Euler integral algorithm is applied. In this program the measures, such as subbing the step size and limiting drift of the yield surface, are taken to ensure the error to be controlled within the tolerance. The comparison between the numerical simulation results and laboratory experimental results show that the constitutive model can reasonably simulate the whole stress history, accordingly, the effectiveness of the algorithm program is verified.