A numerical study on frictional shakedown in large-scale three-dimensional conforming elastic contacts