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

set log on \`-xxhb?e  
set logfile mgan.log [g=4'4EZc  
new `dl^)
4J  
plot block group a\B?J  
plot show `nc=@" 1  
generate zone radcyl p0 = (0,0,0) p1 = (20,0,0) & CE| *&G  
p2 = (0,1,0) p3 = (0,0,20) p4 = (20,1,0) p5 = (0,1,20) & f Avh!g  
p6 = (20,0,20) p7 = (20,1,20) p8 = (2.4,0,0) p9 = (0,0,2.4) p10 = (2.4,1,0) & e#WASHZN  
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 UA4MtTp`  
gen zone radtunnel p0 = (0,0,0) p1 = (0,0,-20) & nS Vr,wU  
p2 = (0,1,0) p3 = (20,0,0) & zqRps8=  
size 10 5 10 10 dimension 2.2,2.4,2.2 fill group tunnel2 wX!0KxR/Z  
gen zone reflect normal (-1,0,0) origin =(0,0,0) e{^lD.E  
group section1 range y 0 1 group tunnel any group tunnel2 any 6!=q+sw/X  
save model.sav azRp4~2?  
pau -f?,%6(1  
;德鲁克-普拉格模型的材料参数 7$*x&We  
model dr rV*Ri~Vx  
prop bulk 1.41e9 shear 0.97e9  kshear 1.603e7 qvol 0.917 6.|[;>Km  
;给定边界条件 EQ"+G[j~x  
fix x range x -20.1 -19.9 ="AaC!E,W  
fix x range x 19.9 20.1 5E}]
U,$  
fix z range z -20.1 -19.9 sn'E}.uhXH  
fix y range y -0.1 0.1 X'xnJtk  
fix y range y 0.9 1.1 e YiqTWn:  
;初始应力赋值 ~{{7y]3M-  
set grav 0 0 -10 WP^wNi ~>  
ini density 2750 1DH P5q  
ini szz -17.6e6 grad 0 0 27500 3,Iu!KB  
ini sxx -17.6e6 grad 0 0 27500 ]7q|) S\  
ini syy -17.6e6 grad 0 0 27500 y6H`FFqK  
hist unbal Oz<#s{Z  
solve $9W,1wg  
save chushihua1_maogan.sav )d{fDwrx1  
ini xvel 0 yvel 0 zvel 0 <uo@k'   
ini xdis 0 ydis 0 zdis 0 yEkwdx5!(  
mod null range group section1 hdH3Jb_hl(  
solve tK6=F63e  
save kaiwa4.sav AMK(-=  
;第一排 FU)=+m  
sel cable id=1 begin 0 0 2.4 end 0 0 4.9 nseg 25     ih: XC  
sel cable id=1 begin 0.92 0 2.2 end 1.88 0 4.53 nseg 25 fW=eB'Sl  
sel cable id=1 begin 1.7 0 1.7 end 3.46 0 3.46 nseg 25 f$--y|=  
sel cable id=1 begin 2.2 0 0.92 end 4.53 0 1.88 nseg 25 oS<*\!&D  
sel cable id=1 begin 2.4 0 0 end 4.9 0 0 nseg 25 YmHu8H_Q  
sel cable id=1 begin 2.4 0 -0.8 end 4.9 0 -0.8 nseg 25 Sj viH  
sel cable id=1 begin 2.4 0 -1.6 end 4.9 0 -1.6 nseg 25 ^bLFY9hSC  
sel cable id=1 begin -0.92 0 2.2 end -1.88 0 4.53 nseg 25 |!CAxE0d$B  
sel cable id=1 begin -1.7 0 1.7 end -3.46 0 3.46 nseg 25 ,Z8)DC=  
sel cable id=1 begin -2.2 0 0.92 end -4.53 0 1.87 nseg 25 ROO@EQ#`Z  
sel cable id=1 begin -2.4 0 0 end -4.9 0 0 nseg 25 K :>O X  
sel cable id=1 begin -2.4 0 -0.8 end -4.9 0 -0.8 nseg 25 '{)Jhl47   
sel cable id=1 begin -2.4 0 -1.6 end -4.9 0 -1.6 nseg 25 M5N#xgR  
sel cable pro emod 2.06e11 xcarea 1.57e-3 gr_per 1. & ^3QJv{)Q  
yten 25e4 gr_k 17.5e6 gr_c 20e4 range id 1 t"vkd  
solve ,hp8b$  
save maogan4.save   HI[Pf%${  
.kB!',v\  
h&rZR
`g  
请大家帮忙看看，我这个程序为什么运行了，支护没效果，谢谢！！