自己做的巷道支护例子，有些不懂，请各位高手提点建议，万分感谢！
1、命令流是否存在问题？需做哪些改进？
2、锚杆支护后变形和不加锚杆时的变形差不多，原因何在？
再次表示感谢！
new
plot add block group
plot show
; 产生网格体
;巷道上部混凝土衬砌
gen zon cshell p0 0 0 0 p1 3.1 0 0 p2 0 50 0 p3 0 0 3.1 &
dim 2.95 2.95 2.95 2.95 size 2 50 10 &
group 'concrete liner1'
;巷道上部
gen zon cyl p0 0 0 0 p1 2.95 0 0 p2 0 50 0 p3 0 0 2.95 &
size 5 50 10 &
group tunnela
;; 巷道下部
gen zone brick p0 0 0 -1.6 p1 add 2.95 0 0 p2 add 0 50 0 &
p3 add 0 0 1.6 size 5 50 6 &
group tunnelb
;;巷道两帮混凝土衬砌
gen zone brick p0 2.95 0 -1.6 p1 add 0.15 0 0 p2 add 0 50 0 &
p3 add 0 0 1.6 size 2 50 6 &
group 'concrete liner2'
; 围岩
;;
gen zone brick p0 0 0 -9 p1 8.571 0 -9 p2 add 0 50 0 p3 0 0 -1.6 &
p4 8.571 50 -9 p5 0 50 -1.6 p6 2.95 0 -1.6 p7 2.95 50 -1.6 &
size 5 50 8 ratio 1 1 0.769 &
group '围岩1'
gen zone brick p0 8.571 0 -9 p1 12 0 -9 p2 add 0 50 0 p3 2.95 0 -1.6 &
p4 12 50 -9 p5 2.95 50 -1.6 p6 3.1 0 -1.6 p7 3.1 50 -1.6 &
size 2 50 8 ratio 1 1 0.769 &
group '围岩2'
gen z radc p0 0 0 0 p1 12 0 0 p2 0 50 0 p3 0 0 9 &
dim 3.1 3.1 3.1 3.1 size 5 50 10 8 rat 1 1 1 1.3 &
group '围岩3'
gen zone brick p0 3.1 0 -1.6 p1 12 0 -9 p2 add 0 50 0 p3 3.1 0 0 &
p4 12 50 -9 p5 3.1 50 0 p6 12 0 0 p7 12 50 0 &
size 8 50 6 ratio 1.3 1 1 &
group '围岩4'
gen zone brick p0 0 0 9 p1 12 0 9 p2 add 0 50 0 p3 add 0 0 10 &
size 5 50 4 ratio 1 1 1.2 &
group '围岩上1'
gen zone brick p0 12 0 9 p1 25 0 9 p2 add 0 50 0 p3 add 0 0 10 &
size 3 50 4 ratio 1.2 1 1.2 &
group '围岩上2'
gen zone brick p0 12 0 0 p1 25 0 0 p2 add 0 50 0 p3 add 0 0 9 &
size 3 50 5 ratio 1.2 1 1 &
group '围岩右1'
gen zone brick p0 12 0 -9 p1 25 0 -9 p2 add 0 50 0 p3 add 0 0 9 &
size 3 50 6 ratio 1.2 1 1 &
group '围岩右2'
gen zone brick p0 0 0 -21 p1 12 0 -21 p2 add 0 50 0 p3 add 0 0 12 &
size 7 50 4 ratio 1 1 0.769 &
group '围岩下1'
gen zone brick p0 12 0 -21 p1 25 0 -21 p2 add 0 50 0 p3 add 0 0 12 &
size 3 50 4 ratio 1.2 1 0.769 &
group '围岩下2'
gen merge 0.02
;; 指定 Mohr-Coulomb 材料模型
mo mohr
pro den 2200 bulk 2.68e10 she 7e9 fric 38 coh 3.8e6 ten 4.86e6 &
range z -21 -9
pro den 2300 bulk 1.57e10 she 1.078e10 fric 40 coh 3e6 ten 2.31e6 &
range z -9 9
pro den 2100 bulk 2.56e10 she 8.52e9 fric 34 coh 4e6 ten 1.29e6 &
range z 9 19
;; 边界条件
fix x range x -.1 .1
fix x range x 24.9 25.1
fix y range y -.1 .1
fix y range y 49.9 50.1
fix z range z -9.1 -8.9
;; 初始应力状态
set grav 0 0 -10
apply szz -23e6 range z 18.9 19.1
ini sxx -25e6
ini syy -30e6
ini szz -17e6
hist unbal
plot cont dis out on
so
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
ini sta 0
sav 'maogana.sav
;;;开挖
mo null range group tunnela
mo null range group tunnelb
;;安装锚杆长度为2400 800*1000 共51排
def place_cables
loop n (1,51)
y_d = float(n - 1)
command
sel cable id n beg 0.000,y_d,2.95 end 0.000,y_d,5.35 nseg 3
sel cable id n beg 0.793,y_d,2.841 end 1.438,y_d,5.153 nseg 3
sel cable id n beg 1.527,y_d,2.524 end 2.77,y_d,4.577 nseg 3
sel cable id n beg 2.149,y_d,2.021 end 3.898,y_d,3.664 nseg 3
sel cable id n beg 2.613,y_d,1.369 end 4.739,y_d,2.482 nseg 3
sel cable id n beg 2.885,y_d,0.616 end 5.232,y_d,1.117 nseg 3
sel cable id n beg 2.950,y_d,-0.182 end 5.340,y_d,-0.33 nseg 3
sel cable id n beg 2.95,y_d,-1.3 end 5.268,y_d,-1.921 nseg 3
sel cable id n pro emod 2.06e11 xcarea 1.57e-3 gr_per 1. &
yten 25e4 gr_k 17.5e6 gr_c 20e4
end_command
end_loop
end
place_cables
plot grid sel geom red
求解
监控变量
his nstep=5
his unb
hist gp zdisp 0 0 2.95
hist gp xdisp 0 4 2.95
hist gp zdisp 0 8 2.95
hist gp zdisp 0 12 2.95
hist gp zdisp 0 16 2.95
hist gp xdisp 0 20 2.95
hist gp zdisp 0 24 2.95
hist gp zdisp 0 28 2.95
hist gp zdisp 0 32 2.95
hist gp zdisp 0 36 2.95
hist gp zdisp 0 40 2.95
hist gp zdisp 0 44 2.95
hist gp zdisp 0 48 2.95
hist gp zdisp 0 50 2.95
plot cont dis out on
求解
set mech force 50
;设置大变形
set large
solve
save 'maoganb'.sav
