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

见了一个分离式隧道的模型，模拟自重应力场竖向应力成这样的了。不知道怎么回事，竖向位移也很莫名其妙的。 T*C F5S  
程序如下。 #)48dW!n  
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 UH+#Nel+!  
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 r4X\/  
gen zon reflect dip 90 dd 270 orig 7.9 0 0 R^$EnrY(<  
*V+,X  
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 xI#rnx*  
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 7h41E#  
gen zon reflect dip 90 dd 270 orig 7.9 0 0 range group liningyx any group inrockyx any group outrockyx any 3,*A VcQA  
]BRwJ2< x  
m
MWhUr  
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 jA~omX2A  
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 >o&%via}  
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 1P'L<z  
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 $nIE;idk  
t,0}}9%?  
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 s[/d}S@ >  
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 OUO'w6m!  
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 "s_S!;w@  
QGGBI Ku   
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 U#w0E G  
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 iYw1{U  
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 55] MRv  
e.XD5~Ax  
gen zon reflect dip 90 dd 270 orig 0 0 0 kJNg>SN*@#  
plot block group %Q.M& U  
SJ;{ Hg  
set gravity 0 0 -10 #}~?8/h!  
fix z range z -24.01,-23.99 R*k;4*1u  
fix x range x -65.01,-64.99 )S`Yl;oL  
fix x range x 64.99,65.01 U;u4ey  
fix y range y -0.01,0.01 42,dH
Ydt  
fix y range y 0.99,1.01 H<g8u{ $  
model mohr a~JZc<ze  
ini density 2400 *r9D+}Y(4  
prop bulk=6e9,shear=10e9 fric=45,coh=1.2e6 tension=15e5 4&e<Sc64  
set mech ratio=5e-6 R6z *!W{  
solve n
? =O@yq  
拜托了，在线等。

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

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

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

初始应清零的~ ;1r|Bx
<5  
然后再按照施工步序模拟~

不知道  不应该啊

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

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

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

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