新手金币少,有劳大神了。我想做俯斜式挡土墙的模拟,在初始地应力计算就不收敛了,监测的最大不平衡力呈一直循环波动状态。
下面是我的代码和计算的截图:
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 静力阶段
|
|