Моделирование физических процессов в горном деле : компьютерное моделирование

МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РФ 

ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ АВТОНОМНОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ  
ВЫСШЕГО ОБРАЗОВАНИЯ  
«НАЦИОНАЛЬНЫЙ ИССЛЕДОВАТЕЛЬСКИЙ ТЕХНОЛОГИЧЕСКИЙ УНИВЕРСИТЕТ «МИСиС» 

ГОРНЫЙ ИНСТИТУТ 

 

 
 
 

 

 

 

 
 

 

№ 2984 

Кафедра физических процессов горного производства
и геоконтроля 

А.С. Вознесенский 
М.Н. Красилов 
Я.О. Куткин 

Моделирование физических 
процессов в горном деле 

Компьютерное моделирование 

Практикум 

Рекомендовано редакционно-издательским 
советом университета 

Москва  2018

УДК 550.8:004.94 
 
В64 

Р е ц е н з е н т  

д-р техн. наук, проф. М.Н. Шуплик 

Вознесенский А.С. 

В64  
Моделирование физических процессов в горном деле : ком
пьютерное моделирование : практикум / А.С. Вознесенский, 
М.Н. Красилов, Я.О. Куткин. – М. : Изд. Дом НИТУ «МИСиС», 2018. – 97 с. 

ISBN 978-5-906953-08-7 

Практикум содержит методические рекомендации для выполнения прак
тических занятий по ознакомлению и получению навыков работы с системами математических расчетов и компьютерного моделирования MATLAB, 
MathCAD, COMSOL Multiphysics на примерах задач из области механики 
горных пород и массивов. 

Предназначен для студентов специальности 21.05.04 «Горное дело», а 

также в качестве дополнительного материала для студентов специальности 
21.05.05 «Физические процессы горного или нефтегазового производства». 
Может представлять интерес для студентов других специальностей, научных 
и инженерных работников, чья деятельность связана с компьютерным моделированием в области исследования геосред. 

УДК 550.8:004.94 

 
 
 
ISBN 978-5-906953-08-7 

 А.С. Вознесенский, 

М.Н. Красилов, 
Я.О. Куткин, 2018 

 НИТУ «МИСиС», 2018 

СОДЕРЖАНИЕ 

Предисловие ............................................................................................ 6 
1. Оформление отчетов по лабораторным работам и другой 

документации с помощью процессора Word ................................. 8 
1.1. Цель работы ................................................................................ 8 
1.2. Настройка стилей ........................................................................ 8 
1.3. Ввод текста .................................................................................. 11 
1.4. Составление содержания ........................................................... 13 
1.5. Обновление содержания ............................................................ 14 
1.6. Вставка в текст ссылок на литературу ...................................... 14 
1.7. Задание для выполнения ............................................................ 16 
1.8. Содержание отчета ..................................................................... 17 
1.9. Пример оформления титульного листа .................................... 17 
Литература ......................................................................................... 18 
Контрольные вопросы ....................................................................... 18 

2. Расчет в среде Matlab напряжений и деформаций в массиве 

пород вокруг вертикальной горной выработки ............................. 19 
2.1. Цель работы ................................................................................ 19 
2.2. Общие сведения .......................................................................... 19 
2.3. Выполнение работы ................................................................... 20 

2.3.1. Начало работы .................................................................... 20 
2.3.2. Запуск редактора ................................................................ 21 
2.3.3. Сохранение программы ..................................................... 22 
2.3.4. Ввод программы с ручным вводом начальных 

данных ..................................................................................... 22 

2.3.5. Ввод программы с чтением данных из файла .................. 26 

2.4. Содержание отчета ..................................................................... 31 
Литература ......................................................................................... 31 
Контрольные вопросы ....................................................................... 31 
Варианты исходных данных расчетов по работам № 2, № 4, 

№ 7 (вертикальная выработка) .................................................. 32 

3. Расчет в среде Matlab напряжений в массиве пород вокруг 

горизонтальной горной выработки ................................................. 34 
3.1. Цель работы ................................................................................ 34 
3.2. Общие сведения .......................................................................... 34 
3.3. Выполнение работы ................................................................... 35 

3.3.1. Начало работы .................................................................... 35 
3.3.2. Ввод программы с чтением данных из файла .................. 36 

3.4. Содержание отчета ..................................................................... 37 
Контрольные вопросы ....................................................................... 38 
Варианты исходных данных расчетов по работам № 3, № 5, 

№ 6 (горизонтальная выработка) .............................................. 39 

4. Расчеты и построение графиков в среде Mathcad напряжений 

и деформаций в массиве пород вокруг вертикальной горной 
выработки .......................................................................................... 41 
4.1. Цель работы ................................................................................ 41 
4.2. Общие сведения .......................................................................... 41 
4.3. Выполнение работы ................................................................... 43 

4.3.1. Порядок выполнения работы ............................................ 43 
4.3.2. Запуск системы Mathcad и подготовка к вводу 

программы расчета ................................................................. 43 

4.3.3. Ввод текста программы и ее одновременная отладка..... 43 
4.3.4. Получение численных результатов расчета и 

графиков .................................................................................. 48 

4.3.5. Перенос полученных графиков в Word-файл отчета ...... 50 

4.4. Анализ полученных результатов .............................................. 51 
4.5. Содержание отчета ..................................................................... 51 
Контрольные вопросы ....................................................................... 51 

5. Расчеты и построение графиков в среде Mathcad напряжений 

в массиве пород вокруг горизонтальной горной выработки ........ 52 
5.1. Цель работы ................................................................................ 52 
5.2. Общие сведения .......................................................................... 52 
5.3. Выполнение работы ................................................................... 53 

5.3.1. Текст программы с вводом исходных данных 

вручную ................................................................................... 53 

5.3.2. Текст программы с вводом исходных данных из 

файла ....................................................................................... 55 

5.4. Анализ полученных результатов .............................................. 57 
5.5. Содержание отчета ..................................................................... 57 
Контрольные вопросы ....................................................................... 58 

6. Расчет напряжений и деформаций в массиве горных пород 

вокруг горизонтальной горной выработки путем 
моделирования методом конечных элементов в среде 
COMSOL Multiphysics ..................................................................... 59 
6.1. Цель работы ................................................................................ 59 
6.2. Общие сведения .......................................................................... 59 
6.3. Выполнение работы ................................................................... 59 
6.4. Визуализация результатов моделирования .............................. 67 

6.5. Анализ полученных результатов............................................... 70 
6.6. Содержание отчета ..................................................................... 70 
Контрольные вопросы ....................................................................... 70 

7. Расчет напряжений и деформаций в массиве горных пород 

вокруг вертикальной горной выработки путем 
моделирования методом конечных элементов в среде 
COMSOL Multiphysics ..................................................................... 71 
7.1. Цель работы ................................................................................ 71 
7.2. Общие сведения .......................................................................... 71 
7.3. Выполнение работы ................................................................... 71 
7.4. Анализ полученных результатов............................................... 78 
7.5. Содержание отчета ..................................................................... 79 
Контрольные вопросы ....................................................................... 79 

8. Моделирование сейсмического действия взрыва на карьере ......... 80 

8.1. Цель работы ................................................................................ 80 
8.2. Общая информация и исходные параметры расчета .............. 80 

8.2.1. Геометрия модели .............................................................. 80 
8.2.2. Моделирование взрыва ...................................................... 82 

8.3. Среда моделирования ................................................................. 83 
8.4. Порядок выполнения работы..................................................... 83 
8.5. Выполнение работы ................................................................... 84 

8.5.1. Начало работы .................................................................... 84 
8.5.2. Константы и переменные................................................... 85 
8.5.3. Черчение модели по индивидуальным параметрам ........ 86 
8.5.4. Задание свойств элементов модели .................................. 86 
8.5.5. Граничные условия ............................................................ 86 
8.5.6. Задание параметров сетки и разбивка модели на 

элементы ................................................................................. 87 

8.5.7. Параметры решателя (Solver Parameters) и расчет 

модели ..................................................................................... 88 

8.5.8. Построение сейсмограмм на различных расстояниях 

от скважины ............................................................................ 90 

8.5.9. Анализ полученных результатов ...................................... 91 

8.6. Содержание отчета ..................................................................... 93 
Литература ......................................................................................... 93 
Варианты исходных данных для индивидуальных заданий.......... 94 
Контрольные вопросы ....................................................................... 96 

 

ПРЕДИСЛОВИЕ 

Практикум предназначен для первичного освоения студентами 

навыков работы с современным программным обеспечением, используемым для математических расчетов и компьютерного моделирования объектов горного производства, шахтного и подземного 
строительства. В качестве объектов моделирования выбраны вертикальная и горизонтальная подземные выработки, а также карьер. 
Студентам предлагается рассчитать статические напряжения и деформации вокруг них с помощью различных программных средств – 
Matlab, MathCAD, COMSOL Multiphysics. Расчеты производятся по 
аналитическим выражениям механики горных пород, а также методом конечных элементов в различных программных средах, что позволяет проверить их правильность путем сравнения результатов решения схожих задач между собой. Кроме того, при моделировании 
взрывов на карьере приобретаются навыки моделирования динамических механических процессов. 

В первой работе осуществляется подготовка файла будущего от
чета с помощью процессора Word, настраиваются стили и оформляется титульный лист. 

Во второй и третьей работах в среде Matlab по аналитическим 

формулам производится статический расчет напряжений и деформаций в массиве горных пород вокруг горизонтальной и вертикальной 
выработок. В четвертой и пятой работах аналогичные расчеты осуществляются в системе Mathcad. 

Шестая и седьмая работы посвящены аналогичным расчетам и 

моделированию напряжений и деформаций вокруг вертикальной и 
горизонтальной выработок, осуществляемых методом конечных элементов в среде COMSOL Multiphysics. 

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

Для успешного решения практических задач кроме расчетов в 

статике необходимо моделирование также и динамических процессов, протекающих во времени. В восьмой работе на примере карьера 
осуществляется моделирование сейсмического действия взрыва и 
определяется минимальное безопасное расстояние до взрывной 
скважины. 

Указанные работы рассчитаны на выполнение в течение одного 

семестра и составляют вводную часть курса компьютерного моделирования процессов горного производства. Объем работ уточняется 
преподавателем в зависимости от базовой подготовки студентов и 
количества часов, отведенных на выполнение.  

1. Оформление отчетов по лабораторным работам 

и другой документации с помощью процессора Word 

1.1. Цель работы 

Освоение настройки режимов и интерфейса Word для подготовки 

документов большого объема и подготовка файла отчета по практическим занятиям [1]. 

1.2. Настройка стилей 

В процессоре Word по умолчанию в коллекции стилей установ
лено большое их количество, многие из которых не нужны и только 
затрудняют работу. Для подготовки и редактирования большинства 
документов и отчета по НИР достаточно трех базовых стилей текста 
и трех стилей заголовка. Остальные стили являются вспомогательными. При необходимости могут быть добавлены и другие стили. 

При выполнении работы обратите внимание на следующее. Если 

слева находится жирная двойная линия, выполняйте действия. Там, где 
такой линии нет, просто знакомьтесь с изложенной информацией. 

Результаты работы сохраняйте в файл в папку D:\\Study\, в которой со
здайте еще одну папку с номером года (если она отсутствует), а в ней – 
вашу личную папку, в которой будете хранить результаты своих работ. 

На рис. 1.1 показана часть ленты с коллекцией экспресс-стилей. 

Для удобства работы здесь должны содержаться стили Обычный, 
Основной текст, Центр, Заголовок 1, Заголовок 2, Заголовок 3 в 
указанной последовательности. Другие стили рекомендуется удалить 
через правую кнопку мыши, чтобы они не мешали работе. 

 

Рис. 1.1. Окно экспресс-стилей 

В табл. 1.1 приведены рекомендованные стили и их параметры, 

устанавливаемые при настройке. Эти стили настраиваются один раз, 
сохраняются в файле Normal.dot и появляются всякий раз при создании нового документа Word. 

Таблица 1.1 

Рекомендуемые стили для оформления отчета по ГОСТ 7.32–2001 (автоматическое обновление должно быть снято) 

Стиль 
Раздел «Шрифт» 
Раздел «Абзац» 
Другие 

Название 
стиля 
Основан 
на стиле… 

Стиль следующего абзаца 
Шрифт 
Размер 
шрифта 

Видо
доизизмемененение 

Первая 
строка: 
отступ/ 
выступ 

Отступ: 
слева/ 
справа 

Интервал: 
перед/после 

Межстрочный 
интервал 

Выравнивание по 
… 

Добавить в 

коллекцию 

экспрессстилей 

Обычный 
Н/А 
Обычный 
Times New 

Roman 
14 
Нет 
Нет 
0/0 
0/0 
1,5* 
ширине 
Да 

Основной 
текст 
Обычный Основной текст 
То же 
14 
Нет Отступ 1,25 
0/0 
0/0 
1,5 
ширине 
Да 

Центр** 
Обычный 
Центр 
То же 
14 
Нет 
Нет 
0/0 
0/0 
1,5 
центру 
Да 

Заголовок 1 
Обычный Основной текст 
То же 
14 
Нет Отступ 1,25 
0/0 
18/18 пт 
1,5 
центру 
Да 

Заголовок 2 
Обычный Основной текст 
То же 
14 
Нет Отступ 1,25 
0/0 
18/18 пт 
1,5 
центру 
Да 

Заголовок 3 
Обычный Основной текст 
То же 
14 
Нет Отступ 1,25 
0/0 
18/18 пт 
1,5 
центру 
Да 

Оглавление 1 Обычный 
Обычный 
То же 
14 
Нет 
Нет 
0/0,9 
0/0 
1,5 
ширине 
Нет 

Оглавление 2 Обычный 
Обычный 
То же 
14 
Нет 
Нет 
0,42/0,9 
0/0 
1,5 
ширине 
Нет 

Оглавление 3 Обычный 
Обычный 
То же 
14 
Нет 
Нет 
0,99/0,9 
0/0 
1,5 
ширине 
Нет 

Знак концевой сноски 
Основной 

шрифт 
абзаца 

Н/А 
То же 
14 
Нет 
– 
– 
0/0 
– 
– 
Нет 

Текст концевой сноски 
Обычный Текст концевой 

сноски 
То же 
14 
Нет 
Нет 
0/0 
0/0 
1,5 
ширине 
Нет 

_____________ 
* Можно менять для изменения всех стилей. 
** Отсутствует в стандартных стилях, необходимо создать при настройке. 
Н/А – неактивно. 

Для изменения шаблонов стилей на вкладке ГЛАВНАЯ в разделе 

Стили нажмите маленькую кнопку со стрелкой в правом нижнем 
углу. Появится окно Стили (рис. 1.2). Для настройки стилей, их появления и исчезновения в окне Стили нажмите кнопку Управление 
стилями в нижней части (см. рис. 1.2), появится окно Управление 
стилями, выберите вкладку Рекомендации (рис. 1.3). Установите 
меню Порядок сортировки на По алфавиту. 

Рис. 1.2. Окно Стили 

Кнопка  
Создать стиль 
Кнопка Управление 
стилями 

 

Рис. 1.2. Окно Стили 

В этом режиме устанавливаются или скрываются стили в списке 

рекомендованных стилей в окне Стили. Должны быть показаны стили: Обычный, Основной текст, Центр, Заголовок 1, Заголовок 2, 
Заголовок 3, Оглавление 1, Оглавление 2, Оглавление 3, Знак 
концевой сноски, Текст концевой сноски. 

Здесь же устанавливается приоритет появления стилей в окне экс
пресс-стилей на ленте. Проверьте приоритеты стилей: 

Обычный, Основной текст, Центр – приоритет 1; 
Заголовок 1, Заголовок 2, Заголовок 3 – приоритет 10. 
Для создания нового стиля используйте кнопку Создать стиль 

там же. Стиль Центр не включен по умолчанию, если его нет в списке стилей, то создайте его заново.