基于FLAC3D复合格栅加筋岸坡数值分析.docVIP

基于FLAC3D复合格栅加筋岸坡数值分析 　　摘要：为减少环境影响和尽快恢复生态环境，某滨江路的护岸工程采用低桩承台＋两级格栅加筋陡坡的复合式护岸新型结构。本文采用FLAC3D有限元差分方法，分析了复合格栅加筋岸坡的整体稳定性、低桩承台的受力和变形、格栅加筋岸坡体的受力和变位。并分别对岸坡加筋和不加筋作了对比分析。得出复合岸坡加筋后，桩台受力状况改善、变形显著减少，岸坡的变位也显著减小。计算分析结果为新型结构在工程中的应用提供了依据，也可为类似工程参考。 　　关键词：复合格栅加筋岸坡；有限元差分方法；数值分析；整体稳定性；变形 　　 引言： 　　随着经济建设的发展和需要，城市防洪护岸综合整治工程建设项目也越来越多。三峡工程蓄水后，库区城市建设及经济发展与紧缺的土地资源的矛盾越来越突出，同时，城市防洪护岸工程建设的自然环境条件也更加复杂，人们对工程建设的生态恢复、环境与景观要求也越来越高。要解决这些矛盾，需要用新的工程建设和工程设计理念来研究新的结构和提出新的工程实用技术。 　　城镇防洪的护岸工程按新的要求和设计理念，应尽可能节约土地，减少对河道的占用，减少对环境的影响，尽可能快地恢复生态环境，同时要尽可能形成较多的建设用地，充分考虑工程建设的综合效益。近20年来，在重庆已建的100多公里长的滨江路工程中，针对不同的建设要求和建设条件，滨江路工程采用了路堤式和桥梁型式，路堤护岸工程又分别采用了衡重式挡墙、空箱式挡墙、连拱挡墙、无底肋锚扶壁挡墙、加筋土挡墙、加筋土边坡、抗滑桩柱挡墙等多种支挡结构型式。加筋土挡墙和加筋土边坡对于护岸工程，在技术先进性、工程经济性、施工方便性等方面具有明显的优势。三峡工程蓄水后，原枯季低水位情况不复存在，原在河岸大开挖对河道带来很大影响的施工方法也不能再用。鉴于此，某段滨江路工程根据变化的情况，按新的建设要求和建设条件，提出了采用低桩承台＋两级格栅加筋陡坡的复合式护岸的新型结构型式，能较好地解决工程建设的工期、河道环境保护、生态恢复、城市景观、节约工程投资、增大工程综合效益等各项目标要求。 　　对提出的上述新型复合式加筋土结构，结构的受力和变形计算，现有的规范并未给出相应的现成计算方法。本文在以往研究的基础上，应用FLAC3D，对复合加筋土护岸结构的整体稳定性、低桩承台的受力和变形、加筋土护岸工程的受力和变形进行数值分析，为工程应用提供计算分析依据。 　　一、工程概况 　　1.复合结构型式。某滨江路堤防工程推荐采用的复合式加筋陡坡结构，其代表断面如图1所示。坡脚采用双桩低桩承台结构。两级加筋陡坡，分级高度各为10m，坡比1：1，两级陡坡间留有4m的马道。坡面采用有防冲刷且易于快速恢复河岸植被生态的特殊坡面结构。加筋体部分均采用砂卵石回填，并分层碾压密实，按设计要求及其他工程现场检测，相对密实度均在93%以上（重型标准）。 　　图1 复合构断面图 　　2.岩土体及材料计算参数。结构断面岩土体从下到上依次为中风化岩层、强风化岩层、粘土层和开山土加筋体回填层。桩为Φ1.2m嵌岩灌注桩，。其中材料强度参数如下表1： 　　目前，针对界面的接触问题已经有了很多的接触单元，如Goodman单元、Desai薄层单元。本文采用FLAC3D程序提供的土工格栅接触单元（geogrid）来模拟土体与格栅的界面特性。Geogrid单元的特点是在法线方向上与格栅表面刚性连接，在切线方向相当于格栅表面的弹性滑块。其具体数值见下表2： 　　二、FLAC3D计算原理及求解流程 　　1.基本计算原理。FLAC3D是基于拉格朗日差分方法来求解运动方程和动力方程的一种理论，相对于其他软件，在算法上，FLAC3D即使模拟静态系统，也采用动态运动方程进行求解，这使得FLAC3D在模拟物理上的不稳定过程不存在数值上的障碍。首先，程序将计算区域内模型的介质离散化成若干二维单元，这些单元之间由节点相互联结，当对某一个节点施加荷载作用后，其运动方程可以写成时间步长为Vt的有限差分形式。在某一微小时间段内，作用于该节点的荷载只对周围若干节点有影响。根据单元节点的速度和时段变化，求出单元之间的相对位移，进而求出单元应变，然后再利用单元材料的本构关系求出单元应力。随着时间段的增长，这一过程将扩展到整个区域，从而求出单元之间的不平衡力，然后将不平衡力重新作用到节点上，再进行下一步迭代过程，直到平衡力足够小或节点位移趋于平衡为止。 　　2.求解流程。FLAC3D基本的求解流程如下图2所示： 　　图2 FAC3D一般求解流程 　　三、计算模型 　　由于本模型相对复杂所以采用Ansys软件建立模型，然后通过ANSYS-TO-FLAC3D软件导入到FLAC3D中，再进行计算。考虑到边界条件和网格疏密对模型计算结果的影响，本模型取坡脚