初始应力不平衡

yan0416

大家给看看啊,初始应力怎么平衡不了！！谢谢了[/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

yan0416

高手指点下啊？感激不尽！！！

wangqingniu

差不多了吧，我平时也就这样直接开挖了

sksky

4# 海底冬天

这不是初始化应力吗，szz是自重引起的，sxx和syy是水平应力

yan0416

4# 海底冬天

水平方向的应力，根据测压系数确定的。n=u/1-u，u泊松比。水平应力=自重应力*n

yan0416

为什么y方向会隆起???????搞不明白!!!!

silision

ini szz -13.137e6 grad 0 0 2.5e4
ini sxx -4.379e6 grad 0 0 0.83e4
ini syy -4.379e6 grad 0 0 0.83e4  
apply szz -9.798e6 range z 150
中，ini szz 和apply szz分别代表什么意思？我觉得重复了啊。。
你这样算的收敛很快吧？