初始地应力场的生成

jsaturday

;网格建立
gen zone brick p0 3 0 -9 p1 16.5 0 -9 p2 3 1 -9 p3 3 0 0 size 27 2 18 rat 1 1 1
gen zone brick p0 16.5 0 -9 p1 25 0 -9 p2 16.5 1 -9 p3 16.5 0 0 size 10 2 18 rat 1 1 1

gen zone brick p0 0 0 -14 p1 3 0 -14 p2 0 1 -14 p3 0 0 -9 size 6 2 10 rat 1 1 1
gen zone brick p0 3 0 -14 p1 16.5 0 -14 p2 3 1 -14 p3 3 0 -9 size 27 2 10 rat 1 1 1
gen zone brick p0 16.5 0 -14 p1 25 0 -14 p2 16.5 1 -14 p3 16.5 0 -9 size 10 2 10 rat 1 1 1
gen zone reflect ori 0 0 0 norm 1 0 0

;生成初始应力场

model mohr
prop bulk 500e6 sh 330e6 coh 2e3 fri 25
fix x y z ran z -14.1 -13.9
fix x ran x 24.9 25.1
fix x ran x -25.1 -24.9
fix y ran y 0.9 1.1
fix y ran y -1.1 -0.9

set grav 0 0 -10
ini dens 2000
set large
hist unbal
set mech ratio 1e-4
solve elas
最后显示illegal geometry in zone 4223 是模型失稳还是网格建立的过程中有问题，请过来人帮忙分析一下，谢谢

南湾湖

jsaturday 发表于 2012-12-24 20:29
彭文斌书上的例题就有的，至于为什么就不得而知了

他那个是在开挖之后设置的large

wqq8889

求初始地应力的时候不要用set large，这时候的变形可能达到几米，网格变形过大报错

jsaturday

wqq8889 发表于 2012-12-21 12:51
求初始地应力的时候不要用set large，这时候的变形可能达到几米，网格变形过大报错 ...

可是如果不设置大变形，那么生成的初始地应力场会不会失真，有意义吗？

gongzhenyu

初始应力计算，不应该设置set large吧，最好是采用弹性计算

jsaturday

gongzhenyu 发表于 2012-12-24 20:28
初始应力计算，不应该设置set large吧，最好是采用弹性计算

彭文斌书上的例题就有的，至于为什么就不得而知了

gongzhenyu

我的认为是弹性计算，可以初始化地应力，其目的就是得到地应力，其余的，可以忽略了。

gongzhenyu

如果是规则的网格，最好是用ini szz value grand ……来进行赋值
其实这个没有多大的作用，只是减少计算时间，快速收敛。
如果是不规则的网格，就直接用solve elas进行计算了。
这些计算，只是为了得到一个初始地应力的值，目的就是为了后面的施工工况服务的。

芒草尖之舞

很明显呀，建模出错呀