高面板堆石坝3维真非线性地震反应分析方法.docVIP

PAGE  PAGE 9 高面板堆石坝三维真非线性地震反应分析方法 及模型试验验证基金项目：“九五”国家科技攻关项目(96-221-02-03-02)和国家电力青年基金项目(97Q02和SPQKJ06) 作者简介：赵剑明（1970—），男，山东肥城人，博士/博士后，高级工程师，岩土及土工抗震研究 * 赵剑明 汪闻韶 常亚屏 (中国水利水电科学研究院岩土工程研究所，北京，100044) 摘要 基于所建立的土石料三维粘弹塑性动力本构模型，并采用新型三维各向异性有厚度薄单元来模拟面板和堆石的接触面特性，结合考虑孔隙水压力的消散和扩散的有效应力分析方法，开发了高面板堆石坝地震反应分析的三维真非线性有效应力分析方法。通过典型坝例分析，比较了真非线性分析与等效线性分析结果的差异，表明真非线性方法较真实地反应了结构的地震反应，而且能够直接计算出坝体的残余变形，在理论上更为合理。采用模型坝土石料低应力状态下实测的静力和动力特性参数，对三维模型坝进行了地震反应分析，并与相应的大型振动台模型试验结果进行了对比分析，两者结果相当一致，验证了本研究所提出的计算方法的可靠性。 关键词：高面板堆石坝；动力本构模型；地震反应；真非线性；模型试验 1 前言 混凝土面板堆石坝以其安全性、经济性、适应性等优点受到国内外坝工界的普遍重视。目前相当多的面板堆石坝位于强震区，使得面板堆石坝抗震问题研究的重要性和迫切性越来越突出。面板堆石坝地震反应的二维分析已有了较多的计算方法和相应程序，但由于面板堆石坝多修建在狭谷之中，具有明显的三维效应，仅做二维分析是不够的，因此面板堆石坝地震反应的三维动力分析成为必要；随着计算技术和计算机的快速发展，也使得三维分析成为可能，所以近期面板堆石坝的三维动力分析引起了国内外学者的重视。 地震反应分析方法从基于的本构模型来分可分为两大类，一类是基于等价粘弹性模型的等效线性分析方法，另一类是基于（粘）弹塑性模型的真非线性分析方法。等价粘弹性模型概念明确，使用方便，在参数的确定和应用方面积累了较丰富的试验资料和工程经验，目前应用较多。然而，等价粘弹性模型的缺陷是明显的：等价粘弹性模型在加荷与卸荷时模量相同，因而不能直接计算土体在周期荷载连续作用下发生的残余变形；不能考虑应力路径的影响；不能考虑土的各向异性；较大应变时误差大，偏于不安全。由此可见，基于等价粘弹性模型的等效线性分析方法得到的地震反应并不是土体真实的地震反应，要想得到土体真实的地震反应，应采用基于（粘）弹塑性模型的真非线性分析方法。 2 高面板堆石坝三维真非线性地震反应分析方法 2.1 土石料的三维真非线性动力本构模型 “八五”期间，汪闻韶院士和李万红[1]建立了一个二维的适用于土石坝的粘弹塑性动力本构模型，初步应用于工程分析，取得了一些较好的结果。本研究根据面板坝的特点，结合试验资料，对该模型进行推广和完善，建立了适用于面板坝的三维粘弹塑性动力本构模型，并开发了相应三维真非线性分析方法。 本研究采用的粘弹塑性模型将土视为粘弹塑性变形材料，模型由初始加荷曲线、移动的骨干曲线和开放的滞回圈组成。这种真非线性模型的特点是：（1）与等效线性粘弹性模型相比，能够较好地模拟残余应变，用于动力分析可以直接计算残余变形；在动力分析中可以随时计算切线模量并进行非线性计算，这样得到的动力响应过程能够更好地接近实际情况。（2）与基于Masing准则的非线性模型相比，增加了初始加荷曲线，对剪应力比超过屈服剪应力比时的剪应力应变关系的描述较为合理；滞回圈是开放的；考虑了振动次数和初始剪应力比等对变形规律的影响。模型的数学表达式如下： 初始加荷曲线： (1) 骨干曲线： (2) 滞回圈： (3) 式(2)和式(3)中，在加荷时取（－）、（＋），在卸荷时取（＋）、（－）。 在此非线性动力模型中，骨干曲线和滞回圈的原点不断移动产生残余变形，即有： (4) 以上诸式中：τ和γ为剪应力和剪应变，τmax为极限剪应力，τmax=τf /Rf ；Rf为破坏比；τf为破坏剪应力；φ 为有效内摩擦角；σ 为有效正应力；γ0为骨干曲线和滞回圈原点相应的剪应变或称塑性剪应变；γh是以γ0为零点的剪应变；A、B为模型参数；DRSd为动剪应力比幅值；DRS动剪应力比，DRS=RS-RS0，RS=τ/σ，RS0为初始剪应力比。 模型参数A和B以及γ0可以用剪应力比控制的循环三轴试验来测定[1]，主要受振次、动剪应力比幅值和初始剪应力比的影响比较大。模型参数A和B也可由等效线性粘弹性模型参数换算近似得到，换算原则是使两变形模型的骨干曲线重合和滞回圈包围的面积相等。 2.2 单元型式 堆石坝体及地基主要采用三维八结点六面体等参单元来模拟，在边界不规则处采用六结点五面体三棱柱单元来填充。 混凝土面板可用薄板单元来模拟