A fully coupled numerical model was presented to analyze the process of dissipation of the excess pore water pressure and solid skeleton deformation based on the poroelastic governing equations. The numerical model was verified against a two-dimensional consolidation of homogeneous unsaturated soil. The numerical results show that the non-homogeneity has its own importance not only during the process of dissipation of the excess pore water pressure but also for solid skeleton deformation. As the gradient of soil modulus increases, dissipation rate of excess pore water pressure increases and the amplitude of the soil's deformation decreases in the initial compaction stage and final steady stage.