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

见了一个分离式隧道的模型，模拟自重应力场竖向应力成这样的了。不知道怎么回事，竖向位移也很莫名其妙的。 -
Q JPJ.  
程序如下。 (jj`}Qe3U  
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 @n.n[zb\|  
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 +WMXd.iN,  
gen zon reflect dip 90 dd 270 orig 7.9 0 0 i|AWaG)  
;s3\Z^h4kd  
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 v@q&B|0  
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 |#<z\u }  
gen zon reflect dip 90 dd 270 orig 7.9 0 0 range group liningyx any group inrockyx any group outrockyx any ]d7A|)q  
` V [4  
8Yf*vp>T/x  
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 u7RlxA:  
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 (s&]V49  
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 sP2Uj  
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 qhcx\eD:?  
ZS(%!+M  
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 |&W4Dkn  
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 Z}>F
V~4  
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 e`LkCy[_  
_(8#  
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 vxC];nCC#  
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 oPni4^g 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 4Otq3s34FT  
zaLPPm&f  
gen zon reflect dip 90 dd 270 orig 0 0 0 t^zE^:06  
plot block group *>mjUT}cP  
_xm<zy{`S  
set gravity 0 0 -10 5XB]p|YU~s  
fix z range z -24.01,-23.99 Do\j_  
fix x range x -65.01,-64.99 MMpId Uhr  
fix x range x 64.99,65.01 .Tq8Qdl  
fix y range y -0.01,0.01 '7oCWHq[  
fix y range y 0.99,1.01 8]b;l; W5  
model mohr _&9P&Zf4  
ini density 2400 _E'}8.#{  
prop bulk=6e9,shear=10e9 fric=45,coh=1.2e6 tension=15e5 [TUs^%2@  
set mech ratio=5e-6 V]+y*b.60  
solve Q3&DA1b`  
拜托了，在线等。

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

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

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

初始应清零的~ ,VEE<*'X  
然后再按照施工步序模拟~

不知道  不应该啊

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

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

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

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