为什么边界附近应力场往上翘呢

zhuxuekaihao

做了一个简单的模型，代码如下：
new
ge zone brick size 25 25 5 p0=(0,0,0) &
   p1=(120,0,0) p2=(0,100,0) p3=(0,0,10) ratio 1 1 1 group 1 ;最底部砂岩
ge zone brick size 50 50 10 p0=(0,0,10) &
   p1=(120,0,10) p2=(0,100,10) p3=(0,0,24.5) ratio 1 1 1 group 2 ;放顶煤煤层
ge zone brick size 25 25 2 p0=(0,0,24.5) &
   p1=(120,0,24.5) p2=(0,100,24.5) p3=(0,0,27.5) ratio 1 1 1 group 3 ;泥岩
ge zone brick size 25 25 4 p0=(0,0,27.5) &
   p1=(120,0,27.5) p2=(0,100,27.5) p3=(0,0,35.5) ratio 1 1 1 group 4 ;砂岩
set gravity 0 0 -10
fix x y z range z -0.1 0.1
fix x range x -0.1 0.1
fix x range x 119.9 120.1
fix y range y -0.1 0.1
fix y range y 99.9 100.1
apply szz -1e6 range z 35.4 35.6plot con szz out on magf 1
solve
得到的主应力图如下：

应力分布的数值是正确的，可是为什么除坐标系所在的角之外，其余3角都往上翘啊，难道是边界的问题？

ld1339427940

边界fix应该没有问题吧 貌似没有attach face？

zhuxuekaihao

ld1339427940 发表于 2013-8-20 14:49
边界fix应该没有问题吧 貌似没有attach face？

嗯，那我试试attach吧

zhuxuekaihao

ld1339427940 发表于 2013-8-20 14:49
边界fix应该没有问题吧 貌似没有attach face？

还是不行啊

yangxin546

你用塑性模型算自重平衡就是这样的吧，难道你否认用的不是塑性模型？

zhuxuekaihao

yangxin546 发表于 2013-8-25 16:39
你用塑性模型算自重平衡就是这样的吧，难道你否认用的不是塑性模型？

什么是塑形模型啊？您说的是sovel elastic？