Изложены алгоритмы прогноза напряженно-деформированного состояния неоднородного слоистого физически нелинейного породного массива. Моделирование базируется на применении комплекса вычислительных методов: вариационного, дискретного продолжения по числовому параметру, квазилинеаризации нелинейных краевых задач, конечных разностей, матричной прогонки и общего итерационного процесса.
С развитием вычислительной техники в математическом моделировании наряду с традиционными аналитическими методами все шире применяются численные методы. Среди численных методов решения наиболее широко известны метод конечных разностей, метод конечных элементов, метод граничных элементов, метод граничных интегральных уравнений и другие методы, получившие большое распространение в инженерной и научной среде. К числу наиболее эффективных численных методов относится МКР, широко используемый при решении важных прикладных задач геомеханики.