有限元强度折减在边坡稳定性分析中的应用.docxVIP

有限元强度折减在边坡稳定性分析中的应用 1 元一开始就是再分配 1.1 边坡稳定性验算 元强度折叠法的基本思想与传统的轮廓平衡法一致，可以称为强度储存系数。基本原则是将岩质参数c和c值从sf中排除，并获得一组新的cf和f值。然后继续计算新材料参数，直到计算收敛。在满足一定计算精度的前提下, 取最后一次计算收敛时的安全系数为边坡稳定的安全系数。此折减系数对应的Fos被称为坡体的最小安全系数, 此时岩土体达到临界状态, 发生剪切破坏。 Fos=τ/τf 式中:Fos为边坡安全系数;τ为边坡岩土体的抗剪强度 (τ=c+σntanφ) ;τf为边坡滑面剪应力,τf=cf+σntanφf;cf=c/SRF;φf=arctan|1/SRFtanφ|φf=arctan|1/SRFtanφ|。 1.2 应力张量的计算 本文用折减系数法求解边坡稳定问题时, 采用的是理想弹塑性模型, 其本构模型采用Drucker-Prager准则。 F=αI1+J2??√=kF=αΙ1+J2=k 式中:I1、J2分别为应力张量的第一不变量和应力偏量的第二不变量;α、k是与岩土材料内摩擦角φ和粘聚力c有关的常数, 其中: α=sinφ/3(3+sin2φ)??????????√k=3√ccosφ/3+sin2φ????????√α=sinφ/3(3+sin2φ)k=3ccosφ/3+sin2φ 其屈服面在π平面上为不等角度的六边形的内切圆, 用此准则时, 塑性区最大, 结果稍偏保守。 1.3 同一边坡的破坏荷载 有限元计算中, 采用关联流动法则还是非关联流动法则, 取决于ψ值 (剪胀角) :当ψ=φ时, 为关联流动法则;当ψ≠φ时, 为非关联流动法则。总的来说, 采用非关联流动法则所得破坏荷载比同一类型材料采用关联流动法则所得的破坏荷载小, 如果忽略剪胀角 (ψ=0) , 将会得到较为保守的结果。值得注意的是:当ψ=0时, 正好与郑颖人等提出的广义塑性力学理论相符, 这时对应的塑性面与q轴垂直。相关文献的结果表明:对同一边坡, 不论采用关联流动法则还是非关联流动法则, 计算结果相差不大。因此, 本项目中的所有计算都取ψ=0, 即满足非关联流动法则。 2 工程计算的示例 2.1 数值模拟施工过程 某顺层岩质边坡主要有第四纪残积土、全风化、强风化、弱风化和微新岩层等地层和F2093、F2055、F2059、F2074等断层构成 (见图1) ;根据工程需要, 按表1进行开挖施工, 但开挖支护完成后, 该边坡发生座滑。本文通过有限元计算模拟施工过程, 并由此评价该边坡的稳定性。 计算时, 计算范围为X向取700 m,Y向自地表至高程0 m。计算域内包含了F2093、F2055、F2059、F2074、F2106、F2113、F2078、F2085等断层以及第四纪残积土、全风化、强风化、弱风化和微新岩层等地层。计算域的底部采用铰支约束, 两侧均采用法向约束。边坡的有限元网格剖分模型如图2所示, 网格剖分均采用6节点三角形单元。各类围岩均采用弹塑性Mohr-Coulomb模型, 岩体力学参数计算采用值见表2、表3。计算中考虑了水的影响, 并模拟了支护过程。 2.2 预应力释放型 由图3～6可知, 边坡开挖后, 由于围岩卸荷, 下部应力释放, 主应力减小, 而孔隙压力不变, 所以有效主应力减小。进行支护后, 在边坡支护段附近的应力场主应力都变大, 但离支护较远处应力几乎没有影响。 2.2.1 巷道围岩变形随边坡变形的变化 图7显示:边坡在开挖280台阶以上部分变形不显著, 在开挖260～280高程时, 变形突然加大。开挖后进行支护对边坡上部影响不大, 对边坡下部而言, 将有效减小岩体变形。从图8～10可以看出, 260平台开挖后, 节理间的剪切位移比较大, 整个边坡有沿着在开挖区出露的断层F2085、F2078向下移动的趋势。 2.2.2 塑性区的减少 图11～13显示, 边坡的屈服单元与塑性区、最大剪应变主要集中在坡顶全、强风化带内;支护完成后, 塑性区有所减少。 开挖后, 采用有限元强度折减方法计算边坡的安全系数, 得到座滑前边坡安全系数为0.837, 支护后安全系数为0.943, 较支护前提高了0.106, 但安全系数均小于1, 所以边坡整体不稳定 (见表4) 。 3 基于极限平衡法的边坡稳定分析方法 本文将有限元强度折减法成功应用于某岩质边坡工程的稳定性分析中, 得到了边坡稳定的安全系数, 不仅能有效地模拟边坡的开挖过程, 无需假定滑动面的位置和形状, 而且可得每一步骤中应力应变的演变过程, 程序会自动地寻找坡体的破坏面。工程实践表明, 这一方法是适于边坡稳定分析的一种有效的实用方法。主要结论如下: 1) 有限单元法不需要做任何假定, 计算模型满足岩土体应力应变关系,