Geostudio软件模拟降雨入渗过程中边界条件的探讨

在Geostudio软件中研究边坡稳定性问题,常需要考虑降雨入渗的影响。软件有多种入渗边界条件的设置,如选择不当则会造成模拟的失真。文章通过建立假设模型,讨论了不同降雨强度条件下降雨入渗边界设置的问题,发现降雨入渗边界的设置与降雨强度和土体的饱和渗透系数关系密切;当降雨强度小于饱和渗透系数时,入渗边界应设置为流量边界,即降雨全部渗入土体中;当降雨强度大于饱和渗透系数时,入渗边界应设置为水头边界。

基于COMSOL Multiphysics的非饱和裂隙土降雨入渗特性研究

基于COMSOL Multiphysics软件对非饱和裂隙土降雨入渗特性进行数值模拟研究。通过将裂隙和基质分别离散成有限单元,建立了能充分模拟土中裂隙流、基质流以及裂隙-基质流量交换的离散裂隙-孔隙介质模型。结合"空气单元"的概念,对裂隙土的上边界进行模拟。该方法不仅能描述降雨初期雨水沿裂隙优先入渗的现象,还能描述当降雨量大于裂隙土入渗量时雨水沿地表流走的现象。通过对地表以下2 m深度内低渗含裂隙土体进行模拟,分析了裂隙的几何特征、基质的水力特性、前期水分条件以及降雨强度对非饱和裂隙土降雨入渗过程的影响。结果表明,在非饱和裂隙土中,存在两个主要的渗流过程:一是水沿裂隙优先流动;二是水不断从裂隙吸入基质中,基质吸收水的作用抑制了裂隙中优势流的发展。与裂隙的几何特征相比,基质的水力特性对非饱和裂隙土渗流的影响较大。增大基质的饱和渗透系数可能使由裂隙流主导的渗流过程转变为由基质流主导的渗流过程,而基质的非饱和特性与裂隙土的初始含水率改变了土体的储水能力,从而加速或延缓了降雨入渗至某一深度的时间。降雨强度对土体入渗速率和入渗量均有影响,当超过裂隙土的入渗能力时,多余积水沿地表流走,断面入渗率随时间趋于稳定值。

降雨诱发边坡破坏数值模拟两个关键问题的解决方法

研究降雨因素对边坡稳定性的影响,最有效的方法是室内外试验和数值模拟方法,最为重要的是设置与实际情况接近的降雨入渗边界与动态流量边界下初始孔隙水压力分布问题。为有效地解决这两方面的问题,首先将现有的降雨数值边界施加方法所得模拟结果与室内试验结果进行了对比分析,指出了传统方法的不足之处,并在此基础上提出了新的降雨数值边界条件,即在土体表面施加一层很薄的"空气单元",表述的是一个动态的边界条件。降雨入渗边界施加在"空气单元"表面,而不直接施加在土体单元表面,应用该方法完成的入渗土柱数值模拟结果更加接近现场和试验测得的降雨入渗率,为有效模拟降雨诱发边坡破坏的边界条件提供准确、真实可行的方法。同时,由于相对渗透系数和负孔隙水压力均为饱和度的函数,通过设置流体的抗拉强度来赋予节点负孔隙水压力模拟基质吸力,在数值模拟中每一个计算步根据节点饱和度自动更新负孔隙水压力和相对渗透系数,可实现FLAC软件中非饱和土的有效渗流分析。此外,为解决非饱和区土体施加负孔隙水压力所带来的问题,提出数值模拟过程中施加动态的流量边界,通过模拟自然界中地下水的补给,很好地解决使用FLAC进行数值模拟时动态流量边界条件下初始孔隙水压力的分布问题。