求助：关于深基坑开挖

隗小宝

小弟刚学flac，需要做三维基坑模拟方面的毕业论文。看到帖子：3维基坑数值模拟现场直播（代码已经贴出），完工后的照片也贴出。http://forum.simwe.com/thread-565835-1-3.html。非常需要上面的代码，请问哪位仁兄有上面的命令流，能传我一下吗？非常感谢，我的邮箱是jacklove2001@163.com

yangxin546

你给我QQ,我发给你

yangxin546

用ANSYS写的，贴给你，附件下载地址http://bbs.simwe.com/thread-911988-1-1.html
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

yangxin546

用ANSYS写的，贴给你，附件下载地址http://bbs.simwe.com/thread-911988-1-1.html
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

隗小宝

非常感谢，谢谢谢谢！！！！！！！！