有自由面非稳定渗流分析的改进截止负压法

截止负压法是求解有自由面渗流问题较为有效的算法之一。针对其算法特点,结合渗流本质,指出该算法的不足之处在于只引用了单参数的罚函数,从而造成自由面附近区域的渗流计算不够精确。为此,提出了引用双参数罚函数的改进截止负压法,改进后的算法更能准确地考虑单元的部分饱和非饱和作用,精确地计算单元结点外力或者结点不平衡力,进而提高计算精度和收敛速度。最后,将改进算法用于求解有自由面非稳定渗流问题,编制了相应的有限元分析程序,并通过算例进行验证。