flac地应力平衡完地面有隆起怎么回事啊？ [复制链接]

flac地应力平衡完地面有隆起怎么回事啊？这个高了很久了都不行啊大神给分析下吧

gen zone cylinder p0(0 0 0) p1(2.45 0 0) p2(0 40 0) p3(0 0 2.45) &
size 5 8 6  ratio 0.7 1 1 group soil

gen zone cshell p0(0 45 0) p1(2.5 45 0) p2(0 85 0) &
p3(0 45 2.5) p4(2.5 85 0) p5(0 85 2.5) &
p8(2.45 45 0) p9 (0 45 2.45) p10(2.45 85 0) &
p11(0 85 2.45) &
size 1 16 6 5 dimension 4 4 4 4 &
ratio = 1 1 1 0.7   group soilout

gen zone radcyl p0(0 90 0) p1(40 90 0) p2(0 130 0) p3(0 90 40) &
     p4(40 130 0) p5(0 130 40) p6(40 90 40) p7(40 130 40) &
   p8(2.5 90 0) p9(0 90 2.5) p10(2.5 130 0) p11(0 130 2.5) &
size 5 16 6 30 dimension 4 4 4 4 ratio = 0.7 1 1 1.15  group soils

gen zone radcyl p0(0 -70 0) p1(40 -70 0) p2(0 -10 0) p3(0 -70 40) &
p4(40 -10 0) p5(0 -10 40) p6(40 -70 40) p7(40 -10 40) &
p8(2.5 -70 0) p9(0 -70 2.5) p10(2.5 -10 0) p11(0 -10 2.5) &
size 5 24 6 30 dimension 4 4 4 4 ratio = 0.7 1 1 1.15 group soildown

gen zone cylinder p0(0 -140 0) p1(2.45 -140 0) p2(0 -80 0) p3(0 -140 2.45) &
size 5 24 6  ratio 0.7 1 1 group soild

gen zone cshell p0(0 -210 0) p1(2.5 -210 0) p2(0 -150 0) &
p3(0 -210 2.5) p4(2.5 -150 0) p5(0 -150 2.5) &
p8(2.45 -210 0) p9 (0 -210 2.45) p10(2.45 -150 0) &
p11(0 -150 2.45) &
size 1 24 6 5 dimension 4 4 4 4 &
ratio = 1 1 1 0.7   group soilp

gen zone reflect normal(0 0 -1) origin(2.5 50 0)
gen zone reflect normal(1 0 0) origin(0 0 0)  

ini y add -90 range group soils
ini y add -45 range group soilout
ini y add 10 range group soildown
ini y add 80 range group soild
ini y add 150 range group soilp

gen merge 1e-3
;attach face range y 0.1 -0.1


model mohr                        
prop bulk 1.666e8 shear 7.693e7 coh 20000 fric 12 di 0 ten 1e10  
ini dens 800
plot add con yd ou on
plot add axes


fix y range y -59.9 -60.1
fix x range y -59.9 -60.1
fix z range y -59.9 -60.1

fix x range x -39.9 -40.1
fix x range x 39.9 40.1

fix z range z -39.9 -40.1
fix z range z 39.9 40.1

set grav 0 -10 0

ini syy 0 grad 0 8000 0 range group y -60 40
ini szz 0 grad 0 3429 0 range group y -60 40
ini sxx 0 grad 0 3429 0 range group y -60 40
;plot con ydis
hist nstep = 5
hist unbal
hist gp zdisp 0 0 0
pause
;step 20000
so
save 11.sav

那根本不是隆起. 有微弱差异是模型的网格造成的.
而且, 正常模拟顺序应该是先做弹性平衡, 再开挖, 再计算你的工程.

设置参数的问题把，个人认为