【求助】求助！初始应力场，哪里错了呢？？
本帖最后由 haozhenggong 于 2011-1-5 08:29 编辑初始应力场，埋深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 0p5 0 20 1 p6 15 0 3 p7 15 20 3 group 1
gen zonebrick size 15 20 6 p0 0 0 1 p1 15 0 3 p2 0 20 1 p3 0 0 3 p4 15 20 3p5 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.5e9bulk 3.3e9cohesion 4.0e6 fric 37ten 1.20e6'
macro niyan 'shear 2.07e9 bulk 2.2e9cohesion 1.35e6 fric 30ten 0.5e6'
macro xishayan 'shear 1.60e9 bulk 2.7e9cohesion 2.0e6 fric 35ten 1.0e6'
macro mei 'shear 4.60e9 bulk 1.0e9cohesion 1.2e6 fric 28ten 0.6e6'
prop zhongshayan range group 4
prop xishayan range group 2
prop mei range group 3
prop xishayan range group 1
ini density1430 rangegroup 3
ini density2580 rangegroup 4
ini density2540 rangegroup 2
ini density2540 rangegroup 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 0p5 0 20 1 p6 15 0 3 p7 15 20 3 group 1
gen zonebrick size 15 20 6 p0 0 0 1 p1 15 0 3 p2 0 20 1 p3 0 0 3 p4 15 20 3p5 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.5e9bulk 3.3e9cohesion 4.0e6 fric 37ten 1.20e6'
macro niyan 'shear 2.07e9 bulk 2.2e9cohesion 1.35e6 fric 30ten 0.5e6'
macro xishayan 'shear 1.60e9 bulk 2.7e9cohesion 2.0e6 fric 35ten 1.0e6'
macro mei 'shear 4.60e9 bulk 1.0e9cohesion 1.2e6 fric 28ten 0.6e6'
prop zhongshayan range group 4
prop xishayan range group 2
prop mei range group 3
prop xishayan range group 1
ini density1430 rangegroup 3
ini density2580 rangegroup 4
ini density2540 rangegroup 2
ini density2540 rangegroup 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.5e4range 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
这个应力图更看不懂了
两者为什么不一样呢，还有 我的初始应力场设置的对不对呢？？？高手解答下 个人觉得采用S-B法初始地应力平衡时，还应该加上Z方向的底面约束，即:fix z ran z -.1.1. 本帖最后由 haozhenggong 于 2011-1-8 21:13 编辑
同意楼上hello_kitty ！！！！！！！！ S-B法初始地应力平衡时，还应该加上Z方向的底面约束，即:fix z ran z -.1.1. 你的问题解决了吗？ 这种结果是参数的问题，当你把参数都换成一样的以后，你会发现szz会呈梯度变化；但是我发现原岩竟然能算平衡，这不合理啊，按道理模型会一直往下沉，不会平衡的。 试一试:handshake:handshake:handshake S-B法更适合 楼主 问一下 你的SB法 怎摸改善的？？？谢谢了
页:
[1]