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

见了一个分离式隧道的模型，模拟自重应力场竖向应力成这样的了。不知道怎么回事，竖向位移也很莫名其妙的。 PcC@}3  
程序如下。 R ABw(b  
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 _Z(t**Zh6y  
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 0BDw}E\  
gen zon reflect dip 90 dd 270 orig 7.9 0 0 1dLc/,|  
b$klm6nMvm  
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 (T*$4KGV  
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 k\[(;9sf.  
gen zon reflect dip 90 dd 270 orig 7.9 0 0 range group liningyx any group inrockyx any group outrockyx any 9:p-F+  
&IN%2c  
Aax;0qGbH  
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 O2>c|=#  
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 l~"T>=jq3  
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 5TJd9:\Af  
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 SAdT#0J  
bY#BK_8 :  
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$Vv  
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 Dy.i^`7\  
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 }. &ellNQ  
@$jV"Y  
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 ir<e^a  
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 y.lWyH9  
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 aqEZhMy  
qt}M&=}8Q  
gen zon reflect dip 90 dd 270 orig 0 0 0 41<~_+-@  
plot block group >[Vc$[62  
~)f^y!PMQ  
set gravity 0 0 -10 ;p+'?%Y}  
fix z range z -24.01,-23.99 ./ {79  
fix x range x -65.01,-64.99 To(I<W|{  
fix x range x 64.99,65.01 Kn:Ml4[;  
fix y range y -0.01,0.01 zx%WV@O9  
fix y range y 0.99,1.01 #DgHF*GG+>  
model mohr V<UChD)N`  
ini density 2400 e%cTFwX?n  
prop bulk=6e9,shear=10e9 fric=45,coh=1.2e6 tension=15e5 J'Pyn  
set mech ratio=5e-6 N[@H107`  
solve vS\2zwb}  
拜托了，在线等。

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

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

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

初始应清零的~ SyVGm@  
然后再按照施工步序模拟~

不知道  不应该啊

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

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

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

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