羊曲面板堆石坝应力变形非线性有限元分析.docVIP

羊曲水电站面板堆石坝应力变形非线性有限元分析 任喜平 摘 要：基于邓肯E～B材料本构模型，本文建立了羊曲水电站混凝土面板堆石坝的有限元计算模型，并运用中点增量法，进行了施工期和蓄水期的应力变形非线性有限元分析，获得了不同时期坝体、面板应力变形的分布与变化规律。通过有限元分析与计算，可以看出羊曲水电站面板堆石坝是安全、可靠的。 关键词：混凝土面板堆石坝；应力变形；非线性；有限元分析 1.工程概况 黄河羊曲水电站位于青海省海南州兴海县与贵南县交界处，是一座以发电为主的大型水利水电枢纽工程，是黄河干流龙羊峡水电站上游河段规划梯级电站“茨哈水电站、班多水电站和羊曲水电站水电站”三个规划梯级电站的最下一级。水库正常蓄水位2715.0m，相应库容14.724。坝体从上游到下游依次分粘土铺盖区、盖重区、混凝土面板、垫层区、特殊垫层区、过渡区、主堆石区、下游堆石区。 2.计算原理 2.1.材料本构模型 面板堆石坝的材料主体为垫层、过渡层及堆石体等散粒体材料，在荷载作用下这些材料通常呈现非线性变形特性。目前，国内通常采用非线性邓肯E～B本构模型来模拟这些材料的应力应变关系。该本构模型的基本原理如下[2(4]： 材料的切线弹模为： （1） 式中，S为应力水平，定义为实际主应力差与破坏时主应力差的比值，即： （2） 为破坏比，其值小于1.0，定义为破坏时的主应力差与主应力差渐近值的比值: （3） Ei为初始切线模量，定义为: （4） 式中，K、 （5） 式中，Kb、 （6） 式中，Kur为材料的试验参数。 卸荷判定采用邓肯在1984年的土石坝计算程序中所提出的卸荷准则。定义应力状态函数为： （7） 式中，S为应力水平，pa为大气压力。 若以Fmax表示应力状态历史记录中的最大值，当计算的F＞Fmax时，判定为加荷情况，取用切线弹模Et值。当计算的F0.75Fmax时，认为处于卸荷状态，取用卸荷模量Eur。当0.75 FmaxFFmax时，认为处于过渡状态，按如下线性插值取用卸荷模量，即： （8） 2.2.材料非线性问题的有限元解法 目前，材料非线性问题的基本解法有迭代法和增量法两种，其中按材料切线模量迭代计算过程的不同，增量法又分为基本增量法和中点增量法两种。相比较而言，由于中点增量法具有分级迭代易于收敛、计算精度较高等优点，因此在工程实际中应用更为普遍。中点增量法的基本步骤如下[5]： （1）设第i级荷载增量为 ,根据初始应力 确定Et,i-1、νt，i-1等参数, 并进而形成刚度矩阵 。 （2）加荷载增量之半，即 /2于结构，由线性方程组解位移增量。 （9） （3）根据 ，求 、 ，并进而累加求 、 。 （4）由 或 确定Et,i-1/2等材料模量参数，再形成矩阵 。 （5）在 基础上重新加i级全荷载增量 用下式解出本级位移增量 ： （10） （6）由 求 、 ，进而求 、 ,得到施加i级荷载增量 后的结果。 （7）i＝n(荷载总级数)，就得到全部结果。 3.计算条件 3.1.计算模型 沿坝轴线选取最大剖面所在的坝段建立计算模型。 计算范围：上、下游方向分别取至上、下游坝坡坡脚，沿坝轴线方向取一个面板块的宽度即长为12m的坝段，坝底取至河床基岩面。 边界条件：底部取为固定铰约束，模型左右两侧面施加沿坝轴线方向的约束。 混凝土采用线弹性模型，垫层、过渡层及堆石体材料采用非线性邓肯E～B模型。 选用8节点等参单元进行剖分，共剖分为102个单元、596个节点。计算的几何模型和有限元网格剖分结果分别见图2、图3。 图2 大坝几何模型图 图3 有限元网格图 3.2.计算参数 (1) 水库特征水位 正常高水位：1856.00m 死 水 位：m 下游水位：m (2) 混凝土面板按照线弹性材料计算，计算中取材料密度=2450kg/m3，泊松比=0.167，弹性模量=2.0×