软黏土不排水循环应力应变关系的数值模拟

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