做桩基负摩阻力计算加地应力平衡语句后不收敛

huoshaolei

用hypermesh做的前处理，放到abaqus程序里进行计算，做群桩1/4模型，做好后地应力平衡步没有问题，不加地应力平衡语句（*initial conditions, type=stress, input=qaba.inp ）可以算过去，可是加了地应力平衡语句后就一直算不过去，无法收敛，在附件放了主要的inp文件，因原附件太大，把节点单元的数据都删掉了，请大神帮忙看一下，问题出在哪里了？不尽感激！
模型概述：

1、第一层软土，厚32m；重度17.8 kN/m3，计算时按透水考虑，故取浮重度7.8kN/m3；弹模暂取5倍平均压缩模量，为15MPa，泊松比取0.4；平均内摩擦角10°，黏聚力15kPa；桩土接触面上的摩擦系数取经验值0.20（流塑）。

2、第二层黏土，厚49m；重度19.7 kN/m3，计算时按透水考虑，故取浮重度9.7kN/m3；弹模暂取5倍平均压缩模量，为30MPa，泊松比取0.3；平均内摩擦角20°，黏聚力40kPa；桩土接触面上的摩擦系数取经验值0.25（可塑）。

3、第三层碎石土和强风化岩层，厚9m；重度取22 kN/m3，计算时按透水考虑，故取浮重度12kN/m3；弹模暂取200MPa，泊松比取0.2；平均内摩擦角40°，黏聚力1kPa。

4、再往下为中风化岩层，作为计算边界处理。

5、桩身混凝土，弹性模量3000MPa，泊松比0.2，重度24kN/m3，取浮重度14 kN/m3。

6、顶部回填5m的碎石，重度取22 kN/m3（水上，天然重度），直接按照均布荷载考虑110kPa。

7、取1/4群桩模型。

8、桩长为75m，桩径1.5 m。