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

R8tF/dx>7  
我的这个模拟是基于这样一个问题：采矿中开挖4条同样支护方式的巷道，但是巷道的位置不同（煤柱的宽度不同） tVunh3-  
;在其中使用了几重循环，可能对初学者有参考价值，我也常常看论坛的资料，现在自己贡献一点。 S."7+g7Ar  
;这个文件可以直接运行的。 Kh)FyV  
;建立网格 820m*900m*350m  网格51*30*30=4.59万网格 s|Zx(.EP  
gen zone brick size 10 30  4 & .!i0_Rv5x  
p0   0 0 0  p1 295 0 0 p2 0 900 0   p3 0 0 40  ratio 0.8 1  1 yLW/ -%I#u  
gen zone brick size 8 30  4 & :lK4 db  
p0 295 0 0  p1 335 0 0 p2 295 900 0 p3 295 0 40 " \:ced  
gen zone brick size 15 30  4 & a"QU:<-v  
p0 335 0 0  p1 485 0 0 p2 335 900 0 p3 335 0 40 \OcMiuw  
gen zone brick size 8 30  4 & 84'?um  
p0 485 0 0  p1 525 0 0 p2 485 900 0 p3 485 0 40 R uLvG+  
gen zone brick size 10 30  4 & NY@"&p'Q  
p0 525 0 0  p1 820 0 0 p2 525 900 0 p3 525 0 40 ratio 1.25 1 1 >!}`%pk(  
gen zone brick size 10 30  10 & kg3ppt  
p0   0 0 40  p1 295 0 40 p2 0 900 40   p3 0 0 50  ratio 0.8 1  1 rs:Q%V ^  
gen zone brick size 8 30  10 & A:eG5K}
 
p0 295 0 40  p1 335 0 40 p2 295 900 40 p3 295 0 50 RlsVC_H\  
gen zone brick size 15 30  10 & Zm(dY*z5:J  
p0 335 0 40  p1 485 0 40 p2 335 900 40 p3 335 0 50 A W)a">|  
gen zone brick size 8 30  10 & '&by3y5w-3  
p0 485 0 40  p1 525 0 40 p2 485 900 40 p3 485 0 50 pCC7(Ouo  
gen zone brick size 10 30  10 & Pd~MiyO;K  
p0 525 0 40  p1 820 0 40 p2 525 900 40 p3 525 0 50 ratio 1.25 1 1 J{Tq%\a3  
gen zone brick size 10 30  6 & zE`R,:VI  
p0   0 0 50  p1 295 0 50 p2 0 900 50   p3 0 0 110  ratio 0.8 1  1 $}R$t-  
gen zone brick size 8 30  6 & Pr+~Kif  
p0 295 0 50  p1 335 0 50 p2 295 900 50 p3 295 0 110 B wC+o
v=  
gen zone brick size 15 30  6 & ]TSg!H  
p0 335 0 50  p1 485 0 50 p2 335 900 50 p3 335 0 110 &aPl`"j  
gen zone brick size 8 30  6 & -(>qu.[8=  
p0 485 0 50  p1 525 0 50 p2 485 900 50 p3 485 0 110 kX[fy7rVt  
gen zone brick size 10 30  6 & aV>aiR=  
p0 525 0 50  p1 820 0 50 p2 525 900 50 p3 525 0 110 ratio 1.25 1 1 t};~H\:  
gen zone brick size 10 30  10 & aqMc6N`z  
p0   0 0 110  p1 295 0 110 p2 0 900 110   p3 0 0 350  ratio 0.8 1 1.25 \dq!q=b\  
gen zone brick size 8 30  10 & %y"J8;U  
p0 295 0 110  p1 335 0 110 p2 295 900 110 p3 295 0 350  ratio 1 1 1.25 H{V)g  
gen zone brick size 15 30  10 & XoR>H4xh  
p0 335 0 110  p1 485 0 110 p2 335 900 110 p3 335 0 350  ratio 1 1 1.25 q* !3C  
gen zone brick size 8 30  10 & Ue$zH"w  
p0 485 0 110  p1 525 0 110 p2 485 900 110 p3 485 0 350  ratio 1 1 1.25 )ofm_R'q*  
gen zone brick size 10 30  10 & 8$vK5Dnn8  
p0 525 0 110  p1 820 0 110 p2 525 900 110 p3 525 0 350  ratio 1.25 1 1.25 5]"SGP  
8g*hvPc  
;本构关系、参数 0mexF@  
m m }nud  
prop den=2700  bulk=15e9  shear=60e8   fric=42  coh=2.61e6   range z 0  44 9<s4yZF@x  
prop den=1500  bulk=1e9   shear=4e8    fric=30  coh=1.4e6    range z 44 50   a3@w|KLt  
prop den=2700  bulk=15e9  shear=60e8   fric=42  coh=2.61e6   range z 50 350 0oU=RbC  
;边界、初始条件 ksWSMxm  
fix z range z  -0.1 0.1 x*}bo))hb  
fix x range x  -0.1 0.1 /\=MBUN  
fix x range x  819.9 820.1 zgKY4R{V  
fix y range y  -0.1 0.1 +.Bmkim  
fix y range y  899.9 900.1 R1ktj  
;重力梯度化 V U~Dk);Bv  
set grav 0 0 -10 & ,L9OU  
ini szz -9.45e6  grad 0,0,2.7e4   range z 50 350 ~`eHHgX  
ini szz -8.85e6  grad 0,0,1.5e4   range z 44 50 vR>o}%`  
ini szz -9.378e6 grad 0,0,2.7e4   range z 0 44 $-vo}k%M  
hist unbal *P2[qhP2  
plot set plane nor 0 1 0 origin 0 100 0 XU-m"_t  
plot con szz   plan M>l^%`  
step 3000 &L4 q10-N
 
hist unbal i0pU!`0  
save 4001.sav wW`}VKu  
;没开挖前的的自重平衡，为状态1 Vhg1/EgUr  
A^Cj1:,  
rest 4001.sav zAScRg$:?  
;开挖巷道 SEF6B45}1  
ini state=0 YQ(Po!NI\'  
ini xdis 0 ydis  0  zdis 0 $.K?N@(W  
ini xvel 0 yvel  0  zvel 0 GsV4ZZ  
;以下是自己写的几重循环的嵌套，可能对初学者有用。 Z1,gtl ?  
;开挖四条巷道加锚杆支护的实验 P $>`  
；定义一个函数 aa_2 HJhH-\{@  
def aa_2 c-CYdi@  
;定义一个循环变量 aaa,为了实现在不同位置开挖4条同样参数巷道，巷道位置循环，循环1 sR_xe}-  
loop aaa (1,4) Imym+  
  caseof aaa w+AuMc  
  case 1 B0)]s<<  
   x_d1=295 OXhAha`R  
  case 2 >+9JD%]x]  
   x_d1=520 8
i0  
             ;10m煤柱 h##WA=1QZ  
  case 3 [}o~PN:sT(  
   x_d1=310 F 3RB  
  case 4 0[a}n6XTk  
   x_d1=505 (ku5WWJ  
  end_case (#Mp 5C'X  
;定义几个变量，在下面的循环中要用到 (> "QVxr  
;左帮锚杆的起、终点变量 ]~c+'E`  
x_d11=x_d1-1.4 D`o<,Y  
;右帮锚杆的起、终点变量 V1bh|+o9  
x_d2=x_d1+5 F@vbSFv)/  
x_d21=x_d2+1.4 sQ1jrkm  
;开挖巷道 f_6`tq m%  
       command r
%=[
},JQ  
   model null range  x x_d1 x_d2 y 0 900 z 44 47 L-h$Z
0]_F  
       endcommand jBtj+TL8  
;在一定的范围内给巷道加锚杆，加帮锚杆和顶锚杆 2WCLS{@'  
loop k (250,500) .p@N:)W6  
        y_d=k*0.8-0.8 \e:d)^cbh  
;帮锚杆的参数及循环，每帮加3个锚杆 LtKB v4  
  loop a_2(1,3) m *X7T  
   z_d1=44.0  '3,\@4  
   z_d1=z_d1+0.7*a_2 T] | d5E  
   ;左、右两帮锚杆的定义 y{=NP  
   command `~F5wh~  
    sel cable beg x_d1,y_d,z_d1    end x_d11,y_d,z_d1 nseg 4 c{3rl;Cs  
    sel cable beg x_d2,y_d,z_d1    end x_d21,y_d,z_d1 nseg 4 EESN\_{~.  
   sel cable prop   emod 2e10   ytens 2.5e8   xcarea 1.0    gr_coh 1e10   gr_k 2e10   gr_per 1.0 ^PqF<d6  
   endcommand gA^q^>7  
  end_loop 9y]J/1#  
;顶锚杆的参数及循环，顶锚杆加6根 (.X]F_*sc  
  loop a_3(1,6) FxU'LN<;HY  
   x_d3=x_d1 Wk7WK` >i  
        x_d3=x_d3+0.7*a_3 g1[&c+=U`P  
     ;顶锚杆的定义 'ZHdV,dd  
   command < K%j  
    sel cable beg x_d3,y_d,47  end x_d3,y_d,48.4 nseg 4 80hme+e  
   sel cable prop   emod 2e10   ytens 2.5e8   xcarea 1.0    gr_coh 1e10   gr_k 2e10   gr_per 1.0 trYTs,KV  
   endcommand M<`|CVl  
  end_loop PpOlt.yui  
;结束每一个循环，大家可以根据位置来判断对应关系 t&RruwN_;  
end_loop /9hR  
end_loop rY}B-6qJn  
end 1mW%  
aa_2 S*t%RZ~a  
step 500 D<):ZfUbI  
;开挖工作面、软化顶板 uHM@h{r  
ini state=0 $k0H9_  
model null range  x 300 520 y 300 600 z 44 47 ;采底煤 <Oz66bTze  
step 500 RUXCq`)"<  
model null range  x 300 520 y 300 600 z 47 50 ;放顶煤 Dkw7]9Qm  
prop den=2700  bulk=7.5e9 shear=30e8 fric=42 coh=1.3e6  range  x 300 520 y 300 600 z 50 70 hnxc`VX>g  
step 1000 "yh Pm  
save 401.sav 0i@:KYP  

这是个好东西啊。。。

谢谢分享

谢谢分享  学习了

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