拱坝稳定计算理论方法论述.docVIP

拱坝稳定计算理论方法论述 　　摘要：拱坝在施工期间与运行期间均有不同程度的开裂情况，有些裂缝的出现严重影响了拱坝的安全，甚至是溃坝的前兆，而有些裂缝的出现改变了拱坝的应力结构，释放了超限的应力从而改善了拱坝的应力状况。所以，检查拱坝是否开裂并对开裂坝体进行裂缝成因及修补措施分析的关键是计算坝体应力应变、复核坝体抗滑稳定性。拱坝坝肩坝基抗滑稳定分析有许多方法，常用三种主要的分析方法有刚体极限平衡法、模型试验和有限元法。 关键词：拱坝 稳定 理论 中图分类号： TV642.4 文献标识码： A 文章编号： 随着计算机广泛应用于自然科学和有限元方法的不断成熟，有限元法应用更加广泛，特别是在大、中型水利工程中已经成为不可替代的有效方法。 1刚体极限平衡法 与其他分析坝肩抗滑稳定的方法相比，刚体极限平衡法理论成熟、概念清晰、计算简单，为过去和现阶段的工程所普遍采用，也是目前规范规定采用的方法。 在刚体极限平衡法中，有二维和三维两种计算方法： 二维抗滑稳定分析即局部稳定分析，一般是取任一高程单位高度的拱圈，并将每层拱圈对应的坝肩抗力体岩石沿垂直方向切取单位高度，或沿水平方向切取单位宽度进行计算研究。因此，如果以该法计算的各高程拱圈均能独立维持稳定，则拱坝坝肩的整体稳定当更无问题。反之，少数拱圈的稳定性不足，并不意味着拱坝一定不安全，而应该进一步研究拱坝的整体稳定性是否满足要求。 三维抗滑稳定分析即整体稳定分析，该法是从坝肩被断层、节理、裂隙分割的整个岩石块体群中，人为地取出若干被判定为有可能滑动的块体进行力学分析计算，然后比较其安全系数值，取最小的一个作为坝体的安全系数。该方法力学概念清晰，与二维方法相比，能更全面的体现出拱坝受力条件与地质薄弱面分布的复杂性，因此三维刚体极限平衡法使用更多。 在傅作新教授提出的拱坝的上滑稳定分析的改进方法中，也使用了刚体极限平衡法[21]，总体而言，刚体极限平衡法是一种传统的稳定分析方法，其理论成熟、概念清晰、计算简单，为过去和现阶段的工程所普遍采用，但是该方法计算比较粗略，引入较多假定，采用的岩体物理力学模型属于“刚塑性”，因此，在拱坝坝肩稳定分析方面还有其不足之处[1]。 2有限元法 从结构力学的角度出发，有限元法实质上是将一个连续的实体人为地离散为若干单元，相邻的单元在结点处相互连接。先对每个单元的应力应变特性进行研究，计算单元的刚度矩阵，进而组成整体刚度矩阵，再以静力等效的方式把所有的外荷载转移到各个结点上，组成结点的荷载列阵，最后通过结点上的平衡条件及整体的边界条件建立方程，计算出结点位移，进而求出单元的应变和应力[2-3]。 随着线性及非线性有限元理论的不断完善，有限元法在拱坝坝肩稳定分析中的应用也越来越广泛。我国有关规范中规定：拱座抗滑稳定的数值计算方法以刚体极限平衡法为主，1、2级拱坝或地质情况复杂的拱坝还应辅以有限元法或其他方法进行分析[4]。 采用有限单元法分析拱坝坝肩稳定时，可以将坝体、坝肩作为整体计算，得到坝体内的应力和变形分布，据以分析坝肩的稳定。一般取一定范围的岩体作为拱坝的地基，在此范围内划分坝体及地基岩体单元，要求结点处三维力平衡、三维变位相容。其优越之处在于：拱坝各单元的材料特性可以不同，能更好地符合拱坝材料分区的实际情况；地基内各单元的材料特性可以是非均匀、各向异性的，还可以考虑地质构造的影响；坝体内各单元之间力的平衡和变位相容是对每一个结点的要求，这样能更好地符合坝体内力平衡和变位共容的条件。因此，用有限元方法进行坝肩稳定分析逐渐成为一种发展趋势，通过对坝体及地基岩体的应力、变形分析，进而分析坝肩岩体及坝体结构的稳定性。 根据毕肖普50年代重新定义的滑坡稳定安全系数概念，即每个条块的安全系数等于该条块底部的抗剪强度与土体发挥的剪应力之比，再根据Mhor—Colun准则，最后对整个滑动面进行加权处理，取条块底部的弧长为权函数，得整个滑面的安全系数。求得安全系数值后，取其中的最小值作为坝肩的抗滑稳定安全系数值，其对应的滑动面即为最危险滑动面。 显而易见，采用有限单元法能够充分考虑岩体不同区域材料的力学特性、坝肩岩体构造的复杂性以及岩体变形对坝体结构的影响，并且能够一次算出稳定安全系数。另外，在分析过程中，可以得到危险滑移面上应力、应变及稳定安全系数的变化规律，定量的指出最危险区域的位置，这是前面几种评价安全度方法所不及的[5]。因此，用有限元方法进行坝肩稳定分析逐渐成为一种发展趋势。对于拱坝的上滑稳定问题，也可以用非线性有限单元法进行分析。此法可以给出拱坝开始加载直至破坏全过程的应力和变形规律。为了检验沿坝基面上滑的可能性，可以在拱坝坝基面布置薄层单元，并在计算时校核沿坝基面和横缝面