问题描述:
三种材料,金属,水泥土,混凝土三种材料,在金属套筒内壁抹一层水泥土,用混凝土柱剪切水泥土,问题类似桩土接触问题。下面是我的命令流(分三层材料,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
没法计算,求大神指教问题出在哪里?第一次做有接触的问题,对接触不懂,也可以帮忙讲解下 |