给大家分享一个三软巷道支护的例子

kywenli

建模就不在这儿放了了，就给个计算的过程

kywenli

；平衡
def para
;基岩
p1=0.25
E1=2.1e7
bulk1=E1/(3*(1-2*p1))
shear1=E1/(2*(1+p1))
coh1=1.99e6
fric1=30
tens1=1.99e6
dens1=2300
;泥质砂岩
p2=0.33
E2=8.55e6
bulk2=E2/(3*(1-2*p2))
shear2=E2/(2*(1+p2))
coh2=1.4e6
fric2=27
tens2=1.99e6
dens2=2300
;粘土岩
p3=0.33
E3=5.4e6
bulk3=E3/(3*(1-2*p3))
shear3=E3/(2*(1+p3))
coh3=8.0e4
fric3=25.5
tens3=1.50e6
dens3=2300
;煤层
p4=0.25
E4=1.65e6
bulk4=E4/(3*(1-2*p4))
shear4=E4/(2*(1+p4))
coh4=1.0e5
fric4=17
tens4=0.12e6
dens4=1310

;砂质泥岩
p5=0.26
E5=3.0e6
bulk5=E5/(3*(1-2*p5))
shear5=E5/(2*(1+p5))
coh5=2.73e5
fric5=25
tens5=1.2e6
dens5=2305
;泥岩
p6=0.26
E6=1.0e7
bulk6=E6/(3*(1-2*p6))
shear6=E6/(2*(1+p6))
coh6=2.73e6
fric6=25
tens6=1.2e6
dens6=2305
;细砂岩
p7=0.26
E7=2.1e7
bulk7=E7/(3*(1-2*p7))
shear7=E7/(2*(1+p7))
coh7=2.8e7
fric7=29
tens7=1.2e6
dens7=2305
;覆岩
p8= 0.27
E8=2.5e7
bulk8=E8/(3*(1-2*p8))
shear8=E8/(2*(1+p8))
coh8=2.552e6
fric8=44
tens8=1.58e6
dens8=2531
end
para
model mohr
prop bulk=bulk1 shear=shear1 fric=fric1 coh=coh1 tension=tens1 density=dens1 range group 1
prop bulk=bulk2 shear=shear2 fric=fric2 coh=coh2 tension=tens2 density=dens2 range group 2
prop bulk=bulk3 shear=shear3 fric=fric3 coh=coh3 tension=tens3 density=dens3 range group 3
prop bulk=bulk4 shear=shear4 fric=fric4 coh=coh4 tension=tens4 density=dens4 range group 4
prop bulk=bulk4 shear=shear4 fric=fric4 coh=coh4 tension=tens4 density=dens4 range group 5
prop bulk=bulk5 shear=shear5 fric=fric5 coh=coh5 tension=tens5 density=dens5 range group 6
prop bulk=bulk6 shear=shear6 fric=fric6 coh=coh6 tension=tens6 density=dens6 range group 7
prop bulk=bulk7 shear=shear7 fric=fric7 coh=coh7 tension =tens7 density=dens7 range group 8
prop bulk=bulk8 shear=shear8 fric=fric8 coh=coh8 tension=tens8 density=dens8 range group 9
fix x range x -.1 .1
fix x range x 69.9 70.1
fix y range y -.1 .1
fix y range y 19.9 20.1
fix z range z -.1 .1
set grav 0 0 -10
apply szz -5.25e6 range z 68.69 68.89
pl bl st
solve

kywenli

；开挖
mo null range group 5
set large
step 25

kywenli

；加锚杆
ini xdis=0 ydis=0 zdis=0
hist unbal
;上运
;zuosangen2.3
def place_cables
loop n (1,40)
y_d=float(n)-0.7
command
sel cable id 1 beg 30.00,y_d,24.82   end 28.80,y_d,24.82    nseg 6
sel cable id 1 beg 30.00,y_d,25.72   end 28.80,y_d,25.72    nseg 6
sel cable id 1 beg 30.00,y_d,26.62   end 28.80,y_d,26.62    nseg 6
;yousangen
sel cable id 1 beg 34.00,y_d,26.27   end 35.20,y_d,26.27    nseg 6
sel cable id 1 beg 34.00,y_d,27.17   end 35.20,y_d,27.17    nseg 6
sel cable id 1 beg 34.00,y_d,28.07   end 35.20,y_d,28.07    nseg 6

;shangjiugen
sel cable id 1 beg 33.00,y_d,28.60   end 32.25,y_d,31.40    nseg 6
sel cable id 1 beg 32.00,y_d,28.24   end 31.25,y_d,30.67    nseg 6  
sel cable id 1 beg 31.00,y_d,27.88   end 30.25,y_d,29.95    nseg 6

sel cable id=1 pretension 60e3
sel cable id=1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 &
              gr_coh 10e5 gr_k 2e7
endcommand
endloop
end
place_cables
plot add cont disp out on
plot add sel cable force

solve

liuxueqiang

上传个txt多好

qinggu20082008

下回去好好看看的，一定有收获的

Visual1988

犹抱琵琶半遮面

bgt1208

好人。。。。。。。。

cuhammer

能不能把建模的程序也发上来  谢谢  好好学习一下

江beyond

要遇到不懂得还望楼主能多多赐教！

icanlhs

谢谢楼主！

xqwzhym3

模型也发上来啊  flac建模也是必修的一课呢

superlxystar

怎么没有加初始地应力，这样的结果合理吗

liqiangre

不错,!学习!...................

youpeihai632

好东西~

侯健

挺好的，现下载看看

phymos

全长锚固锚杆

huawen226

谢谢楼主！

bobyzhang

最好把建模也放出来

syi450643708

有没有讲述具体点的啊！！