Discontinuum-based numerical methods can contain the multiple discontinuities in a model and reflect the thermal, hydraulic and mechanical characteristics of discontinuities. Therefore, discontinuum methods can be appropriate to simulate the model which require the detailed analysis of the coupled thermo-hydro-mechanical processes in fractured rock such as geothermal energy, CO2 geo-sequestration, and geological repository of the high-level radioactive waste. TOUGH-3DEC, the three-dimensional discontinuum simulators for the coupled thermo-hydro-mechanical analysis, was developed by linking the integral finite difference method TOUGH2 and the explicit distinct element method 3DEC to describe the coupled thermo-hydro-mechanical processes in both porous media and discontinuity. TOUGH2 handles thermo-hydraulic analysis by the internal simulation module, and 3DEC performs mechanical study based on the constitutive models of porous media and discontinuity with coupling the thermal and hydraulic response from TOUGH2. The thermal and hydraulic couplings are the key processes and should be carefully verified by sufficient cases, so this study performed the thermomechanical and hydro-mechanical simulations which are modelling the analytic solutions including the uniaxial consolidation, fracture static opening, and the heating of a hollow cylinder problems. Each thermo-mechanical and hydro-mechanical verification case is also validated by comparing with the results of the other continuum and discontinuum-based numerical methods. TOUGH-3DEC results follow the analytic solutions and show better accuracy than the continuum-based numerical methods in the static fracture opening problem. The developed TOUGH-3DEC simulator can be expanded to coupled thermo-hydro-mechanical-chemical analysis in fractured rock mass, and the simulator needs to be verified by more complicated coupled processes problems which require in the chemical coupling.