该文从饱和多孔介质中的达西定律和连续性原理出发,系统研究了有潜水面的三维非稳定渗流问题的数值计算方法,在已有的稳定渗流计算方法-截止负压法的基础上,提出了一个精确确定非稳定渗流问题的自由面的有限元数值计算模型.该模型采用一般土工问题中通用的方程和网格,在计算迭代过程当中不必调整网格,使得它能用同一套网格同时联合求解渗流场和应力场问题.第一部分建立数值计算模型.根据达西定律和连续性方程,阐述了一般的复杂非稳定渗透的定解问题,推导出了以压力场为未知函数的弱表达式.第二部分采用罚函数有限元和等参数单元,用Galerkin逼近方法和隐式差分法,导出了有限元方程,建立了固定网格迭代的数值计算模型,编制了通用了的程序,并用大量的算例验证该方法的精确性和通用性,再以一个实际工程的具体运用进一步验证了该方法的优点.结合算例,分析研究了该模型中的若干参数对计算结果的影响.第三部分是有限元数值计算的前处理.根据工程常用的由排水廊道和排水孔组成的排渗体的一般型式,开发了该类排渗体的网格数据自动生成软件,使得既可以优选更多的防渗方案又节约大量时间.