We perform high resolution N-body simulations for f(R) gravity based on a self-adaptive particlemesh code MLAPM. The Chameleon mechanism that recovers General Relativity on small scales is fully taken into account by self-consistently solving the non-linear equation for the scalar field. We independently confirm the previous simulation results, including the matter power spectrum, halo mass function and density profiles, obtained by Oyaizu et al. (Phys.Rev.D 78, 123524, 2008) and Schmidt et al. (Phys.Rev.D 79, 083518, 2009), and extend the resolution up to k � 20 h/Mpc for the measurement of the matter power spectrum. Based on our simulation results, we discuss how the Chameleon mechanism affects the clustering of dark matter and halos on full non-linear scales.