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