<p class="p15">The main innovation of this paper includes two parts. One part is the discrete formulas of Thermo-hydro-mechanical (THM) coupling equations and another part is the discussion of the truncation errors based on the Taylor formula. There are many THM coupling problems in unsaturated soils, which are very important in both theoretical and engineering applications. The numerical computing of coupling equations is increasingly important. Considering the deformation of unsaturated soils skeleton, fluid flow and heat transfer, constitutive relationships of the THM coupled behavior are given. Then, the constitutive equations are derived and a closed problem is formed. The equations are dispersed by difference method and the truncation errors of the discrete formulas are given.</p>