斜坡降雨入渗在Flac中模拟

为了简化分析,在模拟斜坡降雨入渗暂态渗流时,通常没有考虑入渗和产流的耦合过程,通过施加流量边界或零水头边界进行地下水渗流计算,通过在FALC中利用内嵌的fish语言,考虑地表径流和地下渗流的耦合以及雨向影响下斜坡降雨入渗过程,并通过算例分析了考虑和不考虑耦合、对斜坡孔隙水压力分布以及考虑各种雨向作用下对坡面积水深度的影响。     

改进的斜坡降雨入渗与坡面径流耦合算法研究

大气降雨时坡面径流与非饱和入渗的耦合计算,是降雨诱发滑坡及水土流失等研究中的重大问题之一。从坡面径流和降雨入渗控制方程着手,建立了基于有限元方法的耦合方程求解模型,实现了耦合问题的直接求解。该耦合计算模型与求解方法避免了坡面径流和降雨入渗间的迭代计算,也不需计算两者间的流量交换,较大程度地提高了计算效率和计算精度。算例表明,该方法正确可行,计算效率高,稳定性好。     

考虑径流影响的滑坡降雨入渗二维有限元模拟及应用

降雨是滑坡失稳的常见诱发因素之一,现有数值模拟方法由于不能考虑来自基岩边坡的径流对滑体渗流的补给而低估了真实的降雨入渗水量,导致在评价降雨对滑坡稳定性特别是深层滑动影响时存在一定不足。以Richards方程和有限元法为基础,忽略降雨对渗透性极低的滑床和滑带的影响,将滑坡渗流计算域缩小为滑体,以避免因滑体与滑床（滑带）渗透性差异巨大引起的数值计算困难;依据基岩边坡水平长度和滑体降雨入渗边界饱和情况,修正降雨入渗边界,实现了考虑径流补给的滑坡降雨入渗简化数值模拟。算例表明,该方法所得渗流场更加符合实际情况。以三峡库区某滑坡为例,模拟了2006年10月～2009年12月间滑坡在库水升降和降雨条件下的渗流场和安全系数的演化过程;计算结果表明,考虑径流补给时滑坡后部的渗流场饱和区域明显较大,稳定系数降低较多,与位移监测资料显示的佐证较为吻合;若不考虑径流补给,则降雨对滑坡稳定性影响不大。     

降雨条件下坡面径流和入渗耦合的数值模拟

降雨条件下坡面径流和降雨入渗的模拟互为条件,必须采用耦合分析的方法才能较好地解决这一问题。从坡面径流和降雨入渗控制方程着手,建立了用有限元方法求解该问题的耦合方程,并编制了相应的计算程序。为了加快迭代收敛的速度,还对地表饱和的判断提出了一种近似处理的方法。算例表明,所提出的方法可行,能较好地反映实际现象。该方法可为边（滑）坡稳定、水土流失、山地洪水分析等提供计算依据。     

考虑气压势影响的降雨入渗数值模拟研究

利用一维下渗模型的数值模拟方法研究下渗试验过程中气压势与饱和水力传导度K_s的关系。数值分析表明,可通过由下渗过程中表层湿润的饱和度来折减K_s值的途径考虑气压势的影响,从而改进计算方法。下渗试验的模拟计算结果表明,改进后的方法比常用的方法有更高的精度。     

降雨条件下道路边坡地下水渗流分析

降雨是影响道路边坡稳定的主要因素之一,根据岩土饱和-非饱和渗流理论,考虑降雨入渗的影响,利用有限元方法,对强降雨条件下公路边坡地下水渗流场动态进行了数值模拟,得到了边坡地下水压力水头、总水头变化、流速的变化规律,为道路边坡的稳定性分析和滑坡预测提供重要的分析数据。     

降雨入渗规律的实测与分析

关于降雨入渗可降低土坡的稳定性的研究很多,其机制也逐渐趋明确。然而,一场降雨对土坡稳定产生影响的入渗雨量计算问题还没有得到很好地解决。Holten等研究者的一些计算入渗量的方法在土坡的实际应用中难以满足精度要求;而一些有限元的边界处理方法则建立在假定的入渗规律基础上,多少存在不足。研究设计了室内降雨入渗土柱试验,在初步揭示降雨入渗过程和规律基础上,探讨了准确反映降雨入渗量的有限元计算方法,并通过数值试验进一步揭示室内试验未能明确的降雨入渗规律。     

降雨入渗条件下边坡岩体饱和非饱和渗流计算

简要分析了连续介质饱和非饱和渗流数学模型,并讨论了渗流有限元计算中的有关问题,同时研究了降雨入渗机理及模拟方法,在此基础上编写了非饱和渗流程序SUSC。运用该程序模拟了某边坡降雨过程渗流场的变化情况。计算结果表明,在降雨入渗作用下边坡顶部迅速被饱和,随后表面雨水逐渐向边坡深部下渗,形成从坡顶往深部压力水头等值线由高(零)→低→高的封闭现象。随着降雨的继续,边坡顶部负压区进一步缩小,且负压绝对值减小。降雨结束后,由于上部地下水的继续下渗,在边坡的一定范围、一定时间内压力水头继续升高。根据计算结果可知,局部地方的压力水头最高值出现在降雨结束2 d左右,往后整个边坡非饱和区地下水压力水头全面下降,逐渐恢复原状。模拟结果总体可靠,可作为边坡稳定性分析评价及边坡排水加固的参考依据。同时表明,运用上述饱和非饱和渗流模型及降雨模拟方法,计算降雨条件下边坡岩体的渗流场是可行的。     

成都黏土降雨入渗规律研究

成都黏土作为一种典型的膨胀土,具有很强的水敏性,其降雨入渗规律较复杂.结合现场试验和数值计算结果,分析不同降雨条件下的降雨入渗规律.研究表明:降雨入渗深度范围内大部为非饱和区,降雨入渗深度主要受降雨时长、孔隙率和基质吸力控制,降雨入渗饱和深度主要受降雨强度、孔隙率和饱和渗透系数控制,成都地区单日降雨情况下,降雨入渗深度...     

考虑优势流作用的降雨入渗边坡可靠度分析

降雨诱发滑坡是世界上最普遍的地质灾害,而关于降雨入渗对边坡的可靠度分析,优势入渗的影响被长期忽视。使用Comsol Multiphysics对降雨条件下优势入渗做数值求解,利用无限边坡模型计算边坡安全系数;并应用改进Cholesky分解法生成空间相关随机场,使用蒙特卡洛方法分析降雨过程中边坡的可靠度。结合确定性与可靠度计算对比均质入渗与优势入渗在降雨过程中边坡安全性变化：（1）降雨强度较低时,优势入渗安全性较好,而降雨强度较高时,均质入渗更稳定;（2）均质入渗中参数空间变异性是边坡失稳破坏的关键因素,而优势入渗的边坡失稳则由湿润峰快速推进所导致;（3）针对优势入渗模型研究,发现基质域与优势域水力交换强度较大时边坡有更大概率失稳,而较小的水力交换强度可能影响边坡底部的失稳破坏。     

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

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

降雨过程中边坡临界滑动场

降雨入渗过程中孔隙水压力的升高与基质吸力的降低引起边坡稳定性的下降,是导致边坡滑塌的主要诱导因素。利用饱和-非饱和渗流有限元计算得到的孔隙水压力场,基于Fredlund提出的非饱和土抗剪强度理论,对边坡临界滑动场进行改进,提出可以考虑降雨过程的边坡临界滑动场数值模拟方法,能够方便、快速地计算出边坡局部、整体安全系数和相对应的临界滑动面在降雨过程中的变化历程。将该法用于一个典型均质边坡和一个非均质边坡在降雨过程中的稳定性计算,分析降雨持续时间、降雨强度和非饱和强度参数取值等因素对边坡稳定性的影响,并将计算结果与其他方法进行比较,结果表明临界滑动场方法能搜索任意形状最危险滑面,计算的安全系数合理。     

降雨入渗条件下非饱和路堤变形与边坡的稳定数值模拟

对于常常处于非饱和状态的路堤在降雨入渗条件下变形问题的数值模拟研究,由于涉及到非饱和土力学,一直是个热点和难点。现有的强度折减有限元法一般仅对强度指标进行折减,但折减后得到的路堤变形场与实际变形场存在一定差异。论文分析认为,在降雨入渗影响下非饱和路堤不仅强度参数会降低,而且变形参数也会随着应力水平和强度参数的变化而变化。采用基于变模量的弹塑性强度折减法,对不同降雨入渗深度下路堤的变形规律进行数值模拟,并与室内模型试验结果进行比较;同时采用强度折减法和极限平衡法中的Bishop法,分析了不同降雨入渗深度对非饱和路堤稳定性和滑弧变化特征的影响规律。     

降雨入渗条件下软岩边坡稳定性分析

降雨条件下软岩边坡的失稳模式是进行边坡处治的基本依据之一。为研究软岩边坡在降雨条件下的稳定性,提出了一种既能考虑渗流场对边坡稳定性的影响,又能体现岩石软化效应对边坡失稳所带来的不利作用的新方法。运用二维渗流数值计算方法,对降雨条件下的边坡孔隙水压力大小及暂态饱和区面积在空间及时间上的分布进行了模拟。并将渗流场计算结果与暂态饱和区岩石软化试验所得岩石物理力学参数随时间的取值相结合,采用强度折减法分析软岩边坡在降雨入渗条件下的稳定性。对算例边坡的研究表明：降雨入渗条件下软岩边坡的失稳在降雨初期表现为边坡表层局部分层垮塌。随着降雨历时的增长,失稳形式则表现为局部分层垮塌与整体滑移相结合。降雨停止后,边坡负孔隙水压力的消散,对软岩边坡安全系数的继续降低具有一定的延缓作用。     

降雨作用下煤系土路堑边坡稳定分析

基于掌握的降雨气象资料,利用非饱和瞬态渗流有限元程序对开挖后和风化后两种情况下的煤系土路堑边坡进行了降雨入渗坡体渗流场数值计算和分析。依据非饱和土力学强度理论,基于Morgenstern-Price方法,调用不同降雨时间的渗流计算结果,对两种情况下煤系土路堑边坡的瞬态稳定性进行了计算和分析。研究结果表明,水分在风化煤系土路堑边坡体内的渗透迁移对路堑边坡稳定性和最危险滑动面有时空影响;在煤系土路堑边坡降雨入渗和稳定性分析中考虑风化层较为合理。     

A semi‐analytical solution of surface flow and subsurface flow on a slope land under a uniform rainfall excess

A new approach is proposed to analyze the surface flow and subsurface flow passing over a pervious ground under a uniform rainfall excess. The flow field is divided into two regions that are called water layer and soil layer. To figure out the hydraulic behavior of overland flow on an inclined plane under a rainfall event, the simplified Navier–Stokes equations are employed for the surface water flow, and the flow inside the soil layer is porous media flow, which is governed by Biot's (1956, 1962) theory of poroelasticity. The velocity distribution of overland flow is nonzero at the ground surface. The relation between water depth and slope length was developed first. The profile of surface water flow was then found backwards from the downstream end of the flow section by the Runge–Kutta method. After that, the flow velocity and flow discharge of each layer could also be obtained via the water depth. Finally, the variation of fluid shear stress inside the soil layer is also discussed. Copyright © 2011 John Wiley & Sons, Ltd.