论坛 › G0201：FLAC3D/FLAC2D › 隧道开挖（配有程序和云图）

zhoubo8708 发表于 2010-7-8 08:58:33

隧道开挖（配有程序和云图）

以下就是我的程序，请各位看看：存在一个问题，在加了壳单元的情况下，整个地表是向上隆起的，请多多指教！！

gen zone radcyl   p0 0 0 0   p1 6.5 0 0   p20 5 0   p30 0 6.5dim 3.1 3.1 3.1 3.1   size 3 5 12 3ratio 1 1 1 1group tutl fill group suidao
gen zone reflect normal 0 0 -1

gen zone brick p0 6.5 0 -6.5p1 add 18 0 0p2 add 0 5 0 p3 add 0 0 13                        size 18 5 12 group ytutl


gen zone brick p0 0 06.5p1 add 6.5 0 0p2 add 0 5 0 p3 add 0 0 7                     size 6 5 7 group s-tutl-1
gen zone brick p0 6.5 06.5p1 add 18 0 0p2 add 0 5 0 p3 add 0 0 7                     size 18 5 7 group s-tutl-1
gen zone brick p0 00-20p1 add 6.5 0 0p2 add 0 5 0 p3 add 0 0 13.5                        size 6 5 14 group x-tutl-1
gen zone brick p0 6.5 0-20p1 add 18 0 0p2 add 0 5 0 p3 add 0 0 13.5                        size 18 5 14 group x-tutl-1
model mohr
property bu 2.5e7 sh 7e6 fric 30 cohesion 2.1e4 tension 2e4
in de 1900
set gravity 0 0-9.81

fix x range x 24.5
fix x range x 0
fix y range y 0
fix y range y 5
fix z range z -20
set mech ratio=1e-5
solve
save quan-moxing.sav
in xdisp=0 ydisp=0 zdisp=0
def excate_step1
n=1
a1=n-1
a2=n

end
excate_step1
mod null range cyl end1 0 a1 0 end2 0 a2 0 rad 3.1
sel shell id=1 range cyl end1 0 a1 0 end2 0 a2 0 rad 3.1 group suidao not
sel shell id=1 property iso=(3.45e10,0.3) thick=0.35 density 2450   
                                                
set mech ratio=1e-5
solve
save section1.sav

def excate_step2
n=2
a1=n-1
a2=n
end
excate_step2
mod null range cyl end1 0 a1 0 end2 0 a2 0 rad 3.1
sel shell id=1 range cyl end1 0 a1 0 end2 0 a2 0 rad 3.1 group suidao not

sel shell id=1 property iso=(3.45e10,0.3) thick=0.35 density 2450                                             
set mech ratio=1e-5
solve
save section2.sav
def excate_step3
n=3
a1=n-1
a2=n

end
excate_step3
mod null range cyl end1 0 a1 0 end2 0 a2 0 rad 3.1
sel shell id=1 range cyl end1 0 a1 0 end2 0 a2 0 rad 3.1 group suidao not
sel shell id=1 property iso=(3.45e10,0.3) thick=0.35 density 2450   
set mech ratio=1e-5
solve
save section3.sav
def excate_step4
n=4
a1=n-1
a2=n
end
excate_step4
mod null range cyl end1 0 a1 0 end2 0 a2 0 rad 3.1
sel shell id=1 range cyl end1 0 a1 0 end2 0 a2 0 rad 3.1 group suidao not

sel shell id=1 property iso=(3.45e10,0.3) thick=0.35 density 2450                                                   
set mech ratio=1e-5
solve
save section4.sav
def excate_step5
n=5
a1=n-1
a2=n

end
excate_step5
mod null range cyl end1 0 a1 0 end2 0 a2 0 rad 3.1
sel shell id=1 range cyl end1 0 a1 0 end2 0 a2 0 rad 3.1 group suidao not
sel shell id=1 property iso=(3.45e10,0.3) thick=0.35 density 2450                                                   
set mech ratio=1e-5
solve
save section5.sav

zhoubo8708 发表于 2010-7-8 09:05:04

songli269 发表于 2010-8-26 13:07:17

对称结构取一半计算时，初始要控制转动。

liuji2009 发表于 2010-11-10 10:53:05

你这个是几步开挖？什么工法

chicprince 发表于 2010-11-10 12:57:30

边界条件要设置好

silinfe 发表于 2010-12-9 20:51:31

我遇到地表隆起这种情况
如何解决呢？

msh2010 发表于 2011-3-7 10:48:59

本帖最后由 msh2010 于 2011-3-21 10:43 编辑

楼主好想法

love36524 发表于 2011-3-7 17:12:30

我想请问一下，有没有在别的软件里面建模，然后导入flac里面计算的/?

这样编程序也忒麻烦了( ⊙ o ⊙ )啊！

ruanqingsong 发表于 2011-5-25 09:58:42

wwwwwwwwwwwwwwwwwww

cumtdabing 发表于 2011-9-25 16:55:37

:) 可以输出锚杆szz应力么
？

gaochao20wo 发表于 2011-9-26 10:46:04

想学习一下楼主，请问楼主的QQ号多少？我的496123461

ycs008 发表于 2011-11-30 16:43:08

我也遇到了这样的问题

ycs008 发表于 2011-12-1 21:58:47

加我qq394343652，我也遇到你这样的问题，我们一起交流交流

查看完整版本: 隧道开挖（配有程序和云图）