对锚杆、锚索支护不发挥作用的求助及探讨(附命令流)
本帖最后由 star486 于 2009-6-10 22:16 编辑我模拟的锚杆、锚索支护不发挥一点作用， 不知道为什么？？改了好几次参数都不行，而且地应力较大，岩性参数较小，但位移与无支护时无区别，不管如何调岩石参数都这样，望高手指点！
命令流在4#、5# 光说不行的，这样大家想帮忙都没有法！LZ还是传上命令流上来，大家看看哪里出问题嘛 好的，我传上来，呵呵 初始建模如下：
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 103
gen zone brick p0=(0,0,6) p1=(16,0,6) p2=(0,40,6) p3=(0,0,10) size 4 102
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 105
;模型PART2
gen zone brick p0=(16,0,0) p1=(38,0,0) p2=(16,40,0) p3=(16,0,6) size 11 103
gen zone brick p0=(16,0,6) p1=(38,0,6) p2=(16,40,6) p3=(16,0,10) size 44 408
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 103
gen zone brick p0=(38,0,6) p1=(54,0,6) p2=(38,40,6) p3=(38,0,10) size 4 102
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 105
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 沿空巷道开挖过程及加锚杆过程如下：
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 cableid=1 begin 34.2,y_d,14 end 33.2,y_d,15.73 nseg 5
sel cableid=1 begin 35.1,y_d,14 end 35.1,y_d,16 nseg 5
sel cableid=1 begin 36,y_d,14 end 36,y_d,16 nseg 5
sel cableid=1 begin 36.9,y_d,14 end 36.9,y_d,16 nseg 5
sel cableid=1 begin 37.8,y_d,14 end 38.8,y_d,15.73 nseg 5
sel cableid=2 begin 34,y_d,11.6 end 32.53,y_d,10.75 nseg 4
sel cableid=2 begin 34,y_d,12.5 end 32.3,y_d,12.5 nseg 4
sel cableid=2 begin 34,y_d,13.4 end 32.53,y_d,14.25 nseg 4
sel cableid=3 begin 38,y_d,11.6 end 39.47,y_d,10.75 nseg 4
sel cableid=3 begin 38,y_d,12.5 end 39.47,y_d,12.5 nseg 4
sel cableid=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 2e10ytens 154e3 xcarea 0.00031 gr_per 0.09106 &
gr_k 2e7 gr_coh 10e5 gr_fric 35 range id 1
sel cable prop emod 2e10ytens 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 35range id 4
;施加预应力
sel cablepretension 40e3 range id 1
sel cable id=2 pretension 30e3 range id 2 3
;锚索预应力
sel cablepretension 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
propshear 0.06e9 bulk 0.1e9 cohesion 0.8e6 fric 20 tension 0.6e6 density 2000range 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=5beg 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=6beg 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=7beg 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=8beg mbw_x y_0 mbw_z end mew_x y_0 mew_z nseg 12
end_command
end_loop
end
zhihu
;施加预应力
sel cableprop emod 2e10ytens 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 cableprop emod 2e10ytens 99e3 xcarea 0.00033 gr_per 0.09106 &
gr_k 2e7 gr_coh 10e5 gr_fric 35 rangeid 6 7
sel cablepretension 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 cablepretension 150e3 range id 8
solve
save 支护效果.sav
set log off 锚杆已经受力了，而且很多超过了预应力，问题量怎么位移及塑性区等与无支护时几本无差别啊，紧急求助，望各位大牛指正！ 你的预应力锚索肯定不起作用！！因为你没有模拟锚头！
详参：http://forum.simwe.com/viewthread.php?tid=758239&extra=page%3D1%26amp%3Bfilter%3Dtype%26amp%3Btypeid%3D1 但是加上锚杆后发现有力啊，怎么回事啊， 我是用全锚，应该不用模拟锚头吧 望高手指正！！ 能否将命令流打个包啊？ 都是模拟隧道，有没有基坑的？ 都是模拟隧道，有没有基坑的？
页:
[1]