FLAC深埋工程初始地应力场的S-B法生成的改进失败 [复制链接]

深埋工程初始地应力场的S-B法生成的改进   （FLAC/flac3d基础与工程实例 陈育明编著）p101页

原来的命令
new
gen zone brick p0 0 0 0 p1 60 0 0 p2 0 60 0 p3 0 0 90 &
p4 60 60 0 p5 0 60 90 p6 60 0 150 p7 60 60 150 &
size 6 6 10
model elas
pro bulk 10e10 she 10e10
ini den 2500
apply sxx -1e9 grad 0 0 1.1111111e7 range x -.1 .1
apply sxx -1e9 grad 0 0 7.6666666e6 range x 59.9 60.1
apply syy -1e9 grad 0 0 8.3333333e6 range y -.1 .1
apply syy -1e9 grad 0 0 8.3333333e6 range y 59.9 60.1
apply szz -1e8 grad 0 0 8.3333333e5 range z 0 120
set grav 0 0 -10
step 30000
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
plot cont szz

想把这个过程用动画演示出来，但是没有成功，希望高人指点，小弟在这里感激不尽！！！！

new
gen zone brick p0 0 0 0 p1 60 0 0 p2 0 60 0 p3 0 0 90 &
p4 60 60 0 p5 0 60 90 p6 60 0 150 p7 60 60 150 &
size 6 6 10
model elas
pro bulk 10e10 she 10e10
ini den 2500
hist unbal
apply sxx -1e9 grad 0 0 1.1111111e7 range x -.1 .1
apply sxx -1e9 grad 0 0 7.6666666e6 range x 59.9 60.1
apply syy -1e9 grad 0 0 8.3333333e6 range y -.1 .1
apply syy -1e9 grad 0 0 8.3333333e6 range y 59.9 60.1
apply szz -1e8 grad 0 0 8.3333333e5 range z 0 120
set grav 0 0 -10
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
plot set rot 10 0 30
plot con szz ou on
plot add hist 1
set movie avi file 7-7.avi
movie start
step 30000
solve

怎么没有人回答呢？

参数怎么确定？如1.1111等

还是没有人回应·····

失败是怎么个失败法，建议楼主把问题说清楚，最好截图

可以到ITASCA官方论坛，找技术工程师帮您看看；
地址：http://bbs.itasca.cn

看不懂！

我实验了一下,你可以看一下计算过程,到最后时Meach ratio 在1.3e-2左右而系统默认的是1e-5所以会无限计算下去,你应该把最后一行的SOLVE去掉.

非常感谢，我试试

请问例子中的一些参数是如何得到的
apply sxx -1e9 grad 0 0 1.1111111e7 range x -.1 .1
apply sxx -1e9 grad 0 0 7.6666666e6 range x 59.9 60.1
apply syy -1e9 grad 0 0 8.3333333e6 range y -.1 .1
apply syy -1e9 grad 0 0 8.3333333e6 range y 59.9 60.1
apply szz -1e8 grad 0 0 8.3333333e5 range z 0 120
中的1.11111111e7 ,7.6666666e6和8.333333333333e5

还有最后apply szz -1e8 grad 0 0 8.3333333e5 range z 0 120，为什么是range z 0 120 ,不是和之前一样取两个面呢

那如果此时的不平衡力很大，也可以算初始应力场么