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

set log on 1gShV ]2  
set logfile mgan.log o\ow{gh9  
new 7lDaok  
plot block group t+!gzZ  
plot show \kV7NA  
generate zone radcyl p0 = (0,0,0) p1 = (20,0,0) & /-{O\7-D  
p2 = (0,1,0) p3 = (0,0,20) p4 = (20,1,0) p5 = (0,1,20) & ?PE1aB+{:  
p6 = (20,0,20) p7 = (20,1,20) p8 = (2.4,0,0) p9 = (0,0,2.4) p10 = (2.4,1,0) & df4^C->:  
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 ONN{4&7@<  
gen zone radtunnel p0 = (0,0,0) p1 = (0,0,-20) & ;}eEG{`Y  
p2 = (0,1,0) p3 = (20,0,0) & ,rFLpQl  
size 10 5 10 10 dimension 2.2,2.4,2.2 fill group tunnel2 A,lw-(.z4Z  
gen zone reflect normal (-1,0,0) origin =(0,0,0) |<3Q+EB^  
group section1 range y 0 1 group tunnel any group tunnel2 any ss`q{ARb  
save model.sav K;y\[2;}e,  
pau B#GZmv1  
;德鲁克-普拉格模型的材料参数 OpbT63@L  
model dr !qXq y}?w  
prop bulk 1.41e9 shear 0.97e9  kshear 1.603e7 qvol 0.917 ~I\r1Wj;  
;给定边界条件 GQ-e$D@SfB  
fix x range x -20.1 -19.9 O3C)N I\i  
fix x range x 19.9 20.1 [> &+*c  
fix z range z -20.1 -19.9 0Dm`Ek3A7x  
fix y range y -0.1 0.1 ?X_0Iy}1  
fix y range y 0.9 1.1 I0l.KiBm  
;初始应力赋值 hw ;dm  
set grav 0 0 -10 ccN&h  
ini density 2750 *T>#zR{  
ini szz -17.6e6 grad 0 0 27500 /cL9?k;o  
ini sxx -17.6e6 grad 0 0 27500 <Xs@ \  
ini syy -17.6e6 grad 0 0 27500 FJjF*2.  
hist unbal ?%dCU~ z  
solve Gp
^ owr  
save chushihua1_maogan.sav F
0BOhlK  
ini xvel 0 yvel 0 zvel 0 ;h-G3>Il  
ini xdis 0 ydis 0 zdis 0 p#;dLM/EA  
mod null range group section1 _N,KHxsG8B  
solve iTugvb
 
save kaiwa4.sav O5TK&j  
;第一排 U.pr} hq  
sel cable id=1 begin 0 0 2.4 end 0 0 4.9 nseg 25     h"BhTx7E}  
sel cable id=1 begin 0.92 0 2.2 end 1.88 0 4.53 nseg 25 @0UwI
%.  
sel cable id=1 begin 1.7 0 1.7 end 3.46 0 3.46 nseg 25 )1Ma~8Y%r  
sel cable id=1 begin 2.2 0 0.92 end 4.53 0 1.88 nseg 25 W}R=  
sel cable id=1 begin 2.4 0 0 end 4.9 0 0 nseg 25 oX2DFgz  
sel cable id=1 begin 2.4 0 -0.8 end 4.9 0 -0.8 nseg 25 +wz`_i)!  
sel cable id=1 begin 2.4 0 -1.6 end 4.9 0 -1.6 nseg 25 lYZ@a4TA  
sel cable id=1 begin -0.92 0 2.2 end -1.88 0 4.53 nseg 25 Nx4DC  
sel cable id=1 begin -1.7 0 1.7 end -3.46 0 3.46 nseg 25 GrLM${G  
sel cable id=1 begin -2.2 0 0.92 end -4.53 0 1.87 nseg 25 c;21i;&,9  
sel cable id=1 begin -2.4 0 0 end -4.9 0 0 nseg 25 j |'#5H`  
sel cable id=1 begin -2.4 0 -0.8 end -4.9 0 -0.8 nseg 25 `!,\kc1  
sel cable id=1 begin -2.4 0 -1.6 end -4.9 0 -1.6 nseg 25 @%G'U&R{  
sel cable pro emod 2.06e11 xcarea 1.57e-3 gr_per 1. & }GNkB  
yten 25e4 gr_k 17.5e6 gr_c 20e4 range id 1 mU?&\w=v$  
solve ZaRr2Z:!  
save maogan4.save   3\p]esse  
`z$uw  
p~,3A:i  
请大家帮忙看看，我这个程序为什么运行了，支护没效果，谢谢！！