Taking a wide range of the initial conditions, including spatial density distribution and mass function etc, the dynamical evolution of globular clusters in the Milky Way is investigated in detail by means of Monte Carlo simulations. Four dynamic mechanisms are considered: stellar evaporation, stellar evolution, tidal shocks due to both the disk and bulge, and dynamical friction. It is found that stellar evaporation dominates the evolution of low-mass clusters and all four are important for massive ones. For both the power-law and lognormal initial clusters mass functions, we can find the best-fit models which can match the present-day observations with their main features of the mass function almost unchanged after evolution of several Gyr. This implies that it is not possible to determine the initial mass function only based on the observed ones today. Moreover, the dispersion of the modelled mass functions mainly depends on the potential wells of host galaxies with the almost constant peaks, which is consistent with current observations. (C) 2010 COSPAR. Published by Elsevier Ltd. All rights reserved.