Размеры и конфигурация зоны запредельных деформаций вокруг выработки Текст научной статьи по специальности «Энергетика и рациональное природопользование»

д.т.н., проф. Литвинский Г. Г. (ДонГТУ, г. Алчевск, ЛНР, ligag@ya.ru)

РАЗМЕРЫ И КОНФИГУРАЦИЯ ЗОНЫ ЗАПРЕДЕЛЬНЫХ ДЕФОРМАЦИЙ

ВОКРУГ ВЫРАБОТКИ

Работа направлена на повышение достоверности и дальнейшее развитие методов прогноза проявлений горного давления в подземных выработках при образовании зон запредельных деформаций (ЗЗД). Впервые получено порождающее решение краевой задачи о равновесии разрушенных пород вокруг горной выработки с учетом возникающей в ЗЗД радиальной прочностной неоднородности пород. Показано, что прочность разрушенных пород в ЗЗД линейно зависит от радиальной (минимальной) компоненты тензора напряжений. Методом малого параметра решен ряд неосесимметричных задач горной геомеханики, установлены закономерности влияния на конфигурацию и размеры ЗЗД основных горнотехнических факторов. Результаты будут полезны при разработке новых методов прогноза горного давления и пересмотре нормативных материалов по обеспечению устойчивости выработок.

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

УДК 622.831.2

1 Актуальность проблемы и цель исследований

Расчет проявлений горного давления (ПГД) в подготовительных горных выработках является неотъемлемым атрибутом процесса проектирования для обеспечения безопасной и надежной их эксплуатации. Методы расчета ПГД претерпели долгую историю создания, развития и совершенствования, что отражено в многочисленных исследованиях и обширных обзорах отечественных и зарубежных ученых [1-5 и др.].

В настоящее время общепринято считать, что в условиях больших глубин, на которые повсеместно переходят рудники и шахты, из всех возможных многообразных форм ПГД (вывалы, своды естественного равновесия, продольно-поперечный изгиб и выпучивание слоев и др.) начинает преобладать разрушение пород в виде формирования зон запредельного деформирования (ЗЗД) пород вокруг всей выработки или на части ее породного контура.

ЗЗД вокруг выработки возникает, растет и принимает окончательную форму под влиянием многочисленных на нее влияющих неосесимметричных факторов. К та-

ким факторам в первую очередь следует отнести:

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

- анизотропию и неоднородность прочностных свойств пород вокруг выработки;

- массовые силы, действующие на породы внутри ЗЗД;

- отличие формы породного контура выработки от кругового очертания;

- неравномерность отпора (реакции) крепи на различных участках породного контура выработки.

При этом в подавляющем большинстве существующих методов расчета ПГД по умолчанию принималось, что влияющие на ЗЗД факторы являются осесимметрич-ными или несущественными. Однако практика ведения горных работ и многочисленные шахтные инструментальные наблюдения показали, то форма ЗЗД претерпевает существенные искажения и чаще всего заметно отличается от круговой формы. А форма ЗЗД предопределяет амплитуду и изменчивость смещений пород

Науки о земле

в выработку, распределение и интенсивность нагрузки на крепь.

Следовательно, для повышения достоверности и надежности расчетных методов прогнозирования ПГД необходимо установить размеры и форму ЗЗД и оценить влияние самых важных горно-геологических и горнотехнических факторов. Достижению этой цели и посвящена данная работа.

Особенностью проведенного исследования является изучение изменения формы ЗЗД для наиболее общего случая разрушения пород вокруг выработки на основе «порождающего» решения краевой задачи развития ПГД. Под порождающим решением мы понимаем исходное решение осе-симметричной краевой задачи об образовании ЗЗД с учетом радиальной прочностной неоднородности внутри зоны разрушения.

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

Объект исследования — зона разрушения и трещинообразования пород массива вокруг выработки, т. е. ЗЗД.

Предмет исследования — закономерности изменения формы ЗЗД при учете воздействия наиболее существенных не-осесимметричных факторов.

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

Основные задачи исследования:

- для получения порождающего решения дать обоснование расчетной схемы и получить решение осесимметричной краевой задачи об образовании ЗЗД вокруг выработки с учетом радиальной прочностной неоднородности пород в этой зоне;

- выполнить постановку и решение неосе-симметричного класса задач методом малого параметра для учета важнейших, ранее неучтенных неосесимметричных факторов;

- произвести анализ и оценить степень влияния каждого из неосесимметричных факторов на форму ЗЗД.

2 Обоснование расчетной схемы задачи

Попытки разработать методы расчета параметров ЗЗД делались многими зарубежными и отечественными учеными начиная с середины прошлого века (Р. Феннер, А. Лабасс, Л. Уокер, Г. Шпа-келлер, Ю. М. Либерман, Ю. З. Заславский и др.), подробный анализ их исследований можно найти во многих обзорах [1-5 и др.].

На ранних стадиях развития горной промышленности, на сравнительно малых глубинах Н ведения горных работ, разрушение пород в окрестности горных выработок происходило в основном в виде образования сводов естественного равновесия [1].

С ростом глубины исходное напряжённое состояние массива (рис. 1, а) претерпевает серьезные изменения: вертикальная компонента, равная весу вышележащей толщи пород ог = yH , возрастает пропорционально глубине Н, а боковые (горизонтальные) напряжения ох = оу составляют

лишь часть Л от вертикальных, с глубиной приближаясь к ог.

Соотношение между горизонтальными ах и вертикальными ог напряжениями

в нетронутом горном массиве было названо коэффициентом бокового распора Л

Л = ох/ог < 1. (1)

Своды естественного равновесия (СЕР) возникали в условиях «малых глубин», когда коэффициент бокового распора отвечал неравенству Л<1. Это вызывало появление в кровле и (или) почве выработки растягивающих напряжений, и происходило разрушение пород в виде трещин разрыва (см. рис. 1, б).

Главная особенность «больших» глубин разработки, где зачастую напряженное состояние близко к гидростатическому Л « 1, состоит в том, что разрушение пород вокруг выработок происходит от напряже-

Науки о земле

нии сжатия, приводящих к развитию трещин сдвига. Последние появляются вначале вблизи породного контура выработки, а затем распространяются далее в окру-жающии массив и образуют зоны разрушения и запредельных деформации (ЗЗД).

Большинство ранее предложенных геомеханических моделеи разрушения и деформирования пород вокруг горных выработок базировались примерно на одинаковых предпосылках и расчетных схемах (рис. 1, в). А именно, в нетронутом горном массиве деиствует гидростатическое исходное поле напряжении от веса вышележащих пород

О = Оу = О = Ун, Л~ 1 (2)

где у — плотность горных пород, кг/м3;

Н — глубина расположения горноИ выработки, м.

Вокруг круглоИ выработки при Л«1 возникает осесимметричное плоско деформированное упругое поле напряжении в виде радиальноИ ог и окружноИ (тангенциальной) Од компонент (см. рис. 1, в).

Расчетные формулы, полученные решением краевых осесимметричных задач при разных математических моделях массива пород, систематизированы в таблице 1.

Если постулировать, что порода в ЗЗД после разрушения становится близкой по свойствам к сыпучей среде, то ее прочность равна нулю <с (r) = 0 , и получаем

решение Р. Феннера. При допущении, что прочность внутри ЗЗД постоянна <с (r) = Const, получаем решение

А. Лабасса, а пластичная среда была исследована в известной задаче А. Галина (табл. 1).

Анализ всех этих моделей показывает, что им присущ общий существенный недостаток — в них предполагается, что внутри ЗЗД породы однородны, т. е. после разрушения их прочностные свойства постоянны и не изменяются вдоль радиуса по мере удаления от контура выработки до границы ЗЗД. Такое несоответствие расчетных схем реальному поведению пород в ЗЗД приводит к необъяснимо большому разбросу расчетных радиусов зон разрушения и смещений пород для различных принятых ранее моделей их поведения.

На самом деле внутри ЗЗД прочность разрушенных пород весьма изменчива: она минимальна на контуре выработки и постепенно возрастает по мере приближения к внешней границе ЗЗД.

У//////////////////////////////////////////^

a)

z ▲ У

Н

Ox

Рисунок 1 — Расчетные схемы ПГД вокруг выработки: а) напряжения в массиве горных пород; б) свод естественного равновесия СЕР; в) зона запредельных деформаций ЗЗД

Науки о земле

Таблица 1

Параметры ЗЗД вокруг горной выработки для разных моделей горных пород

Модель среды в ЗНД Напряжения в ЗНД а < r < r* Радиус ЗНД r* и смещения пород и в ЗНД

1. Сыпучая ^ = ов - А°г = 0 (Р. Феннер) f У-1 " =91 a J "в = ^ [ ^ J r* = а r u = — 2G 2 p ( A +1) q _ (p - "*)(" 1 (A-1) xA+1 *)

2. Идеально-пластичная ^ = Ов -Ог - 2к = 0 r "r = q + k ln — а r "в = q + k ln—+ 2k а p - k - q r* = а exp- 2 - k * 2 p - "r r* U =~----- 2G r

3. Идеально-связная ^ = Ов - Аог -Осж = 0 (А. Лабасс) "r =(q+c^) -c "в = A (q+c ^ ) - c 1 2 ( p + c ) (A-1) r* = а |_( A +1)( q + c )J * / \A+1 u = P20* Г(rJ

4. Связная, после разрушения — сыпучая ^ = Ов - Аог = 0 (проф. Ю. М. Либерман) f JA-1 " =q l aJ / \A-1 "в = A ( ^ ) r* = а u= 2 p - "сж ( A +1) q ?-"* r f r* 2G r l r 1 (A-1) j A+1

Примечание: а — радиус выработки, р = уН — давление пород; q — отпор крепи; А = (1 + Sinp) / (1 - SmpУ; с, G — сцепление и модуль сдвига пород; г*, аг — радиус границы ЗЗД и радиальное напряжение на ней; и — смещения пород на контуре выработки.

Поэтому целесообразно получить в качестве исходной модели поведения пород в ЗЗД новое решение, где этот принципиальный недостаток будет устранен. Только тогда можно принять это решение в качестве исходного — порождающего. Оно может служить основой для последовательного отказа от принятых в ней идеали-заций, главной из которых следует считать осесимметричность. Остальные расчетные схемы, где не учитывается радиальная неоднородность прочности частично разрушенных пород в ЗЗД, будут охватываться предложенным порождающим решением как частные случаи.

3 Получение порождающего решения задачи об образовании ЗЗД

Пусть горная выработка круглой формы радиуса г0 расположена на глубине Н. Исходное напряженное состояние массива является гидростатическим, тогда все компоненты тензора напряжений равны давлению р вышележащей толщи пород. Так, для произвольно расположенной в массиве декартовой системы координат х, у, z для напряжений будет соблюдаться равенство (рис. 1, б)

Ох = Оу = О = P, (3)

Науки о земле

где р — давление вышележащей толщи горных пород p = уН, Па ;

у — усредненный по всей глубине Н объемный вес пород, у « 25 кН/м3.

На достаточно удаленном расстоянии от забоя протяженной выработки (не менее 3...4 ее радиусов Го) массив находится в плоско деформированном состоянии. Если направить ось z по вертикали, ось у по горизонтали и ось х вдоль оси выработки, то в окрестности выработки деформации пород удовлетворяют равенствам sx = Const,

sy , s7 = Var .

У "

Введем полярную систему координат r,в с началом r = 0 в центре выработки (см. рис. 1, в). Задача состоит в определении напряженно-деформированного состояния пород при образовании зоны запредельных деформаций вокруг выработки, на контуре которой установлена податливая крепь с постоянным отпором, равным q = Const.

До образования ЗЗД вокруг горной выработки действует упругое поле напряжений по известному из теории упругости закону (в полярной системе координат r, в ) для области r0 < r <<х

< = p+(p - q)f—1 ; = 0; (4)

где ог , Од и тгд — радиальное, тангенциальное (окружное) и касательное напряжения вокруг выработки.

При осевой симметрии нормальные напряжения являются главными, причем тангенциальное Од = 01 является максимальным, а радиальное ог = 03 — минимальным.

Если прочность горных пород ос на

контуре горной выработки окажется меньше действующих напряжений (4), возникает ЗЗД, где породы полностью или частично разрушены, т. е. возникает прочностная неоднородность.

Для решения задачи о распределении напряжений ог, Од в кольцевой ЗЗД вокруг круглой горной выработки радиуса г0, в которой установлена крепь постоянного сопротивления q, необходимо решить систему из двух уравнений:

а) уравнение равновесия в полярной системе координат г, в при осевой симметрии

dor or -Oq

= 0;

lrQ

= 0; (ro < r < rz)

(5)

где rz — радиус ЗЗД вокруг выработки,

подлежащий определению;

б)уравнение состояния (прочности) разрушенной горной породы внутри ЗЗД при (r0 < r < rz ), которое принимаем в виде линейного паспорта прочности

ав = kpar + о-с, (6)

где ос — показатель прочности пород в ЗЗД на одноосное сжатие с учетом масштабного фактора;

kp — коэффициент угла внутреннего

трения пород, равный

kp =(1 + Sin р) / (1 - Sin р) ,0 < kp < 3;

где р — угол внутреннего трения пород, град., 0 < р < 30°.

На уравнении (6), несмотря на его кажущуюся простоту, следует остановиться подробнее ввиду его высокой важности. В отличие от всех ранее выполненных постановок подобных задач (см. табл. 1), когда прочность пород внутри ЗЗД всегда считали постоянной ос (r) = Const, мы отказываемся от такого упрощения.

На основании изучения поведения пород при запредельном деформировании, пришло понимание, что на самом деле прочность породы в ЗЗД ос изменяется вдоль

радиуса от минимальной на контуре горной выработки до исходной (ненарушенной) на границе между ЗЗД и упругой зоной.

r

r

Науки о земле

Вся сложность математической модели о развитии проявлений горного давления вокруг горной выработки состоит в правильном отражении этой важнейшей закономерности.

Теоретическими и экспериментальными исследованиями запредельного деформирования и разрушения материалов и горных пород при их объемном нагружении в стабилометрах доказано, что изменение их прочности на ниспадающей ветви диаграммы полностью контролируется боковым подпором, т. е. минимальной компонентой объемного нагружения [6, 7 и др.].

Такая закономерность хорошо прослеживается на рисунке 2, где приведены результаты испытаний образцов угля при их многокомпонентном нагружении в стабиломет-ре. Графики запредельного поведения угля в координатах «деформация — напряжение» - о*! при разных значениях бокового давления О3 показывают почти прямо пропорциональную зависимость о*1 от уровня минимального напряжения 03 (рис. 2).

Поэтому вполне обоснованно можно полагать, что остаточная прочность пород ос (г) вдоль радиуса г в ЗЗД пропорциональна минимальной компоненте напряженного состояния, т. е. радиальному напряжению ог . В первом, достаточно приемлемом приближении вполне допустимо принять эту зависимость линейной

(7)

"c (r) = kc"r

"в = kp"r + "c (r) = kp"r + kc"r = = (kp + kc )"r .

Рисунок 2 Диаграмма испытаний угля при трехосном сжатии (по данным ВНИМИ [6])

После подстановки (8) в исходное дифференциальное уравнение равновесия (5) и его интегрирования с учетом граничного условия (при г = г0 ог = д) получим значения напряжений внутри ЗЗД

"r = q

( r JM f r V

, "в = (m +1) q

l r0 )

l r0 )

(9)

he

= 0,

(r0 < r < rz ),

где кс — коэффициент сопротивления пород горному давлению, зависящий от условий разрушения пород внутри ЗЗД и определяемый из граничных условий развития области разрушения пород вокруг горной выработки.

Подставляя (7) в (6), получим искомое уравнение условия прочности пород в ЗЗД

где ц = кр + кс -1 — коэффициент, учитывающий сопротивление разрушению и суммирующий влияние коэффициентов кр,

кс, он показывает степень сопротивления породы развитию ЗЗД вокруг выработки.

Чтобы найти граничный радиус ЗЗД г2 , воспользуемся условием равенства радиальных ог и тангенциальных Ов напряжений по обе стороны границы ЗЗД г = г2,

которые следует определить по уравнениям (4) и (9). После несложных преобразований получим искомую зависимость

rz = r0

2 p

(m + 2) q

1/ m

=r

2 p -"c

(kp+1) q

1/M

(10)

Для определения значения коэффициента сопротивления кс воспользуемся до-

Науки о земле

полнительным условием, требующим, чтобы на границе с упругой областью прочность пород была равна исходному своему значению в массиве, т. е. ос (г2 ) = ос = kcoг (г2 ) . Используя выражение (10), получим для коэффициента сопротивления kc и коэффициента разрушения ц удобные формулы для определения по исходным данным

kc =( kp + 1) U = 2

Р '2p — Стс (kp — 1) P + °c

(11)

2 P — °c

r, Y

о = q

2p 2 p — о c

\ ro;

U + 2

kp +1 (12)

пРи (r = rz ).

ного радиуса выработки Г), получим искомое распределение упругих напряжений за пределами ЗЗД

= P

1 + U

ч2

Г

U + 2 тгв =0; (rz ^r .

(13)

Из (11) можно установить пределы возможного изменения установленных коэффициентов. Так, эти пределы равны: 0 < kc < да для коэффициента сопротивления и kp +1 < ц < да для коэффициента устойчивости где нижний предел получен для сыпучей среды ос = 0, а верхний (при отсутствии ПГД в выработке) для обоих случаев одинаков и достигается при 2 р -ос = 0, т. е. когда ЗЗД не образуется.

Чтобы получить распределение напряжений вокруг выработки за пределами ЗЗД при г > г2, используем известное решение теории упругости о напряжениях в окрестности подкрепленных круглых отверстий.

За контур отверстия мы принимаем границу ЗЗД гг, на которой действует радиальное напряжение огг, а его определим из ранее полученной зависимости (9) при г = г2 . Подставляя (10) в (9), получим для

радиального подпора на границе гг

Это распределение напряжений при возрастании радиуса г ^ да асимптотически стремится к уровню исходного напряженного состояния массива пород ог = Од = р. На границе ЗЗД сохраняется непрерывность радиальных и тангенциальных напряжений и смещений.

Более того, оказывается, что радиальное напряжение Огг на границе ЗЗД не зависит

от размеров выработки и отпора q установленной в ней крепи. А если из уравнения (9) мы найдем такое значение отпора крепи, при котором ЗЗД не образуется, т. е. г2 = г0, то это будет не что иное, как требуемая величина отпора (грузонесущей способности) идеально жесткой крепи, совпадающая с (12)

4 = О, = = Х-О, (г, = г) ,(14)

1

U + 2

kp +1

где — минимальное значение отпора жесткой крепи, устанавливаемой в горной выработке, которое предотвращает возникновение ЗЗД.

Зависимость (14) дает верхнюю теоретическую границу для выбора отпора 4 жесткой крепи в горной выработке, которая прямо зависит только от исходного поля напряжений и от коэффициента устойчивости /. Если порода не обладает прочностью, т. е. является идеально сыпучей средой, тогда kc = 0 и коэффициент устойчивости принимает минимальное значение

Подставляя в (4) вместо отпора крепи 4 значение Огг из (12) и гг вместо началь-

U = kp — 1 =

2Sin p 1 — Sin p

z

Нижнюю границу отпора для податливой крепи можно считать минимальной грузонесущей способностью qm, при возникновении ЗЗД она всегда должна быть больше нуля qm > 0. При образовании ЗЗД горная выработка без крепи существовать не может и будет полностью разрушена, т. е. произойдет «завал» выработки.

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

Важным новым отличием решения является учет прочностной неоднородности при разрушении пород: в пределах ЗЗД прочность частично разрушенных пород (разупрочнения) изменяется прямо пропорционально радиальным напряжениям. Такая предложенная закономерность органически вытекает из поведения пород на ниспадающей ветви диаграммы нагруже-ния, т. е. имеет убедительное экспериментальное обоснование.

Это позволило значительно улучшить учет особенностей реальных процессов разрушения пород в ЗЗД и стать достаточно достоверным базисом для дальнейших построений расчетных методов прогноза проявлений горного давления в выработках.

4 Исследование закономерностей влияния неосесимметричных факторов на ПГД

Основная идея решения неосесимет-ричного класса задач горной геомеханики состояла в следующем. Исследование закономерностей влияния неосесимметрич-ных факторов на форму и размеры зоны неупругих деформаций производилось асимптотическим методом малого параметра путем последовательного отказа от принятых в порождающем решении (7-9) идеализаций расчетной схемы, которые, хотя и являются общепринятыми в подобного рода решениях, заметно снижают достоверность результатов расчетов.

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

агг + г~1(тгв в + аг -ав) = -£гуътв,

*тв,т + Г ~1(ав,в + 2тг0) = -£rУcos0,

где запятой обозначено дифференцирование, а £г — малый параметр при учете

сил гравитации.

Напряжения в пластической области представляли рядами по степеням малого параметра

а1} = 1 £ка¥}, Ш = г ,в) (15)

к=0

и выражали через функцию напряжений

а? > = г-и^1 + г -2и > = и<к>; ¿Цк =-(г\, (16)

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

Решения этих дифференциальных уравнений получали в бесконечных рядах, комплексные коэффициенты которых определяли из граничных условий задачи. Затем задачу сводили к отысканию границы ЗЗД и функций Н. И. Мусхелишвили, голоморфных во всей плоскости вне границы и непрерывных вплоть до нее. Граничные значения этих функций определяли по функции в пластической области.

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

Науки о земле

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

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

5 -л да

7 = о(д) = г*(д + X X 8(к])д1~к), (17)

j=1

к=0

где д — аффикс единичной окружности;

/ — мнимая единица;

£] — малый параметр, соответствующий ]-му фактору;

0] — угол отклонения оси симметрии ]-го

фактора от вертикали;

Зк]) — коэффициент искажения формы

ЗЗД для ]-го фактора.

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

Общее качественное представление о влиянии изученных факторов на форму ЗЗД можно получить по рисунку 3, а количественный анализ показал следующее (см. табл. 2):

1. Негидростатичность исходного поля напряжений была учтена с помощью решения упругопластической задачи по методу, предложенному проф. А. Г. Прото-сеней (1972 г.). ЗЗД принимает форму эллипса, вытянутого по направлению действия минимального главного напряжения.

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

плексный ряд Фурье по окружному направлению. Коэффициенты разложения в ряд Фурье можно получить из полярных диаграмм прочности породы, построенных, например, с помощью портативного прибора для экспресс-испытаний горных пород ППЭИ-2м.

Форма ЗЗД представляет собой цилиндрическую поверхность с сечением, большая ось которой вытянута перпендикулярно напластованию горных пород. Степень «отклика» границы ЗЗД на возмущающее влияние прочностной анизотропии быстро затухает с возрастанием номера гармоники при разложении функции прочности в комплексный ряд Фурье. Поэтому в практических расчетах есть смысл рассматривать гармоники не выше третьего или четвертого порядка.

При номере гармоники к = 1 прочностная анизотропия пород приобретает особый физический смысл прочностной неоднородности, под которой понимается изменение средневзвешенной прочности пород вокруг выработки по полярному углу. Решение позволяет оценить амплитуду конгруэнтного сдвига ЗЗД в сторону пород с пониженной прочностью. Степень влияния прочностной неоднородности весьма ощутима.

Рисунок 3 Изменения конфигурации ЗЗД под влиянием различных неосесимметричных факторов

Науки о земле

3. При значительном развитии ЗЗД вокруг выработки собственный вес среды в пластической области приводит при малых значениях отпора крепи к сдвигу границы ЗЗД вверх по отвесной линии. Количественная оценка этого эффекта методом малого параметра для критической ситуации (обрушение пород) не вполне пригодна. Поэтому получена оценка сверху, позволяющая вычислить эффективный отпор крепи в кровле (-) и в почве (+) выработки

Исследование формы зоны запредельн

радиуса а, закрепленной крепью с конструктивным отпором q

V • а

Я* = q + —л . (18)

¡ —1

Зависимость (18) дает очень важную в практическом отношении оценку минимальной грузонесущей способности крепи, устанавливаемой в выработке.

Таблица 2 х деформаций (ЗЗД) вокруг выработки

Возмущающий фактор и его аналитическое представление Коэффициент формы границы ЗЗД

Разнокомпонентность (негидростатичность) исходного поля напряжений массива горных пород с коэффициентом бокового распора Лф 1 ^ = 1—Л(2+и) 1 + л

Анизотропию и неоднородность прочностных свойств пород вокруг горной выработки учтем через переменную удельную силу сцепления 0{в) = и~А 1[1 + е^ (0)]; §Р =—2 и — А +1 х к и(и + 1)(и + 2) к2 х(Мк + 2ReNk, г—к)

X f (в) = 2 свке1кв —X

Массовые силы в зоне запредельных деформаций учитываем в виде дифференциального уравнения °г,Г + Г ~1(тгв,в + — °в) = —^^пв 5(3) = V а (и +1) (и + 2) ; 1 2и • д*г* (и — 2) (2и — 1)' V а я* = q ± . и — 1

Переменный отпор крепи по контуру горной выработки выразим в виде нормальных и касательных к контуру усилий X X N = 2 /кек ; Т = 2 сркеЛв ; —X —X $к4) = — 2/kRe(Mk, Л—2) ир к

/к = /—к; п = —п к

Некруговая (произвольная) форма контура выработки как функция текущего угла 9

X Го = 1 + е5/(в) = 1 + е5 2 щеЛв ; —X ^ = 4 ¡^^(М^ Г ~Лк) и + 2 к

м>к = ж"1 | /(вуквс1в —ж

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

а) граница ЗЗД ближе к контуру выработки там, где вектор поверхностных усилий больше;

б) чем больше радиус ЗЗД, тем меньше влияние неравномерности отпора;

в) чем выше порядковый номер гармоники неравномерности отпора, тем меньше ее влияние на форму ЗЗД, что позволяет ограничиваться первыми тремя гармониками разложения вектора в ряд Фурье.

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

Исходя из особенности поведения решения, можно учесть возмущение контура выработки гармониками не выше четвертого порядка, чего вполне достаточно для практических расчетов. Область пригодности решения ограничена случаем границы ЗЗД и контура выработки.

Выполняя обобщение результатов исследований в этом направлении, можно заключить, что изученные факторы существенно различаются по своему влиянию на ЗЗД. Наиболее «дальнодействующими», а потому и важнейшими, являются разно-компонентность исходного поля напряжения массива и прочностная анизотропия горных пород.

Неравномерность отпора крепи и некруговой контур выработки заметно сказываются на форме ЗЗД лишь при незначительном ее развитии, и их влияние быстро (по степенному закону) затухает по мере роста радиуса ЗЗД, что позволяет назвать их факторами «ближнего действия». Массовые силы в ЗЗД следует учитывать при незначительном конструктивном отпоре крепи и достаточном развитии ЗЗД.

Полученные 1...5 закономерности позволяют с учетом конкретных горногеологических и горнотехнических условий нахождения горной выработки прогнозировать не только качественные, но и количественные показатели и особенности ПГД в выработках, что обеспечит надежность и долговечность их эксплуатации, повысит безопасность работ и снизит неоправданные издержки на их крепление, поддержание и ремонт, повысит эффективность подземных работ. Важным направлением дальнейших исследований становится пересмотр существующих и разработка новых отраслевых нормативных документов с учетом полученных научных результатов по прогнозированию ПГД в выработках, расположенных на больших глубинах разработки и находящихся в сложных горно-геологических условиях.

Выводы

Проведенные исследования позволили получить следующие новые научные результаты:

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

2. Главным недостатком существующих модельных представлений и принятых идеализаций можно считать принимаемую по умолчанию гипотезу о неизменности прочностных свойств пород внутри ЗЗД и отсутствие учета влияния неосесиммет-ричных факторов на размеры и конфигурацию ЗЗД.

3. На основании изучения поведения пород при запредельном деформировании впервые сформулирована и доказана закономерность о линейной зависимости в ЗЗД прочности от радиальной компоненты тензора напряжений.

4. Впервые решена краевая задача о равновесии ЗЗД вокруг выработки с учетом

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

5. На основе полученного порождающего решения упругопластической задачи впервые сформулирован и решен с применением метода малого параметра ряд не-осесимметричных задач горной геомеханики, установлены закономерности влияния на конфигурацию и размеры ЗЗД основных горнотехнических факторов: не-гидростатичности исходного поля напряжений, прочностной анизотропии и неоднородности горных пород, некруговой форма выработки, сил гравитации пород внутри ЗЗД.

Библиографический список

6. Изучено влияние на конфигурацию ЗЗД различных неосесимметричных факторов, произведена их сравнительная оценка и сопоставление. Показано, что наиболее «дальнодействующими» являются негидростатичность исходного поля напряжений и прочностная неоднородность пород, менее значимо проявляются некруговая форма выработки и вес разрушенных пород в ЗЗД, а влияние отпора крепи тем больше, чем он ближе к критическому своему значению.

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

1. Протодьяконов, М. М. Давление горных пород и рудничное крепление [Текст] / М. М. Протодьяконов. — М. — Л. : ОГИЗ, 1931. — 186 с.

2. Вопросы теории горного давления [Текст] : сб. переводов / под ред. проф. П. И. Городецкого, А. А. Борисова. — М. : Госгортехиздат, 1961. — 300 с.

3. Jovanovic, Р. Proektovanje i proracun podgrade horizontalnih podzemnih prostorija [Text]. T. 2 / Petar Jovanovic. — Beograd: Rudarsko-geoloski fakultet Univerziteta, 1994. — 316 str.

4. Заславский, Ю. З. Проявления горного давления в глубоких шахтах Донбасса [Текст] / Ю. З. Заславский. — М. : Недра, 1963. — 238 с.

5. Tajdus, А. Geomechanika w budovnictwie podziemnym. Projectowanie i budowa tuneli [Text] / Antoni Tajdus, Marek Gala, Krzysztof Tajdus. — Krakow : Akademija AGH, 2012. — 762 s.

6. Ставрогин, А. Н. Механика деформирования и разрушения горных пород [Текст] / А. Н. Ставрогин, А. Г. Протосеня. — М. : Недра, 1992. — 224 с.

7. Литвинский, Г. Г. Аналитическая теория прочности горных пород и массивов [Текст]. — Донецк : Норд-Пресс, 2008. — 207 с.

© Литвинский Г. Г.

Рекомендована к печати д.т.н., проф., зав. каф. РПМ ДонНТУПетренко Ю. А.,

к.т.н., доц. каф. СГ ДонГТИ Шульгиным П. Н.

Статья поступила в редакцию 02.11.2020. Prof. Litvinsky G. G. (DonSTI, Alchevsk, LPR)

DIMENSIONS AND CONFIGURATION OF THE ZONE OF EXTREME DEFORMATIONS AROUND THE WORKING

The article is aimed at increasing the reliability and further development of methods for predicting the manifestations of rock pressure in underground workings during the formation of the extreme deformation zones (EDZ). For the _ first time, a generating solution of boundary problem of the equilibrium destroyed rocks around mine workings is obtained, considering the radial strength inhomogeneity of rocks that occurs in EDZ. It is shown that the strength of destroyed rocks in EDZ linearly depends on

the radial (minimum) component of the stress tensor. By the small parameter method there has been resolved a number of nonaxisymmetric problems of geomechanics, as well the regularities of influence on configuration and size of major EDZ mining factors were determined. The results will be useful in developing new methods for predicting rock pressure and reviewing regulatory materials for ensuring the stability of workings.

Key words: mine working, rock pressure, extreme deformations, calculation methods, rock stability, destroy zone configuration, strength heterogeneity, nonaxisymmetric factors.