abaqus土层地震反应分析中土体模型怎么解决？

我在进行abaqus土层地震反应分析时，建立的土体模型要着重体现土的非线性特征，就是应力随应变是非线性变化的（弹性模量E会折减，应该吧，我也不清楚），有谁知道该怎么弄吗？abaqus中自带的库伦摩尔模型是理想线弹塑性的。有人知道该怎么弄把土体模型弄成非线性的么？

或者有umat材料子程序也可以。

知道的大侠教教我。要是涉及到私人研究成果，有偿求教也行，并保证有关资料只用于本人毕设中。