为了使广义西原模型可以描述冻土的各个变形阶段,用非线性牛顿体替代线性牛顿体进行改进,采用类比的方法将冻土单轴应力状态下的本构方程推广到三维状态;在ABAQUS中利用二次开发平台,编写了改进广义西原模型的UMAT子程序,并在单轴、三轴蠕变条件下进行检验。单轴蠕变的数值解与解析解计算结果十分吻合,两淮地区深部冻结粘土三轴蠕变试验模拟值与实验值相符。表明改进的广义西原模型可以很好地描述冻土蠕变变形特征,包括加速蠕变阶段,UMAT子程序可以用于冻结法施工工程数值模拟。