对锚杆、锚索支护不发挥作用的求助及探讨(附命令流)

star486

我模拟的锚杆、锚索支护不发挥一点作用， 不知道为什么？？改了好几次参数都不行，而且地应力较大，岩性参数较小，但位移与无支护时无区别，不管如何调岩石参数都这样，望高手指点！
命令流在4#、5#

jihezhilv

光说不行的，这样大家想帮忙都没有法！LZ还是传上命令流上来，大家看看哪里出问题嘛

star486

好的，我传上来，呵呵

star486

初始建模如下：

new
plot block state
plot show


;模型PART1

gen zone brick p0=(0,0,0) p1=(16,0,0) p2=(0,40,0) p3=(0,0,6) size 4 10  3

gen zone brick p0=(0,0,6) p1=(16,0,6) p2=(0,40,6) p3=(0,0,10) size 4 10  2

gen zone brick p0=(0,0,10) p1=(16,0,10) p2=(0,40,10) p3=(0,0,11) size 4 20 2

gen zone brick p0=(0,0,11) p1=(16,0,11) p2=(0,40,11) p3=(0,0,14) size 4 40 6

gen zone brick p0=(0,0,14) p1=(16,0,14) p2=(0,40,14) p3=(0,0,16) size 4 20 2
gen zone brick p0=(0,0,16) p1=(16,0,16) p2=(0,40,16) p3=(0,0,20) size 4 10 2

gen zone brick p0=(0,0,20) p1=(16,0,20) p2=(0,40,20) p3=(0,0,30) size 4 10  5

;模型PART2

gen zone brick p0=(16,0,0) p1=(38,0,0) p2=(16,40,0) p3=(16,0,6) size 11 10  3

gen zone brick p0=(16,0,6) p1=(38,0,6) p2=(16,40,6) p3=(16,0,10) size 44 40  8

gen zone brick p0=(16,0,10) p1=(38,0,10) p2=(16,40,10) p3=(16,0,11) size 44 80 2

gen zone brick p0=(16,0,11) p1=(38,0,11) p2=(16,40,11) p3=(16,0,14) size 44 80 6

gen zone brick p0=(16,0,14) p1=(38,0,14) p2=(16,40,14) p3=(16,0,16) size 44 80 4

gen zone brick p0=(16,0,16) p1=(38,0,16) p2=(16,40,16) p3=(16,0,20) size 44 40 8

gen zone brick p0=(16,0,20) p1=(38,0,20) p2=(16,40,20) p3=(16,0,30) size 11 10 5

;模型PART3


gen zone brick p0=(38,0,0) p1=(54,0,0) p2=(38,40,0) p3=(38,0,6) size 4 10  3

gen zone brick p0=(38,0,6) p1=(54,0,6) p2=(38,40,6) p3=(38,0,10) size 4 10  2

gen zone brick p0=(38,0,10) p1=(54,0,10) p2=(38,40,10) p3=(38,0,11) size 4 20 2

gen zone brick p0=(38,0,11) p1=(54,0,11) p2=(38,40,11) p3=(38,0,14) size 4 40 6

gen zone brick p0=(38,0,14) p1=(54,0,14) p2=(38,40,14) p3=(38,0,16) size 4 20 2

gen zone brick p0=(38,0,16) p1=(54,0,16) p2=(38,40,16) p3=(38,0,20) size 4 10 2

gen zone brick p0=(38,0,20) p1=(54,0,20) p2=(38,40,20) p3=(38,0,30) size 4 10  5


save 初始网格.sav
restore 初始网格.sav

attach face

model mohr

;赋材料性质
macro 砂泥岩互层 'shear 0.5e9 bulk 0.8e9 cohesion 5e6 fric 40 tension 4e6 density 2500'
macro 砂质泥岩   'shear 0.25e9 bulk 0.4e9 cohesion 3.2e6 fric 35 tension 2e6 density 2500'
macro 煤         'shear 0.18e9 bulk 0.3e9 cohesion 2.5e6 fric 30 tension 1.6e6 density 1350'
macro 粉砂岩     'shear 0.3e9 bulk 0.4e9 cohesion 3.5e6 fric 36 tension 2e6 density 2600'
macro 细砂岩     'shear 0.5e9 bulk 1.0e9 cohesion 6e6 fric 43 tension 5e6 density 2700'


prop 砂泥岩互层 range z 0 10
prop 砂质泥岩   range z 10 11
prop 煤         range z 11 14
prop 粉砂岩     range z 14 16
prop 细砂岩     range z 16 30
   

fix x range x -0.1 0.1
fix x range x 53.9 54.1
fix y range y -.1 .1
fix y range y 39.9 40.1
fix z range z -.1 .1

apply szz -16.25e6 range z 29.9 30.1
set grav 0 0 -9.8
initial szz -17e6 grad 0 0 2.5e4
initial syy -10e6
initial sxx -10e6

solve
save 初始平衡.sav

star486

沿空巷道开挖过程及加锚杆过程如下：
set log on
set logfile zhihu.txt
restore 初始平衡.sav
plot show
plot block state



;上区段巷道开挖
model null range x 34 38 y 0 40 z 11 14
step 50

; 安装锚杆

def place_cable1
   loop n (1,45)
     y_d=n * 0.9 - 0.9
     command

      sel cable  id=1 begin 34.2,y_d,14   end 33.2,y_d,15.73 nseg 5
      sel cable  id=1 begin 35.1,y_d,14   end 35.1,y_d,16    nseg 5
      sel cable  id=1 begin 36,  y_d,14     end 36,y_d,16    nseg 5
      sel cable  id=1 begin 36.9,y_d,14   end 36.9,y_d,16    nseg 5
      sel cable  id=1 begin 37.8,y_d,14   end 38.8,y_d,15.73 nseg 5


      sel cable  id=2 begin 34,y_d,11.6 end 32.53,y_d,10.75 nseg 4
      sel cable  id=2 begin 34,y_d,12.5 end 32.3,y_d,12.5 nseg 4
      sel cable  id=2 begin 34,y_d,13.4 end 32.53,y_d,14.25 nseg 4


      sel cable  id=3 begin 38,y_d,11.6 end 39.47,y_d,10.75 nseg 4
      sel cable  id=3 begin 38,y_d,12.5 end 39.47,y_d,12.5 nseg 4
      sel cable  id=3 begin 38,y_d,13.4 end 39.47,y_d,14.25 nseg 4


sel beam id=1 begin 34.2 y_d 14 end 37.8 y_d 14 nseg=4
sel beam id=2 begin 34 y_d 11.6 end 34 y_d 13.4 nseg=4
sel beam id=3 begin 38 y_d 11.6 end 38 y_d 13.4 nseg=4

sel beam prop emod 2e11 nu 0.15 xcarea=2.04e-3 xcj 0 xciy 72e-6 xciz 72e-6 range id 1 3
      endcommand
   endloop
end




def place_cable2
   loop n (1,15)
     a_d=n * 2.7-2.7+0.45
     command

;          顶板
      sel cable id=4 begin 36,a_d,14   end 36,a_d,20 nseg 12
      
     endcommand
   endloop
end
place_cable1
place_cable2

sel cable   prop emod 2e10  ytens 154e3 xcarea 0.00031 gr_per 0.09106 &
                         gr_k 2e7 gr_coh 10e5 gr_fric 35 range id 1
sel cable   prop emod 2e10  ytens 99e3 xcarea 0.00031 gr_per 0.09106 &
                         gr_k 2e7 gr_coh 10e5 gr_fric 35 range id 2 3

                     
sel cable prop emod 2e12 ytension 300e3 xcarea 0.0002  &
                        gr_per 0.09106 gr_k 2e7 gr_coh 15e5 gr_fric 35  range id 4
;施加预应力

        sel cable  pretension 40e3 range id 1

         sel cable id=2 pretension 30e3 range id 2 3




;锚索预应力
        
        sel cable  pretension 150e3 range id 4


step 1000
save 前一面巷道挖.sav


;上区段工作面开采
model null range x 38 54 y 0 40 z 11 14
step 200


model mohr range x 34 54 y 0 40 z 11 14
prop  shear 0.06e9 bulk 0.1e9 cohesion 0.8e6 fric 20 tension 0.6e6 density 2000  range x 34 54 y 0 40 z 11 14
solve
save 条件准备.sav
ini xdisp=0
ini ydisp=0
ini zdisp=0
;监测

hist id=1 gp zdisp 27 4 14.2

hist id=2 gp zdisp 27 4 10.8

hist id=3 gp xdisp 24.8 4 12.5

hist id=4 gp xdisp 28.8 4 12.5



hist id=5 gp zdisp 27 19 14.2

hist id=6 gp zdisp 27 19 10.8

hist id=7 gp xdisp 24.8 19 12.5

hist id=8 gp xdisp 28.8 19 12.5


def dingbanshuzhu
   array bw(2,5) ew(2,5)

   bw(1,1)=25.2
   bw(1,2)=26.1
   bw(1,3)=27
   bw(1,4)=27.9
   bw(1,5)=28.8
  

   bw(2,1)=14
   bw(2,2)=14
   bw(2,3)=14
   bw(2,4)=14
   bw(2,5)=14

   ew(1,1)=24.2
   ew(1,2)=26.1
   ew(1,3)=27
   ew(1,4)=27.9
   ew(1,5)=29.8

   ew(2,1)=15.7
   ew(2,2)=16
   ew(2,3)=16
   ew(2,4)=16
   ew(2,5)=15.7
end

def yankongshuzhu
   array ybw(2,3) yew(2,3)

   ybw(1,1)=29
   ybw(1,2)=29
   ybw(1,3)=29
   

   ybw(2,1)=11.6
   ybw(2,2)=12.5
   ybw(2,3)=13.4
  

   yew(1,1)=30.47
   yew(1,2)=30.7
   yew(1,3)=30.47
   


   yew(2,1)=10.75
   yew(2,2)=12.5
   yew(2,3)=14.25
   
end

def shitishuzhu
   array sbw(2,3) sew(2,3)
   ;x起
   sbw(1,1)=25
   sbw(1,2)=25
   sbw(1,3)=25
   

   sbw(2,1)=11.6
   sbw(2,2)=12.5
   sbw(2,3)=13.4
  

   sew(1,1)=23.5
   sew(1,2)=23.3
   sew(1,3)=23.5
   


   sew(2,1)=10.75
   sew(2,2)=12.5
   sew(2,3)=14.25
   
end
def maosuoshuzhu
   array mbw(2,1) mew(2,1)

   mbw(1,1)=27
   

   mbw(2,1)=14
   

   mew(1,1)=27
   

   mew(2,1)=20

end
   dingbanshuzhu
   yankongshuzhu
   shitishuzhu
   maosuoshuzhu


def dingbanzuobiao
   
   bw_x=bw(1,index)
   bw_z=bw(2,index)
   ew_x=ew(1,index)
   ew_z=ew(2,index)
   
end
def shitizuobiao
   sbw_x=sbw(1,mm)
   sbw_z=sbw(2,mm)
   sew_x=sew(1,mm)
   sew_z=sew(2,mm)

end
def yankongzuobiao
   ybw_x=ybw(1,ii)
   ybw_z=ybw(2,ii)
   yew_x=yew(1,ii)
   yew_z=yew(2,ii)
end
def maosuozuobiao
   mbw_x=mbw(1,1)
   mbw_z=mbw(2,1)
   mew_x=mew(1,1)
   mew_z=mew(2,1)
end


model null range x 25 29 y 0 40 z 11 14
step 50
;zhihu
def zhihu

   loop n(1,45)
      y_d=n*0.9-0.9
   

       loop index(1,5)
          dingbanzuobiao
          command
               sel cable id=5  beg bw_x y_d bw_z end ew_x y_d ew_z nseg 5
          endcommand
        end_loop


        loop mm(1,3)
        shitizuobiao
         command
             sel cable id=6  beg sbw_x y_d sbw_z end sew_x y_d sew_z nseg 4
            
          endcommand
         endloop


         loop ii(1,3)
        yankongzuobiao
              command
               sel cable id=7  beg ybw_x y_d ybw_z end yew_x y_d yew_z nseg 4

              end_command
        end_loop   
   ;钢带
  command
     sel beam id=5 begin 25.2 y_d 14 end 28.8 y_d 14 nseg=4
     sel beam id=6 begin 25 y_d 11.6 end 25 y_d 13.4 nseg=4
     sel beam id=7 begin 29 y_d 11.6 end 29 y_d 13.4 nseg=4

     sel beam prop emod 2e11 nu 0.15 xcarea=2.04e-3 xcj 0 xciy 72e-6 xciz 72e-6 range id 5 7
  endcommand
end_loop
  
  ;加锚索
  loop n(1,15)
   y_0=0.45+2.7*n-2.7
     maosuozuobiao
          command
             sel cable id=8  beg mbw_x y_0 mbw_z end mew_x y_0 mew_z nseg 12
          end_command
   end_loop
  
end
zhihu
;施加预应力



        sel cable  prop emod 2e10  ytens 154e3 xcarea 0.00031 gr_per 0.09106 &
                         gr_k 2e7 gr_coh 10e5 gr_fric 35 range id=5

        sel cable pretension 40e3 range id=5







         sel cable  prop emod 2e10  ytens 99e3 xcarea 0.00033 gr_per 0.09106 &
                         gr_k 2e7 gr_coh 10e5 gr_fric 35 range  id 6 7
        
        sel cable  pretension 30e3 range id 6 7


   
      
;锚索预应力


         sel cable prop emod 2e12 ytension 300e3 xcarea 0.0002  &
                        gr_per 0.09106 gr_k 2e7 gr_coh 15e5 gr_fric 35 range id 8

        
        sel cable  pretension 150e3 range id 8


solve
save 支护效果.sav

set log off

star486

锚杆已经受力了，而且很多超过了预应力，问题量怎么位移及塑性区等与无支护时几本无差别啊，紧急求助，望各位大牛指正！

jihezhilv

你的预应力锚索肯定不起作用！！因为你没有模拟锚头！
详参：http://forum.simwe.com/viewthrea ... %26amp%3Btypeid%3D1

star486

但是加上锚杆后发现有力啊，怎么回事啊， 我是用全锚，应该不用模拟锚头吧

star486

望高手指正！！

ether2008

能否将命令流打个包啊？

888438

都是模拟隧道，有没有基坑的？

li_2001

都是模拟隧道，有没有基坑的？