陕南蓝商高速公路程家院高边坡稳定性的数值模拟分析.docVIP

陕南蓝商高速公路程家院高边坡稳定性的数值模拟分析 　　摘 要：针对陕南蓝商高速公路程家院岩质高边坡的地质和工程地质条件，利用FLAC3D有限差分软件建立三维工程地质模型，对该高边坡的六级开挖过程进行了数值模拟，通过分析边坡开挖过程以及边坡形成后的应力、应变规律以及塑性区的分布特征，研究了边坡开挖过程对边坡稳定性的影响。同时运用强度折减法计算了边坡的安全系数，并结合数值模拟结果对该边坡的稳定性进行了综合评价，以此为依据，提出了该边坡治理的技术方案。 　　关键词：高速公路 高边坡 数值模拟 边坡稳定性 　　一、引言 　　随着我国国民经济的发展，高等级公路修建在我国发展很快，但高速公路建设工程中所遇到的岩质边坡稳定性问题也相应地增多。能否对边坡的变形、位移、破坏和灾变进行准确的分析计算和预测预报，直接关系到人们生命财产、物资财产和环境资源的安全。边坡稳定性分析研究自然而然地成为这项工作中的核心内容[1-3]。 　　目前边坡稳定评价方法除传统的极限平衡分析法外，数值模拟分析也是主要方法之一。 　　二、工程地质条件及建模 　　1.程家院高边坡概述 　　陕南蓝商高速公路是一条山区高速公路，地形地质条件较为复杂，路堑高边坡发育众多。程家院高边坡位于高速公路K59+300～K59+530路线右侧，走向93°，平均坡面倾角45°左右，坡体上缓下陡，除在坡体上部较缓的地区上覆4～10cm厚的土层外，主体由碎裂状强风化岩石组成，坡面产状大体一致的节理有三组（55°∠30°，5°∠25°，35°∠12°），利用工程地质类比法分析，坡体中上部的土体及全风化岩体由于节理、裂隙作用，容易相切和追踪形成倾向临空不良结构面，边坡设计开挖等级为六级，一级坡率为1?0. 5，二级为1?0.75，三四级坡率1?1，五六级坡率为1?1.25，边坡开挖后，岩石风化加重，边坡极易产生坍塌、错落和滑移等现象。 　　2.数值模型的建立 　　根据边坡地形地质条件，建立的设计边坡计算模型。计算模型沿x向边坡倾向宽度为150m，沿y向边坡走向长度为250m，z向边坡垂直高度150m。边坡几何模型开挖面的生成均在ANSYS中完成，网格划分后保存单元和节点几何信息，然后通过接口程序转化为FLAC3D的前处理数据格式。在FLAC3D中导入这些数据之后生成的网格模型（图1），共7220个节点，14231个单元。计算模型除坡面为设为自由边界外，模型底部（z=0）设为固定约束边界，模型四周设为单向约束。在初始条件中，不考虑构造应力，仅考虑自重应力产生的初始应力场。岩土体强度参数的选取以工程地质岩石力学研究成果为依据。模型为近似理想的弹塑性模型，采用摩尔库仑准则。岩土体物理力学参数的具体取值见表1。 　　表1 岩土体力学参数表 　　Table 1 Mechanical parameter of rock and soil 　　层 　　序 岩 性 容重 　　（kg/m3） 弹性 　　模量 　　（MPa） 泊松比 抗拉 　　强度 　　（MPa） 内聚力 　　（MPa） 摩擦角 　　（o） 　　1 上部强风化 　　岩石 2300 5000 0.25 0.6 0.4 25 　　2 下部中风化 　　岩石 2300 6000 0.25 0.7 0.6 25 　　三、数值模拟结果的分析 　　1.开挖过程的对比分析 　　从图2中可以看出，随着边坡的开挖边坡，在边坡的中下部产生较大的沿坡向的位移，由第一步开挖的3mm，到开挖最后一级中达到了29mm。这主要是由于随着开挖的进行，由上往下边坡的坡率逐渐加大，最后一级边坡的坡率达到了63°，对整体边坡的坡率改动较大，且边坡主要是破裂的全风化岩体很容易发生向下的滑塌，从位移图2可看出开挖过程对边坡的稳定性影响较大。 　　（a）～（d）分别代表开挖第一、三、五、六步水平方向位移 　　图2 边坡开挖过程水平方向位移对比图 　　2.开挖后边坡的稳定性分析 　　2.1位移分布规律 　　图3、图4分别为边坡的整体位移云图、铅垂方向位移云图。从整体位移云图来看，位移最大的部分集中在边坡后部；铅垂方向位移云图与整体位移云图相似，且数值相近，且最大铅垂位移也出现在边坡后部。这表明边坡在开挖后，边坡后部的位移以“沉降”模式为主。图2（b）中最大水平位移则出现在边坡的中下部，并以此为中心，水平位移呈“同心圆”状逐渐减小向四周扩散，这表明，相对其它地段，该地段岩体水平向变形突出，向水平方向方向隆起。 　　2.2应力分布规律 　　（1）图5中（a）、（b）分别为边坡的第一主应力和第三主应力云图（拉应力为正，压应力为负），从边坡主应力云图来看，边坡由于开挖卸荷在坡体两侧出现了范围很小拉应力区，其余地方均以压应力为主，若边坡发生