我没有翻译,但是有单元问题求助!在建立由二维拉伸成三维时,老是有问题,哪位高手帮忙解决一下,十分感谢!我的qq是765010043.
/clear
/prep7
et,1,solid45
et,2,mesh200,6
!土体材料参数
mp,ex,1,3.94e6
mp,prxy,1,0.35
mp,dens,1,1828
mp,ex,2,20.6e6
mp,prxy,2,0.30
mp,dens,2,2160
mp,ex,3,500e6
mp,prxy,3,0.33
mp,dens,3,2160
!~管片材料参数
mp,ex,4,27.6e9
mp,prxy,4,0.2
mp,dens,4,2500
mp,ex,5,1e9
mp,prxy,5,0.2
mp,dens,5,2100
!建立平面内模型
local,11,0,0,0,0
local,12,1,0,0,0
csys,11
wpcsys,-1
cyl4,,,,,2.7,90
cyl4,0,0,2.7,0,3,90
cyl4,0,0,3,0,3.2,90
rectng,4.5,0,4.5
aovlap,all
nummrg,all
numcpm,all
rectng,4.5,31.5,0,4.5
nummrg,all
numcpm,all
!划分网格
mshkey,1
type,2
lesize,1,,,6
lesize,2,,,6
lesize,3,,,6
amesh,1
lesize,4,,,6
lesize,8,,,2
lesize,9,,,2
amesh,2
lesize,5,,,6
lesize,10,,,1
lesize,11,,,1
amesh,3
lesize,12,,,3
lesize,13,,,3
lesize,6,,,3
lesize,7,,,3
lesize,14,,,8,2
lesize,16,,,8,0.5
amap,4,7,6,8,10
amap,5,9,8,11,12
!利用对称性
arsym,y,all
nummrg,all
numcpm,all
allsel,all
!建立下方土层
rectng,0,4.5,-4.5,-26
rectng,4.5,31.5,-4.5,-26
nummrg,all
numcpm,all
lesize,28,,,3
lesize,29,,,5,0.5
lesize,30,,,5,2
lesize,32,,,5,0.5
lesize,31,,,8,2
amesh,11
amesh,12
!建立上方土层
rectng,0,4.5,4.5,15
rectng,4.5,31.5,4.5,15
nummrg,all
numcpm,all
lesize,34,,,3
lesize,35,,,4
lesize,33,,,4
lesize,36,,,4
lesize,37,,,8,0.5
amesh,13
amesh,14
nummrg,all
numcpm,all
save
!利用对称性建立平面模型
allsel,all
arsym,x,all
nummrg,all
numcpm,all
save
!建立体模型
*do,i,1,20
k,,0,0,-3*i
*enddo
l,3,38
*do,i,1,19
l,37+i,38+i
*enddo
*do,i,1,20
lesize,64+i,,,1
*enddo
type,1
mat,1
esys,11
extopt,aclear,1
csys,12
*do,i,1,20
asel,s,loc,z,-3*(i-1)
vdrag,all,,,,,,64+i
*enddo
allsel,all
nummrg,all
numcpm,all
allsel,all
save
!分层赋予参数
csys,0
nsel,s,loc,y,-10.5,5
esln,s
mpchg,2,all
nsel,s,loc,y,-11,-26
esln,s
mpchg,3,all
!加载求解
Finish
/solu
csys,0
/view,1,-1,0.5,1
/replot
eplot
save
nsubst,10
nlgeom,on
nropt,full
pred,on
lnsrch,on
outres,all,last
nsel,s,loc,x,-31.5
d,all,ux
nsel,s,loc,x,31.5
d,all,ux
nsel,s,loc,y,-26
d,all,uy
nsel,s,loc,z,0
nsel,a,loc,z,-60
d,all,uz
allsel,all
acel,0,10,0
step=0
save
!求解自重应力场
time,1
solve
save
finish
/post1
PLNSOL,S,Y,0,1
PLNSOL,U,Y,0,1 |