jsaturday 发表于 2012-12-20 21:42:55

初始地应力场的生成

;网格建立
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 是模型失稳还是网格建立的过程中有问题,请过来人帮忙分析一下,谢谢

南湾湖 发表于 2013-5-28 21:22:12

jsaturday 发表于 2012-12-24 20:29 static/image/common/back.gif
彭文斌书上的例题就有的,至于为什么就不得而知了

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

wqq8889 发表于 2012-12-21 12:51:56

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

jsaturday 发表于 2012-12-21 15:32:19

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

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

gongzhenyu 发表于 2012-12-24 20:28:07

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

jsaturday 发表于 2012-12-24 20:29:41

gongzhenyu 发表于 2012-12-24 20:28 static/image/common/back.gif
初始应力计算,不应该设置set large吧,最好是采用弹性计算

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

gongzhenyu 发表于 2012-12-25 14:53:38

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

gongzhenyu 发表于 2012-12-26 11:26:32

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

芒草尖之舞 发表于 2012-12-28 19:29:16

很明显呀,建模出错呀
页: [1]
查看完整版本: 初始地应力场的生成