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

set log on %IE;'aa }  
set logfile mgan.log XA(.O|VZ  
new +> d;%K  
plot block group ZHUW1:qs  
plot show SSTn|  
generate zone radcyl p0 = (0,0,0) p1 = (20,0,0) & GN Ewq$  
p2 = (0,1,0) p3 = (0,0,20) p4 = (20,1,0) p5 = (0,1,20) & B/qN1D]U.  
p6 = (20,0,20) p7 = (20,1,20) p8 = (2.4,0,0) p9 = (0,0,2.4) p10 = (2.4,1,0) & <RpTk*Yo^=  
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 U;*t5l  
gen zone radtunnel p0 = (0,0,0) p1 = (0,0,-20) & )siWc_Z4  
p2 = (0,1,0) p3 = (20,0,0) & 10.ZBfn  
size 10 5 10 10 dimension 2.2,2.4,2.2 fill group tunnel2 a7uL{*ZR  
gen zone reflect normal (-1,0,0) origin =(0,0,0) !`4ie  
group section1 range y 0 1 group tunnel any group tunnel2 any Hz[1c4)'F  
save model.sav a(*"r:/lD  
pau $tXW/  
;德鲁克-普拉格模型的材料参数 G<k.d"<  
model dr LFqY2,#i  
prop bulk 1.41e9 shear 0.97e9  kshear 1.603e7 qvol 0.917 nqx0#_K-E  
;给定边界条件 fAz4>_4  
fix x range x -20.1 -19.9 JiO8EIM  
fix x range x 19.9 20.1 q*,HN(&l?  
fix z range z -20.1 -19.9 WNb2"W  
fix y range y -0.1 0.1 n/ :#:  
fix y range y 0.9 1.1 Wd+G)Mu_=  
;初始应力赋值 R,["w98a  
set grav 0 0 -10 ZMr[:,Jp  
ini density 2750 6d6Dk>(V  
ini szz -17.6e6 grad 0 0 27500 ZXkrFA |  
ini sxx -17.6e6 grad 0 0 27500 Gc*=n*@^K  
ini syy -17.6e6 grad 0 0 27500 |pbetA4&  
hist unbal .Ln98#ZR  
solve ;K|K]c  
save chushihua1_maogan.sav QVkji7)ZT  
ini xvel 0 yvel 0 zvel 0 #G?#ot2o  
ini xdis 0 ydis 0 zdis 0 BNA1"@9q  
mod null range group section1 g.blDOmlc  
solve lq]8zm<\)]  
save kaiwa4.sav __s'/6u  
;第一排 ?(Tin80=r  
sel cable id=1 begin 0 0 2.4 end 0 0 4.9 nseg 25     )F\
kGe  
sel cable id=1 begin 0.92 0 2.2 end 1.88 0 4.53 nseg 25 >?)_, KL  
sel cable id=1 begin 1.7 0 1.7 end 3.46 0 3.46 nseg 25 pDJN}XtjT  
sel cable id=1 begin 2.2 0 0.92 end 4.53 0 1.88 nseg 25 ,quUGS  
sel cable id=1 begin 2.4 0 0 end 4.9 0 0 nseg 25 e-UWbn'~  
sel cable id=1 begin 2.4 0 -0.8 end 4.9 0 -0.8 nseg 25  W
/u(9  
sel cable id=1 begin 2.4 0 -1.6 end 4.9 0 -1.6 nseg 25 \CE+P5  
sel cable id=1 begin -0.92 0 2.2 end -1.88 0 4.53 nseg 25 {Y=k`t,  
sel cable id=1 begin -1.7 0 1.7 end -3.46 0 3.46 nseg 25 #`=>Mza  
sel cable id=1 begin -2.2 0 0.92 end -4.53 0 1.87 nseg 25 0-aaLC~Z>  
sel cable id=1 begin -2.4 0 0 end -4.9 0 0 nseg 25 :?= 1aiS  
sel cable id=1 begin -2.4 0 -0.8 end -4.9 0 -0.8 nseg 25 B0p;Zh  
sel cable id=1 begin -2.4 0 -1.6 end -4.9 0 -1.6 nseg 25 5D?{dA:Rq  
sel cable pro emod 2.06e11 xcarea 1.57e-3 gr_per 1. & C+ Y;D:  
yten 25e4 gr_k 17.5e6 gr_c 20e4 range id 1 mA%}ijR6y  
solve [`E_/95  
save maogan4.save   Mae2L2vc  
bivo7_  
y!
:vX6l  
请大家帮忙看看，我这个程序为什么运行了，支护没效果，谢谢！！