想咨询一下有关flac3d水位下降时的边界问题 [复制链接]

想咨询一下有关flac3d水位下降时的边界问题：水位上升的时候用fish编制了一个简单的函数，就是按照渗流时间（fltime）不断释放，定义以及固定更新边界孔压值和施加库水压力，这个在计算水位上升的时候效果还行，就是太慢，要调整流体模量。虽说这样算出来的就不是真实时间对应的渗流场，但是跟seep对比了下，也可以接受。内部水是可以渗进去的。但是同样的办法在计算水位下降的时候似乎就不太对了！！在计算时间很长以后，虽然边界的孔压值变了，但是坡体内部的渗流场水位依然很高。似乎根本就没向外排水，只有边界上的孔压值变化了，根本就不可接受。跟seep计算结果差很多（倒不是说seep的结果多么的正确，主要是我想通过一个参照来说明我对流体模量的修改造成的影响到什么样的范围）。

不知道哪儿错了，还望指正。

命令流如下：

def waterraise
pp_1=0+fltime*7.6e3/36000
pp_1_2=pp_1-2e3
p1=-1*pp_1
p2=-1*pp_1_2
line=pp_1/1e4
if line<0.2
command
apply remove nstress
free pp ran plane norm 1 0 1.5 ori 0 0 3.46 z 0 line

ini pp pp_1 grad 0 0 -1e4  ran plane norm 1 0 1.5 ori 0 0 3.46 z 0 line
fix pp  ran plane norm 1 0 1.5 ori 0 0 3.46 z 0 line
apply nstress p1 grad 0 0 -1e4 ran plane norm 1 0 1.5 ori 0 0 3.46 z 0 line
endcommand
else
command

;中间差不多，就是对边界分段赋值，命令太长，就不都贴出来了

command

endif
end

set fishcall 10 waterraise

solve age 3600

在水位上升和下降中，思路和程序上是一样的，但是为什么结果差那么多呢？