求助，FLAC程序，巷道锚杆支护 [复制链接]

set log on i#W*'  
set logfile mgan.log 5~+XZA#2  
new /L v1$~  
plot block group cp6WMHLj   
plot show ^Ga&
}-  
generate zone radcyl p0 = (0,0,0) p1 = (20,0,0) & !aNh!  
p2 = (0,1,0) p3 = (0,0,20) p4 = (20,1,0) p5 = (0,1,20) & sJZ!sznn  
p6 = (20,0,20) p7 = (20,1,20) p8 = (2.4,0,0) p9 = (0,0,2.4) p10 = (2.4,1,0) & G[Tl
%w  
p11 = (0,1,2.4) size 10 5 10 10 dimension 2.4 2.4 2.4 2.4 ratio 1 1 1 1 fill group tunnel s_ bR]G  
gen zone radtunnel p0 = (0,0,0) p1 = (0,0,-20) & \
+nGOvM  
p2 = (0,1,0) p3 = (20,0,0) & ~R]35Cp-#  
size 10 5 10 10 dimension 2.2,2.4,2.2 fill group tunnel2 4(cJ^]wb^  
gen zone reflect normal (-1,0,0) origin =(0,0,0) n~gLPHY  
group section1 range y 0 1 group tunnel any group tunnel2 any XOy2lJ/  
save model.sav OP! R[27>  
pau L{8;Ud_2r  
;德鲁克-普拉格模型的材料参数 vF4]ux&  
model dr H?
`g!cX  
prop bulk 1.41e9 shear 0.97e9  kshear 1.603e7 qvol 0.917 s B 20/F  
;给定边界条件 h#qN+qt}  
fix x range x -20.1 -19.9 1n=_y o  
fix x range x 19.9 20.1 a~]bD  
fix z range z -20.1 -19.9 9$1)k;ChP/  
fix y range y -0.1 0.1 ^,8R,S\}$  
fix y range y 0.9 1.1 e3[:D5  
;初始应力赋值 4n0xE[-  
set grav 0 0 -10 K@u&(}  
ini density 2750 r"{<%e  
ini szz -17.6e6 grad 0 0 27500 -YGbfd<wq  
ini sxx -17.6e6 grad 0 0 27500 6^)}PX= *  
ini syy -17.6e6 grad 0 0 27500 v;EQ, NL  
hist unbal b-^p1{A0zW  
solve `>lzlEhKV  
save chushihua1_maogan.sav jd-glE,Y/  
ini xvel 0 yvel 0 zvel 0 ZcE_f>KV  
ini xdis 0 ydis 0 zdis 0 )?aaBaN$  
mod null range group section1 %O[N}_XHEh  
solve h9s >LY  
save kaiwa4.sav O;9?(:_  
;第一排 EbY,N:LK  
sel cable id=1 begin 0 0 2.4 end 0 0 4.9 nseg 25     NjuiD].  
sel cable id=1 begin 0.92 0 2.2 end 1.88 0 4.53 nseg 25 ZWUP^
V  
sel cable id=1 begin 1.7 0 1.7 end 3.46 0 3.46 nseg 25 s5&@Cxzl  
sel cable id=1 begin 2.2 0 0.92 end 4.53 0 1.88 nseg 25 MCTsi:V>+  
sel cable id=1 begin 2.4 0 0 end 4.9 0 0 nseg 25 IE2"rQT  
sel cable id=1 begin 2.4 0 -0.8 end 4.9 0 -0.8 nseg 25 v-b0\_  
sel cable id=1 begin 2.4 0 -1.6 end 4.9 0 -1.6 nseg 25 B%s7bS  
sel cable id=1 begin -0.92 0 2.2 end -1.88 0 4.53 nseg 25 ]< l6s  
sel cable id=1 begin -1.7 0 1.7 end -3.46 0 3.46 nseg 25 L1{GL #qV  
sel cable id=1 begin -2.2 0 0.92 end -4.53 0 1.87 nseg 25 6xK[34~6  
sel cable id=1 begin -2.4 0 0 end -4.9 0 0 nseg 25 uQ1@b-e`5  
sel cable id=1 begin -2.4 0 -0.8 end -4.9 0 -0.8 nseg 25 7- d.ZG  
sel cable id=1 begin -2.4 0 -1.6 end -4.9 0 -1.6 nseg 25 }T2xXbU  
sel cable pro emod 2.06e11 xcarea 1.57e-3 gr_per 1. & dihjpI_  
yten 25e4 gr_k 17.5e6 gr_c 20e4 range id 1 Sjj>#}U  
solve l|2D/K5  
save maogan4.save   YY((#"o;l  
}L=/A7Nk>  
d*8 $>GA  
请大家帮忙看看，我这个程序为什么运行了，支护没效果，谢谢！！