纵缝对重力坝地震发硬影响的研究.docVIP

纵缝对重力坝地震反应影响的研究 涂 劲，候顺载，陈厚群(中国水利水电科学研究院 工程抗震研究中心) 摘 要：本文采用显式有限元结合人工透谢边界进行波动分析，配合动接触力模型考虑纵缝接触非线性问题，用静动组合方法考虑静动荷载的耦合作用，形成新的重力坝地震反应分析的计算体系，对一典型的含纵缝重力坝进行了地震反应分析。在考虑地基辐射阻尼的条件下得出了纵缝的存在及纵缝初始宽度、填充情况、缝面摩擦系数等因素对坝体地震反应影响的一系列结论。 关键词：纵缝；重力坝；人工透射边界；动接触力法收稿日期：1999-11-22 作者简介：涂劲(1973-)，女，福建长汀人，博士、工程师，主要研究方向：水工结构抗震研究。 　　重力坝是水电建设中最常见的坝型，其特点是断面尺寸巨大，与地基接触面广，坝内应力较低、较均匀，耐久性好，所以对抵抗地震、渗漏、洪水及战争破坏等意外荷载的能力都较强。重力坝是大体积结构，为了防止开裂和满足施工要求，需要分缝分块浇筑，沿坝轴线方向设置的横缝可以通过计算中以单个坝段为研究对象加以考虑，而当由于浇 筑能力限制和施工速度要求而采用柱状块方法浇筑时，坝的横断面上还将有纵缝存在，通过纵缝灌浆并不能使坝体完全恢复整体作用，纵缝将削弱坝体刚度，并产生不利的应力分布。文献[1，2]分别以点面接触的物理模型和将接触点对的接触条件引入基本方程的有限元混合法对含纵缝的重力坝进行了地震反应分析，但计算过程都较为繁琐，且对于边界条件分别以刚性地基和无质量地基来简化，未考虑地基的辐射阻尼作用。本文以显式有限元结合人工透射边界的波动分析方法，配合动接触力模型考虑纵缝接触非线性问题，用静动组合方法考虑静动荷载的耦合作用，形成新的重力坝地震反应分析的计算体系，对一典型的含纵缝重力坝进行地震反应分析。 1 计算方法 1.1 无限地基中的波动模拟方法 本文采用廖振鹏[3]建议的多次透射人工边界方法模拟地震波向无穷远处辐射，并与显式有限元的内点计算相结合，形成了完全解耦的近场波动数值模拟方法。在有限元分析中，首先将结构和感兴趣的部分地基切割出来进行离散化，成为近场计算区，在近场计算区的地基外附加一定厚度线性的人工边界区，人工边界区的外部节点称为人工边界节点，人工边界以内的所有节点统称为内点。计算模型如图1. 　 　　人工边界节点的位移包括自由场与散射场两部分，自由场含入射场和均匀无限半空间的自由面反射场，是已知的，而散射场位移可由透射人工边界条件给出，N阶透射公式为： (1) 式中：CNn=N!/(N-n)!n!. 　　式(1)具有解耦特征，为建立一个完全解耦的近场波动计算模型，内点采用集中质量有限元法解耦，显式有限元积分格式如下(角标ij表示向量在i节点j自由度方向的分量)： (2) ijt+Δt=2/Δt(uijt+Δt-utij-tij (3) 　　显式有限元结合人工透射边界的波动模拟方法具体计算步骤为：(1)由内点计算公式计算有限元计算区内各节点t+Δt时刻的位移反应；(2)计算t+Δt时刻人工边界区各点自由场位移；(3)对人工边界区，将t+Δt时刻以及相关前几时刻的总位移减去相应自由场位移，获得多次透射公式外推t+Δt时刻散射场位移所需人工边界区节点前几时刻散射场位移；(4)用透射公式计算t+Δt时刻人工边界节点的散射场位移；(5)将人工边界节点t+Δt时刻的散射场位移与自由场位移相加，得到总位移，进入下一时步。 1.2 纵缝开合的动接触力计算方法 在内点计算中，对于接触点，将uijt+Δt分为三部分： uijt+Δt=uij-t+Δt+Δuijt+Δt+Δvijt+Δt (7) 其中，uij-t+Δt可以由前一时刻的运动状态直接得到。 (8) 而后两项由接缝的动接触状态确定。 Δuijt+Δt=Δt2/2miNtij=Ntij/Mi (9) 式中：Mi=2mi/Δt2. Δvijt+Δt=Δt2/2miτtij=τtij/Mi (10) Tu(ui′t+Δt-uit+Δt)=0 (11) 令Δ1i=Ti(i′t+Δt-it+Δt，由式(7)和式(11)可得： Nti=MiMi′/Mi+Mi′iΔ1i (12) 令Δ2i=Tt[(i′t+Δt-it+Δt-ti′-ti)],在切向根据静摩擦定律可得： τti=MiMi′/Mi+Mi′iΔ2i (13) |τti|=μ|Nti| (14) 式中：μ为动摩擦系数；τti的符号由Δ2i确定。 　　最后由式(10)计算Δvit+Δt、Δvi′t+Δt。 　　由以上计算求得it+Δt、Δuit+Δt、Δvtt+Δt，从而得到接触点总位移uit+Δt，这就是动接触力模型计算格式，它与显式有限元的内点计算配套，且不存在人为选扔触刚度的问题，不会发生接触面的相互侵入现象。