Строится математическая модель массива с учетом его внутренней структуры, анизотропии, разупрочнения, аккумулирования и высвобождения упругой энергии. Эффективные упругие модули зернистой микронеоднородной среды определяются методом граничных интегральных уравнений с использованием тензора фундаментальных решений третьего рода. Диаграмма межзеренного проскальзывания учитывает стадию разупрочнения. Получено численное решение задачи о деформировании целика горного массива