Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Расчет крепи капитальных горных выработок

..pdf
Скачиваний:
8
Добавлен:
12.11.2023
Размер:
16.09 Mб
Скачать

Из полученной 10. М. Либерманом зависимости можно составить следующее выражение для нагрузки на крепь:

Р = ( у Я -

(1 - s i n ф) ( - £ ) “

(9.1)

Выражение для смещений поверхности выработки имеет следущий вид (если угол внутреннего трения в зоне разрушения и упругой области одинаков):

 

 

R

 

 

<*сж

л

 

V * -

Асж

 

 

UR =

 

Sill <P ( Y# +

 

2

(9.2)

 

2G

2

•) |^(1 — sin ф)----- -

 

 

 

 

/ L'

"

Р

 

 

где показатель А принимает значения:

 

 

 

 

д __ 1— sin ф

— при условии несжимаемости материала в зоне разру­

 

sin ф

шения;

 

 

 

 

 

 

 

___

 

 

 

 

 

 

А =

— при условии

пластического

потенциала

(ассоции­

 

sill ф

 

рованный закон

течения).

 

 

 

 

 

 

 

 

 

 

 

10. М. Либерман отмечает, что при условии пластического потен­ циала получается сильная зависимость смещения пород от реакции крепи, которая на практике не наблюдается. Поэтому, принимая во внимание равную обоснованность обоих вариантов зависимости (9.2), он рекомендует отдать предпочтение условию несжимаемости.

Решение 10. М. Либермана обобщено М. Т. Алимжановым на слу­ чай, когда порода в зоне разрушения обладает некоторым сцепле­ нием К г <С К, где К — сцепление в нетронутом массиве [4]. В этом случае справедливы выражения:

Р =

[ ( Y# -

“^р") (1 —sin ф

) К г ctgср] ( - £ ) “ _ * ! ctgq>; (9.3)

 

 

 

1

п

/

 

sin ср

a N[(Y tf--^ )(l-sin < P )+ * iC tg q >

« п = - ^ - 81п Ф( у Я + - ^ . ) V

2 PJ+ Ki ctg <р

Смещения получены при условии пластического потенциала. Аналогичная модель исследована Д. В. Хоббсом [220], который

воспользовался эмпирическим условием пластичности в зоне разру­

шения:

 

<*1 = Яо| + а8.

(9.4)

Здесь S , Ъ — эмпирические постоянные (для аргиллита В =

49,5;

Ь = 0,61).

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

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

тангенциальных напряжений (рис. 22), вызванный резким измене­

нием свойств материала при переходе через эту границу [252]. Принципиальная возможность учета непрерывного изменения

свойств материала (сцепления) в зоне разрушения показана в работе В. С. Сажина [153], хотя сам автор рассматривает изменение сце­ пления как следствие взрывных работ при проходке выработки.

В работе А. Г. Протосени [141] решение В. С. Сажина обобщено на случай неравнокомпонентного поля напряжений на бесконечности при условии возникновения пластической области по всему контуру

 

 

отверстия.

Неоднород­

 

 

ность

сцепления

задана

 

 

в виде:

 

 

 

 

 

 

К (г) =

* ( 00) - * ,

( 4 ) "

 

 

 

( К ^

0;

 

 

 

 

п = 1,

2,

6).

 

(9.5)

Рис. 22. Распределение напряжений

н массиве вокруг

TJ>

у

п р

л

выработки при наличии зоны

разрушении

Jn.aK И В

раООТв JD. L .

Ь а -

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

Представляет интерес предложение Е. Б. Дружко [63], который выделил между зоной полного разрушения и упругим массивом промежуточную зону — зону постепенного перехода породы в не­ нарушенное состояние. Условие пластичности в этой зоне имеет вид:

 

 

 

<*6 -

Ог =

д7-Н'е °сж+ 00"

(9-6 )

где R'c — граница между

зоной полного

разрушения,

в которой

 

К =

0, и зоной постепенного перехода.

 

Нагрузка на крепь определяется соотношением

 

Р =

3К cos ср

, 2уЯ

К etg ср

 

(ж)“ ^

а —1

2“I- ct

 

 

 

 

 

'“ - «

(

ж - 1)

 

 

 

 

 

 

 

Некоторые

соображения,

связывающие

изменение

механиче­

ских свойств пород с величиной деформации пород в зоне разруше­ ния, изложены в работе А. П. Максимова и др. [114].

В работе Г. Л. Фисенко [172] высказан ряд критических замеча­ ний по поводу решения 10 . М. Либермана. В предлагаемом решении

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

ния пересекается с огибающей для нетронутого массива. Г. Л. Фисенко отмечает, что в ряде случаев ломаную огибающую можно с не­ большой погрешностью заменить прямолинейной. В этом случае расчетные формулы соответствуют упругопластической однородной модели, в частности — зависимости (8.3), но при этом физический смысл входящих величин несколько иной.

При определении смещений пород предлагается учитывать объем­ ное расширение их в зоне разрушения.

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

Критерии замены упругой или упругопластической модели рас­ сматриваемой моделью соответствуют критериям первого (для хруп­ ких пород) или второго предельного состояния (см. гл. I).

При значительных смещениях пород и развития зоны разрушения создаются условия для явлений вывалообразования, соответству­ ющих жесткопластической модели взаимодействия пород и крепи. На это обстоятельство обратил внимание Ю. М. Либерман, который полагал, что давление на крепь не может быть меньше веса пород в зоне разрушения. На основании указанного условия им получено нелинейное уравнение

(9.8)

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

чивающим применение рассматриваемой модели, является образова­ ние сползающего объема. Следовательно, минимальная нагрузка на крепь определяется формулой (7.19), а максимальное смещение пород — следующим выражением, получение которого не вызывает затруднений:

“ я - - k Ф № + Т ) [ * < Р - ^ Г 1 О - т ) ] “ <»•*»

§ 10. ВЯЗКОУПРУГАЯ МОДЕЛЬ ВЗАИМОДЕЙСТВИЯ ПОРОД И КРЕПИ

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

03

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

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

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

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

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

Теория наследственной ползучести горных пород

Эта теория основывается на математическом аппарате (раздел теории интегральных уравнений), разработанном Л. Больцманом (1878 г.) и развитом в дальнейшем В. Вольтерра (1913—1931 гг.)

* Не следует путать со структурными реологическими моделями, которые представляют собой схематическое изображение взаимного расположения вяз-; ких, упругих л других элементов и являются наглядной иллюстрацией уравне­ ния состояния данной конкретной среды (модели Максвелла, Фойгта, Шведова, Бингама [209, 216]).

U

для описания явлений линейной ползучести. Впоследствии Ю. Н. Работнов показал, что задачу теории линейной наследственной ползучести можно формально рассматривать как задачу теории упругости, в которой вместо упругих постоянных необходимо при­ менить временные интегральные операторы с ядром ползу­ чести [65, 142]. Указанное положение, названное Ю. Н. Работновым принципом Вольтерра, позволяет решать задачи механики горных пород, в которых граничные условия и объемные силы могут при­ ниматься не зависящими от времени.

Уравнение состояния среды в теории линейной наследственной ползучести представляется интегральным уравнением Больцмана— Вольтерра:

_1_

L (t, x)G(x)dx

( 1 0 . 1 )

Е

 

 

где L (£, т) — ядро ползучести (наследственности).

При постоянном напряжении это уравнение можно представить

в виде:

 

«(*) = -§-(1 + £ * ) = -§-.

(Ю.2)

t

где L* = JL (tf, т) dx — интегральный оператор типа Вольтерра;

о

Евременной оператор.

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

Получение интегральных операторов и особенно их функций при произвольном ядре ползучести связано с серьезными математи­ ческими трудностями. Ю. Н. Работновым [142] предложена в каче­ стве ядра ползучести экспоненциальная функция дробного порядка

LU

г-рп (I - т)п(1-а)

 

Эа (—Р, t — T) = (f — т)-“2

(10.3)

Г [(я-М) U “~а)]

п=0

 

Свойство экспоненциальной функции дробного порядка, определя­ емое теоремой умножения,

Э*а (*)3* (у) = э « ---- (* + V) (10.4)

позволяет упростить преобразование интегральных операторов, которые определяются по таблицам дробно-экспоиенциальной функции [143]. Однако для числовых расчетов с использованием ЭВМ таблицами воспользоваться нельзя. В этом случае применяется аппроксимация М. И. Розовского [145]:

Эа (—Р) 1 -jp [1 — ехр {—Putf1'*)],

(10.5)

65

где

^ = (1_а)1-а; р= 8Г (1 — а).

Возможность применения теории наследственной ползучести к ре­ шению задач механики горных пород впервые экспериментально обосновал Ж. С. Ержанов [65]. Он показал, что деформирование ряда горных пород до определенного уровня нагружения соответ­ ствует закону линейной наследственной ползучести с ядром в виде степенной функции (Абелево ядро):

L ( t , т) = 8(* — т)~а,

(10.6)

которое является частным случаем дробно-экспоненциальной функ­ ции (10.3). Используются и другие ядра [39, 196].

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

боту) [1, 32], неровности породного

контура сечения выработки

п различные формы сечения [66, 67],

слоистость пород, выража­

ющаяся в анизотропии свойств массива при различной ориентировке

выработки относительно элементов залегания

[50, 66],

физическая

и геометрическая нелинейность

[50, 68] и т. д.

 

Исследования показали, что

упруговязкая

модель

взаимодей­

ствия пород и крепи близка упругой. Отличие заключается в том, что упруговязкая модель учитывает так называемое «последействие», т. е. деформирование пород во времени после снятия или приложения нагрузки. В задачах линейной ползучести при незакрепленных выработках поля напряжений в упруговязкой и упругой средах точно совпадают. При наличии жесткой крепи напряжения в упруго­ вязкой среде релаксируют во времени. Указанные особенности опре­ деляют условия применимости упруговязких моделей взаимодей­ ствия пород и крепи к реальным объектам.

При решении задач теории ползучести горных пород важное значение имеет отмеченный выше (см. § 6) принцип И. В. Родина,

согласно которому поле напряжений в массиве вокруг выработки делится условно на «природное» и «снимаемое». Ползучесть горных пород проявляется только под действием снимаемого поля напря­ жений.

Функции ползучести и их применение для практических расчетов *

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

сдругой — главным образом в расшифровке операторных выражений

*Обоснование функций ползучести выполнено при участии А. М. Линь-

кова.

66

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

В качестве примера рассмотрим функцию ползучести при исполь­ зовании Абелева ядра (10.6). Из выражения (10.2), подставляя ядро ползучести в интегральный оператор, получим выражение для вре­ менной функции модуля деформации

 

Et

Е

(10.7)

 

1 + Ф ’

где Ф =

Ы1-*

 

 

 

 

1 —а

 

 

Здесь и далее а и 6 — параметры ползучести. Далее, пользуясь

выводом об отсутствии объемного последействия, сделанным

10. Н. Работновым на основании экспериментальных данных

[142],

получим уравнение

 

 

 

 

1 — 2р*

_

1 — 2|LI

const,

(10.8)

Et

Е

 

 

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

^ = 0 .5 — ^ - .

(10.9)

Другие временные операторы, содержащие Е и |л, преобразуются, как обычные алгебраические выражения с использованием исходных временных функций Et и [it. Например, выражение для временной функции модуля сдвига получим, подставив в известную зависимость

Gt

 

Et

(10: 10)

2(1

+14)

 

 

выражения (10.7) и (10.9), тогда

Gt

G

(10.11)

 

Ф

1 +

1,5

 

1 + JLI

 

Для сравнения приведем выражения для временных операторов рассмотренных величин, которые получены с использованием основ­ ных свойств «Э-операторов и аппроксимации (10 .5):

Ё ~ Е е х р [—к>6Г(1 —а)

[I = 0,5 + (fi — 0,5) exp [—w6r (1 —а) ti-“];

(10.12)

67

 

 

 

 

 

 

Т а б л и ц а 18

 

6

/, сут

Е

*1

G

м

 

E t

 

 

M t

 

 

 

 

 

 

 

 

 

 

 

0,6

0,0094

0,01

1,00

1

1,00

1,00

 

 

1

0,89

1

0,87

1,05

 

 

100

0,80

1

0,75

1,00

0,7

0,1

0,01

0,94

1

0,92

 

 

1

1,00

1

0,88

 

 

100

0,83

1

1,25

 

0,01

0,01

1,00

1

1,00

0,90

 

 

1

0,97

1

0,96

1,00

 

 

100

0,90

1

0,90

1,03

 

0,001

0,01

1

1

1

 

 

1

1

1

1

 

 

100

1

1

1

0,8

0,0094

0,01

1

1

1

1

 

 

1

1

1

1

1

 

 

100

1

1

1

1

Втабл. 18 сопоставлены временные функции (10.7), (10.9), (10.11),

атакже более сложная функция

M t

M l + Urf)

(10.13)

Et (1

—М/)

 

 

с точными значениями временных операторов, вычисленными после их расшифровки по таблицам [143]. Сходимость является вполне удовлетворительной, тем более, что аппроксимация (10.5) дает зна­ чительно меньшую точность.

Заметим, кстати, что в работе [65] в формуле временного опера­

тора для М имеются опечатки, поэтому приведем ее в исправленном виде:

2

(*+-!£)[‘+т£5г*41^

(10.14)

 

 

 

Предлагаемый математический аппарат распространяется и на иррациональные временные операторы, используемые при анализе анизотропной среды. Расшифровка иррациональных операторов встречает серьезные затруднения [77], поэтому при решении задач наследственной ползучести с учетом анизотропии пород вводятся иногда недостаточно обоснованные упрощения. Например, пред-

68

полагается, что ползучесть сказывается лишь на модуле деформации материала при сжатии его под углом 45° к плоскости изотропии. Остальные характеристики среды полагаются постоянными, не зависящими от времени. Но даже и в этом случае расшифровка интегральных операторов оказывается не простой. В работах [50, 75] расшифровка наиболее распространенного временного оператора приводится на шести страницах. Этот оператор имеет следующий вид:

где

 

А — 2к -f- AM,

 

(10.15)

 

 

 

 

 

 

 

т =

(r<>

2fi2 (1

 

 

1-У?

 

 

 

 

G2 — временной оператор модуля

сдвига

для плоскостей, нор­

 

мальных к плоскости изотропии:

 

 

 

Go

1

 

1

 

 

 

 

 

Ei

 

 

 

 

E2

 

E i i E о,

p 2 — модули упругости при сжатии в плоскости изотро­

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

коэффициенты поперечной деформации; ЕАЪ— временной оператор модуля деформации при сжатии под углом 45° к плоскости изо­ тропии.

Окончательное выражение иррациональной функции операторов

(10.15) имеет вид:

 

А

и

(10.16)

 

2^-ар

где i/А — некоторая функция упругих констант; х, Р — параметры ползучести.

До преобразования иррациональной функции (10.15) ее можно представить в виде:

 

 

= УК 1/1 +

хЗа (р).

(10.17)

В табл. 19

дано

сравнение расчетных значений функций (10.16)

и (10.17)

при а

= 0,7 и параметрах

ползучести, взятых

из

 

 

 

 

 

Таблица

10

Порода

 

Р

-А- при /, сут

 

 

10

100

 

 

 

 

 

 

Песчаник

 

0,719

0,068

1,010

1,003

 

Песчанистый сланец

5.062

0,071

1,009

1,002

 

69

работы [50]. Сравнение говорит в пользу функции (10.17), к тому же выражение (10.16) является приближенным и справедливо при (3^“а > 1 .

Для упрощения получения временных функций при ядре ползу-

чести

(10 .6) построены номограммы для определения функций пол­

 

 

 

 

 

 

зучести Ф (рис. 23). Функ­

 

 

 

г 0,25

ции

ползучести

можно

так­

0,02 -

 

 

 

же получать непосредственно

 

 

 

 

0 ,0 5 -

 

 

 

-0 ,5

из экспериментов. Например*

0,05 -.

5

 

 

0,1 1

0,001

- /

при

одноосном

сжатии

по­

г

стоянными

усилиями

функ­

0,2'-

-

0,002

 

- 0,005

 

ция ползучести определяется

0 ,5 -

- 0,005

Н

по

формуле

 

 

 

 

1

0,01

 

 

 

 

 

 

 

 

 

1

 

 

 

75

 

 

 

 

 

 

 

 

 

2 4

 

 

 

Ф =

 

 

 

1 ,

(10.18)

 

 

 

-10

п

 

e<i)

3

 

 

 

 

 

b J

 

 

 

£ 12

 

 

 

 

 

 

 

 

 

5

V

 

 

 

 

где

£Q\

 

г ) 1 — продольные

 

 

 

t,cym

деформации — мгновенная и

0,03-

 

 

 

 

 

 

-0 ,5

в момент времени t при на­

0,05-.

 

 

 

 

пряжениях

о t.

 

Функцию

0,1

^

 

Я

 

Ч

ползучести, поданным экспе­

0, 2 -

 

0,001

риментов,

удобно

представ­

0,5 ±

 

г

 

лять

графически

(рис.

24),

 

-

0,002

 

1

-

 

-о,ооь

- 2

тогда

отпадет необходимость

 

L0,01

в аппроксимации эксперимен­

2 1

 

 

5-.

 

 

 

■3

тальной кривой.

 

 

 

 

 

 

 

Рассмотрим

 

несколько

 

 

 

 

 

10 \

 

 

 

5

примеров

использования

1bJ

 

 

 

функций

ползучести.

 

 

 

 

 

 

 

 

Учет

статического давле­

 

Ср

 

 

 

t,cym

ния подземных вод

на крепь

 

 

 

 

ствола.

В

массиве,

облада­

0,05-

 

 

 

 

 

 

Г 5

ющем линейной

наследствен­

0,1\

 

 

 

 

 

 

 

 

ной ползучестью, образуется

0,

2-

 

 

 

г 10

вертикальная

выработка

0,5\

 

 

 

 

 

 

15

круглого

сечения,

причем

11

 

 

 

 

 

 

20

 

 

 

 

 

 

 

 

 

Z-.

 

 

 

\-30

 

 

 

 

 

 

 

 

 

з -

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

W

 

 

 

 

 

 

 

 

 

10 -

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

50

 

 

 

 

 

 

 

 

 

2Q-.

 

 

 

£-7У

 

 

 

 

 

 

 

 

 

30 z

 

 

 

 

 

 

 

 

 

 

 

 

я н

 

 

 

100

 

 

 

 

 

 

 

 

 

Рис. 23. Номограммы для определения

функции

 

 

 

 

 

 

 

 

 

ползучести Ф при Абелевом ядре в зависимости от

 

 

 

 

 

 

 

 

 

 

 

 

 

времени:

 

 

 

 

 

 

 

 

 

 

а I =

12 ч; б — t — Ь сут‘, в i =

100 сут

 

 

 

 

 

 

 

 

 

70

Соседние файлы в папке книги