【求助】求助！初始应力场，哪里错了呢？？

haozhenggong

初始应力场，埋深800m  按25KN/m 算 侧压比0.5，常规方法：
new
gen zone brick size 15 20 3 p0 0 0 0 p1 15 0 0 p2 0 20 0 p3 0 0 1 p4 15 20 0  p5 0 20 1 p6 15 0 3 p7 15 20 3 group 1
gen zone  brick size 15 20 6 p0 0 0 1 p1 15 0 3 p2 0 20 1 p3 0 0 3 p4 15 20 3  p5 0 20 3 p6 15 0 5 p7 15 20 5 group 2
gen zone brick size 15 20 1 p0 0 0 3 p1 15 0 5 p2 0 20 3 p3 0 0 4 p4 15 20 5 p5 0 20 4 p6 15 0 6 p7 15 20 6 group 3
gen zone brick size 15 20 5 p0 0 0 4 p1 15 0 6 p2 0 20 4 p3 0 0 8 p4 15 20 6 p5 0 20 8 p6 15 0 8 p7 15 20 8 group 4
plot add block group yellow white black blue
Plot set rotation 20 0 30
Plot set center auto
Plot add axes
Plot show
model mohr
macro zhongshayan       'shear 2.5e9  bulk 3.3e9  cohesion 4.0e6   fric 37  ten 1.20e6'
macro niyan             'shear 2.07e9 bulk 2.2e9  cohesion 1.35e6   fric 30  ten 0.5e6'
macro xishayan          'shear 1.60e9 bulk 2.7e9  cohesion 2.0e6   fric 35  ten 1.0e6'
macro mei               'shear 4.60e9 bulk 1.0e9  cohesion 1.2e6     fric 28  ten 0.6e6'
prop zhongshayan          range   group 4
prop xishayan             range   group 2
prop mei                  range    group 3        
prop xishayan             range    group 1

ini density  1430     range  group 3
ini density  2580     range  group 4
ini density  2540     range  group 2
ini density  2540     range  group 1
ini sxx -7.2e6 grad 0 0 1.125e4 range x -.1 .1   
ini sxx -7.2e6 grad 0 0 1.125e4 range x 14.9 15.1
ini syy -7.2e6 grad 0 0 1.125e4 range y -.1 .1
ini syy -7.2e6 grad 0 0 1.125e4 range y 19.9 20.1
ini szz  -16e6   range z -.1 .1
;确定边界条件
fix x y z                  range z -.1  .1
fix x                     range x -0.1  -0.1
fix x                    range x 14.9 15.1
fix y                     range y -0.1 0.1
fix y                    range y 19.9 20.1

set grav 0 0 -10
solve
ini xdis=0 ydis=0 zdis=0
ini xvel 0 yvel 0 zvel 0
plo cont szz  
得到如下图  感觉不对：



S-B方法：：：
new
gen zone brick size 15 20 3 p0 0 0 0 p1 15 0 0 p2 0 20 0 p3 0 0 1 p4 15 20 0  p5 0 20 1 p6 15 0 3 p7 15 20 3 group 1
gen zone  brick size 15 20 6 p0 0 0 1 p1 15 0 3 p2 0 20 1 p3 0 0 3 p4 15 20 3  p5 0 20 3 p6 15 0 5 p7 15 20 5 group 2
gen zone brick size 15 20 1 p0 0 0 3 p1 15 0 5 p2 0 20 3 p3 0 0 4 p4 15 20 5 p5 0 20 4 p6 15 0 6 p7 15 20 6 group 3
gen zone brick size 15 20 5 p0 0 0 4 p1 15 0 6 p2 0 20 4 p3 0 0 8 p4 15 20 6 p5 0 20 8 p6 15 0 8 p7 15 20 8 group 4
plot add block group yellow white black blue
Plot set rotation 20 0 30
Plot set center auto
Plot add axes
Plot show
model mohr
macro zhongshayan       'shear 2.5e9  bulk 3.3e9  cohesion 4.0e6   fric 37  ten 1.20e6'
macro niyan             'shear 2.07e9 bulk 2.2e9  cohesion 1.35e6   fric 30  ten 0.5e6'
macro xishayan          'shear 1.60e9 bulk 2.7e9  cohesion 2.0e6   fric 35  ten 1.0e6'
macro mei               'shear 4.60e9 bulk 1.0e9  cohesion 1.2e6     fric 28  ten 0.6e6'
prop zhongshayan          range   group 4
prop xishayan             range   group 2
prop mei                  range    group 3        
prop xishayan             range    group 1

ini density  1430     range  group 3
ini density  2580     range  group 4
ini density  2540     range  group 2
ini density  2540     range  group 1

app sxx -7.2e6 grad 0 0 1.125e4 range x -.1 .1
app sxx -7.2e6 grad 0 0 1.125e4 range x 14.9 15.1
app syy -7.2e6 grad 0 0 1.125e4 range y -.1 .1
app syy -7.2e6 grad 0 0 1.125e4 range y 19.9 20.1
apply szz -16e6 grad 0 0 2.5e4  range z 0 8
set grav 0 0 -10
solve
ini xdis=0 ydis=0 zdis=0
ini xvel 0 yvel 0 zvel 0
plo cont szz

这个应力图更看不懂了

两者为什么不一样呢，还有 我的初始应力场设置的对不对呢？？？高手解答下

hello_kitty

个人觉得采用S-B法初始地应力平衡时，还应该加上Z方向的底面约束，即:fix z ran z -.1  .1.

haozhenggong

同意楼上hello_kitty      ！！！！！！！！ S-B法初始地应力平衡时，还应该加上Z方向的底面约束，即:fix z ran z -.1  .1.

咬钩鱼

你的问题解决了吗？

咬钩鱼

这种结果是参数的问题，当你把参数都换成一样的以后，你会发现szz会呈梯度变化；但是我发现原岩竟然能算平衡，这不合理啊，按道理模型会一直往下沉，不会平衡的。

止戈之殇

试一试:handshake:handshake

yanling205

S-B法更适合

炫光光

楼主 问一下 你的SB法 怎摸改善的？？？谢谢了