请教高手桩土接触的问题!!!

我的毕业论文要用到ANSYS做桩和土的接触分析,可是我是一个学俄语的,所以很吃力,下面是我自己写的桩突接触的程序,但是总是的不到结果,不知道是什么原因,希望高手们能帮我分析一下,小弟在这里先谢谢了!!!!!!!!



/filname,zhuangtu
/title，xuebao
!/config,nres,2000000
/units,si !国际单位制
fe=0.2 !摩擦系数
fk=5 !法向接触刚度
fd=3 !切向接触刚度
ftoln=0.1 !初始渗透因子

/PREP7
et,1,plane42
ET,2,SOLID45
ET,3,TARGE170
ET,4,CONTA173

KEYOPT,4,12,2 !!!设置摩擦类型，法向不分开，但是允许滑动

MP,EX,1,2.5E+010 !桩的参数
MP,PRXY,1,0.2
MP,DENS,1,2500

MP,EX,2,3.8E+006 ！第一层土的参数
MP,PRXY,2,0.3
TB,DP,2,,,
TBMODIF,1,1,15000
TBMODIF,1,2,21.5
TBMODIF,1,3,
MP,DENS,2,1857

MP,EX,3,2.5E+006 ！第二层土的参数
MP,PRXY,3,0.4
TB,DP,3,,,
TBMODIF,1,1,16800
TBMODIF,1,2,19.2
TBMODIF,1,3,
MP,DENS,3,1690

MP,DENS,4,1890
MP,EX,4,1.0E+010
MP,PRXY,4,0.23
TB,DP,4,,,
TBMODIF,1,1,0
TBMODIF,1,2,30
TBMODIF,1,3,

MP,EX,5,8.8E+006
MP,PRXY,5,0.3
TB,DP,5,,,
TBMODIF,1,1,25000
TBMODIF,1,2,23
TBMODIF,1,3,
MP,DENS,5,1800

MP,EX,6,1.2E+007 ！第五层土的参数
MP,PRXY,6,0.21
TB,DP,6,,,
TBMODIF,1,1,0
TBMODIF,1,2,35
TBMODIF,1,3,
MP,DENS,6,1900

!!!!!!!!!!!!!!!!!!!!!!建模

block,0,2,0,-0.4,0,-2
block,0,2,-0.4,-1.6,0,-2
block,0,2,-1.6,-5.6,0,-2
block,0,2,-5.6,-10.6,0,-2
block,0,2,-10.6,-15,0,-2
block,0,2,-15,-30,0,-2
wpro,,90.000000,
CYLIND, ,0.35,0,30,270,360,
allsel
FLST,2,6,6,ORDE,2
FITEM,2,1
FITEM,2,-6
VSBV,P51X,7
! vplot
CYLIND, 0.2,0.35,0.4,1.6,270,360,
CYLIND, 0.2,0.35,1.6,5.6,270,360,
CYLIND, 0.2,0.35,5.6,10.6,270,360,
CYLIND, 0.2,0.35,10.6,15,270,360,
CYLIND, 0.2,0.35,15,30,270,360,
allsel
vglue,all
nummrg,all
numcmp,all
CYLIND, ,0.35,0,0.4,270,360,
CYLIND, ,0.2,0.4,1.6,270,360,
CYLIND, ,0.2,1.6,5.6,270,360,
CYLIND, ,0.2,5.6,10.6,270,360,
CYLIND, ,0.2,10.6,15,270,360,
CYLIND, ,0.2,15,30,270,360,
wpro,,-90.000000,

WPSTYLE,,,,,,,,0

vadd,12,13,14,15,16
vglue,17,5
NUMCMP,LINE
NUMCMP,VOLU
！NUMCMP,KP

!!!!!!!!!!!!!!!!!!!!!!!!!土体网格的划分

lsel,s,,,40,67,27 ！！！！vsel,s,,, aslv,all, lsla,all lplot
lsel,a,,,39,65,26
lsel,s,,,66,68,2
lsel,a,,,1
lsel,a,,,3,6,3
lsel,s,,,9,38,29
allsel
accat,1,6
lccat,2,11
lccat,10,4
lsel,s,,,10,11,1
lsel,a,,,2,4,2
lesize,all,,,5
lsel,s,,,40,67,27
lsel,a,,,39,65,26
lesize,all,,,10,0.3
lsel,s,,,66,68,2
lsel,a,,,1
lsel,a,,,3,6,3
lesize,all,,,4
lsel,s,,,9,38,29
lesize,all,,,10
type,2
mat,2
vmesh,6

allsel
accat,34,35
lccat,10,4
lccat,13,18
lsel,s,,,10,18,8
lsel,a,,,4,13,9
! lplot
lesize,all,,,5
lsel,s,,,67,72,5
lsel,a,,,65,73,8
! lplot
lesize,all,,,10,0.3
lsel,s,,,17,19,2
lsel,a,,,69
lsel,a,,,70,71,1
lesize,all,,,12
lsel,s,,,9,12,3
lesize,all,,,10
type,2
mat,2
vmesh,7

allsel
accat,40,39
lccat,20,23
lsel,s,,,20,23,3
lesize,all,,,5
lsel,s,,,77,78,1
lesize,all,,,10,0.3
lsel,s,,,49,50,1
lsel,a,,,74,76,1
lesize,all,,,40
lsel,s,,,25
lesize,all,,,10
type,2
mat,3
vmesh,8

allsel
accat,45,44
lccat,28,24
lsel,s,,,24,28,4
lesize,all,,,5
lsel,s,,,82,83,1
lesize,all,,,10,0.3
lsel,s,,,53,54,1
lsel,a,,,80,81,1
lsel,a,,,79
lesize,all,,,50
lsel,s,,,33
lesize,all,,,10
type,2
mat,4
vmesh,9

allsel
accat,49,50
lccat,30,31
lsel,s,,,30,31,1
lesize,all,,,5
lsel,a,,,87,88,1
lesize,all,,,10,0.3
lsel,s,,,57,58,1
lsel,a,,,84,86,1,
lesize,all,,,44
lsel,s,,,41
lesize,all,,,10
type,2
mat,5
vmesh,10

allsel
accat,54,55
lccat,32,34
lsel,s,,,32,34,2
lesize,all,,,5
lsel,a,,,92,93,1
lesize,all,,,10,0.3
lsel,s,,,89,91,1
lsel,a,,,61,62,1,
lesize,all,,,150
lsel,s,,,45
lesize,all,,,10
type,2
mat,6
vmesh,11

allsel
lsel,s,,,7,9,2
lsel,a,,,12,15,3
lesize,all,,,10
lsel,s,,,21,22,1
lsel,a,,,17,19,2
lesize,all,,,12
lsel,s,,,14,16,2
lsel,a,,,5,8,3
lesize,all,,,3
type,2
mat,2
vmesh,1

allsel
lsel,s,,,25,27,2
lesize,all,,,10
lsel,s,,,49,52,1
lesize,all,,,40
lsel,s,,,26,29,3
lesize,all,,,3
type,2
mat,3
vmesh,2
allsel
lsel,s,,,33,36,3
lesize,all,,,10
lsel,s,,,55,56,1
lesize,all,,,50
lsel,s,,,35,37,2
lesize,all,,,3
type,2
mat,4
vmesh,3
allsel
lsel,s,,,41,43,2
lesize,all,,,10
lsel,s,,,59,60,1
lesize,all,,,44
lsel,s,,,42,44,2
lesize,all,,,3
type,2
mat,5
vmesh,4

allsel
lsel,s,,,45,47,2
lesize,all,,,10
lsel,s,,,46,48,2
lesize,all,,,3
type,2
mat,6
vmesh,5

TYPE,2
MAT, 6
ESIZE,0.08
VSEL,s,,,12
VMESH,12

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!以下是桩的网格划分

TYPE,2
MAT,1
ESIZE,0.05, ,
MSHAPE,1,3D
MSHKEY,0 !!!! 用的是自由划分
VSEL,s,,,13
vplot
VMESH,13

!!!!!!!!!!下面是用鼠标操作的去除网格划分时连接的面和线

/PREP7
CM,_Y,AREA
ASEL,R,ACCA
*SET,_Z1,ARINQR(0,13)
*IF,_z1,ne,0,then
ADELE,ALL
*SET,_Z2,ARINQR(0,13)
*SET,_Z3,_Z1-_Z2
*IF,_Z3,NE,0,THEN
CMSEL,S,_Y
CMDELE,_Y
*ENDIF
*ELSE
CMSEL,S,_Y
CMDELE,_Y
*ENDIF
!*
CM,_Y,LINE
*SET,_Z1,LSINQR(0,13)
*IF,_z1,ne,0,then
LSEL,R,LCCA
*SET,_Z2,LSINQR(0,13)
LDELE,ALL
*SET,_Z3,_Z1-_Z2
*IF,_Z3,NE,0,THEN
CMSEL,S,_Y
CMDELE,_Y
*ENDIF
*ELSE
CMSEL,S,_Y
CMDELE,_Y
*ENDIF

allsel
keyopt,4,9,0
keyopt,4,12,2

mp,mu,1,0.2
mp,mu,2,0.2
mp,mu,3,0.2
mp,mu,4,0.2
mp,mu,5,0.2
mp,mu,6,0.3

!!!!!!做六个接触对，每层土一个，其中第一层土作了两个，桩的底面单独作了一个接触对

R,10, , ,5,0.1, , !应该设置最大接触摩擦,需要计算
RMORE, , ,1000000, , ,
real,10
vsel,s,,,12
vplot
asel,s,,,61,89,28
aplot
type,3
nsla,s,1
esurf,all
esel,s,type,,3
！eplot
allsel
vsel,s,,,1,6,5
vplot
asel,s,,,2,32,30
aplot
type,4
nsla,s,1
esurf,all
esel,s,type,,4
！eplot
allsel

r,11, , ,5,0.1, , !应该设置最大接触摩擦,需要计算
RMORE, , ,1000000, , ,
real,11
asel,s,,,66
aplot
type,3
nsla,s,1
esurf,all
esel,s,type,,3
eplot
allsel
asel,s,,,5
aplot
type,4
nsla,s,1
esurf,all
esel,s,type,,4
eplot
allsel

r,12, , ,5,0.1, , !应该设置最大接触摩擦,需要计算
RMORE, , ,1000000, , ,
real,12
asel,s,,,90
aplot
type,3
nsla,s,1
esurf,all
esel,s,type,,3
eplot
allsel
asel,s,,,16
aplot
type,4
nsla,s,1
esurf,all
esel,s,type,,4
！eplot
allsel

r,13, , ,5,0.1, , !应该设置最大接触摩擦,需要计算
RMORE, , ,1000000, , ,
real,13
asel,s,,,95
aplot
type,3
nsla,s,1
esurf,all
esel,s,type,,3
eplot
allsel
asel,s,,,20
aplot
type,4
nsla,s,1
esurf,all
esel,s,type,,4
！eplot
allsel

r,14, , ,5,0.1, , !应该设置最大接触摩擦,需要计算
RMORE, , ,1000000, , ,
real,14
asel,s,,,98
aplot
type,3
nsla,s,1
esurf,all
esel,s,type,,3
！eplot
allsel
asel,s,,,24
aplot
type,4
nsla,s,1
esurf,all
esel,s,type,,4
！eplot
allsel

r,15, , ,5,0.1, , !应该设置最大接触摩擦,需要计算
RMORE, , ,1000000, , ,
real,15
asel,s,,,80
aplot
type,3
nsla,s,1
esurf,all
esel,s,type,,3
！eplot
allsel
asel,s,,,60
aplot
type,4
nsla,s,1
esurf,all
esel,s,type,,4
！eplot
allsel

/solu
asel,s,loc,x,0
aplot
da,all,symm
asel,s,loc,x,2
aplot
da,all,all
asel,s,loc,z,0
aplot
da,all,symm
asel,s,loc,z,-2
da,all,all
asel,s,loc,y,-30
da,all,all
allsel

ACEL,0,9.8,0
csys,5
nsel,s,loc,z,0
nsel,r,loc,x,0.01,0.34
esln,s
eplot
csys,0
sfe,all,3,pres,,1000

nropt,full,,off ！！ 全牛顿-拉普拉斯解选项
antype,static ！！ 分析类型
!nlgeom,on ！！ 打开大变形
time,1
autot,on ! ！ 自动载荷步打开
!nsub,400
pred,on ！！ 预测求解器
outres,all,all ！！ 记录所有的结果文件
allsel
solv
finish