桩土的接触问题，求教！急！

soarsun

我做的是个桩土共同作用的分析，桩与土接触的地方设置接触单元，但是接触单元设置后始终没有起作用，求指教！急！
我的命令流：
/prep7
Et,1,solid45
桩单元，土单元
注意接触单元号码
et,2,170
et,3,173
keyopt,3,9,0
keyopt,3,12,2
Mp,ex,1,2.0e9
桩的弹性模量
Mp,nuxy,1,0.25
桩的泊松比
Mp,dens,1,2000
桩的密度
Mp,ex,2,2.0e6
第一层土的弹性模量
Mp,nuxy,2,0.3
第一层土的泊松比
Mp,dens,2,1600
第一层土的密度!!
Tb,dp,2
Tbdata,1,19000，15，20
!!
mp,mu,3,0.2
r,1
r,2
r,3
建立断面模型，
k,1,0,0,0
k,2,5,0,0
k,3,5,0,-5
k,4,0,0,-5
k,5,0,-42,0
k,6,5,-42,0
k,7,5,-42,-5
k,8,0,-42,-5
v,1,2,3,4,5,6,7,8
wprot,0,-90,0
生成圆柱体
cylind,,1,0,-42,0,90
wprot,0,90,0
vsbv,1,2,sepo,delete
土体减去桩体
nummrg,all
合并了减去桩后的土体！！！！
numcmp,all
wprot,0,-90,0
再次生成圆柱体（桩）
cylind,,1,0,-21,0,90
wprot,0,90,0
wprot,0,-90,0
再次生成圆柱体（土）
wpoff,0,0,-21
cylind,,1,0,-21,0,90
wpoff,0,0,21
wprot,0,90,0
vsel,s,,,1,3,2
vplot
vglue,all
numcmp,all


allsel
vsel,s,,,3
vplot
vatt,2,2,1
lsel,s,,,1,10,9
lesize,all,0.5
lsel,s,,,4,7,3
lesize,all,0.5
lsel,s,,,2
lesize,all,0.5
vsweep,3

vsel,s,,,2
vplot
vatt,2,2,1
lsel,s,,,24,25,1
lesize,all,,,4
lsel,s,,,26
lesize,all,0.5
vsweep,2

vsel,s,,,1
vplot
vatt,1,1,1
lsel,s,,,18,19,1
lesize,all,,,4
lsel,s,,,22
lesize,all,0.5
vsweep,1
Allsel


Vsel,s,,,1                          刚性接触面单元
Vplot
Asel,s,,,5,7,2
Aplot
Nsla,s,1
Esln,s,0
Real,3
type,2
mat,3
Esurf,all
Esel,s,type,,2
Eplot
/psymb,esys,1

Allsel


vsel,s,,,2,3,1                        柔性接触面单元
vplot
Asel,s,,,10,15,5
Aplot
Nsla,s,1
Esln,s,0
Real,3
Type,3
Esurf,all
Esel,s,type,,3
Eplot
/psymb,esys,1
allsel
finish
/solu
加载求解
asel,s,loc,x,0
aplot
da,all,symm
asel,s,loc,x,5
aplot
da,all,all
asel,s,loc,z,0
aplot
da,all,symm
asel,s,loc,z,-5
aplot
da,all,all
asel,s,loc,y,-42
aplot
da,all,all
allsel
acel,0,9.8,0





asel,s,,,6
sfa,all,1,pres,500000



allsel
antype,static
设置分析选项并求解
nlgeom,on
time,1
nsub,100
outres,all,all
allsel
solv
finish

jiangsuwubing

可以试着用接触向导去创建接触对