岩质边坡持续降雨工况下的稳定性分析.docVIP

岩质边坡持续降雨工况下的稳定性分析 　　摘 要：伴随着国民经济的日益发展，交通运输、矿山开采、水利及国防设施等项目的建设，工程中遇到岩质边坡的问题越来越多，其不良影响引起了人们的关注，而边坡稳定性评价则是边坡工程的核心问题，是岩土工程领域的一个重要课题，该文针对一岩质边坡采用ABAQUS分析方法对其的稳定情况进行分析，首先根据该边坡工程地质剖面图及其边坡物质组成和工程地质性质建立边坡变形演化计算模型，再在持续降雨工况下进行数值模拟分析，最后得出可靠结论，为该边坡的除险加固提供依据。 　　关键词：岩质边坡 ABAQUS 有限元数值模拟 　　中图分类号：TV223 文献标识码：A 文章编号：1672-3791（2015）02（a）0049-02 　　边坡滑坡是地质灾害中的一种，是山体斜坡自然或人为的变形现象，能给人类造成巨大的生命财产损失，而边坡稳定性分析是判断边坡是否失稳、是否需要加固及采取防护措施的重要依据，是边坡工程中最重要和基本的问题。 　　要降低滑坡所造成的危害，则需要对滑坡进行深入的研究，首要任务就是正确的评价滑坡的稳定性。 　　随着计算机的迅速发展，数值模拟方法也得到迅猛的发展，目前边坡稳定性分析中也大量采用该方法，该文主要采用ABAQUS有限元数值模拟来对某水电站一级厂房后边坡进行稳定性分析与评价。 　　1 ABAQUS分析方法 　　1.1 计算机模型建立 　　ABAQUS分析方法可根据有限元基本方程、Mohr-Coulomb模型及库仑摩擦模型针对边坡物质组成和工程地质性质进行建模[1]。 　　（1）有限元基本方程。 　　有限元分析是在地质体力学形态本构模型和工程地质模型的基础上求解地质体应力及稳定状态过程的数值模拟方法，以达到对地质体变形破坏进行描述的目的。有限元基本方程包括平衡方程、物理方程和几何方程。 　　①平衡方程。 　　②物理方程。 　　根据塑性增量理论，各向同性材料的应力应变关系可以写成式4。 　　③几何方程。 　　（2）Mohr-Coulomb模型。 　　Mohr-Coulomb屈服准则[2]假定：作用在某一点的剪应力等于该点的抗剪强度时，该点发生破坏，剪切强度与作用在该面的正应力呈线性关系。 　　（3）库仑摩擦模型。 　　ABAQUS中的库仑摩擦（Coulomb friction）模型[3]可以用来描述滑坡体物质分异面处之间的相互作用。库仑摩擦模型用摩擦系数来量化接触面之间的摩擦行为。如果两表面间的接触压力为，则接触面之间的摩擦力的上限值为，当接触面上的剪力达到该上限值时，接触面就会发生相对移动。ABAQUS通过自动选择一个“弹性滑动”的罚刚度因子，表示在粘结的接触面之间允许发生少量的相对位移，这个允许的“弹性滑动”是非常小比例的特征单元尺寸。 　　1.2 选取边界条件 　　边界条件是数值计算结果精确程度的重要影响因素，它的选取将直接影响到数值计算结果的可靠性。边界条件可分为固定边界条件，黏性边界条件，远置边界条件等[4]，具体的选择根据工程情况反复调试以达到选择的合理性。 　　1.3 单元选取 　　模型单元选用ABAQUS中的修正6节点三角形二次平面应变单元（CPE6M），此单元精度较高，且能模拟任意的几何形状。在模型中，滑动带等比较关心其应力应变状态的部位布置较密的单元，而在滑床布置较疏的单元。 　　1.4 物理力学参数 　　对于Mohr-Coulomb模型，需要输入的参数为弹性模量，泊松比，内摩擦角，凝聚力，密度，剪胀角及初始塑性应变等。 　　2 工程实例 　　某水电站一级厂房后边坡发育在两冲沟之间，平面上呈北东向狭长型展布，山体三面临空，冲沟发育，各条冲沟向上延伸上至尖峰并汇集于尖峰之下，由尖峰山脊向上下游的延伸，在地貌上形成圈椅状背形态，各条冲沟座落于圈椅内，形成“爪状”地貌形态。山坡和冲沟内有大量有崩积堆积层，厚度20m，部分冲沟下雨时很少见到流水现象。山坡上直径约15cm以上的部分松树、彬树在树根段可见马刀形弯曲，具有“马刀树”特征。 　　2013年的一次持续降雨，使该边坡的地下水位急剧上升，滑面的抗剪强度降低、滑体的容重增大、静水压力升高，又由于原来厂房基坑和坡脚的开挖破坏了坡体原有的稳定，边坡出现了明显的位移，现就该边坡持续降雨工况下进行数值模拟分析。 　　2.1 建立计算模型 　　根据该边坡工程地质剖面图及其边坡物质组成和工程地质性质建立边坡变形演化计算模型，如图1所示。将边坡划分为熔结凝灰岩滑床，夹碎块石粉质粘土和熔结凝灰岩破碎带组成的滑体和泥砾层组成的滑动带。 　　2.2 滑动面抗剪参数反演 　　运用不平衡推力法反算模型的物理力学参数，计算结果如表1所示。 　　2.3 持续降雨状态下数值分析结果