A Monte Carlo approach has been developed for simulations of the angular and energy distributions for heavy evaporation residues (ER) produced in heavy ion fusion-evaporation reactions. The approach uses statistical model approximations of the HIVAP code for the calculations of initial angular and energy distributions inside a target, which are determined by neutron evaporation from an excited compound nucleus. Further step in the simulation of transmission of ER heavy atoms through a target layer is performed with the TRIM code that gives final angle and energy distributions at the exit from the target. Both the simulations (neutron evaporation and transmission through solid media) have been separately considered and good agreement has been obtained between the results of simulations and available experimental data. Some applications of the approach have been also considered.