论坛 › G0201：FLAC3D/FLAC2D › 锚杆支护模拟后，如何进行后处理

hwaust 发表于 2009-10-5 15:47:07

锚杆支护模拟后，如何进行后处理

我做了一个锚杆支护的东西，想调出锚杆受力和变形情况，请问哪位兄弟帮我一下，给出调出锚杆应力和应变的命令流？还有衬砌结构的应力和变形怎么调？？？？？

lookcity 发表于 2009-10-5 16:53:07

如图

hwaust 发表于 2009-10-7 09:41:30

谢谢楼上兄弟

55709673 发表于 2009-10-7 13:49:00

2楼的兄弟，你这是哪个版本啊？

lookcity 发表于 2009-10-7 14:10:11

2楼的兄弟，你这是哪个版本啊？
55709673 发表于 2009-10-7 13:49 http://forum.simwe.com/images/common/back.gif
3.00-261

hwaust 发表于 2009-10-7 15:16:39

二楼兄弟，我按照您讲的方法怎么不能显示衬砌结构的内力和位移图呀，有没有命令流显示呀？？

lookcity 发表于 2009-10-7 16:05:49

显示锚杆轴力：plot add sel cable force

衬砌：
plot add sel l coup str n
plot add sel l coup str s

hwaust 发表于 2009-10-14 14:10:20

7# lookcity
你好，帮我看看有没有问题？我求解后，发现位移图不连续，还有怎样在程序中添加喷射混凝土衬砌支护？谢谢帮忙。
new
plot add block group
plot show
; 产生网格体
;巷道上部混凝土衬砌
gen zon cshell p0 0 0 0 p1 3.1 0 0 p2 0 50 0 p3 0 0 3.1 dim 2.95 2.95 2.95 2.95 size 2 50 10 group 'concrete liner1'
;巷道上部
gen zon cyl p0 0 0 0 p1 2.95 0 0 p2 0 50 0 p3 0 0 2.95size 5 50 10 group tunnela
;; 巷道下部
gen zone brick p0 0 0 -1.6 p1 add 2.95 0 0 p2 add 0 50 0 p3 add 0 0 1.6 size 5 50 6 group tunnelb
;;巷道两帮混凝土衬砌
gen zone brick p0 2.95 0 -1.6 p1 add 0.15 0 0 p2 add 0 50 0 p3 add 0 0 1.6 size 2 50 6group 'concrete liner2'

; 围岩
;;
gen zone brick p0 0 0 -9p1 8.571 0 -9p2 add 0 50 0 p3 0 0 -1.6 p4 8.571 50 -9 p5 0 50 -1.6 p6 2.95 0 -1.6 p7 2.95 50 -1.6 size 5 50 8 ratio 1 1 0.769 group '围岩1'

gen zone brick p0 8.571 0 -9p1 12 0 -9p2 add 0 50 0 p3 2.95 0 -1.6p4 12 50 -9 p5 2.95 50 -1.6 p6 3.1 0 -1.6 p7 3.1 50 -1.6 size 2 50 8 ratio 1 1 0.769 group '围岩2'

gen z radc p0 0 0 0 p1 12 0 0 p2 0 50 0 p3 0 0 9 dim 3.1 3.1 3.1 3.1 size 5 50 10 8 rat 1 1 1 1.3 group '围岩3'

gen zone brick p0 3.1 0 -1.6 p1 12 0 -9p2 add 0 50 0 p3 3.1 0 0 p4 12 50 -9 p5 3.1 50 0 p6 12 0 0 p7 12 50 0 size 8 50 6 ratio 1.3 1 1 group '围岩4'

gen zone brick p0 0 0 9 p1 12 0 9p2 add 0 50 0 p3 add 0 0 10size 5 50 4 ratio 1 1 1.2 group '围岩上1'
gen zone brick p0 12 0 9 p1 25 0 9p2 add 0 50 0 p3 add 0 0 10 size 3 50 4 ratio 1.2 1 1.2 group '围岩上2'

gen zone brick p0 12 0 0 p1 25 0 0p2 add 0 50 0 p3 add 0 0 9 size 3 50 5 ratio 1.2 1 1 group '围岩右1'

gen zone brick p0 12 0 -9 p1 25 0 -9p2 add 0 50 0 p3 add 0 0 9 size 3 50 6 ratio 1.2 1 1 group '围岩右2'

gen zone brick p0 0 0 -21 p1 12 0 -21p2 add 0 50 0 p3 add 0 0 12 size 7 50 4 ratio 1 1 0.769 group '围岩下1'
gen zone brick p0 12 0 -21 p1 25 0 -21p2 add 0 50 0 p3 add 0 0 12 size 3 50 4 ratio 1.2 1 0.769 group '围岩下2'

gen merge 0.02

;; 指定 Mohr-Coulomb 材料模型
mo mohr
pro den 2200 bulk 2.68e10 she 7e9 fric 38 coh 3.8e6 ten 4.86e6range z -21 -9
pro den 2300 bulk 1.57e10 she 1.078e10 fric 40 coh 3e6 ten 2.31e6range z -9 9
pro den 2100 bulk 2.56e10 she 8.52e9 fric 34 coh 4e6 ten 1.29e6 range z 9 19
;; 边界条件
fix x range x -.1 .1
fix x range x 24.9 25.1
fix y range y -.1 .1
fix y range y 49.9 50.1
fix z range z -9.1 -8.9
;; 初始应力状态
set grav 0 0 -10
apply szz -23e6 range z 18.9 19.1
ini sxx -25e6
ini syy -30e6
ini szz -17e6
hist unbal
plot cont dis out on
so
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
ini sta 0
sav 'maogana.sav
;;;开挖
mo null range group tunnela
mo null range group tunnelb

;;安装锚杆长度为2400 800*1000 共51排
def place_cables
    loop n (1,51)
         y_d = float(n - 1)
      command
            sel cable id n beg 0.000,y_d,2.95 end 0.000,y_d,5.35 nseg 8
            sel cable id n beg 0.793,y_d,2.841 end 1.438,y_d,5.153 nseg 8
            sel cable id n beg 1.527,y_d,2.524 end 2.77,y_d,4.577 nseg 8
            sel cable id n beg 2.149,y_d,2.021 end 3.898,y_d,3.664 nseg 8
            sel cable id n beg 2.613,y_d,1.369 end 4.739,y_d,2.482 nseg 8
            sel cable id n beg 2.885,y_d,0.616 end 5.232,y_d,1.117 nseg 8
            sel cable id n beg 2.950,y_d,-0.182 end 5.340,y_d,-0.33 nseg 8
            sel cable id n beg 2.95,y_d,-1.3 end 5.268,y_d,-1.921 nseg 8            
sel cable id n pro emod 2.06e11 xcarea 1.57e-3 gr_per 1. yten 25e4 gr_k 17.5e6 gr_c 20e4
      end_command
   end_loop
end
place_cables


plot grid sel geom red
;
; 求解
;监控变量
his nstep=5
his unb
hist gp zdisp 0 0 2.95
hist gp xdisp 0 4 2.95
hist gp zdisp 0 8 2.95
hist gp zdisp 0 12 2.95
hist gp zdisp 0 16 2.95
hist gp xdisp 0 20 2.95
hist gp zdisp 0 24 2.95
hist gp zdisp 0 28 2.95
hist gp zdisp 0 32 2.95
hist gp zdisp 0 36 2.95
hist gp zdisp 0 40 2.95
hist gp zdisp 0 44 2.95
hist gp zdisp 0 48 2.95
hist gp zdisp 0 50 2.95
plot cont dis out on
;求解
set mech force 50
;设置大变形
set large
solve
save 'maoganb'.sav

lookcity 发表于 2009-10-14 14:41:12

纵向50m一次性开挖与实际不符
喷层用sel liner range 指定一个范围就可以了然后指定属性

hwaust 发表于 2009-10-16 17:12:15

好的，谢谢

hwaust 发表于 2009-10-20 17:42:24

楼上兄弟，我用sel liner range 进行范围设置怎么求解老是出错，你能不能帮我把命令流补充完整，谢谢？
还有sel liner range 和sel shell range 定义对结果有没有影响？急。。。。。。

hwaust 发表于 2009-10-30 10:08:05

楼上兄弟帮帮忙

liaozhongzy 发表于 2010-5-8 17:45:14

看过感觉还是很有用的，谢帖了

coastjia 发表于 2010-5-11 08:59:49

学习一下，以后也许会用到

fanggaokui 发表于 2010-5-20 19:50:59

类似的问题迷茫中

haixing123 发表于 2010-6-2 15:46:19

正在迷茫中

daiguafan 发表于 2010-6-2 21:34:54

还不是很会啊

beikeren 发表于 2010-6-4 09:56:53

大家都是高手啊，学习一下:victory:

lforever 发表于 2012-1-3 17:08:29

学习下1

lt20070169 发表于 2012-10-6 19:30:58

想请教哈前辈锚杆的位置怎么定，就是cable beginend 的x y z怎么定。。谢谢

查看完整版本: 锚杆支护模拟后，如何进行后处理