做了一个桩土接触分析，大家讨论讨论

/filname,pile-soil contact 11_12_2
/prep7
et,1,plane42
et,2,solid45 !这里面用不上
keyopt,1,3,1
mp,ex,1,2.5e10 !桩的弹性模量
mp,nuxy,1,0.2 !桩的泊松比
mp,dens,1,2500 !桩的密度
mp,ex,2,2.5e8 !土的弹性模量
mp,nuxy,2,0.4 !土的泊松比
mp,dens,2,2000 !土的密度
tb,dp,2
tbdata,1,19,32,30 !粘聚力c为19,摩擦角为32度,膨胀角为30

RECTNG,0,1,0,8, !面1
RECTNG,1,10,0,7, !面2
RECTNG,0,10,0,-8, !面3
/pnum,area,1
/pnum,line,1
asel,s,,,2,3,1
aglue,all
numcmp,all
allsel
aplot

lsel,s,,,1,3,2 !划分面1
lesize,all,,,2
lsel,s,,,2,4,2
lesize,all,,,16
type,2
mat,1
amesh,1

lsel,s,,,6,8,2 !划分面2
lesize,all,0.5
lsel,s,,,5
lesize,all,0.5
lsel,s,,,7
lesize,all,0.5
type,2
mat,2
amesh,2

lsel,s,,,11 !划分面3
lesize,all,,,2
lsel,s,,,5
lesize,all,0.5
lsel,s,,,12
lesize,all,0.5
lsel,s,,,10
lesize,all,0.5
lsel,s,,,9
lesize,all,0.5
type,2
mat,2
amesh,3
allsel

!合并面3，面2重复单元
asel,s,,,2
asel,a,,,3
aplot
ESLV,S!选择当前有效体上的单元
numcmp,all
allsel
aplot

!接触单元设置
et,3,169
et,4,171
keyopt,4,9,0
keyopt,4,12,2
r,10
mp,mu,2,0.2

asel,s,,,1
nsla,s,1
nsel,r,loc,x,1
cm,target,line
r,10
type,3
esurf

allsel
asel,s,,,2
nsla,s,1
nsel,r,loc,x,1
cm,contact,line
r,10
type,4
esurf
allsel


!底端接触
allsel
asel,s,,,1
nsla,s,1
nsel,r,loc,y,0
cm,target,line
r,10
type,3
esurf

allsel
asel,s,,,3
nsla,s,1
nsel,r,loc,y,0
cm,contact,line
r,10
type,4
esurf
allsel,all

!设置边界条件
/solu
asel,s,,,2
nsla,s,1
dl,6,2,ux
asel,s,,,3
nsla,s,1
dl,12,3,ux
dl,9,3,uy
allsel

!施加重力加速度及荷载
!acel,0,9.8,0,
asel,s,,,1
!lplot
sfl,3,pres,600

!设置分析选项并求解
antype,static
nlgeom,on
time,1
auto,on
nropt,full
nsub,100,10000,100
outres,all,all
allsel
solve
finish

算了一下，结果云图不对，桩侧都没有反应，加了轴对称约束结果也一样，希望大家提供修改意见