fini
/clear
/filname,dig
/title,element birth and death-dig
/prep7
et,1,42,,,1
mp,ex,1,2.675e9 !土的材料特性
mp,nuxy,1,0.3375
mp,dens,1,1000
mp,ex,2,2.85e10 !混凝土的材料特性
mp,nuxy,2,0.17
mp,dens,2,2548
rectng,0,7.5,0,4,
rectng,0,7.5,4,4+11/3,
rectng,0,7.5,4+11/3,4+11*2/3,
rectng,0,7.5,4+11*2/3,15,
rectng,4,4.5,0,15,
aptn,all !分割
/pnum,line,1
/pnum,area,1
/number,0 !不同的颜色和数字标记选择的项目
lplot
lsel,s,loc,x,2 !选择中心坐标为x=2的实体
lesize,all,,,4,,,,,1
lsel,s,loc,x,4.25
lesize,all,,,2,,,,,1
lsel,s,loc,x,6
lesize,all,,,3,,,,,1
lsel,s,loc,x,0
lsel,a,loc,x,4
lsel,a,loc,x,4.5
lsel,a,loc,x,7.5
lesize,all,,,3,,,,,1
alls
mshape,0,2d !2维四边形网格
mshkey,1 !映射网格
amesh,all
finish
save
!计算挖掘前的地基状态
/solu
antype,0 !静力学选项
nlgeom,on !打开大变形开关
nropt,full
outres,all,all
time,1
lsel,s,loc,y,0
dl,all,,uy,
sfgrad,pres,0,y,15,-3e4 !土压力荷载斜率，从Y=15开始沿Y轴施加每单位-3E4荷载
lsel,s,loc,x,7.5
sfl,all,pres,0, !周围土压力
acel,,9.8
allsel,all
sbctran !将几何实体上的荷载和边界条件转化为节点的荷载
/psf,pres,norm,2,0, !以箭头方式显示垂直面的方向的可见的荷载
eplot
save
solve
!计算地基开挖和添加混凝土围衬过程基坑分三次开挖，每次开挖后立即浇筑混凝土围衬
!第一次开挖和浇筑混凝土围衬
asel,s,,,9
asel,a,,,15
allsel,below,area !选择面9和面15面以下的实体，即面9和面15中的面，线，点实体全部选择
eplot
cm,e_death1,elem
ekill,all
allsel,all
time,2
solve
asel,s,,,9
allsel,below,area
eplot
mpchg,2,all
cm,e_birth1,elem
ealive,all
allsel,all
time,3
solve
!第二次开挖和浇筑混凝土围衬
asel,s,,,17
asel,a,,,13
allsel,below,area
eplot
cm,e_death2,elem
ekill,all
allsel,all
time,4
solve
asel,s,,,17
allsel,below,area
eplot
mpchg,2,all
cm,e_birth2,elem
ealive,all
allsel,all
time,5
solve
!第三次开挖和浇筑混凝土围衬
asel,s,,,16
asel,a,,,11
allsel,below,area
eplot
cm,e_death3,elem
ekill,all
allselall
time,6
solve
asel,s,,,16
allsel,below,area
eplot
mpchg,2,all
cm,e_birth3,elem
ealive,all
allsel,all
time,7
solve
fini
/post1
/expand,18,axis,,,10 !以Z为轴，转动XY平面，每10度1转，转18次。
!————制作总位移动画过程--------------
/show,dig_sum,grph
set,,,,,,,1 !读入数据组1，即读入开挖前（time=1)的数据组，
plnsol,u,sum,0,1 !以1倍的缩放因子显示位移的矢量和（即X,Y,Z三方向总位移），且仅显示结构变形图
cmsel,u,e_death1 !反选被开挖的e_birth1部分
set,,,,,,,2 !读入数据组2，即读入开挖第一部分后的数据组(time=2)
plnsol,u,sum,0,1
cmsel,a,e_birth1
set,,,,,,,3
plnsol,u,sum,0,1
cmsel,u,e_death2
set,,,,,,,4
plnsol,u,sum,0,1
cmsel,a,e_birth2
set,,,,,,,5
plnsol,u,sum,0,1
cmsel,u,e_death3
set,,,,,,,6
plnsol,u,sum,0,1
cmsel,a,e_birth3
set,,,,,,,7
plnsol,u,sum,0,1
/show,term !定向到上次指定的图形设备
!----------------制作环向应力动画过程——-----------------------
alls
/show,dig_SZ,grph
set,,,,,,,1
plnsol,s,z,0,1
cmsel,u,e_death1
set,,,,,,,2
plnsol,s,z,0,1
cmsel,a,e_birth1
set,,,,,,,3
plnsol,s,z,0,1
cmsel,u,e_death2
set,,,,,,,4
plnsol,s,z,0,1
cmsel,a,e_birth2
set,,,,,,,5
plnsol,s,z,0,1
cmsel,u,e_death3
set,,,,,,,6
plnsol,s,z,0,1
cmsel,a,e_birth3
set,,,,,,,7
plnsol,s,z,0,1
/show,term
