渗流场和应力场耦合应力变形数值分析.pdfVIP

水电站设计 第28卷第3期 D H P S 2012年9月 渗流场和应力场耦合应力变形数值分析 郝春华 ，龚 霞 (1．重庆市交通规划勘察设计院，重庆 401121；2．林同梭国际工程咨询 (中国)有限公司，重庆 401121) 摘 要：土坝的实际运行和理论研究均表明，当大坝建成蓄水后，上游库水位的抬高将使水体在坝体和坝基中产生从上游到下游 的流动，形成渗流场，从而对土坝原有的应力场产生影响，且这种影响是相互的，其实质是渗流作用下渗流场与应力场的耦合作 用。本文结合具体算例阐述坝体在耦合作用下的变形情况。 关键词：土石坝；渗流场；应力场；耦合分析；应力；变形；有限元法 中图法分类号：TV223．6 文献标识码：B 文章编号：1003—9805(2012)03—0062—04 1 前 言 2 耦合模型 众所周知，在土石坝的施工过程中，不论土石坝 2．1 土体平衡方程 如何被压实，在施工期间及初期蓄水时，总是要发生 Fredlund和 Rahardjo曾对非饱和土壤介质的土 沉降变形的，因此控制好坝体的沉降是至关重要的， 体二维本构模型进行过描述，建立了土体本构关系。 尤其是在施工和蓄水同时进行的情况下坝体的沉降 相应地，对于饱和区，孔隙气压等于大气压，水体充 变化对坝体的安全具有很大影响。目前渗流场和应 满介质孔隙，可以建立如下本构关系： 力场耦合分析越来越受到关注，两者是不可分割的， {△ }=[D]{△ }+{m}Au (1) 将两者联系起来分析，能更好地反映出在坝体填筑 式中 [D]——线弹性矩阵，与弹性模量及泊松比 过程中应力、应变的动态演化。 有关 ； 渗流场和应力场的耦合分析在裂隙岩体中早已 {m}——各向同性单位张量； 成为热点问题，目前，己有较多的研究成果。对于土 △“——孑L隙水压力增量。 体而言，最有代表性的是基于Biot固结理论的渗流 据虚位移原理，将式(1)离散 ，可得如下形式： 与应力的耦合模型。但该模型无论是以各向同性土 ∑[B][D][曰]{△6}+∑{m}(N){△“}=∑F (2) 体作为分析对象，还是以各向异性饱和软土作为分 设[K]=[][D][]，[ ]={m}(Ⅳ)， 析对象，分析原理都是注重研究外荷载下孔隙水压 简写之 ，即为 ： 力与有效应力 (或总应力)及相应的变形之间的关 []{△6}[L]+[L]{△ }={△F} (3) 系。这种分析方法其实与渗透力作用无关，所探讨 式中 []——劲度矩阵； 的渗流场实质是孔隙水压力分布场 ，并不是渗流作 [Ld]——耦合矩阵。 用下所形成的渗透力场。也有学者采用耦合分析理 2．2 渗流连续性方程 论探讨渗流作用下渗流场与应力场的相互影响，但 根据研究区岩土体渗流的特点，可 以将其概化 研究成果或是针对连续介质 (岩体)的渗流行为建 为非均质各向同性渗流模型。根据达西定律，可以 立数学模型 ，或是针对均质土坝建立耦合模型 ，较 建立如下渗流方程： 少分析非均质土坝的渗流场与应力场相互作用对土 00w + + ： 0 (4) Ox