不同支护应力场模拟
请教:FLAC如何才能有效实现不同支护强度应力场的区别和规律 newplo 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
以上是我写的一个简单的模型,但是改变预紧力以后也体现不出煤岩体周围应力场的区别,希望能得到高手的指点,不胜感谢! 怎么就没有人给点帮助啊。。。。。:'( 是所有人都不知道怎么实现。还有高人不愿意说啊? 又做了一些相关条件的模拟,区别还是不大,,,希望有懂的高人出来指点一二 支护应力场只存在大牛的论文里 我也正在学,求份程序,602498845@qq.com 好赞啊
页:
[1]