24小时热门版块排行榜    

查看: 9221  |  回复: 1

晓浪哥

新虫 (初入文坛)

[求助] 求助flac3d不收敛问题 已有1人参与

新手金币少,有劳大神了。我想做俯斜式挡土墙的模拟,在初始地应力计算就不收敛了,监测的最大不平衡力呈一直循环波动状态。
下面是我的代码和计算的截图:
new
gen zone brick &
p0 0 0 0 p1 8 0 0 p2 0 10 0 p3 0 0 1 size 8 10 2 group 岩石地基
gen zone brick &
p0 8 0 0 p1 11.46 0 0 p2 8 10 0 p3 8 0 1 size 4 10 2 group 岩石地基
gen zone brick &
p0 11.46 0 0 p1 30 0 0 p2 11.46 10 0 p3 11.46 0 1 size 20 10 2 group 岩石地基
gen zone brick &
p0 8 0 1 p1 11.46 0 1 p2 8 10 1 p3 9 0 5  p4 11.46 10 1 p5 9 10 5 p6 10 0 5  p7 10 10 5 size 4 20 8 group 挡土墙
gen zone brick &
p0 11.46 0 1 p1 30 0 1 p2 11.46 10 1  p3 11.46 0 5 size 20 20 8  group 墙后填土
gen zone wedge &
p0 11.46 0 5 p1 10 0 5 p2 11.46 10 5  p3 11.46 0 1 p4 10 10 5  p5 11.46 10 1 size 2 20 8 group 墙后填土

;添加接触面
gen separate 墙后填土
interface 1 wrap 墙后填土 挡土墙
interface 1 prop kn 4e8 ks 4e8 coh 0.5e3 fric 18

plot show
;静力计算
config dyn
set dyn off
model elas
set grav 0, 0, -9.8
pro bulk 4e10 she 2.4e10 dens 2700 range group 岩石地基
pro bulk 4e10 she 2.4e10 dens 2700 range group 挡土墙
pro bulk 4e10 she 2.4e10 dens 2700 range group 墙后填土
hist unbal
;set mech force 1.09e4
solve (这个solve就不能收敛了)
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
model elas range group 墙后填土
pro bulk 0.514e8 shear 0.231e7 coh 0.134e5 fric 30 ten 1e10 dens 1810
model elas  range group 岩石地基
pro bulk 0.18e11 shear 0.125e11 dens 2400
model elas  range group 挡土墙
pro bulk 0.11e11 shear 0.11e11 dens 2300
fix x range x -0.1 0.1
fix x range x 29.9 30.1
fix y range y -0.1 0.1
fix y range y 9.9 10.1
fix x y z range z -0.1 0.1
;hist unbal
;set large
;set mech ratio 1.07e-3
solve
save 静力阶段

求助flac3d不收敛问题
`VL5Z)LA0E}VP{1$V5AZJQ4.png
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

哒哒已注册

禁虫 (初入文坛)

本帖内容被屏蔽

2楼2018-06-27 14:52:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 晓浪哥 的主题更新
信息提示
请填处理意见