Abstract：Bathe algorithm converts the locating of the phreatic surface to a nonlinear constitutive problem. And the implication of Signorini condition can simulate a seepage face as a head-fixed boundary through iterative calculation. Because an unsaturated seepage problem is also a nonlinear flow problem，the implication of Bathe algorithm and Signorini condition makes it possible to model both saturated and unsaturated seepage problems with a unified method by a minimal modification to an ordinary finite element method，and avoiding solving a variational inequality system. The followings are discussed mainly：(1) the improvements of Bathe algorithm in converge and its generalized form in a three-dimensional model；(2) the implement of the switching algorithm of Signorini condition to solving seepage problems；and (3) the under-relaxation scheme to improve the mass conservation and converge properties when an unsteady unsaturated problem is solved. Finally，some numerical examples are solved to evaluate the applicability of the proposed method；and the results are compared with those available in the literature.
[1] LAM L，FREDLUND D G，BARBOUR S L. Transient seepage model for saturated-unsaturated soil systems：a geotechnical engineering approach[J]. Canadian Geotechnical Journal，1987，24(4)：565–580.
[2] HUYAKORN P S，THOMAS S D，THOMPSON B M. Techniques for making finite elements competitive in modeling flow in variably saturated porous media[J]. Water Resources Research，1984，20(8)：1 099–1 115.
[3] BATHE K J，KHOSHGOFTAAR M R. Finite element free surface seepage analysis without mesh iteration[J]. International Journal for Numerical and Analytical Methods in Geomechanics，1979，3(1)：13–22.
[4] DESAI C S，LI G C. A residual flow procedure and application for free surface flow in porous media[J]. Advances in Water Resources，1983，6(1)：27–35.
[5] 张有天，陈 平，王 镭. 有自由面渗流分析的初流量法[J]. 水利学报，1988，19(8)：18–26.(ZHANG Youtian，CHEN Ping，WANG Lei. Initial flow method for seepage analysis with free surface[J]. Journal of Hydraulic Engineering，1988，19(8)：18–26.(in Chinese))
[6] LACY S J，PREVOST J H. Flow through porous media：a procedure for locating the free surface[J]. International Journal for Numerical and Analytical Methods in Geomechanics，1987，11(6)：585–601.
[7] 王 媛. 求解有自由面渗流问题的初流量法的改进[J]. 水利学报，1998，29(3)：69–74.(WANG Yuan. The modified initial flow method for 3D unconfined seepage computation[J]. Journal of Hydraulic Engineering，1998，29(3)：68–74.(in Chinese))
[8] LARABI A，SMEDT F D E. Numerical solution of 3D groundwater flow involving free boundaries by a fixed finite element method[J]. Journal of Hydrology，1997，201(1/4)：161–182.
[9] 郑 宏，戴会超，刘德富. 改进的有自由面渗流问题的Bathe算法. 岩土力学，2005，26(4)：505–512.(ZHENG Hong，DAI Huichao，LIU Defu. Improved Bathe?s algorithm for seepage problems with free surfaces[J]. Rock and Soil Mechanics，2005，26(4)：505–512.(in Chinese))
[10] ZHENG H，LIU D F，LEE C F，et al. A new formulation of Signorini?s type for seepage problems with free surfaces[J]. International Journal for Numerical Methods in Engineering，2005，64(1)：1–16.
[11] AITCHISON J M，POOLE M W. A Numerical Algorithm for the solution of signorini problems[J]. Journal of Computational and Applied Mathematics，1998，94(1)：55–67.
[12] GENUCHTEN M T V. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils[J]. Science Society of America Journal，1980，44(5)：891–898.
[13] FREDLUND D G，ANQING X，SHANGYAN H. Predicting the permeability function for unsaturated soils using the soil-water characteristic curve[J]. Canadian Geotechnical Journal，1994，31(4)：533–546.
[14] GARDNER W R. Some steady-state solutions of the unsaturated moisture flow equation with application to evaporation from a water table[J]. Soil Science，1958，85(4)：228–232.
[15] PHOON K K，TAN T S，CHONG P C. Numerical simulation of richards equation in partially saturated porous media：under-relaxation and mass balance[J]. Geotechnical and Geological Engineering，2007，25(5)：525–541.
[16] KAZEMZADEH-PARSI M J，DANESHMAND F. Unconfined seepage analysis in earth dams using smoothed fixed grid finite element method[J]. International Journal for Numerical and Analytical Methods in Geomechanics，2011，36(6)：780–797.
[17] BARDET J P，TOBITA T. A practical method for solving free-surface seepage problems[J]. Computers and Geotechnics，2002，29(6)：451–475.
[18] BORJA R I，KISHNANI S S. On the solution of elliptic free-boundary problems Via Newton?s method[J]. Computer Methods in Applied Mechanics and Engineering，1991，88(3)：341–361.
[19] LEE K K，LEAP D I. Simulation of a free-surface and seepage face using Boundary-fitted coordinate system method[J]. Journal of Hydrology，1997，196(1/4)：297–309.
[20] STRACK O D L，BARNES R J，VERRUIJT A. Vertically integrated flows，discharge potential，and the dupuit-forchheimer approximation[J]. Ground Water，2006，44(1)：72–75.
[21] SRIVASTAVA R，YEH T C J. Analytical solutions for one- dimensional，transient infiltration toward the water table in homogeneous and layered soils[J]. Water Resources Research，1991，27(5)：753–762.
[22] CELIA M A，BOULOUTAS E T，ZARBA R L. A general mass- conservative numerical solution for the unsaturated flow equation[J]. Water Resources Research，1990，26(7)：1 483–1 496.
[23] RULON J J，RODWAY R，FREEZE R A. The development of multiple seepage faces on layered slopes[J]. Water Resources Research，1985，21(11)：1 625–1 636.
2013, Vol. 32 
