李炳阳，杨旭亮

(1.中国电建集团北京勘测设计研究院有限公司，北京 100024；2.青海民族大学土木与交通工程学院，青海 西宁 810007)

渗流区域内的渗透压力、渗透坡降、渗透流速及渗流量的求解，是土质地基闸基渗流计算的基本内容，对水闸防渗排水设计起着重要支撑作用。在实际工程中，闸基渗流区域的边界条件通常较为复杂，解析求解闸基渗流运动的拉普拉斯方程难度很大，因而常采用一些近似且实用的方法，如改进阻力系数法和流网法[1]。改进阻力系数法系根据地下轮廓特点将整个渗流区域分成几个典型流段，求解各个典型流段的阻力系数，进而算出任一流段的水头损失、渗流压力及其他渗流要素，由于不需要绘制流网，且计算精度满足工程应用要求，目前已经被设计人员广泛采用。但当地下轮廓线复杂或闸下为多层地基土时，改进阻力系数法存在计算过程繁琐、费时和适用性问题。规范[2]指出，复杂土质地基上的重要水闸，应采用数值计算法进行计算。在水闸渗流场数值分析中建立饱和—非饱和渗流模型，利用有限元软件分析两种工况下水闸底板的扬压力和水力坡降，研究闸底区域渗流流场的分布规律[3]；在水闸闸基渗流分析中采用改进阻力系数法和Autobank软件进行计算，对比分析了数值方法计算结果的可靠性[4-5]。利用slide软件对渗流作用下的土质边坡、尾矿坝等进行稳定性计算[6- 8]，表明该有限元软件可以较好地模拟地下水渗流作用。……