Ansys 桩土水平作用分析自由度超限

我刚刚入手ANSYS做个小项目，关于试验的。该项目是关于桩土共同作用的。不考虑桩的竖向外力（施加自重即可），只考虑水平外力。给桩施加水平荷载，看桩身变形。土体根据勘测数据建立了分成模型。但是我在运行后，总是出现DOF EXCEEDED，MAX....,LIMIT 0.0000。然后运算没有结果，error很多，运算时间也比较长。有没有高手帮我看看。多谢多谢！
finish
/clear
/prep7
*afun,deg
et,1,solid45
et,2,solid65
et,3,plane42
keyopt,2,7,1
!定义土体模型材料参数
mp,ex,1,1.8e+6
mp,nuxy,1,0.4999
mp,dens,1,1640
tb,dp,1
tbdata,1,7.7e3,2.2,0
mp,ex,2,6.61e+6
mp,nuxy,2,0.4999
mp,dens,2,1980
tb,dp,2
tbdata,2,16.6e3,10.4,0
mp,ex,3,5.3e+6
mp,nuxy,3,0.4999
mp,dens,3,1960
tb,dp,3
tbdata,3,36.5e3,9.4,0
mp,ex,4,6.4e+6
mp,nuxy,4,0.4999
mp,dens,4,1910
tb,dp,4
tbdata,4,61.7e3,13.5,0
!定义桩模型材料参数
r,2,5,0.0128,0,90
MP,EX,5,2.0E11
mp,nuxy,5,0.3
mp,ex,6,3e10
mp,nuxy,6,0.2
mp,dens,6,2500
tb,dp,6
tbdata,6,0.35,0.9,1.43e6,14.3e6,,,,,,
!建立几何模型
block,0,18,0,18,6,14
block,0,18,0,18,1,6
block,0,18,0,18,0,1
/view,1,1,1,1
cyl4,6,9,0.2,,0.2,,15 !4
cyl4,12,9,0.2,,0.2,,15 !5
vsbv,1,4,sepo,,keep
vsbv,2,4,sepo,,keep
vsbv,3,4,sepo,,keep
vsbv,6,5,sepo,,keep !3 上
vsbv,1,5,sepo,,keep !6 中
vsbv,2,5,sepo,,keep !1 下
!网格划分
!分割体
wpstyl,,,,,,,,1
wprota,,,-90
wpoffs,,,-6
vsbw,3
vsbw,6
vsbw,1
wpoffs,,,-3
vsbw,7
vsbw,8
vsbw,9
wpoffs,,,-3
vsbw,1
vsbw,7
vsbw,8
wpstyl,,,,,,,,0

esize,1
vglue,2,10,9,13,3,11,14,1
vglue,18,20,17,1,6,12,15,7
!对土体分配属性
vsel,s,volu,,2,8,6
vsel,a,volu,,16,19,3
vatt,1,,1
vsel,s,volu,,9,11,2
vsel,a,volu,,1,3,2
vatt,2,,1
vsel,s,volu,,6,10,4
vsel,a,volu,,13,14,1
vatt,3,,1
!对面映射划分网格
amap,38,14,12,10,16
amap,56,16,10,12,14
amap,85,104,108,102,106
amap,131,106,102,108,104
!土体扫掠
vsel,s,volu,,2,9,7
vsel,a,volu,,6
vsweep,all,38,61
vsel,s,volu,,8,11,3
vsel,a,volu,,10
vsweep,all,56,80
vsel,s,volu,,3,16,13
vsel,a,volu,,13
vsweep,all,85,106
vsel,s,volu,,1,19,18
vsel,a,volu,,14
vsweep,all,131,100
!对桩1进行切割
wpstyl,,,,,,,,1
vsel,s,,,4,5,1
vsbw,5
wpoffs,0,0,6
vsbw,4
kwpave,29
wprota,,-90
vsel,s,,,5,7,2
vsel,a,,,12,15,3
vsbw,all
wpoffs,0,-0.1
wprota,,135
vsel,s,,,4,17,13
vsel,a,,,23,24,1
vsbw,17
wprota,0,-90
vsbw,4
wprota,,45
wpoffs,,,0.2
wprota,,45
vsbw,23
wprota,,-90
vsbw,24
wprota,,45
vadd,23,5,12,4 !合并桩中间体
aadd,118,23,20,90
aadd,4,117,21,19
vglue,7,17,15,24,25
!对桩2进行切割
kwpave,37
vsel,s,,,20,22,1
vsel,a,,,18
wpoffs,,,0.1
wprota,,45
vsbw,20
wprota,,90
vsbw,18
wprota,,-135
wpoffs,,,0.2
wprota,,-45
vsbw,22
wprota,,90
vsbw,21
vadd,18,4,22,12 !合并桩中间体
aadd,99,19,115,103
aadd,4,98,102,114
vglue,23,26,20,21,5
wpstyl,,,,,,,,0
!对桩1进行网格划分
vsel,s,,,23,26,1
vsel,a,,,15,17,2
vsel,a,,,5,7,2
vsel,a,,,20,21,1
vatt,6,2,2 !对桩分配属性
esize,0.025
amap,53,46,29,30,48
amap,135,48,30,31,52
amap,132,52,31,32,50
amap,122,50,32,29,46
amesh,22
esize,0.4
vsweep,25,135,134
vsweep,7,53,52
vsweep,17,132,125
vsweep,15,122,97
vsweep,24,22,20
!对桩2进行网格划分
esize,0.025
amap,138,58,38,39,69
amap,141,69,39,40,60
amap,121,60,40,37,57
amap,95,57,37,38,58
amesh,21
esize,0.4
vsweep,23,138,137
vsweep,5,95,90
vsweep,26,141,140
vsweep,20,121,120
vsweep,21,21,19

!接触建模
tshap,cyli
et,4,170 !target170目标面
et,5,173 !contac173接触面
keyopt,5,9,0
keyopt,5,12,1 !考虑只受水平力，且变形量不大，因此不考虑接触面滑动
r,4,0.2
type,4
asel,s,,,91,94,3
asel,a,,,25,26,1
asel,a,,,65,66,1
asel,a,,,34,47,13
nsla,s,1
esurf
type,5
asel,s,,,9,10,1
asel,a,,,35,36,1
asel,a,,,77,129,52
asel,a,,,130,180,50
asel,a,,,50,70,20
asel,a,,,51,71,20
asel,a,,,29,39,10
asel,a,,,6,13,7
asel,a,,,31,32,1
asel,a,,,58,63,5
asel,a,,,75,76,1
asel,a,,,88,89,1
nsla,s,1
esurf
/solu
asel,s,loc,y,18
asel,a,loc,y,0
asel,a,loc,x,18
asel,a,loc,x,0
da,all,ux,0
da,all,uy,0
da,all,uz,0
da,all,rotx,0
da,all,roty,0
da,all,rotz,0
asel,s,,,100,106,6
asel,a,,,61,80,19
da,all,uy,0
da,all,rotx,0
da,all,roty,0
da,all,rotz,0
asel,s,,,20,52,32
asel,a,,,125,134,9
asel,a,,,90,97,7
asel,a,,,137,140,3
asel,a,,,19,120,101
da,all,uz,0
da,all,rotx,0
da,all,roty,0
da,all,rotz,0
allsel
acel,,,9.8 !设置重力加速度
fk,29,fx,-1000
fk,39,fx,1000 !施加荷载
!设置分析选项
antype,static
nlgeom,on !大变形开关打开
nsubst,50
outres,all,all
solve
finish