钉结护面板作用下土石坝地震稳定性分析.docxVIP

钉结护面板作用下土石坝地震稳定性分析 0 高土石坝地震稳定性分析 中国已经建造了近100米的高坝，主要位于高地震烈度地区。随着中国经济的快速发展和西部大发展战略的实施，大量高坝也在建设，高度也在不断增加。其中，双江口主墙堆石坝高312m，双河口、古水、马集坝高近300m。如果这些高坝在地震中受损，其结果将是一场灾难。因此，有效的抗洪措施是确保水库地震稳定性的必要条件。 1 坝体稳定分析 振动台模型试验结果表明:土石坝由于坝体结构对地震波的放大效应, 坝体上部加速度反应往往较大, 若大坝遭遇强震而发生破坏, 则破坏将首先从坝顶部开始, 坝顶部堆石体将最先失去平衡而产生松动、滑动乃至坍塌等.汶川地震中, 紫坪铺面板堆石坝下游坝坡也发生了局部的滑动.因此, 在地震区修建高土石坝时应特别重视坝顶区土体的稳定.“八五”期间孔宪京根据大型振动台模型试验结果, 提出了钉结护面板加固方案 (图1、2) , 并采用二维弹塑性有限元分析方法分析了钉结护面板加固技术对坝体沉降及坝坡永久变形的抑制效果, 计算结果显示, 采用钉结护面板抗震方案可以使坝顶区的永久位移明显减小.钉结护面板技术已经逐渐被工程单位认可, 并在新疆吉林台面板坝和糯扎渡心墙堆石坝等工程中进行了应用. 钉结护面板技术作为一种有效的抗震措施, 目前多采用工程类比进行设计, 缺乏计算方法, 不便于工程应用.稳定安全系数和基于Newmark法的滑移量作为坝坡稳定性的综合评价指标, 很容易为工程技术人员理解, 但迄今为止, 尚未见到有关钉结护面板加固方案对坝坡稳定性影响的研究成果.因此, 本文试图建立考虑钉结护面板作用的动力有限元稳定性计算方法, 在此基础上, 结合300 m级的两河口心墙堆石坝进行钉结护面板对强震区高土石坝坝坡稳定的影响研究. 2 考虑到钢筋作用的安全因素和滑动位移计算方法 2.1 有限元模型建立 土石坝动力稳定安全系数的计算方法主要有拟静力极限平衡法和动力有限元时程法.考虑地震过程中坝体应力的瞬时变化, 计算出每一时刻坝坡抗滑稳定安全系数, 称之为动力有限元时程法.动力有限元时程法可以考虑岩土材料的不均匀性以及其非线性的应力-应变特性, 从合理性而言, 其优于拟静力极限平衡法, 但目前广泛采用的动力有限元时程法没有考虑加筋的作用.本文在常规动力有限元时程法的基础上, 考虑了钢筋的作用, 如图3所示.即将堆石体和护面板用平面等参单元离散, 钢筋采用杆单元离散, 建立加固后大坝的有限元模型, 分别计算出大坝的震前应力和地震时每一瞬时的动应力.假定破坏时钢筋拉力沿滑弧的切向作用, 根据单元的静、动应力叠加结果可对大坝进行稳定计算, 其安全系数为 Fs=n∑i=1(ci+σitanφi)li+m∑j=1Fjn∑i=1τili+m∑j=1fj(1)Fs=∑i=1n(ci+σitanφi)li+∑j=1mFj∑i=1nτili+∑j=1mfj(1) 式中:ci、φi分别为第i单元土体的凝聚力和内摩擦角,li为第i单元滑弧面的长度,fj为滑弧穿过的第j个钢筋的拉力,Fj为滑弧穿过的第j个钢筋的抗拉力,σi、τi分别为第i单元滑弧面上法向应力和切向应力, 由下式表示: σn=σx+σy2-σx-σy2cos2α-τxysin2α(2)τ=σx-σy2sin2α-τxycos2α(3)σn=σx+σy2?σx?σy2cos2α?τxysin2α(2)τ=σx?σy2sin2α?τxycos2α(3) 式中:σx=σsx+σdx,σy=σsy+σdy,τxy=τsxy+τdxy;σsx为单元的静水平应力,σdx为单元的动水平应力,σsy为单元的静竖向应力,σdy为单元的动竖向应力,τsxy为单元的静剪应力,τdxy为单元的动剪应力. 稳定计算中, 在每一个时刻均采用枚举法根据单元应力自动搜索最危险滑弧, 可以考虑最危险滑弧随时间的变化, 这种方法比通常用拟静力法确定最危险滑弧更为合理和精确. 2.2 滑动体滑动角加速度 对于任意滑弧 (如图3所示) , 可通过下式计算滑块绕圆心的滑动角加速度: 当某时刻出现瞬时滑动时, 滑弧的滑动量为 在整个时间段里可能出现多次瞬时滑动, 则累计滑动量为 D=n∑i=1Di(7)D=∑i=1nDi(7) 3 计算模型 3.1 大坝网格、加筋 本文采用了两河口水电站心墙堆石坝作为计算算例, 大坝最大坝高为295 m.图4为采用钉结护面板加固后的大坝网格图, 图5为大坝加筋和护面板示意图.其中钢筋采用杆单元模拟, 沿坝高方向和轴向间距均为4 m. 3.2 等效静力计算 静力计算时, 筑坝材料采用邓肯E-B模型, 计算参数见表1.动力计算时, 筑坝材料采用沈珠江等改进的等效线性模型, 计算参数见表2.钢筋和护面板均采用线弹性模型, 见表3.静力计算考虑了大坝填