We present N-body simulations of globular clusters including gravitational field of the Galaxy, in order to study effects of tidal field systematically on the shape of outer parts of globular clusters using NBODY6. The Galaxy is assumed to be composed of central bulge and outer halo. We investigate the cluster of multi-mass models with a power-law initial mass function (IMF) starting with different initial masses, initial number of particles, different slopes of the IMF and different orbits of the cluster. We have examined the general evolution of the clusters, the shape of outer parts of the clusters, density profiles and the direction of tidal tails. The density profiles appear to become somewhat shallower just outside the tidal boundary consistent with some observed data. The position angle of the tidal tail depends on the location in the Galaxy as well as the direction of the motion of clusters. We found that the clusters become more elongated at the apogalacticon than at the perigalacticon. The tidal tails may be used to trace the orbital paths of globular clusters.