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

set log on BLMcvK\9  
set logfile mgan.log ~xvQ?c?-  
new OtQKDpJq  
plot block group O8dD
oP\F2  
plot show zq1&MXR)l  
generate zone radcyl p0 = (0,0,0) p1 = (20,0,0) & 8|7Tk[X1j  
p2 = (0,1,0) p3 = (0,0,20) p4 = (20,1,0) p5 = (0,1,20) & "OK(<x]3;>  
p6 = (20,0,20) p7 = (20,1,20) p8 = (2.4,0,0) p9 = (0,0,2.4) p10 = (2.4,1,0) & }[i35f[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 $V5Ol6@2  
gen zone radtunnel p0 = (0,0,0) p1 = (0,0,-20) & cs@5K$v  
p2 = (0,1,0) p3 = (20,0,0) & zWjGGTP~3&  
size 10 5 10 10 dimension 2.2,2.4,2.2 fill group tunnel2
D4
o?  
gen zone reflect normal (-1,0,0) origin =(0,0,0) +;;fw |/  
group section1 range y 0 1 group tunnel any group tunnel2 any Eda sGCo  
save model.sav :B=`^>RK  
pau ]/d2*#  
;德鲁克-普拉格模型的材料参数 df8rf8B-  
model dr ?0_i{BvN  
prop bulk 1.41e9 shear 0.97e9  kshear 1.603e7 qvol 0.917 2ix_,yTO  
;给定边界条件 jv_sRV  
fix x range x -20.1 -19.9 k5}i^^.  
fix x range x 19.9 20.1 +m>Kb edl  
fix z range z -20.1 -19.9 OQ[>s(`*{  
fix y range y -0.1 0.1 V?Nl%M[b  
fix y range y 0.9 1.1 lbv, jS  
;初始应力赋值 ]7xAL7x  
set grav 0 0 -10 (iWNvVGS  
ini density 2750 Gl1`Nx0  
ini szz -17.6e6 grad 0 0 27500 wlKL|N  
ini sxx -17.6e6 grad 0 0 27500 W%1fm/G0  
ini syy -17.6e6 grad 0 0 27500 Ho"FB|e  
hist unbal AHZ6  
solve &6}vvgz  
save chushihua1_maogan.sav K5w22L^=+  
ini xvel 0 yvel 0 zvel 0 ^i:%;oeG  
ini xdis 0 ydis 0 zdis 0 Hs-NP#I  
mod null range group section1 +d'h20  
solve RbrvY  
save kaiwa4.sav }#yRaIp  
;第一排 XUSfOf(  
sel cable id=1 begin 0 0 2.4 end 0 0 4.9 nseg 25     #Q$4EQB  
sel cable id=1 begin 0.92 0 2.2 end 1.88 0 4.53 nseg 25 PK+][.6H  
sel cable id=1 begin 1.7 0 1.7 end 3.46 0 3.46 nseg 25 P, S9gG9  
sel cable id=1 begin 2.2 0 0.92 end 4.53 0 1.88 nseg 25 9qqEr
~  
sel cable id=1 begin 2.4 0 0 end 4.9 0 0 nseg 25 8{d`N|k  
sel cable id=1 begin 2.4 0 -0.8 end 4.9 0 -0.8 nseg 25 ~i  &K,  
sel cable id=1 begin 2.4 0 -1.6 end 4.9 0 -1.6 nseg 25 .R-:vU880  
sel cable id=1 begin -0.92 0 2.2 end -1.88 0 4.53 nseg 25 H!45w;,I  
sel cable id=1 begin -1.7 0 1.7 end -3.46 0 3.46 nseg 25 SH`
"o  
sel cable id=1 begin -2.2 0 0.92 end -4.53 0 1.87 nseg 25 wh~~g qi9  
sel cable id=1 begin -2.4 0 0 end -4.9 0 0 nseg 25 *`l>1)B>  
sel cable id=1 begin -2.4 0 -0.8 end -4.9 0 -0.8 nseg 25 "
&2D6  
sel cable id=1 begin -2.4 0 -1.6 end -4.9 0 -1.6 nseg 25 hw1ZTD:Y  
sel cable pro emod 2.06e11 xcarea 1.57e-3 gr_per 1. & wT.V3G  
yten 25e4 gr_k 17.5e6 gr_c 20e4 range id 1  |u$AzI  
solve M#|TQa
N  
save maogan4.save   t*D[Q$v  
aa".d[*1  
[e&$4l IS  
请大家帮忙看看，我这个程序为什么运行了，支护没效果，谢谢！！