To describe the behaviours of geomaterials such as time-dependency, anisotropy and destructuration, multiple hardening parameters and laws are generally needed for application in advanced elasto-viscoplastic models. Time integration with stress updating is a key step in the application of elasto-viscoplastic models to engineering practice. However, the robustness of time integration algorithms for such complicated models has rarely been studied, creating difficulties in selecting and improving algorithms. This paper focuses on use of three typical implicit time integration algorithms—Katona, Stolle and cutting plane—for integration of an advanced elasto-viscoplastic model. First, all selected algorithms are improved to fit the characteristics of the advanced model with multiple hardening parameters and are combined with adaptive substepping procedures to enhance their performance. Then a step-changed undrained triaxial test is simulated at the integration point level, on the basis of which variations in iteration number and relative error of stresses with step size are investigated and compared. Furthermore, the advanced model using different algorithms is implemented into finite element code, with global convergence and calculation time investigated and compared for two boundary value problems: a biaxial test and an embankment. All comparisons demonstrate that the modified cutting plane algorithm with substepping is the most robust and efficient one, followed by the modified Stolle with substepping and the modified Katona with substepping, for an advanced model with multiple hardening parameters.

为了描述诸如时间相关性，各向异性和破坏性之类的岩土材料的行为，在先进的弹粘塑性模型中应用通常需要多个硬化参数和定律。将时间与应力更新集成在一起是将弹粘塑性模型应用到工程实践中的关键步骤。但是，很少研究这种复杂模型的时间积分算法的鲁棒性，这给选择和改进算法带来了困难。本文重点介绍了如何使用三种典型的隐式时间积分算法（Katona，Stolle和切割平面）对高级弹塑性粘塑性模型进行积分。第一，所有选定的算法均经过改进，以适应具有多个强化参数的高级模型的特征，并与自适应子步长程序相结合以提高其性能。然后，在积分点水平上模拟了变步不排水三轴试验，在此基础上，研究并比较了迭代次数和应力相对误差随步长的变化。此外，使用不同算法的高级模型被实现为有限元代码，并研究了全局收敛性和计算时间，并比较了两个边界值问题：双轴试验和路堤。所有比较结果都表明，经过改进的带子步距的切割平面算法是最可靠，最有效的算法，