初始地应力场无法收敛
大家看看我这命令流有什么问题 为什么计算就是收敛不了 顶起来 没人帮忙吗 本帖最后由 jihezhilv 于 2010-8-22 12:27 编辑抽空给你看了一下：问题还真不少，呵呵
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 -25p5 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.91e6shear=1.5e6coh=20.0e3 fric=20 tension=0 range group clay1
prop bulk=4.44e6shear=1.48e6 coh=5e3 fric=10 tension=0 range group clay2
prop bulk=16.67e6shear=5.56e6 coh=10e3 fric=25 tension=0 range group clay3
;这里clay3中的coh你原来是10e10时候不对吧？我直接给你改为10e3
;pause
hist id=1unbal
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 Sini dens 1800 ; 既然三个岩性的密度都一样，为什么还分别赋值呢
这个效果是一样的，没有关系，分开也可以 4# zsq
感谢啊 呵呵 3# jihezhilv
多谢指导 看了帖子，我知道边界条件原来这么重要，谢谢1 学习了。。。。。。。。。。。 学习学习。。 如果修改后还不平衡又是什么原因呢？大胆猜测一下:$
页:
[1]