根据修正的SPV(Stavropoulou-Papanastasiou-Vardoulakis)管涌控制方程,采用Newton-Raphson求解法和传导矩阵调整法,对二维带自由面的均质土石坝坝体渗流侵蚀问题进行了有限元数值分析。数值分析表明:在具有渗流自由面的均质土石坝坝体中,溢出面和自由面附近区域的水力梯度较大,比其他区域更容易大于管涌临界水力梯度。管涌侵蚀区(孔隙率明显增大)和细颗粒最大浓度区首先出现在溢出面处,然后沿着渗流自由面附近区域向上游扩展,最终达到自由面附近的进水面,形成孔隙率明显增大的侵蚀贯通通道。管涌侵蚀增大堤坝的渗透系数,它对堤坝渗透速度的影响较大,但对孔压分布和自由面位置的影响较小。