一个圆形基坑开挖实例

xiangruide

我做了一个圆形基坑开挖的例子，算出来后水平方向的位移应该很正常，但是竖直方向的位移始终都是向上的，基坑底隆起正常但周围也是隆起感觉有点不对，请大家帮看看问题出在哪儿。
       基坑开挖的半径是15m。深度为30m，y方向设为竖直方向，x、z方向为水平方向。我把原文件和三个方向的位移图都贴出。
new
;简单试验模型
gen zone radcyl p0 0 0 0 p1 30 0 0 p2 0 60 0 p3 0 0 30 dim 16 16 16 16 size 3 30 12 7 group 原岩
gen zone cshell p0 0 0 0 p1 16 0 0 p2 0 60 0 p3 0 0 16 size 2 30 12 3 dim 15 15 15 15 group 连续墙1
gen zone cshell p0 0 0 0 p1 15 0 0 p2 0 60 0 p3 0 0 15 size 2 30 12 7 dim 14 14 14 14 group 连续墙2 fill group 开挖部分
gen zone brick  p0 30 0 0 p1 45 0 0 p2 30 60 0 p3 30 0 30 size 3 30 6 group 原岩
gen zone brick  p0 0 0 30 p1 30 0 30 p2 0 60 30 p3 0 0 45 size 6 30 3 group 原岩
gen zone brick  p0 30 0 30 p1 45 0 30 p2 30 60 30 p3 30 0 45 size 3 30 3 group 原岩
gen zone reflect normal -1 0 0 origin 0 0 0
gen zone reflect normal 0 0 -1 origin 0 0 0
plot grid
plot add sur yellow
model mohr
def derive1
s_mod1=y_mod1/(2.0*(1.0+p_ratio1))
b_mod1=y_mod1/(3.0*(1.0-2.0*p_ratio1))
end
set y_mod1=10e6 p_ratio1=0.3
derive1
def derive2
s_mod2=y_mod2/(2.0*(1.0+p_ratio2))
b_mod2=y_mod2/(3.0*(1.0-2.0*p_ratio2))
end
set y_mod2=11.9e6 p_ratio2=0.3  
derive2
def derive3
s_mod3=y_mod3/(2.0*(1.0+p_ratio3))
b_mod3=y_mod3/(3.0*(1.0-2.0*p_ratio3))
end
set y_mod3=50e6 p_ratio3=0.27  
derive3
def derive4
s_mod4=y_mod4/(2.0*(1.0+p_ratio4))
b_mod4=y_mod4/(3.0*(1.0-2.0*p_ratio4))
end
set y_mod4=13.1e6 p_ratio4=0.29  
derive4
prop bulk b_mod1  shear s_mod1  coh 30000  fric 10 range y 48 60
prop bulk b_mod2  shear s_mod2  coh 40000  fric 25 range y 38 48
prop bulk b_mod3  shear s_mod3  coh 45000  fric 35 range y 32 38
prop bulk b_mod4  shear s_mod4  coh 40000  fric 28 range y 0 32
ini dens 1930 range y 48 60
ini dens 2020 range y 38 48
ini dens 2100 range y 32 38
ini dens 2030 range y 0 32
set gravity 0 -10 0
fix x range x -45.2 -44.8
fix x range x 44.8 45.2
fix z range z -45.2 -44.8
fix z range z 44.8 45.2
fix x y z range y -0.1 0.1
solve
save gra.sav
;开挖
ini xdisp 0 ydisp 0 zdisp 0
hist gp zdisp 18 60 0
m n range group 开挖部分 y 56 60
m n range group 连续墙2 y 56 60
plot con ydis
plot add hist 1
step 300
m n range group 开挖部分 y 52 56
m n range group 连续墙2 y 52 56
m e range group 连续墙2 y 56 60
prop bulk 34.5e9 shear 14.78e9 range group 连续墙2 y 56 60
step 300
m n range group 开挖部分 y 48 52
m n range group 连续墙2 y 48 52
m e range group 连续墙2 y 52 56
prop bulk 34.5e9 shear 14.78e9 range group 连续墙2 y 52 56
step 300
m n range group 开挖部分 y 44 48
m n range group 连续墙2 y 44 48
m e range group 连续墙2 y 48 52
prop bulk 34.5e9 shear 14.78e9 range group 连续墙2 y 48 52
step 300
m n range group 开挖部分 y 40 44
m n range group 连续墙2 y 40 44
m e range group 连续墙2 y 44 48
prop bulk 34.5e9 shear 14.78e9 range group 连续墙2 y 44 48
step 300
m n range group 开挖部分 y 36 40
m n range group 连续墙2 y 36 40
m e range group 连续墙2 y 40 44
prop bulk 34.5e9 shear 14.78e9 range group 连续墙2 y 40 44
step 300
m n range group 开挖部分 y 30 36
m n range group 连续墙2 y 30 36
m e range group 连续墙2 y 36 40
prop bulk 34.5e9 shear 14.78e9 range group 连续墙2 y 36 40
step 300
m e range group 连续墙2 y 30 36
prop bulk 34.5e9 shear 14.78e9 range group 连续墙2 y 30 36
solve
save 1.sav

xiangruide

这是水平方向上的两个位移图：

xiangruide

下边是竖直方向上的位移图：

                                                   

横之爱

不知道你注意了没有，在你最后那张竖直方向位移图上，基坑外产生都是向上的位移，与坑内隆起位移方向一致，并且靠近基坑的位移比远离基坑的坑后隆起位移要更大一点，常理看来，坑后应该是沉降的，这是什么原因造成的？能否解释？

xiangruide

就是找不出原因想请各位大侠帮忙的，我认为也是不对的。

笑沧海

发现flac算软土总是隆起；认为由于开挖卸荷坑底隆起是正常的，通常开挖10m,隆起1cm；但是周围的土体实际应该是沉降，由于模拟过程中很多因素未能考虑，比如地下水影响，时间效应，地质缺陷、应力释放程度、地连墙刚度等等，计算出来的结果不理想，我在模拟隧道开挖时也有此困惑；

emily127

我现在做的基坑开挖模拟也是这个问题，墙后土体怎么是隆起的呢？而且在连续墙周围的隆起很大！！
不知道有没有方法能够消除这个隆起

qiudh

这好像是flac软件本身的问题了，我在做边坡开挖的时候也会出现这个问题，在开挖上面平面的上也出现的位移向上的情况。

xiangruide

看来遇到这样问题的人很多，可是怎么就没有人帮着解决呢？我听一位大侠说在表面加力来解决，可好像与实际情况不一样，并且这个力加多大也不知道。希望成功模拟的大侠们多指点。

zhongnan

我也遇见过这一问题，后来在表面加力来解决就可以了

笑沧海

好像不应该单纯为了计算结果，随意地施加力或者改变参数之类的吧，呵呵；

xiangruide

我觉得笑沧海说的有道理，如果为了结果而随便加力或者改变参数之类的那还模拟它做什么。

lxrswf

我在竖井开挖时也有些问题

hondu5

我做基坑模拟也遇到这样的问题,真不知道怎么回事,后来我通过改变各个参数,发现土体的粘聚力影响比较大,在开挖到一定的深度,才出现沉降,而且一沉降就不可收拾.很快就达到很大.弄得看起来那些支护好象作用都不大似的.实在是让人琢磨不透,还望各位大狭帮忙解决,或是问培训这个软件的那些老师啊!期待..................

sy1884

研究研究，现在都是讨论抬升问题的啊！

lsjing393

我做的基坑开挖也遇到这个问题