各位大虾帮忙看下这个桩土接触的问题！

babecheer

各位大虾，小弟做的是一个桶形基础（桩）的桩土接触分析，在桶的一角施加了一个竖直向下的集中荷载，现在比较疑惑的是在基础附近发生了较大的位移，而基础附近的土体应力却没有什么变化，而应力则主要集中在基础上，觉得这个结果很不合理，因该是程序哪个地方出了问题，可是检查了好几遍都没找到问题所在，请各位高手指点一下,附命令流：
fini
/cle
/title,bucket   
/filname,bucket
*SET,bucket_ex,2e5                                           ！桶的参数
*SET,bucket_possion,0.125   
*SET,bucket_dens,7.86e-9
*SET,bucket_bithickness,2   
*SET,bucket_dithickness,10  
*SET,bucket_inradius,50
*SET,bucket_lowheight,100   
*SET,soil_ex,1.75                                     ！土体参数
*SET,soil_possion,0.4   
*SET,soil_dens,1.6e-9   
*SET,soil_c,0.019   
*SET,soil_fai,15
*SET,soil_pengzhang,30  
*SET,soil_radius,600
*SET,soil_height,1000   
/prep7  
et,1,plane42
et,2,solid45
mp,ex,1,bucket_ex   
mp,nuxy,1,bucket_possion
mp,dens,1,bucket_dens   
mp,ex,2,soil_ex
mp,nuxy,2,soil_possion  
mp,dens,2,soil_dens
tb,dp,2
tbdata,1,soil_c,soil_fai,soil_pengzhang
wprota,,-90
cyl4,,,bucket_inradius,,,180,soil_height
cyl4,,,bucket_inradius+bucket_bithickness,,bucket_inradius,180,soil_height  
cyl4,,,bucket_inradius+bucket_bithickness,,soil_radius,180,soil_height  
wpoffs,,,soil_height-bucket_dithickness-bucket_lowheight
vsel,s,,,1,3,2  
vsbw,all
vsel,s,,,2  
vsbw,all,sepo   
alls
numcmp,all  
wpoffs,,,bucket_lowheight   
vsel,s,,,4  
vsbw,all,sepo   
WPCSYS,-1,0
alls
numcmp,all  
vsel,s,,,2,6,4  
vglue,2,6   
vatt,1,,2   
vsel,s,,,1  
vsel,a,,,3,5,1  
vsel,a,,,7  
vglue,all   
vatt,2,,2   
alls
numcmp,all  
vsel,s,,,2  
lesize,13,,,20  
lesize,19,,,20  
lesize,37,,,4   
lesize,36,,,4   
lesize,21,,,12  
lesize,2,,,12   
lesize,31,,,4   
lesize,32,,,4   
amesh,14
vsweep,all  
vsel,s,,,3  
lesize,36,,,4   
lesize,37,,,4   
lesize,38,,,10  
lesize,42,,,10  
lesize,43,,,10  
lesize,47,,,4   
lesize,45,,,4   
lesize,46,,,12  
lesize,21,,,12  
amesh,19
vsweep,all  
vsel,s,,,4  
lesize,48,,,4   
lesize,49,,,4   
lesize,39,,,4   
lesize,41,,,4   
lesize,44,,,1   
lesize,50,,,1   
lesize,51,,,1   
lesize,26,,,12  
lesize,40,,,12  
amesh,24
vsweep,all  
vsel,s,,,1  
vsel,s,,,1  
lesize,19,,,20  
lesize,10,,,20  
lesize,22,,,1   
lesize,1,,,1
lesize,21,,,12  
lesize,12,,,12  
amesh,5
vsweep,all  
vsel,s,,,5  
lesize,34,,,1   
lesize,27,,,1   
lesize,51,,,1   
lesize,53,,,10  
lesize,28,,,11  
lesize,26,,,12  
lesize,24,,,12  
amesh,28
vsweep,all  
vsel,s,,,6  
vsel,s,,,6  
lesize,12,,,12  
lesize,14,,,20  
lesize,10,,,20  
lesize,56,,,12,0.8  
lesize,54,,,12,0.8  
lesize,16,,,12  
lesize,9,,,12   
amesh,30
vsweep,all  
vsel,s,,,7  
vsel,s,,,7  
lesize,17,,,11  
lesize,59,,,11  
lesize,8,,,12,0.8   
lesize,56,,,12,0.8  
lesize,7,,,12   
lesize,5,,,12   
lesize,12,,,12  
amesh,35
vsweep,all  
save
alls
FLST,2,35,5,ORDE,2  
FITEM,2,1   
FITEM,2,-35
ACLEAR,P51X
vsel,s,,,4,5,1  
eslv,s  
nslv,s,1
eplot   
nummrg,all  
numcmp,all  
vsel,s,,,1,3,1  
vsel,a,,,6,7,1  
eslv,s  
nslv,s,1
eplot   
nummrg,all  
numcmp,all  
alls
et,3,170                              ！设置接触面
et,4,173
keyopt,4,9,0
keyopt,4,12,2   
r,10
real,10
mp,mu,2,0.2
vsel,s,,,4,5,1  
vplot   
asel,s,,,11,15,4
asel,a,,,17,26,9
aplot   
cm,bucket_target,area   
type,3  
nsla,s,1
esurf,all   
/VIEW,1,1,1,1   
/ANG,1  
/REP,FAST   
esel,s,type,,3  
eplot   
alls
asel,s,,,7,18,11
asel,a,,,21,34,13   
aplot   
cm,soil_contact,area
type,4  
nsla,s,1
esurf,all   
esel,s,type,,4  
eplot   
alls
gplot   
allsel  
nsel,s,,,1148   
f,all,fy,-30                             ！在一个节点上施加集中力
FINISH  
/SOL
WPCSYS,-1,0
csys,0  
asel,s,loc,z,0  
da,all,symm
asel,s,loc,y,0  
da,all,all  
csys,5  
asel,s,loc,x,soil_radius
da,all,all  
csys,0  
antype,static   
allsel  
acel,0,9800,0   
nlgeom,on   
time,1  
nsubst,10   
outres,all,all  
alls
solve

yuansu

你后来找到原因了吗？我也有类似的问题啊