初始应力计算不收敛
请高手指点,谢谢!new
gen zone brick p0 0 0 0 p1 372 0 0 p2 0 200 0 p3 0 0 8 size 62 20 2 group zsg1;黏土层
gen zone brick p0 0 0 8 p1 372 0 8 p2 0 100 8 p3 0 0 50 size 62 10 7 group zsg2 ;细砂土层
gen zone brick p0 0 100 8 p1 372 100 8 p2 0 200 8 p3 0 100 50 size 62 10 7 group zsg3 ;砂壤土层
gen zone brick p0 200 0 50 p1 236 0 50 p2 200 200 50 p3 230 0 60 p4 236 200 50 &
p5 230 200 60 p6 236 0 60 p7 236 200 60 size 6 20 5 group zsg4;堤身壤土层
gen zone brick p0 236 0 50 p1 272 0 50 p2 236 200 50 p3 236 0 60 p4 272 200 50 &
p5 236 200 60 p6 242 0 60 p7 242 200 60 size 6 20 5 group zsg5;堤身壤土层
gen zone brick p0 100 0 50 p1 200 0 50 p2 100 200 50 p3 124 0 58 size 25 20 4 group zsg6;淤背砂土
attach face rang x 100 200 z 49.9 50.1
plot
add surface white
show
;pause
model mohr range group zsg1
prop den=1980.0 bulk=2.29e6 shear=9.94e5 coh=92650.0 tens=853 friction=24.15 dilation=20
model mohr range group zsg2
prop den=1900.0 bulk=2.09e8 shear=1.25e8 coh=0.0 tens=0.0 friction=34.77 dilation=20
model mohr range group zsg3
prop den=1900.0 bulk=2.19e7 shear=1.20e7 coh=1627.0 tens=520.0 friction=31.73 dilation=20
model mohr range group zsg4
prop den=1710.0 bulk=5.09e6 shear=1.30e6 coh=27000.0 tens=853.0 friction=30.74 dilation=20
model mohr range group zsg5
prop den=1710.0 bulk=5.09e6 shear=1.30e6 coh=27000.0 tens=853.0 friction=30.74 dilation=20
model null range group zsg6
set gravity 0,0,-9.81
fix x y z range z -0.1 0.1
fix x range x 371.9 372.1
fix x range x -0.1 0.1
fix y range y -0.1 0.1
fix y range y 199.9 200.1
step 5000
;set mechanical force 50
;solve
;initial xdisplacement=0 ydisplacement=0 zdisplacement=0
;model mohr range group zsg6
;prop den=1880.0 bulk=2.0e8 shear=3.0e7 coh=733.0 tens=1.0e2 friction=28.0 dilation=20
;step 200
;set mechanical force 50
;solve
plot create stress1
plot add contour smax plane behind shade on
plot add axes
plot show
plot create stress2
plot add contour smin plane behind shade on
plot add axes
plot show
plot create stress3
plot add contour smid plane behind shade on
plot add axes
plot show
save cstress.sav
set log on
pr gp pos
pr zone gp
pr zone str
set log off 我感觉你的参数问题很大啊，怎么粘土层的c值怎么那么大啊；另外细沙图层的变形模量也很大。
最好按照flac的标准去写命令流，比如 ini dens 我做过的粘土是63900 我试了一下 你的几何模型有问题
原帖由 zsg6537 于 2008-5-14 08:12 发表 http://forum.simwe.com/images/common/back.gif
请高手指点,谢谢!
new
gen zone brick p0 0 0 0 p1 372 0 0 p2 0 200 0 p3 0 0 8 size 62 20 2 group zsg1;黏土层
gen zone brick p0 0 0 8 p1 372 0 8 p2 0 100 8 p3 0 0 50 size 62 10 7 group zsg2 ;细砂土层 ... 我真无聊，把你的问题修改好了
原帖由 chaoshawnli 于 2008-12-5 16:47 发表 http://forum.simwe.com/images/common/back.gif
我试了一下 你的几何模型有问题
好好学学 谢谢版主 谢谢版主 谢谢 谢谢版主谢谢 几何模型没有问题，关键是你计算时没有用solve elastic，而是用的step5000
页:
[1]