水利工程论文-坝基深层软岩对土石坝防渗系统的影响摘要：本文采用邓肯模型对西霞院土石坝的变形进行三维非线性有限元分析计算，分析坝基深层存在软岩对土石坝防渗系统的影响。关键词：土石坝软岩三维非线性有限元变形影响土石坝对地质条件要求相对较低，在深覆盖层上直接修筑土石坝可以减少覆盖层开挖，节省工程投资。目前对建造在深厚覆盖层上的土石坝模拟计算中，坝基影响范围只取至覆盖层底部，覆盖层以下的基岩作为刚性体考虑，但实际工程情况并不完全相同。当覆盖层以下是变形模量很小的软岩时，坝体以及防渗结构可能产生更大的沉降和变形,这将影响防渗系统的可靠性。由于国内很少有这方面的工程实例介绍，故本文采用西霞院河床段土石坝作为计算模型，应用邓肯模型对坝体、坝基进行三维有限元计算，分析坝基深层存在软岩对土石坝坝体变形的影响以及对防渗系统的影响。1工程概况西霞院工程建于小浪底大坝下游16公里处的黄河干流上，坝址位于吉利区和孟津县白鹤镇之间，属于低山丘陵区。工程的开发任务是以反调节为主，结合发电，兼顾灌溉、供水综合利用。工程布置采用土石坝，坝长约2600m，坝体采用砂卵石坝壳复合土工膜防渗斜墙，坝基防渗采用砼防渗墙。坝址地形呈宽阔"U"型河谷，宽约3km。河床两侧分布高漫滩及Ⅱ级阶地。坝址处两岸河漫滩表层为2m～7m的新近沉积的砂壤土、砂层等，结构松散，靠近两岸Ⅱ级阶地其上部多分布有粉质壤土。河漫滩下部为砂砾石层，一般厚20m～28m，中等密实状态。砂砾石层下部基岩由上第三纪砂岩与粘土岩互层组成，饱和抗压强度平均值为4.3MPa，属极软岩类岩体。2计算模型及材料参数2.1计算模型西霞院土石坝坝顶高程139.0m，最大坝高21.0m，坝顶宽8米，上游边坡1：2.75，下游边坡1：2坝，上游设壤土铺盖与斜墙相连进行防渗。计算工况：上游水位134.75m，下游水位126.23m。取100m河床段作深部软岩对坝体变形的影响分析，河床段剖面示意图见图1。图1河床段最大剖面图（单位：m）2.2计算参数选择坝基由表部松散粉细砂、中部砂砾石层及下部上第三系基岩地层组成。表部松散粉细砂，厚3.0m，粉细砂上部为砂壤土，厚1.0m，下部为砂层，厚2.0m，结构疏松；地层中部为砂砾石层，厚20m；地层下部为软岩，厚30m，弹性模量只有500Mpa（属于极软岩）。坝体及坝基E～B模型计算参数见表1。表1E～B模型计算参数参数材料土粒比重KN/m3CKN/m3Φ(O)KnRfKbmKurnur坝体压实砂砾2.6621.00325500.450.802200.511000.45地基砂砾石Q42.6621.00325000.480.822000.410000.48地基砂砾石Q32.6821.00326000.480.812400.512000.48砼防渗墙C1524.0E=22000MPa,ν=0.167软岩17.1E=500MPa,ν=0.3接触面参数CKN/m3Φ(O)K1n‘Rf’防渗墙与砂层018300001.00.9防渗墙与砂砾层022360001.00.92.3单元剖分及计算坐标空间模型坐标系建立如下：顺河向为X轴，指向下游为正；竖直向为Y轴，竖直向上为正；坝轴线为Z轴，指向河谷中心（右岸）为正，指向左岸为负。土坝的有限元计算网格见图3.2。整个结构共剖分单元5208个，结点总数6315个。图2有限元计算网格图2.4荷载分级情况在有限元模拟分析中，结合施工顺序共分十级加荷。第一级为坝基砂卵石层；第二级为河槽段截流围堰体；第三级为围堰防渗墙及第一级坝身；第四及至第七级为坝身；第八级分别为复合土工膜及上游面混凝土护坡；第九级为上游校核洪水134.75m及相应的下游水位126.23m。3计算结果及分析3.1对坝体变形的敏感性分析考虑软岩情况与不考虑软岩情况下坝体的最大水平位移值及沉降值见表2。在两种情况下水平位移、最大沉降的绝对值均不大，考虑软岩与不考虑软岩的水平位移最大值变化微小，但考虑软岩影响情况下最大沉降值比不考虑软岩影响情况大2.5cm，坝体最大沉降值增大约5.5%。从水平位移最大值和沉降最大值以及坝体水平位移及沉降变形图可以看出：软岩的存在增加坝体变形，但对土石坝本身来说影响不大。表2考虑软岩与不考虑软岩坝体最大位移比较表位移工况向上游最大位移值(cm)向下游最大位移值(cm)最大沉降值(cm)考虑软岩影响-10.622.6-45.1不考虑软岩影响-10.222.4-42.63.2对坝坡面复合土工膜小主应变敏感性分析上述讨论的变形值都是...