求助：降雨入渗的边界条件转变

aiter

小弟在做一个降雨入渗问题，边界条件是一大障碍。降雨初期，有多少水入渗多少，这时候边界条件是流量已知的第二类边界条件。降雨一定时间后，土体入渗能力下降，降水一部分不能入渗，这时候边界条件成为水头已知的第一类边界条件。请问这个边界条件转变怎么实现呢？谢谢谢谢····

jjliu

你用的是什么软件啊？商业软件？geo－seep不可以计算边界条件转化的计算的。自己编写软件是可以的。你可以商业软件中可以这样处理，计算在第一种边界条件下计算稳定后的状态作为第二种边界条件计算的初始条件。我也是作非饱和渗流计算的，有机会可互相讨论。

acamus

作降雨产流的有限元计算就行了,得出坡面的入渗比, 然后作为边界条件进行饱和-非饱和计算. 因为地面径流水深值很小, 压力值可以不考虑.

changqing

可以在坡面上定义一虚拟土层，其渗透系数为下面土层的100~1000倍，虚拟土层内的渗流相当于表面径流。
模型定义为非稳定渗流，土坡逐渐饱和，下渗的雨水逐渐变少，不能渗入的雨水沿虚拟土层渗透。

jjliu

版主的办法挺好的，但是非稳定渗流的初始条件如何确定啊？

geoguotao

嘿嘿，好像我师兄也在研究这个问题，我导师说这个没有很好的办法

aiter

jjliu wrote:
你用的是什么软件啊？商业软件？geo－seep不可以计算边界条件转化的计算的。自己编写软件是可以的。你可以商业软件中可以这样处理，计算在第一种边界条件下计算稳定后的状态作为第二种边界条件计算的初始条件。我也是作非饱和渗流计算的，有机会可互相讨论。

这只是权宜之际，我试想解决这个问题，哎。3DFEMwater可以解决这个边界条件转化问题，但是我是用abaqus，好像还没发现可以改。

aiter

changqing wrote:
可以在坡面上定义一虚拟土层，其渗透系数为下面土层的100~1000倍，虚拟土层内的渗流相当于表面径流。
模型定义为非稳定渗流，土坡逐渐饱和，下渗的雨水逐渐变少，不能渗入的雨水沿虚拟土层渗透。

谢谢，我开始也是用这个方法。一维问题好像可以。但是对于斜坡，好像水都流到最低点，不能坡面入渗了，就是雨水优先从渗透系数大的土层渗走了。

aiter

jjliu wrote:
版主的办法挺好的，但是非稳定渗流的初始条件如何确定啊？

我目前采用是土水曲线等都和下层土一样，就是改变渗透系数。后来设置各向异性的渗透系数都没搞定。