关于桩土场初应力

我用平面应变单元模拟桩和土的，而且加了2D面面接触单元。不过在加初应力的时候遇到了问题。我开始
是把桩单元改成土的材料属性，加重力加速度。可是感觉应力有点奇怪，和全部只是土的不一样，也许是
还有接触单元的原因吧，相当于土之间有了结构面。可是把接触单元杀死，合并节点后就不可逆了。不知
道怎么加初始应力呀。求大神解救。