降雨入渗参数对粗粒土路堤暂态饱和区影响的数值模拟.pdf

第42卷第3期河海大学学报(自然科学版)Vo1．42No．32014年5月JournalofHohaiUniversity(NaturalSciences)May2014DOI：10．3876／j．issn．1000—1980．2014．03．013降雨入渗参数对粗粒土路堤暂态饱和区影响的数值模拟曾铃，付宏渊，周功科(1．长沙理工大学土木与建筑学院，湖南长沙410076；2．长沙理工大学交通运输工程学院，湖南长沙410076)摘要：为探讨降雨入渗条件下粗粒土路堤边坡暂态饱和区的形成条件，分析降雨入渗参数对暂态饱和区发展规律的影响，基于室内外试验数据与实际气象数据，拟定不同的数值计算方案，采用降雨入渗理论对降雨入渗参数影响下粗粒土路堤边坡雨水入渗过程进行数值模拟。结果表明：降雨强度g与饱和渗透系数k的比值是决定路堤内雨水入渗方式及暂态饱和区形成的关键因素，当q／k≥1．5时，路堤边坡表面及坡脚水平面将形成一定深度的暂态饱和区，暂态饱和区的形成增大了路堤发生失稳的可能；降雨入渗仅会引起路堤边坡坡面以下水位大幅度上升，对路面中部以下地下水位影响有限；路堤边坡暂态饱和区的扩展速度主要取决于g与k的大小，而暂态饱和区形成的时间与扩展面积则受到降雨强度、降雨历时与饱和含水率的共同影响。关键词：道路工程；暂态饱和区；降雨入渗；粗粒土路堤；数值模拟中图分类号：U416．12文献标志码：A文章编号：】000—1980(2014)03—0250—07Numericalsimulationofefectsofrainfallinfiltrationparametersontransientsaturatedareasofcoarse—grainedsoilembankmentZENGLing，FUHongyuan，ZHOUGongke(1．SchoolofCivilEngineeringandArchitecture，ChangshaUniversityofScienceandTechnology，Changsha410076，China；2．SchoolofTrafficandTransportationEngineering，ChangshaUniversityofScienceandTechnology，Changsha410076，China)Abstract：Inordertoexploretheformationconditionsoftransientsaturatedareasofacoarse．grainedsoilembankmentslopewithrainfallinfihration．andtoinvestigatetheeffectsofrainfallinfihrationparametersOilthedevelopmentofthetransientsaturatedareas，wedesignedvariousnumericalcalculationschemesusingfieldandlaboratorydataandactualmeteorologicaldata．andconductednumericalsimulationoftherainfallinfiltrationprocessofthecoarse—grainedsoilembankmentslopeundertheinfluenceofrainfallinfihrationparametersbasedonthetheoryofrainfallinfihration．Theresultsshowthattherainfallinfiltrationdirectionandformationoftransientsaturatedareasweredeterminedbytheratioofrainfallintensity(q)tosaturatedpermeabilitycoeficient(kt)．Whenq／kt≥1．5，transientsaturatedareaswithacertaindepthformedattheembankmentslopesurfaceandslopetoe．Theformationofthetransientsaturatedareasincreasedthepossibilityofembankmentfailure．Therainfallinfihrationonlycausedthewaterlevelundertheembankmentslopesurfacetoincreasesignificantly．buthadlimitedinfluenceonthegroundwaterlevelbelowthecentralroad．Theexpandingspeedofthetransientsaturatedareaswasmainlydeterminedbyqandkt，whiletheformationtimeandexpansionareaofthetransientsaturatedareaswereinfluencedbyq，rainfallduration，andsaturatedwatercontent．Keywords：roadengineering；transientsaturatedareas；rainfallinfiltration；coarse—grainedsoilembankment；numericalsimulation收稿日期：2013—03—11基金项目：国家自然科学基金(51078042，51278067)；湖南省科学技术厅科技计划重点项目(2013FJ2008)作者简介：曾铃(1986一)，男，重庆江津人，博士研究生，主要从事公路工程研究。E-mail：zlbingqing3@126．conl 第3期曾铃，等降雨入渗参数对粗粒土路堤暂态饱和区影响的数值模拟路堤的稳定是确保公路服务能力的关键。对于公路路堤而言，其稳定性受降雨、填料性质、施工工艺及车辆荷载等因素的影响，在众多因素中，降雨入渗导致的路堤填料含水率增加、基质吸力降低是路堤变形与失稳的主要原因[1-4]。降雨入渗造成的路堤失稳容易导致路面结构应力集中或过早破损，从而影响公路的服务能力，降低行车安全性。因此，了解降雨条件下路堤稳定性及内部渗流场变化规律对路堤的设计与施工具有指导意义，许多学者对此进行了深入研究，并得出许多可靠的结论。陈晓斌等认为降雨人渗对路堤变形影响明显，入渗作用扩大了路堤塑性区范围，使路堤稳定性降低。李佳等认为水在非饱和黄土中的渗流过程是暂态饱和区不断下移扩大的过程。付宏渊等认为暂态饱和区的形成时间由降雨强度决定，而暂态饱和区扩展面积大小则受降雨历时和降雨强度共同影响，暂态饱和区的形成及扩展范围与初始地下水位线密切相关。邹维列等认为新建路堤随着荷载、降雨次数和降雨量的增大，堤坡侧向出现破坏的可能性最大；对于上覆荷载小的堤坡部位，若采用具有膨胀性的黏土填筑，压实度不宜过高；路堤坡脚附近含水率与基质吸力受降雨量影响最大。2010年，Tasi等研究了降雨入渗条件下非饱和土体密度与抗剪强度参数的关系，认为在密度增大与抗剪强度参数降低的情况下，极易引起路堤边坡失稳。Oh等¨。。通过对压实土路堤边坡的研究，得出降雨引起路堤表层含水率增大是造成非饱区浅层破坏的主要原因。尽管人们对降雨条件下边坡的渗流过程及暂态饱和区形成特点有了一定的认识，但对于填方路堤而言，不同地域的气候条件及所用的填筑材料均存在差异，其渗流过程比边坡复杂。目前，针对路堤降雨入渗的研究主要侧重于降雨对路堤边坡稳定性的影响及边坡内部渗流过程两方面，而对于暂态饱和区的形成条件，以及不同降雨人渗因素影响下暂态饱和区形成时间与空间变化规律的研究较少。因此，有必要开展相关研究。笔者以某高速公路路堤断面为研究对象，结合室内外土工试验数据及实际降雨资料，通过数值模拟研究降雨条件下粗粒土路堤边坡暂态饱和区的形成条件，以及在降雨强度、土体饱和渗透系数、土体饱和含水率变化下粗粒土路堤暂态饱和区面积在时间及空间上的分布规律，为进一步分析暂态饱和区对粗粒土路堤稳定性的影响提供参考。降雨入渗基本理论1．1降雨入渗过程分析经过碾压密实后的公路路堤一般处于非饱和状态¨，降雨过程中路堤非饱和区雨水的人渗与饱和区地下水的运动互相联系，雨水在路堤内部的入渗过程属于典型的饱和一非饱和渗流过程。根据Coleman和Bodman的研究¨引，当均质土体表面有积水入渗时，典型含水率分布剖面从上往下可分为4个区，即表层有一厚度较薄的饱和区，以下是含水率变化较大的过渡区，其下是含水率分布较均匀的传导区，最下层是湿润程度逐渐减小的湿润区，湿润区的前缘称为湿润锋。随着降雨时间的持续，传导区会不断地向路堤内部发展，湿润区与湿润锋会向路堤内部移动，路堤内部的含水率分布曲线逐渐平缓。笔者基于这一认识，对暂态饱和区的形成条件及变化规律进行研究。1．2降雨入渗计算模型由于降雨入渗的影响因素众多，且人渗过程较为复杂，为了便于全面揭示降雨人渗过程，人们提出了多种渗流模型lI。。渗流模型主要分为2类：一类是基于长期现场观测数据建立的半经验模型；另一类是基于渗流理论(如Richards方程)，通过对实际问题进行一定的假设和限定建立的模型，这类模型以Green．Ampt模型为代表。对于实际工程中的土坡，其内部土壤与边界条件较为复杂，无法运用Green—Ampt模型直接求其解析解，只能根据土壤水运动的基本方程和初始条件，采用有限元或有限差分法求出近似解。本文研究的问题属于典型的二维饱和一非饱和渗流问题，其控制方程为()+(，)+仲=mwpwgoh(1)式中：——岩体裂隙中的总水头；后|j}—、Y方向的渗透系数；——源汇项；m——比水容量；p——水的密度；g——重力加速度；￡——时间。渗透系数k为fk饱和区，^、，=1【()0≤<1zJ 河海大学学报(自然科学版)第42卷式中：Jj}——饱和渗透系数；(0)——非饱和土的渗透系数；体积含水率。其中非饱和区渗透系数随含水率与基质吸力的变化而变化，本文采用vanGenuchten模型培对其进行拟合。2数值计算模型及方案确定2．1数值计算模型公路路堤计算模型的路面宽为60m，填方高度为12m。路堤中使用土一石(崩解完成后的炭质泥岩)混填粗粒土填料，由于路堤两侧对称，故选取半剖面作为研究对象。路堤典型剖面二维有限元网格如图1所示。计算模型单元数为2997，节点数为3114，为了保证计算精度，通过设置辅助线将网格划分为四边形单元。渗流边界条件：模型底部、两侧铅直位置、路面设置32为不透水边界。当程序运行时，自动判断降雨强度q与土体渗透性的关系。如果g小于路堤坡面、坡脚水平面g20表层土体的渗透性，按流量边界处理，大小为降雨强度；嚣，16一如果q大于路堤坡面、坡脚水平面表层土体的渗透性，一8：部分雨水沿路堤坡面与坡脚水平面流失，会在坡面形成405lO1520253O354o4550556O6570一薄层水膜，此时可按定水头边界处理。由于水膜很薄，宽度／m计算中取水头值等于地表高程。图2为公路路堤降雨人图1路堤典型剖面二维有限元网格渗形象化过程。Fig．1Two·dimensionalfinitedeme~modelof渗流初始条件：根据路堤坡面及坡脚多个位置的钻embankmenttypicalsection孔取样资料，设置如图1所示的初始地下水位进行稳态流计算，以此结果作为边坡初始渗流场。渗流水力学模型：采用各向异性达西渗流模型。渗流计算水力特性及参数：预崩解炭质泥岩粒径多位于0．075—10mm之间，属于中粗粒组，并且崩解后的炭质泥岩泥化成分较大，在压实度为93％的炭质泥岩一土混填情况下，路堤填土具有很高的密实性。通过对现场获取的多组原状土样进行室内饱和渗流试验，测得压实图2路堤降雨入渗形象化过程度为93％的炭质泥岩路堤填料饱和渗透系数为Fig．2Rainfallinfiltrationprocessofembankme~9．5×10～一7．36~10一m／s。本文选取饱和体积含水率为0．15，饱和渗透系数为3．OxlO～m／s，由于路堤非饱和区基质吸力的存在，路堤内部非饱和区渗透系数为非饱和渗透系数，采用vanGenuchten模型对非饱和区渗透系数进行拟合，渗透系数及体积含水率随负孔隙水压力P的变化见图3、图4。此外，由于路堤填筑碾压过程造成了其内部土体的各向异性，取各向异性渗透系数比／：0．5。p／MPa图3渗透系数与负孔隙水压力关系曲线图4体积含水率与负孔隙水压力关系曲线Fig．3RelationshipbetweenpermeabilitycoefficientandFig．4Relationshipbetweenvolumetricwatercontentandnegativeporewaterpressurenegativeporewaterpressure 第3期曾铃，等降雨入渗参数对粗粒土路堤暂态饱和区影响的数值模拟2．2数值分析方案为全面了解粗粒土路堤在降雨条件下的渗流特性，获得不同参数影响下路堤内部渗流过程的变化规律，笔者通过参考当地气象资料，并结合岩土体室内外试验数据，对降雨条件下炭质泥岩粗粒土路堤边坡在降雨强度、土体饱和渗透系数以及土体饱和含水率单因素变化情况下路堤边坡暂态饱和区面积在时间及空间上的分布规律进行探讨。主要数值计算模拟方案见表1。3降雨条件下暂态饱和区形成条件分析对表1中数值计算方案进行初步计算可知，q／k是决定路堤边坡暂态饱和区出现与否的关键控制因素。因此，为了对降雨条件下暂态饱和区的形成与发展有更深入的了解，在方案3的基础上，保持其q不变，不断改变k的大小，使q／k∈[1，4]。经大量试算可知，降雨过程中q／k的大小是决定路堤内部雨水入渗方式与暂态饱和区形成的关键因素。当q／k。≥1．5时(图5(a))，降雨持续一段时间后，路堤边坡坡面及坡脚水平面位置会逐渐出现一定深度的暂态饱和区，暂态饱和区的面积詈随可『疆薹着降雨历时的增长而逐渐扩大，同时坡脚附近的地下水位也逐渐升高，并在坡脚附近率先与地下水位面连通，连通后路堤边坡坡面以下水位出现大幅度上升；当q／k<1．5时(图5(b))，路堤边坡随着降雨的持续，在其坡面与坡脚水平面并没有暂态饱和区形成，而是通过雨水的持续入渗直接引起坡脚附近地下水位上升。从图5可以看出，坡脚附近地下水位最易受雨水人渗影响，坡脚以下土体也最先趋于饱和，饱和区范围由坡脚附近迅速向周围土体扩散，并逐渐覆盖坡脚水平面以下区域，同时由坡脚处沿着坡面向上延伸。在整个降雨过程中，雨水入渗仅会引起路堤边坡坡面以下水位大幅度上升，对道路路面以下地下水位影响有限。路堤边坡表面(a)方粟3(b)方粟7图5暂态饱和区及地下水位随时间的变化Fig．5TransientsaturatedareasandgroundwaterlevelVS．time暂态饱和区的形成及其对路堤渗流场的影响是导致路堤失稳的因素之一，通过在路堤边坡表面附近不同位置设置特征点A、B、C(图1)，研究路堤在降雨条件下孔隙水压力与体积含水率的变化规律(图6，限于篇幅，只列出方案3与方案7的变化过程对比曲线)。从图6可以看出，方案3中特征点A、B、C的孔隙水压力完全丧失所需的时间分别为28h、39h、44h，孔隙水压力出现的最大值分别为598．23kPa、414．18kPa、216．34kPa，土体饱和所需时间分别为28h、35h、39h；方案7中特征点A、B、C的孔隙水压力完全丧失所需的时间分别为30h、37h、49h，孔隙水压力出现的最大值分别为323．55kPa、213．36kPa、98．98kPa，土体饱和所需时间分别为30h、38h、49h。由此可知，暂态饱和区的出现，缩短了路堤内部非饱和区基质吸力完全丧失 河海大学学报(自然科学版)第42卷以及土体饱和所需的时间，增大了孔隙水压力的变化幅度，尤其在坡脚处的特征点A孔隙水压力变化最显著，其原因是人渗雨水在重力作用下迅速向坡脚汇集。依据加拿大学者Fredlund基于摩尔一库伦公式提出的非饱和土体抗剪强度理论可以得出_1：暂态饱和区的出现将导致路堤内部孔隙水压力升高、土体密度增大，从而引起土体抗剪强度减小、路堤边坡下滑力增大，造成路堤边坡发生失稳的可能增大。由于坡脚附近孔隙水压力最大、体积含水率最高，因此边坡发生整体滑移的剪出口位置最有可能出现在坡脚附近。一方案3特征点A一方案3特征点C一方案7特征点B一方案3特征点B一方案7特征点A一方案7特征点Ct／ht／h(a)孔隙水压力(b)体积含水率图6降雨过程中孑L隙水压力与体积含水率对比曲线Fig．6Comparisonofporewaterpressureandvolumetricwatercontentduringrainfall3．1不同降雨强度影响下路堤边坡暂态饱和区变化规律分析为探讨q对路堤边坡暂态饱和区变化规律的影响，分别对第1组所列工况进行计算，并选取不同时刻路堤内部暂态饱和区空间分布及暂态饱和区面积随时间的变化进行对比，如图7所示(图中编号均为表1中的方案编号)。从图7可以看出，当q=1．85xlOm／s(方案1)时，经过36h降雨，路堤内部依然没有饱和区出现，而q=6．Ox10m／s的方案4中，暂态饱和区已经发展到路堤的中上部。由图7(b)可知，当降雨持续72h后，各方案暂态饱和区持续增大，但大小关系与36h保持一致。由此可以推断：q越大，雨水在路堤内部扩展的速度越快，在相同降雨时间内雨水的人渗量也越大。由计算可知，第1组中4种工况暂态饱和区出现所需的时间(坡脚附近地下水位线上升所需时间)分别为40h、30h、17h、12h，所形成的最终面积分别为98．35m、227．62m、326．28m、406．95m。由此可知，在降雨时间相同的情况下，路堤内部由雨水人渗所形成的暂态饱和区面积与q成正比；q越大，路堤内部暂态饱和区出现所需的时间越短，最终所形成的面积越大。因此，在路堤的设计阶段，为确保路堤的长期稳定，应充分考虑当地气候条件对路堤运营期稳定性的影响，做好路堤边坡防排水设施的设计。图7不同降雨强度影响下暂态饱和区及地下水位对比Fig．7TransientsaturatedareasVS．groundwaterlevelunderinfluenceofva~ousrainfallintensities3．2不同饱和渗透系数影响下路堤边坡暂态饱和区变化规律分析图8为不同饱和渗透系数条件下不同时刻路堤边坡内暂态饱和区空间分布对比(图中编号均为表1中的方案编号)。从图8可以看出，降雨持时分别为36h与72h时，在不同饱和渗透系数影响下，路堤内部水位线的变化并不明显。第Ⅱ组的4种工况中暂态饱和区出现所需的时间分别为8h、17h、18h、18h，所形成的最终面积分别为342．93m、326．28m、324．1m、322．15m，即在相同的降雨强度与降雨时间情况下，改变路堤土体的饱和渗透系数，对路堤内部暂态饱和区出现所需的时间与扩展范围影响较小。3．3饱和含水率变化对路堤边坡暂态饱和区变化规律的影响为分析不同饱和含水率条件下路堤边坡暂态饱和区的变化规律，分别对第Ⅲ组所列工况进行计算，由此 第3期曾铃，等降雨入渗参数对粗粒土路堤暂态饱和区影响的数值模拟边坡表面图8不同饱和渗透系数影响下暂态饱和区及地下水位对比Fig．8TransientsaturatedareasVS．groundwaterlevelunderinfluenceofvarioussaturatedpermeabilitycoeficients得出路堤边坡内部暂态饱和区空间分布，如图9所示(图中编号均为表1中的方案编号)。由图9可以看出，在降雨时间与初始含水率均相同的情况下，4种工况中暂态饱和区的空间分布与出现的时间差异较大。4种工况暂态饱和区出现所需的时间分别为11h、17h、22h、27h，所形成的最终面积分别为440．5m、326．3m、246．1m、197．5m。以上分析可知：路堤填料的饱和含水率越低，达到饱和所需的时间越短，雨水在路堤内部的入渗速度越快，对降雨条件下路堤的稳定性越不利。坡表面(a)36h图9不同饱和含水率影响下暂态饱和区及地下水位对比F-g．9Transientsaturatedareasvs．groundwaterlevelunderinfluenceofvarioussaturatedwatercontents4结论a．q／k是决定路堤内部雨水人渗方式及暂态饱和区形成的关键因素。当q／k≥1．5时，在路堤边坡坡面逐渐会出现一定深度的暂态饱和区，随着降雨的持续，暂态饱和区在坡脚附近率先与地下水位线联通，并引起路堤边坡坡面以下水位出现大幅度上升；当q／k<1．5时，则是通过雨水的持续人渗直接引起坡脚附近地下水位线上升。b．降雨条件下，坡脚附近地下水最易受雨水人渗的影响，坡脚以下土体最先趋于饱和，饱和区范围由坡脚附近迅速向路堤截面横向扩散，并逐渐覆盖坡脚水平面以下区域，同时由坡脚处沿着坡面向上延伸。在整个降雨阶段，雨水入渗仅会引起路堤边坡坡面以下水位线大幅度上升，对路堤中部地下水位影响有限。c．降雨强度、土体饱和渗透系数以及土体饱和含水率对雨水人渗过程的影响程度不同。路堤边坡暂态饱和区范围的扩展速度主要取决于降雨强度与饱和含水率的大小，暂态饱和区出现的时间与扩展范围主要受降雨强度、降雨历时与饱和含水率的共同影响，受饱和渗透系数的变化影响较小。在路堤的设计与施工过程中，除充分考虑当地的气候条件与路堤填料对边坡稳定性的影响以外，路堤边坡与坡脚附近水平面的防排水措施不可忽视。参考文献：[1]徐海涛，梁超．降雨条件下顺层边坡稳定性分析[J]．水利水电科技进展，2013，33(5)：73-76．(XUHaitao，LIANGChao．Stabilityanalysisofbeddedslopesunderrainfallconditions[J]．AdvancesinScienceandTechnologyofWaterResources，2013，33(5)：73—76．(inChinese))[2]李亮，刘兴旺，赵炼恒．降雨入渗对路基稳定性影响因素分析[J]．铁道科学与工程学报，2007，4(2)：19-23．(LILing，LIU 河海大学学报(自然科学版)第42卷Xingwang，ZHAOLianheng．Factorsanalysisofsubgradestabilityduetorainfallinfihration[J]．JournalofRailwayScienceandEngineering，2007，4(2)：19-23．(inChinese))[3]刘俊新，刘育田，胡启军．非饱和地表径流一渗流和流固体耦合条件下降雨入渗对路堤边坡稳定性研究[J]．岩土力学，2010，31(3)：903-910．(LIUJunxin，LIUYutian，HUQijun．Stabilityofembankmentslopesubjectedtorainfallinfiltrationconsideringbothrunof-undergroundseepageandfluid—solidcoupling[J]．RockandSoilMechanics，2010，31(3)：903-910．(inChinese))[4]沈水进，孙红月，尚岳全，等．降雨作用下路堤边坡的冲刷一渗透耦合分析[J]．岩石力学与工程学报，201l，30(12)：2456—2462．(SHENShuijin，SUNHongyue，SHANGYuequan，eta1．Scouring—penetrationcouplinganalysisofembankmentslopeunderrainfallaction[J]．ChineseJournalofRockMechanicsandEngineering，2011，30(12)：2456-2462．(inChinese))[5]陈晓斌，徐望国，刘小平．降雨入渗对粗粒土路堤变形与稳定性的影响[J]．中南大学学报：自然科学版，2011，42(3)：766—771．(CHENXiaobin，XUWangguo，LIUXiaoping．Effectofrainfallinfiltrationonadditionalsettlementandstabilitydebasementofgranularsoilfillingsembankment[J]．JournalofCentralSouthUniversity：ScienceandTechnology，2011，42(3)：766-771．(inChinese))[6]李佳，高广运，黄雪峰．非饱和原状黄土边坡浸水试验研究[J]．岩石力学与工程学报，2011，30(5)：1043-1048．(LIJia，GAOGuangyun，HUANGXuefeng．Experimentalresearchonimmersionforunsaturatedintactloessslope[J]．ChineseJournalofRockMechanicsandEngineering，2011，30(5)：1043—1048．(inChinese))[7]付宏渊，曾铃，蒋中明，等．降雨条件下公路边坡暂态饱和区发展规律[J]．中国公路学报，2012，25(3)：59-64．(FUHongyuan，ZENGLing，JIANGZhongming，eta1．Developinglawoftransientsaturatedareasofhighwayslopeunderrainfallconditions[J]．ChinaJournalofHighwayandTransport，2012，25(3)：59-64．(inChinese))[8]邹维列，李聪，汪建峰，等．非饱和路堤对加载和降雨入渗响应的模型试验研究[J]．岩土工程学报，2009，31(10)：1512—1519．(ZOUWeilie，LICong，WANGJianfeng，eta1．Modeltestsonresponsesofunsaturatedroad·embankmentstoloadingandrainfallinfiltration[J]．ChineseJournalofGeoteehnicalEngineering，2009，31(10)：1512—1519．(inChinese))[9]TASITL，CHENHF．Effectsofsaturationonshallowlandslidestriggeredbyrainfall[J]．EnvironmentalEarthSciences，2010，59(6)：1285—1295．[10]OHWT，VANAPALLISK．Influenceofraininfiltrationonthestabilityofcompactedsoilslopes[J]．ComputersandGeotechnics，2010，37(5)：649-657．[11]刘新喜，夏元友，蔡俊杰，等．降雨入渗下强风化软岩高填方路堤边坡稳定性研究[J]．岩土力学，2007，28(8)：1705—1709．(LIUXinxi，XIAYuanyou，CAIJunjie，eta1．Studyonstabilityonhighfilledembankmentslopeofhighlyweatheredsoftrockunderrainfallinfiltration[J]．RockandSoilMechanics，2007，28(8)：1705—1708．(inChinese))[12]王瑞钢，闫澍旺，邓卫东．降雨作用下高填土路堤边坡的渗流稳定分析[J]．中国公路学报，2004，17(4)：25—30．(WANGRuigang，YANShuwang，DENGWeidong．Analysisofseepagestabilityofhill—filledembankmentslopeduetorainfallinfiltration[J]．ChinaJournalofHighwayandTransport，2004，17(4)：25—30．(inChinese))[13]张蔚榛．地下水与土壤水动力学[M]．北京：中国水利水电出版社，1996．[14]SHAHROKHABADIS，TOUFIGHMM．Thesolutionofunconfinedseepageproblemusingnaturalelementmethod(NEM)coupledwithgeneticalgorithm(GA)[J]．AppliedMathematicalModelling，2013，31：2775-2786．[15]GREENWH，AMPTGA．Studiesonphysics：flowofairandwaterthroughsoils[J]．JAgrSci，1911(4)：1—24[16]RUSHTONKR，YOUNGSEG．Drainageofrechargetosymmetricallylocateddownstreamboundarieswithspecialreferencetoseepagefaces[J]．JournalofHydrology，2010，380：94—103．[17]雷志栋，杨诗秀，谢森传．土壤水动力学[M]．北京：清华大学出版社，1998．[18]vanGENUCHTENMT．Aclosedformequationforpredictionthehydraulicconductivityofunsaturatedsoils[J]．SoilSciSocAmJ，1980，44：892—898．[19]刘建华，查旭东，付宏渊，等．考虑降雨入渗条件下岩质边坡稳定性分析[J]．公路交通科技，2009，26(10)：33—37．(LIUJianhua，ZHAXudong，FUHongyuan，eta1．Analysisofrockslopestabilityunderrainfallinfiltrationcondition[J]．JournalofHighwayandTransportationandResearchandDevelopment，2009，26(10)：33—37．(inChinese))