We introduce a new multigrid method to study the lattice statics model arising from nanoindentation. A constrained Cauchy-Born elasticity model is used as the coarse-grid operator. This method accelerates the relaxation process and considerably reduces the computational cost. In particular, it saves quite a bit when dislocations nucleate and move, as demonstrated by the simulation results.