我建个模型 巷道开挖后 为什么 输入plot block state shear 只有空模型没有塑性破坏区域啊希望那个大神帮忙解决一下 谢谢啦
这是我的命令：建立模型
new
;泥岩
gen zone brick p0 0 0 0 p1 45 0 0 p2 0 20 0 p3 0 0 10 &
p4 45 20 0 p5 0 20 10 p6 45 0 10 p7 45 20 10 &
size 45 20 10
group '泥岩'
;煤
gen zone brick p0 0 0 10 p1 45 0 10 p2 0 20 10 p3 0 0 20 &
p4 45 20 10 p5 0 20 20 p6 45 0 20 p7 45 20 20 &
size 90 40 20
group '煤' range group '泥岩' not
;粗砂岩
gen zone brick p0 0 0 20 p1 45 0 20 p2 0 20 20 p3 0 0 30 &
p4 45 20 20 p5 0 20 30 p6 45 0 30 p7 45 20 30 &
size 90 40 20
group '粗砂岩' range group '泥岩' not group '煤' not
;中砂岩
gen zone brick p0 0 0 30 p1 45 0 30 p2 0 20 30 p3 0 0 33 &
p4 45 20 30 p5 0 20 33 p6 45 0 33 p7 45 20 33 &
size 45 20 3
group '中砂岩' range group '泥岩' not group '煤' not group '粗砂岩' not
;粗砂岩
gen zone brick p0 0 0 33 p1 45 0 33 p2 0 20 33 p3 0 0 40 &
p4 45 20 33 p5 0 20 40 p6 45 0 40 p7 45 20 40 &
size 45 20 7
group '粗砂岩' range group '泥岩' not group '煤' not group '粗砂岩' not &
group '中砂岩' not
;固定边界
fix x range x -0.1 0.1
fix x range x 45.1 44.9
fix y range y -0.1 0.1
fix y range y 19.9 20.1
fix z range z -0.1 0.1
;plot gpfix red
;设置重力加速度
set grav 0 0 -10
;参数模型
model mohr
;泥岩
prop density 2580 bulk 4.89E+09 shear 1.87E+09 coh 4.26E+06 fric 40.78 tens 5.62E+06 &
range group '泥岩'
;煤
prop density 1350 bulk 2.66E+09 shear 6.32E+08 coh 3.51E+06 fric 38.55 tens 2.10E+06 &
range group '煤'
;粗砂岩
prop density 2360 bulk 6.09E+09 shear 2.18E+09 coh 6.56E+06 fric 50.2 tens 2.58E+06 &
range group '粗砂岩'
;中砂岩
prop density 2380 bulk 3.02E+09 shear 1.90E+09 coh 6.48E+06 fric 46.74 tens 2.58E+06 &
range group '中砂岩'
;粗砂岩
prop density 2450 bulk 6.09E+09 shear 2.18E+09 coh 6.56E+06 fric 50.2 tens 5.62E+06 &
range group '粗砂岩'
;施加应力
apply szz -3.6E6 range z 39.9 40.1
;原岩应力
ini szz -3.6E6 grad 0 0 2.2E4 range z 0 40
ini syy -1.55E6 grad 0 0 9.5E3 range z 0 40 ;侧压力系数k=v/(1-v)=0.43
ini sxx -1.55E6 grad 0 0 9.5E3 range z 0 40
model null range x 20 25 y 0 20 z 10 14
;plot bl pro bulk ;可选，显示和查看分组体积模量
solve elas ;按弹性求解法生成初始应力场
save 无支护.sav
