Рассматривается численное моделирование напряженно-деформированного состояния участков литосферы в некоторых областях месторождений полезных ископаемых. В качестве объектов исследования выбраны двумерные структурные модели Вилюйской и Тунгусской синеклиз на основе профиля «Кимберлит-1981». Для описания процессов деформирования применена модель упругопластической среды с неассоциированным законом течения в форме Друкера – Прагера – Николаевского, которая была адаптирована с использованием модели «jelly sandwich». Локализация пластической деформации, области положительных горизонтальных напряжений соответствуют местонахождениям полезных ископаемых.