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

i;Cs,Esnf  
我的这个模拟是基于这样一个问题：采矿中开挖4条同样支护方式的巷道，但是巷道的位置不同（煤柱的宽度不同） {6{y"8  
;在其中使用了几重循环，可能对初学者有参考价值，我也常常看论坛的资料，现在自己贡献一点。 +13h*  
;这个文件可以直接运行的。 KiW4>@tY  
;建立网格 820m*900m*350m  网格51*30*30=4.59万网格 wI.i\S  
gen zone brick size 10 30  4 & e~R; 2bk  
p0   0 0 0  p1 295 0 0 p2 0 900 0   p3 0 0 40  ratio 0.8 1  1 \Zc$X^}vN  
gen zone brick size 8 30  4 & ?&;d#z*4  
p0 295 0 0  p1 335 0 0 p2 295 900 0 p3 295 0 40 Q|QVm,m  
gen zone brick size 15 30  4 & KilgeN
:  
p0 335 0 0  p1 485 0 0 p2 335 900 0 p3 335 0 40 9rXbv4{  
gen zone brick size 8 30  4 & CvfXm  
p0 485 0 0  p1 525 0 0 p2 485 900 0 p3 485 0 40 w}+#w8hu  
gen zone brick size 10 30  4 & QK _1!t3  
p0 525 0 0  p1 820 0 0 p2 525 900 0 p3 525 0 40 ratio 1.25 1 1 hoY.2 B_  
gen zone brick size 10 30  10 & %x@ D i`
;  
p0   0 0 40  p1 295 0 40 p2 0 900 40   p3 0 0 50  ratio 0.8 1  1 ah<1&UG,  
gen zone brick size 8 30  10 & >dKK [E/[d  
p0 295 0 40  p1 335 0 40 p2 295 900 40 p3 295 0 50 w1HE^ /  
gen zone brick size 15 30  10 & b~DtaGh  
p0 335 0 40  p1 485 0 40 p2 335 900 40 p3 335 0 50 rt">xVl  
gen zone brick size 8 30  10 & 9ZvBsG)  
p0 485 0 40  p1 525 0 40 p2 485 900 40 p3 485 0 50 7pMl:\  
gen zone brick size 10 30  10 & fm$eJu  
p0 525 0 40  p1 820 0 40 p2 525 900 40 p3 525 0 50 ratio 1.25 1 1 )
Y%>t  
gen zone brick size 10 30  6 & vzV,} S*c  
p0   0 0 50  p1 295 0 50 p2 0 900 50   p3 0 0 110  ratio 0.8 1  1 n,sf$9"  
gen zone brick size 8 30  6 & n][/c_]q  
p0 295 0 50  p1 335 0 50 p2 295 900 50 p3 295 0 110 K$OxeJP?F  
gen zone brick size 15 30  6 & (t&]u7Atr  
p0 335 0 50  p1 485 0 50 p2 335 900 50 p3 335 0 110 -c-af%xD  
gen zone brick size 8 30  6 & j.FA
!4L  
p0 485 0 50  p1 525 0 50 p2 485 900 50 p3 485 0 110 x
g=}MoX  
gen zone brick size 10 30  6 & (ii6w d<*  
p0 525 0 50  p1 820 0 50 p2 525 900 50 p3 525 0 110 ratio 1.25 1 1 :$k':0 n  
gen zone brick size 10 30  10 & x,$N!X  
p0   0 0 110  p1 295 0 110 p2 0 900 110   p3 0 0 350  ratio 0.8 1 1.25 .N2yn`
  
gen zone brick size 8 30  10 & )006\W|t9  
p0 295 0 110  p1 335 0 110 p2 295 900 110 p3 295 0 350  ratio 1 1 1.25 :#D?b.=  
gen zone brick size 15 30  10 &
 aEUC  
p0 335 0 110  p1 485 0 110 p2 335 900 110 p3 335 0 350  ratio 1 1 1.25 Vp8t8X1`  
gen zone brick size 8 30  10 & Fe 3*pUt  
p0 485 0 110  p1 525 0 110 p2 485 900 110 p3 485 0 350  ratio 1 1 1.25 V.zKjoky@  
gen zone brick size 10 30  10 & T=r-6eN  
p0 525 0 110  p1 820 0 110 p2 525 900 110 p3 525 0 350  ratio 1.25 1 1.25 @sQ^6FK0G  
r=GF*i[3  
;本构关系、参数 I)#=#eI*:  
m m q/y4HT,x  
prop den=2700  bulk=15e9  shear=60e8   fric=42  coh=2.61e6   range z 0  44
iEx
.BQ
+  
prop den=1500  bulk=1e9   shear=4e8    fric=30  coh=1.4e6    range z 44 50   fxfzi{}uj  
prop den=2700  bulk=15e9  shear=60e8   fric=42  coh=2.61e6   range z 50 350 &:}e`u@5|  
;边界、初始条件 r@C2zF7  
fix z range z  -0.1 0.1 L9tjHC]
 
fix x range x  -0.1 0.1 P^m+SAAB  
fix x range x  819.9 820.1 gXr"],OM;  
fix y range y  -0.1 0.1 +./c=o/v  
fix y range y  899.9 900.1 ' >(])Oq,  
;重力梯度化 XMhDx  
set grav 0 0 -10 HQHFD0hv  
ini szz -9.45e6  grad 0,0,2.7e4   range z 50 350 Ow7NOhw  
ini szz -8.85e6  grad 0,0,1.5e4   range z 44 50 1d/-SxhZ  
ini szz -9.378e6 grad 0,0,2.7e4   range z 0 44 RC
7|@a  
hist unbal K&FGTS,  
plot set plane nor 0 1 0 origin 0 100 0 i9Fg  
plot con szz   plan +5J"G/f  
step 3000 p@epl|IZp  
hist unbal 'J^ M`/  
save 4001.sav 50!/%  
;没开挖前的的自重平衡，为状态1 $r`K4g  
n]M1'yU  
rest 4001.sav pR_cI]{=SA  
;开挖巷道 \b{Aj,6,  
ini state=0 FTM(y CN  
ini xdis 0 ydis  0  zdis 0 6k6M&a  
ini xvel 0 yvel  0  zvel 0 yMdEH-?/  
;以下是自己写的几重循环的嵌套，可能对初学者有用。 D|-^}I4  
;开挖四条巷道加锚杆支护的实验 `$og]Dn;  
；定义一个函数 aa_2 x._IP,vRx^  
def aa_2 FR~YO|4?  
;定义一个循环变量 aaa,为了实现在不同位置开挖4条同样参数巷道，巷道位置循环，循环1 F) w.q  
loop aaa (1,4)  2|'v[  
  caseof aaa <p@c%e,_  
  case 1 a*LT<N  
   x_d1=295 XL[/)lX{  
  case 2 YnnpgR.  
   x_d1=520 y,qP$5xiq  
             ;10m煤柱 gcYx-gA}  
  case 3 fR_ jYP1  
   x_d1=310 m87,N~DP  
  case 4 _&S?uz m  
   x_d1=505 k=w;jX&;`  
  end_case ;>^oe:@  
;定义几个变量，在下面的循环中要用到 Bvzu{B%  
;左帮锚杆的起、终点变量  V|=PaO  
x_d11=x_d1-1.4 p- 5)J&  
;右帮锚杆的起、终点变量 B$~oZ'4v  
x_d2=x_d1+5 {\-rZb==F2  
x_d21=x_d2+1.4 .C^1.)  
;开挖巷道 !NWz  
       command &
`>[4D*  
   model null range  x x_d1 x_d2 y 0 900 z 44 47 B;9"=0  
       endcommand \s<7!NAE4  
;在一定的范围内给巷道加锚杆，加帮锚杆和顶锚杆 ,#3}TDC  
loop k (250,500) :}d`$2Dz  
        y_d=k*0.8-0.8 Ol,;BZHc\  
;帮锚杆的参数及循环，每帮加3个锚杆 [8*Ovd  
  loop a_2(1,3) |8I
#`  
   z_d1=44.0 cBf9-k  
   z_d1=z_d1+0.7*a_2 8r '  
   ;左、右两帮锚杆的定义 gfE<XrG  
   command .DSn H6O  
    sel cable beg x_d1,y_d,z_d1    end x_d11,y_d,z_d1 nseg 4 (;utiupW  
    sel cable beg x_d2,y_d,z_d1    end x_d21,y_d,z_d1 nseg 4 M99#\0=/  
   sel cable prop   emod 2e10   ytens 2.5e8   xcarea 1.0    gr_coh 1e10   gr_k 2e10   gr_per 1.0 d,=Kv  
   endcommand i`o}*`//  
  end_loop ""Ul6hRgv  
;顶锚杆的参数及循环，顶锚杆加6根 ?DcRD)X  
  loop a_3(1,6) n^|;J*rD  
   x_d3=x_d1 8 vNgePn  
        x_d3=x_d3+0.7*a_3 <hi@$.u_Q^  
     ;顶锚杆的定义 ~;
St,Fw<<  
   command {MmHR  
    sel cable beg x_d3,y_d,47  end x_d3,y_d,48.4 nseg 4 TR!^wB<F  
   sel cable prop   emod 2e10   ytens 2.5e8   xcarea 1.0    gr_coh 1e10   gr_k 2e10   gr_per 1.0 \)'nxFKqV  
   endcommand 1);$#Dlt k  
  end_loop `|K,E  
;结束每一个循环，大家可以根据位置来判断对应关系 5 e:Urv77  
end_loop b?Wg|D  
end_loop )6|7L)Dk  
end WqN=D5  
aa_2 `(A6uakd  
step 500 \m-fLX  
;开挖工作面、软化顶板 hoxn!x$?  
ini state=0 ~~:w^(s9  
model null range  x 300 520 y 300 600 z 44 47 ;采底煤 {zoUU  
step 500 [% chN/  
model null range  x 300 520 y 300 600 z 47 50 ;放顶煤 8r[ZGUV  
prop den=2700  bulk=7.5e9 shear=30e8 fric=42 coh=1.3e6  range  x 300 520 y 300 600 z 50 70 }Ictnb  
step 1000 2m.RM&TdB  
save 401.sav j]7|5mC78  

这是个好东西啊。。。

谢谢分享

谢谢分享  学习了

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