小女子大菜鸟一个，初学ansys,现学现卖，想做个接触分析，是混凝土和石砌体之间的剪应力，命令流如下，算来算去都不收敛，着急的很啊，请高手帮我看下是怎么回事？有好心人也可以加我qq8376330指点一下（可以有偿），谢谢好心人！
finish
/clear
/filname,modeltest19
/prep7
et,1,65
et,2,65
mp,ex,1,11143.5
mp,prxy,1,0.2
mp,mu,1,0.7
*SET,fc1,11.73
*SET,ft1,1.04
tb,concr,1 !定义混凝土
tbdata,,0.5,0.95,ft
tb,miso,1,,11
tbpt,,0.0002,fc1*0.19
tbpt,,0.0004,fc1*0.36
tbpt,,0.0006,fc1*0.51
tbpt,,0.0008,fc1*0.64
tbpt,,0.001,fc1*0.75
tbpt,,0.0012,fc1*0.84
tbpt,,0.0014,fc1*0.91
tbpt,,0.0016,fc1*0.96
tbpt,,0.0018,fc1*0.99
tbpt,,0.002,fc1
tbpt,,0.0033,fc1*0.85
mp,ex,2,2600 ！定义砌体
mp,prxy,2,0.25
mp,mu,2,0.7
*SET,fc2,15.89
*SET,ft2,1.09
tb,concr,2
tbdata,,0.5,0.95,ft2
tb,miso,2,,5
tbpt,,0.001,2.6
tbpt,,0.0015,3.68
tbpt,,0.0045,5.78
tbpt,,0.0072,0.47
tbpt,,0.02,4.17
tbpt,,0.0072,6.47
k,1 ！建立模型并划分网格
k,2,350
k,3,350,100
k,4,200,125
k,5,200,350
k,6,0,350
a,1,2,3,4,5,6,
k,7,200,150
k,8,350,150
k,9,350,500
k,10,50,500
k,11,50,400
k,12,200,375
a,7,8,9,10,11,12
voffst,1,200
voffst,2,200
mshkey,1
esize,25
vsweep,all
!定义接触单元
et,3,170
et,4,174
keyopt,4,9,1
keyopt,4,12,2
r,2,,,0.01,0.4,0.1
rmore,,,,,,,
rmore,300000
asel,s,,,7
nsla,s,7
nsel,r,loc,y,150,350
real,2
type,3
esurf
allsel,all
asel,s,,,16
nsla,s,16
nsel,r,loc,y,150,350
real,2
type,4
esurf
allsel,all
finish
/solu
allsel,all ！施加约束
asel,s,loc,y,0
nsla,s,1
d,all,all
nsel,s,loc,y,500
d,all,uy,-0.04 ！位移加载并求解
allsel,all
nlgeom,on
time,1
autots,on
nsubst,100,1000,100
outres,all,all
solve