Iterative solvers are of increasing interest in geomechanics with the move towards 3D finite element modelling. Potentially, these methods can lead to reduced computational complexity as, unlike direct methods, they do not require the full system matrix to be assembled. In general, however, iterative solvers have not been widely adopted in geomechanics due to problems with convergence. This pap...