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

set log on |$Td-M^)  
set logfile mgan.log fI6F};I5}T  
new ^-DK<jZ^  
plot block group '?\Hm'8  
plot show !)GPI?{^5  
generate zone radcyl p0 = (0,0,0) p1 = (20,0,0) & xed$z  
p2 = (0,1,0) p3 = (0,0,20) p4 = (20,1,0) p5 = (0,1,20) & DGcd|>q  
p6 = (20,0,20) p7 = (20,1,20) p8 = (2.4,0,0) p9 = (0,0,2.4) p10 = (2.4,1,0) & bkb}M)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 V5=Injs*  
gen zone radtunnel p0 = (0,0,0) p1 = (0,0,-20) & {+!_; zzZ  
p2 = (0,1,0) p3 = (20,0,0) & <R2bz1!h.  
size 10 5 10 10 dimension 2.2,2.4,2.2 fill group tunnel2 q;rU}hAzG0  
gen zone reflect normal (-1,0,0) origin =(0,0,0) OA+W$
  
group section1 range y 0 1 group tunnel any group tunnel2 any ^VA)vLj@  
save model.sav d/e9LK  
pau t;'__">:q  
;德鲁克-普拉格模型的材料参数 7{6wNc  
model dr _v-sb(* J  
prop bulk 1.41e9 shear 0.97e9  kshear 1.603e7 qvol 0.917 il:""x7^y  
;给定边界条件 jsuQR  
fix x range x -20.1 -19.9 N3,EF1%  
fix x range x 19.9 20.1 IySlu^a  
fix z range z -20.1 -19.9 xaPTTa  
fix y range y -0.1 0.1 =uHTpHR  
fix y range y 0.9 1.1 1*XqwBV  
;初始应力赋值 7Ev~yY;N  
set grav 0 0 -10 H]cCyuCdH  
ini density 2750 d%WF
gf}  
ini szz -17.6e6 grad 0 0 27500 1UOFTI2S|  
ini sxx -17.6e6 grad 0 0 27500 mWZVO,t$  
ini syy -17.6e6 grad 0 0 27500 Gb"PMai  
hist unbal  A/9 wr  
solve 9rhz#w  
save chushihua1_maogan.sav 7JbN WN  
ini xvel 0 yvel 0 zvel 0 bp }~{]:b  
ini xdis 0 ydis 0 zdis 0 W*P/~U=  
mod null range group section1 17-K~ybc  
solve ,\VNs'j  
save kaiwa4.sav m
V-MJ$3r  
;第一排 3 Tt8#B  
sel cable id=1 begin 0 0 2.4 end 0 0 4.9 nseg 25    
[=Xvp z  
sel cable id=1 begin 0.92 0 2.2 end 1.88 0 4.53 nseg 25 k7j;'6  
sel cable id=1 begin 1.7 0 1.7 end 3.46 0 3.46 nseg 25 W_?S^>?l/  
sel cable id=1 begin 2.2 0 0.92 end 4.53 0 1.88 nseg 25 NS\'o )J  
sel cable id=1 begin 2.4 0 0 end 4.9 0 0 nseg 25 0'gJSrgNI  
sel cable id=1 begin 2.4 0 -0.8 end 4.9 0 -0.8 nseg 25 kM.zX|_  
sel cable id=1 begin 2.4 0 -1.6 end 4.9 0 -1.6 nseg 25 )9}z^+TH  
sel cable id=1 begin -0.92 0 2.2 end -1.88 0 4.53 nseg 25 /Z^+K  
sel cable id=1 begin -1.7 0 1.7 end -3.46 0 3.46 nseg 25 }RXm=ArN  
sel cable id=1 begin -2.2 0 0.92 end -4.53 0 1.87 nseg 25 Q~jUZ-qN  
sel cable id=1 begin -2.4 0 0 end -4.9 0 0 nseg 25 dme_I
vt  
sel cable id=1 begin -2.4 0 -0.8 end -4.9 0 -0.8 nseg 25 #6JG#!W  
sel cable id=1 begin -2.4 0 -1.6 end -4.9 0 -1.6 nseg 25 b(wiJ&t  
sel cable pro emod 2.06e11 xcarea 1.57e-3 gr_per 1. & /gxwp:&lY  
yten 25e4 gr_k 17.5e6 gr_c 20e4 range id 1 'i}Q R~pe  
solve L;Nm"[`  
save maogan4.save   [xHK^JP 8F  
C3|M\[*fp  
8ED6C"6  
请大家帮忙看看，我这个程序为什么运行了，支护没效果，谢谢！！