岳城水库土坝稳定渗流有限元剖析.docVIP

岳城水库土坝稳定渗流有限元剖析 　　【摘 要】土坝渗流问题是引起土坝破坏的主要原因之一。准确计算土坝渗流，确定渗流自由面和渗流等势线是非常必要的。本文利用有限元的方法，对均质土坝的稳定渗流进行了分析，并运用于岳城水库大坝，为今后该大坝的安全监控及防渗除险加固提供了依据。 　　【关键词】土坝；渗流；有限元 　　Yue City Reservoir Dam seepage finite element analysis 　　Du Er-xia,Liu Shao-xuan,Zhao Shao-ping,Wang Ya-ru 　　(Hebei University of Architecture and Civil Engineering Baoding Hebei 071002) 　　【Abstract】The seepage through earth dams is one of the main reason of dam destruction. It is necessary to calculate the free surface and the amounts of seepage . The steady seepage through earth dams is analyzed with the finite element method and the seepage through Yuecheng earth dam is analyzed in this paper. Based on the result of the calculation, many reasonable suggestions are posed to help the management and reinforcement of Yuecheng reservoir earth dam. 　　【Key words】Earth dam;Steady seepage;Finite element method 　　我国是世界上兴修水利最早，拥有水利设施数量最多的国家。土坝岸堤的渗流破坏一种形式为管涌另一种为整个土体的滑坡，从计算分析和统计数据可知堤坝岸坡的稳定性渗流的破坏作用很大，垮坝失事有30%到40%源于渗流。本研究通过对土坝的渗流分析模型的比较，选择适合的渗流的计算模型，并针对该模型运用有限元的方法编制分析二维渗流的计算程序，迭代计算出稳定渗流自由面，并进一步研究其对堤坝岸坡的稳定的影响。 　　1. 计算原理 　　考虑土坝中土和水的压缩性，二维非均质各项异性土体渗流，根据Darcy定理及水流连续性条件，当渗流系数主方向与坐标方向一致时，非稳定渗流可以用如下微分方程表示：?功?x(k??x??H??x)+?功?z(k??z??H??z)=0（1） 　　式中 k??x 、k??z 是 以x、z轴为主轴方向的渗透系数； H是水头函数。给出 　　公式（1）的边界条件通常有三类，即： 　　水头边界：H=H??0（2） 　　流量边界条件：k??x??H??x+k??z??H??zl??z+q=0（3） 　　不透水边界： 　　??H??n=0 （4） 　　式中： H??0为已知水头，l??x、 l??z 是边界表面的法向余弦； q为单位边界表面的流入或流出量，n为边界的外法线方向。对于均质土的单一流场来说，计算时不管渗透系数为何值，都能求得相同的流网；但对于不同土层分布的渗流场来说，就得考虑各土层的渗透系数。 　　2. 有限元程序的编制 　　本文用有限元法求解二维稳定渗流的问题，选用三角形单元为基本单元形式，3个角点i 、j 、 m的次序应按逆时针方向循环。水头插值函数采用线性插值函数。特别指出，单元划分时应使研究的水头值在单元内部接近一个常数，在急变区域应划分小单元，变化缓慢区可划分大单元；而且单元边长不宜过于悬殊，邻边上也不能出现间断的结点。 　　根据变分原理，可用有限元法讨论稳定渗流的微分方程所对应的泛函： 　　 I(h) =∫Q［12k??x(??h??x)??2+12k??z(??h??z)??2］dxdz（5） 　　由上述原理可建立单元方程并汇总为总的方程组，形成一系列计算公式而编制程序求解。将所有单元的泛函微分式汇总起来，取极值，就得到整个流场的泛函对结点值求导数的方程组。将此方程组写成总的矩阵形式为： 　　［K］｛h｝=｛f｝（6） 　　3. 计算实例 　　河北省岳城水库土坝为粘壤土填筑的大致均质土坝，坝高50m，坝顶宽10m，坝底宽600m，自由面变化的范围内土质的给水度μ=0.047 。水库上游水位为130m ，下游水位 105m。其渗透系数分别为第三纪地层为8.5