不同支护应力场模拟

xuyoulin

请教：FLAC如何才能有效实现不同支护强度应力场的区别和规律

xuyoulin

new
plo blo gro
gen zon brick p0 0 0 0 p1 100 0 0 p2 0 50 0 p3 0 0 6 size 50 25 3 group meiceng
gen zon brick p0 0 0 6 p1 100 0 6 p2 0 50 6 p3 0 0 11 size 50 25 5 group xishayandingban
gen zon brick p0 0 0 11 p1 100 0 11 p2 0 50 11 p3 0 0 50 size 50 25 20 group shayandingban
gen zon brick p0 0 0 -20 p1 100 0 -20 p2 0 50 -20 p3 0 0 0 size 50 25 10 group diban

;建模完成
   
     save moxing.sav
     rest moxing.sav

plo con szz

m m

;煤层
pro density 1550 bulk 1.6e9 she 0.9e9 fric 30 coh 2.7e6 ten 0.85e6 range group meiceng

;底板
pro density 2520 bulk 2e9 she 1.2e9 fric 33 coh 8e6 ten 4e6 range group diban

;顶板
pro density 2700 bulk 10e9 she 6e9 fric 30 coh 14e6 ten 5.2e6 range group xishayandingban
pro density 2700 bulk 10e9 she 6e9 fric 30 coh 14e6 ten 5.2e6 range group shayandingban
      ;setting initial conditions
;
        ini xdisp 0 ydisp 0 zdisp 0
        ini xvel 0 yvel 0 zvel 0
        ini state 0


;setting boundary conditions
      fix x range x -0.01 0.01
              fix x range x 99.99 100.01
                 fix y range y -0.01 0.01
                    fix y range y 49.99 50.01
                       fix z range z -20.01 -19.99

;设置重力参数
   set grav 0 0 -10

   apply szz -13.045e6 range z 49.99 50.01
;H=520
    ;assign initial stress state
        ini szz -13.045e6
            ini sxx -4.97e6
               ini syy -9.15e6
      solve
             save ini.sav
             rest ini.sav

plot con szz
model null range x 20 24 y 0 35 z 0 4
    def bolts
             k=0.5
                j=0.9
                   k2=35
                loop while k<k2
                                             
command
;
               sel cable id=1 begin 18 k 0.2 end 20 k 0.2 nseg=6     ;左侧5根
               sel cable id=1 begin 18 k 1.1 end 20 k 1.1 nseg=6
                        sel cable id=1 begin 18 k 2 end 20 k 2 nseg=6
                        sel cable id=1 begin 18 k 2.9 end 20 k 2.9 nseg=6
                        sel cable id=1 begin 18 k 3.8 end 20 k 3.8 nseg=6

                        sel cable id=1 begin 20.2 k 4 end 20.2 k 6.2 nseg=6   ;顶板5根
               sel cable id=1 begin 21.1 k 4 end 21.1 k 6.2 nseg=6
               sel cable id=1 begin 22 k 4 end 22 k 6.2 nseg=6
               sel cable id=1 begin 22.9 k 4 end 22.9 k 6.2 nseg=6
               sel cable id=1 begin 23.8 k 4 end 23.8 k 6.2 nseg=6
               
               sel cable id=1 begin 24 k 0.2 end 26 k 0.2 nseg=6     ;右侧5根
               sel cable id=1 begin 24 k 1.1 end 26 k 1.1 nseg=6
               sel cable id=1 begin 24 k 2 end 26 k 2 nseg=6
               sel cable id=1 begin 24 k 2.9 end 26 k 2.9 nseg=6
               sel cable id=1 begin 24 k 3.8 end 26 k 3.8 nseg=6

               sel cable prop emod 1e11 ytens 8.39e4 xcarea 2.54e-4 &
                              gr_coh 1.3e5 gr_k 7e9 gr_per 1.0 range id=1
               sel cable pretension 50e4 range id=1
               end_command
               k=k+j
               end_loop
                                                   
               end
               bolts
solve
   save 掘进50KN修改.sav

以上是我写的一个简单的模型，但是改变预紧力以后也体现不出煤岩体周围应力场的区别，希望能得到高手的指点，不胜感谢！

xuyoulin

怎么就没有人给点帮助啊。。。。。

xuyoulin

是所有人都不知道怎么实现。还有高人不愿意说啊？

xuyoulin

又做了一些相关条件的模拟，区别还是不大，，，希望有懂的高人出来指点一二

wangjun1988

支护应力场只存在大牛的论文里

jackjia1988

我也正在学，求份程序，602498845@qq.com

mahoubao

好  赞啊