求ABAQUS接触剪应力极限值设置

霜晨

本人在做饱和粉土中锚的抗拔力仿真，遇到问题，abaqus里定义接触属性中，摩擦力Tangential Behavior采用惩罚系数法，但是这个方法里摩擦力剪应力极限shear stress limit只能设置为定值。然而饱和粉土不是粘土，内摩擦角不可忽略，所以剪应力极限是个随土深度变化而变化的值。
求大神帮忙，我怎么设置才能使剪应力极限值是个变值