初始地应力场的生成
;网格建立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 static/image/common/back.gif
彭文斌书上的例题就有的,至于为什么就不得而知了
他那个是在开挖之后设置的large 求初始地应力的时候不要用set large,这时候的变形可能达到几米,网格变形过大报错 wqq8889 发表于 2012-12-21 12:51 static/image/common/back.gif
求初始地应力的时候不要用set large,这时候的变形可能达到几米,网格变形过大报错 ...
可是如果不设置大变形,那么生成的初始地应力场会不会失真,有意义吗? 初始应力计算,不应该设置set large吧,最好是采用弹性计算
gongzhenyu 发表于 2012-12-24 20:28 static/image/common/back.gif
初始应力计算,不应该设置set large吧,最好是采用弹性计算
彭文斌书上的例题就有的,至于为什么就不得而知了 我的认为是弹性计算,可以初始化地应力,其目的就是得到地应力,其余的,可以忽略了。 如果是规则的网格,最好是用ini szz value grand ……来进行赋值
其实这个没有多大的作用,只是减少计算时间,快速收敛。
如果是不规则的网格,就直接用solve elas进行计算了。
这些计算,只是为了得到一个初始地应力的值,目的就是为了后面的施工工况服务的。 很明显呀,建模出错呀
页:
[1]