This paper introduces the model of pressure solution for granular aggregates established by Taron et al. into the finite element code for analysis of thermo-hydro-mechanical (T-H-M) coupling in porous media developed by the first author. Using the Mohr-Coulomb yield criterion, a hypothetical disposal model for nuclear waste located in a saturated quartz aggregate rock mass with a laboratory scale is simulated. Two computation cases, the elastic analysis and the elastoplastic analysis are designed. Then the corresponding numerical simulation for a disposal period of 4 years is carried out. The states of temperatures, solute concentrations in the intergranular fluid film and at the pore space, migration and precipitation masses, porosities and permeabilities, pore pressures, flow velocities and stresses and plastic zones in the rock mass are investigated. It is shown that, because of the stress adjustment and the increased molecular diffusivity in the elastoplastic analysis, there are obvious changes of solution, migration and precipitation of aggregate medium in the plastic zones, and the seepage field (including pressures and flow velocities of pore water) and the stress field are markedly influenced. But the differences of solution, migration and precipitation of aggregate medium in the elastic zones for two cases are not significant.