基于光滑有限单元法的土石坝无压渗流场数值模拟

获得自由面准确位置是无压渗流分析中最被关注的问题.采用光滑有限元法对土石坝无压渗流场进行数值模拟,对不同模型的自由面位置进行求解,研究了渗透异常体对自由面和渗流参数的影响.根据梯度光滑技术,建立了无压渗流问题的光滑有限元模型,该方法的优势在于将单元面积分优化为沿单元边界的线积分,简化了边界相交单元的内部积分过程,降低了单元畸形对计算结果的影响.通过将矩形均质坝和梯形均质坝的经典算例计算结果和其他数值方法的结果进行对比,说明了本文提出的方法的可行性.进一步研究了不同的渗透率异常体对水头分布、渗流速度和流体压力等渗流参数的具体影响,这些结果为利用渗流参数分布的变化反演渗透异常体位置,以及进行土石坝多场耦合提供了理论依据和可行策略.