通过UMAT实现基于DP屈服准则的改进西原模型的三维粘弹塑性（蠕变）本构模型

在岩石工程中，大量的失稳现象都与岩石的蠕变特性相关。传统西原模型是目前可以比较好地描述岩石蠕变过程曲线的元件模型，但是，西原模型使用的元件为黏弹、黏塑性元件（如图1），难以描述岩石屈服破坏后进入加速阶段的蠕变变形。滑坡预报，特别是临滑预报在地质灾害防治领域具有重要意义。

通过编写abaqus UMAT子程序，可得到如下结果：

（1）应力状态较小时，仅发生弹性应变和粘弹性应变，最后随时间趋于稳定值。

（2）单元屈服时，发生粘弹塑性应变，应变随加载时长逐渐增加，但尚未达到触发应变，曲线呈现两阶段特性。

（3）随着加载时长的增加，应变进一步增加，超越触发应变后，进入快速蠕变阶段，应变快速增加，曲线呈现三阶段蠕变特性。

参考文献：

[1] 齐亚静, 姜清辉, 王志俭, 等. 改进西原模型的三维蠕变本构方程及其参数辨识[J]. 岩石力学与工程学报, 2012, 31(2): 347-355.

[2] 沈才华, 张兵, 王媛, 等. 基于DP屈服准则的西原本构模型及其运用[J]. 地下空间与工程学报, 2016, 12(2): 402-407.

最后，有需要欢迎通过微信公众号联系我们。

微信公众号：320科技工作室。