Robust multigrid methods for nonsmooth coefficient elliptic linear systems