In the THM modeling of multiphase medium, the coupling effects of skeleton, suction, and temperature have been integrated via the concept of state surfaces of void ratio and degree of saturation. Based on proposed formulation, a fully coupled numerical model for the behavior of soil deformation, water flow, air flow, heat flow in unsaturated soil has been developed and integrated in a finite element code θ-Stock by the first author. This program is conceived with this idea that it will be able to analyze the response of a soil in different states of humidity to mechanical, thermal loading, and also damage phenomena. Damage model is dedicated to unsaturated brittle rocks. It mixes phenomenological and micromechanical concepts and is formulated based on the use of independent state variables. The expression of the liquid permeability is modified in order to represent the influence of fracturing on interstitial fluid flows. The final matrix form of established field equations of the proposed model for unsaturated case has been encoded for this particular purpose, in a finite element program which had been developed for dry and saturated soils previously.