A finite element algorithm for consideration of contact constraints is presented. It is characterized by introducing the geometric constraints, resulting from contact conditions, directly into the algebraic system of equations for the incremental displacements of an incremental iterative solution procedure. The usefulness of the proposed algorithm for efficient solutions of contact problems involving large displacements and large strains is demonstrated in the numerical investigation.