水平面渗流有限元计算中减压井点处理论文.docVIP

　　水平面渗流有限元计算中减压井点处理论文 摘要：为解决水平面二维有限元方法计算减压井时井点附近奇异区问题，在李祖贻等针对特定网格划分推导出修正井水位法的成果基础上，推导适应一般网格划分的修正公式。通过特例分析计算，验证修正公式精度和适应性，减压井计算受单元尺寸、尺寸差异影响很小，流量和井点以外节点水头值满足精度要求。针对有限元计算时涉及的井阻力情况，给出实用的修正公式处理方法。 关键词：有限元 减压井 渗流 中断面法 20 (1华南理工大学 土木工程系，广东 广州510641；2广东省水利水电科学研究院.freel点水头ha，单元流量qe按中断面法计算1： (2) 式中：Δ为单元面积；bi,bj,bm,ci,cj,cm为单元节点对边向y及x轴投影长度，可在一般有限元书中找到。注意到图2中流量定义与井出水相反，按井点习惯出水为正，反号后由式(2)推得： (3) 设流向井点的水均匀分布，总流量q与qe的关系为q=nqe=2π=qe/φ，整理可得 (4) 按修正的要求，同一流量下，半径为r圈节点距井点ma，节点数mn，两节点对应圆心角φ/m ，坐标旋转为图3所示情况下，φ角对应区域的m+1个节点坐标为(macosi/φm,masini/φm(i=0,.freel)，节点水头hm。该圈与m-1圈间共(2m-1)n个单元，其中mn个单元有两个节点在m圈上，(m-1)n个单元只有一个节