Запропоновано неявні двоточкові диференціально-різницеві паралельні ітераційні алгоритми декомпозиції області для розв'язування задач про контакт багатьох пружних тіл. На основі скінченноелементних апроксимацій здійснено іх програмну реалізацію для випадку плоских контактних задач. Досліджено вплив ітераційних параметрів на швидкість збіжності розроблених алгоритмів. Проведено порівняння числової ефективності двоточкових та одноточкових ітераційних методів декомпозиції області.