复杂地形条件下高面板堆石坝应力变形的三维有限元分析.docxVIP

复杂地形条件下高面板堆石坝应力变形的三维有限元分析 自20世纪80年代以来，中国的混凝土板堆填技术取得了很大进步，无论水库数量还是水库高度，都处于世界领先地位。然而，目前混凝土板堆填的设计和施工水平仍然取决于工程经验，尤其是在特殊地形和地质条件下的高板堆填石坝。目前还缺乏一套完整的工程经验。理论研究成果已成为项目的重要依据。 某工程可行性研究阶段的推荐坝型为混凝土面板堆石坝.大坝高107 m,上游坝坡为1∶1.4,下游坝坡按变坡设计,综合坝坡1∶1.487.大坝坝址处岸坡陡峻且起伏变化大,岸坡平均自然坡度约为1∶0.2.在坝轴线处,两岸大部分岸坡几乎直立且上部岸坡突然向两岸拓宽,河谷大致呈对称的复式矩形断面.在面板部位,左岸岸坡几乎直立,右岸为一斜陡坡,且在上部有一水平台地.整个大坝所处河谷,从上游到下游均有多处较大起伏或突变.目前国内外尚未见如此复杂地形条件下的已建面板坝.为此,本文运用三维非线性有限元法分析了这种“恶劣”地形条件下高面板堆石坝的应力变形特性,并研究了复杂地形边界对大坝应力变形可能产生的不利影响及相应对策. 1 堆磨料应力应变模型修正 堆石料的变形不仅随荷载大小变化,而且还与加载的应力路径有关,应力应变呈明显的非线性关系.除此之外,堆石料应力应变关系还具有压硬性、剪胀性和各向异性等特征.目前,堆石料常用的本构模型主要有弹性非线性模型和弹塑性模型2类.其中弹性非线性模型中的邓肯双曲线模型因使用简便,参数确定经验丰富,计算结果与实测结果较为符合,在工程界得到了普遍应用.文献建议在三维计算时应对该模型进行修正,即用平均主应力p代替小主应力σ3,广义剪应力q代替主应力差σ1-σ3,以考虑中主应力σ2的影响.本文采用按上述方法修正的邓肯E-B模型来描述堆石料的应力应变关系. 混凝土面板与坝体垫层料性质相差悬殊,因此可在两者之间设置无厚度的Goodman接触单元,以模拟两者之间的相互作用.此外,为了反映面板垂直缝和周边缝中止水材料的力学作用,计算中用六面体连接单元进行模拟.本文计算中大坝各分区堆石料的材料参数、接触面参数及止水材料连接单元模型参数均借用天生桥一级面板坝工程的相关参数. 混凝土材料采用线弹性模型,材料常数取为:ρ=2.45 g/cm3,E=20 GPa,μ=0.167. 2 三维非线性开采的计算与分析 2.1 单元剖分和加载 计算域边界选取时完全模拟了现有的复杂地形条件.整个计算域共剖分结点2 455个,单元1 871个,其中,面板单元和接触单元均为137个,面板垂直缝单元126个,周边缝单元27个,剖分的网格如图1所示.计算过程详细模拟了坝体的分期分断面施工、面板的分期浇筑和库水位的分期蓄降,加载次序如图2所示. 2.2 大坝堆石体应力 竣工期大坝堆石体向上下游最大水平位移分别为19.2cm和27.7cm,最大沉降为67.6cm,最大沉降点大约位于坝高2/3处,沉降率为0.64%,这些数值均在堆石坝变形的合理范围之内.坝体向上游的水平位移明显小于向下游的水平位移,这主要是由于先期填筑的一期“子坝体”受库水预压后变形模量有所增大且对后期坝体向上游的变形有一定限制作用,以及坝体次堆石料较主堆石料“软”的缘故.也正因如此,坝体的同一沉降值,其等值线范围也偏向于下游次堆石区.蓄水后,堆石体水平位移随库水位的增加而逐渐向下游发展.向上游的水平“鼓出”位移被库水逐渐压回,向下游的水平位移也有一定的增加.向上下游最大水平位移分别变为7.3cm和29.2cm,其中向上游的位移改变量明显比向下游的位移改变量要大,这是由于库水通过面板主要作用在上游坝体所致.同样,库水引起的坝体沉降变位增量也主要集中在上游主堆石区内,但由于上游主堆石模量较大,坝体竖向位移因蓄水引起的增量并不大,最大沉降仅增至70.0cm.图3给出了大坝沿坝轴线纵剖面堆石体在竣工期和蓄水期的沉降等值线.由图3可见,大坝堆石体的沉降变形规律较好,等值线分布与河谷轮廓相适应,最值区域位于坝体的中上部.两岸直立岸坡对堆石体变形影响不是很明显,但河床底部的突起对沉降等值线分布影响较大,等值线上凸明显,说明该部位以上堆石体与两侧堆石体变形有较大的沉降梯度,这对坝体受力不利. 大坝堆石体大小主应力分布及变化规律总体较好,大小主应力极值之比也符合堆石体应力分布的一般规律.竣工期大坝上下游侧的主应力等值线总体上接近平行于坝坡,说明堆石体应力主要与“土柱”质量(高度)相关.水库蓄水后,上游堆石体主应力等值线较竣工期有明显“上抬”,同一数值等值线的包含范围也有所增加,这主要是水压力作用的结果.计算表明,堆石体应力水平小于1.0,个别部位尤其是边界突变处,应力水平接近于1.0,说明堆石体在这些部位可能会出现破坏区.因此,施工前最好将突变的地形进行削坡和填塞处理,使建基面能够平缓变