强震区堆石坝应力变形的DDA法分析.pdfVIP

第47卷 第2期 人 民 长 江 Vo1．47，No．2 2016年 1月 Yangtze River Jan．， 2016 文章编号：1001—4179(2016)02—0072—06 强震区堆石坝应力变形的DDA法分析 刘 勇，王 运 生，陈 怡 西，唐 起 (成都理工大学 地质灾害防治与地质环境保护 国家重点实验室，四川 成都 610059) 摘要 ：基于非连续变形分析方法 (DDA)，针对堆石坝 中堆积颗粒 的应变特性，对各部分程序建立或改进 了接 口，实现了堆石料的随机分布，并且能够满足级配要求；同时也实现了堆积颗粒在筑坝中逐层堆积、逐层加载 的过程。探索了DDA中弹簧刚度和步位移比的最佳取值范围。在此基础上，分析 了堆石坝施工、蓄水及地震 作用下的应力变形特征。结果表明：施工过程中坝体的沉降量最大；蓄水后，坝体两侧水平位移变化明显，出 现 向下游偏移 的趋势；在地震波作用下，坝体水平 向和竖向加速度峰值随高程增加而增大。符合一般堆石坝规 律，表明改进后的DDA程序具有较好适用性。 关 键 词 ：DDA；试验模型；应力；变形 ；地震波；加速度 ；堆石坝 中图法分类号：TV641．1 文献标志码 ：A DOI：10．16232／j．cnki．1001—4179．2016．02．017 自20世纪80年代以来 ，堆石坝由于其结构简单 、 连续介质提供 了很好的手段 。本文结合地震波的 施工速度快及适应性强等优势在我国发展很快 。 传播特点，针对 DDA程序的不足之处，对程序进行了 目前在建和已建成的堆石坝 已超过 200多座，其 中有 部分改进、拓展，并将其运用于堆石坝 的应力变形分 相当数量的堆石坝建在地震烈度较强的地区，如四川 析。 岷江紫坪铺堆石坝等。近年来，随着我国对西部水力 1 DDA方法基本理论 资源的大力开发，堆石坝作为一种推广坝型，其抗震问 题更加突出，尤其是在汶川8．0级地震、芦山7．0地震 DDA方法是平行于有限元法的另一种数值分析 后 ，对堆石坝的抗震稳定更为关注。但我 国堆石坝的 方法 ，与有限元法区别在于它以计算不连续面的滑移 、 发展起步较晚，实际震害和记录资料较少 ，特别是缺乏 开裂和旋转等大位移的静力和动力问题为主要优势。 高地震烈度区堆石坝的地震损伤资料。 目前，国内外 该方法有如下理论基础 一 ：计算几何模型 由被节理 在采用有限元法及等效线性模式等动力数值模拟的方 等结构面切割而形成的块体单元组成；块体单元的运 法对堆石坝开展动力计算方面取得了长足的进展 。 动和变形 由切 向应变 、法 向应变 、位移和转动组成 ， 有限元法及等效线性模型方法是基于连续介质发 DDA方法 以各块体单元的位移量为未知量 ；块体运动 展起来的数值分析方法，而堆石料作为堆石坝的筑坝 过程 中各单元之间可以接触，也可以分离，但必须满足 材料 ，是 由坚硬颗粒组成的散粒集合 ，属于典型的非连 库仑准则和块体间无嵌入 、无拉伸条件 ；结合边界条 续介质，相应堆石坝的分析理论上应该用非连续介质 件，通过最小势能原理把块体间的接触问题和块体本 的方法。1986年 ，石根华提 出了非连续变形分析方法 身的变形问题统一到矩阵中求解，得出总体平衡方程 ； (DDA) ，它不要求相邻块体单元的顶点重合，分析 采用基于时步的隐式方法求解方程组 ，并得 出块体系 过程和实际过程很接近 ，为研究如块状岩体之类的非 统平衡时的位移场、应力场等，以及各块体在运动过程 收稿 日期 ：2015—06—03