aibbin 发表于 2009-8-5 16:07:22

边坡三维建模!怎么计算都不出结果?

哪位兄弟帮忙看看,我是从书上找的例子,书上的例子能计算(边坡取单位宽度,相当于二维)。我把边坡宽度加大了(三维的),可是就是没有结果,帮忙看看哪里出毛病了,谢谢了!

gen zone brick p0 0 0 0 p1 20 0 0 p2 0 50 0 p3 0 0 10 size 10 10 10 ratio 0.97 1 1 group ceng11
gen zone brick p0 20 0 0 p1 30 0 0 p2 20 50 0 p3 20 0 10 size 3 10 10 group ceng12
gen zone brick p0 0 0 10 p1 20 0 10 p2 0 50 10 p3 0 0 20 p4 20 50 10 &
               p5 0 50 20 p6 15 0 18 p7 15 50 18 size 10 10 10 ratio 0.97 1 1 group ceng2
gen zone brick p0 0 0 20 p1 15 0 18 p2 0 50 20 p3 0 0 28 p4 15 50 18 &
               p5 0 50 28 p6 10 0 24 p7 10 50 24 size 10 10 10 ratio 0.97 1 1 group ceng3
fix x y z range z -0.1 0.1
fix xrange x -0.1 0.1
fix xrange x 29.9 30.1
fix y
model elas
prop density 2000 bulk 3e9 shear 1e9 range group ceng11
prop density 2000 bulk 3e9 shear 1e9 range group ceng12
prop density 2000 bulk 3e9 shear 1e9 range group ceng2
prop density 2000 bulk 3e9 shear 1e9 range group ceng3
set gravity 0 0 -10.0
solve
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
model mohr
prop dila 20 bulk 1.0e8 shear 3.0e7 coh 12380 tens 1.0e6 fric 20 range group ceng11
prop dila 20 bulk 1.0e8 shear 3.0e7 coh 12380 tens 1.0e6 fric 20 range group ceng12
prop dila 21 bulk 1.1e8 shear 2.9e7 coh 11111 tens 1.2e6 fric 21 range group ceng2
prop dila 19 bulk 1.1e8 shear 3.1e7 coh 12222 tens 1.1e6 fric 18 range group ceng3
ini density 2000 range group ceng11
ini density 2000 range group ceng12
ini density 2100 range group ceng2
ini density 2200 range group ceng3
solve fos file slope3Dfos.sav associated

haiqsim 发表于 2009-8-8 18:05:25

我觉得你的既然是三维计算,边界条件还用“fix y”的话,是否不大合适啊?
人家本来是按二维的算的,所以用“fix y”,你这个三维的,相当于也把所有的网格节点都限制y方移动,个人感觉不大妥当。
改为约束两侧吧。

再者,算这个安全系数,是非常非常非常非常慢的,你得等着。等4、5个小时也是可能的啊。
页: [1]
查看完整版本: 边坡三维建模!怎么计算都不出结果?