CC BY
49
-------------------------------- © П.В. Макаров, И.Ю. Смолин,
Е.П. Евтушенко, А.А. Трубицын,
Н.В. Трубицына, С.П. Ворошилов,
С.И. Голоскоков, Д.В. Ботвенко,
2009
УДК 622.8
П.В. Макаров, И.Ю. Смолин, Е.П. Евтушенко,
А.А. Трубицын, Н.В. Трубицына, С.П. Ворошилов,
С.И. Голоскоков, Д.В. Ботвенко
ВОЗМОЖНОСТИ ПРИМЕНЕНИЯ ЭВОЛЮЦИОННОГО ПОДХОДА ПРИ МОДЕЛИРОВАНИИ ПОВЕДЕНИЯ ГЕОСРЕДЫ
Представлено решение задачи о механическом поведении геосреды вокруг выработки с точки зрения эволюционного подхода.
Ключевые слова: моделирование, геосреда, горная выработка.
Т} азличные подземные полости антропогенного происхож--мГ дения (горные выработки, подземные сооружения) могут представлять серьезную угрозу, как для наземных сооружений, так и для безопасного ведения горных работ. Поэтому одной из актуальных проблем является задача управления кровлей в забоях угольных шахт. В последнее десятилетие эта проблема стала особенно острой в связи со значительным увеличением нагрузок на забои. Высокие скорости их подвигания приводят к тому, что горный массив вокруг выработки и, особенно, в породах кровли оказывается в существенно неравновесном и неустойчивом состоянии [1]. Это приводит к зависанию протяженных участков кровли. Значительное возрастание шагов первой посадки кровли (до 100 м и более) и последующих посадок (до нескольких десятков метров) грозит катастрофическими неуправляемыми обрушениями кровли и представляет серьезную угрозу безопасности при ведении горных работ. Наряду с уже разработанными технологическими способами и приемами управления кровлей вблизи очистного забоя [2, 3], понимание основных закономерностей эволюции горного массива вокруг выработки и возможных сценариев развития событий в сверхбыстрых, катастрофических режимах: обрушений кровли, различных горных обвалов, формирование макроскопических раз-
рывов и т.д., представляет собой актуальную фундаментальную задачу геомеханики, решение которой должно стать научной основой для многих горных технологий.
Механическое поведение кровли над выработанным пространством, и всего горного массива вокруг выработки, включая все значимые элементы геосреды, авторы рассматривают с эволюционной точки зрения [4-6]. Согласно развиваемому подходу, геосреда в поле сил тяжести, а также при других возможных силовых воздействиях, эволюционирует по законам общей теории эволюции, которая к настоящему времени формируется в рамках нелинейной динамики или синергетики [4-6]. Однако привлечение уже ставших традиционными методов и подходов нелинейной динамики, которые основываются на анализе общих свойств решений базовых уравнений синергетики, оказалось практически мало эффективным для выяснения сценариев эволюции конкретных физических и природных систем. И это понятно, базовые уравнения синергетики не являются математическими моделями изучаемых систем, а схожесть общих закономерностей любого эволюционного процесса ничего нам не говорит о конкретном эволюционном сценарии исследуемой системы. Следовательно, для каждого физического или природного процесса должна быть создана оригинальная эволюционная модель, которая детально описывает изучаемый процесс.
В работах [4-6] показано, что если задача моделирования физического процесса сформулирована как эволюционная, то в силу общности характерных свойств решений систем нелинейных уравнений и решения этой задачи, а в более общем случае уравнений математической физики, в частности уравнений механики сплошных сред, выражающих законы сохранения, вместе с нелинейными определяющими уравнениями первой и второй группы (эволюционными уравнениями), проявляют все свойства решений нелинейных динамических систем, а также качественные закономерности, достаточно подробно изученные при анализе общих решений базовых уравнений синергетики. Это, прежде всего, иерархичность и многомасштабность изучаемых физических и природных процессов, масштабная инвариантность, самоподобие или самоафинность, пространственная локализация распределений параметров (деформаций и повреждений при развитии деформационных процессов), пороговость физических и природных явлений, когда, например, при переходе на новый масштаб разрушения должен быть преодо-
лен характерный энергетический порог; неравномерность процессов эволюции во времени - наличие медленной квазистационарной стадии и последующая сверхбыстрая эволюция в режиме с обострением и т.д. [4-9]. Также было показано, что решения соответствующих уравнений механики деформируемого твердого тела демонстрируют все или некоторые из перечисленных особенностей только в том случае, когда задача сформулирована как эволюционная и граничные условия, например, условия нагружения, обеспечивают для управляющих параметров достижение ими пороговых значений, при превышении которых в системе наблюдаются изменения пространственно временной симметрии в распределениях параметров, т.е. происходит самоорганизация - образование новых нестационарных диссипативных структур [4, 5]. В работе [1] предложено решение задачи об обрушении кровли как задачи эволюционной, подобраны скорости накопления повреждений на рассматриваемых пространственно временных масштабах, что обеспечило приемлемое согласие с натурными данными. В настоящей работе развивается эта концепция.
Основные задачи настоящего исследования носят также в основном методический характер, хотя параметры, обеспечивающие эволюцию геосреды в поле сил тяжести, подбирались по возможности в соответствии с данными по посадкам кровли для реальных забоев Кузбасса [1]. Основное внимание в работе уделяется общей феноменологии явления без привязки к разнообразным физическим процессам и механизмам, могущим оказать кардинальное влияние на реализацию конкретного сценария эволюции или смену эволюционных сценариев. Например, обводненность среды, наличие разломов, слоистость, другая исходная структурная неоднородность и т.д. могут быть в дальнейшем учтены при решении конкретных прикладных задач в рамках развиваемой концепции. Отметим также следующий фундаментальный аспект, кардинально меняющий сценарий разрушения при решении задачи с позиций развиваемого подхода в отличие от традиционного усредненного макроскопического описания. С позиций традиционной макроскопической точки зрения сплошная деформируемая среда представляет собой один заранее выбранный (продиктованный решаемой задачей) уровень или масштаб макроскопического феноменологического описания вне зависимости от физического содержания и механизмов, которые могут быть отражены в уравнениях состояния. Эволюционная
точка зрения изначально предполагает наличие в среде иерархии масштабов [1, 4, 5], которые и выявляются в среде по мере развития деформационных процессов, формируя в нагружаемом материале иерархии нестационарных диссипативных структур, например, систем трещин разных порядков в хрупком материале или иерархию сравнительно слабодеформированных блоков, разделенных системами сопряженных полос локализованной деформации [4, 5]. Понятно, что при численном моделировании минимальные размеры таких структурных элементов ограничены возможностями счета (дискретными шагами по пространству). Другими словами, эволюционный подход описывает процессы самоорганизации - формирование в нагружаемой геосреде новых нестационарных структур, их взаимодействия и трансформации. Эти структуры не могут быть навязаны среде извне, а генетически присущи самой нелинейной среде. По этой причине нелинейные определяющие уравнения включают две группы: первая группа определяющих уравнений выражает связи между макропараметрами; вторая группа эволюционных уравнений представляет собой кинетические уравнения, задающие генерацию параметров в зависимости от событий микроуровня. Таким образом, в нагружаемой геосреде, в том числе и первоначально структурно-неоднородной, формируется неоднородное распределение параметров. Эти неоднородности образуются в среде в силу внутренних причин - нелинейных свойств среды, а не в силу специально подобранных условий нагружения. Это и есть по И. Пригожину [7] диссипативные структуры.
Следуя работам [1, 4, 5] выпишем систему уравнений механики твердого деформируемого тела, которая вместе с определяющими или эволюционными уравнениями является в общем случае нелинейной и открытой. Она описывает эволюцию любых твердых тел, в том числе и нагружаемой геосреды. Связь этой открытой нелинейной системы с внешним миром осуществляется через управляющие параметры, которые совместно с граничными условиями (их также в ряде случаев можно причислить к управляющим параметрам) определяют условия нагружения исследуемого объема геосреды [1, 4, 5].
Такая эволюционирующая открытая система будет способна к самоорганизации только в том случае, если в ней имеются объемные, распределенные в пространстве источники, способные генерировать и перераспределять неоднородности в распределении па-
раметров, фактически формировать новые и трансформировать уже имеющиеся структуры: тепло (тепловые структуры [7]), дефекты и повреждения в деформируемых твердых телах [4] и т.д. Диссипативные нестационарные структуры, способные к развитию и приводящие к усложнению и совершенствованию системы, т.е. ее адаптации к внешним воздействиям, формируются источниками [8], в то время как на стоках возможно формирование только стационарных структур [4, 8]. По этой причине деформируемое твердое тело - это идеальная среда, способная к самоорганизации, так как в ней под приложенными воздействиями генерируются различные дефекты и повреждения разных масштабов [4]. Движущими силами эволюции геосреды являются сила тяжести и внутреннее тепло. Для решения задачи эволюции горного массива вокруг выработки ограничимся только гравитацией.
Таким образом, полная система нелинейных уравнений при лагранжевом подходе к описанию движения сплошной среды включает:
уравнения, выражающие законы сохранения д<5н да
длр£^
РV = Р0*0 > РУг + РГ> ’ РЕ = °г}8г] ; (1)
эволюционные определяющие уравнения первой группы 6у = М0Т ~ 0Р )5у + 2М*(8у ~ 8у ) > (2)
эволюционные уравнения второй группы
0у=а 0У+с(0)’ іу=^ ■ б* • ^ ) ■ (3)
Здесь р о , р — начальное и текущее значение плотности материала; У0, V — начальное и текущее значение объёма некоторой частицы материала; Х1 — координаты в декартовой системе координат наблюдателя; Уг — компоненты вектора скорости перемещений; ог]- — компоненты тензора напряжений; ^ — компоненты
1 ( ду,- ду,л
— компоненты тензора
вектора массовых сил; 8, = —
■Т =
= 2
+ —] удх] дх у
скорости деформации; Е — внутренняя энергия единицы начально-
го объёма; qi — компоненты вектора теплового потока; точка над
символом означает материальную производную по времени. Используется разложение полной скорости деформации на упругую и
•Т ■ Е ■ Р тт
неупругую составляющие 8- = 8- + 8- . Для учета независимости
от жесткого вращательного движения производная по времени от напряжений записана в форме коротационной производной Яуман-■ J ■
на Оу = Су +С^(йк- -ск-®^ ,
1 ( дуг ду]л
Кдх] дхг
компоненты тензора скорости враще-
ния (вихря). А, В, С, F — некоторые функции, определяемые при конкретном задании кинетик (3). Использованы также следующие обозначения
0Т = 8Т, 0Р =8у, 8* = £ НУ =Я еУеУ , еу = 8у - 30РЪ]
Р
где в- - компоненты девиатора тензора скорости пластических деформаций, S - — компоненты девиатора тензора напряжений, Р
— среднее давление, 5- — символ Кронеккера.
Эволюционные уравнения первой группы (2) записаны в релаксационной форме и выполняют несколько важнейших функций в развитии эволюционного процесса:
1. Упругие напряжения растут в каждой локальной точке нагружаемой среды пропорционально росту полных скоростей де-
• Т
формации 8- , т.е. в соответствии со скоростью нагружения среды.
Такое нагружение может быть осуществлено либо непосредствен-
• Т
ным заданием полной скорости деформирования 8- , либо косвенно через приложенные массовые и поверхностные силы и перераспределения параметров в соответствии с развитием течения среды, определяемом законами сохранения (1). Именно в этом смысле
скорости полных деформаций 8- выполняют роль управляющих
параметров, согласно терминологии Г. Николиса и И. Пригожина [7]. В любом случае, внутри расчетной области они сложным образом зависят от параметров течения среды.
2. Релаксируют напряжения в соответствии с развитием неупругих реакций в нагружаемой среде 8Р . Следовательно, релаксационные определяющие уравнения первой группы обеспечивают установление в каждой точке деформируемой среды динамического равновесия между внешними воздействиями и откликом среды на нагружение.
3. Через релаксационные определяющие уравнения первой группы реализуются отрицательные и положительные обратные связи. Отрицательная обратная связь реализуется через закон релаксации и стабилизирует деформационный процесс, выравнивая неоднородности. Действительно, чем выше напряжения, тем выше скорость генерации неупругих сдвигов и повреждений, а значит и выше скорость релаксации (2). С уменьшением напряжений уменьшается и генерация повреждений и сдвигов. Процесс стабилизируется в состоянии динамического равновесия в данной частице среды, из которого его может вывести, например, изменение
полной скорости деформации 8- , как управляющего параметра.
Понятно, что эта стабилизация определяется не только скоростью генерации неупругих деформаций и повреждений, а также механизмами упрочнения и деградации физико-механических характеристик среды, о чем будет сказано далее. Положительная обратная связь не только разгоняет автокаталитический процесс деградации прочностных характеристик среды - локализация неупругих деформаций (повреждений) снижает прочностные характеристики среды в этих локальных областях, что в свою очередь, усиливает в них процессы локализации, - но и выполняет созидающую роль, формируя в среде новые диссипативные деформационные структуры [4, 8].
Эволюционные уравнения второй группы (3) отражают реакцию среды на нагружение. В то время как полные скорости дефорТ
мации 8- , так или иначе, определяются условиями нагружения,
8р и 0Р зависят только от внутренних свойств среды и есть результат работы распределенных объемных источников, генерирую-
щих в среде дефекты. Кинетические уравнения (3) интегрально отражают масштабы микроуровня и определяют возможности среды перераспределять и диссипировать подводимую к ней энергию. В этом смысле они отвечают за способность нагружаемой среды к формированию в ней диссипативных нестационарных структур, т.е. за самоорганизацию. В этом и заключается созидательная роль положительной обратной связи [8]. С другой стороны, эти уравнения описывают процессы накопления в среде дефектов и повреждений, и тем самым полностью определяют деградацию физико-механических свойств среды.
Входя в эволюционные уравнения первой группы (2) 8р и 0Р участвуют в игре между положительными и отрицательными обратными связями, определяя эволюционные сценарии.
Модель нагружаемой геосреды обсуждена в работах [1, 4, 5] и учитывает внутреннее трение, дилатансию, накопление повреждений и деградацию прочностных характеристик геосреды. Предельная поверхность перехода от упругой реакции среды к неупругому поведению определена по модели В.Н. Николаевского [9]
у^ = Г, 0Р = /Р = 2Л/21/2 = 3А8Р^, (4)
где а - коэффициент внутреннего трения, Л - коэффициент дила-тансии. Изменения прочностных параметров модели описываются следующими уравнениями [1, 5]:
Г = Го(1 + А(в) - Ве (в))(1 - В(с)), А (в) = 2h —,
в *
Вв(в) = 2/1 в] = А(в)-в*, О(с) = |<С- С")' Л
^ в * ) в * ■> (С*)*
для а > а0 (5)
Здесь к - параметр модели, отвечающий за упрочнение, в = 2/р 12 = 38 ^ - накопленная пластическая деформация (параметр Одквиста), в* - критическая деформация, после достижения которой преобладает деградация материала, о - эффективное напряжение, а о0, о*, (* и п = 2 - параметры модели, определяющие пороговое напряжение, начиная с которого начинают накапливаться повреждения, предельное напряжение и скорость процесса накопления повреждений, В(о) — поврежденность среды.
Усечение предельной поверхности в области растяжения также определяется с учетом повреждений и отражает различие прочностных свойств геосреды при сжатии и растяжении. Фактически среда всегда разрушается в области растягивающих напряжений, абсолютная величина которых невелика и уменьшается с ростом повреждённости среды:
Р*= Ро* (1 - В(а)). (6)
Таким образом, принимается, что разрушение и раскрытие трещин происходит при наличии растягивающих напряжений, которые присутствуют в локальных областях неоднородной среды практически при любом виде нагрузок.
Выбор степенных зависимостей в выражениях (5) был продиктован соображениями, связанными с изучением локализации и образованием тепловых структур и глубоким замечанием С.П. Курдюмова о том, что не любая нелинейность приводит к возникновению структур, т.е. к самоорганизации. Сложные нестационарные структуры формируются в среде со степенными источниками. При других видах нелинейностей задача часто вырождается при приближении к моменту обострения [8]. Степенные зависимости в законах накопления повреждений следуют также из идеи самоподобия (или самоафинности) [10, 11].
Реальное время эволюционного процесса Т существенно больше возможного времени t непосредственно счета [1]. В данном случае Т задается характерными временами подвигания забоя. За время АТ забой продвигается на один дискретный пространственный интервал Ах. Шаг по времени процесса АТ выбирается из условия
АТ = —, (7)
с
где п = 3, Ь - размер расчетной области, с - скорость звука в геосреде. Условие (7) обеспечивает установление в геосреде практически стационарного распределения параметров после малого возмущения, вызванного увеличением выработанного пространства. Скорости накопления повреждений в геосреде выбираются по тестовым расчетам в соответствии с известными данными по реальным временам посадки кровли [1], которые составляют от нескольких дней до десятков дней.
Одной из фундаментальных проблем геомеханики, которая пока не нашла должного объяснения, является проблема изменения отклика геосреды от вязкой реакции до хрупкого поведения в зависимости от условий нагружения [12, 13]. Такое поведение геосреды является вполне объяснимым в рамках развиваемой эволюционной концепции. В зависимости от сценария эволюции, который определяется как нелинейными свойствами среды, скоростью генерации повреждений, так и балансом положительных (разгон автокаталитического процесса) и отрицательных (стабилизация деформационного процесса) обратных связей, среда может реагировать от типично вязкопластического поведения до хрупкого разрушения [4]. На рис. 1 представлены результаты расчета развития повреждений в кровле для разных скоростей генерации повреждений. Изменение только одного параметра * в законе генерации (5) полностью изменило эволюционный сценарий. На рис. 1, а мы видим развитое вязкопластическое течение, а на рис. 1, б - типичное фрактальное дерево ветвящихся трещин. В полосах локализованной неупругой деформации (повреждений) (рис. 1, а) накопленные неупругие деформации более чем на порядок превышают деформации при хрупком растрескивании на рис. 1, б. В нагруженной геосреде сформировались диссипативные структуры - блоки.
В зависимости от условий нагружения (скорости подвигания забоя) меняется и сценарий обрушения кровли. При малых скоростях подвигания забоя (рис. 2) можно наблюдать равновесный режим [1], когда обрушения кровли происходят вслед за продвижением забоя. При высоких скоростях продвижения забоя наблюдаются протяженные участки зависшей кровли, реализуется существенно неравновесный режим, при котором поврежденность среды в кровле достигает критических значений только тогда, когда забой успевает уйти на значительные расстояния (рис 3).
Рис. 1. Распределение повреждений вокруг горной выработки для вязкого (а) и хрупкого (б) разрушения в кровле
Подобные неравновесные режимы чрезвычайно опасны и грозят катастрофическими обрушениями больших участков кровли. Спусковыми механизмами подобных катастрофических обрушений могут быть различные причины (обводнение, наличие газов, быстрый рост исходных трещин и разломов в горном массиве). Все эти причины могут перевести медленную стадию эволюции повреждений в кровле в сверхбыстрые стадии в режимах с обострением [4, 8]. Эти вопросы требуют специального пристального изучения и в настоящей работе не рассматриваются. Мы рассматриваем тот случай, когда исключительно за счет нелинейных свойств среды без каких либо провокаций
10 суток 12 суток 14 суток
Рис. 2. Стационарный режим обрушения кровли. Кровля обрушается вслед за продвижением забоя. Длина расчетной области 80 м,, выработки - 50 м
12 суток 14 суток 16 суток
Рис. 3. Неравновесный режим обрушения с образование протяженных участков зависшей кровли. Длина расчетной области 80 м, выработки - 50 м
Рис. 4. Первая и последующие посадки кровли для случая, глубина выработки 200 м, длина выработки 300 м
извне медленная стадия эволюции повреждений в кровле сменяется сверхбыстрым режимом с обострением.
На рис. 4 показаны результаты моделирования первой посадки кровли и нескольких последующих. Расчеты выполнены для типичной ситуации, наблюдаемой в реальных условиях [1]. Шаг последующих обрушения приблизительно в два раза меньше, чем протяженность первой посадки. Каждая последующая посадка меньше предыдущей, так как увеличение пролета выработанного пространства ускоряет процесс накопления повреждений в кровле. На всех рис. 1-4 мы видим также развитие повреждений в почве и ее выпучивание в выработанное пространство. В любой точке, где функция поврежденности достигает единицы, и среда оказывается разрушенной полностью, т.е. материал исчерпывает все возможности сопротивляться разрушению и его прочность становится нулевой, реализуется типичный эволюционный сценарий — предшествующая длительная квазистационарная стадия медленного накопления повреждений и медленной деградации прочности сменяется на режим сверхбыстрой эволюции. Понятно, что для разных точек среды продолжительность этих стадий различна в зависимости от эволюции напряженно-деформированного состояния среды в этой конкретной точке. Подобное поведение материалов и сред хорошо известно из экспериментов [13].
Так как эволюция геосреды происходит практически в стесненных условиях, обусловленных горным давлением и смягчаемых только наличием выработанного пространства, роль дилатансион-ных процессов в развитии повреждений в горном массиве огромна. Роль дилатансии в эволюции напряженно-деформированного состояния в геосредах, а также в связи с развитием различных ката-
строфических событий, широко обсуждается в литературе [13, 15, 16].
Подробному рассмотрению этой сложной проблемы будет посвящена отдельная работа. Здесь остановимся только на нескольких вопросах принципиального характера. Ещё М.М. Протодьяко-новым [15] было замечено, что в геосреде начиная с некоторой глубины, где горное давление создает напряженное состояние в среде, близкое к предельным величинам для этого материала, в геосреде за счет дилатансионных процессов формируются арки, принимающие на себя давление вышележащих пород, и с дальнейшим ростом глубины выработки горное давление перестает расти. Наши расчеты также демонстрируют формирование в поврежденной среде подобных арок (рис. 1-4). Пока мы говорим только о качественной стороне проблемы, так как для выяснения количественных характеристик необходимы экспериментальные данные о величинах дилатансии для конкретных геосред в зависимости от приложенных нагрузок. Тестовые расчеты показали, что картина разрушения кардинально меняется в зависимости от принятой модели нарастания дилатансии. При большой дилатан-сии возможна ситуация, когда большинство периферийных трещин в кровле запирается, их рост замирает, и развиваются несколько или одна магистральная трещина в центральной части кровли.
Понимание того факта, что разрушение есть многомасштабный иерархический процесс, а нагружаемая среда неоднородна во всей иерархии масштабов, хотя бы потому, что в ней всегда присутствует исходная поврежденность, все ставит на свои места. По мере накопления повреждений в каждой точке геосреды средняя ее прочность только падает от некоторого исходного значения до нуля, и в этом случае значение функции поврежденности D(a) = 1 (рис. 5). При среднем макроскопическом напряжении в точке А, равном оА, локальные напряжения на более низких масштабах могут существенно превышать и средние напряжения, и величину прочности материала, что и приводит к появлению повреждений разных масштабов, падению средней макроскопической прочности и росту среднего напряжения в этой точке за счет дилатансии. Кривые оА и оМ сближаются. Когда оА = оМ, то только в этой локальной точке А, а не во всем образце, будет иметь место
Рис. 5. Зависимость усредненной макроскопической прочности гесреды от накопленных повреждений: оМ - макроскопическая прочность геосреды; оА - среднее напряжение в некоторой точке А частично поврежденной геосреды; ое - локальные напряжения на промежуточных масштабах
локальное макроскопическое разрушение. В этом случае величина поврежденности D(o) = 1 в локальном макроскопическом объеме А. В рассматриваемой эволюционной модели все нижележащие уровни учтены усреднено через функцию поврежденности D(o). В уравнении (5) для D(o) при малой величине о0, меньшой макроскопического предела текучести У, модель позволяет накапливать повреждения и дилатансию уже на макроскопически упругой стадии, но не упругой на малых масштабах. Такой подход на взгляд авторов позволяет качественно правильно учитывать дилатансию и описывать поведение хрупкой геосреды в стесненных условиях деформирования.
Таким образом, в работе задача о механическом поведении геосреды вокруг выработки решается с точки зрения развиваемого авторами эволюционного подхода. Расчеты эволюции повреждений в кровле над выработанным пространством показали, что в кровле
формируется иерархия блоков разных масштабов, разделенных полосами локализованных повреждений. Показано также, что в зависимости от конкуренции отрицательных обратных связей, стабилизирующих деформационный процесс и сглаживающих неоднородности в распределениях параметров, и положительных обратных связей, разгоняющих автокаталитический процесс деградации нагружаемой среды, сценарий эволюции геосреды может меняться от типичного вязкопластического течения до сугубо хрупкой реакции. Отмечена важная роль дилатансии в развитии деформации в стесненных условиях, показано, что дилатан-сионные процессы приводят к образованию в поврежденной среде арочных форм. Дана критика концепции дилатансионного упрочнения.
------------------------------------------- СПИСОК ЛИТЕРАТУРЫ
1. Макаров П.В., Смолин И.Ю., Евтушенко Е.П., Трубицын А.А., Трубицына Н.В., Ворошилов С.П. Моделирование обрушения кровли над выработанным пространством // Физ. мезомех. - 2008. - Т. 11. - № 1. - С. 44-50.
2. Якоби О. Практика управления горным давлением. - М.: Недра, 1987.
- 566 с.
3. Никишин Д.Ю. Обоснование эффективных способов предотвращения обрушений пород кровли в лавах надрабатываемых слоев: дис. ... канд. техн. наук. - С.-Петербург, 2007 - 153 с.
4. Макаров П.В. Математическая теория эволюции нагружаемых твердых тел и сред // Физ. мезомех. - 2008. - Т. 11. - № 3. - С. 19-35.
5. Макаров П.В., Смолин И.Ю., Стефанов Ю.П., Кузнецов П.В., Трубицын А.А., Трубицына Н.В., Ворошилов С.П., Ворошилов Я.С. Нелинейная механика геоматериалов и геосред. - Новосибирск: Академич. изд-во «Гео», 2007.
- 235 с.
6. Макаров П.В. Эволюционная природа блочной организации геоматериалов и геосред. Универсальный критерий фрактальной делимости // Геология и геофизика. - 2007. - Т. 48. - № 7. - С. 724-746.
7. Николис Г., Пригожин И. Познание сложного. - М. Едиториал УРСС, 2003. - 344 с.
8. Курдюмов С.П., Князева Е.Н. Законы эволюции и самоорганизации сложных систем. - М.: Наука, 1994. - 236 с.
9. Гарагаш И.А., Николаевский В.Н. Неассоциированные законы течения и локализации пластической деформации // Успехи механики. - 1989. - Т. 12.
- № 1. - С. 131-183.
10. Малинецкий Г.Г., Потапов А.Б. Современные проблемы нелинейной динамики. - М. Едиториал УРСС, 2002. - 356 с.
11. Карпинтери А., Лачидонья Дж., Пуцци С. Прогноз развития трещин в полномасштабных конструкциях на основе анализа показателя Ь и статистики Юла // Физ. мезомех. - 2008. - Т. 11. - № 3. - С. 75-87.
12. Наймарк О.Б. Структурно-скейлинговые переходы и автомодельные закономерности развития землетрясений // Физ. мезомех. - 2008. - Т. 11. - №
2. - С. 89-106.
13. Гольдин С.В. Дилатансия, переупаковка и землетрясения // Физика Земли. - 2004. -№ 10. - С. 37-54.
14. Си Дж. Мезомеханика, понятие сегментации и мультискейлинговый подход: нано-микро-макро // Физ. мезомех. - 2008. - Т. 11. - № 3. - С. 5-18.
15. Протодьяконов М.М. Давление горных пород и рудничное крепление.
Ч. 1. Давление горных пород. Изд. 3-е. - М.: Госгориздат, 1933. - 127 с.
16. Шемякин Е.И. Механика горного массива // Научные сообщения Института горного дела им. А.А. Скочинского. -2007. - Вып. 333. - С. 137-152.
Makarov P. V., Smolin I.J., Evtushenko E.P., Trubitsyn A.A.,
Trubitsyna N. V., Voroshilov S.P., Goloskokov S.I.,
Botvenko D. V.
POSSIBILITIES OF EVOLUTIONARY APPROACH APPLICATION AT MODELLING OF GEO-ENVIRONMENT BEHAVIOUR
In this work it is presented the decision of a problem on mechanical behaviour of geo-environment round development from the point of view of the evolutionary approach.
Key words: modelling, geo-environment, mining development.
— Коротко об авторах --------------------------------------------------
Макаров П.В. - доктор физико-математических наук, проф., ИФПМ СО РАН, г. Томск,
Смолин И.Ю. - кандидат физико-математических наук, ИФПМ СО РАН, г. Томск,
Евтушенко Е.П. - Томский государственный университет, г.Томск, Трубицын А.А. - доктор технических наук, профессор, НЦ ВостНИИ, г. Кемерово,
Трубицына Н.В. - доктор технических наук, НЦ ВостНИИ, г. Кемерово, Ворошилов С.П. - кандидат физико-математических наук, НП «Куз-басс-ЦОТ», г. Кемерово,
Голоскоков С.И. - кандидат технических наук, НЦ ВостНИИ, г. Кемерово,
Ботвенко Д.В. - кандидат технических наук, НЦ ВостНИИ, г. Кемерово.
