摘要：对程星磊等提出的总应力形式增量弹塑性本构模型进行二次开发，以模拟复杂应力状态下软黏土的响应。通过Newton-Raphson算法,对材料非线性问题进行迭代求解;针对本构模型中应力反向等一系列关键性问 题，应用欧拉切线算法编写了有限元程序，并结合子增量方法提高了计算精度。预测了软黏土在轴对称应力状 态下的响应，得到了应力应变关系曲线，将其与单元预测结果进行比较,二者趋于一致，从而验证了有限元程序 编写的合理性。利用该程序模拟三轴不固结不排水试验，模拟结果与试验吻合良好，表明该本构模型有限元程 序可以反映轴对称应力状态下软黏土的不排水应力应变特性,可应用于更加复杂边值问题的模拟计算。