土石坝心墙水力劈裂的参数分析

土石坝心墙水力劈裂的参数分析

力裂是指由于水流的提高，在岩体和土壤之间发生和扩大裂缝的物理现象。土石坝心墙的生产力是土石坝设计中最受关注的问题之一。土石坝蓄水初期集中渗漏的主要原因之一也是水库破坏（内部侵蚀或排水）的重要因素之一。事实上,在水坝工程中土石坝的建设数量一直居于首位.我国己建成的水坝绝大多数为土石坝.目前,已建成坝高100m以上的土石坝17座,另有超过24座高100m以上的土石坝正在建设和设计论证之中,如两河口水电站土质心墙堆石坝(坝高295m)、长河水电站心墙堆石坝(坝高240m)、狮子坪水电站砾石土心墙堆石坝(坝高136m)、糯扎渡水电站砾质粘土心墙堆石坝(坝高261.5m)等.对中小型土石坝,经验比较丰富,有关理论也比较成熟.但在高土石坝设计建设中仍有许多关键的岩土工程技术问题亟待解决,心墙的水力劈裂问题就属其一.本文通过对水力劈裂发生影响较大的参数分析做数值分析,得出控制水力劈裂发生的坝壳与心墙的弹模比和心墙的泊松比的关系曲线.1不稳定区域内考虑墙内部分荷载的e-模型土体本质上是非线性材料,应遵循非线性规律.但影响土体应力应变关系的因素很多,如心墙和坝壳料的弹性模量、泊松比、心墙的坝坡、密度等,若同时考虑这些因素,不仅使解答困难,也难以分辨和估量其中某一因素的影响程度.很显然,要分析某一因素的影响,就必须假定其它因素不变,当我们考虑泊松比ν对水力劈裂的影响时,就假定弹模E以及其它因素不变;同样,当我们考虑弹模E的影响时,就假定泊松比ν以及其它因素不变.这样考虑后,土的非线性问题可简化为线性问题,讨论邓肯双曲线的E-ν模型也就成了讨论线弹性形式的E-ν模型.一般心墙材料的压缩性比坝壳的高,因此心墙竖向变形能力比坝壳大,心墙与坝壳刚度差异使得坝壳限制了心墙的这种变形,这样就产生了两者之间的“拱效应”,从而导致心墙承受的部分荷载向坝壳传递,应力也相应发生重分布,结果心墙应力下降.当心墙与坝壳刚度差异愈大,拱效应作用就愈明显,心墙应力就愈小,心墙就有可能产生水力劈裂.研究表明,心墙的泊松比和坝壳与心墙的弹模比会产生较强的“拱效应”,对水力劈裂的发生有着重要影响.因此,本文主要针对这两个参数做进一步分析.建立两个不同坝高,不同心墙坡度的大坝模型,计算中也采用不同的计算参数,将两者结果进行比较分析.第一个大坝模型的坝高为160m,心墙坡度为0.4,正常蓄水位150m,坝壳弹模Es为300MPa,坝壳泊松比νs为0.35.第二个大坝模型的坝高为200m,心墙坡度为0.6,正常蓄水位180m,坝壳弹模Es为150MPa,坝壳泊松比νs为0.35,如图1和2所示.因为心墙泊松比νc常用的取值范围为0.25到0.45,所以本文的讨论也只针对该范围,由于在设计时坝壳与心墙的弹模比一般不会太大,所以本文只在坝壳与心墙的弹模比Re为1～8的范围内讨论.2心墙泊松比与水库和心墙之间的弹模型比的值2.1坝壳与心墙弹性模量er的变化在讨论心墙泊松比νc变化时,假设坝壳的泊松比νs不变,且取值为0.35,并假设坝壳的弹性模量Es=300MPa不变,逐渐降低心墙的弹性模量,来控制坝壳与心墙弹模比eR的变化,计算得到心墙上游面孔压和竖向应力σy随深度的变化曲线,如图3所示.2.2坝壳与心墙弹模比的确定在讨论心墙泊松比νc变化时,同样假设坝壳的泊松比νs不变,且取值为0.35,并假设坝壳的弹性模量Es=150MPa不变,逐渐降低心墙的弹性模量,来控制坝壳与心墙弹模比eR的变化,计算得到心墙上游面孔压和竖向应力σy随深度的变化曲线,如图4所示由计算参数可知,坝体1的坝壳弹模Es为300MPa,坝体2的坝壳弹模Es为150MPa,在相同的坝壳与心墙的弹模比时,坝壳与心墙的弹模差是不同的,但对比坝体1和坝体2的结果,两个坝体发生水力劈裂的临界弹模比Re,cr相同,可以看出影响水力劈裂发生的是坝壳与心墙弹模比的大小,而不是坝壳与心墙弹模差的大小,说明本章选取坝壳与心墙弹模比作为研究的对象是正确的.2.3心墙泊松比与坝壳与心墙的临界弹模比关系比较图3和图4,可以看出,对于两个不同坝高,不同心墙坡度,不同心墙厚度的大坝得出了一致的结果,那就是:当心墙泊松比νc取0.25时,坝壳与心墙的临界弹模比Re,cr为3,当心墙泊松比νc取0.3时,坝壳与心墙的临界弹模比Re,cr为4,当心墙泊松比νc取0.35时,坝壳与心墙的临界弹模比Re,cr为5,当心墙泊松比νc取0.4时,坝壳与心墙的临界弹模比Re,cr为7,当心墙泊松比νc取0.45时,坝壳与心墙的临界弹模比Re,cr大于8.由上面的结论,可绘出心墙泊松比与坝壳与心墙的临界弹模比的对应关系表格,如表1所示由表1,根据拉格朗日插值多项式,可过表中四个点,拟合出一个三次多项式,如式(1),所示假定心墙泊松比为横坐标,坝壳与心墙的临界弹模比为纵坐标,由式(1),可绘出心墙泊松比与坝壳与心墙的临界弹模比关系图,如图5所示接下来,在两座大坝心墙泊松比νc都取0.45,其他参数的选取与之前计算参数相同的情况下再分别计算,得到坝壳与心墙的弹模比eR为10,11,12时心墙上游面竖向应力σy随深度的变化曲线,如图6所示虽然图5中的拟合曲线,在描述心墙泊松比取值为0～0.2时,结果较保守,但由于设计时,心墙泊松比一般取值为0.25～0.45之间,通过两个不同模型计算结果对比,得出一致的结果,且与拟合曲线结果的相符,所以该拟合曲线在心墙泊松比取值为0.25～0.45之间时,能较好的表达心墙泊松比与坝壳和心墙的临界弹模比的关系.3坝壳与心墙的临界弹模比与临界弹模比1)由于两个坝体的坝壳弹性模量不同,在心墙泊松比一定时,通过计算得出相同的坝壳与心墙的临界弹模比,说明了影响水力劈裂发生的是坝壳与心墙弹模比的大小,而不是坝壳与心墙弹模差的大小.2)通过对两个坝体计算结果的整理,求出了心墙泊松比与坝壳心墙的临界弹模比关系表达式:同时,在这里值得注意的时,本文在求心墙泊松比和坝壳与心墙的临界弹模比关系表达式时,是假定其它参数不变的,实际上,坝壳的密度和心墙上游的倾角对水力劈裂的影响也是较大的,所以在设计时,要综合考虑所有因素,将水力劈裂的可能性降到最低.由图3可以看出,当心墙泊松比νc取0.25时,坝壳与心墙的临界弹模比Re,cr为3,当心墙泊松比νc取0.3时,坝壳与心墙的临界弹模比Re,cr为4,当心墙泊松比νc取0.35时,坝壳与心墙的临界弹模比Re,cr为5,当心墙泊松比νc取0.4时,坝壳与心墙的临界弹模比Re,cr为7,当心墙泊松比νc取0.45时,即使坝壳与心墙的临界弹模比eR为8,心墙竖向应力仍远大于心墙上游面水压力,也就意味着不会发生水力劈裂.由图4可以看出,当心墙泊松比νc取0.25时,坝壳与心墙的临界弹模比Re,cr为3,当心墙泊松比νc取0.3时,坝壳与心墙的临界弹模比Re,cr为4,当心墙泊松比νc取0.35时,坝壳与心墙的临界弹模比Re,cr为5,当心墙泊松比νc取0.4时,坝壳与心墙的临界弹模比Re,cr为7,当心墙泊松比νc取0.45时,即使坝壳与心墙的临界弹模比eR为8,心墙竖向应力仍