求助隧道求解自重应力场出错

luolibra

new
;1/4上半部分
gen zon radcyl p0 0 0 0 p1 11.06 0 0 p2 0 10 0 p3 0  0 8.4 &
                        size 5 2 10 4 dim 5.8 5.8 5.8 5.8 rat 1 1 1 1.2 group outrockys
gen zon cshell p0 0 0 0 p1 5.8 0 0 p2 0 10 0 p3 0 0 5.8  &
                        size 1 2 10 4 dim 5.4 5.4 5.4 5.4 rat 1 1 1 1 group liningys fill group inrockys
gen zon reflect dip 90 dd 270 orig 0 0 0
;1/4下半部分
gen zon radcyl p0 0 0 0 p1 0 0 -6.9 p2 0 10 0 p3 11.06 0 0 &
                        size 1 2 10 4 dim 3.9 5.8 3.9 5.8 rat 1 1 1 1.2 group outrockyx
gen zon cshell p0 0 0 0 p1 0 0 -3.9 p2 0 10 0 p3 5.8 0 0 &
                        size 1 2 10 4 dim 3.5 5.4 3.5 5.4 rat 1 1 1 1 group liningyx fill group inrockyx
gen zon reflect dip 90 dd 270 orig 0 0 0 range group liningyx any &
                        group inrockyx any group outrockyx any
;隧道中地层网格
gen zon brick p0 -11.06 0 -6.9 p1 -44.24 0 -6.9 p2 -11.06 10 -6.9 p3 -11.06 0 0 &
                        size 6 2 5 group outrock1
gen zon brick p0 11.06 0 -6.9 p1 44.2 0 -6.9 p2 11.06 10 -6.9 p3 11.06 0 0 size 6 2 5 group outrock1
gen zon brick p0 -11.06 0 8.4 p1 -44.24 0 8.4 p2 -11.06 10 8.4 p3 -11.06 0 0 &
                        size 6 2 5 group outrock2
gen zon brick p0 11.06 0 8.4 p1 44.24 0 8.4 p2 11.06 10 8.4 p3 11.06 0 0 size 6 2 5 group outrock2
;隧道上侧土体网格
gen zon brick p0 -44.24 0 8.4 p1 -11.06 0 8.4 p2 -44.24 10 8.4 p3 -44.24 0 28.4 &
                        size 6 2 10 group outrock3
gen zon brick p0 -11.06 0 8.4 p1 11.06 0 8.4 p2 -11.06 10 8.4 p3 -11.06 0 28.4 &
                        size 10 2 10 group outrock3
gen zon brick p0 11.06 0 8.4 p1 44.24 0 8.4 p2 11.06 10 8.4 p3 11.06 0 28.4 &
                        size 6 2 10 group outrock3
;隧道下侧土体网格
gen zon brick p0 -44.24 0 -6.9 p1 -11.06 0 -6.9 p2 -44.24 10 -6.9 p3 -44.24 0 -32.1 &
                        size 6 2 10 group outrock4
gen zon brick p0 -11.06 0 -6.9 p1 11.06 0 -6.9 p2 -11.06 10 -6.9 p3 -11.06 0 -32.1 &
                        size 10 2 10 group outrock4
gen zon brick p0 11.06 0 -6.9 p1 44.24 0 -6.9 p2 11.06 10 -6.9 p3 11.06 0 -32.1 &
                        size 6 2 10 group outrock4


;重力加速度
set gravity 0 0 -10

;给定边界条件
fix z range z -32.09 -32.11
fix x range x -44.25 -44.23
fix x range x 44.23 44.25
fix y range y -0.01 0.01
fix y range y 9.99 10.01

;求解自重应力场
model mohr
ini density 2200
new
;1/4上半部分
gen zon radcyl p0 0 0 0 p1 11.06 0 0 p2 0 10 0 p3 0  0 8.4 &
                        size 5 2 10 4 dim 5.8 5.8 5.8 5.8 rat 1 1 1 1.2 group outrockys
gen zon cshell p0 0 0 0 p1 5.8 0 0 p2 0 10 0 p3 0 0 5.8  &
                        size 1 2 10 4 dim 5.4 5.4 5.4 5.4 rat 1 1 1 1 group liningys fill group inrockys
gen zon reflect dip 90 dd 270 orig 0 0 0
;1/4下半部分
gen zon radcyl p0 0 0 0 p1 0 0 -6.9 p2 0 10 0 p3 11.06 0 0 &
                        size 1 2 10 4 dim 3.9 5.8 3.9 5.8 rat 1 1 1 1.2 group outrockyx
gen zon cshell p0 0 0 0 p1 0 0 -3.9 p2 0 10 0 p3 5.8 0 0 &
                        size 1 2 10 4 dim 3.5 5.4 3.5 5.4 rat 1 1 1 1 group liningyx fill group inrockyx
gen zon reflect dip 90 dd 270 orig 0 0 0 range group liningyx any &
                        group inrockyx any group outrockyx any
;隧道中地层网格
gen zon brick p0 -11.06 0 -6.9 p1 -44.24 0 -6.9 p2 -11.06 10 -6.9 p3 -11.06 0 0 &
                        size 6 2 5 group outrock1
gen zon brick p0 11.06 0 -6.9 p1 44.2 0 -6.9 p2 11.06 10 -6.9 p3 11.06 0 0 size 6 2 5 group outrock1
gen zon brick p0 -11.06 0 8.4 p1 -44.24 0 8.4 p2 -11.06 10 8.4 p3 -11.06 0 0 &
                        size 6 2 5 group outrock2
gen zon brick p0 11.06 0 8.4 p1 44.24 0 8.4 p2 11.06 10 8.4 p3 11.06 0 0 size 6 2 5 group outrock2
;隧道上侧土体网格
gen zon brick p0 -44.24 0 8.4 p1 -11.06 0 8.4 p2 -44.24 10 8.4 p3 -44.24 0 28.4 &
                        size 6 2 10 group outrock3
gen zon brick p0 -11.06 0 8.4 p1 11.06 0 8.4 p2 -11.06 10 8.4 p3 -11.06 0 28.4 &
                        size 10 2 10 group outrock3
gen zon brick p0 11.06 0 8.4 p1 44.24 0 8.4 p2 11.06 10 8.4 p3 11.06 0 28.4 &
                        size 6 2 10 group outrock3
;隧道下侧土体网格
gen zon brick p0 -44.24 0 -6.9 p1 -11.06 0 -6.9 p2 -44.24 10 -6.9 p3 -44.24 0 -32.1 &
                        size 6 2 10 group outrock4
gen zon brick p0 -11.06 0 -6.9 p1 11.06 0 -6.9 p2 -11.06 10 -6.9 p3 -11.06 0 -32.1 &
                        size 10 2 10 group outrock4
gen zon brick p0 11.06 0 -6.9 p1 44.24 0 -6.9 p2 11.06 10 -6.9 p3 11.06 0 -32.1 &
                        size 6 2 10 group outrock4


;重力加速度
set gravity 0 0 -10

;给定边界条件
fix z range z -32.09 -32.11
fix x range x -44.25 -44.23
fix x range x 44.23 44.25
fix y range y -0.01 0.01
fix y range y 9.99 10.01

;求解自重应力场
model mohr
ini density 2200


错误是Zone **** has illegal geometry.
请帮忙解决一下，谢谢。

程小杰

这个报错是Flac3D里面计算的时候第一次平衡造成了网格单元的流行破坏造成的，Flac3D里面的网格单元是不可以发生沙漏破坏的，什么是沙漏破坏，就是本来不想交的两个网格线不可以产生新的交叉，这个容易理解吧。需要解决的方法是第一次平衡可以直接用应力施加，地应力，可以不做第一次solve，直接用重力梯度施加，然后ratio设置xy两个方向的应力比，默认0.5即可。还有一个方法，第一次solve可以用弹性模型，还有一个解决方法，网格大一点，均衡一点。

luolibra

没人帮忙解决一下吗？

chenpf09

估计是建立模型出的问题，你在仔细合一下，模型的建立情况。