由于软土具有蠕变特性,而商业化软件提供的岩土本构模型往往不能满足软土地区实际工程数值分析的需求。将一维广义西原模型的应力-应变关系推广为三维应力状态的弹-粘弹-粘塑本构方程,采用D-P屈服准则推导出用于数值计算的粘弹性柔度矩阵和粘塑性柔度矩阵,并成功添加到ADINA有限元程序中。然后利用曲线拟合法确定流变模型参数,并利用二次开发的粘弹塑模型,对软土的三轴蠕变进行数值模拟,通过数值模拟和三轴蠕变试验的结果比较,验证了模型的正确性和可靠性。文中二次开发采用的方法具有普遍性,可以为其他粘弹塑模型的开发提供借鉴和参考。