初始地应力场无法收敛

冬日的阳光

大家看看我这命令流有什么问题 为什么计算就是收敛不了

冬日的阳光

顶起来 没人帮忙吗

jihezhilv

抽空给你看了一下：问题还真不少，呵呵
new
;set log on
;==========set group clay======================
;这么简单的模型直接建立就行了，乍还用镜像啊
gen zone brick p0 0 0 -25 p1 25 0 -25 p2 0 25 -25 p3 0 0  -12 &
               p4 25 25 -25  p5 0 25 -12 p6 25 0 -12 &
               p7 25 25  -12 ratio 1.1 1.1 1 size 5 5 2 group clay1
gen zone brick p0 0 0 -37 p1 25 0 -37 p2 0 25 -37 p3 0 0 -25 &
              p4 25 25 -37 p5 0 25 -25 p6 25 0 -25 &
              p7 25 25 -25 ratio 1.1 1.1 1 size 5 5 2 group clay2
gen zone brick p0 0 0 -47 p1 25 0 -47 p2 0 25 -47 p3 0 0 -37 &
              p4 25 25 -47 p5 0 25 -37 p6 25 0 -37 &
              p7 25 25 -37 ratio 1.1 1.1 1 size 5 5 2 group clay3
              
gen zone reflect dip 90 dd 90
gen zone reflect dip 90 dd 180
;plot block group
fix x y z range z -47.1 -46.9
fix x range x 24.9 25.1    ;这里你弄成了24.1
fix x range x -24.9 -25.1
fix y range y 24.9 25.1
fix y range y -24.9 -25.1
model mohr  ;ran z -47 -12
;prop bulk=3.91e6 coh=10e10 tension=1e10 shear=1.5e6 range group clay1
;ini dens 1800 range group clay1
;prop bulk=4.44e6 coh=10e10 tension=1e10 shear=1.48e6 range group clay2
;ini dens 1800 range group clay2
;prop bulk=16.67e6 coh=10e10 tension=1e10 shear=5.56e6 range group clay3
ini dens 1800     ; 既然三个岩性的密度都一样，为什么还分别赋值呢
set grav 0 0 -10
;solve
prop bulk=3.91e6  shear=1.5e6  coh=20.0e3 fric=20 tension=0 range group clay1
prop bulk=4.44e6  shear=1.48e6 coh=5e3 fric=10 tension=0 range group clay2
prop bulk=16.67e6  shear=5.56e6 coh=10e3 fric=25 tension=0 range group clay3
;这里clay3中的coh你原来是10e10时候不对吧？我直接给你改为10e3
;pause
hist id=1  unbal
pl set back white
pl add con szz
pl add hist 1
pl show
solve elastic     ;这是最为简便求解初始地应力的方法，弹性求解法。与你的更改强度求解法的效果是一样。
;plot block group
;solve
save gp-3-130.sav

zsq

Sini dens 1800     ; 既然三个岩性的密度都一样，为什么还分别赋值呢
这个效果是一样的，没有关系，分开也可以

冬日的阳光

4# zsq
感谢啊 呵呵

冬日的阳光

3# jihezhilv
多谢指导

zhuxuekaihao

看了帖子，我知道边界条件原来这么重要，谢谢1

af2723

学习了。。。。。。。。。。。

pillar1990

学习学习。。

Gitar扶摇直上

如果修改后还不平衡又是什么原因呢？大胆猜测一下