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

水平面渗流有限元计算中减压井点 处理 摘要 为解决水平面二维有限元方 法计算减压井时井点附近奇异区问题， 在李 祖贻等针对特定网格划分推导出修正井水 位法的成果基础上，推导适应一般网格划分 的修正公式。通过特例分析计算，验证修正 公式精度和适应性，减压井计算受单元尺寸、 尺寸差异影响很小，流量和井点以外节点水 头值满足精度要求。针对有限元计算时涉及 的井阻力情况，给出实用的修正公式处理方 法。 关键词有限元 减压井渗流中断面法 20减压井是堤防防渗加固的一种常用工 程措施，需要定量分析它的出水量和减压效 果。多数堤防强透水层水平成层分布，满足 缓变渗流条件[1],可以用水平面二维有限 元方法。井的公式大多在缓变渗流条件下推 导出的，因此在水平面有限元计算中，可将 减压井设为一个节点，以出水量和井内外水 位差的协调，将井的公式与其结合，解决计 算问题。因井点附近是奇异区，水头分布为 对数旋转面，无法用有限个平面或低阶曲面 拟合，必须作特殊处理。 李祖贻等［1,2］以在井周划分为4个相同的 等腰直角三角形单元的特殊情况 （图1），推 出修正井水位法及修正井周单元渗透系数 法。设hw为井水位，h0为井节点计算水头， ha为井周节点计算水头。井水位修正的要求 是：给定修正量 △ h,当h0满足h0=hw+A h 时，计算得井的出水量Q和井周节点水头ha 与解析解相同。由此得 式中：q=Q