高填方加筋黄土路堤稳定性的有限元分析

第27卷第2期地理科学Vol.27No.22007年04月SCIENTIAGEOGRAPHICASINICAApr.,2007高填方加筋黄土路堤稳定性的有限元分析111,2王家鼎,谢婉丽,骆凤涛(1.西北大学地质系大陆动力学国家重点实验室,陕西西安710069;2.中交通力公路勘察设计工程有限公司,陕西西安710075)摘要:采用有限元法对国内最高(62m)的山西北张沟加筋黄土高填路堤稳定性进行了数值模拟,对其最大主应力、最小主应力、剪应力和塑性破坏区及其变形进行了计算,然后将其计算结果与未加筋路堤相比较,得出高填方加筋黄土路堤及地基变形特性和应力特性的变化规律。结果表明,由于加筋的作用,降低了黄土高路堤的变形,提高了其稳定性。这与现场沉降观测结果基本一致。关键词:加筋黄土;高填方路堤;土工格栅;有限元分析中图分类号:TU432文献标识码:A文章编号:1000-0690(2007)02-0268-05随着西部开发的不断深入,黄土高原地区公m以上的高路堤几乎尚未发现。另一方面,公路建路、铁路等基础建设取得了长足的发展,大大地推设中由于现有理论假设的局限,高路堤稳定分析方动了地方经济建设的可持续发展,特别是近年来黄法并不完全适用于冲沟地形条件,沉降变形分析与土高原地区高速公路、高速铁路的建设极大地缓解实际存在较大差别,这种理论研究滞后于工程实践了黄土高原的交通压力。然而,黄土高原地貌单元的现象在相当大的程度上限制了这种结构形式的十分复杂,沟壑纵横,加之黄土疏松,在外荷载作用推广应用。[1,2]下极易改变其性质,诱发次生灾害。为此,高基于以上原因,为研究高速公路是否可采用高速公路以何种方式穿过黄土沟谷,既经济,又安全,填加筋黄土路堤跨越黄土深沟,本文拟采用有限元是个值得研究的问题。法对山西北张沟62m的高填方加筋黄土路堤稳定加筋路堤是近几十年发展起来的介于挡墙与性进行数值模拟计算,并结合现场工程地质调查研一般填方边坡之间的一种新型加筋土结构,是在路究工作,将两者进行对比分析,为公路设计和施工堤和地基组成的土工结构物的适当位置加入具有提供基础资料。抗拉性能的土工合成材料组成的一种复合体。由1工程概况于土与加筋材料之间的摩擦作用,使加筋土形成了稳定的整体复合结构,从而提高了土体强度。土工北张沟高填方路堤位于山西省霍州市三教乡格栅作为一种筋材,其使用始于1980年。国内外村东约300m处的黄土深沟内,起终点桩号分别为[3～8]研究表明:由于土工格栅具有高强度、低延伸K102+008.095和K102+186.020,全长177.07率、高稳定性等特点,所以采用其加固高填路堤具m。该段原设计为一座45m+80m+45m连续刚有路堤稳定性好、边坡安全度高、沉降及差异沉降构桥,墩高33m,柱式台,桩基础。小和承载力强等诸多优点,因而在土建工程各行业1.1原设计方案变更理由得到广泛的应用,然而,高填方加筋黄土路堤稳定为了保证工期,节约工程费用,经实际踏勘认性的研究并不多见。在山区修建高等级公路时,为为:该处黄土沟形狭窄,沟底平缓,地质条件良好,满足公路线型的要求,路堤填筑高度的加大已愈来高路堤方案基本可行,我们提出改桥为高路堤方愈受到重视。如何改进和完善处理措施,确保路基案,原因如下:成形后的稳定,目前国内还处于研究探索阶段,601)该处地貌类型为黄土塬—冲积扇区,线路收稿日期:2005-09-23;修订日期:2006-01-11基金项目:国家自然科学基金项目(40572157,40502028)、高等学校博士学科点专项科研基金(20050697013)资助。作者简介:王家鼎(1962-),男,陕西蓝田人,博士,教授,博士生导师,从事工程地质研究。E2mail:wangjd@nwu.edu.cn 2期王家鼎等:高填方加筋黄土路堤稳定性的有限元分析269与冲沟斜交,设计高架桥通过时,桥长177m,而实接,平台向内设2%的反坡。边坡坡面与主沟两侧际冲沟顶垂直宽度仅为100m左右,沟底垂直宽度坡面交界处设急流槽;急流槽与路基两侧排水沟相仅为10～25m。接,直设计至沟底。急流槽为0.6m×0.6m的矩2)沟深62m,两侧山坡较陡,由于桥梁跨度形断面,其结构与主线相同;坡面采用拱形骨架护有限(如再增大跨径,将会造成边孔浪费),桥墩承坡进行防护,并植草进行绿化。台难于座落在两侧陡坡上,因此,必然导致高桥墩,4)沟底排水设计。沿沟底设盲沟,其尺寸为增加施工难度;桥墩施工平台形成时挖方量大;承2.0m×2.0m,以解决沟中及地下水的排泄。为防台上临空面过高(最高达45m),对桩基安全施工止盲沟堵塞,用土工布与土体隔离。带来一定隐患。2计算方案的选定3)桥位上游汇水面积很小,汇水总面积不超2过3km,采用高填方路堤,只设置1孔2.5m的泻用有限元的方法从理论上研究加筋黄土高路水隧洞即可满足排洪要求。堤,就是选用合适的物理力学模型来模拟加筋路4)桥位处土质较好,均为粘性黄土,是良好堤,然后将计算的结果与未加筋路堤相比较,从而的填筑路堤填料,而且北张沟北岸(祁县方向)约1探索加筋高路堤的工作状态和加筋机理。43km范围有20×10m黄土废方,可利用此黄土废2.1有限元方法进行加筋黄土高填路堤分析特点方填筑路基,以减少借方数量。有限元法是一种十分有效的数值分析方法。1.2高填方路堤设计方案进行加筋黄土高填路堤计算有如下几个特点:1)路基边坡坡率。距路面标高8m范围内(1)可以考虑黄土体的非线性弹塑性本构关的第一级边坡坡率采用1∶1.5,第一级平台宽度采系及变形对应力的影响;(2)能适应各种复杂的用2.0m。第二级边坡高度采用12m,坡率采用边界条件;(3)可以模拟路堤与地基、土体与加筋1∶1.75,第二级平台宽度采用3.0m。第三、四级材料之间的共同作用;(4)能够模拟高路堤的失边坡高度均采用12m,坡率采用1∶2.0,第三级平稳过程及其滑移面形状;(5)可以模拟高路堤施工台宽度采用4.0m,第四级平台宽度采用5.0m。的逐级加载过程。以下各级边坡高度、坡率及平台宽度与第四级相同2.2高填方加筋黄土路堤有限元分析过程(图1)。1)高填方加筋黄土路堤的几何参数与施工过程的模拟。本文选用与路线相交跨度最大的断面作为计算对象,即图1所示的高填方加筋路堤———地基系统。由于受冲沟地形的影响,路堤底部有1.5°的角度,路堤中心高度为62m,路堤顶宽为34.54m,路堤底宽为405.4m,地基深度约为一倍图1高填方加筋路堤示意图的路堤高度,平均深度为60m,地基宽度为路堤宽Fig.1Sketchofhighreinforcedembankment度向两边各延伸30m,即525.4m。路堤的施工步骤为逐层填土、铺筋、碾压,直到路堤顶面。有限元2)土工格栅的选择与布置。本次选用的土分析中,把一层填土当作是一次加荷,由此计算出工格栅主要材料是聚丙烯。其主要物理力学性质当前施工范围内格栅与土体内各处的变形和应力,2如下:单位面积质量500g/m;纵向尺寸100mm,上一层填土引起的变形和应力与先前值迭加,直至横向尺寸20mm;卷幅宽1.1m,卷长50m;最大抗整个结构完成。这样就实现了整个施工顺序的模拉强度30kN/m;延伸率10%。从高程700m开拟。始,止高程742m共铺设21层土工格栅,每2m设2)数值模型的建立。根据工程实际情况,上计一层(图1)。述高填方加筋土路堤所用材料沿道路轴向相对均3)排水及防护。每级平台上均设平台截水匀,而且路堤有足够的长度,故将高填方加筋黄土沟,尺寸为0.4m×0.4m,平台及水沟均用25cm路堤作为平面应变问题处理;对于该黄土冲沟地基厚浆砌片石进行加固防护,截水沟与两侧急流槽相和路基填土基本参数,我们开展了大量现场原位试 270地理科学27卷验、室内土工试验及黄土地区勘察经验等选取(表1),土单元采用平面四节点等参单元;由于土工格栅的模量远远大于黄土的模量,因此把土工格栅视为线弹性的杆单元;土与土工格栅的接触面采用双曲线古德曼四边形无厚单元,其参数是根据按一定比例缩小的黄土加筋模型承载模拟试验来确定的(表2);采用弹塑性模型和摩尔库仑屈服准则进行计算。表1地基和路堤填土的计算参数Table1Calculatedparametersofhigh图3最大主应力等值云图reinforcedembankmentandfoundationFig.3Continuouscontoursofσ1土层E(kPa)μγ(kN/m3)C(kPa)φ(°)填土1104000.2517.84325.72)最小主应力:最小主应力的等值云图见图地基310000.2719.23524.84,依然是在基础底面应力值最大。在填土和格栅的接触部位,尤其在整个加筋区的上下面处,出现表2接触面计算参数Table1Calculatedparametersofboundary应力集中现象,说明由于填土与加筋之间本身力学betweensoilandgeogrid性能的差异,导致此现象的发生。KsKnγCφVm类型Rf(kN/m)(kN/m)(kN/m3)(kPa)(°)(m)筋-土5110002112150.0213)网格划分和边界条件。网格划分如图2所示,共有659个单元、723个节点。模型中的左右边界条件为水平约束,底部约束条件为双向约束。图4最小主应力等值云图Fig.4Continuouscontoursofσ33)剪应力:剪应力的等值云图见图5。加筋后剪应力出现在路堤两侧一定范围内,说明在设计图2网格划分和边界条件土工格栅的位置时,不一定要在整个土体内布置土Fig.2Finiteelementmeshandboundarycondition工格栅,只要在发生塑性区的范围布置格栅即可满4)分析方法。为了探索加筋高路堤的工作足工程要求。状态和加筋机理,采用岩土二维有限元应用程序分别对加筋与未加筋黄土高路堤模拟分析。3计算结果分析根据计算结果,分别绘出加筋与未加筋高路堤的应力、变形和屈服图,并进行如下的分析。3.1应力分析1)最大主应力:最大主应力等值云图见图3。图5剪应力等值云图总体来看,它随着填土深度的增加而增加,在地基Fig.5Continuouscontoursofτ最底部为最大值;但在加筋部位,和为未加筋的路3.2塑性区堤相比较,值有所下降。说明土工格栅和填土的相塑性区的等值云图见图6。与未加筋的路堤互作用,可以使路堤中的降低。 2期王家鼎等:高填方加筋黄土路堤稳定性的有限元分析271相比,加筋后塑性区的深度和范围明显变小,未加使路堤中的降低充分利用格栅与填土自身的强度,筋路堤的塑性区在地基底部,加筋后塑性区发生在有效提高黄土路堤的强度和稳定性,减少路堤和地加筋区和未加筋区交界处,且范围减小很多。基的变形。将这种加筋土技术应用于黄土山区公路高填方路堤,具有可建筑高度大(现已可达62m)、造价低、节约用地、施工速度快、提高道路使用品质的优点,是一种有效的加固方法。2)加筋后剪应力出现在路堤两侧一定范围内,说明在设计要注意格栅的抗拉强度,同时设计土工格栅的位置时,不一定要在整个黄土体内布置图6变形后的塑性区土工格栅,只要在发生塑性区的范围布置格栅即可Fig.6Continuouscontoursofplasticstrain满足工程要求。为公路建设降低造价,缩短工期,3.3变形分析各项技术指标更加合理以及公路的设计方法与施高填方加筋黄土路堤的变形见图7。加筋路工方法等提供理论依据,使公路的修建创造更大的堤的最大沉降变形出现在路堤中心处,与未加筋路经济与社会效益,具有良好的推广应用前景。堤相比,加筋路堤中心沉降量(0.82m)明显小于3)土工格栅的使用可以减小路堤沉降和路未加筋路堤沉降量(1.9m)。由此可见,土工格栅堤横断面差异,提高黄土路堤的整体性和连贯性。对减小路堤沉降十分有效。同时,路堤横断面差异4)采用有限元方法对加筋黄土高路堤进行(即路肩和路堤中心沉降比)也明显减小,和未加了数值模拟,其结果较理想,并且该结果与现场沉筋情况相比,弯盆明显变浅。这是由于随着变形的降观测结果一致。发展,填土受剪胀作用,使得格栅与填土的相互作参考文献:用加强,提高了加筋土的整体性和刚度,有效地减小了地基中的附加应力,使路堤下的应力扩散更趋[1]王家鼎,张倬元.黄土自重湿陷变形的脉动液化机理[J].地均匀,说明土工格栅对调整地基的不均匀沉降有显理科学,1999,19(3):271～276.[2]王家鼎,惠泱河.黄土地区灌溉诱发滑坡群的系统研究[J].著效果。加筋后路堤的侧向位移比未加筋情况下地理科学,2002,22(3):305～310.的小,加筋后路堤的侧向位移发生在地表附近,沿[3]朱湘,黄晓明.有限元方法分析影响加筋路堤效果的几个深度明显减小。说明土工格栅的铺设可以减小最因素[J].土木工程学报,2002,35(6):86～92.大侧向位移值,且改变侧向位移随深度的分布,使[4]JewellRA.SoilReinforcementwithGeo-textiles[M].Lon2侧向位移集中在地表附近。don:CIRIA,1996.[5]RoweRK,SkinnerGD.Numericalanalysisofgeosyntheticre2inforcedretainingwallconstructedonalayeredsoilfoundation[J].GeotextilesandGeomembranes,2001,19:387-412.[6]Fannin,OSigrudsson.Fieldobservationsonstabilizationofun2pavedroadswithgeo-synthetics[J].JGeotech.Eng.ASCE,1996,(7):544-553.图7变形后的网格图[7]朱湘,黄晓明.加筋路提的室内模拟试验和现场沉降观测Fig.7Deformedmesh[J].岩土工程学报,2002,24(3):386～388.[8]JonesCJFP,LawsonCR,AyresDJ.Geotextilereinforced现场埋设沉降观测标7个,并及时进行了沉piledembankments[A].DenHoedt(ed).Geotextile,Geo-降观测计算,观测数据计算可靠,施工期沉降观测membranesandRelatedProducts[m].Balkema,Rotterdam,分别为:地基沉降52～58cm,路基下部沉降64～1990.IEBN90:155-160.70cm,路基中部沉降29～31cm,和计算结果比较[9]黄晓明,朱湘.公路土工合成材料应用原理[M].北京:人民交通出版社,2001.一致。[10]钱家欢,殷宗泽(主编).土工原理与计算(第二版)[M].北京:中国水利水电出版社,1996.4结论1)土工格栅在高填黄土路堤中的使用,可以 272地理科学27卷FiniteElementAnalysisonStabilityofHighReinforcedEmbankmentinLoessValley112WANGJia2Ding,XIEWan2Li,LUOFeng2Tao(1.GeologicalDepartmentofNorthwestUniversity,StateKeyLaboratoryoftheContinentalDynamics,Xi’an,Shaanxi710069;2.ChinaTransportationTong2LiHighwayEngineeringCo.,Ltd.,Xi’an,Shaanxi710075)Abstract:NumericalsimulationanalysesofthehighestreinforcedembankmentinBeizhangloessvalley,HuozhouTown,ShaanxiProvince,China,werecarriedoutwithfiniteelementmethod.Itisanalyzedforthemaximumandimumstrainsandstresses,shearstress,plasticzoneofthehighestembankment.Thecharacteris2ticsofstrainandstressofhighreinforcedembankmentandfoundationaregotten,basedonthecomparisonbe2tweenreinforcementandnon2reinforcement.Theresultsshowthatthegeo2gridnotonlyenhancethestabilityofhighembankmentandfoundationbutalsoreduceitsdeformation.TheFEManalysisresultsagreewiththefieldtests.Keywords:reinforcedloess;highembankment;geo2grid;analysisoffiniteelement中国城市群第三版已发行自1992年和2001年《中国城市群》第一版、第二版相继问世以来,中国的经济社会得到了迅速发展,城市化进程明显加快,城市化问题已列入国家的重大发展战略。国家的十一五发展规划纲要又明确指出:城市群的发展将作为中国推进城镇化的主体形态。姚士谋、陈振光等的城市群著作第三版于2006年12月修改后再版就是顺应了这个重要的历史时刻。自90年代中后期,中国已有20多个省区在建设部与各省建设厅的主持下,做了不少城市群规划研究,中国城市群将近80万字的著作,根据中国各地大量的规划实践,结合国内外的研究动态与经验,比较系统地分析了中国6大城市群、7个城镇密集地区的形成发育条件、历史与自然基础,提出了有关城市群的概念、结构特征、空间布局、类型和指标,分析了不同类型、不同规模、不同功能的城市群区域特征与问题。同时运用科学发展观,深入研究中国城市化的问题与趋势,提出了城市化健康发展的新思路,为中国城镇化与城镇发展提出了重要的理论依据,也是中国系统地、综合地研究城市群问题的重要著作,起到了重要的指导作用。本书曾于1998年获得中科院自然科学二等奖以及2004年华东地区优秀图书奖。本书第三版出版,周干峙院士、吴传钧院士及崔功豪教授为本书作了序言,并热情地推荐本书。本书共分为10章,内容丰富、观点鲜明、资料翔实、图文并茂,并附有500多条专业英文生词、短语,欢迎广大读者购买本书,定价:58元。订购地址:合肥市金寨路96号,中国科学技术大学出版社,邮编:230026,联系人:高哲峰。