面板堆石坝面板地震反应分析ChinaInstituteofWater.docVIP

面板堆石坝面板地震反应分析-China Institute of Water 紫坪铺曲板堆石坝地震永久变形分析 赵剑明常亚屏陈宁 （屮国水利水电科学研究院岩土所，北京， 100044） 摘要 针对紫坪铺面板堆石坝工程，根据堆石料的大型动三轴试验结果，确 定了坝料的残余体应变和残余剪应变模式，特别是非饱和料的残余体应变 模式。在所建立的面板坝三维非线性动力反应有限元法基础上，结合孔隙 水压力扩散和消散的计算，建立了一套同时计入残余体应变和残余剪应变 的面板堆石坝地震永久变形的计算方法，并对紫坪铺面板坝工程进行了地 震永久变形的计算，分析了坝体地震永久变形的量值和分布情况，为大坝 的抗震设计提供了冇力的技术依据。 关键词 面板堆石坝，地震反应，残余体应变，地震永久变形 1前言 紫坪铺水库工程位于四川省成都市西北60多公里的岷江上游，是一 座以灌溉和供水为主，兼有发电、防洪、环境保护、旅游等综合效益的大 型水利枢纽工程，是都江堰灌区的调节水源工程。本工程位置在都江堰市 麻溪乡，上游与岷江干流映秀湾电站尾水衔接，下游距都江堰市约9公里。 该工程挡水建筑物为钢筋混凝土面板堆石坝，最大坝高156m,坝顶 高程884.00m,坝顶长度663.77m,上游坝坡1:1.4,下游坝坡1:1.5~1:1.4, IE常蓄水位877.0m,设计洪水位871.2m,校核洪水位883.1m,总库容11.12 亿m3。坝基覆盖层由漂卵砾石夹砂层透镜体组成，厚10^23m, 1 结构松散，具一定架空结构，但仍具有一定抗压缩变形能力，抗剪强 度也较高，渗透系数在5.79?10-11.57?10cm/So水库工程定为一等工程， 主要建筑物均为I级建筑物。设计洪水频率为千年一遇，校核洪水为可能 最大洪水。本工程区基本地震烈度为VII度，面板坝按VIII度地震烈度设防。 受四川省水利水电勘测设计研究院委托，中国水利水电科学研究院岩 土工程研究所通过室内大型动三轴试验，详细研究了紫坪铺水库曲板堆石 坝坝料的动力工程性质，其中重点研究了坝料的地震残余变形特性，结合 在所建立的面板坝三维非线性动力反应有限元法，对大坝遭遇8度地震时 的地震永久变形进行了分析。 2面板坝地震永久变形的三维非线性有限元分析方法 2.1地震永久变形研究概况 作为大坝抗震安全性评价的的重要指标，面板堆石坝的地震永久变形 计算是面板堆石坝抗震研究中一个重要课题。土石坝地震永久变形的计算 方法除了利用弹塑性模型直接计算残余变形的真非线性分析方法外，主要 有两大类⑴：一类是滑动体位移分析法，一类是整体变形分析法。滑动体 位移分析法的基本岀发点是假定土石坝的残余变形主要是由地震时坝坡 及地基发生瞬态失稳时滑移体产生位移造成的，主要基于Newmark的刚休 滑动面假设和屈服加速度概念。整体变形分析法的基本假定是把坝体及地 基作为连续介质来处理。这类方法一般是先进行地震反应分析，求岀坝体 及地基的反应，然后利用材料动力特性的试验结果，加以简化求出坝体残 余变形。目前，基于Serff和Seed等人提出的应变势概念的整体变形计算 方法有着广泛的应用。 在现有的土石坝地震永久变形的计算方法中，多数只考虑了土料的残 余剪应变，而没有计入土料的残余体应变。面板堆石坝中堆石体具有由于 棱角破碎而引起体积收缩的特性，残余体积变形不宜忽略。冃前对饱和土 料的残余体积变形研究较多，而对非饱和土料的残余体积变形研究则很少, 而面板堆石坝在运用期间大部分堆石体处于非饱和状态，因此为了较准确 地计算血板堆右坝的地震残余变形，基于试验资料建立非饱和堆石料的残 余体应变模式是非常重要的。 基于应变势概念的整体变形计算方法的基本过程是首先对大坝进行 地震反应分析，确定大坝在地震作用下的应力应变状态，然后利用所建立 的残余应变模式计算残余应变，并根据应变势的概念计算大坝的地震永久 变形。因此要准确计算面板堆石坝的地震永久变形的两个关键问题是：建 立可靠的面板坝地震反应分析方法和堆石体残余应变计算模式。2.2面板 坝地震反应分析方法2.2.1本构模型 本文计算中土石料的静应力应变关系采用了邓肯E-B非线性模型。 为了反映土的动应力应变关系的两个基本特征：非线性和滞后性，并 表示出剪切模量和阻尼比与动剪应变幅的关系，本文采用了中国水利水电 科学研究院大型动三轴试验结果。根据试验测得动剪切模量比G/Gmax和 动阻尼比X与动剪应变丫的关系曲线，采用参考剪应变Y r= t则得到较为 单一的max/Gmax归一后， G/Gmax-y/yr关系曲线和入~丫/丫「关系曲线，图1是本文 ?2?2 所研究的面板堆石坝工程的主堆石料的相应曲线。动力计算时输入相 应关系曲线的控制数据，根据应力应变值进行内插和外延取值，用于计算。 根据试验结果，最大动剪切模量Gma