Abstract：In this paper，firstly the momentum equilibrium equation，the continuity equation for water，the energy balance equation and the elastoplastic matrix are established；secondly by using the Galerkin method and isoparametric elements and 4-node elements with equal thickness for porous media and discontinuities respectively，all the governing equations are discretized in the space domain and in the time domain respectively；then a two dimensional FEM code for analysis of coupled thermo-hydro-mechanical elastoplastic processes in saturated-unsaturated porous media with discontinuities is developed tentatively. Through the numerical computation for an assumed coupled thermo-hydro-mechanical phenomena in nuclear waste repository considering the case without discontinuity and the case with discontinuities，it is shown：because of the existence of discontinuities with very large permeability，the time at which the buffer material becomes saturation is put off；and the displacement，stress and fluid flow fields are very different compared with the case without discontinuity.
