地下洞室群围岩稳定三维非线性有元分析.docVIP

PAGE PAGE 609 地下洞室群围岩稳定三维非线性有限元分析 余卫平 汪小刚 中国水利水电科学研究院，北京，100044 摘要：大型地下洞室群的稳定性是地下水电站面临的一个关键问题，三维非线性有限元方法常常用于评价地下洞室群的稳定性。本文介绍了基于增量变塑性刚度迭代法的岩体开挖的弹塑性有限元分析方法，提出了可以考虑抗剪作用的隐式锚杆单元和可以考虑预应力效果的隐式锚索单元，并引入可以量化的各种破坏区体积作为地下洞室群围岩稳定性评估指标。从而为大型地下洞室群围岩稳定三维非线性有限元分析提供了完整的思路和方法。工程应用实例表明，这一思路和方法是切实可行并且有效的。 关键词： 弹塑性有限元 围岩稳定 隐式锚杆（锚索） 单元 1 引言 目前水利水电领域地下洞室群的规模越来越大，其布局和所处的地质环境也越来越复杂。例如我国拟建和正在建设的溪洛渡、小湾、龙滩、水布垭、糯扎渡等水电站地下厂房洞室群，其规模都是空前的，洞室群的围岩稳定对工程安全至关重要，有必要对洞室群的围岩稳定性进行全面的分析和评价。有限元方法已经成为地下洞室群围岩稳定性分析的有力手段。许多文献在基于有限元的地下洞室群围岩稳定性的分析方法、锚固支护的模拟以及稳定性评价指标等方面进行了许多探索，目前这一领域的研究仍然十分活跃[1～3]。但运用有限元方法分析评价复杂地质条件下大型地下洞室群的围岩稳定性，并经受工程实践的检验，仍需要广大岩土工程科技人员的艰辛探索。首先，大型地下洞室群面临的地质环境复杂，涉及地应力、温度和地下水等复杂的地质环境。其次，大型地下洞室群面临复杂的工程地质条件，地层构造千差万别，不同的地层岩性、地质构造具有不同的本构关系和力学特性，其数值模拟方法需要人们去探索。最后，大型地下洞室群涉及复杂的支护结构，如锚杆、锚索、混凝土喷层以及钢拱架等，这些支护结构的有限元模拟方法也需要人们去探索。 本文尝试从岩体施工开挖的弹塑性有限元分析方法、锚固支护的数值模拟方法以及岩体开挖施工扰动条件下围岩破坏指标的确定等几个方面来探讨基于三维非线性有限元的大型地下洞室群围岩稳定分析方法，并通过工程实例的计算分析情况来验证本文提出思路和方法的可行性和有效性。 2 岩体开挖的弹塑性分析方法 2.1 弹塑性本构关系和迭代计算方法 本文采用增量变塑性刚度法进行迭代计算[1]。该方法是将洞室群的开挖荷载分解为弹性荷载和塑性荷载两部分。对弹性荷载采用一次施加于结构上进行计算，对塑性荷载则分级进行加载计算。对于每一级塑性荷载，可按式(1)进行迭代计算： (1) 式中：和为结构的弹性刚度矩阵和塑性刚度矩阵；和分别为第级塑性荷载和位移增量。迭代计算时，屈服函数采用Zienkiewicz-Pande的双曲线屈服准则[4，5]： (2) 式(2)中各变量的说明参阅文献[4，5]。 当按式(2)判断单元进入塑性状态后，应先将单元应力状态沿屈服面的法向拉回到屈服面，再根据屈服面上的应力状态修正塑性刚度矩阵代入式(1)进行迭代计算。 对破坏状态判别和迭代完毕后，由于有限元数值计算的误差，单元的应力状态不一定位于屈服面上。为保证迭代计算的收敛性，使屈服单元的应力状态沿屈服面滑动，对每一级增量荷载迭代计算完毕后，可以将单元应力状态沿垂直于屈服面的方向将应力状态拉回到屈服面上来。所修正的应力值按下式计算[6]： (3) 式中：为修正前应力状态对应的屈服函数。 将应力状态修正完后再进入下一步增量荷载的迭代计算。 2.2 分期开挖施工优化评估方法 对不同方案和施工方法进行优化评估时，应综合分析围岩的变形、应力分布规律与围岩开裂破坏区、塑性破坏区以及回弹破坏区的大小，并考虑锚杆（索）应力情况，最终对围岩稳定情况作出评价。以下介绍这几种破坏区的判别准则[7]。 开裂区（TD），是指围岩受张力而产生的张裂破坏。本文采用岩体的临界应变作为判断岩体单元是否产生张裂破坏的准则。即认为围岩的最大应变值大于围岩的临界应变时，围岩产生开裂。 塑性区（PA），是指洞室开挖后围岩应力状态达到临界应力的屈服岩体区域。通常可根据有限单元计算的应力状态按所选择的屈服函数F进行判断。当单元的屈服函数时，则认为该单元岩体已经产生屈服破坏。 卸荷区（RE），是指在分期开挖过程中前期开挖已经进入塑性和开裂破坏的单元，在后期开挖过程中岩体应力产生卸载，返回到屈服面以内的岩体。由于塑性破坏是不可逆的，卸荷区也是一种破坏区。 在计算分析时可统计各种破坏区的体积，以进行各方案的比较分析。 3 锚固支护计算分析方法 图1 锚杆单元示意图 为了便于岩体单元的网格剖分，采用隐式的考虑抗剪（销栓作用）的杆单元来模拟系统锚杆，用考虑预应效果的隐式杆单元来模拟锚索的加固作用[8]。将锚杆单元隐埋在岩体