考虑分缝的重力坝的非线性反应分析.pdfVIP

了f诤．： 考虑分缝的重力坝的非线性反应分析 周迎新扇晶 (大连理工大学土建学院) 匠摘要j在常规的重力坝设计计算中．常将埂体作为一个整体．实际上，在每个坝段都有纵缝存在．这些 纵缝常常是坝体中的薄弱面，对结构的整体性削弱作用是非常明显的。困此，正确有效地模拟分缝对大坝地 震反应的影响是大坝抗震安全的重要问题。奉文采用点面接触的物理模型，接触本构关系采用库伦摩擦定律， 利用修正的势能原理导出系统的接触控制方程。根据遮一思想过程，编制了有限元分析程序，并对一个实际 工程进行了分析。结果表明．本文的方法可以较好地模拟坝体在地震作用下的反应，是一种合理有效的方法。 关键词：分缝结构：点面接触模型；地震反应：非线性分析 (一) 一般约定与假设 在此采用的是点面接触的物理模型”¨”。如图l，将其中一物体定义为主动接触体，另 ～物体定义为被动接触体。分别简称主动体和被动体。图1只显示了两个物体相接触的情况． 同样可以推广到多个物体柏接触的情形。在建立接触控制方程时，采用如下基本条件与假设： (1)主动体上的任何节点不允许在被动体内，即不允许有材料的嵌入。(2)两接触体在接触 面上的接触力大小相等，方向相反。(3)仍采用库仑摩擦定律，但是静与动摩擦系数可采用 不同的值。 (=)接触边界的协调条件 对坐标和位移均采用的是线性插值．局部坐标的原点在一个接触段的中点．主动体和 被动体分别在接触面的两侧，两物体接触面单元的法线的标量积为负。可能出现三种接触状 态：(1)粘式：两物体接触的点无相对运动，变形前后接触点的局部坐标值相同。(2)滑移： 两物体相接触，但变形后接触点间沿接触面切面有相对运动，沿接触面法线方向两接触面坐 标相同。(3)开式：两物体表面某些部位并未接触，此时接触约束释放。 在一个接触面对中，不同部位会有不同的接触状态，故有不同的变形一致条件。主要 讨论粘式条件，滑移条件仅是对粘式条件作了简单的修正，开式状态无接触变形的协调条件， 不必讨论。从图2可见，(卜I)次迭代结束时主动体节点t嵌入了被动体中．故有了材料重叠 避。’，违反了接触边界变形协调条件，应进行修正。对于粘式接触，修正的结果是在(卜【)次 迭代后被动体表面被接触点c应与主动体接触点^共同运动到r”处，即点C与点}在(，)次迭 代结束时位置重合。节点k接触力记为““砟“’。在图2中．可得到如下的几何关系： 一霹m=}+“世。，“—彬，’r (I) 一群“：一罐也-+叫一， (2) 秽m=矽卜’《m～“a罐叫 (3) t一，卢∥，：如．砖r L4) 636‘：“ }L} (5) fl(c-D=彩f“毋虬H6竹u]=衫肝…霹哪一，+n鹆u)_一驾w] 可以得到粘式情况下在接触边界上变形一致条件： ’+“群“=…q“1-J+叫“、(6) 对于滑移情况，沿被动段法线方向主动节点七和被动点f的坐标应相同，而切线方向有相对 运动，其约束为：一…《“=ni(”《“一m雄“) (7) 其中，n，为被动段J的法向向量。即仅有一个自由度受变形一致条件约束。可将式(6)，(7) 合并写成： (、j“4《w=(了●“《“”一…q：w) 式中cJ为一个参数t粘式接触时t=f，】，滑移接触时t、，=q。开式情况无变形一致条件。 (三) 约束变分原理 考虑接触边界条件时的虚功原理的表达式为系统的总虚功取极值 踟=』q曲。dV一』f西～dV—It,s,，dS=()． 9) I’