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

У статті обґрунтовано необхідність використання тривимірної моделі шаруватого трансверсально-ізотропного масиву гірських порід для вивчення та уточнення закономірностей розподілу напружено-деформованого стану навколо очисного вибою. Обґрунтовано параметри тривимірного моделювання масиву гірських пор...

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Datum:2009
1. Verfasser: Сидельников, А.А.
Format: Artikel
Sprache:Russian
Veröffentlicht: Інститут геотехнічної механіки імені М.С. Полякова НАН України 2009
Schriftenreihe:Геотехническая механика
Online Zugang:http://dspace.nbuv.gov.ua/handle/123456789/32984
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
Назва журналу:Digital Library of Periodicals of National Academy of Sciences of Ukraine
Zitieren:Обоснование параметров объемного моделирования массива горных пород вокруг очистной и подготовительных выработок / А.А. Сидельников // Геотехническая механика: Межвед. сб. науч. тр. — Днепропетровск: ИГТМ НАНУ, 2009. — Вип. 82. — Бібліогр.: 9 назв. — рос.

Institution

Digital Library of Periodicals of National Academy of Sciences of Ukraine
id irk-123456789-32984
record_format dspace
spelling irk-123456789-329842012-05-27T12:27:54Z Обоснование параметров объемного моделирования массива горных пород вокруг очистной и подготовительных выработок Сидельников, А.А. У статті обґрунтовано необхідність використання тривимірної моделі шаруватого трансверсально-ізотропного масиву гірських порід для вивчення та уточнення закономірностей розподілу напружено-деформованого стану навколо очисного вибою. Обґрунтовано параметри тривимірного моделювання масиву гірських порід для умов шахт Західного Донбасу. The necessity of use a three-dimensional model of a layered transversally-isotropic rock mass for studying and specification stress and strain distribution laws around longwalls and development faces is substantiated in the article. Parameters of the three-dimensional simulation of a rock mass for Western Donbas conditions are substantiated. 2009 Article Обоснование параметров объемного моделирования массива горных пород вокруг очистной и подготовительных выработок / А.А. Сидельников // Геотехническая механика: Межвед. сб. науч. тр. — Днепропетровск: ИГТМ НАНУ, 2009. — Вип. 82. — Бібліогр.: 9 назв. — рос. http://dspace.nbuv.gov.ua/handle/123456789/32984 622.063.4 ru Геотехническая механика Інститут геотехнічної механіки імені М.С. Полякова НАН України
institution Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection DSpace DC
language Russian
description У статті обґрунтовано необхідність використання тривимірної моделі шаруватого трансверсально-ізотропного масиву гірських порід для вивчення та уточнення закономірностей розподілу напружено-деформованого стану навколо очисного вибою. Обґрунтовано параметри тривимірного моделювання масиву гірських порід для умов шахт Західного Донбасу.
format Article
author Сидельников, А.А.
spellingShingle Сидельников, А.А.
Обоснование параметров объемного моделирования массива горных пород вокруг очистной и подготовительных выработок
Геотехническая механика
author_facet Сидельников, А.А.
author_sort Сидельников, А.А.
title Обоснование параметров объемного моделирования массива горных пород вокруг очистной и подготовительных выработок
title_short Обоснование параметров объемного моделирования массива горных пород вокруг очистной и подготовительных выработок
title_full Обоснование параметров объемного моделирования массива горных пород вокруг очистной и подготовительных выработок
title_fullStr Обоснование параметров объемного моделирования массива горных пород вокруг очистной и подготовительных выработок
title_full_unstemmed Обоснование параметров объемного моделирования массива горных пород вокруг очистной и подготовительных выработок
title_sort обоснование параметров объемного моделирования массива горных пород вокруг очистной и подготовительных выработок
publisher Інститут геотехнічної механіки імені М.С. Полякова НАН України
publishDate 2009
url http://dspace.nbuv.gov.ua/handle/123456789/32984
citation_txt Обоснование параметров объемного моделирования массива горных пород вокруг очистной и подготовительных выработок / А.А. Сидельников // Геотехническая механика: Межвед. сб. науч. тр. — Днепропетровск: ИГТМ НАНУ, 2009. — Вип. 82. — Бібліогр.: 9 назв. — рос.
series Геотехническая механика
work_keys_str_mv AT sidelʹnikovaa obosnovanieparametrovobʺemnogomodelirovaniâmassivagornyhporodvokrugočistnojipodgotovitelʹnyhvyrabotok
first_indexed 2025-07-03T13:27:48Z
last_indexed 2025-07-03T13:27:48Z
_version_ 1836632525896679424
fulltext "Геотехническая механика" УДК622.063.4 А.А. Сидельников, асп. (Национальный горный університет) ОБОСНОВАНИЕ ПАРАМЕТРОВ ОБЪЕМНОГО МОДЕЛИРОВАНИЯ МАССИВА ГОРНЫХ ПОРОД ВОКРУГ ОЧИСТНОЙ И ПОДГОТОВИТЕЛЬНЫХ ВЫРАБОТОК У статті обґрунтовано необхідність використання тривимірної моделі шаруватого тра- нсверсально-ізотропного масиву гірських порід для вивчення та уточнення закономірнос- тей розподілу напружено-деформованого стану навколо очисного вибою. Обґрунтовано параметри тривимірного моделювання масиву гірських порід для умов шахт Західного Донбасу. SUBSTANTIATION PARAMETERS OF THE ROCK MASS 3-D SIMULATION AROUND THE LONGWALLS AND DEVELOPMENT FACES The necessity of use a three-dimensional model of a layered transversally-isotropic rock mass for studying and specification stress and strain distribution laws around longwalls and develop- ment faces is substantiated in the article. Parameters of the three-dimensional simulation of a rock mass for Western Donbas conditions are substantiated. Введение. На сегодняшний день основной особенностью ведения горных работ в украинской части донецкого угольного бассейна является то, что средняя глубина разработки превышает 700 м, максимальная при этом дости- гает 1200 м. В связи с этим, определился широкий круг горнотехнических проблем связанных с негативными проявлениями горного давления в подго- товительных выработках и очистных забоях, которые, в свою очередь, суще- ственно снижают добычу угля, увеличивают себестоимость и значительно повышают уровень производственного травматизма. Применение как актив- ных, так и пассивных способов управления горным давлением в выработках требует знания закономерностей распределения напряженно- деформированного состояния массива горных пород вдоль лавы и вдоль вы- емочного столба. Суть вопроса. На современном уровне развития горной науки и методов изучения НДС вокруг очистных и подготовительных выработок основным и наиболее эффективным способом прогнозирования состояния массива гор- ных пород является метод моделирования геомеханических процессов, кото- рые происходят вследствие ведения горных работ. Данный вопрос исследова- ли при помощи физических, математических, феноменологических, компью- терных и других моделей такие ученые как А.А. Борисов, А.Н. Динник, М.А. Иофис, Г.Н. Кузнецов, Л.В. Новикова, А.В. Савостьянов, В.Д. Слесарев, М.С. Четверик, А.Н. Шашенко. Однако любая модель – это детерминирован- ное представление того или иного реального объекта или процесса и чем меньше будет использовано упрощений при построении модели, тем более она будет соответствовать реальности. Обзор современных подходов к моде- лированию геомеханических процессов, протекающих в массиве горных по- Выпуск № 82 род, позволяет выделить основные принципиальные недостатки, которые в той или иной мере снижают практическую ценность полученных результатов и сужают область эффективного применения того или иного подхода. Среди них можно выделить: 1. Моделирование производится для плоской постановки задачи. Реальные процессы протекают в горном массиве, который имеет три измерения, а зна- чит, горные породы находятся не в плосконапряженном, а в плоскодеформи- рованном состоянии. 2. Горный массив моделируют сплошной средой. Реальный горный массив представляет собой совокупность литологических разностей, которые нахо- дятся в постоянном взаимодействии между собой. 3. Горные породы представляются изотропным материалом с усреднен- ными физико-механическими свойствами. Реальные горные породы в общем случае проявляют анизотропные свойства, но с большой точностью могут быть описаны моделью трансверсально-изотропной среды, которая характе- ризуется пятью независимыми упругими константами. 4. Отсутствие учета влияния крепи на распределение НДС вокруг вырабо- ток. Реальная крепь обладает определенной жесткостью и оказывает влияние на распределение зон повышенного и пониженного горного давления. 5. Моделирование производится для ограниченной части горного массива, который непосредственно примыкает к моделируемому объекту (горной вы- работке). В реальных условиях в формировании НДС участвует вся надуголь- ная толща горных пород, которая вмещает отрабатываемый пласт. В сложившихся экономических условиях функционирования, горнодобы- вающие предприятия предъявляют все более и более жесткие требования к точности выполняемых расчетов параметров НДС массива необходимых для применения эффективных способов управления состоянием массива горных пород, что является неотъемлемой частью использования высокоэффектив- ных технологий ведения горных работ и современной высокопродуктивной горной техники. На современном этапе развития горной науки и компьютерных техноло- гий наибольший интерес представляют теоретические исследования с приме- нением математических моделей геомеханических процессов, протекающих в массиве горных пород. Цель статьи. Обосновать параметры трехмерной модели слоистого трансверсально-изотропного массива горных пород для изучения закономер- ностей распределения напряженно-деформированного состояния вокруг очи- стной выработки. Основная часть. Обоснование геометрических параметров модели. Геометрические параметры модели зависят, прежде всего, от размеров зо- ны влияния очистных работ на массив горных пород. Для правильного опре- деления этих параметров необходимо использовать данные многолетних маркшейдерских наблюдений и рекомендации нормативно-правовых доку- ментов, которые регламентируют порядок определения этой зоны. В настоя- "Геотехническая механика" щее время в Украине юридическую силу имеют "Правила охраны сооруже- ний и природных объектов от вредного влияния подземных горных разрабо- ток на угольных месторождениях" [1], согласно которым границы влияния очистной выработки на земную поверхность (мульды сдвижения) определя- ются граничными углами. Различают граничные углы по простиранию δ0, по падению β0 и по восстанию γ0 пласта. Значения граничных углов условий За- падного Донбасса согласно [1] принимаются δ0 = γ0 = β0 = 65°. Руководствуясь [1] можно определить геометрические размеры модели в плане. Вертикальные размеры модели принимаются из условия моделирова- ния всей толщи горных пород, которые непосредственно участвуют в форми- ровании НДС массива вокруг очистной выработки, т.е. всех литологических разностей от пласта до земной поверхности и не менее пятидесятикратной мощности вынимаемого слоя вглубь подстилающего угольный пласт массива. Согласно [1], зона влияния очистной выработки на породы почвы пласта рас- пространяется не более чем на 50 вынимаемых мощностей. На рис. 1 и 2 представлены расчетные схемы для определения геометриче- ских параметров модели в плане. Рис. 1 – Главное сечение мульды сдвижения по падению пласта Рис. 2 – Главное сечение мульды сдвижения по простиранию пласта Выпуск № 82 Исходя из расчетной схемы, представленной на рис 1, была получена формула для определения горизонтального размера модели по направлению вкрест подвигания очистного забоя: , (1) где Lм – горизонтальный размер модели по направлению вкрест подвигания очистного забоя, м; Lл – длина лавы, м; γ0, β0 – граничные углы сдвижения по падению пласта, град; α – угол падения пласта, град; Н – расстояние от сере- дины лавы до земной поверхности, м. Исходя из расчетной схемы, представленной на рис. 2, была получена формула для определения горизонтального размера модели по направлению подвигания очистного забоя: , (2) где Sм – горизонтальный размер модели по направлению подвигания очистно- го забоя, м; Sл – подвигание очистного забоя, м; δ0 – граничный угол сдвиже- ния по простиранию пласта, град; Н – расстояние от середины лавы до зем- ной поверхности, м. Вертикальные размеры модели массива горных пород должны удовлетво- рять выражению Z ≥ H + 50m (3) где m – вынимаемая мощность пласта, м. Для условий Западного Донбасса с учетом того, что δ0 = γ0 = β0 = 65° и принимая α ≈ 0°, формулы 1 и 2 могут быть представлены в виде: (4) (5) Диапазон глубин пластов принятых к разработке в ОАО "Павлоградуголь" колеблется в пределах 230 – 700 м, длины механизированных лав изменяются от 180 до 220 м. Физико-механические характеристики горных пород вмещающих угольный пласт. Горные породы имеют весьма сложные физико-механические свойства, зави- сящие от их вещественного и гранулометрического состава, условий залегания и происходивших во время их формирования геологических процессов. В общем случае физико-механические характеристики по различным на- правлениям в массиве пород различны, т.е. горные породы проявляют анизо- "Геотехническая механика" тропию свойств, и в связи с этим изменяется характер распределения напря- жений вокруг выработок. Слоистость, напластование, перемежаемость пород с различными механи- ческими свойствами могут обуславливать механическую анизотропность массива, несмотря на то, что каждая из пород, слагающих толщу, является по данным испытаний образцов механически изотропной. Эта анизотропность связана с существенным различием механических свойств на контактах меж- ду слоями и пластами слагающих толщу пород [2]. В настоящее время разработана теория упругости анизотропной среды, позволяющая вычислять напряжения вокруг отверстий различной конфигура- ции. Известен ряд работ о распределении напряжений вокруг отверстий кру- говой, эллиптической, прямоугольной и других форм при различных случаях анизотропии. Однако лишь отдельные теоретические решения доведены не- посредственно до числа с учетом реальных свойств деформирующейся среды, в частности для горных пород. Важное практическое значение имеет частный случай анизотропии, когда массив горных пород может быть представлен упругой средой, в каждой точке которой имеется одно главное направление и бесконечное множество главных направлений в плоскости, нормальной к первому. Подобный тип анизотропии на- зывают трансверсально-изотропным. Он характерен для тонкослоистых осадоч- ных, метаморфических пород, слюд, пегматитовых тел и др. При этом деформи- рование может быть охарактеризовано с помощью пяти независимых констант, определяющих связь между напряжениями и деформациями [3]. Учитывая то, что для большинства горных пород наблюдается практиче- ски линейная зависимость между приращениями деформаций и приращения- ми напряжений, уравнения обобщенного закона Гука для упругого трансвер- сально-изотропного массива можно записать в виде: Здесь Е1 и Е2 – модули упругости в плоскости слоев и в направлении, нор- мальном к ним. Коэффициент поперечных деформаций μ1 характеризует от- ношение продольных и поперечных деформаций в плоскости слоев, а μ2 – от- ношение соответствующих деформаций в плоскости, нормальной к напласто- ванию. Модули сдвига G1 и G2 являются коэффициентами пропорционально- сти между касательными напряжениями τ и сдвиговыми деформациями γ со- ответственно в плоскостях слоев и в плоскости, нормальной к ним. При этом, модуль G1 может быть выражен через Е1 и μ1 с помощью известного соотно- шения: , (6) Выпуск № 82 Модуль G2 для анизотропного тела является величиной независимой и не может быть точно определен через другие постоянные. Вместе с тем для мо- дуля G2 в работе [4] была получена приближенная формула, которая позволят с достаточной для практических расчетов точностью определить модуль G2 для подавляющего большинства горных пород: , (7) Следует заметить, что на возможную связь упругих постоянных в виде (7) для некоторых ортотропных тел указывал в 1863 г. Б. Сен-Венан, в 1935 г. еѐ использовал в расчетах по оценке напряженно-деформированного состояния скальных пород под сооружениями К. Вольф, в 1954 г. – К.В. Руппенейт для оценки влияния анизотропии пород на концентрацию напряжений на контуре эллиптического отверстия. Расчеты по формуле (7) показали хорошую схо- димость экспериментальных и расчетных значений для G2 [4]. Коэффициент анизотропии Ка для пород Западного Донбасса по данным В.Т. Глушко [5] изменяется в пределах 1 – 5 (табл. 1). Таблица 1 – Коэффициент анизотропии для пород Западного Донбасса Значение Ка Количество проб, % для каждого литологического типа Песчаник Алевролит Аргиллит Ка < 1,2 78,4 13,55 6,20 1 < Ка < 1,5 7,2 44,50 18,75 1,5 < Ка < 2,0 7,2 24,45 17,90 Ка > 2,0 7,2 17,65 53,55 Коэффициент Ка связан с модулями упругости выражением [2] . Для усредненных условий Западного Донбасса все физико-механические свойства горных пород необходимые для построения математической модели уп- ругого трансверсально-изотропного горного массива были сведены в табл. 2. Параметры материала моделирующего крепь очистных и подготови- тельных выработок. Для получения реальной картины распределения НДС массива вокруг очистной выработки необходимо учитывать влияние работы крепи как очист- ных, так и подготовительных выработок с вмещающими горными породами. Однако соотношение геометрических параметров крепи и рассматриваемой области горного массива (см. формулы 1-5) не позволяет моделировать крепь в состоянии «как есть». Поэтому влияние реальной крепи, состоящей из сово- купности отдельных модулей (секции механизированной крепи для очистного забоя и крепежные рамы для подготовительных выработок), в модели пред- "Геотехническая механика" ложено заменить распределенной нагрузкой, которая прилагается непосред- ственно в месте контакта крепи с массивом горных пород. Таблица 2 – Основные физико-механические свойства горных пород Западного Донбасса* Е1 Е2 μ1 μ2 G1 G2 σсж 1 σр 1 σсж 2 σр 2 ρ ГПа ГПа ГПа ГПа МПа МПа МПа МПа кг/м 3 Уголь 39,00 32,00 0,27 0,21 16,18 14,34 27,00 6,40 35,00 4,00 1320 Алевролит 19,70 10,00 0,25 0,14 8,62 5,58 21,00 4,20 28,00 2,00 2510 Аргиллит 20,00 8,00 0,21 0,09 9,16 5,05 17,00 3,20 21,30 1,50 2370 Сланец песчанистый 25,20 16,80 0,20 0,15 10,99 8,57 20,00 4,30 27,40 2,80 2540 Сланец глинистый 11,60 7,10 0,27 0,18 4,90 3,59 14,00 2,40 17,00 1,20 2290 Песчаник 25,60 20,50 0,26 0,23 10,39 9,05 23,00 4,50 32,50 3,00 2600 *Индекс 1 соответствует плоскости изотропии. Индекс 2 соответствует нормали к плоскости изотропии. Моделирование крепи подготовительной выработки. В качестве крепи подготовительных выработок на шахтах Украины нашли широкое применение арочные податливые крепи, изготовленные из стали го- рячекатаной профильной типа СВП. Для условий Западного Донбасса было научно обосновано применение разновидности трехзвенной арочной подат- ливой крепи, которая получила название КШПУ [6]. Сопротивление одной рамы крепи КШПУ из СВП-27 не превышает 400 кН, при этом максимальная возможная конструктивная податливость составляет 300 мм. В узле податли- вости используются силы трения между металлическими элементами, стяну- тыми зажимными накладками, поэтому крепь оказывает постоянное сопро- тивление опускающимся породам кровли [7]. С достаточной для практиче- ских расчетов точностью величина давления крепи на кровлю может быть найдена из следующего выражения: где Р – сопротивление одной арки крепи, кН; а – шаг установки крепи, рам/м; b – ширина выработки, м. Для крепи КШПУ из СВП-27 F≈100 кПа. Моделирование крепи очистной выработки. В настоящее время на шахтах Украины добыча угля производится механизи- рованными лавами, которые оснащены забойными комплексами. В своем боль- шинстве механизированные комплексы представлены следующими типами КД- 80, МКД-90, МКД-90Т, КДМ и др. Основной составной частью каждого ком- плекса является гидрофицированная крепь. К силовым элементам крепи относят- Выпуск № 82 ся гидравлические стойки, которых каждая секция мехкрепи насчитывает две или четыре в зависимости от конструкции. Основной режим работы гидравлической стойки – это режим постоянного сопротивления. В номинальном режиме работы каждая секция механизированной крепи КД-80 постоянно оказывает давление на кровлю порядка 2,5 МПа, на почву – 3,0 МПа. Давление на вмещающие пласт по- роды передается через жесткую плиту основания и перекрытия, которая имеет размеры для крепи КД-80 1,3 м×1,5 м (площадь около 2 м 2 ) и расположена на рас- стоянии 2,5 м от плоскости забоя. При построении модели для корректного учета работы механизированной крепи необходимо задавать вышеуказанные силовые и геометрические пара- метры. Моделирование трещиноватости горного массива. Хотя в земной коре встречаются массивы практически нетрещиноватых по- род, однако, как правило, скальные и полускальные породы разбиты тремя или более системами трещин. На глубине эти трещины могут быть плотно сомкнуты, и в их плоскостях может сохраняться значительное сцепление – до 20-30 % от сцепления в сплошном теле. Трещины могут быть раскрытые, заполненные ми- лонитом или зияющие. В зонах тектонического дробления густота трещин может возрастать вплоть до щебнеобразного состояния породы. Трещиноватый массив можно рассматривать как некогда монолитный и подвергнутый запредельному деформированию до определенной степени, предельный график которого нахо- дится где-то в промежутке между предельным графиком монолита и графиком остаточной прочности. В зависимости от степени трещиноватости сцепление массива будет составлять некоторую долю λ от сцепления образца. Показатель λ называют коэффициентом структурного ослабления. В табл. 2 приведены ориентировочные значения коэффициента структур- ного ослабления для различных пород [8]. Таблица 2 – Значение коэффициента структурного ослабления Характеристика породы Коэффициент структурного ослабления при прочности пород на сжатие в образце (МПа) < 2 2 – 10 10 – 40 > 40 Без четко видимой трещиноватости 0,90 0,70 0,60 0,50 Плотного сложения с нормальносе- кущей трещиноватостью 0,50 0,40 0,35 0,30 Плотного сложения с прерывистыми кососекущими трещинами 0,30 0,25 0,20 0,15 С нарушенной структурой 0,10 0,08 0,06 0,03 Малотрещиноватые массивы с высокими значениями коэффициента структурного ослабления (λ > 0,25) сохраняют основные свойства монолит- ных образцов, в частности хрупкость, разупрочнение в запредельной области деформирования, склонность к динамическому разрушению [8]. Системы основных трещин характерных для рассматриваемой области масси- ва горных пород могут быть смоделированы секущими плоскостями, которые пе- "Геотехническая механика" ресекают моделируемые горные породы в зависимости от ориентации трещин в натуре. Причем на контакте поверхностей задаются условия взаимодействия по- род. Как правило, такое взаимодействие учитывается коэффициентом трения по- коя для пары "порода-порода", но в зависимости от реальных условий могут быть рассмотрены различные варианты сопряжений с определенной контактной проч- ностью. По такой же схеме учитываются условия взаимодействия на контакте слоев слагающих моделируемый массив горных пород. Коэффициенты трения для пары "порода-порода" детально изучены Л.И. Бароном и представлены в [9]. На основании [9] для горных пород ха- рактерных в условиях Западного Донбасса была составлена табл. 3 Таблица 3 – Коэффициенты трения покоя для пары "порода-порода" Горная порода fп Горная порода fп Каменный уголь 0,40 – 0,45 Сланец песчанистый 0,38 – 0,41 Алевролит 0,35 – 0,39 Сланец глинистый 0,36 – 0,40 Аргиллит 0,33 – 0,37 Песчаник 0,38 – 0,42 Как видно из табл. 3, коэффициент трения покоя для всех типов пород с достаточной для практических расчетов точностью может быть принят fп = 0,35÷0,40. Выводы. Трехмерная модель слоистого трансверсально-изотропного мас- сива горных пород, учитывающая влияние силовых характеристик крепи, по- зволит уточнить закономерности распределения напряженно-дефор- мированного состояния массива вокруг очистной выработки. Обоснованные параметры моделирования позволяют наиболее адекватно отобразить процес- сы, происходящие в горном массиве, что позволит уточнить параметры зон опорного горного давления и зон разгрузки. Данные полученные в результате моделирования могут быть использованы для выбора и обоснования эффек- тивного способа управления состоянием массива горных пород вокруг очист- ной и подготовительных выработок. СПИСОК ЛИТЕРАТУРЫ 1. Правила охраны сооружений и природных объектов от вредного влияния подземных горных разработок на угольных месторождениях /Министерство угольной промышленности СССР. – М., Недра, 1981.- 288 с. 2. Лехницкий С.Г. Теория упругости анизотропного тела. М.: Наука,1977, - 415 с. 3. Турчанинов И.А., Иофис М.А., Каспарьян Э.В. Основы механики горных пород. Л., "Недра", 1977. 503 с. 4. Батугин С.А. Анизотропия массива горных пород. Новосибирск, Наука, Сиб. отделение, - 1988, 82 с. 5. Глушко В.Т. Проявление горного давления в глубоких шахтах. – К.: "Наукова думка", 1971. – 196 с. 6. Инструкция по поддержанию горных выработок на шахтах Западного Донбасса. – С.- Петербург – Павлоград: ВНИМИ, 1994. 7. Якоби О. Практика управления горным давлением. – М.: Недра, 1987. – 566 с. 8. Фадеев А.Б. Метод конечных элементов в геомеханике. М.: Недра, 1987. – 221 с. 9. Барон Л.И. Характеристики трения горных пород. – М.: "Наука", 1967. – 208 с. Рекомендовано до публікації д.т.н. С.Ф. Власовим 30.06.09