适用于膨胀性非饱和土的边界面模型的数值实现

基于完全隐式的Euler向后积分算法,对一个适用于膨胀性非饱和土的边界面模型提出一种应力积分方法。模型中的塑性变形具有双重机制即加载–湿陷机制和微观结构机制,并考虑力学行为与持水滞回间的耦合作用。依据当前应力状态、力学与持水模型中各塑性机制的发生情况,将应力积分过程中的当前增量步分为多种情况,针对每一情况利用Euler向后算法进行试算,并依据判定准则选取结果。其次,利用该模型分别对已有的吸力循环试验和三轴压缩试验进行了验证,以此反映模型的预测能力和算法的运行能力。最后,对比不同步长的应变控制试验路径下的计算结果,以此验证应力积分方法的稳定性和精度。