无压渗流分析的复合单元法

针对无压渗流问题，在传导矩阵调整法的基础上，对穿越自由面的单元用复合单元来处理，每个复合单元具有两套结点水头，一套用于饱和区的水头插值；另一套则用于无水区。通过变分原理，推出了控制方程并在程序中实现。应用复合单元法可以较好地解决自由面的插值问题，消除了穿过自由面单元的水上部分对自由面附近结点水头的贡献，同时对逸出面边界条件的处理也更加方便和准确。算例表明了该方法的有效性。