求助！！初始应力平衡之后，达到如下的结果 ，请高手指点一下 [复制链接]

代码

初始网格和初始平衡图分别如下

能详细说明一下问题吗？就云图看，可能是网格划分的问题。

可能是网格的连接问题

我也觉得可能是网格划分有问题，但是，我看别人做的也是网格疏密不同，怎么没有出现这种问题呢？  还有，会不会和参数的选择有关系？  谢谢 ！！！命令如下，请多多指教，谢谢你~~~

网格都对齐了呀~~，也attach了. 程序如下，请多多指教~~

new
plot block state
plot show

;老底，细砂岩17m
gen zone brick size 10 30  3 &
p0   0 0 0  p1 40 0 0 p2 0 600 0   p3 0 0 15  ratio 1 1 0.85
gen zone brick size 5 30  3 &
p0   40 0 0  p1 70 0 0 p2 40 600 0   p3 40 0 15  ratio 1 1 0.85
gen zone brick size 5 30  3 &
p0   70 0 0  p1 90 0 0 p2 70 600 0   p3 70 0 15  ratio 1 1 0.85
gen zone brick size 28 30  3 &
p0   90 0 0  p1 230 0 0 p2 90 600 0   p3 90 0 15  ratio 1 1 0.85


;老底，细砂岩8m
gen zone brick size 6 30  2 &
p0   0 0 15  p1 24 0 15 p2 0 600 15   p3 0 0 25  
gen zone brick size 8 30  4 &
p0   24 0 15  p1 40 0 15 p2 24 600 15   p3 24 0 25  
gen zone brick size 15 30  4 &
p0   40 0 15  p1 70 0 15 p2 40 600 15   p3 40 0 25
gen zone brick size 10 30  4 &
p0   70 0 15  p1 90 0 15 p2 70 600 15   p3 70 0 25  
gen zone brick size 28 30  2 &
p0   90 0 15  p1 230 0 15 p2 90 600 15   p3 90 0 25  

;直接底，泥岩3m
gen zone brick size 6 30  1 &
p0   0 0 25  p1 24 0 25 p2 0 600 25   p3 0 0 28
gen zone brick size 8 30  1 &
p0   24 0 25  p1 40 0 25 p2 24 600 25   p3 24 0 28
gen zone brick size 30 30  3 &
p0   40 0 25  p1 70 0 25 p2 40 600 25   p3 40 0 28
gen zone brick size 10 30  1 &
p0   70 0 25  p1 90 0 25 p2 70 600 25   p3 70 0 28
gen zone brick size 28 30  1 &
p0   90 0 25  p1 230 0 25 p2 90 600 25   p3 90 0 28

;3下煤层，6m
gen zone brick size 6 30  2 &
p0   0 0 28  p1 24 0 28 p2 0 600 28   p3 0 0 34
gen zone brick size 8 30  2 &
p0   24 0 28  p1 40 0 28 p2 24 600 28   p3 24 0 34  
gen zone brick size 30 30  6 &
p0   40 0 28  p1 70 0 28 p2 40 600 28   p3 40 0 34
gen zone brick size 10 30  2 &
p0   70 0 28  p1 90 0 28 p2 70 600 28   p3 70 0 34
gen zone brick size 28 30  2 &
p0   90 0 28  p1 230 0 28 p2 90 600 28   p3 90 0 34

;伪顶，砂质泥岩2m
gen zone brick size 6 30  1 &
p0   0 0 34  p1 24 0 34 p2 0 600 34   p3 0 0 36
gen zone brick size 8 30  1 &
p0   24 0 34  p1 40 0 34 p2 24 600 34   p3 24 0 36
gen zone brick size 30 30  2 &
p0   40 0 34  p1 70 0 34 p2 40 600 34   p3 40 0 36
gen zone brick size 10 30  1 &
p0   70 0 34  p1 90 0 34 p2 70 600 34   p3 70 0 36
gen zone brick size 28 30  1 &
p0   90 0 34  p1 230 0 34 p2 90 600 34   p3 90 0 36

;直接顶，中砂岩6m
gen zone brick size 6 30  2 &
p0   0 0 36  p1 24 0 36 p2 0 600 36   p3 0 0 42
gen zone brick size 8 30  2 &
p0   24 0 36  p1 40 0 36 p2 24 600 36   p3 24 0 42
gen zone brick size 30 30  6 &
p0   40 0 36  p1 70 0 36 p2 40 600 36   p3 40 0 42
gen zone brick size 10 30  2 &
p0   70 0 36  p1 90 0 36 p2 70 600 36   p3 70 0 42
gen zone brick size 28 30  2 &
p0   90 0 36  p1 230 0 36 p2 90 600 36   p3 90 0 42

;基本顶，细砂岩8m
gen zone brick size 6 30  2 &
p0   0 0 42  p1 24 0 42 p2 0 600 42   p3 0 0 50
gen zone brick size 8 30  4 &
p0   24 0 42  p1 40 0 42 p2 24 600 42   p3 24 0 50
gen zone brick size 15 30  4 &
p0   40 0 42  p1 70 0 42 p2 40 600 42   p3 40 0 50
gen zone brick size 10 30  4 &
p0   70 0 42  p1 90 0 42 p2 70 600 42   p3 70 0 50
gen zone brick size 28 30  2 &
p0   90 0 42  p1 230 0 42 p2 90 600 42   p3 90 0 50

;基本顶，粉砂岩40m
gen zone brick size 10 30  7 &
p0   0 0 50  p1 40 0 50 p2 0 600 50   p3 0 0 90  ratio 1 1 1.15
gen zone brick size 5 30  7 &
p0   40 0 50  p1 70 0 50 p2 40 600 50   p3 40 0 90  ratio 1 1 1.15
gen zone brick size 5 30  7 &
p0   70 0 50  p1 90 0 50 p2 70 600 50   p3 70 0 90  ratio 1 1 1.15
gen zone brick size 28 30  7 &
p0   90 0 50  p1 230 0 50 p2 90 600 50   p3 90 0 90  ratio 1 1 1.15

attach face
save 初始网格.sav

restore 初始网格.sav
set log on
set logfile Pingheng.TXT

plot block state
plot show

;本构关系、参数
model mohr

macro 砂质泥岩   'shear 0.5e9 bulk 0.7e9 cohesion 1.5e6 fric 28 tension 1.2e6 density 2500'
macro 煤         'shear 0.35e9 bulk 0.5e9 cohesion 1.3e6 fric 28 tension 0.9e6 density 1350'
macro 中砂岩     'shear 1.5e9 bulk 1.8e9 cohesion 3.5e6 fric 30 tension 2e6 density 2500'
macro 细砂岩     'shear 2.0e9 bulk 4.0e9 cohesion 6e6 fric 38 tension 5e6 density 2500'
macro 粉砂岩     'shear 2.0e9 bulk 4.0e9 cohesion 6e6 fric 38 tension 5e6 density 2500'

prop 粉砂岩        range z 50 90
prop 细砂岩        range z 42 50
prop 中砂岩        range z 36 42
prop 砂质泥岩      range z 34 36
prop 煤            range z 28 34
prop 砂质泥岩      range z 25 28
prop 细砂岩        range z 0 25

;边界条件
fix z range z  -0.1 0.1
fix x range x  -0.1 0.1
fix x range x  229.9 230.1
fix y range y  -0.1 0.1
fix y range y  599.9 600.1

;在上面z=90处施加力
apply szz -21.5e6 range z 89.9 90.1

;重力梯度化
set grav 0 0 -9.8
ini szz -21.5e6 grad 0 0 2.5e4
ini syy -10.55e6 grad 0 0 1.25e4
ini sxx -10.55e6 grad 0 0 1.25e4

hist unbal
plot set plane nor 0 1 0 origin 0 100 0
plot con szz   plan

solve
save 初始平衡.sav
;没开挖前的的自重平衡
set log off

图片:未命名.JPG

图片:未命名2.JPG

刚才将你的命令设置分组运行后的结果图，现有一疑问，你为什么要将中间那块的网格划分那么密呢？还有感觉你建模的命令好多都是重复性的，为什么不直接用一个gen  zon 来建立每层呢？而是每层都用好几个gen zon ? 这样不麻烦吗？

煤层的密度与各层不同，你设置的初始应力梯度在煤层处必然有突变啊 怎么能在整个底层中均匀变化呢？

网格细化的部分肯定有细部结构

谢谢帮忙哈~我现在要建立的模型是，加密的那个部分是巷道开挖的位置，右边的部分是工作面回采的位置。我主要就是想研究的就是工作面回采对巷道的影响。所以，我想巷道周边（包括顶底板、两帮）都要加密啊，不知道是不是需要这样？
我用的gen  zon 多是因为每一层不同块的网格疏密不同呀~~~