Математическое и программное обеспечение для численного моделирования геомеханического состояния углепородного массива
Предлагается подход для построения нелинейной математической модели напряженнодеформированного состояния углепородного массива, в которой учитывается интегральное воздействие на горные породы природных и техногенных сил. Для исследования предложенной нелинейной модели разработан комплекс проблемно-ориентированных программ и приведены результаты вычислительных экспериментов. Автор заявляет об отсутствии конфликта интересов.
Mathematical and software for numerical modeling of the geomechanical state of the carbon mass.pdf Геологическая среда - сложный природный объект. Методология описания сложных систем рассмотрена Г. Хакеном, философия которого предусматривает три уровня описания: микро-, мезо-и макроскопический [1]. Поэтому автором предлагается рассматривать геосреду как иерархическую нелинейную систему, декомпозицию которой, в зависимости от задач исследования, необходимо проводить согласно данным трем уровням. Из геосреды на макроуровне выделяется исследуемый участок геомассива. На микроуровне возникает необходимость введения параметров, которые позволили бы описывать наиболее значимые особенности поведения на мезоуровне - это свойства горных пород и действующие на них силы. На мезо- и макроуровне связь между параметрами обеспечивается посредством уравнений, в которых учитываются не сами процессы на микроуровне, а их проявление и влияние на геосреду при воздействии природных и техногенных сил. Такой подход позволяет рассматривать общие сценарии поведения геомассива, включающего систему горных выработок, а также является основой для исследования посредством математического моделирования его геомеханического состояния [1, 2]. Рис. 1. Схема построения математических моделей напряженно-деформированного состояния Fig. 1. Scheme for constructing mathematical models of the stress-strain state Если ставится задача исследования напряженно-деформированного состояния геомассива при воздействии природных и техногенных сил, то можно ограничиться мезоуровнем, считая его представительным. Геомассив по отношению к техногенным воздействиям находится в двух состояниях -природном и нарушенном. В результате горных работ происходит существенное изменение природного состояния горного массива: локализация деформаций и формирование в области нарушения целостности горных пород разгруженных зон и концентраторов напряжений [3]. При таких условиях наблюдаются необратимые деформации, и математическая модель напряженно-деформированного состояния геомассива на мезоуровне должна строиться на основе законов, соответствующих нелинейному характеру деформирования горных пород [4, 5] (рис. 1). 1. Определяющие соотношения математической модели Установлено, что выемка угольного пласта приводит к смещениям части углепородной толщи. На рис. 2 приведены зоны деформирования пород, которые характерны для геомассива, включающего очистные выработки [6]. По результатам натурных исследований выявлено, что в природном состоянии, формирующемся в геологическом масштабе времени, горные породы сохраняют сплошность, и для определения напряженно-деформированного состояния ненарушенной части геомассива D\\ можно применять соотношения теории упругости [3, 7]. 53 Математическое моделирование /Mathematical modeling Рис. 2. Зоны деформирования Fig. 2. Deformation zones Исследователями отмечено, что при техногенных воздействиях породы частично ослаблены трещинами и обладают пластическими свойствами [8]. Следовательно, для зон влияния горных выработок ^2,1, D2,2, D3,i и Ds,2, выделенных на рис. 2 штриховкой, характерна нелинейная зависимость деформаций от напряжений. В соответствии с теорией сдвижения горных пород зоны деформирования D2,i, D2,2, D3,i, D3,2 определяются контурами /1,1, Гі,4, Г2,і, Г2,4. Контурами Г12, Гі,з, /2,2, Г%з задана область очистного выработанного пространства. Расчетная область D представлена в виде объединения областей D1, D2,1, D2,2, D3,1, D3,2: D = D1 ^ D1,1 ^ D1,2 ^ D2,1 ^ D2,2. Для определения напряженно-деформированного состояния геомассива решена следующая краевая задача: в области D найти горизонтальные u = u(x, у) и вертикальные перемещения v = v(x, y), удовлетворяющие системе дифференциальных уравнений о) V (
Ключевые слова
нелинейная математическая модель,
напряженно-деформированное состояние,
массив горных породАвторы
Цветков Андрей Борисович | Сибирский государственный индустриальный университет | доцент, доктор технических наук, профессор кафедры прикладной математики и информатики Института информационных технологий и автоматизированных систем | atsvet@mail.ru |
Всего: 1
Ссылки
Хакен Г. Информация и самоорганизация. Макроскопический подход к сложным системам. М. : Мир, 1991. 240 с.
Ершов Л.В., Максимов В.А. Математические основы физики горных пород. М. : МГИ, 1968. 293 с.
Борисов А.А. Механика горных пород и массивов. М. : Недра, 1980. 360 с.
Цветков А.Б., Павлова Л.Д., Фрянов В.Н. Сравнительная оценка математических моделей геомеханического состояния массива горных пород // Краевые задачи и математическое моделирование : темат. сб. науч. ст. Новокузнецк, 2014. С. 314-317.
Карташов Ю.М., Матвеев Б.В., Михеев Г.В. Прочность и деформируемость горных пород. М. : Недра, 1979. 269 с.
Турчанинов И.А., Иофис М.А., Каспарьян Э.В. Основы механики горных пород. Л. : Недра, 1989. 488 с.
Грицко Г.И., Власенко Б.В., Посохов Г.Е. Прогнозирование и расчет проявлений горного давления. Новосибирск : Наука, 1980. 159 с.
Ставрогин А.Н., Протосеня А.Г. Пластичность горных пород. М. : Недра, 1979. 301 с
Цветков А.Б., Павлова Л.Д. Адаптация алгоритма метода конечных элементов к системе символьной математики // Наукоемкие технологии разработки и использования минеральных ресурсов : Междунар. науч.-практ. конф. : сб. науч. ст. Новокузнецк, 2015. С. 121-125.
Tsvetkov A.B., Pavlova L.D., Fryanov V.N. Numerical simulation of geomechanical state of coal massif in the vicinity of underground workings in the superimposed seams // IOP Conference Series: Earth and Environmental Science. 2017. V. 84. Art. 012005. URL: http://iopscience.iop.org/article/10.1088/1755-1315/84/1/012005/pdf (дата обращения: 22.03.2022).
Цветков А.Б., Павлова Л.Д., Фрянов В.Н. Свидетельство о государственной регистрации программы для ЭВМ № 2020618419. Программа численного исследования математической модели деформирования геомассива с учетом разномодульности горных пород // Реестр программ для ЭВМ. Дата регистрации 27.07.2020.
Штумпф Г.Г. Физико-технические свойства горных пород и углей Кузнецкого бассейна. М. : Недра, 1994. 447 с.