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

zhoubo8708

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

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

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


gen zone brick p0 0 0  6.5  p1 add 6.5 0 0  p2 add 0 5 0 p3 add 0 0 7                       size 6 5 7 group s-tutl-1
gen zone brick p0 6.5 0  6.5  p1 add 18 0 0  p2 add 0 5 0 p3 add 0 0 7                       size 18 5 7 group s-tutl-1
gen zone brick p0 0  0  -20  p1 add 6.5 0 0  p2 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  -20  p1 add 18 0 0  p2 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

songli269

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

liuji2009

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

chicprince

边界条件要设置好

silinfe

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

msh2010

楼主好想法

love36524

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

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

ruanqingsong

wwwwwwwwwwwwwwwwwww

cumtdabing

可以输出锚杆szz应力么
？

gaochao20wo

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

ycs008

我也遇到了这样的问题

ycs008

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