In this paper, we show the use of parallel computing in the finite element analysis of thermo-hydro-mechanical problems. Our approach is based on iterative solvers, partitioning of vectors and use of Schwarz type overlapping domain decomposition preconditioners. This technique is applied to standard, time-stepping and mixed finite element analysis, respectively.