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

книги / Механика горных ударов и выбросов

..pdf
Скачиваний:
11
Добавлен:
19.11.2023
Размер:
26.32 Mб
Скачать

Общие представления о полных диаграммах, приведенные выше, слу­ жат для выделения характерных уча­

стков деформирования горных

пзрод

и качественного изучения их

особен­

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

имеют горизонтальные участки (рис. 5), то частичное описание дается теорией предельного состояния (деформации в ней не ис­ пользуются и статически неопределимые задачи не рассматрива­ ются [54, 56, 61]). Эта теория может привлекаться и при рассмо­

трении

участков остаточной прочности (горизонтальные

участка

за максимумом нагрузки на рис. 3).

 

 

 

 

 

 

 

ют

Для развития динамических явлений важнейшее значение име­

и

наклонные, падающие участки

запредельных

кривых. Их

а

 

 

 

описание начато

сравни­

 

 

<Vбр

тельно недавно и в про­

 

 

 

 

 

 

стейшей

форме представ­

 

 

 

 

лено

в

следующем

под­

 

 

 

 

разделе.

Здесь

же

оста­

 

 

 

 

новимся

кратко

на

пове­

6,

 

 

 

дении пород при действии

 

 

 

 

растягивающих

напряже­

 

 

 

 

ний, а также при совме­

 

 

 

 

стном действии сжатия и

 

 

 

 

растяжения

по двум

вза­

 

 

 

 

имно

перпендикулярным

РИС. 6 .К оценке прочности при совместном

направлениям.

Эти

све­

дения

 

необходимы

для

действии растягивающих (аз) и сжимающих

 

(ffi)

напряжений:

теории

выбросов

 

угля

а —схема

нагружения; б —'диаграмма критических

(породы) и газа.

горных

сочетаний <т( и

ôj, отвечающих разрушению

Особенностью

прочность

 

пород является их низкая

на растяжение ор по сравнению с прочностью на сжа­

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

рушения [24, 53, 69]. Приложение сжимающей нагрузки в на­ правлении, перпендикулярном к растяжению (рис. 6,а), как пра­ вило, уменьшает отрывающую силу. Сочетания критических на­ пряжений в целом отвечают диаграмме, изображенной на рис. 6,6, и могут быть аппроксимированы формулой

-ÏÜ-— ïi—= 1 •

(1.5)

ff0

 

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

1.2. ЗАПРЕДЕЛЬНЫЕ ДИАГРАММЫ

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

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

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

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

максимумов o i= /(o з)> отвечающими паспорту прочности, можно записать в виде

a i = F (а3, si) = / (а3)+М 0ф

(а3, etm— ei),

(1.6)

где М0— неотрицательный множитель

с размерностью

напряже­

ний; <р(а3, Eim—? i) — безразмерная функция, определяемая видом падающих участков; eim— Деформация в направлении первой глав­

ной оси напряжений при

боковом сжатии <т3

в момент

перехода

к падающему участку.

поскольку при e i= e im точка принадле­

Ясно, что ср(а3, 0 )= 0 ,

жит паспорту прочности.

Функция ф неотрицательна,

так как

<П, f(o3) < 0, а с ростом

|sim—ei| значения

|oi| не возрастают.

По той же причине дф/д(б1т— 8 i)^ 0 . При

достаточно

больших

|Elm — ei| все запредельные кривые становятся

горизонтальными.

При этом

 

д(р!д (eim— 81) =дф /до3= 0 .

(1.7)

Эти же равенства приближенно выполняются и вблизи от то­ чек максимумов. При а3= 0 и больших |eim— &i | имеем ai— —а*, т. е. ф (0, оо) = (ао—а*)/М.

Случаю М0= 0 отвечает обычная теория предельного состоя­ ния. Кроме того, в начале и в конце плоской запредельной дефор­ мации имеет место совпадение характеристик системы уравнений теории предельного состояния с характеристиками теории запре­ дельного состояния. Чтобы убедиться в этом, достаточно присое­ динить к (1.6) уравнения равновесия и учесть, что в начале и в конце запредельной деформации из (1.7) следует dci=daz{dfldaz). Этот факт поясняет наблюдаемое иногда замечательное согласие углов, под которыми пересекаются поверхности разрушения, с уг­ лами пересечения характеристик, вычисленных по теории пре­ дельного состояния.

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

 

е3— G (а3, ei),

 

(1.8)

которое

получается при обработке замеров поперечных деформа­

ций в

опытах, проводимых в условиях

плоской

деформации.

В случае, если опыты проведены при условии <02=03

(е2= е 3), то

соответствующие поперечные деформации

при построении G—

= С (о3, 81) следует удвоить. Если нее имеются данные об измере­ нии объема 0 , то искомую зависимость можно получить, полагая

83= @ — 81.

Неизвестными при решении задачи являются ai, о3, ei, е3 и угол #i, образуемый первым главным направлением с осью Ох. Присоединение (1.6) и (1.8) к уравнению сплошности и двум уравнениям равновесия, в которых компоненты тензоров выражены через главные значения и угол fti, дает систему пяти уравнений

 

 

 

e3 =

G=

 

 

 

E

 

 

 

 

 

при

 

 

4 - [ ‘ +

P8(l + 4

)

 

 

 

 

 

 

 

 

 

при -|- +

|3l-£i <

k l < k

l + PL^ (1.12)

E [Н - Р 8 ( Ч - - § - ; *. — s ( 1 +

T ' ) ( s» + P 'ir + 'è ')—

- « . +

« . +

P #

при |«.|>|».| +

Plx

:

 

 

 

 

 

при

 

+

I *3

 

ф =

 

 

 

(1.13)

 

 

 

 

 

 

 

 

 

 

4 - (° з + ^ )

ПРИ 18,1 > 1 г + Р т г

Помимо

уже

встречавшихся

величин,

здесь

 

ô — абсолютная

величина отношения приращения поперечной необратимой дефор­ мации к приращению осевой необратимой деформации на падаю­ щем участке запредельной кривой; ô = 2ôo, где выражения ôo че­ рез величины, получаемые в опытах на одноосное сжатие, даются

формулой

(1.4); р = Г—sîn^"^ Р — угол тРенияПервым строчкам

в (1.11),

(1.12) соответствуют упругие Деформации, вторым — де­

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

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

Линейность зависимостей (1.11), (1.12) существенно облегчает решение задач и, в частности, использование метода конечных элементов, когда известны направления главных осей. В ряде слу­ чаев решение получается в квадратурах (ем. ниже). Теории пре­ дельного состояния для плоской задачи отвечает присоединение к двум уравнениям равновесия второго из соотношений (1.11) при М — 0. В случаях, когда остаточной прочностью можно пренебречь, в (1.11), (1.12) полагается а * = 0.

Для решения задач и моделирования удобно .нормировать на­ пряжения на прочность при одноосном сжатии ст0, а деформации на отношение прочности к модулю упругости aofE. Отмечая нор­ мированные величины штрихами, имеем:

= о ,Ч ;

o'2= a j a 0; о,3= о3/о0; О* ---®

8', =

е'2 = е8£/о0; е'3 = «,£/о.; в р= в р£,/о,; (1.14)

F'(o'„ e,,) = F(o0o'3, s\oJE)joa\ Ф' (o'„ е',) = Ф ( e / 3, е\а0/£ ) /в.;

G'(s'3, 0 = С(з0о'., е\о*(Е) Efa0,

где F , Ф', G' — нормированные функции, которые получаются из функций F, Ф, G, входящих в формулы (1.6), (1.9), (1.10) и опре­ деляемых экспериментально. Для кусочно-линейной аппроксима­ ции (см. рис. 7), описываемой формулами (1.11) — (1.13), зависи­ мости F', Ф', G' для плоской деформации имеют вид:

 

 

 

л

 

|«',1<Ч -З К 1

P (I +

*) -

*«', - 1 -

1 + Р13', |<

|е',| < |е'.»|+р|о"1

- o ' . +

fb'.

 

 

 

K H « 'J + f t > '.l

ф, | 0

 

К,1<1 + РК,|

(1-15)

 

W , + F')/2

|е', |>1+р|а',|

 

 

 

 

 

К И 1 + Р К 1

[ 1 + Р « ( Я + 1 ) ] а ',- 8 ( 1 + 1 ) ( .'1+ 1 )

 

[1 —j—pô (Я —

1)] о '3 — 5 (Я 4 -

1) (е'* - f- (Ь'з + 1 ) —

В случае, когда опыты проведены при условии 02= 03(62= 83),

функция G' представляет удвоенную нормированную

поперечную

деформацию (е/2= е 'з), a F'

сохраняет

смысл

нормированного

осевого напряжения о'ь

 

 

и G1 получаются из гра­

Графики нормированных функций F'

фиков исходных функций F и G изменением масштабов координат-

 

ных осей. Для кусочно-линейной

 

аппроксимации

нормированные

 

диаграммы

изображены на рис.

 

8.

При

моделировании

после

 

нормировки

нужно

 

подобрать

 

материал, у

которого

аналогич­

 

ным

образом

нормированные

 

диаграммы

близки

к

получен­

 

ным. Лишь такой материал, стро­

 

го говоря, можно считать экви­

 

валентным, так как вид запре­

 

дельных кривых, как будет пока­

 

зано в дальнейшем, влияет на

 

устойчивость выработок.

урав­

 

Использование полных

 

нений для запредельных деформа­

 

ций возможно только для с т а т и ­

РИС. 8. Кусочно-линейные нор­

ч е с к и х

задач,

поскольку

запи­

сать

падающую

диаграмму

при

мированные диаграммы для пло­

ской деформации

динамических

нагрузках

прин­

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

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

«Строгий» ответ на этот вопрос, если не затрагивать структуры материала и условий проведения опытов, такое: ни в какой мере; запредельная деформация — это мистика; говорить о ней бессмыс­ ленно. В этом состоит так называемый динамический парадокс.

Однако экспериментатор, проводя опыты в некотором диапазо­ не скоростей нагружения, вовсе не ощущает физической «бес­ смысленности» падающего участка. На жестких прессах диаграм­ мы хорошо воспроизводятся и изменяются сравнительно слабо при скоростях от 10-8 до 1 м/мин. Такие скорости часто реализу­ ются на практике, и соответствующие диаграммы полезны и важ­ ны для практических целей: было бы наивно думать, что простой перенос образца из шахты в лабораторию влияет на его поведе­ ние в заданных условиях нагружения *.

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

приводится в 4.5.

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

* Масштабный эффект здесь не обсуждается.

2— 133

Вывод о локализации разрушения, начиная с некоторых ско­ ростей нагружения, подтверждается экспериментальными данны­ ми об изменении характера * деформации [8]. Начиная с неко­ торой критической скорости наблюдается резкая локализация дробления вблизи от места приложения нагрузки; однородной деформации в образце нет. Значительно раньше существование критических скоростей было обнаружено шахтньши наблюдения­ ми [46]. Из них следует, что влияние скорости внедрения в ударо­ опасный пласт чрезвычайно велико и приближение ее к некоторым предельным значениям сопровождается как увеличением числа упругих импульсов на единицу подвигания добывающей машины, так и резким возрастанием их интенсивности. Деформирование при этом локализуется около места внедрения и принимает очень бурный характер. Дело обстоит так, как если бы скорость рас­ пространения необратимой деформации в глубь пласта была ог­

раничена.

Это дает основание говорить о «соответствии» или «несоот­ ветствии» скорости нагружения скорости выхода пласта из-под нагрузки [14, 46, 47].

Теоретическое определение критических скоростей представля­ ет очень трудную, интересную и пока .не решенную задачу. Их практическое нахождение выполняется по наблюдениям за измене­ ниями в характере деформаций, которые описаны выше. При ис­ пытаниях на прессах с сервоконтролем достижение критической скорости может отмечаться по переходу от диаграмм первого ти­ па (см. рис. 3) к диаграммам второго типа (см. рис. 4). Послед­ ние не могут быть получены при поступательном движении за­ хватов пресса, требуют возвратного движения, т. е. изменения скорости не только по (величине, но и по направлению. В шахт­ ных условиях такой процес практически не реализуется, и дости­ жение скорости, при которой совершается переход к диаграммам второго типа, вне зависимости от жесткости нагружения сразу ве­ дет к бурному дроблению породы. При скоростях, превышающих критическое значение, модуль спада можно считать неограниченно большим = оо).

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

мулах, описывающих свойства

пород,

величина

сто, полученная в

лабораторных опытах, заменена величиной

о0£> учитывающей

коэффициент длительной прочности

и коэффициент структурно­

го ослабления £с:

 

 

 

Д0;

У»/®О*

 

 

* В отличие от характера разрушения максимальные нагрузки сравнитель но слабо зависят от скорости нагружения [8, 57].

При учете фактора времени будем также во всех формулах, вместо модуля упругости Е и коэффициента Пуассона v использо­ вать 'Переменные модули Eu v* [}5], сохраняя, однако, неизменным отношение К ‘модуля спада к модулю упругости 'к=М1Е (в на­ стоящее время отсутствуют надежные экспериментальные данные- о зависимости этого отношения от скорости нагружения). Разу­ меется, мгновенные значения Et и v* совпадают с £ и -v. При нагружении со скоростью, превышающей критическое значение,, модуль спада считается равным бесконечности. Для того чтобы неусложнять обозначений, ниже будем применять прежние символы (а0, Е, v, М), имея при этом в виду сделанные замечания о по­ правках к значениям, полученным в лабораторных экспериментах при обычных прессовых скоростях нагружения. Случаи, когда экс­ периментальные данные лабораторных опытов используются непо­ средственно, без поправок, оговариваются.

1.3. ВЛИЯНИЕ ЗАПРЕДЕЛЬНЫХ СВОЙСТВ

Рассмотрение задач для областей/в которых действуют сжи­ мающие напряжения, сводится к решению полной системы урав­ нений, описанной в предыдущем подразделе при тех или иных граничных условиях. Соответствующие задачи .могут быть решены численными методами (например, методом конечных элементов). В настоящее время уже разработаны первые программы, которые позволяют выполнять расчеты для достаточно произвольных свойств пород. Тем не менее для теории важно не только уметь численно решать разнообразные задачи, но и располагать замк­ нутыми аналитическими решениями, позволяющими в легкообо­ зримой и компактной форме рельефно выделить особенности и 'упростить исследование основных закономерностей, свойственных, той или иной модели массива горных пород.

Зачастую подобные решения оказываются вполне приемлемы­ ми и с точки зрения практического использования расчетов, по­ скольку сложность и известная неопределенность геологического строения массива и ряд трудно учитываемых технологических фак­ торов нередко делают иллюзорными результаты более «точных» расчетных схем. Поэтому, например, современная практика расче­ та устойчивости капитальных и подготовительных выработок су­ щественно опирается на рассмотрение простейшей, решаемой в. замкнутой форме задачи об осесимметричной деформации пород вокруг цилиндрической выработки [10]. Учет факторов, остаю­ щихся за рамками схемы, таких, как неравномерность нагруже­ ния контура крепи, неоднородность строения пород, наличие на­ рушений, негидростатическое состояние на бесконечности и так далее, осуществляется с помощью введения поправок к базовой схеме. К этому следует добавить, что влияние формы контура и нагрузок на нем быстро затухает по-мере удаления от выработки и замена квадратного, трапециевидного или иного сечения на рав­ ное по площади круглое сечение сравнительно мало сказывается на результатах ten .

2 *

19

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

 

 

ческого

явления

при

выходе ее

ч

 

из строя.

Аналогичным

образом

/

получаются

результаты

и

для

сферической

выработки.

 

Оста­

 

 

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

 

 

дробно

не

будем,

ограничившись

 

 

замечанием,

что

в

 

дополнение

 

 

к приводимым ниже данным, ре­

 

 

зультаты для сферической

выра­

 

 

ботки

свидетельствуют

о

неко­

 

\

тором

повышении

устойчивости

/

у забоя и позволяют учесть изме­

 

нения

в

нагрузках,

 

происходя­

 

 

щие по

мере его

подвигания.

РИС. 9. Схема к задаче о цилиндри­

В целом же случай сферической

симметрию не отличается в каче­

ческой выработке

или скважине

ственном

отношении

от

случая

 

 

осевой

симметрии.

 

 

 

 

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

ния динамических

явлений, осуществляемого с помощью бурения

скважин.

 

 

 

круглого

сечения е радиусом R

при

Рассматривается выработка

гидростатическом

сжатии

давлением

q

на

бесконечности

(рис. 9). Уравнения равновесия

и совместности

деформаций име­

ют вид

 

 

 

 

 

 

 

 

± 1 .+

г

= 0 ,

* ^ + '- 1 ^ = 0

(1.16)

 

dr

'

 

dr '

г

 

 

где

г — расстояние

от

центра

выработки

до

рассматриваемой

точки; оь <тз — нормальные напряжения на площадках, соответст­

венно параллельных радиусу и перпендикулярных к нему; еь ез — относительные деформации в направлениях, соответственно перпен­

дикулярных

и параллельных радиусу: ез=dujdr, г\=и(г\ и

смещения по

радиусу (они отрицательны, так как происходят по

направлению к центру выработки).

Вдали от выработки а3(оо) = —Ц\ на ее контуре действует дав­ ление р, т. е. о з (Я )= —р. В случае свободного контура сг3(i/?)=0,