我建的桩土模型,算出来结果不理想,还请各位高手指点一二。我的模型不设接触面,采用线弹性模型,桩土接触地方采用捆绑法,即节点共用,现在结果出现问题,请各位大神帮忙解决,不胜感激!
命令流如下:
/title,Single pile
/nopr
keyw,pr_struc,1
/prep7
et,1,plane42
et,2,solid45
mp,ex,1,2e10
mp,nuxy,1,0.2
mp,dens,1,2500
mp,ex,2,2e8
mp,nuxy,2,0.3
mp,dens,1,2000
cyl4,,,0.5,,,90
cyl4,,,10,,,90
aptn,all
lsel,s,,,1,3,1
lesize,all,,,4
lsel,s,,,7,8
lesize,all,,,10,5
lsel,s,,,4
lesize,all,,,4
lsel,all
meshape,0,2d
mshkey,1
amesh,all
extopt,aclear,1
type,2
mat,1
esize,,30
extopt,aclear,1
voffst,1,15
esize,,2
voffst,1,-1
mat,2
esize,,30
extopt,aclear,1
voffst,3,15
nummrg,all
wpoff,0,0,15
cyl4,,,0.5,,,90
cyl4,,,10,,,90
aptn,14,16
lsel,s,,,25,26,1
lesize,all,,,4
lsel,s,,,21
lesize,all,,,4
lsel,s,,,30,31
lesize,all,,,10 ,5
lsel,s,,,27
lesize,all,,,4
lsel,all
meshape,0,2d
mshkey,1
amesh,14,17,3
esize,,30
extopt,aclear,1
voffst,14,15
voffst,17,15
mat,2
vsel,s,,,1,2,1
vplot
eslv,s
nslv,s,1
eplot
nummrg,all
numcmp,all
vsel,s,,,3,5,1
vplot
eslv,s
nslv,s,1
eplot
nummrg,all
numcmp,all
wpoff,0,0,-15
allsel
eplot
vsymm,x,all
vsymm,y,all
nummrg,all
numcmp,all
csys,1
asel,s,loc,x,10
aplot
da,all,ux
da,all,uy
da,all,uz
asel,s,loc,z,30
aplot
da,all,ux
da,all,uy
da,all,uz
acel,0,0,9.8
asel,s,loc,z,-1
aplot
/psf,pres,norm,3
sfa,all,1,pres,200000
allsel,all
finish
/solu
antype,static
nlgeom,on
autots,on
solve
finish
/post1
/pnum,area,1
/view,1,1,1,1
esel,s,mat,,1
plnsol,s,z
plnsol,u,z
esel,s,mat,,1
path,p1,2
ppath,1,,0,0,0
ppath,2,,0,0,15
pdef,sz,s,z
plpath,sz
pdef,uz,u,z
plpath,uz
|