fini
/cle
/nopr
keyw,pr_struc,1
/prep7
et,1,shell43
et,2,combin14
mp,ex,1,2.8e10
mp,nuxy,1,0.2
mp,dens,1,2500
r,1,0.5
r,2,30e6 !是否需要换算,用的是地基反应模量,单位N/m3
!!
k,1
k,2,5
k,3,5,5
k,4,,5
lstr,1,2
lstr,2,3
lstr,2,3
lstr,3,4
lstr,1,4
a,1,2,3,4
VEXT,1,,,,,10
vdele,1
adele,1,2,1
!!
lsel,s,line,,9,12,1
lesize,all,0.5
lsel,s,line,,1,8,1
lesize,all,0.5
allsel,all
type,1
mat,1
real,1
amesh,all
/pnum,elem,1
!!
nsel,s,loc,y,0
ngen,2,1000,all,,,,-1
nplot
allsel,all
/rep
!!
type,2
mat,1
real,2
*do,i,1,231,1
e,i,i+1000
*enddo
eplot
!!
/sol
nsel,s,loc,y,-1 !约束弹簧底端的所有自由度
d,all,all
allsel,all
nsel,s,loc,y,0 !我约束了箱涵底面的ux,uz两个方向,感觉多余了
d,all,ux
d,all,uz
allsel,all
!!
acel,0,10
sfa,5,,pres,,-63060 !箱涵顶面荷载,暂时按均布考虑,不过这个值
!是否包含顶板的混凝土自重?
allsel,all
sfgrad,pres,0,y,0,3720
sfa,4,,pres,-40000 !对这边施加的方向有点疑惑,怎么看到施加的效果呢?
sfgrad,pres,0,y,0,3720
sfa,6,,pres,-23400
sfa,3,,pres,-85560 !底板的反力,加了弹簧的话,是否还要加反力呢?
/psf,pres,norm,2,0,1
!!
nropt,full
allsel
outres,all,all
solve
fini
/post1
pldisp,1
etable,,smisc,4
pletab,smis4,-1,0
etable,,smisc,7
pletab,smis7,-1,0 !弯距的显示结果好像有很大问题
etable,,smisc,1
pletab,smis1,1,0
etable,,smisc,2
etable,,smisc,7
etable,,smisc,8
plls,smis1,smis2,1,0 !轴力
plls,smis27,smis8,1 !剪力
finish !下面考虑的是杀死所有弹簧单元分析的结果
/solu
allsel,all
esel,s,type,,2
ekill,all
allsel,all
/post1
pldisp,1
etable,,smisc,7
pletab,smis7,-1,0 ! 感觉没多大差别
etable,,smisc,1
pletab,smis1,1,0
etable,,smisc,2
etable,,smisc,7
etable,,smisc,8
plls,smis1,smis2,1,0 !轴力
plls,smis27,smis8,1 !剪力 |