В работе предложен сопряженный алгоритм, представляющий собой комбинацию метода конечных элементов (МКЭ) и метода блочных элементов (МБЭ), для моделирования механического поведения массивов горных пород в окрестности глубоких подземных сооружений. МКЭ используется для расчета напряженно-деформированного состояния (НДС) вмещающего массива горных пород и выявления зон предельного состояния (ПС) в окрестности подземного сооружения посредством использования комплексного критерия предельного состояния. В таких зонах в дальнейшем применяется МБЭ для моделирования механического состояния области массива дискретной структуры с использованием деформируемых блочных элементов, что позволяет точно описывать локализованные зоны нарушения сплошности в массиве (разрушения, сдвигов и вывалов породных масс). Эффективность алгоритма подтверждена численным решением двух классических задач геомеханики: задачи об устойчивости одиночной выработки и задачи моделирования обрушения породных масс при ведении горных работ лавами применительно к месторождениям калийных солей Беларуси. Верификация результатов моделирования данными натурных замеров показала погрешность 5–17 % в количественных показателях при адекватном качественном повторении исследуемых геомеханических процессов. Среди преимуществ разработанного алгоритма учет неоднородности породных массивов (за счет использования методов механики дискретных сред), экономия вычислительных и временных ресурсов при проведении численных расчетов.