请高手帮忙看一下，为什么结果会不收敛呢
/CLEAR,NOSTART
/filname,pile-soil contact
/title，pile
/units,si !国际单位制
fe=0.3 !摩擦系数
fk1=10 !法向接触刚度
fk2=1 !法向接触刚度
fd=3 !切向接触刚度
ftoln=1 !初始渗透因子
!进入前处理器
/prep7
et,1,solid45
et,2,170 !目标单元
et,3,173 !三维四节点四边形接触单元，174三维8节点四边形接触单元
keyopt,3,4,1 !4，定义接触面
keyopt,3,5,1
keyopt,3,9,0 !消除初始渗透
keyopt,3,12,2 !不分开接触，允许滑动，法向不分开
!定义材料属性************************************************************
mp,mu,1,0.2 !桩体材料属性
mp,ex,1,2.5e10
mp,nuxy,1,0.2
mp,dens,1,2500
mp,ex,2,2.0e8 !土层2材料属性
mp,nuxy,2,0.4
mp,dens,2,2000
mp,ex,3,1.8e9 !岩层3材料属性
mp,nuxy,3,0.29
mp,dens,3,2600
!定义参数*****************************************************************
tb,concr,1 !桩体的参数
tbdata,,0.7,0.9,2.6,-1
tb,dp,2 !土层2的D-P参数
tbdata,1,19,32,30
tb,dp,3 !岩层3的D-P参数
tbdata,1,0,35,0
!建立模型*****************************************************************
cylind,0,0.1,0,7,0,90
block,0,0.350,0,0.300,0,1
block,0,0.350,0,0.300,1,7
!布尔操作*****************************************************************
VSEL,S,VOLU,,1,3
VSBV,2,1,SEPO,,KEEP
VSBV,3,1,SEPO,,KEEP
numcmp,all
!网格划分*****************************************************************
allsel,all
vsel,s,,,1 !桩体附属性，划分网格
aslv,s
lsla,s
lplot
lesize,1,,,4
lesize,2,,,4
lesize,3,,,4
lesize,4,,,4
lesize,5,,,4
lesize,6,,,4
lesize,7,,,14
lesize,8,,,14
lesize,9,,,14
type,1
mat,1
mshape,0,3d
mshkey,1
vmesh,all
vsel,s,,,2 !土层2附属性，划分网格
vplot
aslv,s
lsla,s
lplot
lesize,18,,,2
lesize,38,,,2
lesize,17,,,2
lesize,39,,,2
lesize,13,,,4
lesize,10,,,4
lesize,14,,,4,4
lesize,21,,,4,4
lesize,37,,,4,4
lesize,15,,,4,4
lesize,20,,,12
lesize,19,,,12
lesize,16,,,12
lesize,11,,,12
lesize,12,,,12
ACCAT,8,10
LCCAT,18,38
LCCAT,17,39
type,1
mat,2
mshape,0,3d
mshkey,1
vmesh,all
LDELE,40,41
ADELE,20
vsel,s,,,3 !岩层3附属性，划分网格
vplot
aslv,s
lsla,s
lplot
lesize,29,,,2
lesize,35,,,2
lesize,36,,,2
lesize,30,,,2
lesize,22,,,4
lesize,25,,,4
lesize,33,,,4,4
lesize,34,,,4,4
lesize,27,,,4,4
lesize,26,,,4,4
lesize,24,,,2
lesize,23,,,2
lesize,28,,,2
lesize,31,,,2
lesize,32,,,2
ACCAT,16,14
LCCAT,35,29
LCCAT,36,30
type,1
mat,3
mshape,0,3d
mshkey,1
vmesh,all
LDELE,40,41
ADELE,20
!用接触向导建立接触对*************************************************************
ALLSEL,ALL
VPLOT
!桩侧，建立目标面
r,1,,,fk1,ftoln, , , !定义单元实常数
rmore,,,,,,5
r,2,,,fk2,ftoln, , ,
rmore,,,,,,5
mp,mu,6,0.2
allsel,all
vsel,s,,,1
asel,s,,,3
aplot
cm,target1,area !由选择实体生成一个原件
type,2
real,1
mat,6
nsla,s,1
esurf,all
!建立接触面
vsel,s,,,2,3
aslv,s
asel,s,,,6
asel,a,,,12
/pnum,area,1
aplot
cm,contact1,area
type,3
real,1
mat,6
nsla,s,1 !选择依附面上的节点 1表示所有节点 0表示仅面上节点
esurf,all !在选择自由面上生成重叠单元
!滑动面，建立目标面*************************
allsel
vsel,s,,,3
asel,s,,,18
aplot
cm,target2,area !由选择实体生成一个原件
type,2
real,2
mat,6
nsla,s,1 !选择依附面上的节点 1表示所有节点 0表示仅面上节点
esurf,all !在选择自由面上生成重叠单元
! 滑动面，建立接触面 ***********************
allsel,all
vsel,s,,,2
aslv,s
asel,s,,,11
/pnum,area,1
aplot
cm,contact2,area
type,3
real,2
mat,6
nsla,s,1
esurf,all
ALLSEL,ALL
!镜像****************************************************************************
VSYMM,X,1,3
VSYMM,Y,1,6
NUMMRG,node
NUMMRG,elem
NUMMRG,all
NUMCMP,ALL
!施加约束
/solu
allsel,all
da,28,all
da,56,all
da,33,all
da,60,all
da,48,all
da,16,all
da,43,all
da,10,all
!设置非线性选项
acel,0,0,-9.8 !定义重力加速度
neqit,50, !迭代次数
nropt,Modi ! 用改变的牛顿-拉普拉斯方程求解
nlgeom,on ! 打开大变形选项
autots,on ! 打开自动时间步
lnsrch,on !自动线性搜索
pred,on !打开预测求解器
outpr,all,all
OUTRES,all,all
allsel,all
cncheck,detail !检查接触状态
cncheck,summary
/solu
CNCHECK,POST
FINISH
!施加荷载
/solu
allsel,all
sfa,8,,pres,100
sfa,27,,pres,100
solve
