非饱与土坡入渗规律的数值模拟.pdfVIP

第六届全国岩土力学数值分析与解析方法讨论会，中国，广州 1998年8月4～6日 非饱和土坡入渗规律的数值模拟 陈善雄 陈守义 (中国科学院武汉岩土力学研究所武汉430071 提要本文采用求解非饱和土水分流动二维问题的积分有限差分程序MFUS2t模拟r降雨条件下非饱和士坡水分人渗 问题。模拟中，有意识地调整水分扩散牢协，导水率凰．得到了有关物性参数时土壤永分人渗影响的规律f同时．又有意识 调整降雨强度I．取得降雨强度对土坡人渗影响的规律；再次，模拟了不同类型土的水分人渗规律。这些规律性成果对于认 识降雨引起的土坡失稳的机理及规律有重要意义。 关键词非饱和土倒；芒萝，数值模拟 。士一。数 一ri／之41．弓3 ’f !．． 1前 言、 一：、芝’毛 }，参一一 土坡水分的入渗与蒸发对于土坡稳定性起着十分重要的作用。已往的工程实例及研究表明，降雨往往 是引起黪胀土边坡、残积土边坡失稳的主要促发困索。地处热带和亚热带的膨胀土、残积土往往是非饱和 的，由于吸力的存在，土体抗剪强度较高．土坡是稳定的，一旦发生降雨，雨水将渗入土体，而致土体饱和度 增加，吸力锐减并引起抗剪强度大幅度下降。当持续降雨的历时和强度超过一定限度时，则可导致土坡失 稳。因此，土坡水分人渗对于土坡的稳定性具有重要意义，要搞清气候变化对土坡稳定性的影响及预测预 报气候影响下的土坡稳定性，首先，必须搞清土坡水分人渗问题，把握其规律。 非饱和土中水分运动由于其参数的高度非线性．求解十分困难，往往只能求其数值廨。以往已有许多 研究者研究过非饱和土水分运动的有限元、有限差分、积分有限差分解法。但对边坡的二维入渗问题研究 并不多见，本文运用求解非饱和土水分流动二维问题的积分有限差分程序MFUS203，模拟分析非饱和土 坡的人渗问题，在计算中有意调整某些条件(包括降雨强度、物性参数)，寻求非饱和水分人渗受这些条件 变化影响的规律。 2计算模型 为考察降雨条件下，土坡水分入渗的情况，选取图1所示截面土坡进行模拟计算。假设该土坡为一均 匀土坡，饱和体积含水率0，一o．40．DAB边界为相对隔水边界，DCB边界受降雨影响，边坡内初始体积台 0 1000 图I土坡剖面尺寸(图中尺寸单位：cm) 488 非饱和土坡人渗规律的数值模拟 胨善雄陈守义 水率研均为0．08。 3水分扩散率变化对土坡水分人渗的影响 6 6 在进行水分扩散率变化对土坡水分人渗影响分析巾，保持导水率不变，K一一10S m／s(s为饱和 度，以下同)，溺整D一时保持曲线形状不变，即函数型式不变。算例中参数变化情况见袭1，袭中，瓿一 103 5”7。trl2／s。 表1￡b变化情况 情 形 l 2 3 4 5 岛 10鼠 13， 南珥 击甄 志礴 降雨强度均取为I一2．5×i0～cm／min。舆型的上坡体积含水率分布见图2。 图2情况2总降雨量为40cm时(序1600min)土垃体积含水率分布 体积含水率沿深度的分布。 从图中可以看出土坡水分人渗强烈受水分扩散率的影响，从情形1至情形5，即D。由太变小，坡顶体 积含水率增加，曲线愈徒·受水分人渗影响的深度变小，即D。越大在相同降雨强度和降雨量情况下，受降 雨影响的深度愈大。 4导水率变化对土坡水分入渗的影响 在进行导水率变化对土坡水分入渗影响分析中，保持水分扩散率不变，D一=103”“9 m。