flac3d数值模拟巷道锚杆模型

new

;断层左
gen zone brick size 30 50 55 &
p0 0 0 0 p1 50 0 0 p2 0 50 0 p3 0 0 55 p4 50 50 0 p5 0 50 55 p6 30 0 55 p7 30 50 55 rat 1 1 1 group 断层左

;断层右
gen zone brick size 50 50 55 &
p0 52 0 0 p1 90 0 0 p2 52 50 0 p3 32 0 55 p4 90 50 0 p5 32 50 55 p6 90 0 55 p7 90 50 55 rat 1 1 1 group 断层右
;查看计算模型
plot block group
group 灰岩 range z 0 8.5 group 断层左
group 底2 range z 8.5 23.5 group 断层左
group 灰岩 range z 23.5 25 group 断层左
group 底1 range z 25 38 group 断层左
group 煤 range z 38 44.5 group 断层左
group 顶1 range z 44.5 48 group 断层左
group 顶2 range z 48 55 group 断层左

group 灰岩 range z 0 5.5 group 断层右
group 底2 range z 5.5 20.5 group 断层右
group 灰岩 range z 20.5 22 group 断层右
group 底1 range z 22 35 group 断层右
group 煤 range z 35 41.5 group 断层右
group 顶1 range z 41.5 45 group 断层右
group 顶2 range z 45 55 group 断层右
;断层破碎带
gen zone brick size 4 50 110 &
p0 50 0 0 p1 52 0 0 p2 50 50 0 p3 30 0 55 p4 52 50 0 p5 30 50 55 p6 32 0 55 p7 32 50 55 rat 1 1 1 group 断层带

attach face
;查看计算模型
plot block group
save moxing.sav
restore moxing.sav
mo mohr
prop bulk 2.5e9 shear 1.87e9 cohesion 3.2e6 tension 1e6 friction 30 range group 灰岩
prop bulk 1.4e9 shear 0.87e9 cohesion 2.5e6 tension 0.3e6 friction 26 range group 底2
prop bulk 2.5e9 shear 1.87e9 cohesion 3.2e6 tension 1e6 friction 30 range group 灰岩
prop bulk 1.4e9 shear 0.87e9 cohesion 1.7e6 tension 0.3e6 friction 26 range group 底1
prop bulk 0.8e9 shear 0.56e9 cohesion 0.8e6 tension 0.1e6 friction 25 range group 煤
prop bulk 2.6e9 shear 2.5e9 cohesion 1.7e6 tension 0.3e6 friction 30 range group 顶1
prop bulk 2.6e9 shear 2.5e9 cohesion 3e6 tension 2e6 friction 30 range group 顶2
prop bulk 3.5e8 shear 3.2e8 cohesion 0.01e6 tension 0.02e6 friction 15 range group 断层带


ini dens 2200 range group 灰岩
ini dens 2060 range group 底2
ini dens 2200 range group 灰岩
ini dens 2060 range group 底1
ini dens 1220 range group 煤
ini dens 2200 range group 顶1
ini dens 2230 range group 顶2
ini dens 1000 range group 断层带



;应力梯度化
set grav 0 0 -10

ini szz -10.5e6 grad 0,0,2.5e4
ini syy -8.5e6 grad 0,0,2e4
ini sxx -8.5e6 grad 0,0,2e4

apply szz -7.5e6 range x 0 41 y 0 50 z 24.9 25.1
apply szz -7.5e6 range x 42 90 y 0 50 z 21.9 22.1

fix x range x -0.1 0.1
fix x range x 89.9 90.1
fix y range y -0.1 0.1
fix y range y 49.9 50.1
fix z range z -0.1 0.1
hist unbal
solve
save 初始应力.sav


;监测开挖过程中
hist id=2 gp zdis 20 25 38
hist id=2 gp zdis 30 25 38
hist id=2 gp zdis

40 25 38
hist id=2 gp zdis 50 25 38
hist id=2 gp zdis 60 25 38
hist id=2 gp zdis 70 25 38
hist id=2 gp zdis 80 25 38
;
;监测断层内
hist gp zdis 37 25 38.5
hist gp zdis 38.8 25 33.5
hist gp zdis 40.65 25 28.5

ini xdis=0 ydis=0 zdis=0
;定义一个函数 aa_2
def aa_2

x_d1=22

;定义几个变量，在下面的循环中要用到
;左帮锚杆的起、终点变量
x_d11=x_d1-0.9
;右帮锚杆的起、终点变量
x_d2=x_d1+6
x_d21=x_d2+0.9
;巷道的起、终点变量
y_d1=10
y_d2=y_d1+10

;开挖巷道

command
model null range x 0 90 y x_d1 x_d2 z 38.5 42.5
endcommand

;在一定的范围内给巷道加锚杆，加帮锚杆和顶锚杆
loop k (10,80)
y_d=k*0.8-0.8
;帮锚杆的参数及循环，每帮加3个锚杆
loop a_2(1,3)
z_d1=38.5
z_d1=z_d1+0.8*a_2
;左、右两帮锚杆的定义
command
sel cable beg y_d,x_d1,z_d1 end y_d,x_d11,z_d1 nseg 4
sel cable beg y_d,x_d2,z_d1 end y_d,x_d21,z_d1 nseg 4
sel cable prop emod 2e10 ytens 2.5e8 xcarea 1.0 gr_coh 1e10 gr_k 2e10 gr_per 1.0
endcommand
end_loop
;顶锚杆的参数及循环，顶锚杆加6根
loop a_3(1,6)
x_d3=x_d1
x_d3=x_d3+0.8*a_3
;顶锚杆的定义
command
sel cable beg y_d,x_d3,42.5 end y_d,x_d3,43.7 nseg 4
sel cable prop emod 2e10 ytens 2.5e8 xcarea 1.0 gr_coh 1e10 gr_k 2e10 gr_per 1.0
endcommand
end_loop
;结束每一个循环
end_loop
end
aa_2

set large

step 3000

save 巷道.sav