When a landslide takes place, it is believed that a shear band of loose granular media acts as a lubricant between the descending block of soil and the basis on repose. The mechanism involved is known as softening: the granular skeleton looses its stiffness and the shear stress on the block is lost. In the hypothesis of Habib, the friction between grains heats the pore water, increasing its pressure and reducing the effective stress by the Terzagi criterion. Vardoulakis had constructed models on this hypothesis including thermal diffusion and Darcy's law, plus a double dependence of the friction angle on the displacement and the velocity of the rolling block. Hereby we present a discrete element simulation of the process on a tilted shear band between two soil blocks: one bottom at rest and one upper at move. Soil blocks are assumed with uniform permeability and thermal conductivity. The shear band is modeled as a set of Voronoi polygons with elastic, frictional and damping forces between them. Pore water acts with hydrostatic pressure on the grains and on the upper and lower blocks, with a thermodynamic response that is reproduced by the Steam Tables provided by the International Association for the Properties of Water and Steam (IAPWS 97 report). At each time step, the forces on all grains are computed and all translational and rotational movements are integrated. Then, the heat is computed as the work done by all dissipative forces, distributing between water and grains according to their thermal capacities and increasing water temperature and pressure. Finally, this water pressure pushes the grains apart, reducing the shear stress on the upper block and speeding up the landslide. By this simulation procedure we obtain temperature increments on 10 C° that are strong enough to produce softening. Although the model is in two dimensions, it provides new insights on the study of catastrophic landslides evolutions.