Modeling submillimeter galaxies in cosmological simulations: contribution to the cosmic star formation density and predictions for future surveys