the energy balance equation and the elastoplastic matrix are established. Secondly by using the Galerkin method all the governing equations are discretized in the space domain and time domain respectively
then a two dimensional FEM code for coupled thermo-hydro-mechanical elastoplastic analysis of saturated porous media is developed tentatively. Through two computational examples the use try of the code is made.