请教自重应力平衡后应力云图不规律的问题，诚挚感谢

cumt安全

命令流很简要，如下所示：
new
plo blo gro
gen zon bri p0 0 0 0  p1 40 0 0  p2 0 20 0  p3 0 0 1  size 80 40 2  gro 1
gen zon bri p0 0 0 1  p1 40 0 1  p2 0 20 1  p3 0 0 10 size 80 40 18 gro 2
gen zon bri p0 0 0 10 p1 40 0 10 p2 0 20 10 p3 0 0 16 size 80 40 12 gro 3
gen zon bri p0 0 0 16 p1 40 0 16 p2 0 20 16 p3 0 0 29 size 80 40 26 gro 4
gen zon bri p0 0 0 29 p1 40 0 29 p2 0 20 29 p3 0 0 55 size 80 40 52 gro 5


model mohr


prop bulk 5.2e9  shear 3.2e9   coh 6.3e6  fric 40   tens 1.5e6    ran z 0 1
prop bulk 2.8e9  shear 1.3e9   coh 2.5e6  fric 32   tens 1.0e6    ran z 1 10
prop bulk 1.6e9  shear 6.0e8   coh 8.0e5  fric 24   tens 7.0e5    ran z 10 16
prop bulk 2.8e9  shear 1.3e9   coh 2.5e6  fric 32   tens 1.0e6    ran z  16 29
prop bulk 5.2e9  shear 3.2e9   coh 6.3e6  fric 40   tens 1.5e6    ran z  29 55


ini dens 2500 ran z 0 1
ini dens 2300 ran z 1 10
ini dens 1400 ran z 10 16
ini dens 2300 ran z 16 29
ini dens 2500 ran z 29 55


fix x ran x -0.1 0.1
fix x ran x 39.9 40.1
fix y ran y -0.1 0.1
fix y ran y 19.9 20.1
fix z ran z -0.1 0.1


set grav 0 0 -9.8
app nstress -10e6 ran x 0 40 y 0 20 z 55
hist unbal
solve
save 自重应力平衡计算结果.sav
plot con szz ou on
plot add hist 1


rest 自重应力平衡计算结果.sav
ini xvel 0 yvel 0 zvel 0
initial xdisp = 0  ydisp = 0  zdisp = 0


;巷道开挖
model null ran x 11 15 y 0 20 z 11 15
plot block group
solve
save 巷道开挖.sav
plot con  zdis ou on
plot con  szz  ou on


自重应力平衡后 应力云图如下:

cumt安全

新人，求点拨 thx

zhaijun

自重应力平衡计算时，应该分别进行弹塑性分析吧，云图的原因不知道是不是这个原因造成的，你试着算一下看看

cumt安全

zhaijun 发表于 2015-3-24 14:49
自重应力平衡计算时，应该分别进行弹塑性分析吧，云图的原因不知道是不是这个原因造成的，你试着算一下看看 ...

好的，我试试，谢谢

caikuangren

即使把粘聚力和抗拉强度设置的很大或者用solve elastic求解，还会存在类似问题。岩层承压会形变而在边界处又限制了这些变形，故在边界处局部应力异常，是正常现象，不影响你后续开挖的计算。

cumt安全

caikuangren 发表于 2015-3-24 16:40
即使把粘聚力和抗拉强度设置的很大或者用solve elastic求解，还会存在类似问题。岩层承压会形变而在边界处 ...

好的，谢谢师兄，我是安全学院研一的，还请多多指教

zhouff1990

自重平衡采用弹性模型计算试试，将刚度调大。平衡后再改成摩尔库伦模型计算。

cumt安全

zhouff1990 发表于 2015-3-28 09:04
自重平衡采用弹性模型计算试试，将刚度调大。平衡后再改成摩尔库伦模型计算。 ...

好的，谢谢