Thanks to the Gaia mission, it will be possible to determine the masses of approximately hundreds of large main belt asteroids with very good precision. We currently have diameter estimates for all of them that can be used to compute their volume and hence their density. However, some of those diameters are still based on simple thermal models, which can occasionally lead to volume uncertainties as high as 20-30%. The aim of this paper is to determine the 3D shape models and compute the volumes for 13 main belt asteroids that were selected from those targets for which Gaia will provide the mass with an accuracy of better than 10%. We used the genetic Shaping Asteroids with Genetic Evolution (SAGE) algorithm to fit disk-integrated, dense photometric lightcurves and obtain detailed asteroid shape models. These models were scaled by fitting them to available stellar occultation and/or thermal infrared observations. We determine the spin and shape models for 13 main belt asteroids using the SAGE algorithm. Occultation fitting enables us to confirm main shape features and the spin state, while thermophysical modeling leads to more precise diameters as well as estimates of thermal inertia values. We calculated the volume of our sample of main-belt asteroids for which the Gaia satellite will provide precise mass determinations. From our volumes, it will then be possible to more accurately compute the bulk density, which is a fundamental physical property needed to understand the formation and evolution processes of small solar system bodies.