This project’s aim is to analyze theoretical properties of phase space iterative solvers in the context of linear elasticity. A theoretical proof about the existence of a solution for the algorithm under consideration is given. A quantitative estimate of the convergence rate of the method is also provided, starting from the case of linear constitutive laws and moving to more complex cases with non-linearities.