【请各位大牛、斑竹帮忙看看：水位变动条件下的边坡数值模拟】

51746170

本人在搞水位变动条件下的边坡数值模拟，在论坛中找到一个参考例子，仿照该例子写了相关命令流，可是，计算发现，存在问题很大。

       我用的命令流如下：
new
restore 1103.sav
gen merge 1e-5
attach face all
;-----流体参数及边界-----
config fluid
model fl_iso                 ;均质流体模型
prop perm 4e-8 poro 0.1 range group 2  
prop perm 4e-8 poro 0.1 range group 3
;m fl_n range group 1
prop perm 0 poro 0.01 range group 1
ini fdens 1e3
ini fmod 2e8
ini pp 80e4 range z 65 y 0 330  ;145m高水位
fix pp range z 65 y 0 330
fix x range x -0.1 0.1
fix x range x 1499.9 1500.1
fix y range y -0.1 0.1
fix y range y 1559.9 1560.1
fix x y z range z -0.1 0.1
set grav 0 0 -10
;-----流固计算---------
set fluid pcut on           ;消除负孔隙水压力
set fluid on mech off       ;打开流体计算过程，静力计算关闭。
set flow on
plot cont pp out on
plot set rotation 20 0 30
solve
save flow.sav
      主要问题有：
      1.孔隙水压力沿滑坡体左右两侧的小河严重的超过了所设的水位145m，而两侧小河上游的实际高程超过300m，可以看出，计算才运行了2000多步，还没有计算结束 ；
       2.滑坡体中间出现了零星的孔隙水压力，不知道是不是节点的问题；
       3.模型左上角及模型背面也出现了孔隙水压力，而基岩我是设成了不透水的，不知道是怎么搞的？
       我发现，论坛内对水位的讨论比较少，本人在渗流的理论和模拟方面都不是很好，还望各大牛们从模拟的理论和FISH语言方面给予指教！虽然说我是还没划分网格，但是，我划分网格后试算了一下，还是存在如上的问题！
    图1中，group 1 为基岩，group 2 为滑坡体，group 3为土体，其中group 2 和group 3 参数相同。
    图2为采用m fl_n range group 1指定基岩为不透水，图3为指定基岩渗透系数为0，其空隙水压力稍有区别。
    附件中还有在论坛中找到的参考命令流，也贴上来了，不过忘记是哪个帖子了，在此表示感谢。
    不知道是不是这个命令流不能模拟考虑地形地貌的情况？我把命令流和模型都贴上了，还望各大牛们给看看，不甚感激啊 ！

[ 本帖最后由 51746170 于 2008-11-12 13:10 编辑 ]

51746170

过去一天了，有人看，没人回 自己又研究了一天，有点进展，不过还有问题。

把set grav 0 0 -10去掉，ini pp 70e4 后面加上grad 0 0 -1e4 。可以解决水位沿河向上乱跑的问题，中间零星的孔隙水压力也没有啦。不过产生了两个另外的问题。见图。
问题如下：

1.孔隙水压力最大值一直显示的是1.5e5，也就是15m高的水位，但是图中显示的水位实际已经升到了220多米了？:L
2.本人设定的水位为135m，而实际升到了220m，不知道问题出在哪里？

我想问题是不是出在我的理解上，是不是对原来的那个参考例子理解有问题？

继续期待各位大牛、斑竹、高手支招！ 万分期待，多谢了！

51746170

又一天过去了...
从来就没有什么救世主，也不能靠神仙皇帝——只能靠自己...
自力更生吧。

bluecmycmy

估计还是要自己来了，好多东西都要开发

silinfe

4# bluecmycmy

请问楼主模型是有什么软件建立的啊?谢谢！