论坛 › G0201：FLAC3D/FLAC2D › 渗流初始无法平衡

巧遇 发表于 2008-5-12 16:04:49

渗流初始无法平衡

请大家帮我看看，模拟基坑开挖，渗流场的变化，
(1)但是我渗流初始平衡都无法达到啊？？ 要如何修改
(2) 我开挖后的渗流和力学设置对吗？
谢谢了。

new
config fluid
gen zone brick size 20 1 20 p0 -40 0 -30 p1 0 0 -30 p2 -40 1 -30 p3 -40 0 0
gen zone brick size 10 1 20 p0 0 0 -30 p1 20 0 -30p2 0 1 -30 p3 0 0 0
gen zone brick size 20 1 20 p0 20 0 -30 p1 60 0 -30 p2 20 1 -30 p3 20 0 0
group base
gen zone brick size 30 1 10 p0 0 0 0 p1 60 0 0 p2 0 1 0 p3 0 0 20 group base
gen zone brick size 20 1 10 p0 -40 0 0 p1 0 0 0 p2 -40 1 0 p3 -40 0 20 group exc
model fl_iso
prop perm 1e-19 poro 0.3
ini fdensity 1e3
ini fmod 2.0e9 ftens -1e-3
;ini sat 0
;ini sat 1 range z -30.0 12.0
ini pp 1.2e5 grad 0 0 -1.0e4 range z -30.0 12.0
fix pp range x -40.1 -39.9 z -30.0 12.0
fix pp range x 59.9 60.1 z -30.0 12.0
fix pp range z -30.1 -29.9
; --- mechanical model ---
model elas
prop bul 3.9e9 shea 2.8e9
;model null range group excavate
ini density 1.8e3
;ini density 1.5e3 range group wall
fix x range x -40.1 -39.9
fix x range x 59.9 60.1
fix y
fix x y z range z -30.1 -29.9
set grav 0 0 -10
plot cont pp out on
;plot add flow red
;set fl biot on
set fluid on
;set flow on
set mech off
solve
ini xdis 0 ydis 0 zdis 0
ini xvel 0 yvel 0 zvel 0
set flow on
set mech on
model null range group exc
model fl_null range group exc
fix pp 0 range x -40 0 z 0
solve age 1e5
;save stateflow.sav

巧遇 发表于 2008-5-12 20:14:37

好急啊 ，帮忙啊 。。。

巧遇 发表于 2008-5-13 08:46:08

修改了一下，以下程序可以了。


new
config fluid

gen zone brick size 20 1 20 p0 -40 0 -30 p1 0 0 -30 p2 -40 1 -30 p3 -40 0 0
gen zone brick size 10 1 20 p0 0 0 -30 p1 20 0 -30p2 0 1 -30 p3 0 0 0
gen zone brick size 20 1 20 p0 20 0 -30 p1 60 0 -30 p2 20 1 -30 p3 20 0 0
group base

gen zone brick size 30 1 10 p0 0 0 0 p1 60 0 0 p2 0 1 0 p3 0 0 20 group base
gen zone brick size 20 1 10 p0 -40 0 0 p1 0 0 0 p2 -40 1 0 p3 -40 0 20 group exc

model fl_iso
      prop perm 1e-19 poro 0.3
ini fdensity 1e3
ini fmod 2.0e9 ftens -1e-3

;ini sat 0
;ini sat 1 range z -30.0 12.0
ini pp 1.2e5 grad 0 0 -1.0e4 range z -30.0 12.0
fix pp range x -40.1 -39.9 z -30.0 12.0
fix pp range x 59.9 60.1 z -30.0 12.0
fix pp range z -30.1 -29.9

; --- mechanical model ---
model elas
prop bul 3.9e9 shea 2.8e9
;model null range group excavate
ini density 1.8e3
;ini density 1.5e3 range group wall
fix x range x -40.1 -39.9
fix x range x 59.9 60.1
fix y
fix x y z range z -30.1 -29.9

set grav 0 0 -10

plot cont pp out on

;plot add flow red
;set fl biot on
set fluid on
;set flow on
;set mech off   ; 计算时首先要达到力学初始平衡
;set fluid substep 10
;set mech substep 100 auto

solve
ini xdis 0 ydis 0 zdis 0
ini xvel 0 yvel 0 zvel 0

set flow on
set mech on
model null range group exc
model fl_null range group exc
fix pp 0 range x -40 0 z 0
fix pp 0 range x 0 z 0 20; 边界的孔压可设为零

solve
;save stateflow.sav

[ 本帖最后由 巧遇 于 2008-5-13 08:52 编辑 ]

巧遇 发表于 2008-5-13 09:32:07

初始平衡时若不约束水平孔压，效果更好。

巧遇 发表于 2008-5-13 09:35:57

若不约束底部渗流边界，则初始平衡就有问题

巧遇 发表于 2008-5-13 09:43:50

设置了水力梯度的范围后，在水位面也不需要约束孔压为零了，
要不然在计算过程中，会出现负孔压间断出现，跳跃出现的现象
不知道是否正确。？

巧遇 发表于 2008-5-13 09:46:02

经过几次试验，学到了，在计算之前要进行力学平衡，不要先关掉力学计算。


设置了水力梯度的范围后，在模型底部要约束孔压，
但若没有设置水力梯度，则无需在底部约束孔压，而让程序自动计算其他边界的孔压

[ 本帖最后由 巧遇 于 2008-5-13 10:07 编辑 ]

lovemm 发表于 2008-5-13 22:54:31

回复 2# 的帖子

:lol ,我以前没有接触流体，但还是顶:victory:

zzy211321 发表于 2008-6-25 09:43:53

不错，只是还不是很懂，正在学习中

zzy211321 发表于 2008-6-25 10:23:06

为什么在设置地应力平衡时没“set fluid off ”呢？不约束水平孔压是什么意思？等待答复中，谢谢

070542 发表于 2008-6-25 10:58:27

请问一下我用的FLAC 3D为什么识别不了 set flow on7 和

minhrac 发表于 2008-7-16 14:48:53

原帖由 070542 于 2008-6-25 10:58 发表 http://www.simwe.com/forum/images/common/back.gif
请问一下我用的FLAC 3D为什么识别不了 set flow on7 和

可能你没有指定config fluid

data_306 发表于 2008-9-5 10:36:01

巧遇同志，你的初始平衡 ：是不是没有计算力学平衡啊？只做了个流体平衡？
;set mech off   ; 计算时首先要达到力学初始平衡
你关掉了力学平衡

songjinglei 发表于 2008-9-5 15:51:44

不错，只是还不是很懂，正在学习中

syx2006 发表于 2009-2-20 13:32:57

你的计算情况没有考虑基坑开挖过程有效应力消散过程

ddkk1463 发表于 2009-3-26 16:21:36

谢谢搂主,我也正在学习中

yantumi 发表于 2009-5-5 20:56:31

学习了，谢谢。

leoathit 发表于 2009-8-26 16:49:26

有收获，学习

binggehaoyang 发表于 2011-11-6 16:06:21

我觉得你的水压赋得有错误    pp 1.2 grad 0 0 -1.0e4 range z -30.0 12.0这儿 应该是4.2e5

binggehaoyang 发表于 2011-11-6 16:26:00

巧遇 发表于 2008-5-13 08:46 static/image/common/back.gif
修改了一下，以下程序可以了。




你没发现他的边界条件X =-40   60 以及底面都设成了透水边界把孔压固定了吗   我感觉是不是有点错误

查看完整版本: 渗流初始无法平衡