求助，初始应力怎么平衡不了！！ [复制链接]

大家给看看啊,初始应力怎么平衡不了！！谢谢了[/quote]！
gen zone b size 40 10 1  p0 0 0 0      p1 400 0 0       p2 0 200 0      p3 0 0 1.21 group 泥岩1
gen zone b size 40 10 2  p0 0 0 1.21   p1 400 0 1.21    p2 0 200 1.21   p3 0 0 7.07 group 灰岩1
gen zone b size 40 10 3  p0 0 0 7.07   p1 400 0 7.07    p2 0 200 7.07   p3 0 0 13.47 group 粉砂岩1
gen zone b size 40 10 1  p0 0 0 13.47  p1 400 0 13.47   p2 0 200 13.47  p3 0 0 14.96 group 灰岩2
gen zone b size 40 10 4  p0 0 0 14.96  p1 400 0  14.96  p2 0 200  14.96 p3 0 0 31.4  group 粉砂岩2
gen zone b size 40 10 1  p0 0 0 31.4   p1 400 0 31.4    p2 0 200 31.4   p3 0 0 32.46 group 泥岩2
gen zone b size 40 10 3  p0 0 0 32.46  p1 400 0 32.46   p2 0 200 32.46  p3 0 0 39.73 group 粉砂岩3
gen zone b size 40 10 6  p0 0 0 39.73  p1 400 0 39.73   p2 0 200 39.73  p3 0 0 51.93 group 细砂岩1
gen zone b size 40 10 7  p0 0 0 51.93  p1 400 0 51.93   p2 0 200 51.93  p3 0 0 65 group 粉砂岩4
gen zone b size 40 10 3  p0 0 0 65     p1 400 0 65      p2 0 200 65     p3 0 0 67.93 group 10煤
gen zone b size 40 10 4  p0 0 0 67.93  p1 400 0 67.93   p2 0 200 67.93  p3 0 0 71.97 group 泥岩3
gen zone b size 40 10 1  p0 0 0 71.97  p1 400 0 71.97   p2 0 200 71.97  p3 0 0 72.35 group 天然焦
gen zone b size 40 10 2  p0 0 0 72.35  p1 400 0 72.35   p2 0 200 72.35  p3 0 0 75.15 group 岩浆岩
gen zone b size 40 10 2  p0 0 0 75.15  p1 400 0 75.15   p2 0 200 75.15  p3 0 0 78.85 group 细砂岩2
gen zone b size 40 10 25 p0 0 0 78.85  p1 400 0 78.85   p2 0 200 78.85  p3 0 0 132.6 group 泥岩4
gen zone b size 40 10 4  p0 0 0 132.6  p1 400 0 132.6   p2 0 200 132.6  p3 0 0 136.79 group 细砂岩3
gen zone b size 40 10 3  p0 0 0 136.79 p1 400 0 136.79  p2 0 200 136.79 p3 0 0 140.24 group 泥岩5
gen zone b size 40 10 3  p0 0 0 140.24 p1 400 0 140.24  p2 0 200 140.24 p3 0 0 146.29 group 粉砂岩5
gen zone b size 40 10 4  p0 0 0 146.29 p1 400 0 146.29  p2 0 200 146.29 p3 0 0 150 group    泥岩6
m e
pro bu 1e10 she 1e10
in den 2.5e3
fix z range z 0
fix y range y 0
fix y range y 200
fix x range x 0
fix x range x 400
set grav 0 0 -10
ini szz -13.137e6 grad 0 0 2.5e4
in sxx -4.379e6 grad 0 0 0.83e4
in syy -4.379e6 grad 0 0 0.83e4  
apply szz -9.798e6 range z 150
so

看不出程序有什么问题啊，同等高手来解疑

你的各向同性材料模型还需要初始平衡吗？

new
generate zone brick p0 0 0 0 &
p1 30 0 0 &
p2 0 30 0  &
p3 0 0 10 &
p4 30 30 0 p5 0 30 10 &
p6 30 0 10 p7 30 30 10 &
size 30 30 10
gen zone brick p0 0 0 10  &
p1 20 0 10  p2 0 30 10 p3 0 0 40 p4 20 30 10 &
p5 0 30 40  p6 10 0 40 p7 10 30 40 &
size 20 30 30
plot add surface yellow
plot add axes
plot show
model mohr
property bulk 1e8 shear 0.3e8 friction 35 &
cohesion 1e10 tension 1e10
set gravity 0 0 -9.81
initial density 2000
fix  range z 0
fix y range y 0
fix y range y 30
history gp zdisp 12 15 6
set mechanical force 100
solve
我的也是啊 求教

路过学习的

设置重力后先求解