请教一下flac3d渗流的浸润线计算 [复制链接]

如图，一个宽12 m，高16 m，厚1 m 的模型，上游水头为12 m，下游水头为4 m，假定饱和状态下的渗透系数为1×10-7，将模型剖分成192 个
网格（12×16×1），如图1 所示。定义上、下游水头边界条件，全部区域都定义饱和状态的渗透系数进行计算，找出孔隙水压力为0 的等值线，即为求
出的初始浸润线，如图2 所示，该线以上的区域为非饱和区，以下的区域为饱和区。

图片:QQ截图20150322200904.jpg

图片:QQ截图20150322200914.jpg

我的目的是得到图2，求出初始浸润线
下面是我写的命令流


new
config fluid
set mech off fluid on
gen zone brick &
p0 0 0 0 p1 12 0 0 p2 0 1 0 p3 0 0 16 size 12 1 16

model fl_iso
prop perm 1e-7
ini fdensity 1e3
ini fmod 2.0e9 ftens -1e3 ;ftens不能为正值


app pp 120e3 grad 0 0 -10e3  ran x -0.1  0.1  z 0 12
app pp  40e3 grad 0 0 -10e3  ran x 11.9 12.1  z 0 4

plot con pp ou on
solve
但是我算出来的结果却惨不忍睹。。

图片:QQ截图20150322201002.jpg

大家帮我看看哪里出错了啊 ！！！

找到原因了 没有加重力，另外ftens的值太小了。。

swallow_wmj:找到原因了 没有加重力，另外ftens的值太小了。。 (2015-03-23 09:09) 

大侠，第三步的命令流有没有研究出来啊...弄半天都不对啊 不知道哪里错

new
config fluid
set mech off fluid on
gen zone brick &
p0 0 0 0 p1 12 0 0 p2 0 1 0 p3 0 0 16 size 12 1 16

model fl_iso

ini fden=1e3
ini fmod=2.0e8 ftens=-1e3 ;ftens不能为正值

prop perm 1e-5 poro=0.3
ini sat=0

set grav 0 0 -10

plot create ppFv
plot add cont pp out on sk magf 20
plot add flow red
plot show ppFv


ini pp 120e3 grad 0 0 -10e3  ran x -0.1  0.1  z 0 12
ini pp  40e3 grad 0 0 -10e3  ran x 11.9 12.1  z 0 4
fix pp ran x -0.1  0.1  z 0 12
fix pp ran x 11.9 12.1  z 0 4

solve

图片:渗流.png

大神，如果在斜面上fix pp得咋弄？

ini   fix  只能计算出饱和区