We present results from extensive Monte Carlo simulations of polymer models where each lattice site can be visited by up to K monomers and no restriction is imposed on the number of bonds on each lattice edge. These multiple monomer per site (MMS) models are investigated on the square and cubic lattices, for K = 2 and K = 3, by associating Boltzmann weights ω0 = 1, ω1 = e β 1 and ω2 = e β 2 to sites visited by 1, 2 and 3 monomers, respectively. Two versions of the MMS models are considered for which immediate reversals of the walks are allowed (RA) or forbidden (RF). In contrast to previous simulations of these models, we find the same thermodynamic behavior for both RA and RF versions. In three-dimensions, the phase diagrams -in space β2 × β1 -are featured by coil and globule phases separated by a line of Θ points, as thoroughly demonstrated by the metric νt, crossover φt and entropic γt exponents. The existence of the Θ-lines is also confirmed by the second virial coefficient. This shows that no discontinuous collapse transition exists in these models, in contrast to previous claims based on a weak bimodality observed in some distributions, which indeed exists in a narrow region very close to the Θ-line when β1 < 0. Interestingly, in two-dimensions, only a crossover is found between the coil and globule phases.