[交流] 【求助】桩土接触分析

大家帮我看一下这个命令流，!ex8.31 桩土接触分析
finish
/clear
/prep7
apile=1.0 !桩平面尺寸
hpil1=2.0 !桩地面以上高度
hpil2=8.0 !桩入土深度
asoil=7.0 !土体平面尺寸
hsoil=24.0 !土体深度

et,1,solid45
et,2,170
et,3,173
keyopt,3,9,0
keyopt,3,12,2

mp,ex,1,2.6e10 !桩的材料性质
mp,prxy,1,0.167
mp,dens,1,2500

mp,ex,2,2.6e8 !土材料性质
mp,prxy,2,0.42
mp,dens,2,1900
tb,dp,2
tbdata,,19000,31,29

mp,mu,3,0.2 !摩擦系数

r,1
r,2
r,3
!1/4模型
cylind,0,1,0,-hpil2,270,360
cylind,1,7,0,-8,270,360
wpoff,,,-8
cylind,0,1,0,-16,270,360
cylind,1,7,0,-16,270,360

vsel,s,,,2,4,1
vplot
vglue,2,3,4

vsel,s,,,5
aslv,s
lsla,s
lplot
lesize,22,0.05
lesize,23,0.05
lesize,24,,,20
lesize,30,0.05
lesize,45,0.05
lesize,46,0.05
lesize,11,,,20
lesize,43,0.05
lesize,44,0.05
type,1
mat,1
mshape,0,3d
mshkey,1
vmesh,5

vsel,s,,,6
aslv,s
lsla,s
lplot
lesize,50,0.05
lesize,12,0.05
lesize,10,0.05
lesize,48,0.05
lesize,47,0.05
lesize,49,0.05
lesize,13,,,20
lesize,34,,,20
type,1
mat,1
mshape,0,3d
mshkey,1
vmesh,6

vsel,s,,,1
aslv,s
lsla,s
lplot
lesize,1,0.05
lesize,2,0.05
lesize,5,0.05
lesize,6,0.05
lesize,7,0.05
lesize,8,0.05
lesize,9,0.05
lesize,3,,,20
lesize,4,,,20
type,1
mat,2
mshape,0,3d
mshkey,1
vmesh,1

vsel,s,,,2
aslv,s
lsla,s
lplot
lesize,20,0.05
lesize,21,0.05
lesize,18,0.05
lesize,19,0.05
lesize,15,0.05
lesize,17,0.05
lesize,16,,,20
lesize,14,,,20
type,1
mat,1
mshape,0,3d
mshkey,1
vmesh,2

这里到了网格划分，可是当弄上接触单元之后，就出现了附件中的错误，高手帮我看一下是网格划分的太密吗？谢谢

回复此楼