基于能量型砂土的软硬化弹塑性本构模型,结合Abaqus软件提供的用户材料子程序umat接口,提出了调整步长及限制屈服面滑移的向前欧拉显式积分算法,并将其嵌入到Abaqus软件中,实现了本构模型的程序化.数值模拟结果与室内试验结果的比较表明,提出的模型可以较好地模拟砂土加载过程中的应力-应变关系,从而也验证了算法程序的有效性.
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.