В роботi запропоновано наближений метод розв'язання осесиметричних задач пружних тiл з адгезiєю. Цей метод дозволяє обчислювати розмiр контактної плями, величину сили взаємодiї мiж тiлами а також розподiл контактного тиску без iстотного обмеження форми початкового зазору мiж тiлами. Отже вiн є бiльш унiверсальним, нiж iснуючi аналiтичнi теорiї, iнструментом для дослiдження таких явищ, як адгезiйна мiцнiсть спряжених тiл, а також пов'язана з втратою стiйкостi дисипацiя енергiї пiд час циклiчного контактного навантаження. Був запропонований варiацiйний принцип, з якого можливо побудувати наближений розв'язок. Побудована дискретизацiя та отриманi нелiнiйнi рiвняння для визначення невiдомого радiусу кругової плями контакту та вузлових значень контактного тиску. На вiдмiну вiд iснуючих чисельних методiв, в яких область контакту в ходi уточнення змiнюється за рахунок додавання або вилучення окремих граничних елементiв цiлком, запропонований пiдхiд дозволяє обчислювати поступову, неперервну змiну площi контакту. Для отримання повної кривої адгезiйної взаємодiї включно з дiлянками, де система втрачає стiйкiсть, запропонована реалiзацiя методу подовження. При цьому задля досягнення збiжностi контроль за змiнними навантаження здiйснювався з урахуванням змiни площi плями контакту на додачу до приросту значень зближення та сили. Здiйснено оцiнку точностi наближеного методу в порiвняннi з вiдомими аналiтичними розв'язками, яка показала лiнiйну збiжнiсть iз зменшенням розмiру сiтки за обчислюваною силою та площею плями контакту. Запропонований пiдхiд планується розповсюдити на випадок довiльної форми тiл, що контактують, що дозволить дослiдити вплив випадкової шорсткостi поверхонь на їхнi адгезiйнi властивостi Ключовi слова: адгезiйний контакт, метод граничних елементiв, варiацiйний принцип Калькера, хвиляста шорсткiсть, метод неперервного подовження