[实验技术]利用ANSYS进行煤层开采 [复制链接]

finish
/clear
/prep7
et,1,42
mp,ex,1,2e10
mp,prxy,1,0.3
mp,dens,1,2500
mp,ex,2,2e10
mp,prxy,2,0.3
mp,dens,2,2500
rectng,0,40,0,40
rectng,10,30,10,12
aovlap,all
numcmp,all
/pnum,area,1
/pnum,line,1
aplot
lsel,s,,,1,7,2
lsel,a,,,2,8,2
lesize,all,2
mat,1
amesh,2
mat,2
amesh,1
/solu
nsel,s,loc,x,0
nsel,a,loc,x,40
d,all,ux,0
nsel,s,loc,y,0
d,all,uy,0
allsel
acel,,10
deltim,0.02,0.01,0.2
autots,on
pred,on
lnsrch,on
nlgeom,on
nropt,full
cnvtol,f,,0.02,2,0.5
solve
/solu
esel,s,mat,,2
ekill,all
esel,all
esel,s,live
eplot,all
solve
/post1
plnsol,u,y,0,1

以上是我设计的一个煤层开采的ansys程序，但不知哪里出现了问题，采空区不下沉，希望大家帮我分析一下原因……非常感谢

这是两步模拟工况啊，首先模拟自重应力场，之后将拟采部分冻结，再进行变形场分析

后开挖
finish
/clear
/prep7
et,1,42
mp,ex,1,2e10
mp,prxy,1,0.3
mp,dens,1,2500
mp,ex,2,2e10
mp,prxy,2,0.3
mp,dens,2,2500
rectng,0,40,0,40
rectng,10,30,10,12
aovlap,all
numcmp,all
/pnum,area,1
/pnum,line,1
aplot
lsel,s,,,1,7,2
lsel,a,,,2,8,2
lesize,all,2
mat,1
amesh,2
mat,2
amesh,1
/solu
nsel,s,loc,x,0
nsel,a,loc,x,40
d,all,ux,0
nsel,s,loc,y,0
d,all,uy,0
allsel
acel,,10
time,1
deltim,0.02,0.01,0.2
autots,on
pred,on
lnsrch,on
nlgeom,on
nropt,full
cnvtol,f,,0.02,2,0.5
solve
finish
/solu
allsel
esel,s,mat,,2
ekill,all
esel,all
esel,s,live
eplot
solve



先开挖
finish
/clear
/prep7
et,1,42
mp,ex,1,2e10
mp,prxy,1,0.3
mp,dens,1,2500
rectng,0,40,0,40
rectng,10,30,10,12
asba,1,2
/pnum,area,1
/pnum,line,1
aplot
lsel,s,,,1,7,2
lsel,a,,,2,8,2
lesize,all,2
mat,1
amesh,3
/solu
nsel,s,loc,x,0
nsel,a,loc,x,40
d,all,ux,0
nsel,s,loc,y,0
d,all,uy,0
allsel
acel,,10
solve

两个模型的本质是不一样的，开挖的顺序不一样，但为什么处理出来的结果，也就是Y方向的唯一变化量是一样的呢？谢谢了……