通过GDS应力三轴仪对湛江黏土在恒定围压不同偏应力下的三轴蠕变试验,得到几组不同偏应力下轴向应变与时间关系数据,鉴于传统线性元件组合模型无法精确描述黏土蠕变的非线性问题,提出采用分数阶导数元件分别替换传统Burger模型中两个牛顿粘壶线性元件,建立起含两个分数阶导数元件的非线性蠕变模型,采用C++语言将该模型在FLAC3D中进行二次开发,通过建立一个三轴数值算例,将不同偏应力下数值计算结果与试验数据进行对比,结果表明提出的非线性蠕变模型相比Burger模型更适合描述湛江黏土的蠕变过程。