книги / Механика горных ударов и выбросов
..pdfОбщие представления о полных диаграммах, приведенные выше, слу жат для выделения характерных уча
стков деформирования горных |
пзрод |
и качественного изучения их |
особен |
ностей. Однако решение задач теории динамических явлений и других за дач горной геомеханики требует не только качественного познания меха нических свойств, но и количественно го описания различных участков пол ных диаграмм, получаемых в экспери ментах с горными породами в услови ях сжатия. Такое описание допредель ных участков осуществляется с по мощью теории упругости и вариантов теории пластичности. Если диаграммы
имеют горизонтальные участки (рис. 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,