不同防渗形式下大坝基坑渗流场数值模拟.docVIP

不同防渗形式下大坝基坑渗流场数值模拟摘要：文章通过某U形河道基坑渗流场的二维有限元模拟，对比分析了在不同的防渗形式下基坑渗流场实际复杂情况，比较理想地反映出基坑渗流场中等势线、浸润线以及流量等渗流要素在不同条件下的变化规律，为采取经济合理的基坑防渗排水措施提供了一定的理论依据。算均质土石围堰的渗流量时，基岩渗流起主导作用，可以只考虑基岩的流量，围堰的渗流量可以忽略不计。垂直心墙防渗与斜墙与粘土铺盖组合式防渗是很好的选择。 关键词：基坑，渗流场，有限元，防渗，土石围堰 0引言 在水利工程中，基础工程都处于河道周围或者河道中，地下水位较高，在这些工程中进行基坑开挖时若不采取有效的排水措施，容易导致流砂，管涌，严重时会导致基坑塌方、邻近地层掏空下陷等事故。因此进行渗流场分析是非常必要的，对渗流场进行合理的数值模拟有着十分重要的意义。60年代以来，数值方法中的有限单元法随着计算机的迅速发展而逐渐发展起来。由于它能够考虑复杂的边界条件及不同介质，应用于基坑工程，能够较好地模拟开挖工况及土与结构共同作用的问题，而且能够同时计算土体与支护结构的变形和应力，理论上更加合理，是一种很有前景的基坑工程分析方法。 文章根据有限单元法来分析某U形河道基坑渗流场的二维有限元模拟，对比分析了在不同的防渗形式下基坑渗流场实际复杂情况，比较理想地反映基坑渗流场中等势线、浸润线以及流量等渗流要素在不同条件下的变化规律，为采取经济合理的基坑防渗排水措施提供一定的理论依据。 1有限元法计算分析过程 有限单元法以剖分离散和分块插值为指导思想。其基本方法是将连续的求解区域离散化为一组有限个且按一定方式相互联接在一起的单元组合体利用每一个单元内假设的近似函数来分片地表达全求解域上待求的未知场函数。 当不考虑土和水的压缩性，符合达西定律的二向非均质各向同性土体稳定渗流基本方程为： (1) 式中是水头函数：是空间坐标。 通过连续区域的离散化，插值函数的选择，可推导出特性矩阵方程，然后组合各单元的渗透矩阵集合，可形成整个渗流区的总渗透矩阵方程，基坑渗流有限元方程可表示为： [K]{P}= {F}(2) 式中[K]为总渗透矩阵；{P}为未知水头列向量；{F}为自由项列向量。 2算例分析 2.1有限元计算模型 假定在宽100m的U形河道上拟建一土石坝，设计为中坝。现需在上游和下游修建围堰，以便进行基坑开挖施工。考虑到充分利用当地废弃的土石方，施工方便等问题，拟采用土石围堰。认为地基和围堰为各向同性。 围堰的上游设计施工挡水位为7m，下游设计施工挡水位为6.5m，土石围堰上下游坡比挡水一侧设为1：1.8，渗流出口侧的坡比为1：1.6，上游围堰高为9m，顶宽7m，下游围堰高为8m，顶宽7m。取地基深度为10m，上下游延长各40m，基坑长50m。当设防渗墙时，厚度取0.8m。计算模型如图1所示 图1基坑模型 2.2拟采用的计算方案 分4种不同工况进行计算，计算方案如表1所示，各种防渗材料的渗透系数如表2所示。 2.3计算方法和过程 首先根据有限单元法的计算过程，对模型进行网格剖分，由于高出水深的围堰不渗流，所以高出水位部位不需要剖分，剖分采用三角形单元。然后用Fortran语言编制了有限元程序对不同防渗形式下基坑各渗流要素进行计算分析。 2.4计算成果整理和分析 (1)水同时从上游和下游向基坑内渗，如图2至图5所示，分别给出了4种不同方案 下，基坑模型的水头分布图，由图2可以看出，在均质土石围堰中，随着地基渗透系增大，等水头分布线在围堰与地基交界处转折越来越大，可以看出如果不同渗透系数的两个区域存在交界面，他们的渗透系数相差越大，水头变化越强烈，对建筑物的稳定越不利，所以当两个不同区域的渗透系数相差大时，要采取相应的防护措施。 （2）结合图2和图6可以看出随着中分化基岩渗透系数得增加，通过围堰的浸润线逐渐变陡，渗流溢出点越来越低，再结合表3与图8，可以分析出，随着地基渗透系数与均质土石围堰渗透系数比值的增大，渗流逐渐主要转向基岩，所以在计算均质土石围堰的渗流量时，基岩渗流起主导作用，可以只考虑基岩的流量，围堰的渗流量可以忽略不计。（图2-6中，横坐标表示距离，纵坐标表示高度，单位：m。坐标系中基岩高十米，围堰上游水头表示的是自由水面与基岩面的高差，浸润线数值表示自由水面在坐标系中的数值。图7中，横坐标表示总水头，纵坐标表示距离，单位：m。图8中，横坐标表示渗透系数值，纵坐标表示流量变化。单位：。） （3）由图3～5，再结合图7，可以看出，方案2和方案4都可以明显的降低渗流的浸润线，基本截断渗流，防渗效果比较好；可以