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

set log on f4r)g2Zb[  
set logfile mgan.log Xhm)K3RA*T  
new RoeLf Ow  
plot block group /eNDv(g)M  
plot show qASV\ <n  
generate zone radcyl p0 = (0,0,0) p1 = (20,0,0) &  njg\y  
p2 = (0,1,0) p3 = (0,0,20) p4 = (20,1,0) p5 = (0,1,20) & rhA>;9\  
p6 = (20,0,20) p7 = (20,1,20) p8 = (2.4,0,0) p9 = (0,0,2.4) p10 = (2.4,1,0) & I18<brZJ  
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 tA]Y=U+Q  
gen zone radtunnel p0 = (0,0,0) p1 = (0,0,-20) & i
!dQ Sdf  
p2 = (0,1,0) p3 = (20,0,0) & ".Sa[A;~  
size 10 5 10 10 dimension 2.2,2.4,2.2 fill group tunnel2 T
xXX}6  
gen zone reflect normal (-1,0,0) origin =(0,0,0) ;J&9l >  
group section1 range y 0 1 group tunnel any group tunnel2 any _omz74   
save model.sav Ul%D}(,  
pau \2NT7^H#  
;德鲁克-普拉格模型的材料参数 P*.0kR1n  
model dr Y[Kpd[)[v  
prop bulk 1.41e9 shear 0.97e9  kshear 1.603e7 qvol 0.917 ]d -U
  
;给定边界条件 `}|$eF&  
fix x range x -20.1 -19.9 fs6% M]u  
fix x range x 19.9 20.1 ]Wdnr1d~8  
fix z range z -20.1 -19.9 o`<ps$yT  
fix y range y -0.1 0.1 r>3y8
7  
fix y range y 0.9 1.1 1@{qPmf^  
;初始应力赋值 ewOR
b  
set grav 0 0 -10 _1kcz]]F  
ini density 2750 jRYW3a_7  
ini szz -17.6e6 grad 0 0 27500 Xgyi}~AoaU  
ini sxx -17.6e6 grad 0 0 27500 z]bcg$m  
ini syy -17.6e6 grad 0 0 27500 Gfy9?sa  
hist unbal ?)L
X4GY  
solve 7o4B1YD  
save chushihua1_maogan.sav pA?2UZ  
ini xvel 0 yvel 0 zvel 0 +je{%,*  
ini xdis 0 ydis 0 zdis 0 @]xHt&j  
mod null range group section1 J{h?=vK  
solve CwQRHi  
save kaiwa4.sav S^*ME*DDz  
;第一排 3KN>t)A#  
sel cable id=1 begin 0 0 2.4 end 0 0 4.9 nseg 25     ]B>g~t5J  
sel cable id=1 begin 0.92 0 2.2 end 1.88 0 4.53 nseg 25 &JHqUVs^  
sel cable id=1 begin 1.7 0 1.7 end 3.46 0 3.46 nseg 25 2V)qnMxAZJ  
sel cable id=1 begin 2.2 0 0.92 end 4.53 0 1.88 nseg 25  j2%?-(U  
sel cable id=1 begin 2.4 0 0 end 4.9 0 0 nseg 25 i*2l4  
sel cable id=1 begin 2.4 0 -0.8 end 4.9 0 -0.8 nseg 25 ~fR-cXj"  
sel cable id=1 begin 2.4 0 -1.6 end 4.9 0 -1.6 nseg 25 /XZ\Yy=  
sel cable id=1 begin -0.92 0 2.2 end -1.88 0 4.53 nseg 25 ?fmW'vs  
sel cable id=1 begin -1.7 0 1.7 end -3.46 0 3.46 nseg 25 Ze-M
B0w  
sel cable id=1 begin -2.2 0 0.92 end -4.53 0 1.87 nseg 25 r"\g6<RP  
sel cable id=1 begin -2.4 0 0 end -4.9 0 0 nseg 25 {u{8QKeC  
sel cable id=1 begin -2.4 0 -0.8 end -4.9 0 -0.8 nseg 25 ZtH
{2j0  
sel cable id=1 begin -2.4 0 -1.6 end -4.9 0 -1.6 nseg 25 qF57T>v|  
sel cable pro emod 2.06e11 xcarea 1.57e-3 gr_per 1. & )9'Zb`n  
yten 25e4 gr_k 17.5e6 gr_c 20e4 range id 1 3~6,fTMz{  
solve d gRTV<vM  
save maogan4.save   4VrL@c @  
P[<EFjE  
UcxMA%Pw7$  
请大家帮忙看看，我这个程序为什么运行了，支护没效果，谢谢！！