高边坡稳定性分析.docxVIP

高边坡稳定性分析 1 项目总结 1.1 高边坡岩体结构 工程区（图1）位于大巴山以北、鄂西山脚下的附近。4号公路d区的高坡度是在4号公路和2号公路之间挖掘的斜坡。坡度为91，呈顺坡，坡度高度为19.726.0m，坡度角为40-55。这是岩石的坡度。斜坡的岩性主要为细至中层状粘土岩，夹细泥灰岩。 场区岩体裂隙较发育,据现场调察,该高边坡岩体发育三组裂隙,其中控制裂隙为:90°～120°∠ 35°～50°,面较平直～平直,宽0.1～2.5cm,最宽达5cm,内有方解石及灰绿色泥钙质物充填.间距0.1～0.3m,延伸长3～12m。目前边坡坡度大于裂隙倾角,坡高超过15m,边坡不稳定,最不利结构面为1-1′剖面ZK2孔9.7m处发育的35°裂隙。 1.2 层序地层填土为新堆填材料本构模型 区内地层岩性分三类:第四系人工填土层、泥岩及泥灰岩,且以薄～中层状泥岩为主,夹薄层泥灰岩。第四系人工填土为新近堆填,分布不均,力学性差,工程意义小,且场区泥灰岩呈薄层状夹于泥岩中,厚度薄,1.8～4.8m,工程意义不大,故未对它们进行试验。本次勘察取中风化泥岩若干组,按《岩土工程勘察规范》(以下简称《报告》)计算成果见表1。 2 模型分析 2.1 超静定问题的计算方法 对于边坡稳定分析,目前国内采用的方法主要是条分法,它们的共同点是:假定土体为理想塑性材料,并把土条作为一个刚体,按极限平衡的原则进行受力分析,完全不考虑土体本身的应力～应变关系。各种方法最大的不同之处仅仅在于相邻土条间的内力如何假定,即如何将超静定问题化为静定问题。而有限元法与传统的条分法相比,它不用假定条间力就能使土体内各部份满足平衡条件,边坡失稳自然地发生在土体抗剪强度低于实际剪应力的区域。 本工程是个静力平面应变问题,采用非线性分析。分析的重点是如何在模型中反映节理裂隙的影响,这方面的文献较多,但总的来看,没有公认一致的计算方法。众所周知,坡体失稳主要原因是因为土体抗剪强度低于滑动力而发生的剪切破坏。影响抗剪强度的因素很多,对于岩体,主要是岩体本身物理力学指标、排水条件及节理面的发育等。根据现场调查及钻探揭示的裂隙情况综合分析判定,该高边坡最不利结构面已经确定(图2),因此我们可以认为岩体中存在一个位置已知的软弱结构面(软弱结构面的土性参数与其它区域不同),它是由于节理裂隙的存在而产生。 2.2 分析单元和方法 根据工程地质勘察报告提供的数据,有利用Ansys有限元计算程序,建立了计算模型。有限元模型如图3所示。坡体倾角为50°,坡高为26m,均按工程地质剖面图1-1′取,横向宽度向两边分别延伸高度的1.5倍,以减少模型左右边界的约束(只约束水平方向)对坡体位移的影响,底部约束为固定支座。软弱结构面倾角为35°,共取五条,间距1m,其宽取控制裂隙的中间值10 mm。只考虑重力的作用,由于场区内无地下水,故本次分析不考虑孔隙水压力。分析目标是确定以软弱结构面为滑动面在载荷作用下的安全系数。 单元选择二维结构实体单元PLANE42,PLANE42单元由四个节点定义,每个结点有两个自由度:节点坐标系X,Y方向的平动。该单元有塑性功能,可用于模拟岩体。岩体本构模型选取了Drucker-Prager准则,即假定岩体为理想塑性材料,不产生硬化,同时,假定岩体的膨胀角为零,即在变形过程中岩体不发生体积改变,也即采用非关联流动法则。为了减少迭代次数,此处采用修正Newton-Raphson法求解。 其收敛准则是以各节点处力的不平衡量的范数(各节点力矢量各项平方和的平方根SRSS)是否等于各节点(等效)加载力的范数的10-3,10-3也称为容差。 2.3 滑移表面土体抗剪强度验算 本次试算共采用了两组参数,一组由设计单位地质勘察报告提供,由于这组参数过小,试算时发现安全系数很小,坡体很不稳定,而现场坡体并未出现失稳征兆,可能是参数取的过于保守,因此拟选用另一组为调整后参数,如表2所示。 以过坡角的软弱结构面为假定滑移面。由于弹性模量的大小对安全系数的求解影响不大,所以此处非节理区的弹性模量取设计指标的约一半,软弱结构面弹模取非节理区的一半。 求解可能滑移面的安全系数时用到Mohr-Coulomb准则(图4),即滑移面上土体抗剪强度与实际剪力的比值。计算公式如下: 滑动面上某一点的应力状态沿坡向的分量: σi=σix+σiy2+σix?σiy2cos2α?τixysin2ατi=σix?σiy2sin2α+τixycos2α(1)Fs=τfτs=∑i=1n(c+σitgφ)Δli∑i=1nτiΔli(2)σi=σix+σiy2+σix-σiy2cos2α-τixysin2ατi=σix-σiy2sin2α+τixycos2α(1)Fs=τfτs=∑i=1n(c+σitgφ)Δli∑i=1nτiΔl