做过基坑的人都来看看啊，地表沉降的问题

mthinking

基坑开挖以后怎么再远离基坑一段距离以后第二次出现沉降较大的位置啊，按说是离基坑越远沉降就慢慢减小的，不会再出现“反弹”啊。。。请各位高手看看

mthinking

这是初始文件:
new
gen zone brick p0 0 0 -12 p1 1.5 0 -12 p2 0 24 -12 p3 0 0 0 size 5 35 20 ratio 1 1.05 1   
gen zone brick p0 1.5 0 -12 p1 0 0 -12 p2 1.5 -8 -12 p3 1.5 0 0 size 5 16 20 ratio 1 1.1 1   

gen zone brick p0 0 0 -16 p1 1.5 0 -16 p2 0 24 -16 p3 0 0 -12 size 5 35 5 ratio 1 1.05 1   
gen zone brick p0 1.5 0 -16 p1 0 0 -16 p2 1.5 -8 -16 p3 1.5 0 -12 size 5 16 5 ratio 1 1.1 1

group s1_1 range z -2 0
group s2_2 range z -2 -5
group s3_3 range z -9 -5
group s4_4 range z -16 -9
range name s1_1 group s1_1
range name s2_2 group s2_2
range name s3_3 group s3_3
range name s4_4 group s4_4

model mohr

prop bulk 2.824e7 shear 1.155e7 coh 13000 tens 1500 fric 18.5 range s1_1
prop bulk 3.000e7 shear 6.429e6 coh 12000 tens 3000 fric 12.0 range s2_2
prop bulk 1.333e7 shear 4.444e6 coh 15000 tens 2000 fric 13.5 range s3_3
prop bulk 4.333e7 shear 2.000e7 coh 18800 tens 2500 fric 15.9 range s4_4


ini dens 1800 range s1_1
ini dens 1750 range s2_2
ini dens 1700 range s3_3
ini dens 1810 range s4_4
set grav 0 0 -10

fix x y z range z -16.1 -15.9
fix x range x -0.1 0.1
fix x range x 1.4 1.6
fix y range y -8.1 -7.9
fix y range y 23.9 24.1

ini szz 0 grad 0 0 18000 range z 0 -2
ini sxx 0 grad 0 0 9692.3 range z 0 -2
ini syy 0 grad 0 0 9692.3 range z 0 -2

ini szz -1000   grad 0 0 17500    range z -2 -5
ini sxx 3948.74 grad 0 0 11666.67 range z -2 -5
ini syy 3948.74 grad 0 0 11666.67 range z -2 -5

ini szz -3500    grad 0 0 17000   range z -5 -9
ini sxx -8615.37 grad 0 0 9153.85 range z -5 -9
ini syy -8615.37 grad 0 0 9153.85 range z -5 -9

ini szz 6400      grad 0 0 18100   range z -9 -16
ini sxx -21185.74 grad 0 0 7757.14 range z -9 -16
ini syy -21185.74 grad 0 0 7757.14 range z -9 -16

plot add con szz ou on  
his unbal
his gp zd 1 0 1
so
pause
ini state 0
ini xd 0 yd 0 zd 0
ini xvel 0 yvel 0 zvel 0

group excav_1 range y -8 0 z -2 0
group excav_2 range y -8 0 z -3.5 -2
group excav_3 range y -8 0 z -5.0 -3.5
group excav_4 range y -8 0 z -6.5 -5.0
group excav_5 range y -8 0 z -8.0 -6.5


range name excav_1 group excav_1
range name excav_2 group excav_2
range name excav_3 group excav_3
range name excav_4 group excav_4
range name excav_5 group excav_5

sav gen_ini+tens.sav

mthinking

没有人吗
抓狂了。。。

eagleming

图形上最好能贴上矢量方向，你那个颜色挺难区分的

haiqsim

LZ算了2w多步啊，怎么算的啊，每开挖一次都solve么？

haiqsim

我是这样试算了下，每开挖一次step 200，
最后又step 500，这样得到如下的结果，看来与每次都solve的结果相差很大啊

mthinking

6# haiqsim w w
我每次开挖都要slove至平衡的。
之前生成趋势应力时求20000步，是因为我拍自重固结不充分而影响后面的位移。但这样做了没有什么效果，那看来就不是这个原因了。

haiqsim

不知LZ如何每次开挖都solve至默认平衡，
我开挖到第二步solve时，位移都达到3m了，
mech.ratio还是不能小于1e-5，偶很奇怪。

mthinking

我是加了支护的，这上面的文件只是初始应力平衡 8# haiqsim

mthinking

这是基坑后缘加长的40m后竖向位移图