采矿模型，开挖后，模型未垮落，请高手解答

litingyyhs

采矿模型，开挖后，模型未垮落，下面是代码，请高手解答

new
ro 0.05
set delc off
set ovtol 1
bl 0,0 0,188 700,188 700,0
cr 0,30 700,30
cr 0,36 700,36
cr 0,38 700,38
cr 0,88 700,88


;底板 细砂岩30m
jregion id 1 0,0 0,30 700,30 700,0 delete
jset 0,0 700,0 0,0 15,0 0,0 range jregion 1
jset 90,0 15,0 15,0 35,0 0,0 range jregion 1
jset 90,0 15,0 15,0 35,0 17.5,15 range jregion 1

;煤层 6m
jregion id 3 0,30 0,36 700,36 700,30 delete
jset 0,0 700,0 0,0 3,0 0,30 range jregion 3
jset 90,0 3,0 3,0 50,0 0,30 range jregion 3
jset 90,0 3,0 3,0 50,0 25,33 range jregion 3  
;泥岩 2m
jregion id 4 0,36 0,38 700,38 700,36 delete
jset 0,0 700,0 0,0 2,0 0,36 range jregion 4
jset 90,0 2,0 2,0 10,0 0,36 range jregion 4
jset 90,0 2,0 2,0 10,0 5,38 range jregion 4
;细砂 50m
jregion id 5 0,38 0,88 700,88 700,38 delete
jset 0,0 700,0 0,0 5,0 0,38 range jregion 5
jset 90,0 5,0 5,0 15,0 0,38 range jregion 5
jset 90,0 5,0 5,0 15,0 7.5,43 range jregion 5
;泥岩细砂岩互层100
jregion id 6 0,88 0,188 700,188 700,88 delete
jset 0,0 700,0 0,0 10,0 0,88 range jregion 6
jset 90,0 10,0 10,0 20,0 0,88 range jregion 6
jset 90,0 10,0 10,0 20,0 10,98 range jregion 6

;
del area 0.1
gen edge 15 range jregion 1
gen edge 3 range jregion 3
gen edge 2 range jregion 4
gen edge 5 range jregion 5
gen edge 10 range jregion 6

;
;change cons=3 range 0,700 0,188
change jcons=2 range 0,700 0,188
;块体材料
;底板 细砂岩
change mat=1 range jregion 1
prop mat=1 dens=2750 bulk=2.94e10,shear=2.58e10,coh=4.25e6,ten=1.55e6 fric=38 ;底板

;一分层
change mat=3 range jregion 3
prop mat=3 dens=1430 bulk=6e10,shear=5.5e10,coh=1.57e6,ten=0.32e6 fric=28 ;6煤
;直接顶
change mat=4 range jregion 4
prop mat=4 dens=2500 bulk=1.55e9,shear=1.38e9,coh=2.32e6,ten=0.48e6 fric=2 ;泥岩
;老顶
change mat=5 range jregion 5
prop mat=5 dens=2610 bulk=3.69e9,shear=1.59e9,coh=2.85e6,ten=1.11e26 fric=2 ;细砂岩
;泥岩细砂岩互层
change mat=6 range jregion 6
prop mat=6 dens=2500 bulk=4.7e9,shear=2.7e9,coh=3.68e6,ten=1.12e6 fric=3 ;泥岩细砂岩互层


;节理材料
change jmat=1 range jregion 1
prop jmat=1 jkn=2e7  jks=1.5e7 jfri=36 jten=3.19e6 jcoh=3.19e6

change jmat=3 range jregion 3
prop jmat=3 jkn=8e7  jks=8e7  jfri=30 jten=0.1e6 jcoh=0.1e6
change jmat=4 range jregion 4
prop jmat=4 jkn=7e7  jks=7e7  jfri=30 jten=3.19e4 jcoh=3.19e4
change jmat=5 range jregion 5
prop jmat=5 jkn=7e7 jks=7e7 jfri=30 jten=3.19e4 jcoh=3.19e4
change jmat=6 range jregion 6
prop jmat=6 jkn=7e7  jks=7e7 jfri=30 jten=3.19e4 jcoh=3.19e4
;

set jcondf 2
set gravity 0,-10
;
insitu stress -5e5,0,-4.165e5 szz -4.5e6
bound stress 0 0 -4.165e5 range -0.01 700.01 187.99 188.01 ;y=188
bound xvel 0 range -0.1,0.1 0,188
bound xvel 0 range 699.9,700.1 0,188
bound yvel 0 range 0,700 -0.1,0.1

hist solve type 1
solve ratio 1e-5
save mx-20.sav

liluanfei

yw5sehserhes

liluanfei

hsjhahsradsrhaszhdf

1002091319

楼主不错嘛

喂喂小毛驴

这里面没有开挖命令吧，只有初始平衡的运行，刚开始学，有没有大佬再深入解释下

GSH123LLL

你这不是开挖 只是原岩应力计算

kylin1282869061

高手,可以学习以下了

yangyuntao

楼主会了吗，求教学