This study presents calculations results of thermal consolidation process of the porous medium with the rheological Kelvin-Voigt skeleton, obtained numerically with the use of Flex.PDE software. The investigated calculation scheme consisted of the porous column filled with a liquid. The vertical load was applied to the top surface of the column through a porous plate allowing the free flow of liquid through this surface. Numerical solution is based on compression of the sample at appropriately defined boundary conditions. The aim of this study was to describe the influence of external load and temperature gradient on the deformation tests progress at different values of three parameters: λ,