求高手看一段flac程序 [复制链接]

见了一个分离式隧道的模型，模拟自重应力场竖向应力成这样的了。不知道怎么回事，竖向位移也很莫名其妙的。 niV=Ijt{5  
程序如下。 +kKfx!  
gen zon radcyl p0 7.9 0 0 p1 14.3 0 0 p2 7.9 1 0 p3 7.9 0 6.4 size 5 2 10 4 dim 5.6 5.6 5.6 5.6 rat 1 1 1 1.2 group outrockys YQ}bG{V  
gen zon cshell p0 7.9 0 0 p1 13.8 0 0 p2 7.9 1 0 p3 7.9 0 5.9 size 1 2 10 4 dim 5.3 5.3 5.3 5.3 rat 1 1 1 1 group liningys fill group inrockys QeVM9br)m  
gen zon reflect dip 90 dd 270 orig 7.9 0 0 NPM}w!  
/a$RJ6t&3  
gen zon radcyl p0 7.9 0 0 p1 7.9 0 -4 p2 7.9 1 0 p3 14.3 0 0 size 5 2 10 4 dim 3.5 5.6 3.5 5.6 rat 1 1 1 1.2 group outrockyx Vee`q.  
gen zon cshell p0 7.9 0 0 p1 7.9 0 -3.5 p2 7.9 1 0 p3 13.8 0 0 size 1 2 10 4 dim 3.2 5.3 3.2 5.3 rat 1 1 1 1 group liningyx fill group inrockyx 6uf+,F  
gen zon reflect dip 90 dd 270 orig 7.9 0 0 range group liningyx any group inrockyx any group outrockyx any  D=nuK25  
e&(Di,%:  
;yO7!
{_  
gen zon brick p0 0 0 -4 p1 1.5 0 -4 p2 0 1 -4 p3 0 0 0 size 3 2 5 group outrock1 jz2W/EE`w  
gen zon brick p0 14.3 0 -4 p1 65 0 -4 p2 14.3 1 -4 p3 14.3 0 0 size 20 2 5 group outrock1 +<P%v k  
gen zon brick p0 0 0 0 p1 1.5 0 0 p2 0 1 0 p3 0 0 6.4 size 3 2 5 group outrock2 438+zU  
gen zon brick p0 14.3 0 0 p1 65 0 0 p2 14.3 1 0 p3 14.3 0 6.4 size 20 2 5 group outrock2 )5U&^tJ  
w%[`'_[  
gen zon brick p0 0 0 6.4 p1 1.5 0 6.4 p2 0 1 6.4 p3 0 0 36 size 3 2 10 group outrock3 T=w5FT  
gen zon brick p0 14.3 0 6.4 p1 65 0 6.4 p2 14.3 1 6.4 p3 14.3 0 36 size 20 2 10 group outrock3 T7=~l)I  
gen zon brick p0 1.5 0 6.4 p1 14.3 0 6.4 p2 1.5 1 6.4 p3 1.5 0 36 size 10 2 10 group outrock3 +[ 9
44n  
TY3WP
$u  
gen zon brick p0 0 0 -24 p1 1.5 0 -24 p2 0 1 -24 p3 0 0 -4 size 3 2 8 group outrock4 =?f\o*J)  
gen zon brick p0 14.3 0 -24 p1 65 0 -24 p2 14.3 1 -24 p3 14.3 0 -4 size 20 2 8 group outrock4 I)Dd"I  
gen zon brick p0 1.5 0 -24 p1 14.3 0 -24 p2 1.5 1 -24 p3 1.5 0 -4 size 10 2 8 group outrock4 3g]Sp/  
n1xN:A  
gen zon reflect dip 90 dd 270 orig 0 0 0 &dkjT8L$  
plot block group fK|F`F2V  
j6KGri  
set gravity 0 0 -10 jnuovM!x~  
fix z range z -24.01,-23.99 fN TPW]  
fix x range x -65.01,-64.99 *a7&v3X  
fix x range x 64.99,65.01 I2=?H<  
fix y range y -0.01,0.01 u@$C i/J*  
fix y range y 0.99,1.01 }*4K]3et$  
model mohr Q1Sf
7)  
ini density 2400 3%<xM/#  
prop bulk=6e9,shear=10e9 fric=45,coh=1.2e6 tension=15e5 Py[
Z9KLX  
set mech ratio=5e-6 JYB<}
;,  
solve Y&k6Xhuao  
拜托了，在线等。

还没有开挖怎么就有隧道出现了呢

估计是网格划分的原因，只要数值差异不大，开挖前全部归零就是了

看下数值是不是差10的负几次的数量级

初始应清零的~ Brh<6Btl  
然后再按照施工步序模拟~

不知道  不应该啊

在初始应力时最好用这个解算solve elas

还没有应力初始化，就开挖了吧

这个是还没有开挖呢，就是模拟自重力场的时候就这样了。

计算精度的设置应该有关系，可以设置到1e-5看看