ansys 接触问题不收敛 求助!!

finish
/clear
/filname,dangenmaoganshiyan1
/title,paper
*afun,deg
/prep7
r0=100
r1=20
l0=6000

et,1,plane42
et,2,solid45
et,3,conta173
keyopt,3,1,0
keyopt,3,12,2
et,4,targe170
!钢筋
mptemp,1,0
mp,ex,1,2.0e11
mp,nuxy,1,0.3
!mp,dens,1,9870
mp,mu,1,0.5
!tb,biso,1 !定义了双线性等向强化材料
!tbdata,1,930e6,
!/xrange,0,0.01
!tbplot,biso,0
!砂浆
mp,ex,2,2.63e10
mp,nuxy,2,0.19
!mp,dens,2,2100
mp,mu,2,0.75
!tb,dp,2,,,
!tbdata,1,530000,37,34
!基岩
mp,ex,3,3.3e10
mp,nuxy,3,0.2
!mp,dens,3,2780
mp,mu,3,0.5
!tb,dp,3,,,
!tbdata,1,900000,39,35

!建模
k,1000,0,0,0
circle,1000,100
circle,1000,20
l,1,5
l,2,6
l,3,7
l,4,8
a,1,2,6,5
a,2,3,7,6
a,3,4,8,7
a,4,1,5,8
lsel,s,line,,1,4,1
lsel,a,line,,5,8,1
lsel,a,line,,9,12,1
lesize,all,,,3
amesh,all
extopt,esize,60,0
type,2
mat,2
extopt,aclear,1
vext,1,4,1,,,6000 !1michang



circle,1000,20
a,17,18,19,20
lsel,s,line,,33,36,1
lesize,all,,,3 !钢筋划分网格
type,1
amesh,21
extopt,esize,60,0
type,2
mat,1
extopt,aclear,1
vext,21,,,,,6000 !6michang

vsel,s,volu,,5,,
vplot
allsel,all
vplot
csys,0 !注意这的坐标系

circle,1000,800
circle,1000,100
l,25,29
l,26,30
l,27,31
l,28,32
a,25,26,30,29
a,26,27,31,30
a,27,28,32,31
a,28,25,29,32
lsel,s,line,,45,48,1
lsel,a,line,,49,52,1
lesize,all,,,3
lsel,s,line,,53,56,1
lesize,all,,,10,1.2
allsel
aplot
asel,s,area,,27,30,1
aplot
type,1
amesh,all

extopt,esize,60,0
type,2
mat,3
extopt,aclear,1
vext,27,30,1,,,6000 !1michang
vsel,s,volu,,6,9,1
vplot
allsel,all
vplot
!接触对的生成
r,1
RMODIF,1,13,530*1000 !定义初始凝聚力
RMODIF,1,9,650*1000 !定义等效剪应力
real,1
vsel,s,mat,,2
vplot
asel,s,area,,8,,
aplot
cm,target,area
type,4
mat,2
nsla,s,1
esurf,all

allsel,all
vsel,s,mat,,1
vplot
asel,s,area,,23,,
aplot
cm,contact,area
type,3
mat,1
nsla,s,1
esurf,all

r,2
RMODIF,2,13,530*1000 !定义初始凝聚力
RMODIF,2,9,650*1000 !定义等效剪应力
real,2
vsel,s,mat,,2
vplot
asel,s,area,,13,,
aplot
cm,target,area
type,4
mat,2
nsla,s,1
esurf,all

allsel,all
vsel,s,mat,,1
vplot
asel,s,area,,24,,
aplot
cm,contact,area
type,3
mat,1
nsla,s,1
esurf,all

r,3
RMODIF,3,13,530*1000 !定义初始凝聚力
RMODIF,3,9,650*1000 !定义等效剪应力
real,3
vsel,s,mat,,2
vplot
asel,s,area,,17,,

aplot
cm,target,area
type,4
mat,2
nsla,s,1
esurf,all

allsel,all
vsel,s,mat,,1
vplot
asel,s,area,,25,,
aplot
cm,contact,area
type,3
mat,1
nsla,s,1
esurf,all

r,4
RMODIF,4,13,530*1000 !定义初始凝聚力
RMODIF,4,9,650*1000 !定义等效剪应力
real,4
vsel,s,mat,,2
vplot
asel,s,area,,20,,
aplot
cm,target,area
type,4
mat,2
nsla,s,1
esurf,all
allsel,all
vsel,s,mat,,1
vplot
asel,s,area,,26,,
aplot
cm,contact,area
type,3
mat,1
nsla,s,1
esurf,all

!外接触
r,5
RMODIF,5,13,900*1000 !定义初始凝聚力
RMODIF,5,9,1000*1000 !定义等效剪应力
real,5
vsel,s,mat,,3
vplot
asel,s,area,,34,,
aplot
cm,target,area
type,4
mat,3
nsla,s,1
esurf,all
allsel,all
vsel,s,mat,,2
vplot
asel,s,area,,6,,
aplot
cm,contact,area
type,3
mat,2
nsla,s,1
esurf,all


r,6
RMODIF,6,13,900*1000 !定义初始凝聚力
RMODIF,6,9,1000*1000 !定义等效剪应力
real,6
vsel,s,mat,,3
vplot
asel,s,area,,39,,
aplot
cm,target,area
type,4
mat,3
nsla,s,1
esurf,all
allsel,all
vsel,s,mat,,2
vplot
asel,s,area,,11,,
aplot
cm,contact,area
type,3
mat,2
nsla,s,1
esurf,all


r,7
RMODIF,7,13,900*1000 !定义初始凝聚力
RMODIF,7,9,1000*1000 !定义等效剪应力
real,7
vsel,s,mat,,3
vplot
asel,s,area,,43,,
aplot
cm,target,area
type,4
mat,3
nsla,s,1
esurf,all
allsel,all
vsel,s,mat,,2
vplot
asel,s,area,,15,,
aplot
cm,contact,area
type,3
mat,2
nsla,s,1
esurf,all

r,8
RMODIF,8,13,900*1000 !定义初始凝聚力
RMODIF,8,9,1000*1000 !定义等效剪应力
real,8
vsel,s,mat,,3
vplot
asel,s,area,,46,,
aplot
cm,target,area
type,4
mat,3
nsla,s,1
esurf,all
allsel,all
vsel,s,mat,,2
vplot
asel,s,area,,19,,
aplot
cm,contact,area
type,3
mat,2
nsla,s,1
esurf,all
allsel,all
eplot

!求解设置
/solu
!设置约束
asel,s,area,,27,30,1
aplot
da,all,all
allsel
asel,s,area,,32,,,
asel,a,area,,37,,,
asel,a,area,,41,,,
asel,a,area,,45,,,
aplot
da,all,symm
da,all,all
alls
dtran
!施加荷载
allsel,all
!acel,0,0,9.8
asel,s,area,,22,,
aplot
f0=500*1000 !加载方式500
a0=3.14*r1**2
sfa,all,,pres,-f0/a0
allsel,all
ftran
lswrite
asel,s,area,,22,,
aplot
f0=500*1000 !加载1000
a0=3.14*r1**2
sfa,all,,pres,-f0/a0
allsel,all
ftran
lswrite
asel,s,area,,22,,
aplot
f0=500*1000 !加载1500
a0=3.14*r1**2
sfa,all,,pres,-f0/a0
allsel,all
ftran
lswrite
!设置求解选项
allsel,all
antype,static
nlgeom,on
!EQSLV, PCGOUT,1.0e-5,0, !选用预条件共轭梯度求解器PCG求解
time,1
nsub,100,500,10
autots,on
neqit,50
cnvtol,f,,0.05,2,2, !cnvtol,f,100,0.05,2,,
nropt,full,,on
outress,all,10
pred,on
lnsrch,on
alls
lssolve,1,3
save


不收敛啊，请问是哪一步出错了呢？