非饱和渗流关键问题课题.ppt

Limetal于1996年对南洋理工大学校内残积土的基质吸力进行了研究，其均值也大约在20kPa左右。 在加拿大西部曾经开挖两个垂直于铁路路堤的深槽，并沿槽壁安置了快拔型张力计。量测结果用基质吸力等值线表示如图所示。 2、数值模拟初始条件 四、降雨数值分析中几个问题 国内外己有的对降雨条件下非饱和渗流影响以及稳定性方面的研究，一般在处理计算入渗量时，通常有这样几种处理方法： （1）以降雨强度作为流量边界的水流计算入渗速率； （2）把降雨强度按一定（不变）比例降低，作为流量边界的水流计算入渗速率； （3）考虑降雨强度与土体渗透性的关系，决定计算入渗速率。当降雨强度高于饱和渗透系数ks时，计算入渗速率按ks取值，否则按降雨强度取值； （4）当降雨强度小于土壤表面的入渗能力时，计算入渗速率取为降雨强度，边界条件为第二类边界条件;当降雨强度大于土壤表面的入渗能力时，入渗的强度就等于土壤的入渗能力，边界条件转换为第一类边界条件。 3、降雨入渗边界的处理 四、降雨数值分析中几个问题 如果土坡中存在较大或较多的竖向裂隙，而土体本身的渗透性不是很大时，在降雨过程分析中如果不考虑竖向裂隙对入渗过程的影响，则计算得到的斜坡入渗量较小。而实际情况是，斜坡中由于存在竖向节理裂隙，降雨会很快沿裂隙形成垂直入渗。在这种情况下，计算结果与实际情况存在较大的差距。因此有必要在入渗过程中考虑竖向裂隙的影响。 对于无积水的入渗过程（第二类边界条件），降雨直接入渗土体中，土中的裂隙中不存在积水，裂隙的存在对整个非饱和渗流场的变化没有明显影响。这种情况下可以忽略掉裂隙的影响，直接采用无裂隙的土体的非饱和入渗模型，引入边界条件求解即可。而对于有积水的入渗过程第一类边界条件，由于降雨强度大于入渗能力，多余的水在地表形成积水或产生地表径流。此时可以认为雨水充满于竖向裂隙中，渗流场由原来的地表入渗变为地表与裂隙的两个侧面共同入渗。由于入渗边界的显著改变，相应的非饱和渗流场也发生明显的变化，所以在研究第一类边界条件下的入渗问题时必须考虑竖向裂隙的影响。 4、降雨入渗分析中裂隙的处理方法 四、降雨数值分析中几个问题 ①视裂隙两侧为边界。这种处理方法把裂隙的两侧看作为土坡的一部分边界。在数值计算时，重新调整网格的部分范围，把裂隙从所研究的空间去除，使其两侧壁变为边界，引入边界条件，就可以实现考虑裂隙对入渗过程影响的目的。由于有限元方法能够考虑边界复杂的边界条件，这种处理裂隙的方法非常适合于用有限元来求解。 ②等效渗透系数法。在非饱和土斜坡的渗流分析中，裂隙的处理方法还可以采用等效渗透系数法。其基本思想是先按不考虑裂隙进行整个区域的离散化划分，在裂隙所在位置按照裂隙的大小把裂隙剖分为一系列小的薄层单元或网格，将裂隙本身的渗透性等效为薄层的渗透性，当然这一薄层单元或网格小区域的竖向渗透性远大于周围土体的渗透性，在计算中给其很大的渗透系数。当降雨入渗时，雨水迅速沿薄层入渗，很快就可充满这一区域，并向周围土体扩散，这种处理也比较符合实际情况，特别是不需要改变非饱和土的渗流方程及边界条件，只需调整某些单元的渗透性即可实现，是一种处理非饱和土中裂隙渗流的简单有效的方法。 4、降雨入渗分析中裂隙的处理方法 四、降雨数值分析中几个问题 汇报提纲 一、非饱和土特性 二、非饱和渗流基本理论 三、非饱和数值求解关键问题 四、降雨数值分析中几个问题 五、应用实例分析 五、应用实例分析 1、计算模型的建立 ①模型两侧：地下水位以上边界按零流量边界处理，地下水位以下为给定水头边界条件。 ②入渗边界：土坡表面及斜坡处即入渗边界，取为流量边界或给定水头边界。当孔隙水压力小于零时为流量边界，反之变为水头边界。如果雨强小于表层土体渗透性，按流量边界处理，大小为降雨强度，如果雨强大于表层土体渗透性，一部分雨水沿坡面流失，会在坡面形成一薄层水膜，此时可按给定水头边界处理。由于水膜很薄，计算中取水头值等于地表高程。 ③模型底面：假设为不透水边界。 五、应用实例分析 2、计算参数与计算工况 五、应用实例分析 3、初始条件 五、应用实例分析 4、结果分析-1 五、应用实例分析 4、结果分析-2 五、应用实例分析 4、结果分析-3 五、应用实例分析 4、结果分析-4 五、应用实例分析 4、结果分析-5 五、应用实例分析 4、结果分析-6 五、应用实例分析 4、结果分析-7 敬请各位批评指正！ 谢谢！ 汇报题目：非饱和渗流关键问题探讨 GeoStudio使用体会交流 汇报提纲 一、非饱和土特性 二、非饱和渗流基本理论 三、非饱和数值求解关键问题 四、降雨数值分析中几个问题 五、应用实例分析 一、非饱和土特性 基质吸力直接影响着非饱和土的力学性质。对于非饱和土体，其基质吸力的大小则