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

книги / Механика деформирования и разрушения горных пород

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

Увеличение бокового давления не приводит к изменению положения J3P на оси абсцисс и, естественно, не влияет на положение Атах. Но оно существенно влияет на количественную сторону процесса. Так, с увеличением ст3 от 0 до 10 МПа число импульсов АЭ сначала уменьшается от 22,5*103 до 2*103, а затем вновь возрастает до 35-103 при ст3=100МПа. При дальнейшем увеличении а 3 значение Nmax уменьшается. Зависимость 7V =/(a3) рассмотрена ниже. Влияние бокового давления на ширину области, где образуются трещины отрыва, обнаружить не удалось.

Приведенные выше результаты относятся к случаю ис­ пытания образцов при контактных условиях I (см. рис. 1.12, ав) на торцах (сухое трение). При этом зона максимального разрыхления расположена в центре деформированного образца, который в данном случае принимает форму «бочки». При испытании образцов с низкомодульными прокладками [кон­ тактные условия II (см. рис. 1.12, г—е)] «бочка» не образуется, и образец сохраняет форму цилиндра от начала до конца процесса деформирования. Качественно данные зависимости совпадают с рассмотренными выше. Отличия носят количест­ венный характер, причем большинство кривых лежат ниже своих аналогов на рис. 1.12, ав.

Указанная закономерность нарушается только в области боковых давлений ст3 ^ 20 МПа. Это вызвано, по-видимому, недостаточной стабильностью в выполнении условий соблюде­ ния равномерности напряженного состояния по высоте образца при малых боковых давлениях. В этой же области боковых давлений наблюдаются и наиболее сильные отличия в характере разрушения образцов с контактными условиями I и II. Так, при сухом трении на торцах образца визуально отмечаются наклонные магистральные трещины. Использование низкомо­ дульных прокладок приводит к образованию отрывных магист­

ральных трещин субвертикальной ориентации.

При давлении

а 3 > 20 МПа различий в

количестве

и углах

наклона линий

Чернова— Людерса для

указанных

контактных условий не

обнаружено.

 

 

 

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

ется к

пластическому.

Из рис. 1.13, в— г видно, что при

давлении

а 3 > 10 МПа

коэффициент остаточной поперечной

деформации ц не превышает 0,5—0,55, а остаточное увеличение объема материала не более 1,3% при достигнутой осевой деформации 12%. Характерная особенность для соли — резкая

Рис. 1.13. Акустическая

эмиссия при испытании образцов

каменной соли.

Контактные условия II

на торцах:

 

а — а 3 = 0; ' — ст3 = 2,5 МПа;

вст3=10М П а; г—а3 = 50МПа

i

зависимость числа генерированных в образце сигналов и ин­

тенсивности

АЭ

от уровня

ст3.

Так, N

при

давлении

а 3 = 50МПа

более

чем на

четыре

порядка

меньше,

чем

в случае, когда а 3 = 5 МПа. Этот факт указывает

на то,

что

при пластическом деформировании материала число генериру­ емых сигналов АЭ в образце незначительно.

Интенсивность АЭ для соли подчиняется тем же закономер­ ностям, что и для мрамора, исключая некоторые особенности, присущие данному материалу. Также, по мере роста а 3, выделяют два характерных максимума интенсивности АЭ. Первый из них тесно коррелирует с областью перехода 9Р через ноль, однако с повышением бокового давления интен­ сивность АЭ в этой области резко снижается.

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

При давлении 5— 10 МПа интенсивность АЭ после перехода через первый максимум значительно снижается. В этой области

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

При давлении более 10 МПа начинает формироваться второй максимум на кривой N / Nmax= f ), пик которого более размыт и находится в области больших осевых деформаций (6^10% ). Причина образования второго максимума, видимо, заключается в обусловленной технологией изготовления неодно­ родности центральной части образцов, которая проявляется при больших значениях гх и сг3. В этом случае в центральной части образца локализуется зона образования микротрещин отрыва.

В пользу данного предположения свидетельствуют резуль­ таты, полученные при деформировании соли с сухим трением на торцах. Понижение степени однородности силового поля в образце увеличивает второй максимум интенсивности АЭ, одновременно уменьшая первый. Как и для мрамора, при повышении бокового давления просматривается тенденция уменьшения значений N и Nmax для контактных условий II. Однако для соли такая тенденция выражена менее отчетливо, что может быть связано, как указывалось выше, с неоднород­ ностью центральной части образцов.

Исследования закономерностей поведения АЭ при сложном напряженном состоянии показали, что характер зависимостей числа сигналов АЭ существенно зависит от соотношения компонент напряженного состояния образца.

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

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

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

1.4. МОДЕЛИРОВАНИЕ ПРЕДЕЛЬНОГО ДЕФОРМИРОВАНИЯ И РАЗРУШЕНИЯ ГОРНЫХ ПОРОД

Вмеханике твердого тела существуют два основных подхода

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

большой интерес представляют модели, осуществляющие ком­ бинированный подход.

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

Статистическая модель деформирования горных пород, объясняющая их механизм разрушения с образованием пло­ щадок сдвига и трещин отрыва, приведена в работах [31, 32].

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

находятся

на конечном расстоянии друг от друга, дал

Т. Екобори

[12]. При построении модели для горных пород

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

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

позволит получить информацию о

поведении напряжений

и перемещений в окрестности особой

точки.

Решим сначала несколько более общую задачу: пусть имеются два клина раствором а и ~Р, разделенные полосой

скольжения, совпадающей

с осью

Ох (рис. 1.14, а).

Однородные граничные

условия

имеют вид:

а п = а т = 0 при 0= а; а п = стт = 0

при

0 = 0;

а т+= а т~=0; [г] = [ а п] = 0 при

0 = 0

(1.30)

или в функциях Колосова — Мусхелишвили (pt , v|/1? /=1; 2:

34

а

6

в

Рис. 1.14. Схемы

расчета:

 

б— расположения контура интегрирования;

в-т-

а—расположения

полосы скольжения;

расположения полосы скольжения

и трещины отрыва

 

 

 

 

 

ср1+гф'1 + 1р,1 =0 при

0 = а;

ф2 + гф'2 + ф2 = 0

при

9 = (3;

 

 

1т(фг+ 2 ф[ + ф,.)= 0

при 0 = 0,

/=1,

2;

 

 

 

Ке(ф1+ 2 ф'1_+\р'1)= Ке(ф2 + 2 ф'2-|-ф2) при

0 = 0;

 

где

1ш(хф1-г ф 1 - ф 1 ) = 1т(хФ2-2ф2-ф2)

ПРИ

е = 0>

0-31)

х = 3—4v; z = x + iy;

v— коэффициент

Пуассона.

 

 

 

Для нахождения собственных чисел и собственных решений

используем преобразование

Меллина

 

 

 

 

 

 

 

 

/ =

ï / z x_1dr,

 

 

 

 

где

.А— параметр.

 

 

0

для аналитических функций

ф,-,

v|Ij

Можно показать [21 ], что

при J = l;

2:

фj = F je-k‘a,

фу.= (7;е “м ,

 

 

 

 

 

 

 

 

 

где Fj=Aj+iBj;

Gj=Cj+iDy

Aj, Bj, Cj,

D}— коэффициенты,

в общем случае

комплексные, зависящие от А.

 

 

 

 

Запишем преобразования Меллина для граничных условий.

 

Из граничного, условия

ф1+ 2ф'1 + ф1=0 при

0 = а

 

 

(zl1 + /5 1)e " iXe-A (/l1- / 5 1)ei*a + 2,+ (C1+/Z)1)e " ,4e= 0.

 

 

Из условия ф ^гф 'х + ф ^ О при 0 = а

 

 

 

 

 

(Лх —ii?! )е -ае —A(/lj +

)e -il,(>’+2)+ (C1 —/.Dj )е -,Хв= 0.

 

Аналогично поступаем с другими граничными условиями. Рассмотрим случай, когда а —Р=2л и полоса скольжения

переходит в разрез, моделирующий трещину (рис. 1.14,6). Мржно показать, что условие существования нетривиаль­

ного решения системы линейных уравнений для нахождения коэффициентов AJt В}, Cj, Dj сводится [21] к нахождению собственных чисел уравнения

[A2 cosAz( 1 —cos* )+ A sin Az sinx—cosAz+ cos 2л A] sin2л A = 0, (1.32)

где л' = 2ос; z= a + P= 2(oc—я).

Из уравнения (1.32) следует, что спектр его решений

распадается на две части:

 

\ =к / 2, Аг= 1; 2,

 

X2cos^z(1—cosA') + Xsin^zsin.Y —cos^z + cos2À,z = 0.

(1.33)

Численное решение уравнения (1.33) методом Ньютона на комплексной плоскости показало, что при всех допустимых значениях параметра X для всех нулевых корней выполняется неравенство |Re?t|>0,5. Следовательно, первым ненулевым решением уравнения (1.33) при всех а является Х= 0,5, и асимптотика решения определяется однородным решением, соответствующим этому значению.

Найдем нетривиальные решения однородной системы линей­ ных уравнений. В общем случае можно зафиксировать какой-

либо коэффициент и

выразить

через него все

остальные [21]:

 

Fx= (1 + Х—À,e2l0t—eL2aot)(À,sin2p + sin2À,|});

 

 

F2 = ( 1+ Х—Хе2,р—e2ap)(^sin2oc + sin2^P);

 

G1=(1 + A.)(1 —X+ Xe~2ia—e2Xai)(X.sin2p + sin2X.p);

G2 = ( 1+ A.) ( 1—}.+ ^ e “ 2,p—e2xp‘)(À.sin2a + sin2^a).

(1.34)

Найдем теперь однородные решения. Из обратного преоб­

разования

Меллина

 

 

 

 

 

 

 

сря = 1/(2тс/) J <р*г

l dX,

 

(1.35)

где г— радиус-вектор,

ц

L

проходит так, как

показано

а контур

на рис. 1.14, б.

что Fj

и

Gj

таковы,

что

интеграл

Можно

показать,

в выражении (1.35) допускает

разложение по

вычетам

 

<P* = Z<PR(k/,

9)г

Xi-

 

 

 

 

?..

 

 

 

 

 

Ограничимся здесь случаем, когда XteR. Тогда можно убедиться, что слагаемое ср*, соответствующее л;, принадлежит

R. Заметим,

что для XteR в общем случае ( p ^ R ,

и здесь

требуется дополнительное

рассмотрение.

j= 1; 2,

Итак, принимая во внимание смысл Fj, Gj,

окончательно

получим:

 

 

 

ФJ = CFJ (X)Z \

 

(1.36)

 

* J = CGJ (X)Z \

(сеЛ; у=1; 2).

 

 

Исследуем теперь первое

собственное решение при Х= 0,5.

Из формул (1.34) и

(1.36)

найдем:

ф1= ф 2

= ф = (0,5е2,с‘- е ia + 0,5)v/z;

Ф, = Ф2 = Ф = ( l,5

-0 ,5 e -2ia- e ia)0,5yfz.

При z = r eR имеем

ср + гф ' + ф= v/r( 1,5 + O,5cos2a —2cosa)e R, следовательно, на линии скольжения касательные напряжения

равны 0. Нормальные напряжения

имеют вид

 

 

а л = —^ ( 1,5+ O,5cos2a —2cosot).

(1.37)

 

2 ф

 

 

Найдем раскрытие берегов трещины. По формуле Коло­

сова — Мусхелишвили

 

 

2G(u +

iv) = y.<p— zip' —ф = (и + 1)ч/ г [ —ехр( —/а/2)+

 

+0,5ехр(/а/2)+0,5ехр(5а//2)];

мг+ ш в = (ы+ /г)ехр( —/а).

Тогда

 

 

 

26(мг+/ив)= (х + 1) , / г [ —ехр (-3 /а/2 )+

 

 

+ 0,5ехр( —/а/2)+0,5ехр(3/а/2)].

 

Отсюда следует

 

 

2GM, = (X + l)N/ r [l,5sin(3/2a) —0,5sin(a/2)].

(1.38)

Аналогично

при 0 = a —2л

 

 

2Gue= —ч/г(х + l)(l,5sin3/2a —0,5sina/2).

(139)

Проанализируем механический смысл формул (1.37)—(1.39).

Пусть осе [л,

Зл/2], тогда

 

 

 

стп = —- (1,5 + 0,5 cos2a —2 cosa) > 0.

 

_

2V^

 

 

Следовательно, на линии скольжения имеет место состояние растяжения. При моделировании растягивающее напряжение принято положительным. Рассмотрим соотношение для пере­ мещения щ

2Смв/[(1 + к )ч/г ] = 1,5 sin3/2a —0,5 sin 1/2а.

Это

выражение___ обращается

в

ноль

при

а = а* = 2(л —a r c s i n 2/3)»250° Можно

показать, что

на ин­

тервале

[я, ос*) величина и0 > 0, т. е.

трещина

раскрывается

и на интервале (а*, Зл/2] ме > 0 берега трещины перекрываются. Аналогично, если Y < 0, т. е. при сжатии на линии скольжения, при а > а* берега трещины пересекаются, при ос > а* расходятся. В горных породах при сжатии на полосах скольжения имеем ап < 0, случай перекрытия берегов трещины механически невоз­ можен, следовательно, возникновение трещины отрыва воз­ можно лишь при ос > а*

Таким образом, максимальный угол между полосой сколь­ жения и трещиной может принимать значение 2л —ос*

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

характеризуются

одинаковыми геометрическими параметрами,

а. именно: L — длина полосы скольжения;

/— длина

трещины;

б— угол среза;

у — угол между линией

среза и

полосой

скольжения (рис. 1.14, в); <р — отклонение направления трещины от вертикали. Первоначально будем полагать, как и в ста­

тистической

модели

[31],

что

трещины располагаются вер­

тикально, т. е. угол

ср = 0.

Полная

работа, затрачиваемая

на

деформацию

тела,

А = Ау + Ап + Ат,

где Аг— работа упругих

сил; Ап— работа

пластических

деформации;

Ат— работа,

за­

трачиваемая

на

образование трещин.

по формуле

 

Работа А

главных напряжений

находится

 

 

 

 

Л = П

I

a.deidn.

 

 

 

 

 

П

i=l,2, 3

 

 

 

Поэтому А может быть вычислена, если известны диаграммы напряжение— деформация.

В предположении линейной упругости материала Ау также может быть вычислена* из зависимости напряжение— дефор­ мация:

Ay= -]I.efE .

 

Работа пластической деформации

 

An = K L n [ u l

(1.40)

где п — число элементов деформации;

[и ]— среднее смещение

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

кого формоизменения. Если принять, что толщина полосы скольжения Л «с 1 и на ней действует условие пластичности Треска с параметром тп, то К к \ п. В качестве тп примем наибольшее касательное напряжение на пределе упругости.

Работа трещинообразования

Ат= 1пр, где

р — плотность

поверхностной энергии, которая

для горных

пород мала.

В частности, это находит свое выражение в том, что их хрупкая прочность — прочность на растяжение — мала. Поэтому можно пренебречь Ат по сравнению с Ап и А .

Рассмотрим остаточное изменение объема 0. Из элемен­

тарных геометрических соображений (см.

рис. 1.14, в)

0р%/[м] sin (л —5—у)п.

(1.41)

В осесимметричном случае остаточное изменение объема

0р = е?-2е$ .

Из оценок выражений (1.40) и (1.41) получим

Ay)/D = АХ/[/sin (л —5—у)].

Отсюда следует

// L = Я / [sin (л —8—у)],

где H = D K / ( A - A y).

Из треугольника на рис. 1.14, в имеем

 

//L = siny/sin(n —5 —у);

 

 

sin (5 + у) sin у = sinа Я.

 

После

преобразований

 

 

cos(2y+ 6) = cos8 —2#sin5.

(1.42)

Из соотношения (1-42), зная Я и 8, определяем угол между

трещиной

и полосой скольжения.

 

Оценим соответствие теоретических результатов с экспе­ риментальными данными при больших и малых значениях параметра c = a 2l a i- Объемное расширение 0Р при достаточно большом значении параметра с мало, тогда мало и Я. Следовательно, cos (2у + 8) = cos 8. Отсюда следует, что 2у + 8 = 8 и у = 0, т. е. угол среза совпадает с углом, под которым расположены линии скольжения. Трещины отрыва при этом не образуются и теоретические и экспериментальные результаты совпадают.

Оценим соответствие расчетного и экспериментального угла

при

малых с

для мрамора I, талькохлорита и мрамора II.

При

с= 0,069

для мрамора

1

[31] имеем

 

 

 

0рж4 • 10-7 м -3 ;

А — Аул 1,705 • 105 Н/м2;

 

 

 

5

= 30°; £ = 5 8 0 -105 Н/м2,

 

 

где

значения

0 и А — Ау отнесены к

единице

объема. Отсюда

следует, что

у = 45°

и (3 = 75°

II

при

с = 0,069 имеем

Р= 71,5°

Аналогично для

мрамора

и для талькохлорита при с= 0

Р«70° При

действии

сжима­

ющих напряжений возникновение трещин отрыва возможно,

как установлено выше,

если

а > а* = 250°

Для мрамора I,

мрамора II й талькохлорита соответственно а = 255°; а = 251,5°;

а = 250° Следовательно,

при

сжатии пород

происходит рас­

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

что оценку выражения (1.36) нельзя удовлетворить ни при

каком

значении угла у.

Это связано с тем, что развитие

полос

скольжения сильно

искажает напряженное состояние

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

1.5. УРАВНЕНИЯ ПЛАСТИЧЕСКОГО ДЕФОРМИРОВАНИЯ ГОРНЫХ ПОРОД

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

На практике получил распространение второй подход, основанный на теоретическом обобщении экспериментальных данных и формулировке уравнений, описывающих главные

черты процесса пластического

деформирования.

Механизм пластического деформирования

горных пород

как неоднородных тел изучен

в работах [31,

32]. Развитие

пластического деформирования на основе статистической мо­ дели определяется числом плоскостей скольжения, по которым происходит сдвиг. С ростом числа плоскостей скольжения пластическая деформация пород увеличивается. Для ее описания предлагается единая кривая пластической деформации [31 ]

ер = ф(т - т«еВс),

где т®, В — постоянные; eÇ — главная пластическая деформация; Ф — известная функция.

Запись ‘единой кривой пластической деформации в коор­

динатах

Дт = т —т°е , в?

вызвана тем, что

предел упругости

горной

породы зависит

от показателя

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

в качестве второй координаты главной пластической дефор­ мации eÇ, а не сдвига eÇ —е§ вызвано тем, что из-за раскрытия

микротрещин зависимость между

eÇ и

cÇ —е*> нельзя пред­

ставить в виде единой кривой для

всех

значений с. Единую

кривую деформации можно описать с помощью

уравнения

г? = ср(т,

с).

(1.43)

С учетом связи между с и т,

ст уравнение (1.43)

примет вид

е? = <р(т,

ст).

(1.44)