沥青混凝土心墙堆石坝的地震反应特性分析.pdfVIP

现代水利水电工程抗震防灾研究与进展 沥青混凝土心墙堆石坝的地震反应特性分析 李守义马成成李炎隆涂幸 (西安理工大学水利水电学院西安710048) 摘要：基于土石坝动力分析中广泛采用的等效非线性本构模型，本文g习-ADINA软件进行了二次开 发，实现了土石坝的动力等效线性化模拟。结合某沥青混凝土心墙堆石坝工程，采用三维非线性 有限元法，进行了地震反应分析，重点对沥青混凝土心墙的绝对加速度反应、动主应力反应及动 剪应力反应的分布规律进行了分析。结果表明：沥青混凝土心墙项部及坝顶附近区域的加速度反 应比较大；由于受地形约束条件影响，沥青混凝土心墙在接近坝基和心墙2／3高度处与岸坡混凝 土垫座相接部位动应力较大，出现了拉应力区域；心墙的动剪应力最大值发生在心墙1／2高度处 与岸坡混凝土垫座相接部位，受地形约束条件影响，心墙靠近两岸混凝土垫座附近出现局部剪应 力集中现象。本文的研究成果为该工程的设计提供了参考依据。 关键词：堆石坝；地震反应；沥青混凝土心墙；等效非线性；有限元分析 0引言 土石坝具有就地取材、造价低、结构简单、抗震能力好等特点，在世界范闱内被广泛采用…。我国地 处环太平洋地震带和欧亚地震带之间，是一个多地震国家，地震区域广阔而分散，而且地震频繁而强烈， 主要发生在中西部及西南地区比划。随着国家西部大开发及大力发展水电战略的实施，在我国西南地区修 建的高土石坝工程大多位于强震区，因此，工程的地震反应特性是备受关注的关键技术问题b1。 目前，土石坝地震反应特性分析主要采用的是拟静力法和地震反应时程分析法。拟静力法既不考虑地震 特性，如振动频率、次数和地震持续时间等，又不考虑坝体材料动力性质和阻尼性质等，故无法反映坝在地 震时的反应特性u61。而时程分析法可以校正由于采用反应谱法振型分解和组合求解结构应力和位移时的 误差；可以计算结构在非弹性阶段的地震反应，对结构进行大震作用下的变形验算，从而了解结构的薄弱 部位，以便采取适当的工程措施；可以计算结构和各结构构件在地展作用下每个时刻的地震反应。因此， 时程分析法的计算结果能更真实地反映结构的地震反应，从而能更准确地揭示结构的薄弱部位，在土石坝 地震反应分析中得到了广泛的应用哺。。本文基于土石坝动力分析的等效非线性本构模型，采用大型有限元 软件ADINA，对某沥青混凝土心墙堆石坝进行了三维地震反应分析，重点研究了沥青混凝土心墙地震反应 特性，为该工程的设计提供了参考依据。 1 土石坝地震反应分析的基本原理 在土石坝地震反应分析中，坝体及坝基可看作由许多质点组成的多自由度体系，经过有限元离散后， 其动力方程表示如下： [M]{艿)+【c】{艿)+[K】{6)=一【M】{疋 (1) 式中：【M卜．表示坝体及坝基的总质量矩阵，【q一表示坝体及坝基的总阻尼矩阵，【明一表示坝体及坝基的 总刚度矩阵；p}一表示坝体及坝基节点加速度列阵，耍酉一表示坝体及坝基节点速度列阵，pj一表示坝 体及坝基节点位移列阵；{疋_}为地震加速度列阵。 质量矩阵【M】可用集中质量法求得，即假定单元的质量集中在节点上；刚度矩阵∞可用常规的有限元 法求得；总阻尼阵【q由各单元阻尼矩阵组成，单元阻尼矩【q。由下式计算： [c】8=％[M】。+屈隧】。 (2) 式中：系数瑾一纠，，6r：,^／w，，其中九为单元的阻尼比，叫1为土体振动基频，可按下式计算： 现代水利水电工程抗震防灾研究与进展 ([K卜602[M】){皖)=0 (3) 式中：{哦)是自振时坝体各单元节点的振幅。 求解时，考虑堆石坝体和覆盖层地基的动力非线性，动剪应力模量6和阻尼比与剪应变y有很大的关 系。剪切模量6采用公式： (4) G=KxPall仃】：：×G／G一 式中：K、，I一相当于6懈时的试验常数；‰一平均有效应力；6／6m毅一为剪应变y的函数。