从建立应力平衡方程、水连续性方程、能量守恒方程入手,利用Galerkin方法,对于岩体和不连续面分别采用等参数单元和具有相同厚度的四节点节理单元,将各控制方程分别在空间域和时间域进行离散,建立了用于分析存在不连续面的饱和–非饱和介质中的热–水–应力耦合弹塑性问题的有限元模型,并开发了相应的计算程序。通过假定的有、无节理的核废料地下处置算例,表现了由于存在渗透系数很大的不连续面,推迟了缓冲材料达到饱和的时间,并对岩体中的位移场、应力场和渗流场产生很大的影响。