Study on performance of Krylov subspace methods at solving large-scale contact problems