降雨工况下渗透力对粗粒土高路堤边坡稳定性的影响.pdf

第3４卷第9期公路交通科技Vol．3４No．92017年9月JournalofＨighwayandTransportationResearchandDevelopmentSep．2017doi:10．3969/j．issn．1002-0268．2017．08．003降雨工况下渗透力对粗粒土高路堤边坡稳定性的影响1，2223何忠明，邓喜，付宏渊，曾铃(1．长沙理工大学特殊环境道路工程湖南省重点实验室，湖南长沙４1011４;2．长沙理工大学交通运输工程学院，湖南长沙４1011４;3．长沙理工大学土木与建筑学院，湖南长沙４1011４)摘要:为分析降雨工况下,渗透力对粗粒土高路堤边坡塑性区分布以及安全系数变化规律的影响,依托现场工程实例,以粗粒土骨架为研究对象,采用geo-studio软件中的seep模块计算边坡降雨过程中及降雨停止后水位线、基质吸力、渗流梯度的变化规律;利用FＬAC3D内置的fish语言将渗流计算结果导入到FＬAC3D中;综合考虑水位线、基质吸力变化,分析了渗透力变化对粗粒土高路堤边坡安全系数以及塑性区分布规律的影响。数值计算结果表明:降雨影响水位线、基质吸力、渗流梯度在边坡空间上的分布;降雨过程中水位线逐渐由坡脚抬升并向粗粒土高路堤边坡内部扩展,降雨停止后坡脚处水位线下降较快且路堤中心水位线高于初始状态水位线高度;水位线抬升过程中基质吸力逐渐消散,降雨停止后基质吸力恢复较缓慢;降雨初期坡脚处x方向渗流梯度逐渐由正值转为负值且负值分布的区域逐渐扩大;随降雨历时增加边坡塑性区由坡脚不断向边坡内部逐渐延伸,降雨停止后边坡塑性区主要发生在坡脚处;降雨0~2４h内边坡安全系数呈逐渐增大的趋势,降雨2４~72h内安全系数下降较快并达到最小值,降雨停止后安全系数逐渐恢复,但安全系数恢复过程具有一定滞后性并且恢复后的数值小于初始值。关键词:道路工程;降雨入渗;FＬAC3D;粗粒土高路堤;渗透力中图分类号:U４16．21文献标识码:A文章编号:1002-0268(2017)09-001４-08InfluenceofSeepageForceonStabilityofCoarseGrainedSoilHighEmbankmentunderRainfallCondition1，2223ＨEZhong-ming，DENGXi，FUＨong-yuan，ZENGＬing(1．ＨunanProvincialKeyＬaboratoryofSpecialEnvironmentRoadEngineering，ChangshaUniversityofScience&Technology，ChangshaＨunan４1011４，China;2．SchoolofTransportationEngineering，ChangshaUniversityofScienceandTechnology，ChangshaＨunan４1011４，China;3．SchoolofCivilEngineeringandArchitectural，ChangshaUniversityofScience&Technology，ChangshaＨunan４1011４，China)Abstract:Inordertoanalyzetheeffectofseepageforceonthedistributionofplasticzoneincoarsegrainedsoilhighembankmentslopeandthechangeofsafetyfactorunderrainfallcondition，basedonfieldengineeringcase，takingframeworkofcoarsegrainedsoilastheresearchobject，thechangesofmatrixsuction，undergroundwatertable，andseepagegradientoftheslopeduringandafterrainfallarecalculatedbyusingtheseepmoduleofgeo-studiosoftware．Then，thecalculationresultofseepageisimputedintoFＬAC3DbyusingfishlanguageofFＬAC3D．Consideringthechangeofundergroundwatertablelineandmatrixsuctioncomprehensively，theinfluenceofseepageforcechangeonthesafetycoefficientofcoarsegrainedsoilhighembankmentslopeandthedistributionofplasticzoneisanalyzed．Thenumericalcalculationresultshowsthat收稿日期:2016-10-2４基金项目:国家自然科学基金项目(515080４2，515080４0，51678073，51578079);特殊环境道路工程湖南省重点实验室开放基金项目(kfj1４0501);湖南省重点研发计划项目(2016SK2023)作者简介:何忠明(1980-)，男，湖南永兴人，博士，教授．(hezhongming４5@126．com) 第9期何忠明，等:降雨工况下渗透力对粗粒土高路堤边坡稳定性的影响15(1)rainfallaffectedthespatialdistributionsofundergroundwatertable，matrixsuctionandtheseepagegradientofslope;(2)theundergroundwatertableliftedgraduallyfromthefootofslopeanddevelopedintothedeepoftheslopeduringrainfall，itdecreasedfastandtheheightoffinalundergroundwatertablenearembankmentcenterishigherthantheinitialstateafterrainfall;(3)thematrixsuctiondissipatedgraduallyduringtheliftofundergroundwatertable，butitrecoveredslowlywhenrainfallstopped;(４)atthebeginningofrainfall，thepositivevaluesofseepagegradientinxdirectionturnedtonegativeandtheareasofnegativeexpandgradually;(5)withthetimeofrainfallincreasing，theplasticzoneofslopeextendedgraduallyfromthefoottotheinsideofslope，andtheplasticzonemainlyoccurredatthefootofslopeafterrainfall;(6)thesafetyfactorofslopeincreasedslowerattheinitial2４h，decreasedfastandreachedtheminimumwhileitlastedduring2４-72h，andgraduallyrestoredafterrainfall，buttherecoveryofsafetyfactorisdelayedandsmallerthantheinitialstate．Keywords:roadengineering;rainfallinfiltration;FＬAC3D;coarsegrainedsoilhighembankment;seepageforce用FＬAC3D中的fish语言将渗透力以体力的方式施0引言加在边坡各节点上，分析水位线升降、基质吸力变粗粒土(粒径在0．075mm≤d≤60mm范围内化、重度变化影响下，渗透力对粗粒土高路堤边坡0的质量占试样总质量的50%以上[1])因其具有孔隙安全系数变化以及塑性区分布规律的影响，可为粗率大、透水性强、抗压强度高、压实密度大等工程粒土高路堤边坡设计提供指导。特性，在湖南高速公路高路堤的填筑施工中已得到1边坡渗流特性及稳定性分析理论[2]广泛应用。但是根据陈晓斌的调研结果可知，在湖南多雨地区，已运营的高速公路粗粒土高路堤边1．1饱和-非饱和渗流理论坡受降雨影响灾变明显。近年来广大学者针对降雨粗粒土高路堤边坡的渗流问题，满足达西定律土体内非恒定渗流，其偏微分方程形式如式(1)工况下，渗透力致使边坡失稳问题，做了一些研究:[10][3]所示:如李广信指出了土骨架和渗透力的定义，为降雨∂∂h∂∂h∂h入渗过程粗粒土骨架受力分析提供了理论依据。陈(kx)+(ky)+Q=mwρwg，(1)∂x∂x∂y∂y∂t[４]祖煜提出并论证了在计算边坡稳定性时，采用水、式中，h为土体中的总水头;k和k为x和y方向的xy土合算，水、土分算两种方法得到的结果较为一致，渗透系数;Q为源汇项;m为比水容量;ρ为水的ww也从侧面反映了以粗粒土骨架为研究对象来分析渗密度;g为重力加速度;t为时间。透力对边坡稳定性影响方法的可行性。连继峰，刘[11]非恒定渗流的初始条件为:[5-8]子振，付宏渊等分析了渗透力、有效应力作用h(x，y，0)=h，(x，y)∈Ω。(2)0下，降雨入渗对边坡稳定性影响。宋子亨，马年有限元方程求解可结合以下边界条件求解:[9-10]祖研究发现基质吸力的存在对土体抗剪强度具h(x，y，t)=h，(x，y)∈S，(3)11有一定提高。综上所述，现有学者关于渗透力对边∂hk|S=q(x，y，t)，(x，y)∈S，(４)22坡稳定性影响的有限元分析方面尚存在一些不足:∂n一方面，重复考虑了孔隙水压力变化和渗透力对边式中，S1为已知水头边界;S2为已知流量边界。坡稳定性的影响;另一方面，单一考虑静孔隙水压在缺乏试验成果的情况下，可采用经典的Van力、基质吸力变化对边坡稳定性的影响，不能完全Genuchten模型，拟合粗粒土高路堤边坡的土水特征[12]反映雨水在边坡内部入渗过程中，渗透力作用方向曲线。模型表达式为:对边坡稳定性的影响，尤其是对于透水性强的粗粒θs-θrθ=θ+，(5)rnm土，渗透力对粗粒土骨架力学性能影响更为明显。[1+(α·h)]11m2因此，针对湖南地区的气候特点，以粗粒土骨æθ-θrö2æθ-θrömK(θ)=Ks{1-[1-]}，架为研究对象，基于饱和-非饱和理论，采用seepèθs-θrøèθs-θrø模拟粗粒土高路堤边坡降雨工况下的渗流特性;应(6) 16公路交通科技第3４卷式中，θ为土壤体积含水率;h为压力水头;θ和θτ=c+σ′tanφ，(9)srf分别为粗粒土饱和体积含水量与残余体积含水量;α式中，c为黏聚力;σ′为剪切滑裂面的有效正应力;和n为经验拟合参数(或曲线形状参数)，而m=φ为内摩擦角。由于式(9)仅考虑了饱和状态下土1-1/n;Ks为饱和渗透系数;K(θ)为渗透系数随坡体的抗剪强度，未考虑非饱和状态下基质吸力对边[15]体内体积含水率的变化关系。坡稳定性的贡献，据此Fredlund等根据非饱和土1．2渗透力计算的双应力变量(σ-μ，μ-μ)理论提出了非饱和aaw渗透力是指粗粒土高路堤边坡在降雨入渗过程抗剪强度公式:中，雨水在坡体内部流动对粗粒土骨架所施加的力;bτ=c+(σ-μ)tanφ+(μ-μ)tanφ，faaw[3]渗透力除了拖曳力外，还有沿流向的巨大压力差。(10)因此仅通过研究降雨入渗前后静孔隙水压力变化，式中，μ为边坡土体的进气值，由于考虑边坡孔隙a不能真正反映雨水入渗过程中，渗透力对边坡稳定b与大气相连，μ=0;φ为基质吸力对内摩擦角的贡a性的影响。降雨入渗过程中雨水对粗粒土骨架作用b献，按经验取值φ=1４°;μ为孔隙水压力。w效果如图1所示。1．4边坡稳定性分析方法在边坡稳定性分析过程中，由于FＬAC3D软件在边坡渗流计算时，软件会将负孔隙水压力全部归零，故在使用FＬAC3D软件分析边坡流固耦合问题时，不能考虑基质吸力对粗粒土高路堤边坡稳定性的影响。因此，为充分反映雨水入渗过程中渗透力对边坡稳定性的影响，以粗粒土骨架为研究对象，采用geo-studio软件的seep模块计算边坡的渗流特性，然后将渗透力导入到FＬAC3D中计算粗粒土高路堤边坡稳定性。根据渗透力定义可知，渗透力是由水头损失产生，故当以粗粒土骨架为研究对象时，可仅考虑渗透力对边坡稳定性影响，不再重复考虑正孔隙水压力变化对边坡稳定性的影响。当以粗粒土骨架为研图1渗透力作用示意图究对象，粗粒土骨架的静力平衡方程如式(11)、Fig．1Schematicdiagramofseepageforce[13](12)所示。边坡稳定性计算的具体流程图如图2渗透力计算如式(7)、(8)表示。由于边所示。坡渗流过程中，内部渗流路径极为复杂且难以确定，G+G′+τ+F=0，(11)因此文中采用geo-studio软件中的seep模块，模拟粗abG′=(γ-γ)V，(12)粒土高路堤边坡渗流过程中x，y两方向的渗流梯bsw式中，G为粗粒土骨架水上部分自重;G′为粗粒土度;再将x，y两个方向渗流梯度乘以相应的水重度ab骨架水下部分浮重;τ为粗粒土骨架抗剪强度;F为(γ)，通过fish语言导入到FＬAC3D中，其中渗透w粗粒土骨架所受渗透力合力;γ为粗粒土饱和重度;力的方向为x，y两个力合力的方向。sγ为水的重度;V为水位线以下粗粒土体积。F=γ·j，(7)ww→→→F=F+F，(8)2边坡计算条件xy式中，F为渗透力;γ为水的重度;j为水力梯度;w2．1计算模型→→Fx，Fy分别为各节点渗透力在x，y两个方向上的粗粒土高路堤边坡依托工程位于湖南省郴州地分力。区，起点桩号为K10+320，边坡高度为20m。边坡1．3非饱和抗剪强度理论分为两级，第1，2级边坡高度分别为8，12m，坡法国著名科学家库仑于1776年根据抗剪强度试率分别为1∶1．5，1∶1．75，降雨入渗过程如图3所示。[1４]验提出土的抗剪强度τf计算公式:利用geo-studio软件中的seep模块建立相应的计算模 第9期何忠明，等:降雨工况下渗透力对粗粒土高路堤边坡稳定性的影响17残余体积含水率为0．08，饱和渗透系数为1．02×-610m/s，并假设粗粒土高路堤填料为各向同性材料。由于通过室内试验测定粗粒土的土水特征曲线具有较大困难，故采用VanGenuchten模型拟合粗粒土渗透系数、体积含水率随基质吸力的变化规律;根据现有研究成果，粗粒土路堤边坡最大基质吸力不超过100kPa，因此渗流分析过程中设置初始最大[16]基质吸力为60kPa。拟合后的VanGenuchte关系图2路堤边坡稳定性计算流程曲线如图5、图6所示。Fig．2Calculationprocessofembankmentslopestability型，为方便观察边坡表层基质吸力变化规律，在边坡坡面以下2m布置5个监测点。计算模型如图４所示，总节点数为512４个，单元个数为3228个。图5x方向渗透系数随基质吸力变化关系Fig．5Relationshipbetweenpermeabilitycoefficientandmatrixsuctioninxdirection图3降雨入渗过程Fig．3Processofrainfallinfiltration图6含水率随基质吸力变化规律Fig．6Relationshipbetweenwatercontentandmatrixsuctioninxdirection通过室内筛分试验、查阅相关文献、物理力学参数试验可知，郴州地区粗粒土的级配曲线变化关图4路堤边坡有限元网格划分系如图7所示，计算选用的粗粒土高路堤边坡的相Fig．4Finiteelementmeshingofembankmentslope[17-18]关参数如表1所示。边界条件:坡脚左侧以及坡面为降雨入渗边界，边坡坡顶为封闭道路，设置为不透水边界，边坡底部设置为不透水边界。2．2计算参数查阅郴州地区气象资料可知，郴州历史上最大连续3日降雨量为357mm。因此，文中考虑极端降雨工况下粗粒土高路堤边坡的稳定性，取降雨强度-6为1．38×10m/s，通过室内试验，郴州地区粗粒图7粗粒土级配曲线土高路堤在标准压实度下饱和体积含水率为0．12，Fig．7Gradingcurvesofcoarsegrainsoil 18公路交通科技第3４卷表1边坡物理力学参数体积含水率随基质吸力的增加逐渐减小，在坡脚处Tab．1Physico-mechanicalparametersofslope基质吸力较小，所以渗透系数较大，雨水入渗量较饱和重度干重度有效黏有效内变形b多，先达到饱和;随着降雨历时的增加，边坡坡面φ/泊松岩性γsat/γd/聚力c/摩擦角模量E/(°)比μ附近逐渐形成暂态饱和区，渗透系数明显大于道路(kN·m-3)(kN·m-3)kPaφ/(°)MPa中心线附近的渗透系数，所以靠近坡面处的水位线粗粒土22193532251４0．35抬升速率明显快于道路中心线附近的水位线抬升由图7可知，研究区内0．075mm≤d≤20mm0速率。约占总量的55%，黏土含量为４5%。3．2边坡降雨入渗基质吸力分布规律3计算结果分析降雨过程中，边坡坡面监测点的基质吸力变化规律如图9所示。3．1边坡降雨入渗水位线分布规律图8为粗粒土高路堤边坡降雨过程中和降雨停止后边坡的水位线变化规律。图9降雨过程中边坡监测点孔隙水压力分布规律Fig．9Distributionofporewaterpressureatmeasuringpointsduringrainfall由图9可知，随降雨历时增加边坡基质吸力(即负孔隙水压力)逐渐消散并转变为正孔隙水压力，降雨停止后边坡内部正孔隙水压力逐渐减小，基质吸力逐渐恢复;此外，初始条件下监测点距离边坡水位线越近，其孔隙水压力越大。产生上述现象的原因为:随着降雨历时增加，粗粒土高路堤边坡内部水位线抬升，边坡坡面附近形成暂态饱和区，孔隙被入渗的雨水填充，最终基质吸力消散，转变为正孔隙水压力并逐渐增大;降雨停止后，雨水从图8高路堤边坡地下水位线分布规律边坡内部排出，基质吸力逐渐恢复。Fig．8Distributionsofundergroundwatertableof3．3边坡降雨入渗渗流速度分布规律highembankment根据渗透力定义可知，雨水入渗过程中渗流速由图8可知，降雨过程中粗粒土高路堤边坡首度矢量大小及分布规律直接反映渗透力分布规律。因先在坡脚处到达饱和，并且随降雨历时增加，粗粒此，为揭示降雨过程中，入渗的雨水在粗粒土高路堤土高路堤边坡水位线逐渐抬升，当降雨历时60h时，边坡中的渗流规律以及渗透力的分布规律，给出降雨边坡水位线高度达到第1级平台坡脚处;此外，靠2４，72h的渗流速度矢量分布，如图10所示。近边坡坡面处的水位线抬升速度明显快于道路中心由图10可知，0~2４h内降雨入渗主要分布在线附近的水位线抬升速率;降雨停止后，边坡水位边坡坡脚以及边坡坡面表层附近，速率分布的面积线逐渐下降，并且靠近边坡坡面处的水位线下降速较小，降雨72h时雨水逐渐向边坡内部转移且降雨率明显快于远离边坡坡面处的水位线下降速率。产过程中部分雨水速度矢量指向坡面。产生上述现象生上述现象的原因为:粗粒土高路堤边坡采用Van的原因为:降雨0~2４h内坡脚先达到饱和，雨水主Genuchten模型拟合渗透系数、体积含水率随基质吸要汇聚在坡脚处，渗透系数较大，所以速度矢量主力的变化关系，在VanGenuchten模型中渗透系数和要分布在坡脚及边坡坡面表层附近;随降雨历时的 第9期何忠明，等:降雨工况下渗透力对粗粒土高路堤边坡稳定性的影响19述现象的原因为:边坡x方向渗流等值线以流入边坡为正，流出边坡为负。由图10(a)可知，降雨初期雨水主要以入渗到边坡内部为主，仅有少量雨水从坡脚处排出，故仅在坡脚附近较小区域出现x方向渗流梯度为负值。由图10(b)可知，当降雨72h时，暂态饱和区逐渐向边坡内部扩展，渗透系数达到饱和渗透系数，雨水持续入渗到边坡内部，故靠近边坡内部x方向渗流梯度分布的区域逐渐扩大;同时，坡脚附近的暂态饱和区在重力作用下，入渗的雨水也不断从坡脚处排出，故坡脚处x方向渗流梯度负值分布区域逐渐扩大。3．5导入渗流计算结果后的三维模型将seep计算的渗流结果转化为渗透力，通过fish将基质吸力、重度修正、渗透力导入到FＬAC3D中。受图10速度矢量分布规律篇幅限制，仅给出了降雨72h对应的基质吸力、重度、Fig．10Distributionofvelocityvector渗透力导入到FＬAC3D模型中的示意图，如图12~图1４增加，水位线逐渐抬升，暂态饱和区的面积逐渐增所示。其中，考虑水位线以下浮力的影响，水位线以下大，渗透系数逐渐增大，此外部分入渗的雨水在重粗粒土的重度采用有效重度。力作用下从坡脚处排出，所以降雨72h时，边坡速度矢量在边坡内部分布的面积较广，部分速度矢量箭头指向坡面。3．4边坡降雨入渗x方向渗流梯度分布规律由图11可知，降雨36h时，x方向渗流梯度等值线仅在坡脚处局部区域出现渗流梯度为负值，其余区域渗流梯度均为正值;降雨72h时，靠近边坡内部的x方向渗流梯度等值线分布区域逐渐扩大，坡脚附近的x方向渗流梯度负值分布区域逐渐扩大，且该区域内的渗流梯度的绝对值逐渐增大。产生上图12导入基质吸力后的三维模型(单位:kPa)Fig．12Three-dimensionalmodelafterimputingmatrixsuction(unit:kPa)-3图13修改密度后的三维模型(单位:kg·m)Fig．13Three-dimensionalmodelafterrevising-3densitty(unit:kg·m)3．6边坡降雨入渗塑性区分布规律为真实反映降雨过程中及降雨停止后塑性区发展规律，图15中塑性区为降雨影响新产生的塑性区，以前图11x方向渗流梯度积累的塑性区未显示。考虑篇幅有限仅列出了降雨工况Fig．11Seepagegradientsinxdirection下，降雨36，72h以及降雨停止后12，４8h的边坡塑 20公路交通科技第3４卷1级路堤边坡坡脚处并不断向第1级边坡坡脚附近延伸，当降雨历时72h时，第2级边坡内形成较大面积的塑性贯通区;降雨停止后，塑性区的面积逐渐减小且发生在坡脚处。分析其原因为:降雨历时36h时边坡水位线从坡脚处逐渐抬升，致使边坡坡脚处粗粒土重度增加，基质吸力减小，导致粗粒土的有效黏聚力减小;此外，边坡降雨过程中雨水首先在坡脚处汇聚，雨水汇聚过程中形成的渗透力使得边坡坡脚处易于发生塑性变形。降雨历时72h，暂态饱和区不断向边坡内部扩展，x方向渗流梯度负值分布区域逐渐扩大，沿滑移方向的渗透力增大，即下滑力增大，故塑性区不断向边坡内部扩展。3．7边坡降雨入渗安全系数分布规律降雨工况下，计算得到的粗粒土高路堤边坡安图14导入渗透力后的三维模型全系数变化规律如图16所示。Fig．14Three-dimensionalmodelafterimpuingseepageforce性区发展变化分布图，如图15所示。图16降雨过程中高路堤安全系数变化规律Fig．16Distributionofsafetyfactorofhighembankmentduringrainfall由图16可知，降雨0~2４h内边坡安全系数先升高到1．8，降雨2４~72h内，边坡安全系数迅速下降，从1．8降到1．４;降雨停止12h内，边坡安全系数持续降低，达到最低值1．31;降雨停止12~４8h内，边坡安全系数逐渐恢复到1．４。究其原因为:降雨0~2４h，地下水位线抬升较少，基质吸力减小较少并且此段时间内，雨水主要从坡脚处入渗到边坡内部，渗流产生的渗透力对边坡具有加固作用，所以安全系数先逐渐上升;2４~72h内，边坡坡面附近形成暂态饱和区，粗粒土渗透能力变强，雨水入渗量增大，水位线抬升速度较快，土体达到饱和，土体重度增大，边坡下滑力增大，基质吸力消散较快，导致粗粒土有效黏聚力减小并且雨水不断从坡脚处排出，渗流梯度负值分布的区域逐渐扩图15高路堤塑性区分布规律大，渗透力方向沿着边坡滑移方向，增大边坡下滑Fig．15Distributionofplasticzoneofhighembankment力，边坡塑性区向边坡内部发展，致使边坡安全系由图15可知，降雨过程中塑性区首先发生在第数下降的速率较快;降雨停止0~12h后，雨水逐渐 第9期何忠明，等:降雨工况下渗透力对粗粒土高路堤边坡稳定性的影响21从边坡坡脚排出，渗流梯度为负值，渗透力方向沿Beijing:ChinaWater&PowerPress，2003．着边坡滑移的方向，致使下滑力增大，边坡安全系[5]连继峰，罗强，蒋良潍，等．顺坡渗流条件下土质边坡浅层稳定分析[J]．岩土工程学报，2015，37数持续降低;降雨停止12~４8h后，水位线下降，(8):1４４0-1４４8．基质吸力逐渐恢复，边坡安全系数逐渐恢复。ＬIANJi-feng，ＬUOQiang，JIANGＬiang-wei，etal．4结论ShallowStabilityAnalysisofSoilSlopesunderSeepageParalleltoSlopeSurface[J]．ChineseJournalof(1)随降雨历时的增加，粗粒土高路堤边坡水GeotechnicalEngineering，2015，37(8):1４４0-1４４8．位线逐渐抬升，靠近边坡坡面的水位线抬升的速率[6]刘子振，言志信．降雨入渗非饱和黏土边坡稳定性的明显快于道路中心线附近的水位线抬升速率;降雨极限平衡条分法研究[J]．岩土力学，2016，37停止后，水位线逐渐下降。(2):350-356．(2)降雨过程中随着水位线的抬升，粗粒土高ＬIUZi-zhen，YANZhi-xin．ＬimitEquilibriumSlice路堤边坡基质吸力逐渐消散并逐渐转变为正孔隙水MethodforUnsaturatedClaySlopeunderRainfallInfiltration[J]．RockandSoilMechanics，2016，37压力，降雨停止后基质吸力逐渐恢复。(2):350-356．(3)随降雨历时增加，粗粒土高路堤边坡x[7]汪丁建，唐辉明，李长冬，等．强降雨作用下堆积层滑方向渗流梯度在坡脚处逐渐出现负值且负值区域坡稳定性分析[J]．岩土力学，2016，37(2):４39-逐渐扩大，因此考虑渗透力来计算边坡稳定性，４４5．能更真实地反映雨水入渗过程中边坡稳定性的变WANGDing-jian，TANGＨui-ming，ＬIChang-dong，et化规律。al．StabilityAnalysisofColluvialＬandslideduetoＨeavy(４)初始条件下粗粒土高路堤边坡塑性区主要Rainfall[J]．RockandSoilMechanics，2016，37(2):发生在坡脚处，随降雨历时的增加，边坡塑性区逐４39-４４5．渐向边坡内部发展;边坡安全系数先升高后逐渐降[8]付宏渊，邱祥，李光裕，等．降雨入渗条件下炭质泥低;降雨停止后边坡塑性区主要发生在坡脚处，安岩路堤动态稳定性[J]．长安大学学报:自然科学版，2017，37(1):33-４2．全系数逐渐恢复，但有一定滞后性。FUＨong-yuan，QIUXiang，ＬIGuang-yu，etal．Dynamic参考文献:StabilityofCarbonaceousMudstoneEmbankmentunderReferences:RainfallInfiltration[J]．JournalofChang’anUniversity:NaturalScienceEdition，2017，37(1):33-４2．[1]JTGD30—2015，公路路基设计规范[S]．[9]宋子亨，杨强，刘耀儒．考虑孔隙水压力作用的岩土JTGD30—2015，SpecificationsforDesignofＨighway体弹塑性模型及其有限元实现[J]．岩土力学，Subgrades[S]．2016，37(增1):500-508．[2]陈晓斌，徐望国，刘小平．降雨入渗对粗粒土路堤变SONGZi-heng，YANGQiang，ＬIUYao-ru．Elastoplastic形与稳定性的影响[J]．中南大学学报:自然科学ModelforGeomaterialConsideringEffectofPoreWater版，2011，４2(3):765-771．PressureandItsFiniteElementsImplementation[J]．CＨENXiao-bin，XUWang-guo，ＬIUXiao-ping．EffectofRockandSoilMechanics，2016，37(S1):500-508．RainfallInfiltrationonAdditionalSettlementandStability[10]马年祖，建德琳，代玉莉，等．山岭重丘区道路改扩DebasementofGranularSoilFillingsEmbankment[J]．建路基持续降雨渗流数值模拟[J]．公路交通科技，JournalofCentralSouthUniversity:NaturalScience2016，33(9):31-37，４5．Edition，2011，４2(3):765-771．MANian-zu，JIANDe-lin，DAIYu-li，etal．Numerical[3]李广信．论土骨架与渗透力[J]．岩土工程学报，SimulationofSubgradePersistentRainfallSeepageof2016，38(8):1522-1528．ReconstructedＨighwayinMountainousandＨeavy-hillyＬIGuang-xin．OnSoilSkeletonandSeepageForce[J]．Area[J]．JournalofＨighwayandTransportationChineseJournalofGeotechnicalEngineering，2016，38ResearchandDevelopment，2016，33(9):31-(8):1522-1528．37，４5．[４]陈祖煜．土质边坡稳定性分析[M]．北京:中国水利[11]付宏渊，曾铃，王桂尧．降雨入渗条件下软岩边坡稳定水电出版社，2003．性分析[J]．岩土力学，2012，33(8):2359-2365．CＨENZu-yu．StabilityAnalysisofSoilSlopes[M]．(下转第29页) 第9期白璐，等:柔性基层沥青路面在动荷载作用下Top-Down开裂模拟研究29SＨAOXian-zhi，SＨAOMin-hua，BIYu-feng，etal．AsphaltPavementunderMeasuredTireGroundPressureofTestingMethodofAsphaltMixturePoissonRatio[J]．ＨeavyVehicle[J]．JournalofTongjiUniversity:NaturalJournalofTongjiUniversity:NaturalScienceEdition，ScienceEdition，2006，3４(1):6４-68．2006，3４(11):1４70-1４7４．[17]DEBEERM，FISＨERC．ContactStressesofPneumatic[15]庄传仪．基于加速加载响应的柔性基层沥青路面设计TiresMeasuredwiththeVehicle-roadSurfacePressure指标与参数研究[D]．西安:长安大学，2012．TransducerArray(VRSPTA)SystemfortheUniversityofZＨUANGChuan-yi．StructuralDesignIndexesandCaliforniaatBerkeley(UCB)andtheNevadaAutomotiveParametersforFlexibleBaseAsphaltPavementBasedonTestCenter(NATC)CR-97/053[R]．Pretoria:AcceleratedＬoadingResponse[D]．Xi’an:Chang’anDivisionofRoadsandTransportTechnology，1997．UniversityPress，2012．[18]OZERＨ，Al-QADIIＬ，DUARTECA．AThree-[16]胡小弟，孙立军．实测重型货车轮载作用下沥青路面dimensionalGeneralisedFiniteElementAnalysisforthe力学响应[J]．同济大学学报:自然科学版，2006，Near-surfaceCrackingProbleminFlexiblePavements3４(1):6４-68．[J]．InternationalJournalofPavementEngineering，ＨUXiao-di，SUNＬi-jun．StressResponseAnalysisof2011，12(４):４07-４19．■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■(上接第21页)2015，36(9):2４78-2４86，2４93．FUＨong-yuan，ZENGＬing，WANGGui-yao．Stability[15]顾慰慈．渗流计算原理及应用[M]．北京:中国建材AnalysisofSoftRockSlopeunderRainfallInfiltration工业出版社，2000．[J]．RockandSoilMechanics，2012，33(8):2359-GUWei-ci．ThePrincipleandApplicationofSeepage2365．Calculation[M]．Beijing:ChinaBuildingMaterials[12]杨金，简文星，杨虎锋，等．三峡库区黄土坡滑坡浸IndustryPress，2000．润线动态变化规律研究[J]．岩土力学，2012，33[16]FREDＬUNDDG，MORGENSTERNNR，WIDGERRA．(3):853-858．TheShearStrengthofUnsaturatedSoils[J]．CanadianYANGJin，JIANWen-xing，YANGＨu-feng，etal．GeotechnicalJournal，1978，15(3):313-321．DynamicVariationRuleofPhreaticＬineinＨuangtupo[17]宋焕宇．粗粒土斜坡高路堤变形性状与稳定性研究ＬandslideinThreeGorgesReservoirArea[J]．Rockand[D]．武汉:华中科技大学，2007．SoilMechanics，2012，33(3):853-858．SONGＨuan-yu．StudyonDeformationCharacteristicsand[13]魏义长，刘作新，康玲玲，等．土壤持水曲线vanStabilityofＨighEmbankmentFilledbyCoarseGrainedGenuchten模型求参的Matlab实现[J]．土壤学报，SoilsonSlopingGround[D]．Wuhan:Ｈuazhong200４，４1(3):380-386．UniversityofScienceandTechnology，2007．WEIYi-chang，ＬIUZuo-xin，KANGＬing-ling，etal．[18]曾铃，付宏渊，何忠明，等．饱和-非饱和渗流条件ParametersEstimationofVanGenuchtenModelforSoil下降雨对粗粒土路堤边坡稳定性的影响[J]．中南大WaterRetentionCurvesUsingMatlab[J]．ACTA学学报:自然科学版，201４，４5(10):361４-PedologicaSinica，200４，４1(3):380-386．3620．[1４]蒋中明，龙芳，熊小虎，等．边坡稳定性分析中的渗ZENGＬing，FUＨong-yuan，ＨEZhong-ming，etal．透力计算方法考证[J]．岩土力学，2015，36(9):ImpactofRainfallonStabilityofGranularSoil2４78-2４86，2４93．EmbankmentSlopeConsideringSaturated-unsaturatedJIANGZhong-ming，ＬONGFang，XIONGXiao-hu，etal．Seepage[J]．JournalofCentralSouthUniversity:StudyofCalculationMethodsofActingForceofSeepageinScienceandTechnologyEdition，201４，４5(10):SlopeStabilityAnalysis[J]．RockandSoilMechanics，361４-3620．