桩土接触类分析报错求助

n3ch

问题描述：
三种材料，金属，水泥土，混凝土三种材料，在金属套筒内壁抹一层水泥土，用混凝土柱剪切水泥土，问题类似桩土接触问题。下面是我的命令流（分三层材料，1/4建模），对混凝土桩施加一个应力或位移力都可以：

!接触分析
finish
/clear
/config,nres,20000
/prep7


et,1,solid45
et,2,solid95
et,3,170
et,4,174
keyopt,4,9,0
keyopt,4,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,ex,3,2.0e11  !钢筒的材料性质
mp,prxy,2,0.28


mp,mu,4,0.2 !桩土摩擦系数
mp,mu,5,0.05!筒土摩擦系数
r,1
r,2
r,3
r,4
r,5

！模型
!桩
cyl4,,,49.5,,50.5,90,120
vatt,1,1,1
esize,5
mshape,0
vmesh,all
！土
cyl4,,,50.5,,52.5,90,120
vatt,2,2,1
esize,1
mshape,0
vmesh,all

！筒
cyl4,,,52.5,,53.5,90,120
vatt,3,3,2
esize,5
mshape,0
vmesh,all
!定义目标面/接触面
asel,s,,,3
nsla,s,1
real,4$type,3$mat,4$esurf,all$allsel,all

asel,s,,,10
nsla,s,a
real,4$type,4$esurf,all$allsel,all

asel,s,,,9
nsla,s,1
real,5$type,3$mat,5$esurf,all$allsel,all

asel,s,,,16
nsla,s,1
real,5$type,4$esurf,all$allsel,all

csys,0
NUMMRG,all
Allsel,all

allsel,all
asel,s,loc,x,0
asel,a,loc,y,0
da,all,symm
/solu
!定义底面约束
local,13,1
csys,13
asel,s,loc,x,52.5
nsla,s
d,all,all
outres,all,all
csys,0
antype,static
ncnv,2
nlgeom,on
autots,off
save
acel,,,-9.8
time,1
nsubst,100
!施加荷载
asel,s,,,2
nsla,s,1
sfa,all,1,pres,0.001
lswrite,1

ALLSEL,all
LSSOLVE,1

finish
没法计算，求大神指教问题出在哪里？第一次做有接触的问题，对接触不懂，也可以帮忙讲解下

htbbzzg

你定义了 5 个实常数：
r,1
r,2
r,3
r,4
r,5
后面定义接触时使用了这些实常数。可是你 5 个定义实常数的命令都没有给出任何实际数据啊？

yangheng8855

钢护筒桩，看看