[资料]锚杆支护程序 [复制链接]

Cb-E<W&2D  
我的这个模拟是基于这样一个问题：采矿中开挖4条同样支护方式的巷道，但是巷道的位置不同（煤柱的宽度不同） ~,1q :Kue  
;在其中使用了几重循环，可能对初学者有参考价值，我也常常看论坛的资料，现在自己贡献一点。 M\6v}kUY  
;这个文件可以直接运行的。 A>2p/iMc  
;建立网格 820m*900m*350m  网格51*30*30=4.59万网格 JU.%;e7  
gen zone brick size 10 30  4 & Bb"4^EOZ,  
p0   0 0 0  p1 295 0 0 p2 0 900 0   p3 0 0 40  ratio 0.8 1  1 $NRb'  
gen zone brick size 8 30  4 & #Kr.!uD  
p0 295 0 0  p1 335 0 0 p2 295 900 0 p3 295 0 40 E\N=p&g$  
gen zone brick size 15 30  4 & j]D =\  
p0 335 0 0  p1 485 0 0 p2 335 900 0 p3 335 0 40 e>Y2q|S85  
gen zone brick size 8 30  4 & ?0%TE\I8  
p0 485 0 0  p1 525 0 0 p2 485 900 0 p3 485 0 40 0l@+xS;  
gen zone brick size 10 30  4 & lM%fgyX  
p0 525 0 0  p1 820 0 0 p2 525 900 0 p3 525 0 40 ratio 1.25 1 1 }]?G"f t K  
gen zone brick size 10 30  10 & gQDK?aQX  
p0   0 0 40  p1 295 0 40 p2 0 900 40   p3 0 0 50  ratio 0.8 1  1 )fL*Ws6  
gen zone brick size 8 30  10 & o+Z9h1z%,  
p0 295 0 40  p1 335 0 40 p2 295 900 40 p3 295 0 50 e;[8GE.   
gen zone brick size 15 30  10 & ,LO-!\L  
p0 335 0 40  p1 485 0 40 p2 335 900 40 p3 335 0 50 I@M^Wu]wW  
gen zone brick size 8 30  10 & mcG$V0D <{  
p0 485 0 40  p1 525 0 40 p2 485 900 40 p3 485 0 50 ]*U')  
gen zone brick size 10 30  10 & %"^XxVJ*  
p0 525 0 40  p1 820 0 40 p2 525 900 40 p3 525 0 50 ratio 1.25 1 1 e.^9&Fk"N  
gen zone brick size 10 30  6 & 6|Q'
\  
p0   0 0 50  p1 295 0 50 p2 0 900 50   p3 0 0 110  ratio 0.8 1  1 ]<LU NxBR  
gen zone brick size 8 30  6 & .)zISa*Xy  
p0 295 0 50  p1 335 0 50 p2 295 900 50 p3 295 0 110 c3t8yifQ  
gen zone brick size 15 30  6 & "?,6{\y,  
p0 335 0 50  p1 485 0 50 p2 335 900 50 p3 335 0 110 (\>'yW{f  
gen zone brick size 8 30  6 & -Lb^O/  
p0 485 0 50  p1 525 0 50 p2 485 900 50 p3 485 0 110 Rw 8o]  
gen zone brick size 10 30  6 & 0M98y!A 5^  
p0 525 0 50  p1 820 0 50 p2 525 900 50 p3 525 0 110 ratio 1.25 1 1 a $%[!vF  
gen zone brick size 10 30  10 & loe>"_`Cq  
p0   0 0 110  p1 295 0 110 p2 0 900 110   p3 0 0 350  ratio 0.8 1 1.25 lM"7 Z  
gen zone brick size 8 30  10 & c`; LF'
!  
p0 295 0 110  p1 335 0 110 p2 295 900 110 p3 295 0 350  ratio 1 1 1.25 d vxEXy  
gen zone brick size 15 30  10 & mh]$g<*m  
p0 335 0 110  p1 485 0 110 p2 335 900 110 p3 335 0 350  ratio 1 1 1.25 r/2:O92E  
gen zone brick size 8 30  10 & `0D1Nh"%k  
p0 485 0 110  p1 525 0 110 p2 485 900 110 p3 485 0 350  ratio 1 1 1.25 7#3)&"j  
gen zone brick size 10 30  10 & 1z:N$O_v  
p0 525 0 110  p1 820 0 110 p2 525 900 110 p3 525 0 350  ratio 1.25 1 1.25 )c !S@Hs  
LL [>Uu?Y  
;本构关系、参数 e6'O,\  
m m Th^#H  
prop den=2700  bulk=15e9  shear=60e8   fric=42  coh=2.61e6   range z 0  44 i8.[d5  
prop den=1500  bulk=1e9   shear=4e8    fric=30  coh=1.4e6    range z 44 50   &Qjl|2  
prop den=2700  bulk=15e9  shear=60e8   fric=42  coh=2.61e6   range z 50 350 -P
&e4sV{  
;边界、初始条件 i`'^ zR(`i  
fix z range z  -0.1 0.1 H-w|JH>g  
fix x range x  -0.1 0.1 RY<b]|  
fix x range x  819.9 820.1 Uk6!Sb  
fix y range y  -0.1 0.1 ^W'[l al.  
fix y range y  899.9 900.1 FJ"9Hs2  
;重力梯度化 hspg-|R  
set grav 0 0 -10 KLW+&.re8  
ini szz -9.45e6  grad 0,0,2.7e4   range z 50 350 eMzCAO  
ini szz -8.85e6  grad 0,0,1.5e4   range z 44 50 &N0|tn  
ini szz -9.378e6 grad 0,0,2.7e4   range z 0 44 v{Vesf  
hist unbal ,ua1xsZl&  
plot set plane nor 0 1 0 origin 0 100 0 $(=1A>40  
plot con szz   plan ]H2aYi$  
step 3000 -~f.>@Wb  
hist unbal Y cpO;md  
save 4001.sav yFsXI0I[p  
;没开挖前的的自重平衡，为状态1 pnJT]?},  
;$< ek(i7  
rest 4001.sav ^ni_%`Ag  
;开挖巷道 4N j?UDa  
ini state=0 )7J>:9h  
ini xdis 0 ydis  0  zdis 0 5zOSb$;  
ini xvel 0 yvel  0  zvel 0 B,,d~\  
;以下是自己写的几重循环的嵌套，可能对初学者有用。 qH"a!  
;开挖四条巷道加锚杆支护的实验 -+|[0hpw  
；定义一个函数 aa_2 n`xh/vGm#  
def aa_2 E2D8s=r  
;定义一个循环变量 aaa,为了实现在不同位置开挖4条同样参数巷道，巷道位置循环，循环1 6QQ oHYtZ  
loop aaa (1,4) <vDm(-i3  
  caseof aaa p]=;t
"  
  case 1 GGtrH~zx
  
   x_d1=295 pSFWNWQ'B  
  case 2 lJ#>Y5Qg  
   x_d1=520 RC%r7K f  
             ;10m煤柱 U$uO%:4%  
  case 3 2zu~#qU[)M  
   x_d1=310 wgrOW]e  
  case 4 ArK9E!`^  
   x_d1=505 Lm#d.AD)  
  end_case kELyD(^P`  
;定义几个变量，在下面的循环中要用到 or`stBx  
;左帮锚杆的起、终点变量 ,c l<74d  
x_d11=x_d1-1.4 [{$0E=&0  
;右帮锚杆的起、终点变量 x$DJ  
x_d2=x_d1+5 V"iLeC  
x_d21=x_d2+1.4 |pSoBA9U  
;开挖巷道 IoOnS)  
       command GJPZ[bo  
   model null range  x x_d1 x_d2 y 0 900 z 44 47 qCN7i&k,  
       endcommand ulJ
YJ+CC!  
;在一定的范围内给巷道加锚杆，加帮锚杆和顶锚杆 e]h'  
loop k (250,500) =]"|x7'!  
        y_d=k*0.8-0.8 =
lQ[%&  
;帮锚杆的参数及循环，每帮加3个锚杆 5AU3s  
  loop a_2(1,3) ;(6lN<iU  
   z_d1=44.0 oW6<7>1M7  
   z_d1=z_d1+0.7*a_2 $t'I*k^N  
   ;左、右两帮锚杆的定义 |Eu~=J7@  
   command [zE
P|  
    sel cable beg x_d1,y_d,z_d1    end x_d11,y_d,z_d1 nseg 4 k
]pD3.QJ  
    sel cable beg x_d2,y_d,z_d1    end x_d21,y_d,z_d1 nseg 4 ;jI"|v{vnS  
   sel cable prop   emod 2e10   ytens 2.5e8   xcarea 1.0    gr_coh 1e10   gr_k 2e10   gr_per 1.0 'U$VOq?!  
   endcommand !Jl0Eu  
  end_loop e8<nPt`C  
;顶锚杆的参数及循环，顶锚杆加6根 [HDO^6U  
  loop a_3(1,6) ! -@!u   
   x_d3=x_d1 >+8I =S  
        x_d3=x_d3+0.7*a_3 r0 C6Ww7u  
     ;顶锚杆的定义 _\PoZ|G4y  
   command Q@VnJ,  
    sel cable beg x_d3,y_d,47  end x_d3,y_d,48.4 nseg 4 u6T?oK9j  
   sel cable prop   emod 2e10   ytens 2.5e8   xcarea 1.0    gr_coh 1e10   gr_k 2e10   gr_per 1.0 >irT|VTf  
   endcommand :/%xK"  
  end_loop !5!$h`g  
;结束每一个循环，大家可以根据位置来判断对应关系 rxeXz<  
end_loop [d>yo_iB  
end_loop RGI6W{\  
end F6VIH(  
aa_2 e/jM+%  
step 500 rd4'y~#S  
;开挖工作面、软化顶板 o{W]mr3D  
ini state=0 ,s&~U<Z  
model null range  x 300 520 y 300 600 z 44 47 ;采底煤 ODA#vAc!  
step 500 @ibPL+~-_  
model null range  x 300 520 y 300 600 z 47 50 ;放顶煤 wJ*
-K-  
prop den=2700  bulk=7.5e9 shear=30e8 fric=42 coh=1.3e6  range  x 300 520 y 300 600 z 50 70 [{LnE:  
step 1000 X.9MOdG70  
save 401.sav de{YgN  

这是个好东西啊。。。

谢谢分享

谢谢分享  学习了

煤层没有分步开采，直接一步到位？