面板堆石坝地震永久变形计算方法研究.pdfVIP

一 ／ 第38卷第5期 大 连 理 工 大 学 学 报 Vo1．38，N0．5 1 3 3 8年3月 Journal of Dalian University of Technology sept．1 9 9 8 面板堆石坝地震永久变形计算方法研究’ 主 邵龙潭， 上 T ，’ (1大连理工大学土木建筑学院，大连 116024) (水利水电科学研究院抗震研究所) ． ． ． TL，二’， ·，5 摘要 针对面板堆右坝提出一种地震永久变形的计算方法，该方法 精动体 ’V f 位移分析法为基本出发点．首先利用有限单元法对面板堆石坝进行静、动力反 应分析{然后采用数学规捌法，求解地震过程每一时间步的最小抗滑稳定安全 系数及相应的滑裂面位置，从而得出整个地震过程中埂体的最危险精裂面}在 最危睦精裂面引入无厚度接触面单元，重新进行有限元动力反应计算，计算出 滑裂面各单元的地震永久变形．结果表明本方法具有合理性和实用性． 关键词．堡 要蛰丝墨毫垒奎墅 分类号 TV315．43，P315．9 随着我国高土石坝建设的发展，土石坝抗震研究从地震反应分析阶段进入到了地震永久 变形阶段．自纽马克(Newmark)首次提出应以大坝地震永久变形作为判别大坝抗震安全度的 标准以来，经过3O多年的研究，目前尚无一个较为通用的模型和定量标准，而工程界采用拟 静力法计算坝体抗滑稳定安全系数作为大坝抗震安全的指标，又存在一些不合理的因素，所 以深入进行土石坝地震永久变形的研究是非常重要的和迫切的． 1 滑动体位移分析法 滑动体位移分析法假定土石坝的地震永久变形主要是由地震时坝坡(及其地基)发生瞬态 失稳时滑动体产生位移造成的．纽马克 认为滑动体本身是刚性的，坝体的永久变形仅出现 在滑动面上．纽马克法首先确定预期滑动体的屈服加速度a，，然后计算地震所引起的预期滑 动体上的平均加速度 (t)，在已知 和aXt)时程曲线后，进行二次积分即可求出滑动体地 震永久变形‰． 滑动体位移分析法的第一个关键问题是如何确定屈服加速度 ．接纽马克公式，堆石体 的屈服加速度岛为 a．= tan(~O一 )·g (1) 式中： 为堆石动摩擦角，0为边坡坡角．而据实测值和模型试验来看，上式所给数值偏小， 因此西特(Seed)。 又根据一些堆石体的滑动破坏特征给出 — tan(尹一 )·g (2) ‘ 景 嚣蓦叠妻悬 曹 ‰ 一。 一。。 李俊杰一男．1363年生．教授，工学博士 第5期 李俊杰等：面板堆石坝地震永久变形计算方法研究 式中： 为滑动体滑坡坡角，见图1． 但按文献[3]的推导，不同形状滑动体的屈服加速度计算公式如下： d —tan(~O一 )·g+Nsin(P一卢)·g／C,cos P cos@一 ) ( ) 式中：在堆石体发生破坏时，_Ⅳ值是变化的，因 而嘶是一个变量，所以选用一恒定的。 值是不 合适的． 第二个关键问题是如何确定滑动体的滑裂面 位置．一种是用拟静力法采用滑块分析来计算求 出最小抗滑稳定安全系数K的面；第二种是采用 有限元法按平均等效动应力值与静应力值叠加按 圉1 堆石边坡动力滑移形式 滑块分析求出最小K的面；第三种可在地震每一 时间步进，行上述第二种过程的求解·因此选用哪种方